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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » OUTCAR里面的total drift是指的什么，怎么减小呢？
81/1返回列表
查看: 4006  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 294求调剂材料与化工专硕 +10 陌の森林 2026-03-18 10/500 2026-03-19 12:35 by peike
[考研] 287求调剂 +3 晨昏线与星海 2026-03-19 4/200 2026-03-19 12:32 by peike
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 材料专硕英一数二306 +5 z1z2z3879 2026-03-18 5/250 2026-03-19 07:43 by BruceLiu320
[考研] 328求调剂，英语六级551，有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 一志愿西南交大，求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 331求调剂（0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 一志愿华中师范071000，325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
信息提示
关闭
请填处理意见
关闭