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

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

你可能感兴趣的:(模拟技术,LAMMPS)