| 查看: 4903 | 回复: 12 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
burningdays铜虫 (小有名气)
|
[求助]
求助VASP计算频率为何总是有虚频？如何消掉？
|
|
小弟初学vasp，想用vasp计算频率，刚开始直接用IBRION=5直接计算，发现有虚频，于是升高精度还是发现有虚频，这让我很郁闷，于是我又装了Phonopy，用DPFT的方法算（PS：这才知道为啥大家都用VASP+PHONOPY+DFPT计算频率，这比直接算快多了。。 ），得到同样的结果，于是我又扩胞（本来是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个虚频，而且数值都很接近，我这前前后后都折腾了快一月了，结果还是这样的。。 ，求各位前辈帮小弟看看我这问题出在哪里，万分感谢！！ |
回复此楼» 猜你喜欢
自荐读博
已经有5人回复
-
求个博导看看
已经有16人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有5人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
需要合成515-64-0，50g，能接单的留言
已经有4人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
带资进组求博导收留
已经有10人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
11楼2018-10-09 19:27:02
2楼2016-12-30 13:47:50
burningdays
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 30.1
- 散金: 10
- 红花: 2
- 帖子: 142
- 在线: 32.6小时
- 虫号: 3437080
- 注册: 2014-09-23
- 专业: 材料物理化学
3楼2016-12-30 14:51:51
burningdays
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 30.1
- 散金: 10
- 红花: 2
- 帖子: 142
- 在线: 32.6小时
- 虫号: 3437080
- 注册: 2014-09-23
- 专业: 材料物理化学
4楼2016-12-30 15:02:34