 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4459 | 回復(fù): 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)計算能量始終不收斂。 希望各位老師，同學(xué)，達人，幫幫忙，看一下，到底問題出在哪？不勝感激。 下面是初始優(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復(fù)制成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，與此同時，磁矩也變的不準(zhǔn)確了。 我的問題是： 1，有人遇到過這樣的情況嗎：優(yōu)化結(jié)構(gòu)后，算單點能量不收斂的。 2，我的INCAR中有什么是不對的嗎，如果需要我貼出什么其他的信息，請在本帖留言。 真心希望問題得到解決，感激不盡。。。。。 [ Last edited by chord7806 on 2012-5-29 at 16:13 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 第一性 |
» 猜你喜歡
085600材料與化工調(diào)劑 324分
已經(jīng)有5人回復(fù)
-
344求調(diào)劑
已經(jīng)有6人回復(fù)
-
266求調(diào)劑
已經(jīng)有9人回復(fù)
-
化學(xué)工程321分求調(diào)劑
已經(jīng)有18人回復(fù)
-
314求調(diào)劑
已經(jīng)有8人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有5人回復(fù)
-
085601專碩，總分342求調(diào)劑，地區(qū)不限
已經(jīng)有5人回復(fù)
-
311求調(diào)劑
已經(jīng)有6人回復(fù)
-
化工學(xué)碩306求調(diào)劑
已經(jīng)有10人回復(fù)
-
286求調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 結(jié)構(gòu)優(yōu)化計算不收斂，是怎么回事？求幫助
已經(jīng)有6人回復(fù)
- 關(guān)于VASP計算收斂的判斷
已經(jīng)有25人回復(fù)
- 關(guān)于VASP靜態(tài)計算
已經(jīng)有4人回復(fù)
- VASP靜態(tài)自洽計算收斂特慢！
已經(jīng)有31人回復(fù)
- VASP NEB 不收斂問題
已經(jīng)有25人回復(fù)
- VASP計算表面吸附 但是能量總是不收斂/。。。
已經(jīng)有9人回復(fù)
- 用VASP計算NEB都十幾天了還不收斂，急�。。�！
已經(jīng)有8人回復(fù)
- vasp中，做截斷能測試時，測試收斂的標(biāo)準(zhǔn)
已經(jīng)有10人回復(fù)
- 用VASP優(yōu)化結(jié)構(gòu)收斂不了，求指教
已經(jīng)有10人回復(fù)
- 做結(jié)構(gòu)優(yōu)化和做靜態(tài)計算能量不一致
已經(jīng)有7人回復(fù)
- 用Dmol優(yōu)化收斂的結(jié)構(gòu)，用vasp不收斂？
已經(jīng)有9人回復(fù)
- vasp4.6收斂的結(jié)構(gòu)在vasp5.2 計算不收斂
已經(jīng)有7人回復(fù)
- 關(guān)于vasp收斂和優(yōu)化
已經(jīng)有24人回復(fù)
- vasp優(yōu)化難收斂，請大家給我看看
已經(jīng)有5人回復(fù)
- VASP弛豫收斂太慢
已經(jīng)有18人回復(fù)
- vasp優(yōu)化力收斂條件
已經(jīng)有6人回復(fù)
- 求助：vasp結(jié)構(gòu)優(yōu)化不能收斂
已經(jīng)有18人回復(fù)
- vasp計算：優(yōu)化能收斂，但同樣的設(shè)置靜態(tài)計算不收斂
已經(jīng)有3人回復(fù)
- VASP優(yōu)化結(jié)構(gòu)遇到的問題。
已經(jīng)有9人回復(fù)
- 【求助】VASP收斂問題
已經(jīng)有15人回復(fù)
- 【求助】結(jié)構(gòu)優(yōu)化不收斂
已經(jīng)有9人回復(fù)
- 【求助】vasp截斷能不收斂
已經(jīng)有7人回復(fù)
- 【求助】VASP結(jié)構(gòu)優(yōu)化
已經(jīng)有7人回復(fù)
- 【求助】VASP結(jié)構(gòu)優(yōu)化不收斂怎么辦？
已經(jīng)有5人回復(fù)
wanghouyuan
銅蟲 (小有名氣)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 224.1
- 紅花: 1
- 帖子: 75
- 在線: 15.2小時
- 蟲號: 1144660
- 注冊: 2010-11-12
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
2樓2012-05-29 16:50:33
chord7806
新蟲 (初入文壇)
- 應(yīng)助: 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
銅蟲 (小有名氣)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 224.1
- 紅花: 1
- 帖子: 75
- 在線: 15.2小時
- 蟲號: 1144660
- 注冊: 2010-11-12
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
4樓2012-05-29 17:29:21
5樓2018-07-07 10:22:39
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 chord7806 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +5 | llllkkkhh 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 344求調(diào)劑 +5 | knight344 2026-03-16 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷 +10 | 曦熙兮 2026-03-15 | 10/500 |
|
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|