在lammps中加恒定剪切力的命令

在 LAMMPS 中，可以使用 fix deform 命令来施加剪切力场。具体步骤如下：

1. 定义一个初始的盒子大小和形状，使用命令：

lattice <style> <a1> <a2> <a3> <xhi> <yhi> <zhi> [align <keyword> <arg>]* [origin <x> <y> <z>]

其中，<xhi>，<yhi>，<zhi> 分别表示盒子的边界。在这个命令中，我们可以设定 x 轴和 y 轴的长度，而 z 轴的长度可以在后面的 fix deform 命令中设定。

2. 使用 create_atoms 命令创建原子，例如：

create_atoms 1 random 100 12345 NULL

其中，1 表示原子类型，random 表示随机分布，100 表示原子数目，12345 表示随机数种子，NULL 表示不使用 Velocities。

3. 使用 pair_style 和 pair_coeff 命令来定义原子之间的相互作用力。

4. 使用 fix npt 命令来模拟恒定温度和压力的正则化运动。

5. 使用 fix deform 命令来施加剪切力场。例如：

fix 1 all deform 1 xy erate 0.001

其中，1 表示 fix 命令的 ID，all 表示对所有原子施加这个力场，deform 表示变形模式，1 表示变形的步骤数，xy 表示 x 轴和 y 轴方向施加剪切力，erate 0.001 表示剪切速率。

注意，fix deform 命令只能在 NPT 模拟中使用。