计算化学公社

 找回密码 Forget password
 注册 Register
Views: 405|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] 在使用vasp进行分子动力学模拟时,使用NVE系综,想把温度控制在300K,可还是从0K

[复制链接 Copy URL]

6

帖子

0

威望

61

eV
积分
67

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 Zfxdo 于 2025-3-18 14:32 编辑

各位老师好,在使用vasp进行分子动力学模拟时(九十多个原子,In2Se3),使用NVE系综,想把温度控制在300K,可还是从0K开始上升,初次使用,还希望多多指导和帮助,以下是我的INCAR内容
# 基础计算控制
ISTART = 0        # 新计算
ICHARG = 2        # 原子电荷初始化
PREC = Accurate   # 高精度模式
ENCUT = 500       # 截断能
LWAVE = F
LCHARG  = F
# 电子自洽收敛
EDIFF = 1E-6      # 电子步收敛标准
NELM = 100        # 最大电子步数
ALGO = Normal     # 标准电子优化算法


# 结构优化参数
MDALGO = 1
ANDERSEN_PROB = 0
POTIM  = 0.5
IBRION = 0        
ISIF = 3          # 优化晶胞形状和原子位置(二维材料可改用ISIF=4固定垂直方向)
NSW = 1000        # 最大离子步数
TEBEG=300         #初始温度控制开关,单位为K
TEEND=300         #弛豫最终温度控制开关,默认值与TEBEG相同
#EDIFFG = -0.02     # 离子步收敛标准(负值表示力收敛,单位eV/?)
#IOPTCELL = 1 1 0 1 1 0 0 0 0
# 交换关联泛函
GGA = PE          # PBE泛函
IVDW = 12         # DFT-D3(BJ)范德华修正(关键!二维层间作用依赖vdW)


# `其他优化
ISMEAR = 0        # Gaussian展宽(金属/半导体均可)
SIGMA = 0.05      # 展宽宽度(eV)
LORBIT = 11       # 输出局域轨道信息(可选)
LREAL = A         # 实空间投影(加速计算)
NCORE  = 8
SYMPREC = 1E-6
#ISPIN = 2
#MGMMOM = 5*0
#ISYM  = 0
LDIPOL = T
IDIPOL = 3
DIPOL  = 0.5 0.5 0.5 "
以下是我目前OSZICAR某一离子步的输出结果” 38 T=     1. E= -.35584486E+03 F= -.35585520E+03 E0= -.35583020E+03  EK= 0.10340E-01 SP= 0.00E+00 SK= 0.00E+00

问题:我想要保持恒温,虽然设置了初始和结束温度为300K,为什么还是从0K开始?我现在用的是NVE系综吗?如何设置可以使用NpT系综?我的目的是要数据集进行机器学习,那在进行分子动力学模拟时选择哪个系综比较合适?

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 02:13 , Processed in 0.156143 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list