计算化学公社

标题: 在使用vasp进行分子动力学模拟时，使用NVE系综，想把温度控制在300K，可还是从0K [打印本页]

作者
Author: Zfxdo 时间: 2025-3-18 14:31
标题: 在使用vasp进行分子动力学模拟时，使用NVE系综，想把温度控制在300K，可还是从0K
本帖最后由 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系综？我的目的是要数据集进行机器学习，那在进行分子动力学模拟时选择哪个系综比较合适？

欢迎光临计算化学公社 (http://bbs.keinsci.com/)