American Chemical Society, Journal of Physical Chemistry C, 5(119), p. 2290-2296, 2015
DOI: 10.1021/jp510328d
Full text: Download
We investigated the initial chemical reactions of BCHMX [cis-1,3,4,6-tetranitrooctahydroimidazo-[4,5-d]imidazole] with the following procedure. First we used density functional theory molecular dynamics simulations (DFT-MD) on the periodic crystal to discover the initial reaction steps. This allowed us to determine the most important reactions through DFT-MD simulations at high temperatures. Then we started with the midpoint of the reaction (unimolecular or bimolecular) from the DFT-MD and carried out higher quality finite cluster DFT calculations to locate the true transition state of the reaction, followed by calculations along the reaction path to determine the initial and final states. We find that for the noncompressed BCHMX the nitro-aci isomerization reaction occurs earlier than the NO2-releasing reaction, while for compressed BCHMX intermolecular hydrogen-transfer and bimolecular NO2-releasing reactions occur earlier than the nitrous acid (HONO)-releasing reaction. At high pressures, the initial reaction involves intermolecular hydrogen transfer rather than intramolecular hydrogen transfer, and the intermolecular hydrogen transfer decreases the reaction barrier for release of NO2 by ∼7 kcal/mol. Thus, the HONO-releasing reaction takes place more easily in compressed BCHMX. We find that this reaction barrier is 10 kcal/mol lower than the unimolecular NO2 release and ∼3 kcal/mol lower than the bimolecular NO2 release. This rationalizes the origin of the higher sensitivity of BCHMX compared to RDX (1,3,5-trinitrohexahydro-1,3,5-triazine) and HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). We suggest changes in BCHMX that might help decrease the sensitivity by avoiding the intermolecular hydrogen-transfer and HONO-releasing reaction.