24小时热门版块排行榜    

查看: 1508  |  回复: 2

蒋士成

捐助贵宾 (小有名气)

[求助] VASP计算出来的BLOCH 波函数的相位问题。 已有1人参与

以下是我发给VASP-support 的邮件
Dear Sir or Madam,
 Recently, I tried to output the wavefunction Un (k, r) ( n is the energy band, k is the crystal momentum) calculated by VASP. As a complex wave function, Un(k, r) should be continuous in respect to k and r, physically. While, I found that when r is fixed, Un(k, r) is not continuous in respect to k. I mean, although the absolute value of Un(k, r) is continuous in respect to k, the phase of Un(k, r) is not (as shown in the figure below). I checked the self-consistency scheme of VASP, I saw that the wave function is initialized using a random number generator, and the convergence criteria is related to energy. I wonder it is the random number generator that introduces random phases into the wavefunction Un(k, r) at last. And, this random phase won’t influence the convergence standard. Can you tell me whether what I think is right? Anyway, I think the phases of the wavefunction are very important physically. So, I want to get some advices from vasp-support about how to deal with the random phases.

大概意思就是BLOCH 波函数Un (k, r) 是个复数,但是VASP计算出来的波函数的相位对晶格动量K 不是连续的,这在物理上是不对的。 但是VASP给我回复的邮件好敷衍,就两句话。
please note that WAVECAR only contains the wavefunctions COEFFICIENTS,
not the wavefuction itself.
the random number generator initiates some starting values of these
coefficients。 it has nothing to do with random phase.
他回复的东西完全没用嘛。有没有人比较精通VASP计算原理的,来说说这个random的phase到底来自于计算的哪些步骤。我觉得可能是来自于初始赋值以及迭代过程中矩阵对角化。 如何才能解决这个随机相位的问题呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

蒋士成

捐助贵宾 (小有名气)

没人?
2楼2017-03-04 07:28:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ra2ghgzh

木虫 (正式写手)

学痴

【答案】应助回帖

感谢参与,应助指数 +1
以下是我抄的:
为求得正确的总电荷密度分布,上式求解 uk(r) 时其k应该取越多越好,毕竟电荷密度是由波函数平方加总而来,每个k都有其贡献。但对于有N个晶胞的材料其总k点数也是N个,此N是10的23次方的数量级。所幸对于相邻近的k点,其所对应的能量本征值与波函数也都很类似,我们因此只要采用所谓k点取样(k point sampling)的技术,在k所定义的范围内尽量均匀公平的取点,如此也能呈现出合乎我们所希望之精确程度的电荷分布、系统总能,以及所有与波函数及其本征能量有关的物理量。
我永远当不了老板
3楼2017-03-04 16:39:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 蒋士成 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 科研孤儿太难了 +6 我4大白菜 2026-08-20 6/300 2026-08-20 17:25 by Dylan.Ke
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 时间戳又变了 +11 wuchongjun 2026-08-20 14/700 2026-08-20 16:26 by 虫友是什么虫
[基金申请] 今天系统多次维护,明天很可能放榜! +9 zju2000 2026-08-16 10/500 2026-08-20 15:56 by Leogzhya
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见