| 9 | 1/1 | 返回列表 |
| 查看: 5317 | 回復(fù): 8 | ||
| 【懸賞金幣】回答本帖問(wèn)題,作者大力2010將贈(zèng)送您 10 個(gè)金幣 | ||
[求助]
VASP結(jié)構(gòu)優(yōu)化以 writing wavefunctions 結(jié)束。 是不是表示不收斂? 已有2人參與
|
||
|
請(qǐng)教大神,我用vasp對(duì)摻雜了 v 和 ta 的異質(zhì)結(jié)進(jìn)行結(jié)構(gòu)優(yōu)化,incar文件時(shí)這樣的: global parameters istart = 1 (read existing wavefunction; if there) ispin = 1 (non-spin polarised dft) # icharg = 11 (non-self-consistent: gga/lda band structures) lreal = .false. (projection operators: automatic) # encut = 400 (cut-off energy for plane wave basis set, in ev) prec = normal (precision level) lwave = .true. (write wavecar or not) lcharg = .true. (write chgcar or not) addgrid= .true. (increase grid; helps gga convergence) # lvtot = .true. (write total electrostatic potential into locpot or not) # lvhar = .true. (write ionic + hartree electrostatic potential into locpot or not) # nelect = (no. of electrons: charged cells; be careful) # lplane = .true. (real space distribution; supercells) # npar = 4 (max is no. nodes; don't set for hybrids) # nwrite = 2 (medium-level output) # kpar = 2 (divides k-grid into separate groups) # ngx = 500 (fft grid mesh density for nice charge/potential plots) # ngy = 500 (fft grid mesh density for nice charge/potential plots) # ngz = 500 (fft grid mesh density for nice charge/potential plots) ivdw = 11 electronic relaxation ismear = 0 (gaussian smearing; metals:1) sigma = 0.05 (smearing value in ev; metals:0.2) nelm = 90 (max electronic scf steps) nelmin = 6 (min electronic scf steps) ediff = 1e-05 (scf energy convergence; in ev) # gga = ps (pbesol exchange-correlation) ionic relaxation nsw = 100 (max ionic steps) ibrion = 2 (algorithm: 0-md; 1-quasi-new; 2-cg) isif = 2 (stress/relaxation: 2-ions, 3-shape/ions/v, 4-shape/ions) ediffg = -2e-02 (ionic convergence; ev/aa) # isym = 2 (symmetry: 0=none; 2=gga; 3=hybrids) 優(yōu)化結(jié)束時(shí),沒(méi)出現(xiàn)“reached required accurancy”,以 writing wavefunctions 結(jié)束: 100 f= -.82977297e+03 e0= -.82954072e+03 d e =-.426037e-01 curvature: -2.30 expect de=-0.187e+00 de for cont linesearch -0.376e-06 trial: gam= 1.72200 g(f)= 0.812e-01 g(s)= 0.000e+00 ort = 0.154e-03 (trialstep = 0.125e+00) search vector abs. value= 0.514e+00 writing wavefunctions |
新蟲(chóng) (著名寫(xiě)手)
新蟲(chóng) (著名寫(xiě)手)
|
以計(jì)算結(jié)果作為初始結(jié)構(gòu),重新計(jì)算。是不是你的步驟比較少,沒(méi)計(jì)算到那一步? 發(fā)自小木蟲(chóng)Android客戶端 |
新蟲(chóng) (初入文壇)
金蟲(chóng) (小有名氣)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué) 0703求調(diào)劑 總分293 一志愿211 +8 | 土土小蟲(chóng) 2026-03-03 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 0703 東華大學(xué) 理學(xué)化學(xué)均過(guò) a 線,總分 281 求調(diào)劑 +3 | 香香 Lu 2026-03-03 | 6/300 |
|
|
[考研] 070300化學(xué) 280 一志愿太原理工 求調(diào)劑 +7 | 拾玖壹 2026-03-04 | 7/350 |
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
[考研] 085601 材料工程 320 +6 | 和樂(lè)瑤 2026-03-03 | 6/300 |
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|
|
[考博] 26申博 +5 | north, 2026-02-28 | 5/250 |
|
|
[考博] 26申博-目前4篇SCI一作 +3 | chen_2024 2026-03-02 | 3/150 |
|
|
[考研] 289求調(diào)劑 +8 | yang婷 2026-03-02 | 10/500 |
|
|
[考研] 0854總分272 +5 | 打小就是老實(shí)人 2026-03-02 | 6/300 |
|
|
[考研] 293求調(diào)劑 +4 | 是樂(lè)渝哇 2026-03-03 | 4/200 |
|
|
[考研] 268求調(diào)劑 +6 | 好運(yùn)連綿不絕 2026-03-02 | 6/300 |
|
|
[考研] 課題組接收材料類(lèi)調(diào)劑研究生 +6 | gaoxiaoniuma 2026-02-28 | 9/450 |
|
|
[考研] 0854復(fù)試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|
|
[考研] 275求調(diào)劑 +7 | 明遠(yuǎn)求學(xué) 2026-03-01 | 7/350 |
|
|
[考博] 博士自薦 +4 | kkluvs 2026-02-28 | 5/250 |
|
|
[考研] 265分求調(diào)劑不調(diào)專(zhuān)業(yè)和學(xué)校有行學(xué)上就 +6 | 禮堂丁真258 2026-02-28 | 9/450 |
|
|
[考研] 284求調(diào)劑 +10 | 天下熯 2026-02-28 | 11/550 |
|
|
[考研] 化工專(zhuān)碩342,一志愿大連理工大學(xué),求調(diào)劑 +6 | kyf化工 2026-02-28 | 7/350 |
|
|
[考研] 311求調(diào)劑 +9 | 南迦720 2026-02-28 | 10/500 |
|