| 查看: 1251 | 回复: 2 | ||
[求助]
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? ![]() |
» 猜你喜欢
2025冷门绝学什么时候出结果
已经有5人回复
Bioresource Technology期刊,第一次返修的时候被退回好几次了
已经有7人回复
真诚求助:手里的省社科项目结项要求主持人一篇中文核心,有什么渠道能发核心吗
已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
孩子确诊有中度注意力缺陷
已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
请问有评职称,把科研教学业绩算分排序的高校吗
已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复

wangyikeco
木虫 (正式写手)
newbie
- 应助: 84 (初中生)
- 金币: 6662
- 帖子: 524
- 在线: 107.1小时
- 虫号: 25598800
- 注册: 2021-03-03
- 性别: GG
- 专业: 凝聚态物性 II :电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wuli_ml: 金币+20, ★★★★★最佳答案, 已解决问的问题 2023-11-25 14:33:07
感谢参与,应助指数 +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
|
大神,增加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














回复此楼