应用案例：

• 铂使用巨正则系综反应力场模拟铂电催化剂的氧化与失活（Small 2019, Angew. 2017）

本教程针对AMS2019.3及之前的版本，AMS2020以后的版本请参考链接：ReaxFF-GCMC：巨正则系综蒙特卡洛模拟

这里以Si表面的氧化为例。

View - Periodic - Show Unit Cell:

点击底部四边形按钮，切出表面：

Edit - Crystal - Generate Super Cell，生成超胞

关于‘μ (kcal/mol)’这一项，上面的例子中，化学势设置为-75kcal/mol。化学势计算公式： [μ(T,P)+k_bT*ln(P/P')-E_d]/2

• μ(T,P)可以通过计算该分子的气相Gibbs自由能得到，参考如何计算分子的频率、红外、零点能、转动能量、转动惯量、熵、焓、热熔与Gibbs自由能。计算的时候，所取温度、压强与系综温度、压强一致即可，如果有实验测得的化学势更好
• k_b为波尔兹曼常数，T为温度
• P'是小分子的分压强。如果只有一种气体，那么k_bT*ln(P/P')＝k_bT*ln1＝0
• E_d是分子的解离能，也就是ADF计算得到的分子的Bonding Energy
• 目前ReaxFF只支持原子的GCMC，不支持分子……

在movie中可以看到Si的氧化结果：

初始帧即为GCMC最后一步的结果。第2帧开始显示逐渐氧化的过程，最后一帧与第一帧相同。