lammps教程：原子组整体漂移怎么办？

大家好，我是小马老师。
在lammps模拟过程中，大家应该会遇到过原子组整体移动的情况，有时候我们不希望原子整体偏移，这时候可以用“fix recenter”命令。
fix recenter语法为：

fix ID group-ID recenter x y z

该命令的作用是使原子组重心固定在坐标点(x,y,z)。
xyz可以设置坐标值，也可以使用固定的参数，如INIT表示原子组初始位置，NULL表示在此方向不固定原子组。
下面以两个Cu球的烧结为例，演示fix recenter命令的使用方法。
左侧的Cu使用如下命令，保持重心不变：

fix 2 left recenter INIT INIT INIT

右侧Cu球不加任何限制。
模拟结果为：

从上图可以看出，左侧Cu球位置基本没变，右侧Cu球移动距离较大。
in文件代码：

boundary    p p p
units       metal
atom_style  atomic
neighbor    2 bin
neigh_modify    every 1 check yes
timestep    0.001
#创建box
region      box block -50 50 -20 20 -20 20 units box
create_box  1 box
#建模
lattice     fcc 3.61
region      ball1 sphere -12.5 0 0 10 units box
create_atoms 1 region ball1
region      ball2 sphere 12.5 0 0 10 units box
create_atoms 1 region ball2
group       left region ball1
#力场设置
pair_style  eam/fs
pair_coeff  * * Cu1.eam.fs Cu
#温度初始化
velocity    all create 300 8999
#升温
dump        1 all atom 100 dump.xyz
fix         1 all nvt temp 300 2000 0.1
fix         2 left recenter INIT INIT INIT
run         20000

更多lammps教程，请关注公众号：lammps加油站