| 3 | 1/1 | 返回列表 |
| 查看: 1332 | 回復(fù): 2 | |||
[求助]
VASP計(jì)算光學(xué)性質(zhì),怎么調(diào)整能量的步長 已有1人參與
|
|
用vasp計(jì)算吸收光譜,INCAR如下: # Basic setup: ISTART = 0 ICHARG = 2 ISPIN = 1 PREC = Accurate ENCUT = 520 ALGO = Normal EDIFF = 1E-05 NELM = 100 ISMEAR = 0 SIGMA = 0.05 ISIF = 2 IBRION = -1 NSW = 0 LCHARG = .TRUE. LWAVE = .TRUE. LELF = .FALSE. LVTOT = .FALSE. LOPTICS = TRUE ! NEDOS = 2000 CSHIFT = 0.1 ! NBANDS = 120 ! 計(jì)算的吸收光譜,部分如下: #Energy xx(%) yy(%) 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.67800000E-01 0.75959367E+00 0.75969136E+00 0.13560000E+00 0.76050169E+00 0.76061024E+00 0.20330000E+00 0.76167577E+00 0.76177341E+00 0.27110000E+00 0.76391894E+00 0.76400573E+00 0.33890000E+00 0.76675176E+00 0.76686025E+00 0.40670000E+00 0.77025761E+00 0.77038779E+00 0.47450000E+00 0.77446585E+00 0.77454177E+00 0.54230000E+00 0.77931500E+00 0.77940176E+00 0.61000000E+00 0.78484194E+00 0.78493951E+00 這個(gè)吸收光譜的能量步長是0.0678eV,怎么調(diào)整這個(gè)步長到0.03eV? |

木蟲 (正式寫手)
newbie
|
在VASP中,吸收光譜的計(jì)算通常是通過在能量范圍內(nèi)進(jìn)行一系列的靜態(tài)計(jì)算來實(shí)現(xiàn)的。為了調(diào)整能量步長,你可以在計(jì)算前增加更多的能量點(diǎn),以確保所需的步長。這可以通過更改 INCAR 文件中的 NEDOS 參數(shù)來實(shí)現(xiàn)。 在你的 INCAR 文件中,NEDOS 參數(shù)的值是2000。你可以通過增加這個(gè)值來增加所計(jì)算的能量點(diǎn)的數(shù)量,從而減小能量步長。這樣,你就可以獲得更密集的吸收光譜數(shù)據(jù)。 例如,將 NEDOS 的值增加到4000: NEDOS = 4000 請記住,增加 NEDOS 的值將增加計(jì)算的時(shí)間,因?yàn)樾枰诟嗟哪芰奎c(diǎn)上進(jìn)行計(jì)算。根據(jù)你的計(jì)算需求和計(jì)算資源,你可能需要權(quán)衡計(jì)算精度和計(jì)算時(shí)間。 調(diào)整 NEDOS 參數(shù)后,重新運(yùn)行你的計(jì)算。在計(jì)算完成后,你將得到一個(gè)具有更密集能量點(diǎn)的吸收光譜。 |
|
大神,增加NEDOS已經(jīng)可以解決,有時(shí)候跑光學(xué)性質(zhì)只有301個(gè)點(diǎn)(NEDOS也設(shè)置了2000)。 能否幫忙看下這個(gè)INCAR,具體啥問題? 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 = 500 (Cut-off energy for plane wave basis set, in eV) # PREC = Accurate LWAVE = .TRUE. (Write WAVECAR or not) LCHARG = .TRUE. (Write CHGCAR or not) Electronic Relaxation ISMEAR = 0 (Gaussian smearing, metals:1) SIGMA = 0.02 (Smearing value in eV, metals:0.2) NELM = 90 (Max electronic SCF steps) NELMIN = 6 (Min electronic SCF steps) EDIFF = 1E-06 (SCF energy convergence, in eV) Ionic Relaxation NSW = 0 (Max ionic steps) IBRION = -1 (Algorithm: 0-MD, 1-Quasi-New, 2-CG) ISIF = 2 (Stress/relaxation: 2-Ions, 3-Shape/Ions/V, 4-Shape/Ions) EDIFFG = -1E-02 (Ionic convergence, eV/AA) IVDW = 11 LORBIT = 11 NPAR = 4 ENDOS = 2000 #EMIN =-20 #EMXA =20 NBANDS = 80 LOPTICS = .TRUE. POTIM = 0.05 AMIX = 0.1 BMIX = 0.001 |

| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中科大材料與化工。353分能調(diào)劑到哪? +5 | MayUxw1 2026-03-08 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-10 | 9/450 |
|
|
[考研] 291分獸醫(yī)專碩 一志愿華中農(nóng) +3 | 種大大太陽 2026-03-05 | 3/150 |
|
|
[考研] 279求調(diào)劑 +3 | 莫xiao 2026-03-10 | 4/200 |
|
|
[考研] 085600 材料與化工 295 求調(diào)劑 +9 | dream…… 2026-03-10 | 9/450 |
|
|
[考博] 26申博求助 +3 | 跳躍餅干 2026-03-10 | 4/200 |
|
|
[考研] 337一志愿華南理工材料求調(diào)劑 +5 | mysdl 2026-03-07 | 5/250 |
|
|
[考研] 環(huán)境調(diào)劑 +5 | 曉看天暮看云 2026-03-09 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 很老實(shí)人 2026-03-09 | 4/200 |
|
|
[考研] 293一志愿華東理工 0817化學(xué)工程與技術(shù) 調(diào)劑 +5 | fjj0912 2026-03-07 | 5/250 |
|
|
[考研] 理學(xué)287分求調(diào)劑 +4 | 看看我. 2026-03-05 | 4/200 |
|
|
[考研] 347求調(diào)劑 +4 | 浮云滿足 2026-03-07 | 4/200 |
|
|
[考研] 材料調(diào)劑 +7 | ounce. 2026-03-05 | 15/750 |
|
|
[考研] 081700學(xué)碩一志愿北京化工大學(xué)數(shù)二英一過六級有競賽求調(diào)劑 +5 | galaxary 2026-03-07 | 7/350 |
|
|
[考研] 一志愿鄭大071000分?jǐn)?shù)282求調(diào)劑 +3 | 研研顏 2026-03-05 | 7/350 |
|
|
[考研] 304求調(diào)劑 +4 | 52hz~~ 2026-03-05 | 5/250 |
|
|
[考研]
|
zbcm_zbcm 2026-03-05 | 6/300 |
|
|
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過來 +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|
|
[考研] 紡織、生物、化學(xué)、材料等專業(yè) +3 | Eember. 2026-03-05 | 7/350 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財(cái) 2026-03-04 | 3/150 |
|