| 查看: 5641 | 回复: 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人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助 用VASP计算频率时有虚频，如何确定虚频是哪个原子或哪几个原子么？
已经有21人回复
- 求助：如何使用VASP计算结构中某个原子或者离子震动频率
已经有2人回复
- vasp计算团簇振动频率 再次求助
已经有0人回复
- VASP计算光学性质，考虑local-field effects(LFE)
已经有4人回复
- 求助声子计算
已经有1人回复
- 金属表面反应过渡态频率计算
已经有1人回复
- 用VASP-DFPT+phonopy计算石墨烯声子谱求助
已经有37人回复
- 求助，VASP如何计算分析空位型扩散的空位振动熵和有效频率
已经有6人回复
- vasp+phonopy计算一个新结构的声子谱有虚频，怎么消掉！！急！！！
已经有4人回复
- 【求助】VASP做频率计算算出的虚频都上2000+cm-1了，请教是不方法错误。
已经有1人回复
- VASP频率计算结果分析细节(求助完毕，分享结果)
已经有11人回复
- 【求助】计算振动
已经有6人回复
- 【求助】请教vasp频率计算中的NSW
已经有7人回复
- 【求助】请教vasp如何做寻找过渡态的频率验证
已经有1人回复
- 【求助】重金悬赏 计算软件的几个问题，悬赏150个金币。
已经有11人回复
- 【求助】请问pyVib2 软件可以显示vasp所计算的频率吗
已经有0人回复
- 【求助】vasp 5.2 光学计算的一些细节
已经有8人回复
- 【求助】VASP有限差分计算力常数和声子振动频率
已经有1人回复
- 【求助完毕】Vasp中Gamma点的振动能频率是验证过度态的充要条件吗?
已经有6人回复
burningdays
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 30.1
- 散金: 10
- 红花: 2
- 帖子: 142
- 在线: 32.6小时
- 虫号: 3437080
- 注册: 2014-09-23
- 专业: 材料物理化学
9楼2017-02-06 09:16:34
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
10