lammps常见设置误区：“刚体”和“固定”命令分不清

在lammps模拟中，有“刚体”和“固定”的设置，在平时咨询过程中，发现还有小伙伴对两者的区别分不大清楚。

在设置”刚体“的时候使用固定的命令，反之也有。

”刚体“和”固定“的设置有相同点，也有不同点。

（1）相同点

被设定为”刚体“以及”固定“的结构，原子之间的间距保存相对固定，在模拟过程中，这个结构不会发生变形。

换一句描述，可以理解为：”刚体“是一个刚体，被“固定”的结构也是一个刚体，这是两者之间相同点。

（2）不同点

第一个不同：“刚体”可以移动，被“固定”的结构静止不动。

比如，在摩擦模拟中，磨粒可以设为一个“刚体”，在摩擦时磨粒自身不发生变形，但是磨粒可以移动以产生磨削效果。

水分子可以设置为一个“刚体”，水分子的键角在模拟时不再发生变化，便于简化运算，加快计算速度。

第二个不同：“刚体"可以具有一定的温度，被”固定“的原子温度为0K。

可以对刚体设置系综，进行控温控压，被“固定”的原子不能具有温度，否则原子就会有移动速度，产生移动。

（3）设置方法

“刚体”设置：

刚体设置主要使用fix rigid、fix rigid/nve、fix rigid/nvt、fix rigid/npt等命令，另外有些命令自带刚体效果，如fix move命令。

“固定”设置：

velocity fixed set 0 0 0

fix 1 fixed setforce 0 0 0

“固定”原子的基本原理是原子的速度设置为0，同时，原子的受力设置为0，两个条件必须同时具备，缺一不可。

如原子的初速度不为0，仅受力被设置为0，原子会沿初始速度进行匀速运动。如果初速度为0，但受力不为0，原子会产生新的速度，也会产生移动。

在模拟时，如果发现被“固定”的原子产生了移动，需要考虑以上两个条件是否同时满足。

以上是对“刚体”和“固定”的简单对比介绍，希望对你有用。

---

公众号：lammps加油站