 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4458 | 回復: 4 | |||
chord7806新蟲 (初入文壇)
|
[求助]
<求助>VASP 結(jié)構(gòu)優(yōu)化收斂，靜態(tài)計算不收斂。（INCAR KPOINTS POSCAR附上） 已有1人參與
|
|
最近計算了一個體系，共40原子。兩種磁性離子。 先進行了結(jié)構(gòu)優(yōu)化。得到了CONTCAR，然后cp CONTCAR POSCAR靜態(tài)計算能量始終不收斂。 希望各位老師，同學，達人，幫幫忙，看一下，到底問題出在哪？不勝感激。 下面是初始優(yōu)化時的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 然后用優(yōu)化的CONTCAR復制成POSCAR做靜態(tài)計算，問題出現(xiàn)了，能量不收斂： 下面是提示內(nèi)容： 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 而且發(fā)現(xiàn)OUTCAR中的 external pressure變的很大：external pressure = -3550.84 kB，與此同時，磁矩也變的不準確了。 我的問題是： 1，有人遇到過這樣的情況嗎：優(yōu)化結(jié)構(gòu)后，算單點能量不收斂的。 2，我的INCAR中有什么是不對的嗎，如果需要我貼出什么其他的信息，請在本帖留言。 真心希望問題得到解決，感激不盡。。。。。 [ Last edited by chord7806 on 2012-5-29 at 16:13 ] |
回復此樓» 收錄本帖的淘帖專輯推薦
 第一性 |
» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有9人回復
-
293求調(diào)劑
已經(jīng)有14人回復
-
0703化學336分求調(diào)劑
已經(jīng)有6人回復
-
304求調(diào)劑
已經(jīng)有9人回復
-
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名
已經(jīng)有3人回復
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復
-
301求調(diào)劑
已經(jīng)有6人回復
-
一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復
-
302求調(diào)劑
已經(jīng)有6人回復
-
26博士申請
已經(jīng)有3人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 結(jié)構(gòu)優(yōu)化計算不收斂，是怎么回事？求幫助
已經(jīng)有6人回復
- 關(guān)于VASP計算收斂的判斷
已經(jīng)有25人回復
- 關(guān)于VASP靜態(tài)計算
已經(jīng)有4人回復
- VASP靜態(tài)自洽計算收斂特慢！
已經(jīng)有31人回復
- VASP NEB 不收斂問題
已經(jīng)有25人回復
- VASP計算表面吸附 但是能量總是不收斂/。。。
已經(jīng)有9人回復
- 用VASP計算NEB都十幾天了還不收斂，急啊�。�！
已經(jīng)有8人回復
- vasp中，做截斷能測試時，測試收斂的標準
已經(jīng)有10人回復
- 用VASP優(yōu)化結(jié)構(gòu)收斂不了，求指教
已經(jīng)有10人回復
- 做結(jié)構(gòu)優(yōu)化和做靜態(tài)計算能量不一致
已經(jīng)有7人回復
- 用Dmol優(yōu)化收斂的結(jié)構(gòu)，用vasp不收斂？
已經(jīng)有9人回復
- vasp4.6收斂的結(jié)構(gòu)在vasp5.2 計算不收斂
已經(jīng)有7人回復
- 關(guān)于vasp收斂和優(yōu)化
已經(jīng)有24人回復
- vasp優(yōu)化難收斂，請大家給我看看
已經(jīng)有5人回復
- VASP弛豫收斂太慢
已經(jīng)有18人回復
- vasp優(yōu)化力收斂條件
已經(jīng)有6人回復
- 求助：vasp結(jié)構(gòu)優(yōu)化不能收斂
已經(jīng)有18人回復
- vasp計算：優(yōu)化能收斂，但同樣的設(shè)置靜態(tài)計算不收斂
已經(jīng)有3人回復
- VASP優(yōu)化結(jié)構(gòu)遇到的問題。
已經(jīng)有9人回復
- 【求助】VASP收斂問題
已經(jīng)有15人回復
- 【求助】結(jié)構(gòu)優(yōu)化不收斂
已經(jīng)有9人回復
- 【求助】vasp截斷能不收斂
已經(jīng)有7人回復
- 【求助】VASP結(jié)構(gòu)優(yōu)化
已經(jīng)有7人回復
- 【求助】VASP結(jié)構(gòu)優(yōu)化不收斂怎么辦？
已經(jīng)有5人回復
wanghouyuan
銅蟲 (小有名氣)
- 應助: 6 (幼兒園)
- 金幣: 224.1
- 紅花: 1
- 帖子: 75
- 在線: 15.2小時
- 蟲號: 1144660
- 注冊: 2010-11-12
- 性別: GG
- 專業(yè): 理論和計算化學
2樓2012-05-29 16:50:33
chord7806
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 108.8
- 帖子: 31
- 在線: 107.3小時
- 蟲號: 724541
- 注冊: 2009-03-17
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
這個是我靜態(tài)計算的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了。。。 多謝關(guān)注。。 |
3樓2012-05-29 17:07:23
wanghouyuan
銅蟲 (小有名氣)
- 應助: 6 (幼兒園)
- 金幣: 224.1
- 紅花: 1
- 帖子: 75
- 在線: 15.2小時
- 蟲號: 1144660
- 注冊: 2010-11-12
- 性別: GG
- 專業(yè): 理論和計算化學
4樓2012-05-29 17:29:21
5樓2018-07-07 10:22:39
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 chord7806 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑（0805） +3 | 煙汐憶海 2026-03-16 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 293求調(diào)劑 +9 | zjl的號 2026-03-16 | 14/700 |
|
|
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +4 | share_joy 2026-03-16 | 4/200 |
|
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學666 2026-03-16 | 6/300 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +7 | 邱gl 2026-03-11 | 8/400 |
|
|
|
[考研] 東南大學364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 0703化學調(diào)劑 +6 | 妮妮ninicgb 2026-03-15 | 9/450 |
|
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
|
[考研] 331求調(diào)劑（0703有機化學 +5 | ZY-05 2026-03-13 | 6/300 |
|
|
|
[考研] 一志愿中科院，化學方向，295求調(diào)劑 +4 | 一氧二氮 2026-03-11 | 4/200 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[考研] 工科調(diào)劑 +4 | Jiang191123！ 2026-03-11 | 4/200 |
|
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
|
[考博] 福州大學楊黃浩課題組招收2026年專業(yè)學位博士研究生，2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|