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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » OUTCAR里面的total drift是指的什么，怎么减小呢？
81/1返回列表
查看: 4139  |  回复: 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一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:37 by ue3ir18jc3
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k37jurhrau 2026-05-16 4/200 2026-05-17 01:35 by ue3ir18jc3
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:57 by ue3ir18jc3
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 x0mp7owy2b 2026-05-15 5/250 2026-05-17 00:29 by xiangfeng
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 7/350 2026-05-16 23:17 by ue3ir18jc3
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 4/200 2026-05-16 23:10 by ue3ir18jc3
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 5/250 2026-05-16 23:05 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 xx7gd5zq4e 2026-05-15 6/300 2026-05-16 22:25 by lindoudou
[基金申请] 重磅！青年科学基金项目（C类）资助增幅预计超过50% +6 水和泥不是水泥 2026-05-13 8/400 2026-05-16 19:50 by Equinoxhua
[有机交流] 求有机合成大神指点三硫酸乙烯酯（CAS：2793408-99-6）的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[公派出国] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 5/250 2026-05-16 16:35 by x28q7dxf75
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 6/300 2026-05-16 16:17 by 0i5p09z61n
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[硕博家园] 考博自荐 +3 科研狗111 2026-05-13 4/200 2026-05-16 11:45 by 科研狗111
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 cjf4bx70cj 2026-05-14 6/300 2026-05-16 11:16 by h3oerqvkv9
[文学芳草园] 裁员滚滚，退居二线 +4 J_wei 2026-05-10 4/200 2026-05-16 10:52 by zh10246
[基金申请] 这年头没有找到涵评专家，还有中面上的可能吗 +11 dd921ww 2026-05-12 13/650 2026-05-16 09:16 by Howard28
[考博] 西南大学考核制博士 +4 lijunjie84 2026-05-11 7/350 2026-05-15 23:20 by 同仁堂教主
[基金申请] 青C资助名额大幅增加！ +11 西葫芦炒鸡蛋 2026-05-13 15/750 2026-05-15 14:36 by ambravo
信息提示
关闭
请填处理意见
关闭