24小时热门版块排行榜    

查看: 1396  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 蒋士成 的主题更新
信息提示
请填处理意见