Molecular dynamics simulations of void defects in the energetic material HMX

Abstract

A molecular dynamics (MD) simulation was carried out to characterize the dynamic evolution of void defects in crystalline octahydro-1, 3, 5, 7-tetranitro-1, 3, 5, 7-tetrazocine (HMX). Different models were constructed with the same concentration of vacancies (10 %) to discuss the size effects of void. Energetic ground state properties were determined by annealing simulations. The void formation energy per molecule removed was found to be 55-63 kcal/mol(-1), and the average binding energy per molecule was between 32 and 34 kcal/mol(-1) according to the change in void size. Voids with larger size had lower formation energy. Local binding energies for molecules directly on the void surface decreased greatly compared to those in defect-free lattice, and then gradually increased until the distance away from the void surface was around 10 Å. Analysis of 1 ns MD simulations revealed that the larger the void size, the easier is void collapse. Mean square displacements (MSDs) showed that HMX molecules that had collapsed into void present liquid structure characteristics. Four unique low-energy conformers were found for HMX molecules in void: two whose conformational geometries corresponded closely to those found in HMX polymorphs and two, additional, lower energy conformers that were not seen in the crystalline phases. The ratio of different conformers changed with the simulated temperature, in that the ratio of α conformer increased with the increase in temperature.

The growth and perfection of β-cyclotetramethylene-tetranitramine (HMX) studied by laboratory and synchrotron X-ray topography

“…[8] The nature of internal defects is still important however, especially when attempting to understand ignition processes at elevated temperatures. [9][10][11]. At these temperatures, the breakdown in mechanical properties and changes within crystalline particulates can affect localised ongoing heating and decomposition processes.…”

Molecular Dynamics Simulation of the Influence of RDX Internal Defects on Sensitivity

“…Then, the RDX crystal can be optimized in a targeted manner. The molecular dynamics simulation was used to simulate many properties before, such as mechanical properties, prediction of crystal morphology, and so on, and the results are all in good agreement with experiments [13][14][15][16][17][18]. The simulation method provides a safe, convenient, fast, and low-cost way to study the influence of internal defects on the sensitivity of RDX crystals.…”

Molecular Dynamics Simulation of β-HMX Crystal Morphology Induced by Polymer Additives

“…A significant trend is that energetic materials need to have higher detonation performance and lower sensitivity. With the in-depth study of high energy and insensitive explosives, the influence of crystal shape, crystal surface state, and crystal internal defects on the performance of explosives has attracted wide attention from researchers [6][7][8][9][10][11][12][13]. Compared to the explosive crystals in the form of flakes and needles, spherical explosives can significantly increase the charge density and improve the detonation performance of the charge.…”