24小时热门版块排行榜    

查看: 677  |  回复: 1

lxbyf

铁虫 (初入文坛)

[交流] 【求助】GW修正中K点的问题已有1人参与

算能带的时候一直用侯博士的gk脚本来生成K点,
现在想在5.2中试一下GW的修正,就在能带计算的时候加入了GW0,但是运行时会提示:
two k-points are equivalent           1          16  this will cause problems in the HF routine
好像是说两个k点一样会造成问题,那应该如何处理,难道不是在计算能带的时候加入gw的参数吗?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gavinliu7390(金币+3):thanks! 2010-05-19 13:41:38
vasp5.2计算能带时,不能像以前采用4.6版本计算能带的方式来进行。
你在vasp上的官方论坛上找找看:
下面的是从上面拷贝过来的,针对HF的计算,可能对GW的计算也适用:
band structures using hybrid functionals have to be calculated the following way:
-) First perform a selfconsistent Hartree-Fock/HSE calculation using a conventional KPOINTS file.
-) Copy the IBZKPT file to KPOINTS, and explicitely add all desired k-points along high-symmetry lines of the BZ that are needed for the bandstructure plot. Add the points at the end of the KPOINTS file, but set the weights of these added k-points to 0.
-) Do not forget to set the number of k-points in KPOINTS correctly (to the number of the k-points used in the standard mesh PLUS the number of the k-points along the lines), such that all k-points are used for the calculations.
-) NKRED can not be used.
-) Perform a second VASP run:
It is recommended to use the Davidson algorithm, since it converges that eigen energies at the new k-points fastest. Since VASP terminates when the total energy is converged to a certain threshold, it is important to force VASP to do a minimum number of steps, so that the orbitals at the new k-points are fully converged (note: since their weight is zero, they do not contribute to the total energy).
This can be done using e.g.

ALGO = N ; NELMIN = 5 ! Davidson, minimum 5 scf-steps
IBRION = 1 ! Use simple charge mixer, since
Pulay might blow up

-) The KS-eigenvalues of the states along the high-symmetry lines are written in OUTCAR, EIGENVAL, vasprun.xml; please cut the k-points required for the bandstructure from one of these files and proceed as usual (using p4vasp or any other graphics tool you usually use to produce bandstructure-plots)

Note: A Hartee-Fock calculation can NOT be continued from an existing CHGCAR file, since the non-local exchange is not determined by the charge density but by the density matrix and/or the KS-orbitals.
2楼2010-05-19 13:41:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lxbyf 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见