在本ReaxFF教程中，解释了从聚合物密度的温度依赖性计算玻璃化转变温度T_g。 详细系统性研究，参考文献：

• Effect of chemical structure on thermo-mechanical properties of epoxy polymers: Comparison of accelerated ReaxFF simulations and experiments, Polymer 159, 354-368 (2018).

本教程包含如下几个步骤：

1. 导入聚合物结构
2. 模拟退火
3. 提取密度-温度性质
4. 计算玻璃化转变温度

本计算对硬件要求较高，建议在集群或工作站上运行。

这里我们使用bond boost法生成的交联环氧聚合物，当然如果有其他结构的*.xyz文件，也可以直接通过AMSinput → File → Import Coordinates导入。

为了精确模拟密度的温度依赖性，我们必须确保聚合物结构尽可能均匀。聚合物结构内部的大孔是不均匀结构的指标，这将导致非常嘈杂和难以分析的结果。模拟退火是一种简单的优化策略，旨在在复杂的高维势能面上产生全局最小能量结构。在退火过程中，MD模拟过程中的温度逐渐升高，希望系统能够克服越来越多的能垒，有效探索能量势能面全局。加热阶段之后是冷却阶段，在此期间，系统被冷却回初始温度。通常，结果结构的能量明显低于起始结构，尽管不能保证它是最佳结构。

参数设置如下：

退火过程中，体系的体积有可能扩大或者缩小，因此使用恒压器（此处使用恒压器还有一个原因是原子的个数足够多，如果体系只有几百元子、一两千原子，恒压器效果就不那么好，原子个数越多，恒压器控压的效果越好）

接下来，我们定义所需的温度程序。为了从计算中提取密度，我们将把温度从298.15提高到598.15，这对于模拟退火来说仍然很低，但对于我们的目的来说已经足够了。我们要实施的温度程序很简单：花费30000步升温25K，随后10000步恒温对密度进行取样，循环此步骤直到达到598.15 K，然后反向冷却至298.15 K。如下图所示意（下图仅仅示意了两步升温、采样，实际需要非常多步升温、采样，详见下文中温度的设置即可知）：

在“Thermostat”中设置温度：

298.15 298.15 323.15 323.15 348.15 348.15 373.15 373.15 398.15 398.15 423.15 423.15 448.15 448.15 473.15 473.15 498.15 498.15 523.15 523.15 548.15 548.15 573.15 573.15 598.15 598.15 573.15 573.15 548.15 548.15 523.15 523.15 498.15 498.15 473.15 473.15 448.15 448.15 423.15 423.15 398.15 398.15 373.15 373.15 348.15 348.15 323.15 323.15 298.15 298.15

“Duration”中设置每个阶段的步数：

10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000 30000 10000

如下图所示：

保存并运行作业，建议使用8~16核。