模擬降溫，有大致兩種辦法，一個是連續（升）降溫，也就是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

例子只提供思路，不保證每個命令的語法正確