|
首页 >> 第一原理 >> 查看话题

求助VASP计算频率为何总是有虚频?如何消掉?

burningdays 举报此信息
小弟初学vasp,想用vasp计算频率,刚开始直接用IBRION=5直接计算,发现有虚频,于是升高精度还是发现有虚频,这让我很郁闷,于是我又装了Phonopy,用DPFT的方法算(PS:这才知道为啥大家都用VASP+PHONOPY+DFPT计算频率,这比直接算快多了。。:cry:),得到同样的结果,于是我又扩胞(本来是3x3x2,扩大到3x3x3),还是有虚频,我都快崩溃了,始终是找不出为啥我的计算总有那么几个虚频,请各位前辈帮忙分析一下,小弟跪谢!!
      
        以下是我的incar,我是先用IBRION=2优化,在改到1优化达到收敛,在用5计算频率。
       第一步优化INCAR:
SYSTEM=BRUCITE3X3X2_OPT1
ENCUT=600
ISTART=0;ICHARGE=2
ISMEAR=0;SIGMA=0.05
PREC=Accurate
LREAL=.FALSE.
IALGO=48
EDIFF=1E-5
EDIFFG=-1E-2
ISIF=3
IBRION=2
LWAVE=.FALSE.
NSW=400
LCHARG=.FALSE.
NPAR=4
NELMIN=6
第二部优化INCAR:
SYSTEM=BRUCITE3X3X2_OPT2
ENCUT=600
ISTART=0;ICHARGE=2
ISMEAR=0;SIGMA=0.05
PREC=Accurate
LREAL=.FALSE.
IALGO=48
EDIFF=1E-7
EDIFFG=-1E-3
ISIF=3
IBRION=1;POTIM=0.2
LWAVE=.FALSE.
NSW=400
LCHARG=.FALSE.
NPAR=4
NELMIN=8
频率计算我是直接在第二步优化基础上改IBRION=5,POTIM=0.015加NFREE=2来做的。
      K点我使用Monkhorst-Pack方法,3x3x3(经过测试的),赝势使用的是PBE标准赝势,晶胞为3x3x2的超胞,计算的晶体为brucite(Mg(OH)2)。
      vasp+phonopy+dfpt的计算INCAR几乎与上面的类似,结果还是有虚频,只是优化的是单胞,扩胞过后用IBRION=8实现的。
    IBRION=5虚频:
266 f/i=    0.005270 THz     0.033113 2PiTHz    0.175794 cm-1     0.021796 meV
267 f/i=    0.007348 THz     0.046168 2PiTHz    0.245099 cm-1     0.030388 meV
268 f/i=    0.070140 THz     0.440705 2PiTHz    2.339633 cm-1     0.290078 meV
269 f/i=    2.263820 THz    14.224001 2PiTHz   75.512910 cm-1     9.362411 meV
270 f/i=    2.267130 THz    14.244800 2PiTHz   75.623325 cm-1     9.376101 meV
提高精度EDIFF=1E-08,EDIFFG=-1E-04之后的虚频:
266 f/i=    0.002587 THz     0.016257 2PiTHz    0.086303 cm-1     0.010700 meV
267 f/i=    0.003011 THz     0.018917 2PiTHz    0.100426 cm-1     0.012451 meV
268 f/i=    0.074536 THz     0.468321 2PiTHz    2.486239 cm-1     0.308255 meV
269 f/i=    2.272819 THz    14.280542 2PiTHz   75.813072 cm-1     9.399626 meV
270 f/i=    2.273241 THz    14.283197 2PiTHz   75.827166 cm-1     9.401374 meV
DPFT虚频:
266 f/i=    0.000825 THz     0.005186 2PiTHz    0.027530 cm-1     0.003413 meV
267 f/i=    0.003373 THz     0.021191 2PiTHz    0.112502 cm-1     0.013948 meV
268 f/i=    0.003373 THz     0.021191 2PiTHz    0.112502 cm-1     0.013948 meV
269 f/i=    2.076583 THz    13.047555 2PiTHz   69.267349 cm-1     8.588060 meV
270 f/i=    2.076583 THz    13.047555 2PiTHz   69.267349 cm-1     8.588060 meV
扩胞到3x3x3之后的虚频:
401 f/i=    0.001677 THz     0.010538 2PiTHz    0.055946 cm-1     0.006936 meV
402 f/i=    0.003328 THz     0.020912 2PiTHz    0.111016 cm-1     0.013764 meV
403 f/i=    0.003328 THz     0.020912 2PiTHz    0.111016 cm-1     0.013764 meV
404 f/i=    2.074640 THz    13.035346 2PiTHz   69.202531 cm-1     8.580023 meV
405 f/i=    2.074640 THz    13.035346 2PiTHz   69.202531 cm-1     8.580023 meV
     总是会有5个虚频,而且数值都很接近,我这前前后后都折腾了快一月了,结果还是这样的。。:cry:,求各位前辈帮小弟看看我这问题出在哪里,万分感谢!!:hand:


展开阅读全文

查看全文请下载APP

用户评论