| 查看: 4443 | 回复: 4 | |||
chord7806新虫 (初入文坛)
|
[求助]
<求助>VASP 结构优化收敛，静态计算不收敛。（INCAR KPOINTS POSCAR附上） 已有1人参与
|
|
最近计算了一个体系，共40原子。两种磁性离子。 先进行了结构优化。得到了CONTCAR，然后cp CONTCAR POSCAR静态计算能量始终不收敛。 希望各位老师，同学，达人，帮帮忙，看一下，到底问题出在哪？不胜感激。 下面是初始优化时的INCAR KPOINTS POSCAR 部分POSCAR： system=****** 1.000000000000000 5.163711736218603 0.000000000000000 0.030191807469367 0.000000000000000 5.546700000000000 0.000000000000000 0.000000000000000 0.000000000000000 14.830600000000000 24 4 4 8 Direct 0.3840999901294709 0.9585000276565552 0.1205499991774557 0.1970999985933304 0.1957000046968459 0.4712499994784595 0.3228000104427338 0.6952999830245971 0.4703499991446734 0.1159000098705291 0.4585000276565551 0.1294500008225444 0.3029000014066697 0.6957000046968459 0.2787500005215405 0.1771999895572662 0.1952999830245970 0.2796500008553267 0.6159000098705292 0.0414999723434448 0.3794500008225443 0.8029000014066697 0.8042999953031542 0.0287500005215406 0.6771999895572662 0.3047000169754029 0.0296500008553267 0.8840999901294708 0.5414999723434448 0.3705499991774557 INCAR: System = ****** PREC = Accurate ISTART = 0 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 500 EDIFF = 1e-6 NELM = 40 GGA = PS NSW = 300 IBRION = 2 ISIF = 3 EDIFFG = -0.001 MAGMOM = 24*0 -3 3 3 -3 -3 3 3 -3 8*0 LDAU = TRUE LDAUU = 0.00 3.00 3.00 0.00 LDAUJ = 0.00 0.00 0.00 0.00 LDAUTYPE = 2 LDAUL = -1 2 2 -1 POTIM = 0.2 LMAXMIX = 4 ISMEAR = 0 SIGMA = 0.01 LORBIT = 11 LCHARG = FALSE LWAVE = FALSE KPOINTS： Automatic 0 Monkhorst-Pack 6 6 2 0 0 0 计算完以后有这样的提示： 73 F= -.30352737E+03 E0= -.30352737E+03 d E =-.220886E-05 mag= 0.0000 curvature: -0.06 expect dE=-0.192E-06 dE for cont linesearch -0.395E-08 ZBRENT: bisectioning opt : 0.4334 next Energy= -303.527366 (dE=-0.221E-05) bond charge predicted N E dE d eps ncg rms rms(c) DAV: 1 -0.303527365963E+03 0.40205E-06 -0.12261E-07 2024 0.516E-03 0.155E-03 DAV: 2 -0.303527366101E+03 -0.13772E-06 -0.10919E-07 2024 0.500E-02 0.155E-03 DAV: 3 -0.303527366167E+03 -0.66240E-07 -0.41386E-08 2024 0.952E-03 0.155E-03 DAV: 4 -0.303527366279E+03 -0.11196E-06 -0.75204E-08 2024 0.228E-02 0.156E-03 DAV: 5 -0.303527366865E+03 -0.58592E-06 -0.15596E-06 2024 0.303E-02 74 F= -.30352737E+03 E0= -.30352737E+03 d E =-.315061E-05 mag= 0.0000 curvature: -0.02 expect dE=-0.165E-07 dE for cont linesearch -0.300E-08 trial: gam= 0.00000 g(F)= 0.285E-06 g(S)= 0.390E-06 ort = 0.815E-06 (trialstep = 0.450E+00) search vector abs. value= 0.675E-06 reached required accuracy - stopping structural energy minimisation OUTCAR里的精度： Total 0.01803 0.01087 0.01049 0.00000 0.00000 0.01127 in kB 0.06868 0.04140 0.03997 0.00000 0.00000 0.04294 external pressure = 0.05 kB Pullay stress = 0.00 kB 和 POSITION TOTAL-FORCE (eV/Angst) ----------------------------------------------------------------------------------- 1.95504 5.28764 1.81560 0.000493 -0.000529 -0.000062 1.04003 1.02590 7.01562 -0.000502 -0.000445 0.000609 1.68657 3.89042 6.97157 -0.000104 0.000333 -0.000392 0.62891 2.52622 1.91875 -0.000493 -0.000529 0.000062 1.57671 3.78732 4.13456 0.000502 -0.000445 -0.000609 0.93018 1.12900 4.17862 0.000104 0.000333 0.000392 3.19679 0.22461 5.64936 -0.000192 0.000558 0.000228 4.13049 4.50077 0.45157 0.000433 -0.000113 -0.000710 3.48222 1.63383 0.49540 0.000351 0.000060 0.000336 4.55507 2.98603 5.55370 0.000192 0.000558 -0.000228 然后用优化的CONTCAR复制成POSCAR做静态计算，问题出现了，能量不收敛： 下面是提示内容： DAV: 57 -0.295043281006E+03 -0.41838E-03 -0.36422E-01 7424 0.295E+01 0.338E+00 DAV: 58 -0.294995226835E+03 0.48054E-01 -0.55161E-01 6560 0.474E+01 0.334E+00 DAV: 59 -0.294998850811E+03 -0.36240E-02 -0.69107E-01 7696 0.140E+01 0.355E+00 DAV: 60 -0.295061875171E+03 -0.63024E-01 -0.12717E-01 6760 0.109E+01 1 F= -.29506188E+03 E0= -.29506186E+03 d E =-.325714E-04 mag= 0.0000 而且发现OUTCAR中的 external pressure变的很大：external pressure = -3550.84 kB，与此同时，磁矩也变的不准确了。 我的问题是： 1，有人遇到过这样的情况吗：优化结构后，算单点能量不收敛的。 2，我的INCAR中有什么是不对的吗，如果需要我贴出什么其他的信息，请在本帖留言。 真心希望问题得到解决，感激不尽。。。。。 [ Last edited by chord7806 on 2012-5-29 at 16:13 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性 |
» 猜你喜欢
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有13人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有7人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有11人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有9人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 结构优化计算不收敛，是怎么回事？求帮助
已经有6人回复
- 关于VASP计算收敛的判断
已经有25人回复
- 关于VASP静态计算
已经有4人回复
- VASP静态自洽计算收敛特慢！
已经有31人回复
- VASP NEB 不收敛问题
已经有25人回复
- VASP计算表面吸附 但是能量总是不收敛/。。。
已经有9人回复
- 用VASP计算NEB都十几天了还不收敛，急啊！！！
已经有8人回复
- vasp中，做截断能测试时，测试收敛的标准
已经有10人回复
- 用VASP优化结构收敛不了，求指教
已经有10人回复
- 做结构优化和做静态计算能量不一致
已经有7人回复
- 用Dmol优化收敛的结构，用vasp不收敛？
已经有9人回复
- vasp4.6收敛的结构在vasp5.2 计算不收敛
已经有7人回复
- 关于vasp收敛和优化
已经有24人回复
- vasp优化难收敛，请大家给我看看
已经有5人回复
- VASP弛豫收敛太慢
已经有18人回复
- vasp优化力收敛条件
已经有6人回复
- 求助：vasp结构优化不能收敛
已经有18人回复
- vasp计算：优化能收敛，但同样的设置静态计算不收敛
已经有3人回复
- VASP优化结构遇到的问题。
已经有9人回复
- 【求助】VASP收敛问题
已经有15人回复
- 【求助】结构优化不收敛
已经有9人回复
- 【求助】vasp截断能不收敛
已经有7人回复
- 【求助】VASP结构优化
已经有7人回复
- 【求助】VASP结构优化不收敛怎么办？
已经有5人回复
wanghouyuan
铜虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 224.1
- 红花: 1
- 帖子: 75
- 在线: 15.2小时
- 虫号: 1144660
- 注册: 2010-11-12
- 性别: GG
- 专业: 理论和计算化学
2楼2012-05-29 16:50:33
chord7806
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 108.8
- 帖子: 31
- 在线: 107.3小时
- 虫号: 724541
- 注册: 2009-03-17
- 专业: 凝聚态物性 II ：电子结构
|
这个是我静态计算的INCAR，你看一下吧： System = ******* PREC = Accurate ISTART = 0 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 500 EDIFF = 1e-6 NELM = 60 GGA = PS NSW = 0 IBRION = -1 ISIF = 3 EDIFFG = -0.001 MAGMOM = 24*0 -3 3 3 -3 -3 3 3 -3 8*0 LDAU = TRUE LDAUU = 0.00 3.00 3.00 0.00 LDAUJ = 0.00 0.00 0.00 0.00 LDAUTYPE = 2 LDAUL = -1 2 2 -1 POTIM = 0.2 LMAXMIX = 4 ISMEAR = 0 SIGMA = 0.01 LORBIT = 11 LCHARG = FALSE LWAVE = FALSE 我就把IBRION由原来的2改成-1了，NSW有300改成0了。。。 多谢关注。。 |
3楼2012-05-29 17:07:23
wanghouyuan
铜虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 224.1
- 红花: 1
- 帖子: 75
- 在线: 15.2小时
- 虫号: 1144660
- 注册: 2010-11-12
- 性别: GG
- 专业: 理论和计算化学
4楼2012-05-29 17:29:21
5楼2018-07-07 10:22:39
