24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 483  |  回复: 0
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 wlscxf 的 10 个金币

wlscxf

铜虫 (小有名气)

[交流] 【求助】vasp计算分子动力学的INCAR文件问题

我用vasp计算金属原子在半导体表面的动力学性质。体系结构先由自洽计算优化得到(K点取7x7x1),然后想看一下低温下体系的动力学性质(这时只取一个k点,为算的快些)。体系的温度从0K逐渐加温到50K,但是计算结果,金属原子在低温下的振荡太剧烈,远远超出了预期,也不符合实际。大家帮忙看一下,是不是INCAR文件设置有问题?INCAR文件内容如下:
general:
   System = 111
   ISTART = 0
# ICHARG = 2         
   PREC = Low
   ENCUT = 300.0
   IBRION = 0
   SMASS = -1
   POTIM = 1.0
   TEBEG = 0.0
   TEEND = 49.2   
   ISYM = 0

   BMIX = 2
   MAXMIX = 50
   APACO = 10.0

# ENAUG = 400.0      
# GGA = PB
# NPAR = 4
   NSW= 1000
   NELM = 100
   NELMIN = 4
   EDIFF = 1E-06
# EDIFFG = -0.01      
# ISIF= 2
# NBANDS = 300  
   ISMEAR = 0
# SIGMA = 0.1
# LWAVE = F
# LCHARG = F
   
# AMIX = 0.2
# BMIX = 0.0001
# AMIX_MAG = 0.8
# BMIX_MAG = 0.0001
   
# Electronic Relaxation 2
# ALGO = V
   IALGO = 38    algorithm
# LDIAG = T    sub-space diagonalisation
   LREAL = A

其中的mix参数设置与否对结果影响不大,LREAL的设置(F or A)对结果稍有影响,但都有振荡。与K点的多少有关系吗?
谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wlscxf 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见