First-principles high-pressure unreacted equation of state and heat of formation of crystal 2,6-diamino-3, 5-dinitropyrazine-1-oxide (LLM-105)

Pressure and temperature effects on the Raman spectra of LLM-105

Currently, only one crystal structure of LLM-105 (2,6-diamino-3,5-dinitropyrazine-1-oxide) (P21/n) has been discovered, and there are still debates on its phase transition point and phase diagram. Based on previous work, we performed crystal structure, Raman spectra, and vibrational properties calculations on LLM-105 crystal. Our results indicate that the crystal structure of LLM-105 remains stable until compressed to 49 GPa, beyond which it may undergo two phase transitions at pressure intervals of 49.0-49.1 GPa and 51.4-51.5 GPa, respectively. Analysis of Raman shift results suggests that these two phase transitions may be reversible, with an intermediate phase possibly acting as a transition phase. Additionally, based on the quasi-harmonic approximation, we fitted the experimental data of LLM-105 lattice expansion state, obtaining the volume at zero pressure and using it for Raman spectra calculations. The results demonstrated the accuracy of this quasi-harmonic approximation method in describing the redshift of Raman peaks during the heating process and the excitation ratio of Raman peaks in different wavenumber ranges.

Review of Applications of Density Functional Theory (DFT) Quantum Mechanical Calculations to Study the High-Pressure Polymorphs of Organic Crystalline Materials

Since its inception, chemistry has been predominated by the use of temperature to generate or change materials, but applications of pressure of more than a few tens of atmospheres for such purposes have been rarely observed. However, pressure is a very effective thermodynamic variable that is increasingly used to generate new materials or alter the properties of existing ones. As computational approaches designed to simulate the solid state are normally tuned using structural data at ambient pressure, applying them to high-pressure issues is a highly challenging test of their validity from a computational standpoint. However, the use of quantum chemical calculations, typically at the level of density functional theory (DFT), has repeatedly been shown to be a great tool that can be used to both predict properties that can be later confirmed by experimenters and to explain, at the molecular level, the observations of high-pressure experiments. This article's main goal is to compile, analyze, and synthesize the findings of works addressing the use of DFT in the context of molecular crystals subjected to high-pressure conditions in order to give a general overview of the possibilities offered by these state-of-the-art calculations.

Ab Initio Crystal Structure Prediction of the Energetic Materials LLM-105, RDX, and HMX

Crystal structure prediction (CSP) is performed for the energetic materials (EMs) LLM-105 and α-RDX, as well as the α and β conformational polymorphs of 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX), using the genetic algorithm (GA) code, GAtor, and its associated random structure generator, Genarris. Genarris and GAtor successfully generate the experimental structures of all targets. GAtor's symmetric crossover scheme, where the space group symmetries of parent structures are treated as genes inherited by offspring, is found to be particularly effective. However, conducting several GA runs with different settings is still important for achieving diverse samplings of the potential energy surface. For LLM-105 and α-RDX, the experimental structure is ranked as the most stable, with all of the dispersion-inclusive density functional theory (DFT) methods used here. For HMX, the α form was persistently ranked as more stable than the β form, in contrast to experimental observations, even when correcting for vibrational contributions and thermal expansion. This may be attributed to insufficient accuracy of dispersion-inclusive DFT methods or to kinetic effects not considered here. In general, the ranking of some putative structures is found to be sensitive to the choice of the DFT functional and the dispersion method. For LLM-105, GAtor generates a putative structure with a layered packing motif, which is desirable thanks to its correlation with low sensitivity. Our results demonstrate that CSP is a useful tool for studying the ubiquitous polymorphism of EMs and shows promise of becoming an integral part of the EM development pipeline.

Pressure-Modulated Dissolution Behavior of LLM-105 Crystals in High-Temperature Water

The exploration of the microstructural evolution and reaction kinetics of energetic materials with high-temperature and high-pressure water contributes to the understanding of their microscopic physicochemical origin, which can provide critical experimental data for the use of energetic materials. As a promising high-energy and insensitive energetic material, LLM-105 has been investigated under extreme conditions such as high pressure and high temperature. However, little information is available about the effect of water on LLM-105 under high pressure and high temperature. In this work, the interaction between LLM-105 and water under HP-HT was investigated in detail. As a result, the dissolving behavior of LLM-105 in water under high pressure and high temperature is related to the initial pressure. When the initial pressure is less than 1 GPa, LLM-105 crystals are dissolved in high-temperature water; when the initial pressure is above 1 GPa, LLM-105 particles are only decomposed in high-temperature water. When the solution is saturated at a high temperature, recrystallization of the LLM-105 sample appears in the solution. High pressure hindered the dissolution process of the sample in HP-HT water because the interaction between the solute and the solvent was weakened by high pressure. The initial pressure is one of the significant parameters that determines whether LLM-105 crystals can be dissolved in high-temperature water. More importantly, water under high pressure and high temperature can not only act as a solvent when dissolving the samples but also act as a catalyst to accelerate the decomposition process. In addition, the HP-HT water reduced the decomposition temperature of the LLM-105 crystal to a large extent. The research in this paper not only provides insights into the interaction between LLM-105 and water but also contributes to the performance of energetic materials under extreme conditions and their practical applications in complex conditions.

Pressure effects on the thermal decomposition of the LLM-105 crystal

Thermal mechanical responses under high temperature and high pressure are basic information to understand the performance of energetic materials. In this work, the pressure effects on the thermal decay of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105) are explored. Up to the initial pressure of 4.6 GPa, the pressure dependent decomposition boundary is built and no phase transition occurs until the decomposition of the LLM-105 crystal. The decomposition temperature is significantly lifted via a weak loading pressure. The experimental measurement confirms the decomposition products, including NO₂, CO₂ and NH₃, which are predicted by the density functional tight-binding molecular dynamics method. The calculation described the details of thermal decay in the initial stages under high pressure. The sudden drop in the shifts of the Raman modes associated with hydrogen bonds under high pressure indicates the strengthening of the intermolecular hydrogen bonds and the occurrence of intermolecular hydrogen transfer prior to crystal decomposition. The simulation supported the existence of intermolecular hydrogen transfer and provided the transfer path and decomposition mechanism. All of these jobs not only contribute significantly to the understanding of thermal decomposition of energetic materials as a function of pressure, but also contribute to the understanding of sensitivity mechanisms and safety issues.

Compression Behavior and Vibrational Properties of New Energetic Material LLM-105 Analyzed Using the Dispersion-Corrected Density Functional Theory

The 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a newly energetic material with an excellent performance and low sensitivity and has attracted considerable attention. On the basis of the dispersion-corrected density functional theory (DFT-D), the high-pressure responses of vibrational properties, in conjunction with structural properties, are used to understand its intermolecular interactions and anisotropic properties under hydrostatic and uniaxial compressions. At ambient and pressure conditions, the DFT-D scheme could reasonably describe the structural parameters of LLM-105. The hydrogen bond network, resembling a parallelogram shape, links two adjacent molecules and contributes to the structure stability under hydrostatic compression. The anisotropy of LLM-105 is pronounced, especially for Raman spectra under uniaxial compression. Specifically, the red-shifts of modes are obtained for [100] and [010] compressions, which are caused by the pressure-induced enhance of the strength of the hydrogen bonds. Importantly, coupling modes and discontinuous Raman shifts are observed along [010] and [001] compressions, which are related to the intramolecular vibrational redistribution and possible structural transformations under uniaxial compressions. Overall, the detailed knowledge of the high-pressure responses of LLM-105 is established from the atomistic level. Uniaxial compression responses provide useful insights for realistic shock conditions.

Predicted Reaction Mechanisms, Product Speciation, Kinetics, and Detonation Properties of the Insensitive Explosive 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a relatively new and promising insensitive high-explosive (IHE) material that remains only partially characterized. IHEs are of interest for a range of applications and from a fundamental science standpoint, as the root causes behind insensitivity are poorly understood. We adopt a multitheory approach based on reactive molecular dynamic simulations performed with density functional theory, density functional tight-binding, and reactive force fields to characterize the reaction pathways, product speciation, reaction kinetics, and detonation performance of LLM-105. We compare and contrast these predictions to 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a prototypical IHE, and 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX), a more sensitive and higher performance material. The combination of different predictive models allows access to processes operative on progressively longer timescales while providing benchmarks for assessing uncertainties in the predictions. We find that the early reaction pathways of LLM-105 decomposition are extremely similar to TATB; they involve intra- and intermolecular hydrogen transfer. Additionally, the detonation performance of LLM-105 falls between that of TATB and HMX. We find agreement between predictive models for first-step reaction pathways but significant differences in final product formations. Predictions of detonation performance result in a wide range of values, and one-step kinetic parameters show the similar reaction rates at high temperatures for three out of four models considered.

Learning to fly: thermochemistry of energetic materials by modified thermogravimetric analysis and highly accurate quantum chemical calculations

The standard state enthalpy of formation and the enthalpy of sublimation are essential thermochemical parameters determining the performance and application prospects of energetic materials (EM). Direct experimental measurements of these properties are complicated by low volatility and high heat release in bomb calorimetry experiments. As a result, the uncertainties in the reported enthalpies of formation for a number of even well-known CHNO-containing compounds might amount up to tens kJ mol-1, while for some novel high-nitrogen molecules they reach even hundreds of kJ mol-1. The present study reports a facile approach to determining the solid-state formation enthalpies comprised of complementary high-level quantum chemical calculations of the gas-phase thermochemistry and advanced thermal analysis techniques yielding sublimation enthalpies. The thermogravimetric procedure for the measurement of sublimation enthalpy was modified by using low external pressures (down to 0.2 Pa). This allows for observing sublimation/vaporization instead of thermal decomposition of the compounds studied. Extensive benchmarking on nonenergetic and energetic compounds reveals the average and maximal absolute errors of the sublimation enthalpies of 3.3 and 11.0 kJ mol-1, respectively. The comparison of the results with those obtained from the widely used Trouton-Williams empirical equation shows that the latter underestimates the sublimation enthalpy up to 140 kJ mol-1. Therefore, we performed a reparametrization of the latter equation with simple chemical descriptors that reduces the mean error down to 30 kJ mol-1. Highly accurate multi-level procedures W2-F12 and/or W1-F12 in conjunction with the atomization energy approach were used to calculate theoretically the gas-phase formation enthalpies. In several cases, the DLPNO-CCSD(T) enthalpies of isodesmic reactions were also employed to obtain the gas-phase thermochemistry for medium-sized important EMs. Combining the obtained thermochemical properties, we determined the solid-state enthalpies of formation for nearly 60 species containing various important explosophoric groups, from common nitroaromatics, nitroethers, and nitramines to novel nitrogen-rich heterocyclic species (e.g., the derivatives of pyrazole, tetrazole, furoxan, etc.). The large-scale benchmarking against the available experimental solid-state enthalpies of formation yielded the maximal inaccuracy of the proposed method of 25 kJ mol-1.

First-Principles-Based Force Field for 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a highly promising energetic material (EM) with high safety. Understanding its microscopic response mechanisms within the external stimulus is meaningful for the design of EMs. In order to comprehend the complicated phenomena, it is necessary to employ molecular simulation methods to investigate the response mechanisms with the force field (FF) at an atomic level. In this work, we developed a tailored FF for LLM-105 based on first-principles calculations. The validity of the FF was evaluated by molecular dynamics simulations. The structural parameters of LLM-105 predicted by FF are in good agreement with the experimental values, such as lattice constant, bond length, bond angle, dihedral angle and center of mass, and so forth. Moreover, the FF possesses good performance to describe the structural response on pressure accurately. In general, our work not only builds a balanced FF in gas and condensed phases, but also provides a useful tool to study the properties about LLM-105 at a large scale, which is helpful to improve the understanding about the balance between energy and safety in EMs.

An Electrospun Preparation of the NC/GAP/Nano-LLM-105 Nanofiber and Its Properties

In this work, an energetic composite fiber, in which 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105) nanoparticles intimately incorporated with a nitrocellulose/glycidyl azide polymer (NC/GAP) fiber, was prepared by the electrospinning method. The morphology and structure of the nanofiber was characterized by scanning electron microscopy (SEM), energy dispersive X-Ray (EDX), fourier transform infrared spectroscopy (IR), X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), and Brunauer-Emmett-Teller (BET). The nanofibers possessed a three-dimensional (3D) net structure and a large specific surface area. Thermal analysis, energetic performance, and sensitivities were investigated, and they were compared with NC/GAP and LLM-105 nanoparticles. The NC/GAP/nano-LLM-105 nanofibers show higher decomposition rates and lower decomposition temperatures. The NC/GAP/nano-LLM-105 decomposed to CO2, CO, H2O, N2O, and a few NO, -CH2O-, and -CH- fragments, in the thermal-infrared spectrometry online (TG-IR) measurement. The NC/GAP/nano-LLM-105 nanofibers demonstrated a higher standard specific impulse (Isp), a higher combustion chamber temperature (Tc), and a higher specialty height (H50). The introduction of nano-LLM-105 in the NC/GAP matrix results in an improvement in energetic performance and safety.

Computational analysis of mesoscale thermomechanical ignition behavior of impacted LLM-105 based explosives

LLM-105 (2,6-diamino-3,5-dinitropyrazine-1-oxide) is an insensitive high explosive crystal which has performance between that of HMX and TATB. An elastoviscoplastic dislocation model is developed for LLM-105 crystal, which accounts for the dislocation evolutions at the crystal interior and crystal wall and strain-rate dependent work hardening. Three different crystal morphology (cubic, icosahedral, rodlike) of LLM-105 based explosive computational models were constructed and subjected to an impact velocity of 200 m s⁻¹ and 500 m s⁻¹. Effects of crystal morphology and initial dislocation density on thermomechanical ignition behavior of LLM-105 based explosives were analyzed. Dislocation density of both crystal interiors and crystal walls in the rodlike LLM-105 based explosive increases slower than that in the cubic and icosahedral explosives. Both the volume averaged and localized stress and dislocation density are the lowest for the rodlike explosive. At the impact velocity of 500 m s⁻¹, a temperature rise due to volumetric work, plasticity work and chemical reaction is sufficiently high to lead to the ignition of the cubic explosive, which shows that the rodlike explosive is the least sensitive among the three explosives. Moreover, with the increase of initial dislocation density, the corresponding volume averaged and localized stress and temperature increase as well. Results presented bridge the macroscale thermomechanical ignition response with the mesoscale deformation mechanisms, which is essential for better understanding the ignition mechanisms and guiding the design of LLM-105 based formulations.

Mesoscale thermomechanical ignition behavior of LLM-105 based explosives was quantified through a developed elastoviscoplastic dislocation model for LLM-105 crystal.

Initial Mechanisms for the Unimolecular Thermal Decomposition of 2,6-Diamino-3,5-dinitropyrazine-1-oxide

The initial channels of thermal decomposition mechanism of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105) molecule were investigated. The results of quantum chemical calculations revealed four candidates involved in the reaction pathway, including the C⁻NO₂ bond homolysis, nitro⁻nitrite rearrangement followed by NO elimination, and H transfer from amino to acyl O and to nitro O with the subsequent OH or HONO elimination, respectively. In view of the further kinetic analysis and ab initio molecular dynamics simulations, the C⁻NO₂ bond homolysis was suggested to be the dominant step that triggered the decomposition of LLM-105 at temperatures above 580 K. Below this temperature, two types of H transfer were considered as the primary reactions, which have advantages including lower barrier and high rate compared to the C⁻NO₂ bond dissociation. It could be affirmed that these two types of H transfer are reversible processes, which could buffer against external thermal stimulation. Therefore, the excellent thermal stability of LLM-105, that is nearly identical to that of 1,3,5-triamino-2,4,6-trinitrobenzene, can be attributed to the reversibility of H transfers at relatively low temperatures. However, subsequent OH or HONO elimination reactions occur with difficulty because of their slow rates and extra energy barriers. Although nitro⁻nitrite rearrangement is theoretically feasible, its rate constant is too small to be observed. This study facilitates the understanding of the essence of thermal stability and detailed decomposition mechanism of LLM-105.

Structural, mechanical properties, and vibrational spectra of LLM-105 under high pressures from a first-principles study

In this work, we report the structure, mechanical properties, and vibrational spectra of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105), an energetic molecular crystal, with a first-principles method based on density functional theory (DFT) using the recentely developped HASEM package. The elastic constants, acoustic velocity, and parameters of equations of state were calculated, and the predicted ordering of stiffness constants is C ₃₃ (38.5 GPa) > C ₁₁ (24.0 GPa) > C ₂₂ (17.7 GPa). We also investigated the structure and equation of state of LLM-105 under hydrostatic pressure up to 100 GPa. The predicted structures are in good agreement with experimental results available from ambient pressure to 20 GPa. Under compressions, the LLM-105 crystal exhibits anisotropic compressibility, with a highly incompressible response along the a-axis and c-axis. It is worth noting that there is a sudden change in the lattice parameters and change rate of volume at ~30 GPa. Based on the intermolecular interaction analysis and vibrational spectra, a phase transition at the hydrostatic pressure of ~30 GPa is predicted.