LAMMPS降(升)温的几个思路

模拟降温，有大致两种办法，一个是连续（升）降温，也就是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

2. 用循环的方法进行准平衡升温抽样

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

例子只提供思路，不保证每个命令的语法正确