| 5 | 1/1 | 返回列表 |
| 查看: 4475 | 回復(fù): 4 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
chord7806新蟲 (初入文壇)
|
[求助]
<求助>VASP 結(jié)構(gòu)優(yōu)化收斂,靜態(tài)計(jì)算不收斂。(INCAR KPOINTS POSCAR附上) 已有1人參與
|
||
|
最近計(jì)算了一個體系,共40原子。兩種磁性離子。 先進(jìn)行了結(jié)構(gòu)優(yōu)化。得到了CONTCAR,然后cp CONTCAR POSCAR靜態(tài)計(jì)算能量始終不收斂。 希望各位老師,同學(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 計(jì)算完以后有這樣的提示: 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)計(jì)算,問題出現(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)后,算單點(diǎn)能量不收斂的。 2,我的INCAR中有什么是不對的嗎,如果需要我貼出什么其他的信息,請?jiān)诒咎粞浴?br /> 真心希望問題得到解決,感激不盡。。。。。 [ Last edited by chord7806 on 2012-5-29 at 16:13 ] |
第一性 |
新蟲 (初入文壇)
|
這個是我靜態(tài)計(jì)算的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)注。。 |
銅蟲 (小有名氣)

銅蟲 (小有名氣)


| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績338 +7 | YYYYYNNNNN 2026-03-26 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 367求調(diào)劑 +5 | 芋泥啵啵… 2026-03-28 | 5/250 |
|
|
[考研] 一志愿:西北大學(xué),英一數(shù)一408-284分求調(diào)劑 +5 | 12.27 2026-03-27 | 5/250 |
|
|
[考研] 085601材料工程找調(diào)劑 +17 | oatmealR 2026-03-29 | 18/900 |
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
[考研] 329求調(diào)劑 +10 | 鈕恩雪 2026-03-25 | 10/500 |
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
[考研] 求調(diào)劑 +7 | 爭取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研]
|
黃粱一夢千年 2026-03-24 | 3/150 |
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|