Fe原子辐照轰击多层石墨烯模拟代码

本文介绍一个Fe原子辐照多层石墨烯的案例。

辐照模拟是指一个或多个原子以较高的速度与另一物体发生碰撞，造成被撞击物体损伤。

因为原子速度较高，因此在模拟时需要合理地设置模拟步长，并根据不同体系在模拟过程中随时调节步长。

主要模拟过程有以下几步：

1. 建模

为防止丢失原子错误，模型采用较大的一个box，边界条件为固定边界。

使用lattice custom命令在box中间创建多层石墨烯。

使用create_atoms命令生成一个Fe原子。

2. 力场设置

本案例力场为混合力场，石墨烯为lammps自带的airebo力场，因为只有一个Fe原子，因此Fe原子力场为LJ力场，参数设为0，C-Fe原子为LJ力场。

3. 辐照过程

选定Fe原子，设定Fe原子能量为10Kev，根据能量公式把能量单位转换为速度单位A/ps

lammps in文件代码如下：

units                metal

dimension            3

boundary            f f f

neighbor            0.3 bin

neigh_modify        delay