《FirePhysChem》：Computational Design and Screening of Fused [5,6,5]-Tricyclic Energetic Materials via Tunable Amino/Nitro Substitution
编辑推荐：
本研究引入了两种新型稠合[5,6,5]-三环骨架结构（吡唑-三嗪-三唑和吡唑-三嗪-四唑），该结构协同结合了吡唑/三嗪的稳定性与三唑/四唑的高能量特性。研究人员通过策略性地引入混合氨基（-NH?）与硝基（-NO?）基团，成功实现了接近黑索金（RDX）水平的优异
本研究引入了两种新型稠合[5,6,5]-三环骨架结构（吡唑-三嗪-三唑和吡唑-三嗪-四唑），该结构协同结合了吡唑/三嗪的稳定性与三唑/四唑的高能量特性。研究人员通过策略性地引入混合氨基（-NH?）与硝基（-NO?）基团，成功实现了接近黑索金（RDX）水平的优异爆轰性能，同时保持了优越的安全参数。完全硝化衍生物表现出卓越的爆轰速度（>9.40 km/s）和爆轰压力（>40.0 GPa），优于三硝基甲苯（TNT）和RDX等传统炸药。研究人员利用一套稳健的参数体系（静电表面势ESP、键解离能BDE、自由空间体积及芳香性指数），系统地阐明了稳定性和感度的电子与结构起源。研究将稠合吡唑-三嗪-四唑骨架确定为下一代高能量密度材料（HEDMs）的特权分子架构。
研究背景与科学问题
含能材料是一类用于火箭推进、军事和民用领域的关键物质，其核心挑战在于如何平衡高能量密度与分子稳定性。传统单环含能化合物如TNT和RDX虽已广泛应用，但在能量输出与安全性之间仍存在固有矛盾。近年来，稠环化合物因其平面共轭骨架能够提供高生成焓、高密度和优异的爆轰性能而受到广泛关注。然而，已报道的稠合[5,6,5]-三环含能材料在结构多样性上仍有局限，难以同时实现高能量与低感度的理想平衡。为解决这一问题，研究人员设计了两种新型稠合[5,6,5]-三环骨架——吡唑-三嗪-三唑（A系列）和吡唑-三嗪-四唑（B系列），旨在通过融合不同杂环的优势，并结合氨基与硝基的可调取代模式，探索高性能与高稳定性兼得的含能分子。
研究内容与意义
研究人员采用高水平量子化学计算方法，对上述两个系列的六种化合物（A2至A4、B2至B4）进行了系统的理论评估。研究发现，完全硝化的衍生物A3和B3展现出超越TNT和RDX的爆轰性能，而混合氨基/硝基取代的A4和B4则在能量与感度之间取得了最佳平衡。该工作不仅提出了新型含能分子骨架设计策略，还通过多维度参数（静电表面势、键解离能、自由空间体积、芳香性指数）建立了结构与性能之间的关联规律，为下一代高能量密度材料（HEDMs）的理性设计提供了重要指导。
关键技术方法概述
所有化合物的几何优化和电子能量计算均采用密度泛函理论（DFT），在B3LYP/6-311+G(2d,p)水平下进行，并通过振动频率分析确认无虚频。气相生成焓（HOF）的计算采用高斯-4（G4）组合技术，固相生成焓则通过减去升华焓获得。晶体密度预测基于B3LYP/6-311++G(d,p)水平的分子表面性质计算。爆轰速度和爆轰压力使用EXPLO5（v7.01.01）程序计算。感度参数包括晶胞自由空间（ΔV）和撞击感度（h??）在B3PW91/6-31G(d,p)水平下评估。键解离能（BDE）在B3LYP/cc-pVTZ水平下计算。静电表面势（ESP）、相互作用区域指示符（IRI）以及芳香性指数（多中心键序MCBO和谐振子芳香性度量HOMA）均使用Multiwfn 3.8程序处理。
研究结果分述
3.1 生成焓（HOF）
生成焓是衡量分子内储存化学能的重要热力学参数。采用G4理论计算的母体骨架A1和B1的气相HOF分别为701.23 kJ/mol和818.69 kJ/mol，B1因含有更高氮含量的四唑环而比A1高出约117 kJ/mol。所有设计化合物的固相HOF均为正值，范围在364.36至724.72 kJ/mol之间。B系列化合物的HOF普遍高于对应的A系列，这主要源于骨架本身的能量贡献差异。在各自系列中，完全氨基化化合物（A2、B2）的HOF最低，完全硝化化合物（A3、B3）的HOF最高。
3.2 物理化学与爆轰性能
氧平衡（OB）、氮含量和密度是决定爆炸能量输出的关键参数。完全硝化的A3和B3的氧平衡分别为-22%和-25%，与RDX相当，远优于TNT（-74%）。它们的密度分别达到1.98 g/cm3和1.97 g/cm3，显著高于TNT（1.65 g/cm3）和RDX（1.80 g/cm3）。混合氨基/硝基取代的A4和B4则表现出中间状态，密度分别为1.82 g/cm3和1.86 g/cm3。爆轰性能方面，完全硝化的A3和B3表现最为突出，爆轰速度分别为9.44 km/s和9.60 km/s，爆轰压力分别为41.83 GPa和40.63 GPa，全面超越TNT（6.94 km/s, 22.0 GPa）和RDX（8.60 km/s, 33.92 GPa）。混合取代的A4和B4的爆轰速度和压力分别为8.13 km/s、24.19 GPa和8.57 km/s、30.96 GPa，介于完全氨基化和完全硝化化合物之间。
3.3 感度与稳定性趋势
感度参数分析显示，完全氨基化的A2和B2具有最小的晶胞自由空间（ΔV分别为43 ?3和38 ?3）和最高的撞击感度值（h??分别为56 cm和55 cm），表明其机械感度最低，安全性最好。完全硝化的A3和B3则相反，ΔV较大（54 ?3和45 ?3），h??较低（27 cm和30 cm），感度较高。平衡参数（ν）的分析进一步印证了这一趋势：A2（ν=0.246）和B2（ν=0.248）最接近理想值0.25，代表最均衡的表面静电势分布和最低感度；而A3（ν=0.112）和B3（ν=0.124）的ν值最小，感度最高。综合感度排序为：A2 ≈ B2（低感度）< A4 < B4 < B3 < A3（高感度）。键解离能（BDE）分析表明，C-NO?键的BDE值均在236-249 kJ/mol范围内，远高于120 kJ/mol的热稳定性阈值，证实了这些化合物的良好热稳定性。相比之下，N-NO?键的BDE较低（100-104 kJ/mol），被认为是分解引发时的潜在触发键。
3.4 静电表面势（ESP）与相互作用区域指示符（IRI）
ESP分析揭示了分子的电荷分布特征。红色区域（正电势）主要集中在中央三嗪环附近，蓝色区域（负电势）分布在三唑/四唑环的氮原子和硝基上。完全硝化的A3和B3的总正ESP值最高（分别为456.45和400.66 kcal/mol），与其较高的感度一致。IRI分析展示了分子间氢键和范德华相互作用的可视化图像，其中B4二聚体显示出更大的连续氢键区域和更均匀分布的范德华相互作用，暗示更强的分子间协同作用网络。
3.5 芳香性指数
多中心键序（MCBO）和谐振子芳香性度量（HOMA）分析表明，B系列的四唑环比A系列的三唑环具有更高的芳香性指数。氨基和硝基取代对芳香性影响较小，说明芳香性主要由稠合杂环骨架本身决定，而非末端取代基。
3.6 与已报道稠合[5,6,5]-三环化合物的性能比较
与文献中12种代表性稠合[5,6,5]-三环含能化合物相比，设计的A3和B3在密度（1.99和1.98 g/cm3）、爆轰速度（9.48和9.64 km/s）和爆轰压力（39.82和41.21 GPa）方面均全面超越所有参比化合物。混合取代的A4和B4也展现出有竞争力的性能，且其C-NO?键BDE值（249和244 kJ/mol）与多数参比化合物相当或更优。
结论总结
本研究设计并系统评估了两个系列的稠合[5,6,5]-三环含能化合物：吡唑-三嗪-三唑（A系列）和吡唑-三嗪-四唑（B系列），并对其进行了氨基和硝基官能化修饰。通过高水平量子化学计算，研究人员深入研究了其分子结构、芳香性、热力学稳定性、感度和爆轰性能。计算结果表明，硝基取代对于获得高爆轰性能至关重要。完全硝化的衍生物A3和B3的氧平衡与RDX相当，密度超过TNT和RDX，爆轰速度和压力均优于这些传统炸药。混合氨基/硝基衍生物A4和B4则是在两个极端之间的折中选择：它们比完全氨基化衍生物具有更好的氧平衡和爆轰性能，同时相比完全硝化化合物保留了更有利的感度特性。MCBO和HOMA芳香性研究表明，稠合杂环骨架对电子离域有显著影响，而氨基和硝基取代基对内在芳香性的影响很小。通过平衡参数、晶格自由空间和撞击感度确定的静电势分析和感度趋势一致表明，完全氨基化衍生物（A2、B2）具有最均衡的电荷分布、最小的晶格自由空间和最低的预测感度，而最敏感的衍生物是完全硝化衍生物（A3、B3）。总之，本研究确定了稠合[5,6,5]-三环骨架，尤其是含四唑的B系列，作为开发高能量密度材料（HEDMs）的特权分子结构。混合氨基/硝基官能化的化合物A4和B4被认为是最有前景的候选物，在高能量性能与最佳稳定性和感度之间展现出了令人瞩目的平衡。
