|
本帖最后由 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系综?我的目的是要数据集进行机器学习,那在进行分子动力学模拟时选择哪个系综比较合适?
|
|