模拟降温,有大致两种办法,一个是连续(升)降温,也就是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
- 用循环的方法进行准平衡升温抽样