求助：quench-annealing dynamics的过程应该如何设置npt

邋遢大王

重复文献中电解液体系的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"