lammps教程：lj/cut力场使用方法

lammps模拟设置中，力场的设置非常重要，力场设置合理才能准确的描述原子之间的受力。

如果两种原子之间没有键连接，可以使用pair_style命令设置这两种原子的力场。

pair_style能够设置的力场类型比较多，比如LJ、MORSE、EMA、MEAM等。

本文主要介绍最常用的力场：lj/cut。

lj/cut力场公式比较简单：

在lammps使用lj/cut只需要设置3个参数：epsilon、sigma、cutoff。

lj/cut有两种写法，第一种是写到in文件中，另一种是写到data文件中，两种写法稍有不同。

1. lj/cut在in文件中使用方法

在in文件中使用lj/cut，需要明确的指定相互作用的原子类型，如：

pair_style lj/cut 10
pair_coeff * * 0.02 3.12
pair_coeff 1 3 0.03 3.22 8.5
上面两句指定全部原子之间的lj/cut力场，截断半径为10，epsilon=0.02、sigma=3.12。