24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1251  |  回复: 2

wuli_ml

银虫 (正式写手)

[求助] VASP计算光学性质,怎么调整能量的步长已有1人参与

用vasp计算吸收光谱,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 !

计算的吸收光谱,部分如下:
  #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

这个吸收光谱的能量步长是0.0678eV,怎么调整这个步长到0.03eV?
回复此楼
于千万年之间,于千万人之中,早一步,晚一步,刚刚好遇到。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyikeco

木虫 (正式写手)

newbie

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wuli_ml: 金币+20, ★★★★★最佳答案, 已解决问的问题 2023-11-25 14:33:07
在VASP中,吸收光谱的计算通常是通过在能量范围内进行一系列的静态计算来实现的。为了调整能量步长,你可以在计算前增加更多的能量点,以确保所需的步长。这可以通过更改 INCAR 文件中的 NEDOS 参数来实现。

在你的 INCAR 文件中,NEDOS 参数的值是2000。你可以通过增加这个值来增加所计算的能量点的数量,从而减小能量步长。这样,你就可以获得更密集的吸收光谱数据。

例如,将 NEDOS 的值增加到4000:
NEDOS = 4000
请记住,增加 NEDOS 的值将增加计算的时间,因为需要在更多的能量点上进行计算。根据你的计算需求和计算资源,你可能需要权衡计算精度和计算时间。

调整 NEDOS 参数后,重新运行你的计算。在计算完成后,你将得到一个具有更密集能量点的吸收光谱。
2楼2023-11-24 22:43:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli_ml

银虫 (正式写手)

引用回帖:
2楼: Originally posted by wangyikeco at 2023-11-24 22:43:50
在VASP中,吸收光谱的计算通常是通过在能量范围内进行一系列的静态计算来实现的。为了调整能量步长,你可以在计算前增加更多的能量点,以确保所需的步长。这可以通过更改 INCAR 文件中的 NEDOS 参数来实现。

在你 ...

大神,增加NEDOS已经可以解决,有时候跑光学性质只有301个点(NEDOS也设置了2000)。
能否帮忙看下这个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楼2023-11-25 14:31:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuli_ml 的主题更新
信息提示
请填处理意见