24小时热门版块排行榜    

查看: 1522  |  回复: 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的回帖

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的回帖
查看全部 3 个回答
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:39 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:03 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 00:51 by 9i7PCSZyMzot
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:03 by 9i7PCSZyMzot
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 4pRyVGCVerAA 2026-09-19 4/200 2026-09-20 00:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:51 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:39 by 9i7PCSZyMzot
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-18 7/350 2026-09-19 21:03 by 9i7PCSZyMzot
[教师之家] 现代”学阀”该如何界定 +3 iaeyuan 2026-09-19 3/150 2026-09-19 20:52 by 没昵称呀
[基金申请] 国社科系统bug了,是不是要放榜了? +7 kynobel 2026-09-16 9/450 2026-09-19 14:52 by 接好运婷
[考博] 申博发邮件 +8 Lmengk 2026-09-15 9/450 2026-09-19 09:27 by 中间体超市
[有机交流] 各位大神,目前国内有哪些比较好用的逆合成软件? 20+3 zapen 2026-09-15 9/450 2026-09-18 18:17 by zapen
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
信息提示
请填处理意见