本文介绍原子平动和振动的设置方法。
lammps命令中,fix move命令功能很强大,可以控制原子的移动和转动,也可以根据用户自定义的公式设置原子的运动方式。
用到的关键词为variable,设定一个随时间变化的variable变量,就可以控制原子按变量进行移动。
语法为:
fix move variable args = v_dx v_dy v_dz v_vx v_vy v_vz
后面6个参数中,前三个控制原子位移,后三个控制原子速度,如果不需要单独控制,则需要设置为NULL。
下面以金刚石球为例,设置球沿x方向平移,并且在平移过程中,沿z轴上下振动。
需要定义两个变量:x和z。
沿x方向移动变量比较简单,随步数增加线性增加:
variable x equal "2*step*dt"
沿z方向变量需要设置一个正弦或余弦函数,实现周期性运动:
variable A equal 20.0
variable T equal 5.0
variable omega equal 3.0