 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4477 | 回復(fù): 4 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
chord7806新蟲 (初入文壇)
|
[求助]
<求助>VASP 結(jié)構(gòu)優(yōu)化收斂，靜態(tài)計算不收斂。（INCAR KPOINTS POSCAR附上） 已有1人參與
|
||
|
最近計算了一個體系，共40原子。兩種磁性離子。 先進(jìn)行了結(jié)構(gòu)優(yōu)化。得到了CONTCAR，然后cp CONTCAR POSCAR靜態(tài)計算能量始終不收斂。 希望各位老師，同學(xué)，達(dá)人，幫幫忙，看一下，到底問題出在哪？不勝感激。 下面是初始優(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ù)此樓» 收錄本帖的淘帖專輯推薦
 第一性 |
» 猜你喜歡
288求調(diào)劑 一志愿哈工大 材料與化工
已經(jīng)有16人回復(fù)
-
0805求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑，一志愿北林食品與營養(yǎng)095500，301分，已過六級，有科研經(jīng)歷
已經(jīng)有3人回復(fù)
-
0856材料與化工調(diào)劑，339
已經(jīng)有5人回復(fù)
-
311（085601）求調(diào)劑
已經(jīng)有10人回復(fù)
-
080200專業(yè)277分，求帶走！
已經(jīng)有4人回復(fù)
-
359求調(diào)劑
已經(jīng)有6人回復(fù)
-
342求調(diào)劑
已經(jīng)有6人回復(fù)
-
070300化學(xué)279求調(diào)劑
已經(jīng)有10人回復(fù)
-
086000生物與醫(yī)藥298調(diào)劑求助
已經(jīng)有5人回復(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é)
4樓2012-05-29 17:29:21
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
5樓2018-07-07 10:22:39
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京科技大學(xué)材料學(xué)碩328分求調(diào)劑 +4 | 1段時間 2026-03-31 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 290分調(diào)劑求助 +10 | 吉祥止止陳 2026-03-25 | 10/500 |
|
|
[考研] 254材料與化工求調(diào)劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +9 | akdhjs 2026-03-31 | 10/500 |
|
|
|
[考研] 362求調(diào)劑 +9 | 西南交材料專碩3 2026-03-31 | 9/450 |
|
|
|
[考研] 343求調(diào)劑 +8 | 愛羈絆 2026-03-28 | 8/400 |
|
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +10 | 26要上岸 2026-03-28 | 10/500 |
|
|
|
[考研] 0856調(diào)劑 +6 | 曲聽筠 2026-03-30 | 6/300 |
|
|
|
[考研] 335求調(diào)劑 +3 | 321* 2026-03-31 | 3/150 |
|
|
|
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +10 | Koxui 2026-03-28 | 10/500 |
|
|
|
[考研]
 105500藥學(xué)求調(diào)劑，一志愿山東大學(xué)藥學(xué)，348分
+3
|
gr哈哈哈 2026-03-28 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
|
[考研] 0856求調(diào)劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
|
[考研]
一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑
+10
|
y7czhao 2026-03-26 | 10/500 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|