| 7 | 1/1 | 返回列表 |
| 查看: 2593 | 回復(fù): 6 | |||
[求助]
求教各位大俠,關(guān)于vasp的運(yùn)行記錄
|
|
如題,以下是vasp的運(yùn)行記錄,想就運(yùn)行記錄請(qǐng)教各位大俠幾個(gè)問(wèn)題: 1. | For optimal performance we recommend that you set | | NPAR = approx SQRT( number of cores) | | This will greatly improve the performance of VASP for DFT. 這段的意思是要設(shè)置NPAR嗎?我沒(méi)有設(shè)置但是我的八個(gè)核仍然還在運(yùn)轉(zhuǎn)啊 2.WARNING: small aliasing (wrap around) errors must be expected 為什么會(huì)出現(xiàn)這個(gè)警告? 3. 為什么經(jīng)過(guò)了一段: N E dE d eps ncg rms rms(c) DAV: 1 0.395219316642E+03 0.39522E+03 -0.18815E+04 5424 0.184E+03 …… DAV: 20 -0.742393695060E+02 -0.51618E-04 -0.48508E-06 3728 0.101E-02 之后,緊接著又來(lái)了一段 N E dE d eps ncg rms rms(c) DAV: 1 -0.742916604512E+02 -0.52343E-01 -0.17052E+01 5424 0.270E+01 0.110E+00 …… DAV: 7 -0.742797664941E+02 0.15331E-04 -0.27305E-04 4544 0.115E-01 此時(shí)電子的能量收斂標(biāo)準(zhǔn)(EDIFF=1E-4)已經(jīng)達(dá)到了,為何還要再?gòu)腄AV1開(kāi)始計(jì)算,而且重復(fù)數(shù)次? 1 F= -.74239370E+02 E0= -.74167302E+02 d E =-.742394E+02 curvature: 0.00 expect dE= 0.000E+00 dE for cont linesearch 0.000E+00 trial: gam= 0.00000 g(F)= 0.308E+00 g(S)= 0.190E-02 ort = 0.000E+00 (trialstep = 0.100E+01) search vector abs. value= 0.310E+00 bond charge predicted 這一段話包含了那些意義呢? 請(qǐng)各位指教 下面是log文件: ----------------------------------------------------------------------------- | | | W W AA RRRRR N N II N N GGGG !!! | | W W A A R R NN N II NN N G G !!! | | W W A A R R N N N II N N N G !!! | | W WW W AAAAAA RRRRR N N N II N N N G GGG ! | | WW WW A A R R N NN II N NN G G | | W W A A R R N N II N N GGGG !!! | | | | For optimal performance we recommend that you set | | NPAR = approx SQRT( number of cores) | | This will greatly improve the performance of VASP for DFT. | | The default NPAR=number of cores might be grossly inefficient | | on modern multi-core architectures or massively parallel machines. | | Unfortunately you need to use the default for hybrid, GW and RPA | | calculations. | | | ----------------------------------------------------------------------------- LDA part: xc-table for Pade appr. of Perdew POSCAR, INCAR and KPOINTS ok, starting setup WARNING: small aliasing (wrap around) errors must be expected FFT: planning ...( 1 ) WAVECAR not read entering main loop N E dE d eps ncg rms rms(c) DAV: 1 0.395219316642E+03 0.39522E+03 -0.18815E+04 5424 0.184E+03 …… DAV: 20 -0.742393695060E+02 -0.51618E-04 -0.48508E-06 3728 0.101E-02 1 F= -.74239370E+02 E0= -.74167302E+02 d E =-.742394E+02 curvature: 0.00 expect dE= 0.000E+00 dE for cont linesearch 0.000E+00 trial: gam= 0.00000 g(F)= 0.308E+00 g(S)= 0.190E-02 ort = 0.000E+00 (trialstep = 0.100E+01) search vector abs. value= 0.310E+00 bond charge predicted N E dE d eps ncg rms rms(c) DAV: 1 -0.742916604512E+02 -0.52343E-01 -0.17052E+01 5424 0.270E+01 0.110E+00 …… DAV: 7 -0.742797664941E+02 0.15331E-04 -0.27305E-04 4544 0.115E-01 2 F= -.74279766E+02 E0= -.74207676E+02 d E =-.403970E-01 trial-energy change: -0.040397 1 .order -0.045188 -0.309599 0.219223 step: 0.5721(harm= 0.5854) dis= 0.01659 next Energy= -74.327040 (dE=-0.877E-01) bond charge predicted N E dE d eps ncg rms rms(c) DAV: 1 -0.743293982141E+02 -0.49616E-01 -0.31287E+00 5424 0.115E+01 0.471E-01 …… DAV: 5 -0.743270296764E+02 0.96815E-04 -0.92136E-04 5512 0.179E-01 3 F= -.74327030E+02 E0= -.74254782E+02 d E =-.876602E-01 curvature: -0.29 expect dE=-0.376E-02 dE for cont linesearch -0.130E-05 trial: gam= 0.03777 g(F)= 0.118E-01 g(S)= 0.104E-02 ort = 0.117E-02 (trialstep = 0.914E+00) search vector abs. value= 0.134E-01 bond charge predicted N E dE d eps ncg rms rms(c) DAV: 1 -0.743306820981E+02 -0.35556E-02 -0.56505E-01 5424 0.517E+00 0.410E-01 …… DAV: 5 -0.743301754215E+02 0.21081E-04 -0.20817E-04 5384 0.846E-02 4 F= -.74330175E+02 E0= -.74257952E+02 d E =-.314575E-02 trial-energy change: -0.003146 1 .order -0.003170 -0.011806 0.005467 step: 0.6250(harm= 0.6250) dis= 0.00454 next Energy= -74.331064 (dE=-0.403E-02) bond charge predicted N E dE d eps ncg rms rms(c) DAV: 1 -0.743310436880E+02 -0.84719E-03 -0.55433E-02 5424 0.162E+00 0.134E-01 …… DAV: 4 -0.743309981039E+02 0.23216E-04 -0.76480E-05 3560 0.623E-02 5 F= -.74330998E+02 E0= -.74258764E+02 d E =-.396843E-02 curvature: -0.32 expect dE=-0.616E-03 dE for cont linesearch -0.130E-06 trial: gam= 0.14599 g(F)= 0.187E-02 g(S)= 0.263E-04 ort =-0.732E-04 (trialstep = 0.857E+00) search vector abs. value= 0.216E-02 bond charge predicted N E dE d eps ncg rms rms(c) DAV: 1 -0.743322307005E+02 -0.12094E-02 -0.58257E-02 5424 0.178E+00 0.105E-01 …… DAV: 4 -0.743322525693E+02 0.37658E-05 -0.28532E-05 3608 0.353E-02 6 F= -.74332253E+02 E0= -.74260047E+02 d E =-.125447E-02 trial-energy change: -0.001254 1 .order -0.001259 -0.001617 -0.000901 step: 1.9360(harm= 1.9360) dis= 0.00553 next Energy= -74.332825 (dE=-0.183E-02) bond charge predicted N E dE d eps ncg rms rms(c) DAV: 1 -0.743327647452E+02 -0.50841E-03 -0.92389E-02 5424 0.224E+00 0.137E-01 …… DAV: 4 -0.743328070278E+02 0.79076E-05 -0.48127E-05 3608 0.481E-02 7 F= -.74332807E+02 E0= -.74260611E+02 d E =-.180892E-02 curvature: -1.11 expect dE=-0.109E-02 dE for cont linesearch -0.649E-07 trial: gam= 0.52253 g(F)= 0.930E-03 g(S)= 0.516E-04 ort =-0.113E-04 (trialstep = 0.107E+01) search vector abs. value= 0.156E-02 bond charge predicted N E dE d eps ncg rms rms(c) DAV: 1 -0.743332034028E+02 -0.38847E-03 -0.77433E-02 5424 0.204E+00 0.144E-01 …… DAV: 4 -0.743331834775E+02 0.23471E-04 -0.78441E-05 3616 0.614E-02 8 F= -.74333183E+02 E0= -.74260994E+02 d E =-.376450E-03 trial-energy change: -0.000376 1 .order -0.000399 -0.001047 0.000248 step: 0.8668(harm= 0.8668) dis= 0.00219 next Energy= -74.333230 (dE=-0.423E-03) bond charge predicted N E dE d eps ncg rms rms(c) DAV: 1 -0.743332115200E+02 -0.45712E-05 -0.30759E-03 5424 0.403E-01 0.272E-02 DAV: 2 -0.743332168524E+02 -0.53324E-05 -0.72029E-05 3688 0.613E-02 9 F= -.74333217E+02 E0= -.74260994E+02 d E =-.409825E-03 curvature: -0.69 expect dE=-0.534E-04 dE for cont linesearch -0.146E-06 trial: gam= 0.11587 g(F)= 0.560E-04 g(S)= 0.210E-04 ort =-0.181E-04 (trialstep = 0.103E+01) search vector abs. value= 0.938E-04 reached required accuracy - stopping structural energy minimisation writing wavefunctions ![]() ![]() |
第一性原理相關(guān)文檔 |

金蟲(chóng) (小有名氣)
|
1.NPAR是并行設(shè)置,多個(gè)核一起計(jì)算的時(shí)候最好設(shè)置成總核數(shù)開(kāi)根號(hào),這樣計(jì)算效率更高,但是不設(shè)也無(wú)妨,而且如果遇到開(kāi)根號(hào)不是整數(shù)的時(shí)候也可以選擇使用別的數(shù),比如32個(gè)核設(shè)成4也是可以算的。 2.起因恐怕是很深層次的量子力學(xué)或者固體物理方面的問(wèn)題,我們也不是很清楚,但是可以通過(guò)INCAR參數(shù)調(diào)整得更精密來(lái)避免,不過(guò)一般情況下可以忽略這點(diǎn)提示,對(duì)結(jié)果影響不會(huì)太大。 3.我不是太明白是什么情況。看底下的LOG文件似乎挺正常的。DAV1…20是表示完成了20個(gè)電子迭代步,然后收斂了,于是有了1F=…也就是說(shuō)完成了第1個(gè)離子迭代步,然后進(jìn)入下一個(gè)循環(huán),也就是第2個(gè)離子步,這里面有包括了7次電子步。具體的其他信息可以參照VASP手冊(cè)里面關(guān)于OSZICAR的描述。 |
|
3.我不是太明白是什么情況?吹紫碌腖OG文件似乎挺正常的。DAV1…20是表示完成了20個(gè)電子迭代步,然后收斂了,于是有了1F=…也就是說(shuō)完成了第1個(gè)離子迭代步,然后進(jìn)入下一個(gè)循環(huán),也就是第2個(gè)離子步,這里面有包括了7次電子步。具體的其他信息可以參照VASP手冊(cè)里面關(guān)于OSZICAR的描述。 前輩,是不是算完第一個(gè)離子迭代步循環(huán),這時(shí)候還沒(méi)有達(dá)到EDIFFG的收斂標(biāo)準(zhǔn),所以要繼續(xù)進(jìn)行計(jì)算?換句話說(shuō),如果EDIFFG的標(biāo)準(zhǔn)越高,那么這種循環(huán)的次數(shù)也就越多?不知道您是否知道我要表達(dá)的意思 |

金蟲(chóng) (小有名氣)
|
應(yīng)該是EDIFF和EDIFFG兩個(gè)參數(shù)吧,默認(rèn)的話前者是后者的10倍。當(dāng)電子迭代到能量差小于EDIFFG的時(shí)候電子迭代結(jié)束,完成第一個(gè)離子步,但是此時(shí)你的體系的離子能量波動(dòng)還很大( d E =-.742394E+02),還遠(yuǎn)沒(méi)有達(dá)到EDIFF要求的標(biāo)準(zhǔn),于是乎再將這個(gè)結(jié)構(gòu)進(jìn)行微擾,再一個(gè)個(gè)電子步地迭代,直到這次又迭代了7次終于達(dá)到了EDIFFG的標(biāo)準(zhǔn),第2個(gè)離子步完成,但是離子步還沒(méi)完(d E =-.403970E-01),所以還要這樣周而復(fù)始地算,直到9 F= -.74333217E+02 E0= -.74260994E+02 d E =-.409825E-03終于達(dá)到了EDIFF的標(biāo)準(zhǔn),于是結(jié)束了 |

金蟲(chóng) (小有名氣)

| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽(yáng)CF 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +5 | S.H1 2026-03-18 | 5/250 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 9/450 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +7 | rare12345 2026-03-18 | 7/350 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么! 2026-03-15 | 18/900 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂(lè) 2026-03-14 | 5/250 |
|
|
[考研] 308 085701 四六級(jí)已過(guò)求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
[論文投稿]
投稿問(wèn)題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|