24小时热门版块排行榜    

查看: 1565  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否退出参与的面上项目解除限项 +11 koalala 2026-08-24 13/650 2026-08-24 10:49 by anjeeshine
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 明天应该可查了!? +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 范进中举一文的中心思想 +5 炎黄贵胄 2026-08-22 6/300 2026-08-24 06:34 by sincosx
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见