| 查看: 2391 | 回复: 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方法的能够给个建议，谢谢大家啦！ |
回复此楼» 猜你喜欢
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
申请26博士
已经有5人回复
-
职称评审没过，求安慰
已经有22人回复
-
垃圾破二本职称评审标准
已经有15人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
3楼2015-07-12 22:33:06
2楼2015-07-12 19:07:54
4楼2015-07-13 10:03:28
5楼2015-07-13 10:20:32