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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » OUTCAR里面的total drift是指的什么，怎么减小呢？
81/1返回列表
查看: 3897  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭