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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Abinit » VASP弛豫收敛太慢
191/1返回列表
查看: 5957  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

guoshukuan

铜虫 (小有名气)

[交流] VASP弛豫收敛太慢

VASP进行计算时第一步先粗算（EDIFF  =  1E-4; EDIFFG = -0.1; IBRION =  2）收敛后在精算（EDIFF  =  1E-5; EDIFFG = -0.01; IBRION =  1）。但第二步精算算了六十多步还没收敛，现在还在算。请高手指点怎样提高计算速度，谢谢。（INCAR附后）   

INCAR:

general:
System = Fe and vacancy in sigma3 GB
ISTART = 0 ; ICHARG=2
PREC=Accurate
ENCUT=320
ISMEAR = 0; SIGMA = 0.1
# IALGO  =  48
ALGO=Fast
GGA=91
ISPIN=2
VOSKOWN = 1
AMIX=0.2
BMIX=0.0001
AMIX_MAG=0.8
BMIX_MAG=0.0001
MAGMOM= 93*0 2*4 1*0
LREAL  =  Auto
# NPAR  =  1
# IMAGES=6
NELMIN = 5  
EDIFF  =  1E-5; EDIFFG = -0.01
NSW  =  60; IBRION =  1
MAXMIX=80
NFREE=10
ISIF  =  2
LWAVE  = .F.
LDAU = .TRUE.
LDAUL = -1  2  -1
LDAUU = 0  6.73  0
LDAUJ = 0  1.0  0
回复此楼

» 收录本帖的淘帖专辑推荐

vasp计算资料 第一性

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2012-02-16 11:32:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

lywiailyw

木虫 (小有名气)
★ ★
guoshukuan(金币+2): 谢谢参与
xueht987: 金币+1, 谢谢回帖，鼓励交流！ 2012-06-12 20:27:59
引用回帖:
10楼: Originally posted by haowenping at 2012-06-11 16:53:07
哦，您说的这个可能是的吧~但是做弛豫的时候加上磁矩这个就没什么必要了，做弛豫的时候是该加上这个ISPIN=2，但是MAGMOM就没必要了，弛豫本来就慢，有的甚至要跑上将近上百个离子步，这样下去每一步都慢，您觉得有 ...

还是有必要的，如果只打开ISPIN，系统默认的MAGMOM是每个原子上都有1μB的磁矩，这样很容易最终跑到铁磁态上面，如果有经验，最好还是指定每个原子的MAGMOM，否则很可能最终跑到一个自旋不合理的态。
另外，对于磁性体系不开ISPIN优化很可能会出现不收敛或者构型扭曲到严重未知状态的情况
赞一下(1人) 回复此楼
11楼2012-06-11 22:24:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)
★
guoshukuan(金币+1): 谢谢参与
引用回帖:
5楼: Originally posted by haowenping at 2012-02-16 19:25:02
就算是总能计算也不是在弛豫里面啊，是在静态里面。这个也不影响总能的！...

弛豫时加不加自旋极化，对晶格常数会有影响，怎么会不影响总能呢？
赞一下(1人) 回复此楼
14楼2012-06-12 20:31:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

haowenping

木虫 (著名写手)
★
guoshukuan(金币+2):谢谢参与
sunyang1988(金币+1): 谢谢交流 2012-02-17 09:05:05
算弛豫的时候可以把磁矩关了吧！
赞一下 回复此楼
2楼2012-02-16 12:14:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoboy185

金虫 (正式写手)
★
guoshukuan(金币+2):谢谢参与
sunyang1988(金币+1): 谢谢交流 2012-02-17 09:05:16
我觉得先升高EDIFF的值优化
然后再降低这个值进行进一步优化，

EDIFFG的值不用改
赞一下 回复此楼
3楼2012-02-16 12:22:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guoshukuan

铜虫 (小有名气)
引用回帖:
楼: Originally posted by haowenping at 2012-02-16 12:14:14:
算弛豫的时候可以把磁矩关了吧！

关了会不会影响总能？
赞一下 回复此楼
4楼2012-02-16 19:22:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)
★
guoshukuan(金币+2):谢谢参与
sunyang1988(金币+1): 谢谢交流 2012-02-17 09:05:38
引用回帖:
4楼: Originally posted by guoshukuan at 2012-02-16 19:22:16:
关了会不会影响总能？

就算是总能计算也不是在弛豫里面啊，是在静态里面。这个也不影响总能的！
赞一下 回复此楼
5楼2012-02-16 19:25:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guoshukuan

铜虫 (小有名气)
引用回帖:
楼: Originally posted by haowenping at 2012-02-16 19:25:02:
就算是总能计算也不是在弛豫里面啊，是在静态里面。这个也不影响总能的！

非常感谢你的回复：
1.我们计算总能就是直接读的弛豫收敛后的free energy。没有做静态计算。

2.手册上说如果允许了原胞形状的弛豫，就不能直接读弛豫后的能量，（Never take the energy obtained at the end of a relaxation run,if you allow for cell shape relaxations.Instead perform one additional static run at the end. ）.但我设ISIF=2。就不需要静态计算了啊。
3.如果MAGMOM= 93*0 2*4 1*0不影响总能，那它的作用在那里。
赞一下 回复此楼
6楼2012-02-17 08:29:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lywiailyw

木虫 (小有名气)
★ ★
guoshukuan(金币+2): 谢谢参与
xueht987: 金币+1, 谢谢回帖，鼓励交流！ 2012-06-12 20:27:03
引用回帖:
5楼: Originally posted by haowenping at 2012-02-16 19:25:02
就算是总能计算也不是在弛豫里面啊，是在静态里面。这个也不影响总能的！...

加磁性做弛豫当然对能量有影响，而且总能计算对于没有垂直电离能之类的情况，结构优化最后一步的能量就是需要的能量，不需要静态计算
赞一下 回复此楼
7楼2012-06-10 15:59:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)
★ ★
guoshukuan(金币+2): 谢谢参与
xueht987: 金币+1, 谢谢回帖，鼓励交流！ 2012-06-12 20:27:19
引用回帖:
7楼: Originally posted by lywiailyw at 2012-06-10 15:59:36
加磁性做弛豫当然对能量有影响，而且总能计算对于没有垂直电离能之类的情况，结构优化最后一步的能量就是需要的能量，不需要静态计算...

可能你说的最后一步优化得到的能量就是需要的能量，但是一般的计算都是为了想得到bands和dos，所以静态计算是必须的，个人认为。
赞一下 回复此楼
8楼2012-06-10 22:23:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lywiailyw

木虫 (小有名气)
★
guoshukuan(金币+2): 谢谢参与
引用回帖:
8楼: Originally posted by haowenping at 2012-06-10 22:23:35
可能你说的最后一步优化得到的能量就是需要的能量，但是一般的计算都是为了想得到bands和dos，所以静态计算是必须的，个人认为。...

算能带和DOS当然要做静态计算，那是另外一码事情了呃。我的意思就是，单考虑能量大小的话，弛豫最后一个离子步的最后一个电子步得到的能量与静态计算得到的能量几乎完全一样；另外，弛豫过程加不加磁性对磁性体系的能量有影响，一般在几十到几百个meV之间，有时候甚至更大。
赞一下 回复此楼
9楼2012-06-11 16:40:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)
★
xueht987: 金币+1, 谢谢回帖，鼓励交流！ 2012-06-12 20:27:43
引用回帖:
9楼: Originally posted by lywiailyw at 2012-06-11 16:40:17
算能带和DOS当然要做静态计算，那是另外一码事情了呃。我的意思就是，单考虑能量大小的话，弛豫最后一个离子步的最后一个电子步得到的能量与静态计算得到的能量几乎完全一样；另外，弛豫过程加不加磁性对磁性体系的 ...

哦，您说的这个可能是的吧~但是做弛豫的时候加上磁矩这个就没什么必要了，做弛豫的时候是该加上这个ISPIN=2，但是MAGMOM就没必要了，弛豫本来就慢，有的甚至要跑上将近上百个离子步，这样下去每一步都慢，您觉得有必要吗？在弛豫的时候加上ISPIN=2，不加MAGMOM对体系能量有何关系呢？我并没有否认在弛豫的时候打开自旋，我只是要表达的意思是没必要看磁矩，个人认为！
赞一下 回复此楼
10楼2012-06-11 16:53:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)
★
xueht987: 金币+1, 谢谢回帖，鼓励交流！ 2012-06-12 20:28:20
引用回帖:
11楼: Originally posted by lywiailyw at 2012-06-11 22:24:46
还是有必要的，如果只打开ISPIN，系统默认的MAGMOM是每个原子上都有1μB的磁矩，这样很容易最终跑到铁磁态上面，如果有经验，最好还是指定每个原子的MAGMOM，否则很可能最终跑到一个自旋不合理的态。
另外，对于磁 ...

这一点我倒真是学习了~我真不知道MAGMOM对弛豫的时候的最终态还会有影响！我再试试看看用加MAGMOM和不加MAGMOM分别得到的结果后看看能量是否变化会很大。
赞一下 回复此楼
12楼2012-06-12 09:42:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lywiailyw

木虫 (小有名气)
★
xueht987: 金币+1, 谢谢回帖，鼓励交流！ 2012-06-12 20:28:32
引用回帖:
12楼: Originally posted by haowenping at 2012-06-12 09:42:20
这一点我倒真是学习了~我真不知道MAGMOM对弛豫的时候的最终态还会有影响！我再试试看看用加MAGMOM和不加MAGMOM分别得到的结果后看看能量是否变化会很大。...

你可以试试，一般是几十到几百个meV
赞一下 回复此楼
13楼2012-06-12 13:58:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lywiailyw

木虫 (小有名气)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
WDD880227: 金币+1, 多谢您的指导 2012-06-14 08:40:41
引用回帖:
14楼: Originally posted by fanarsenal at 2012-06-12 20:31:33
弛豫时加不加自旋极化，对晶格常数会有影响，怎么会不影响总能呢？...

一般VASP弛豫都选择ISIF=2，也就是固定晶格常数了，此时加不加自旋对晶格常数没有影响。ISIF=3不是随便可以用的，有技巧。
对于强磁交换体系，加自旋极化会影响离子的平衡位置，从而影响总能
赞一下 回复此楼
15楼2012-06-13 19:20:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名


★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖仅楼主可见
赞一下
16楼2012-12-16 06:47:49
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

匿名


★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖仅楼主可见
赞一下
17楼2012-12-16 06:50:33
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

lywiailyw

木虫 (小有名气)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fzx2008: 金币+1, 谢谢指教！ 2012-12-16 23:53:11
引用回帖:
17楼: Originally posted by elementAt at 2012-12-16 06:50:33
这样啊，，，那如果在你“确定晶格常数的计算中”， ISPIN是打开的吗？...

不管ISIF=2还是=3，我一般都是先不开ISPIN，等粗优化完成之后再开ISPIN算的
赞一下 回复此楼
18楼2012-12-16 21:40:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

暖暖67

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
弛豫时是收敛的 但是静态算一下结果却不收敛了 怎么
赞一下 回复此楼
19楼2013-06-20 22:19:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 guoshukuan 的主题更新
191/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:59 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:50 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-15 15:27 by CgyNCDVNhGVg
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 15:02 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:50 by VMDqgnlKimtZ
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 14:38 by VMDqgnlKimtZ
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 14:25 by VMDqgnlKimtZ
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:14 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 05:22 by 5BDX0d0WFp7t
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭