What Happens in the First Microseconds: Using Quantum Mechanics to Map How RDX and HMX Begin to Decompose

RDX and HMX are cyclic nitramine compounds. If you have encountered the term "Composition C-4" or "octogen," you have encountered these materials. They are the primary energetic ingredients in most modern military explosives and in high-performance solid rocket propellants. They are exceptionally good at what they do: releasing large amounts of chemical energy very quickly.

Understanding how they burn and decompose is not just academically interesting. It is the foundation for predicting burn rates, designing propellants with desired performance characteristics, modeling the behavior of munitions, and assessing safety under abnormal heating conditions. And for decades, the models used for this purpose were built on a significant gap in knowledge: nobody had mapped, from first principles, what actually happens at the molecular level when these compounds start to break down.

The problem with purely experimental approaches

Thermolysis experiments — heating a small sample and measuring what comes out — can tell you a great deal about the final decomposition products. You can measure that CH₂O and N₂O are the major products, that CO and NO₂ also appear, and how the relative amounts change with temperature and heating rate. These measurements are valuable and reproducible.

What experiments struggle to capture are the initial steps — the first bond that breaks, the first intermediate that forms, the pathway by which the large cyclic molecule begins to unravel. These intermediates are transient: they form and react away on timescales too fast for most detection methods. Yet they are precisely what controls the subsequent chemistry. If you get the first steps wrong, the rest of the mechanism will be wrong in ways that are difficult to diagnose from product measurements alone.

This is the gap that quantum mechanics calculations can fill.

What quantum mechanics calculations actually do

At their core, QM calculations solve (approximately) the Schrödinger equation for the electrons in a molecular system. This gives you the potential energy of the molecule as a function of its geometry — the "energy landscape" that the atoms move on.

To study a decomposition reaction, you map the energy landscape between the reactant configuration and the product configuration. The highest point on this path is the transition state — the energy barrier that must be overcome for the reaction to proceed. From the barrier height and the molecular properties at the transition state, you calculate the rate constant using transition state theory:

$k(T) = \kappa \frac{k_B T}{h} \frac{Q^\ddagger}{Q_R} \exp\!\left(-\frac{\Delta E^\ddagger}{RT}\right)$

where $\Delta E^\ddagger$ is the barrier height, $Q^\ddagger$ and $Q_R$ are partition functions of the transition state and reactant, and $\kappa$ is the tunneling correction (important for hydrogen-transfer reactions at lower temperatures).

The choice of method matters significantly. For the initial nitramine decomposition pathways, the G4(MP2) composite method was used as the high-accuracy benchmark — it gives enthalpies of formation accurate to within ~1 kcal/mol for most organic molecules, which is the level of accuracy needed for reliable rate constants. For the large number of reactions in the full mechanism where G4(MP2) would be prohibitively expensive, the M06-2X density functional was benchmarked against G4(MP2) and found to give accurate barrier heights. This allowed the detailed mechanism to be built at a practical computational cost.

What we found for RDX

RDX is a six-membered ring: three CH₂ groups alternating with three N–NO₂ groups (hence the "hexa" in its chemical name, hexahydro-1,3,5-trinitro-1,3,5-triazine). The central question in the early decomposition is which bond breaks first: the N–NO₂ bond (releasing NO₂), or does the ring itself open first?

The answer is not obvious from the structure alone, and prior literature was divided. The QM calculations mapped both pathways — N–NO₂ homolysis and the ring-opening route — and found that while N–NO₂ homolysis has the lower barrier, the ring-opening reactions play an important and previously underappreciated role in the early decomposition, particularly at the temperatures relevant to liquid-phase and melt-phase conditions.

The calculations also identified the early hydrogen-abstraction reactions that generate key intermediates. These intermediates feed into subsequent reactions that eventually produce the observed final products (CH₂O, N₂O, NO₂, HCN, CO, CO₂) — but understanding how they get there requires knowing the intermediate species, which only the QM pathway analysis reveals.

What we found for HMX

HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) is the eight-membered ring analog of RDX — four CH₂ groups alternating with four N–NO₂ groups. It is denser and more powerful than RDX and is used in advanced propellant formulations where maximum energy density matters.

For HMX, the QM investigation revealed something important that the product measurements had not made clear: autocatalytic decomposition via HONO addition. As NO₂ is released in the early decomposition, it can react with neighboring molecules in the liquid phase (the cage effect — molecules in liquid are surrounded by neighbors and cannot immediately escape) to form HONO. HONO then acts as an autocatalytic species, attacking other HMX molecules and dramatically accelerating the decomposition rate. Hydrogen abstraction by NO₂ was identified as another dominant pathway.

This autocatalytic mechanism explains several features of the experimental data that had been difficult to account for — including the characteristic acceleration in the decomposition rate as the reaction progresses, even at constant temperature.

Why this matters beyond the specific materials

The significance of this work is not only in the RDX and HMX mechanisms themselves. It is in demonstrating a methodology: use quantum mechanics to establish the molecular-level foundation of a reaction mechanism, validate that mechanism against carefully designed experiments, and build upward from there toward a predictive engineering model.

The alternative — fitting global empirical models to experimental data without molecular-level insight — produces models that work within the conditions they were fitted to, but have limited predictive validity outside that range. For propellant design, safety analysis, or the development of new energetic formulations, that limitation is a genuine engineering problem.

The QM-derived mechanisms from this work became the foundation for a comprehensive kinetic model, validated combustion simulations, and a new approach to calculating transport properties from intermolecular potentials — all of which are described in the subsequent work on HMX combustion modeling.

The full work is described in:

• Quantum mechanics investigation of initial reaction pathways and early ring-opening reactions in thermal decomposition of liquid-phase RDX, Combustion and Flame, 2017
• Identification of initial decomposition reactions in liquid-phase HMX using quantum mechanics calculations, Combustion and Flame, 2018