模拟降温，有大致两种办法，一个是连续（升）降温，也就是dt间隔下不停的温度下降（升）；第二种方法是准静态模拟降温，每降低一个温度，跑很长时间；这两种方法大家可以根据自己情况进行模拟。
对于第一种方法，只需要进行fix nvt temp 1000 100 …的设置就行了，降温的快慢，由你run N 的N来决定。
对于第二种方法提供如下思路，也就是使用循环的方法（升温，降温改小参数，+改成-就行啦）：
1.先在最低温度下，跑平衡输出restart.equ低温稳定的结构
units metal
boundary p p p
atom_style atomic
neighbor 0.3 bin
neigh_modify delay 10
lattice bcc 2.87
region box block 0 10 0 10 0 10
create_box 1 box
create_atoms 1 box
pair_style eam/fs
pair_coeff * * fe.eam Fe
mass 1 56
fix 1 all npt temp 373.0 373.0 100.0 iso 1.0 1.0 1000.0thermo 10
timestep 0.001
run 1000000000
write_restart restart.equ
- 用循环的方法进行准平衡升温抽样
units metal
boundary p p p
atom_style atomic
label loopa
variable i loop 25
variable t equal 273.0+100.0*$i
variable v equal vol
read_restart restart.equ
pair_style eam/fs
pair_coeff * * fe.eam Fe
fix 1 all npt temp $t $t 100.0 iso 1.0 1.0 1000.0
thermo 10
timestep 0.001
run 10000000(足够长)
write_restart restart.equ
print “$t $v” append gmm.txt
clear
loop
next i
jump in.script loopa
例子只提供思路，不保证每个命令的语法正确