Molecular dynamics (MD) study to predict performances of a novel hexanitrohexaazaisowurtzitane/1,4-dinitroimidazole (CL-20/1,4-DNI) energetic cocrystal under different temperatures

CONTEXT

To estimate the influence of temperature on properties of 2,4,6,8,10,12-hexanitro- 2,4,6,8,10,12-hexaazaisowurtzitane/1,4-dinitroimidazole (CL-20/1,4-DNI) cocrystal explosive, the supercell crystal of CL-20/1,4-DNI cocrystal model was established. The mechanical properties, sensitivity, and stability of cocrystal model under different temperatures (T = 225 K, 250 K, 275 K, 300 K, 325 K, 350 K) were predicted. Results show that mechanical parameters, including bulk modulus, tensile modulus and shear modulus are the lowest when temperature is 300 K, while Cauchy pressure is the highest, indicating that CL-20/1,4-DNI cocrystal model has better mechanical properties at 300 K. Cohesive energy density (CED) and its components energies decrease monotonically with the increase of temperature, illustrating that the CL-20 and 1,4-DNI molecules are activated and the safety of cocrystal explosive is worsened with the increase of temperature. Cocrystal model has relatively higher binding energy when the temperature is 300 K, implying that the CL-20/1,4-DNI cocrystal explosive is more stable under this condition.

METHODS

The CL-20/1,4-DNI cocrystal model was optimized and the properties were predicted through molecular dynamics (MD) method. The MD simulation was performed with COMPASS force field and the ensemble was set as NPT, external pressure was set as 0.0001 GPa.

Crystal Structure and Noncovalent Interactions of Heterocyclic Energetic Molecules

Nitrogen-rich heterocyclic compounds are important heterocyclic substances with extensive future applications for energetic materials due to their outstanding density and excellent physicochemical properties. However, the weak intermolecular interactions of these compounds are not clear, which severely limits their widespread application. Three nitrogen-rich heterocyclic compounds were chosen to detect their molecular geometry, stacking mode and intermolecular interactions by crystal structure, Hirshfeld surface, RDG and ESP. The results show that all atoms in each molecule are coplanar and that the stacking mode of the three crystals is a planar layer style. A large amount of inter- and intramolecular interaction exists in the three crystals. All principal types of intermolecular contacts in the three crystals are N···H interactions and they account for 40.9%, 38.9% and 32.9%, respectively. Hydrogen bonding, vdW interactions and steric effects in Crystal c are stronger than in Crystals a and b. The negative ESPs all concentrate on the nitrogen atoms in the three molecules. This work is expected to benefit the crystal engineering of heterocyclic energetic materials.

Structure, intermolecular interactions, and dynamic properties of NTO crystals with impurity defects: a computational study

Frontier orbitals distribute in the position of impurity molecules, whose adjacent NTO molecules begin to decompose first.

Computational insights into the formation driving force of CL-20 based solvates and their desolvation process

The O⋯H interactions between CL-20 and solvents are the formation driving force of CL-20 solvates.

Anisotropic response of the co-crystal of CL-20/TNT under shock loading: molecular dynamics simulation

Anisotropic response of the co-crystal of CL-20/TNT under shock loading is studied by analyzing the changes of thermodynamic parameters, product evolution and cluster evolution.

Screening for energetic compounds based on 1,3-dinitrohexahydropyrimidine skeleton and 5-various explosopheres: molecular design and computational study

In this paper, twelve 1,3-dinitrohexahydropyrimidine-based energetic compounds were designed by introducing various explosopheres into hexahydropyrimidine skeleton. Their geometric and electronic structures, heats of formation (HOFs), energetic performance, thermal stability and impact sensitivity were discussed. It is found that the incorporation of electron-withdrawing groups (-NO2, -NHNO2, -N3, -CH(NO2)2, -CF(NO2)2, -C(NO2)3) improves HOFs of the derivatives and all the substituents contribute to enhancing the densities and detonation properties (D, P) of the title compounds. Therein, the substitution of -C(NO2)3 features the best energetic performance with detonation velocity of 9.40 km s-1 and detonation pressure of 40.20 GPa. An analysis of the bond dissociation energies suggests that N-NO2 bond may be the initial site in the thermal decompositions for most of the derivatives. Besides, -ONO2 and -NF2 derivatives stand out with lower impact sensitivity. Characters with striking detonation properties (D = 8.62 km s-1, P = 35.08 GPa; D = 8.81 km s-1, P = 34.88 GPa), good thermal stability, and acceptable impact sensitivity (characteristic height H50 over 34 cm) lead novel compounds 5,5-difluoramine-1,3-dinitrohexahydropyrimidine (K) and 5-fluoro-1,3,5-trinitrohexahydropyrimidine (L) to be very promising energetic materials. This work provides the theoretical molecular design and a reasonable synthetic route of L for further experimental synthesis and testing.

Theoretical Study on CL-20-Based Cocrystal Energetic Compounds in an External Electric Field

An external electric field has great effects on the sensitivity of cocrystal energetic materials. In order to find out the relationship between the external electric field and sensitivity of cocrystals 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane/benzotrifuroxan (CL-20/BTF), 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane/3,4-dinitropyrazole (CL-20/DNP), and 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane/1-methyl-3,5-dinitro-1,2,4-triazole (CL-20/MDNT), density functional theory at B3LYP-D3/6-311+G(d,p) and M062X-D3/ma-def2 TZVPP levels was employed to calculate frontier molecular orbitals, atoms in molecules (AIM) electron density values, bond dissociation energies (BDEs) of the N-NO2 bond, impact sensitivity (H 50), electrostatic potentials (ESPs), and nitro group charges (Q NO2 ) in this work. The results show that a smaller highest occupied molecular orbital-lowest unoccupied molecular orbital gap and the BDEs, as well as H 50, tend to have a larger sensitivity along with the positive directions in the external electric field. Moreover, a smaller local positive ESP (V s max) leads to better stability in the negative electric field. The sensitivity of cocrystal molecules decreases gradually in the negative external electric field with the increase of negative nitro group charges. Finally, the change in the bond lengths, AIM electron density values, and nitro group charges correlate well with the external electric field strengths.

Decomposition mechanism scenarios of CL-20 co-crystals revealed by ReaxFF molecular dynamics: similarities and differences

Understanding the similarities and differences of decomposition mechanisms of CL-20 and its cocrystals is of great interest for practical applications of CL-20 cocrystals. The responses of CL-20 cocrystals to thermal stimulus were investigated using ReaxFF molecular dynamics simulations of two representative cocrystals, CL-20/HMX and CL-20/TNT, under adiabatic conditions and comparing to the baseline system of pure CL-20. The comprehensive chemical details were revealed with the aid of the unique code of VARxMD. The three CL-20-involved reactive systems all exhibit a distinct three-stage character during adiabatic decomposition when using the double peaks of the major intermediate NO₂ amount as the boundary. By taking advantage of the three-stage classification, a clear scenario for the similar stimulus-response of the CL-20 cocrystals can be elucidated for the dominant primary decomposition of CL-20 in stage I and the transition of favored chemical mechanisms from the generation of intermediates/radicals in stage II into their consumption to form stable products in stage III. The similar chemical behaviors are rooted in the dominance of CL-20 chemistry in the initial response of its cocrystals to thermal stimulus. The prolonged reaction zone uncovers the slowed decomposition kinetics of CL-20/HMX and CL-20/TNT, which is associated with the altered kinetics of CL-20 decomposition specifically by N-NO₂ bond scission and CL-20 skeleton decay. The retarded CL-20 decomposition in its cocrystals consequently results in more moderate self-heating and less violent exothermic reactions that agrees with the experimental observations of improved stability and damaged detonation performance of CL-20 cocrystals, particularly for CL-20/TNT. The results obtained in this work suggest that ReaxFF MD simulations can provide useful insight for the modulated chemical properties of varied CL-20 cocrystals.

Seeking a novel energetic co-crystal strategy through the interfacial self-assembly of CL-20 and HMX nanocrystals

Solvent and heat induced self-assembly to CL-20/HMX co-crystals has been investigated. The mechanism towards such process could be concluded to nanoparticle inducing, oriented aggregation, surface integration and co-crystals assembly.

From BTO2− to HBTO− insensitive energetic salt: a route to boost energy

A promising strategy was utilized to boost the detonation performance of insensitive energetic salts.

Preparation of CL-20/TFAZ cocrystals under aqueous conditions: balancing high performance and low sensitivity

CL-20/TFAZ cocrystal explosives were prepared by a self-assembly method under aqueous conditions with both low sensitivity and high detonation performances.