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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » OUTCAR里面的total drift是指的什么，怎么减小呢？
81/1返回列表
查看: 4316  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xiaoxuan1111

铁杆木虫 (正式写手)

[求助] OUTCAR里面的total drift是指的什么，怎么减小呢？ 已有4人参与

用vasp+phonopy计算声子，三个声学支总是在G点处有大约0.5Thz的虚频无法消除。优化结构后得到的OUTCAR力的收敛情况如下：
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      0.98646      3.51321      0.00000         0.000006      0.000006     -0.000000
      4.06704      1.54029      0.00000        -0.000006     -0.000006      0.000000
      1.54029      0.98646      0.00000        -0.000006      0.000006      0.000000
      3.51321      4.06704      0.00000         0.000006     -0.000006     -0.000000
     -0.00000      0.00000      1.77621         0.000000      0.000000      0.000003
     -0.00000      0.00000      4.41310         0.000000      0.000000     -0.000003
      2.52675      2.52675      4.41310         0.000000      0.000000     -0.000003
      2.52675      2.52675      1.77621         0.000000      0.000000      0.000003
     -0.00000      2.52675      3.09466         0.000000      0.000000     -0.000000
      2.52675      0.00000      3.09466         0.000000      0.000000     -0.000000
      0.00000      0.00000      0.00000         0.000000      0.000000     -0.000000
      2.52675      2.52675      0.00000         0.000000      0.000000     -0.000000
-----------------------------------------------------------------------------------
    total drift:                               -0.000130     -0.000142     -0.000628
提高精度优化时，prec也设置为了high， 我感觉total-force应该是足够小了，而total drift一直是这三个值，基本不变。请问我的这个虚频是不是total drift太大导致的，这个值该怎么减小呢？或者说这个虚频该怎么消除呢？
这是我的INCAR
ISTART=0
ICHARG=2
ENCUT=500
IBRION=2
ISIF=3
EDIFF=1E-7
EDIFFG=-0.00001
POTIM=0.01
PREC=High
GGA=PS
NSW=1000
NPAR=2
KPAR=4
NELMIN=8
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
你想要什么，如何得到想要的，现在应该怎么做？
1楼 2016-05-06 16:13:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

armaity

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
xiaoxuan1111: 金币+5, ★有帮助, 提高k点试过了，虚频没什么变化。因为体系比较大，扩胞经常算不动，谢谢了。 2016-05-10 15:11:23
感觉你的情况声子虚频可能不是结构优化引起的，可以考虑计算声子时候的K点提高，或者修改超胞的大小，比如从222到333或者332，根据实际晶格常数调整。
赞一下 回复此楼
2楼2016-05-06 16:25:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saikun

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
xiaoxuan1111: 金币+5, ★有帮助, 谢谢，我先试试吧，毕竟增加k点和扩胞对我这个都不太适用。 2016-05-10 15:13:53
减小total drift, 你可增加参数NGX(YZ) and NGFX(YZ). 但我也觉得声子虚频不是有限total drift引起的。
赞一下 回复此楼
3楼2016-05-06 17:29:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黑暗游侠

铜虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
xiaoxuan1111: 金币+5, ★有帮助, 应该有可能 2016-05-10 15:15:16
声子虚频不是反应的结构不稳定，可能有相变吗？
赞一下 回复此楼
4楼2016-05-09 09:33:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoxuan1111

铁杆木虫 (正式写手)
引用回帖:
4楼: Originally posted by 黑暗游侠 at 2016-05-09 09:33:09
声子虚频不是反应的结构不稳定，可能有相变吗？

声学支0.5Thz的虚频会是结构不稳定产生的吗？
赞一下 回复此楼
你想要什么，如何得到想要的，现在应该怎么做？
5楼2016-05-10 15:16:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awdawdawdd

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
xiaoxuan1111: 金币+5, ★★★★★最佳答案, 非常感谢，要是有具体的操作方法就好了！ 2016-05-16 08:03:40
Gamma点的虚频表明你做unit cell优化时还没有优化完全，应该和supercell的选择关系不是太大

1.如果你用ISIF=3优化的话，有可能体积改变后引起basis set的completeness不够，推荐你优化之后把CONTCAR复制成POSCAR再优化几次，直到ionic step只有一步就结束了为止。

2.还有种可能就是gamma点的虚频是extrapolation的算法导致的，可以考虑施加 acoustic sum rules （asr）
参见http://atztogo.github.io/phonopy ... ?highlight=symmetry 中的symmetry部分

3.如果这些都没有问题有可能是lz的结构不太稳定了，其中一个因素就是symmetry不正确，可以考虑把对称性关掉（ISYM=0）进行relaxation，看看是否有显著地变化。或者lz可以把phonon eigen mode 画出来，然后把相关的原子往这个方向移动一下，然后再进行计算看是否消除imaginary frequency
赞一下(2人) 回复此楼

» 本帖已获得的红花（最新10朵）

xiaoxuan1111
6楼2016-05-16 03:38:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoxuan1111

铁杆木虫 (正式写手)
送红花一朵
引用回帖:
6楼: Originally posted by awdawdawdd at 2016-05-16 03:38:35
Gamma点的虚频表明你做unit cell优化时还没有优化完全，应该和supercell的选择关系不是太大

1.如果你用ISIF=3优化的话，有可能体积改变后引起basis set的completeness不够，推荐你优化之后把CONTCAR复制成POSCAR ...

非常感谢！！这个问题让我愁了好长时间了，我用了第2种方法，在band.conf文件中加了FC_SYMMETRY = 1，虚频变成了-0.0000002584Thz，相当于0了。
赞一下 回复此楼
你想要什么，如何得到想要的，现在应该怎么做？
7楼2016-05-16 08:15:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

big00fish

金虫 (小有名气)
FC_SYMMETRY = 1 管用

发自小木虫IOS客户端
赞一下 回复此楼
8楼2018-05-13 22:35:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaoxuan1111 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 XLGfUIdOmEAO 2026-08-27 3/150 2026-08-28 00:49 by 3RAFyYXYSwGS
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 XLGfUIdOmEAO 2026-08-27 4/200 2026-08-27 22:26 by 3RAFyYXYSwGS
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金不中，共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
关闭
请填处理意见
关闭