A Study of the Shock Sensitivity of Energetic Single Crystals by Large-Scale Ab Initio Molecular Dynamics Simulations

Theoretical investigation of potential energetic material CL-20/TNBP co-crystal explosive based on molecular dynamics method

CONTEXT

The exploration of CL-20 eutectic has been a subject of fervent interest within the realm of high-energy material modification. Through the utilization of density functional and molecular dynamics methods, an investigation into the characteristics of hexanitrohexaazaisowurtzitane (CL-20)/4,4',5,5'-tetranitro-2H,2'H-3,3'-bipyrazole (TNBP)within the molar ratio range of 4:1-1:4 was conducted. This inquiry encompassed the scrutiny of molecular interaction pathways, attachment force, initiating molecular distance, unified energy concentration, and physical characteristics. Furthermore, the EXPLO-5 was harnessed to prognosticate the explosion features and byproducts of unadulterated CL-20, TNBP, and CL-20/TNBP frameworks. The findings delineate a substantial differentiation in the electrostatic charge distribution on the surface between CL-20 and TNBP particles, signifying the preeminence of intermolecular interactions between disparate entities over those within similar entities, thus intimating the plausibility of the eutectic constitution. Remarkably, the identification of maximal attachment force at a molar ratio of 1:1 suggests the heightened likelihood of eutectic formation, propelled primarily by electrostatic and van der Waals forces. The resultant eutectic explosive evinces intermediate reactivity and exemplary mechanical attributes. Moreover, the detonation achievement of the eutectic with a molar proportion of 1:1 straddles that of CL-20 and TNBP, representing a new type of insensitive high-energy material.

METHODS

The testing method employs the Materials Studio software and utilizes the molecular dynamics (MD) method to predict the properties of CL-20/TNBP cocrystals with different ratios and crystal faces. The MD simulation time step is set to 1 fs, and the total MD simulation time is 2 ns. An isothermal-isobaric (NPT) ensemble is used for the 2 ns MD simulation. The COMPASS force field is employed, with the temperature set to 295 K. The prediction of detonation characteristics and products is conducted using the EXPLO-5 software.

Theoretical study of potential energetic material CL-20/DNAN eutectic explosive based on molecular dynamics method

CONTEXT

The exploration of CL-20 eutectic has been a subject of fervent interest within the realm of high-energy material modification. Through the utilization of density functional and molecular dynamics methods, an investigation into the characteristics of hexanitrohexaazaisowurtzitane (CL-20)/2,4-dinitroanisole (DNAN) within the molar ratio range of 9:1-1:9 was conducted. This inquiry encompassed the scrutiny of molecular interaction pathway, attachment force, initiating molecular distance, unified energy concentration, and physical characteristics. Furthermore, EXPLO-5 was harnessed to prognosticate the explosion features and byproducts of unadulterated CL-20, DNAN, and CL-20/DNAN frameworks. The findings delineate a substantial differentiation in the electrostatic charge distribution on the surface between CL-20 and DNAN particles, signifying the preeminence of intermolecular interactions between disparate entities over those within similar entities, thus intimating the plausibility of eutectic constitution. Remarkably, the identification of maximal attachment force at a molar ratio of 4:6 suggests the heightened likelihood of eutectic formation, propelled primarily by electrostatic and van der Waals forces. The resultant eutectic explosive evinces intermediate reactivity and exemplary mechanical attributes. Moreover, the detonation achievement of the eutectic with a molar proportion of 4:6 straddles that of CL-20 and DNAN, representing a new type of insensitive high-energy material.

METHODS

The testing method employs the Materials Studio software and utilizes the molecular dynamics (MD) method to predict the properties of CL-20/DNAN co-crystals with different ratios and crystal faces. The MD simulation time step is set to 1 fs, and the total MD simulation time is 2 ns. An isothermal-isobaric (NPT) ensemble is used for the 2-ns MD simulation. The COMPASS force field is employed, with the temperature set to 295 K. The prediction of detonation characteristics and products is conducted using the EXPLO-5 software.

Insight into the high-temperature oxidation kinetics of acetylene: A first-principles molecular dynamics study

The study on high-temperature oxidation kinetics and kinetic modeling of acetylene (C2H2) has significant importance for its engineering applications. In this paper, the first-principles molecular dynamics method is used to simulate the C2H2 oxidation under high temperatures for the first time. Our results show that there are 38 intermediates and 225 elementary reactions in the process of C2H2 oxidation. The formation mechanisms of "prompt" CO2, as well as gas pollutants CHOCHO and HCOOH are revealed in depth. Four intermediates, CHCHO, CHOCO, CHOCHO and HCOOH, which have significant controversy in current kinetic models, are verified. And a new intermediate, CHOCO2, is discovered. Meanwhile, our simulation also shows that radicals, such as HO2, OH, O, etc. play a key role in promoting the oxidation of intermediates in the early stage of C2H2 oxidation. Combined with quantum chemical calculations, a detailed chemical kinetic model of C2H2/air (FP-C2H2) is built and verified by simulating ignition delay time, species concentration in the flow reactor and premixed laminar flame speed. Our studies provide novel insight for understanding the complex chemical reaction kinetics, and environmental and human health threats from air pollutant formation during C2H2 combustion.

Theoretical investigations on CL-20/ANTA co-crystal explosive via molecular dynamics method

CONTEXT

The study of CL-20 co-crystal has always been a focal point within the field of energetic material modification. In this study, we employed a combination of density functional theory and molecular dynamics simulations to investigate the properties of hexanitrohexaazaisowurtzitane (CL-20)/3-amino-5-nitro-1,2,4-triazole (ANTA) with different molar ratios ranging from 4:1 to 1:4. Additionally, EXPLO-5 software utilized to predict the detonation properties and products of pure CL-20, ANTA, and CL-20/ANTA systems. The results revealed that there was an interaction between CL-20 and ANTA molecules, which had the potential to form a co-crystal. The most likely molar ratio for co-crystal formation was 1:1, and the main driving forces for co-crystal formation were electrostatic force, dispersion force, and van der Waals force. The co-crystal explosive exhibited moderate sensitivity and excellent mechanical properties. Furthermore, the co-crystal detonation performance at a molar ratio of 1:1 was between that of CL-20 and ANTA, representing a new type of insensitive high-energy material.

METHODS

The properties of CL-20/ANTA co-crystal were predicted by molecular dynamics (MD) method under Materials Studio software. For the whole MD simulations, set the temperature at 298 K, and the pressure was 0.0001 GPa. Conducted MD simulation under the NPT ensemble for a total simulation duration of 1 ns. The first 0.5 ns was used for thermodynamic equilibrium, and the last 0.5 ns was used for statistical calculation and analysis. Sampling was recorded every 10 fs during the calculation.

Low-Power Laser Ignition of an Antenna-Type Secondary Energetic Copper Complex: Synthesis, Characterization, Evaluation, and Ignition Mechanism Studies

In recent years, development of new energetic compounds and formulations, suitable for ignition with relatively low-power lasers, is a highly active and competitive field of research. The main goal of these efforts is focused on achieving and providing much safer solutions for various detonator and initiator systems. In this work, we prepared, characterized, and studied thermal and ignition properties of a new laser-ignitable compound, based on the 5,6-bis(ethylnitroamino)-N'2,N'3-dihydroxypyrazine-2,3-bis(carboximidamide) (DS3) proligand. This new energetic proligand was prepared in three steps, starting with 5,6-bis(ethylamino)-pyrazine-2,3-dicarbonitrile. Crystallography studies of the DS3-derived Cu(II) complex (DS4) revealed a unique stacked antenna-type structure of the latter compound. DS4 has an exothermal temperature of 154.5 °C and was calculated to exhibit a velocity of detonation of 6.36 km·s^-1 and a detonation pressure of 15.21 GPa. DS4 showed properties of a secondary explosive, having sensitivity to impact, friction, and electrostatic discharge of 8 J, 360 N, and 12 mJ, respectively. In order to study the mechanism of ignition by a laser (using a diode laser, 915 nm), we conducted a set of experiments that enabled us to characterize a photothermal ignition mechanism. Furthermore, we found that a single pulse, with a time duration of 1 ms and with a total energy of 4.6 mJ, was sufficient for achieving a consistent and full ignition of DS4. Dual-pulse experiments, with variable time intervals between the laser pulses, showed that DS4 undergoes ignition via a photothermal mechanism. Finally, calculating the chemical mechanism of the formation of the complex DS4 and modeling its anhydrous and hydrated crystal structures (density functional theory calculations using Gaussian and HASEM software) allowed us to pinpoint a more precise location of water molecules in experimental crystallographic data. These results suggest that DS4 has potential for further development to a higher technology readiness level and for integration into small-size safe detonator systems as for many civil, aerospace, and defense applications.

Applying machine learning to balance performance and stability of high energy density materials

The long-standing performance-stability contradiction issue of high energy density materials (HEDMs) is of extremely complex and multi-parameter nature. Herein, machine learning was employed to handle 28 feature descriptors and 5 properties of detonation and stability of 153 HEDMs, wherein all 21,648 data used were obtained through high-throughput crystal-level quantum mechanics calculations on supercomputers. Among five models, namely, extreme gradient boosting regression tree (XGBoost), adaptive boosting, random forest, multi-layer perceptron, and kernel ridge regression, were respectively trained and evaluated by stratified sampling and 5-fold cross-validation method. Among them, XGBoost model produced the best scoring metrics in predicting the detonation velocity, detonation pressure, heat of explosion, decomposition temperature, and lattice energy of HEDMs, and XGBoost predictions agreed best with the 1,383 experimental data collected from massive literatures. Feature importance analysis was conducted to obtain data-driven insight into the causality of the performance-stability contradiction and delivered the optimal range of key features for more efficient rational design of advanced HEDMs.

Strategies for Achieving Balance between Detonation Performance and Crystal Stability of High-Energy-Density Materials

Performance-stability contradiction of high-energy-density materials (HEDMs) is a long-standing puzzle in the field of chemistry and material science. Bridging the gap that exists between detonation performance of new HEDMs and their stability remains a formidable challenge. Achieving optimal balance between the two contradictory factors is of a significant demand for deep-well oil and gas drilling, space exploration, and other civil and defense applications. Herein, supercomputers and latest quantitative computational strategies were employed and high-throughput quantum calculations were conducted for 67 reported HEDMs. Based on statistical analysis of large amounts of physico-chemical data, in-crystal interspecies interactions were identified to be the one that provokes the performance-stability contradiction of HEDMs. To design new HEDMs with both good detonation performance and high stability, the proposed systematic and comprehensive strategies must be satisfied, which could promote the development of crystal engineering of HEDMs to an era of theory-guided rational design of materials.

Molecular and Crystal Features of Thermostable Energetic Materials: Guidelines for Architecture of "Bridged" Compounds

Extensive density functional theory (DFT) calculation and data analysis on molecular and crystal level features of 60 reported energetic materials (EMs) allowed us to define key descriptors that are characteristics of these compounds' thermostability. We see these descriptors as reminiscent of "Lipinski's rule of 5", which revolutionized the design of new orally active pharmaceutical molecules. The proposed descriptors for thermostable EMs are of a type of molecular design, location and type of the weakest bond in the energetic molecule, as well as specific ranges of oxygen balance, crystal packing coefficient, Hirshfeld surface hydrogen bonding, and crystal lattice energy. On this basis, we designed three new thermostable EMs containing bridged, 3,5-dinitropyrazole moieties, HL3, HL7, and HL9, which were synthesized, characterized, and evaluated in small-scale field detonation experiments. The best overall performing compound HL7 exhibited an onset decomposition temperature of 341 °C and has a density of 1.865 g cm-3, and the calculated velocity of detonation and maximum detonation pressure were 8517 m s-1 and 30.6 GPa, respectively. Considering HL7's impressive safety parameters [impact sensitivity (IS) = 22 J; friction sensitivity (FS) = 352; and electrostatic discharge sensitivity (ESD) = 1.05 J] and the results of small-scale field detonation experiments, the proposed guidelines should further promote the rational design of novel thermostable EMs, suitable for deep well drilling, space exploration, and other high-value defense and civil applications.