| 查看: 1376 | 回复: 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？  |
回复此楼» 猜你喜欢
297，工科调剂?
已经有4人回复
-
恳请有学校收留
已经有7人回复
-
291求调剂
已经有9人回复
-
300求调剂
已经有11人回复
-
22专硕求调剂
已经有12人回复
-
材料相关专业344求调剂双非工科学校或课题组
已经有25人回复
-
急需调剂
已经有7人回复
-
求调剂
已经有10人回复
-
一志愿华中农业071010，320求调剂
已经有16人回复
-
304求调剂
已经有5人回复
wangyikeco
木虫 (正式写手)
newbie
- 应助: 84 (初中生)
- 金币: 7152.5
- 帖子: 549
- 在线: 112.7小时
- 虫号: 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
