计算化学公社

标题: 求助：quench-annealing dynamics的过程应该如何设置npt [打印本页]

作者
Author: 邋遢大王 时间: 2023-6-20 14:22
标题: 求助：quench-annealing dynamics的过程应该如何设置npt
本帖最后由邋遢大王于 2023-6-20 14:24 编辑

重复文献中电解液体系的MD，其中有一道工序五次quench-annealing dynamics，温度从298 K到 894 K之间循环，但是模拟过程中模型盒子的尺寸大小扩大非常快，且无法回到初始状态。求教是不是应该采用npt系综，以及应该如何设置，模型尺寸的扩张如何解决呢
以下是相关部分lammps代码：
variable a loop 1 5
label loop_start

fix fxnpt all npt temp 894 298 $(100.0*dt) iso 0.98692 0.98692 $(1000.0*dt)
run 800000

next a
jump SELF loop_start
write_data system_after_cycle.data
unfix fxnpt
print "finish cycle"

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