ReaxFF-弹性墙壁：添加外电场、弹性墙壁
LAMMPS中，外加电场时，只考虑电场对原子的“受力”的影响，不考虑对能量（进而对化学反应）的影响。在AMS中，电场的影响不仅仅体现在力上，也体现在能量上。
模拟体系存在外加电场的情况下，如果有原子穿越Cell的边界，就会出现EEM电荷（EEM：ReaxFF计算原子电荷的方法模型）unphysically large相关报错。为了避免穿越Cell边界，因此需要设置一个弹性墙壁，限制原子的活动范围。当然这就变成了一个非周期性的体系，为了避免分子扩散，同理也是需要这个弹性墙壁的。
下面以一个简单的例子，演示如何添加外电场，以及如何添加弹性墙壁。思路如下：首先创建一个常规的盒型体系，去掉周期性（Main中周期性设置为None），然后将其修改为球形，设置外加电场后，设置球形弹性墙壁，球半径略大于模型半径。
本教程适用于AMS2022及其以上版本，2019及以前的版本，请参考：ReaxFF如何外加电场
建模
这里以均相体系H2O和O2混合物为例，由于弹性墙壁是球形（而非旧版的正交盒形），因此我们在创建均相体系之后，可以把模型修改为球形（这样能确保分子扩散后，密度尽量不发生大的变化）。当然非均相体系也是可以的，只是需要注意球半径的设置。
创建均相混合物
首先我们创建H2O和O2的混合物，AMSinput → Edit → Builder：
注意设置好分子种类与个数之后，窗口底部会显示密度，如果没有问题，点击Generate即可生成该体系。
将体系修改为球形
Main → Periodicity设置为None。选中尽量接近坐标原点的某个原子，然后Select → Select within Radius，输入一个半径例如12 Å（建议这个值略小于Cell直径的一半，确保球在Cell内），OK。这样一个半径为12 Å的一个球被选中了，此时球表面有的分子只有部分原子被选中，我们希望完整的分子都被选中，因此我们按Ctrl m键，则相关完整分子都被选中了（也意味着半径略有扩大）。
Select - Invert Selection反选，然后del键删除，这样留下来的，就是刚才半径约为12 Å的一个球：
那么球形弹性墙壁，大致可以设置为12 Å多一些（如果弹性墙壁半径没有包纳全部原子，运行将报错Atoms too close）。本例中，设为17埃才正常运行。
当然如果不介意密度的问题，则模型根本不需要修改为球形，直接用一个大的球，将整个Cell中的原子全部包含进去即可。
设置外加电场
Model → Electric Field，三个数值是电场沿x、y、z方向的分量，例如这里我们设置z方向电场强度为0.01V/Å（注意V/Å是一个非常大的单位1Å的距离将会产生1V的电压差）
该设置，在生成的*.run文件中，体现在System关键词中
System
......省略
ElectrostaticEmbedding
ElectricField 0.0 0.0 0.01
End
......省略
End
设置球形弹性墙壁
照常设置基本分子动力学参数，例如力场、系综、温度、步数、保存轨迹的频率等等，可以参考其他教程。设置完毕后保存作业，在生成的*.run文件中，我们添加如下字段：
EngineAddons
WallPotential
Enabled Yes
Gradient 5.0
Prefactor 1.0
Radius 17
End
End
这里Gradient是描述弹性墙壁的Potential强度变化梯度，单位是1/Å，这个值设置的越小，墙壁越陡峭，可以理解为墙壁的厚度为这个值。因此一般建议这个值设置的很小例如5。这里Prefactor是Potential的指前因子，Radius是球半径，单位为Å，这里我们设置为17 Å，足以将之前的原子全部包围。
这个字段可以添加在*.run文件任何字段之间，例如在MolecularDynamics字段后面（该例子是一个NVT系综）：
MolecularDynamics
NSteps 50000
InitialVelocities
Temperature 298
End
Thermostat
Type NHC
Temperature 298
Tau 5
End
End
EngineAddons
WallPotential
Enabled Yes
Gradient 5.0
Prefactor 1.0
Radius 17
End
End
System
ElectrostaticEmbedding
ElectricField 0.0 0.0 0.01
End
......省略
End
保存*.run文件，并在AMSJobs窗口运行该作业（或在命令行运行该*.run文件）。
结果检验
SCM → Movie：
可以看到经过2 ps以后，分子仍然被限定在一个非常狭小的范围，没有扩散到真空中。