24小时热门版块排行榜    

查看: 4342  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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
回复此楼
你想要什么,如何得到想要的,现在应该怎么做?
已阅   回复此楼   关注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的回帖
查看全部 8 个回答

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的回帖

黑暗游侠

铜虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:40 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:19 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 两块石头 +6 阿美_Lml888 2026-09-03 8/400 2026-09-05 13:52 by 阿美_Lml888
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见