应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求助VASP计算频率为何总是有虚频？如何消掉？
131/2返回列表12下一页
查看: 4871  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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个虚频，而且数值都很接近，我这前前后后都折腾了快一月了，结果还是这样的。。，求各位前辈帮小弟看看我这问题出在哪里，万分感谢！！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2016-12-30 10:30:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

匿名

用户注销 (职业作家)
★ ★
漫天飘雪: 金币+2, 谢谢交流 2017-01-02 10:52:08
本帖仅楼主可见
赞一下
2楼2016-12-30 13:47:50
 已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

FourierLiu

新虫 (初入文坛)
水镁石我也算过的确有虚频，我怀疑是不是开始的结构不对

发自小木虫IOS客户端
赞一下(1人) 回复此楼
7楼2017-01-03 23:51:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

burningdays

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by zyl雨田 at 2016-12-30 13:47:50
小木虫上这么多关于消虚频的帖子，这么多的方法总结。你不找找？

谢谢你的回复！
我在网上查过，之前好几个关于虚频的帖子都是你提的，貌似你最后的解决办法是扩胞吧？但是我这里也扩了胞，结果虚频还是在，可能是我扩胞的方向不对？还是我的优化的设置有问题？
赞一下 回复此楼
3楼2016-12-30 14:51:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

burningdays

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by zyl雨田 at 2016-12-30 13:47:50
小木虫上这么多关于消虚频的帖子，这么多的方法总结。你不找找？

http://muchong.com/html/201406/7591851.html
对，还有一个就是你的这个帖子里面，说可能是对称性的问题，INCAR中加ISYM=0，这个我也试了，也不管用，不过这个我是用单分子水体系试的。我也算了单分子水的频率，也是死活有3个虚频，所以我就怀疑是不是我的优化设置有问题
赞一下 回复此楼
4楼2016-12-30 15:02:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qwangar

金虫 (小有名气)
试试提高优化精度 貌似力精度不够

发自小木虫IOS客户端
赞一下 回复此楼
I know nothing.
5楼2016-12-30 17:10:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

burningdays

铜虫 (小有名气)
引用回帖:
5楼: Originally posted by qwangar at 2016-12-30 17:10:12
试试提高优化精度 貌似力精度不够

我后来把力收敛标准提高到了10-4了  还是不行
赞一下 回复此楼
6楼2016-12-30 19:22:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CarlyQ_Q

金虫 (小有名气)
到底是什么原因，解决了吗？
赞一下 回复此楼
8楼2017-01-11 23:37:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

burningdays

铜虫 (小有名气)
引用回帖:
7楼: Originally posted by FourierLiu at 2017-01-03 23:51:03
水镁石我也算过的确有虚频，我怀疑是不是开始的结构不对

不是结构的问题，是扩胞的问题
赞一下 回复此楼
9楼2017-02-06 09:16:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

burningdays

铜虫 (小有名气)
引用回帖:
8楼: Originally posted by CarlyQ_Q at 2017-01-11 23:37:20
到底是什么原因，解决了吗？

是我的超胞扩的不对，我用3x3x2的胞算出来怎么都有虚频，但是我后来改成2x2x1就没有了，这很奇怪，同样的再把胞扩到4x4x2也是没有虚频，我也不清楚为什么。。
另外，消虚频可以参考这个帖子 dxcharlary前辈的意见http://muchong.com/html/201406/7591851.html
赞一下 回复此楼
10楼2017-02-06 09:24:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 burningdays 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭