**Thermal conductivity tensor of β-HMX as a function of pressure
and temperature from equilibrium molecular dynamics simulations**

A Pereverzev and T Sewell, PROPELLANTS EXPLOSIVES PYROTECHNICS, 48 (2023).

DOI: 10.1002/prep.202200345

We apply the Green-Kubo (G-K) approach to obtain the thermal conductivity tensor of ss-1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (ss- HMX) as a function of pressure and temperature from equilibrium molecular dynamics (MD) simulations. Direct application of the G-K formula exhibits slow convergence of the integrated thermal conductivity values even for long (120 ns) accumulated simulation times. To partially mitigate this slow convergence, we apply a recently implemented numerical procedure that involves physically justified filtering of the MD-calculated G-K heat current and fitting the integrated timedependent thermal conductivity to a physically motivated double-exponential function. The thermal conductivity tensor was determined for pressures 0.1 MPa <= P <= 30 GPa and temperatures 300 K <= T <= 900 K. The thermal conductivity k alpha beta (T,P) increases with increasing pressure, by approximately an order of magnitude over the interval considered, and decreases with increasing temperature. The MD predictions are compared to experimental and other theoretically determined values for the thermal conductivity of ss-HMX. A simple, semi-empirical fitting form is proposed that captures the behavior of kaaoT; PTHORN over the pressure and temperature intervals studied.

Return to Publications page