基于AIMD的均四嗪小分子热分解模拟
一、研究背景
均四嗪（Tetrazine）小分子因其在能源储存、推进剂和爆炸物等领域的应用而受到广泛关注。然而，均四嗪在高温下的热分解机理至今尚不完全清楚。为了深入理解这一过程，本研究使用基于从头算的分子动力学（AIMD）模拟，探讨了均四嗪小分子的热分解反应机理。通过AIMD方法，能够精确模拟均四嗪分子在高温下的反应路径，揭示其热分解过程中的关键中间体和过渡态。
本案例基于xtb软件，采用GFNn-xTB哈密顿量进行AIMD模拟。xtb是一款高效的半经验量子化学程序，在处理中等规模分子体系时兼具精度与速度优势，适合用于模拟键的形成与断裂过程。
二、初始模型构建
首先利用GaussView构建均四嗪分子的初始结构，其分子结构如图1所示。均四嗪为含四个氮原子的六元杂环化合物，具有独特的高氮含量结构。
图1 均四嗪小分子的结构
相关文献表明，均四嗪小分子在高温热解后会生成一个N₂分子和两个HCN分子，其热解产物示意如图2所示。这一反应路径为后续模拟分析提供了重要的参考基准。
图2 均四嗪的热解产物
三、模拟设置
AIMD模拟采用xtb程序执行，具体参数设置如下：模拟温度设置为3500K，以加速热分解过程并使反应在可接受的时间尺度内完成；模拟总时长设置为5ps，时间步长为1fs。1fs为该类模拟体系可接受的最大步长，在涉及键的形成与断裂过程中不宜对键长施加约束，因此shake参数设置为0，即不进行键长约束。
图3 xtb的输入文件
xtb的输入文件配置了温度、热浴、步长等关键参数，确保模拟在给定条件下正确执行。AIMD模拟中，高温设置是为了克服反应能垒，使热分解过程在有限模拟时间内得以观测。
四、模拟结果与分析
经过5ps的AIMD模拟后，可以清晰观察到均四嗪的结构演变及热分解全过程，其时间序列快照如图4所示。
图4 均四嗪的热分解过程
分析模拟轨迹可知：t=1.5ps时，均四嗪分子开始发生显著的扭曲变形，环状结构偏离平面构型；t=2.7ps时，四氮基环中一对相邻氮原子间的键发生断裂，引发环结构的局部破裂，进而触发了热分解连锁反应；t=3.8ps时，反应基本完成，均四嗪成功分解为目标产物，即一个N₂分子和两个HCN分子，与文献预测结果完全一致。
模拟结果充分验证了均四嗪高温热解的反应路径，揭示了环断裂这一关键速率决定步骤，为理解高能材料的热稳定性提供了理论依据。
五、结语
本案例基于xtb软件采用AIMD方法成功实现了均四嗪的热分解过程模拟，完整揭示了从初始结构到N₂和HCN产物的微观反应路径。研究成果不仅为均四嗪分子的热稳定性提供了深入的理论支持，也为能源材料设计与安全性评估提供了重要的参考数据。
