| 查看: 2769 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求助vasp GW方法计算过渡金属错误CHECK_OCCUPATIONS 已有1人参与
|
|||
|
流程参考网上vasp gw计算方法帖子。 原帖地址：http://blog.sina.com.cn/s/blog_7c63531f0101ibx9.html 第一步的INCAR： Start parameter for this run: ISTART=0 ICHARG=2 # ISIF =3 NSW=0 IBRION=-1 # LWAVE=F # LCHARG=F ENCUT = 400.0 EDIFF = 1E-4 #EDIFFG=-1E-2 ISMEAR=0 LREAL=A SIGMA=0.1 #PREC=A POTIM=0.35 #NPAR=4 #NELMIN=4 NBANDS=64 第二步的INCAR： Start parameter for this run: ISTART=1 ICHARG=2 # ISIF =3 NSW=0 IBRION=-1 # LWAVE=F # LCHARG=F ENCUT = 400.0 EDIFF = 1E-4 #EDIFFG=-1E-2 ISMEAR=0 LREAL=A SIGMA=0.1 #PREC=A POTIM=0.35 #NPAR=4 #NELMIN=4 NBANDS=64 ALGO = Exact LOPTICS=.TRUE. （这两步vasp没有提示错误） 第三步INCAR： Start parameter for this run: ISTART=1 ICHARG=1 # ISIF =3 NSW=0 IBRION=-1 # LWAVE=F # LCHARG=F ENCUT = 400.0 EDIFF = 1E-4 #EDIFFG=-1E-2 ISMEAR=0 LREAL=A SIGMA=0.1 #PREC=A POTIM=0.35 #NPAR=4 #NELMIN=4 NBANDS=64 ALGO = GW0 LOPTICS=.TRUE. LSPECTRAL = T NOMEGA = 64 NBANDSGW=64 ENCUTGW=100 PRECFOCK=F LWANNIER90=.TRUE. LRPA=.FALSE. 在以上三步vasp计算过程中保持KPOINTS等其他的文件不变。使用的计算脚本如下： #!/bin/bash mpirun -np 16 vaspwannier9012 >out wait mkdir optic cp INCAR* KPOINTS POTCAR POSCAR WAVE* CHGCAR ./optic cd optic cp INCAR-optic INCAR mpirun -np 16 vaspwannier9012 >out wait mkdir gw cp INCAR-gw KPOINTS POTCAR POSCAR WAVE* CHGCAR ./gw cd gw cp INCAR-gw INCAR mpirun -np 16 vaspwannier9012 >out 问题出在第三步： out文件错误提示如下： NQ= 168 -0.0909 0.4545 0.4000, |.........|.........|.........|.........|.........|.........|.........|......... performed 923544000 updates of chi_q(r,r) total number of BLAS operations 102922.89 Gflops resolving degeneracies of screened two electron integrals calculate QP shifts <psi_nk| G(iteration)W_0 |psi_nk>: iteration 1 writing wavefunctions entering main loop optical routines imaginary and real dielectric function recalculating local potential from charge density CHECK_OCCUPATIONS: ERROR, system is not insulating. CHECK_OCCUPATIONS: ERROR, system is not insulating. CHECK_OCCUPATIONS: ERROR, system is not insulating. 根据提示，我不断增加NBANDS空带的条数，然而始终提示这个错误。 此外，如果把KPOINTS改为如下： KPOINTS 0 G 1 1 1 其他参数保持不变，结果就没问题。但是这时候的结果应该是没有意义的，应该是无法用wannier90算出能带的。恳请各位使用过GW方法的能够给个建议，谢谢大家啦！ |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
国自然面上复盘~欢迎讨论
已经有7人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
基金不中，共勉
已经有5人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有3人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
-
怎么查啊
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- pwscf收敛速度，参数调节与赝势
已经有19人回复
- CESTEP计算第五周期过渡金属氧化物参数设置的一些疑问
已经有11人回复
- 求助quantum-espresso使用问题！
已经有5人回复
- 新手求助一个quantum-espresso的使用问题！
已经有13人回复
- vasp接口wannier90计算问题
已经有24人回复
- 手册中计算能带时k点设置的问题
已经有12人回复
- 新手求助pwscf的进行磁性运算和LDA+U运算
已经有35人回复
- Problem - nscf - c_bands: eigenvalues not converged
已经有20人回复
- vasp中利用GW算法算能带结构
已经有9人回复
- ABINIT和VASP相比，计算GW，哪个更快啊？
已经有13人回复
- vasp跨节点运行出错，mpiexec_node-1 (handle_stdin_input 1089)
已经有5人回复
- ms中吸附能的计算
已经有40人回复
- 第一原理软件的两种磁性计算方法
已经有14人回复
- 【求助】结构优化不收敛
已经有9人回复
- 【求助】VASP计算光学性质时的GW方法（RPA, no local field effects)
已经有5人回复
- 【求助成功】Vasp 5.2.11编译出错
已经有10人回复
- 【求助】对于难收敛的体系
已经有9人回复
- 【求助】pwscf结构优化中结果出的怪异啊
已经有16人回复
- 【求助】VASP计算电子态密度
已经有34人回复
- 【求助】vasp5.2的GW怎么算？
已经有14人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
onion2440: 金币+50, ★★★★★最佳答案, 回复很详细，很有帮助，非常感谢！ 2015-07-13 10:39:16
liliangfang: 金币+1, 谢谢交流 2015-11-22 11:17:29
onion2440: 金币+50, ★★★★★最佳答案, 回复很详细，很有帮助，非常感谢！ 2015-07-13 10:39:16
liliangfang: 金币+1, 谢谢交流 2015-11-22 11:17:29
|
这取决于计算的初始态，如果HSE自洽结果不是金属，就可以算，否则和LDA/GGA并无区别。GW自洽也需要初始态，在这个问题上和G0W0没有区别。 ABINIT可以算，但结果也需要谨慎。 |
6楼2015-07-13 10:32:02
2楼2015-07-12 19:07:54
3楼2015-07-12 22:33:06
4楼2015-07-13 10:03:28
