24小时热门版块排行榜    

查看: 1598  |  回复: 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

禁虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wuli_ml: 金币+20, ★★★★★最佳答案, 已解决问的问题 2023-11-25 14:33:07
本帖内容被屏蔽

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 16:26 by CgyNCDVNhGVg
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 16:15 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-15 15:50 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:50 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-15 14:14 by VMDqgnlKimtZ
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-15 14:02 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:38 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 01:21 by mibUvS8DDCwf
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 23:13 by vZfe6xYu34yj
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
信息提示
请填处理意见