lammps固定原子方法2

大家好，我是小马老师。

本文介绍lammps固定原子的另外一种方法：fix spring/self。

在lammps模拟中，固定原子组最常用的方法是把原子组的速度和受力同时设为0。

使用这种方法，原子组在模拟过程中是完全静止的，并且因为速度为0，所以原子组的温度也为0。

但在有些模拟中，既想固定原子组，又想原子组带有一定的温度。

显然，用传统的固定方式是无法实现的，可以考虑fix spring/self命令。

（1）原理

原子组中的每一个原子都“连接”一个弹簧，原子在力的作用下偏离初始位置时，弹簧被拉长，弹簧产生的拉力把原子拉回到初始位置，所以，原子组被“固定”在初始位置。

（2）语法规则

fix ID group-ID spring/self K dir

K：弹簧系数，弹簧力F