使用动力学蒙特卡罗方法对环四亚甲基四硝胺中热点点火的中尺度模拟

《The Journal of Physical Chemistry A》：Mesoscale Simulations of Hot Spot Ignition in Cyclotetramethylene Tetranitramine Using the Kinetic Monte Carlo Method

　　

使用动力学蒙特卡罗（kMC）方法，计算了二次炸药环四亚甲基四硝胺（HMX）中直径达0.5 μm的圆柱形热点点火临界温度。通过将热比特传输模型与热爆炸起始的有效单步Arrhenius动力学相结合，研究人员能够在介观尺度上观察HMX中持续爆燃反应的均匀成核和生长。kMC模拟捕捉到了量子分子动力学（QMD）模拟在纳米长度尺度上观察到的热爆炸随机起始，该长度尺度需要对爆燃前沿的温度分布进行充分解析。研究人员针对HMX的四组单步Arrhenius动力学评估了爆燃速度和热点临界温度，这些动力学以不同方式考虑了自加热对反应速率加速的影响。结果表明，Manner等人、Burnham和Weese以及Henson等人提出的单步Arrhenius动力学给出了大致相似的热点临界温度，但重新参数化2001年Henson-Smilowitz动力学以包含自加热对反应速率的影响，导致显著不同的行为，热点临界温度比其他模型高出约700 K。讨论了这种行为的原因。

**论文解读：使用动力学蒙特卡罗方法模拟环四亚甲基四硝胺中热点点火的中尺度行为**

**研究背景与动机**

在弱冲击压缩下，凝聚相含能材料内部微小孔隙的塌缩会导致局部高温区域的产生，这些区域即“热点”（hot spot），可能引发化学反应。冲击到爆轰的转变高度依赖于冲击波前沿与塌缩孔隙内化学反应之间的相互作用。爆燃前沿从塌缩孔隙中萌生，通过热传导以极慢的速度（米/秒量级）传播，无法追上超声速的主冲击波（千米/秒量级），但初始放热、体积增大的反应所发射的瞬态压力波能够到达并强化主冲击波。因此，准确预测热点处放热反应的点火与增长，对于理解含能材料的性能与安全性至关重要。

目前存在的主要问题包括：传统基于网格的连续介质模拟需要纳米级分辨率才能准确捕捉爆燃前沿的温度梯度，这直接与分子动力学模拟的尺度重叠；然而，标准网格算法忽略了纳米尺度上热爆炸起始的随机性，而这种随机性在量子分子动力学（QMD）模拟中已被观察到。因此，需要一种能够桥接原子尺度（纳米）与介观尺度（微米），并同时考虑热传输与随机热激活反应的方法。

**研究内容与结论**

本研究应用动力学蒙特卡罗（kMC）方法，模拟了二次炸药环四亚甲基四硝胺（HMX）中直径小于0.5 μm的圆柱形热点的点火与爆燃传播。通过将热比特传输模型与有效单步Arrhenius动力学相结合，研究人员成功在介观尺度上观察到了HMX中持续爆燃反应的均匀成核与生长。kMC模拟捕捉到了与QMD模拟一致的纳米尺度热爆炸随机起始。研究评估了四组不同单步Arrhenius动力学参数（MC24、Burnham–Weese、2001 Henson–Smilowitz、Henson–Smilowitz–FK）下的爆燃速度与热点临界温度。结果表明，前三种动力学模型给出的热点临界温度大致相似，但重新参数化后的Henson–Smilowitz–FK动力学因考虑了自加热效应，导致热点临界温度比其他模型高出约700 K，且爆燃速度显著降低。该论文发表在《The Journal of Physical Chemistry A》。

**关键技术方法**

研究人员主要采用动力学蒙特卡罗（kMC）方法，结合热比特传输模型模拟热传输，使用有效单步Arrhenius动力学描述热爆炸起始，并考虑自加热效应。kMC中的反应速率和热传输速率通过随机采样进行更新。热扩散率α(T)基于温度依赖的热导率λ(T)、密度ρ(T)和热容C_P(T)参数化，这些参数来源于文献中的实验数据和分子动力学模拟结果。四组Arrhenius动力学参数分别来自Manner等人（MC24，基于QMD模拟和ODTX实验）、Burnham和Weese（基于近等温分解实验）、2001年Henson–Smilowitz（基于多种时间爆炸数据）以及重新参数化的Henson–Smilowitz–FK（使用Frank–Kamenetskii绝热诱导时间拟合）。所有模拟均在HMX的β多晶型基础上进行，未考虑多态转变或熔化。

**研究结果**

**3.1 网格收敛性与爆燃前沿结构**
通过一维kMC模拟，研究人员计算了爆燃速度v_d对网格间距h的依赖关系。结果显示，对于所有四组动力学参数，v_d随h减小而急剧增加，说明算法收敛缓慢但存在收敛性。MC24、Burnham–Weese和2001 Henson–Smilowitz动力学在h=1 nm时给出v_d≈10–30 m/s，而Henson–Smilowitz–FK动力学给出的v_d≈3 m/s，显著更慢。爆燃前沿宽度与速度相关：速度最快的MC24动力学前沿最窄（约50 nm），速度最慢的Henson–Smilowitz–FK动力学前沿最宽（约250 nm）。反应区温度在10–20 nm内快速上升，且前沿宽度随网格分辨率提高而变窄。

**3.2 热点临界性在HMX中的应用**
研究人员在二维750×750网格（h=2 nm）上模拟了圆柱形热点（直径d_H≤0.5 μm），背景温度T₀=293 K。通过改变热点初始温度T_H，确定每个d_H下使热点区域至少50%反应的最小T_H，即临界温度T_c。结果表明，T_c在d_H<0.2 μm时快速上升，与Tarver等人的结果定性一致。MC24和Burnham–Weese动力学给出的T_c在整个d_H范围内相差不超过70 K，在d_H=0.5 μm时约为1160 K，比Tarver估计值高出约100 K。2001 Henson–Smilowitz动力学的T_c高100–200 K，而Henson–Smilowitz–FK动力学则高出超过700 K。此外，对于d_H=0.05 μm的热点，MC24、Burnham–Weese和2001 Henson–Smilowitz动力学均观察到点火后熄灭（耗散而无增长），而d_H>0.05 μm的热点则能点火并增长，但前沿速度在进入冷材料后显著降低。Henson–Smilowitz–FK动力学未观察到明确的爆燃前沿，而是呈现均匀热爆炸。

**讨论与结论**

研究讨论了不同动力学模型对热点临界温度的影响。MC24和Burnham–Weese动力学在高温区（T>1000 K）给出相似的诱导时间和T_c。2001 Henson–Smilowitz动力学虽未显式考虑自加热，但可能因误差抵消而与前两者定性一致。而Henson–Smilowitz–FK动力学由于自加热效应的重新参数化，导致反应显著变慢，T_c异常升高。研究指出，小热点的点火与熄灭也可能通过压力波贡献于冲击起爆过程。结论部分翻译如下：研究人员已证明，通过将热传输与kMC框架内的随机热爆炸相耦合，能够有效桥接原子尺度与介观尺度。在微米长度和纳秒时间尺度上对热爆炸和爆燃进行模拟，这是大规模反应性MD模拟难以实现的。kMC形式主义捕捉了随机的、热激活的态-态转变，对应于惰性诱导时间后亚纳米立方体体积材料的热爆炸起始。此外，kMC框架通过依赖于温度的热扩散率α（由参数化热导率λ、密度ρ和热容C_P描述）纳入热传输。经典MD模拟能准确捕捉传导率和密度的温度依赖性，但其热容在所有温度下均等于经典极限。最后，kMC方案提供了一种直接的方法来探索反应速率或能量输出对含能材料对热刺激敏感性的影响。