大家好，我是小马老师。
最近几天有几个粉丝咨询如何把随机选择的辐照PKA原子设置到模型的中心位置，本文介绍一下这个方法。
这种建模，有两个关键点：随机选择和原子居中。
(1)随机选择
随机选择的方法在之前的推文中也介绍过，使用的代码为：
variable tmp equal floor(random(1,686,9090))variable atomid equal ${tmp}print "id=${atomid}"group pka id ${atomid}set group pka type 2
686为data文件或者模型中的最大原子数，可根据自己的模型修改
随机选择pka后的模型如下图所示：
（2）居中显示
居中可使用fix recenter命令：
fix 1 pka recenter 0 0 0后面的参数0 0 0为模型的中心坐标，这个也需要根据自己的模型进行修改
fix recenter命令可以把设置的原子组pka调整到后面三个坐标点为止，同时，整个模型也会自动地进行偏移。
调整后的模型如下图：
蓝色的pka原子已经位于模型的中心位置。
报名一对一培训课程可领完整版《lammps中文教程》，配套学习更轻松！
扫码领取粉丝专属优惠券↓↓↓
