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.

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.

Recent Progress on Synthesis, Characterization, and Performance of Energetic Cocrystals: A Review

In the niche area of energetic materials, a balance between energy and safety is extremely important. To address this "energy-safety contradiction", energetic cocrystals have been introduced. The investigation of the synthesis methods, characteristics, and efficacy of energetic cocrystals is of the utmost importance for optimizing their design and development. This review covers (i) various synthesis methods for energetic cocrystals; (ii) discusses their characteristics such as structural properties, detonation performance, sensitivity analysis, thermal properties, and morphology mapping, along with other properties such as oxygen balance, solubility, and fluorescence; and (iii) performance with respect to energy contents (detonation velocity and pressure) and sensitivity. This is followed by concluding remarks together with future perspectives.

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.

Time for quartet: the stable 3 : 1 cocrystal formulation of FTDO and BTF – a high-energy-density material

A computer simulation of cocrystal structures of [1,2,5]oxadiazolo[3,4-e][1,2,3,4]tetrazine 4,6-dioxide (FTDO) with benzotrifuroxan (BTF) in ratios of (3–1 : 1) was performed. Theoretically and experimentally was shown: a (3 : 1) cocrystal is formed.

Molecular dynamics simulations of a cyclotetramethylene tetra-nitramine/hydrazine 5,5'-bitetrazole-1,1'-diolate cocrystal

An energetic ionic salt (EIS)-based cocrystal formation, cyclotetramethylene tetra-nitramine (HMX)/hydrazine 5,5′-bitetrazole-1,1′-diolate (HA·BTO), is predicted based on molecular dynamics simulations. HA·BTO is a newly-synthesized environmentally friendly energetic ionic salt with good detonation performance and low sensitivity. Calculated powder X-ray diffraction patterns and intermolecular interactions deduce the formation of the new cocrystal structure. Radial distribution function analysis suggests that hydrogen bonds and van der Waals (vdW) forces exist between the H⋯O pairs of HMX and HA·BTO, while the hydrogen bonds between the H of HA·BTO and the O of HMX play a prominent role. The cohesive energy density and mechanical properties are also analyzed. The cohesive energy density of the HMX/HA·BTO cocrystal is larger than that of the composite of HMX and HA·BTO, indicating an improvement in crystal stability by cocrystalization. Compared to both HMX and HA·BTO, HMX/HA·BTO has smaller Young modulus, bulk modulus and shear modulus values, but larger K/G values and a positive Cauchy pressure, suggesting decreased stiffness and improved ductibility. Moreover, the calculated formation energy is −405.79 kJ mol⁻¹ at 298 K, which implies that the proposed cocrystal structure is likely to be synthesized at ambient temperature. In summary, we have predicted an EIS-based cocrystal formation in which the safety and mechanical properties of HMX have been improved via cocrystalization with HA·BTO, and this provides deep insight into the formation mechanism of the EIS-based cocrystal.

An energetic ionic salt-based cocrystal formation, HMX/HA·BTO, is predicted based on molecular dynamics simulations.

Design, preparation, characterization and formation mechanism of a novel kinetic CL-20-based cocrystal

2,4,6,8,10,12-Hexanitrohexaazaisowurtzitane (CL-20)-based cocrystals have gained increasing attention as a means of obtaining insensitive high explosives. However, the design of ideal candidates for these cocrystals remains difficult. This work compares the crystal energies of the CL-20-dinitrobenzene (DNB) and CL-20-2,4,6-trinitrotoluene (TNT) cocrystals with those of the respective pure coformers. The results indicate that the cocrystal formation is driven by the differences in the energies of the cocrystals and the coformers. Furthermore, analysis via Hirshfeld surfaces and two-dimensional fingerprint plots confirms that the O...O, O...H, O...N and C...O interactions were the main force for stabilizing the CL-20-based cocrystal structure. Based on these findings, a novel energetic-energetic cocrystal of CL-20-2,4,6-trinitrophenol (TNP) was designed and prepared by means of a rapid method for solvent removal. The crystal structure was investigated via powder X-ray diffraction methods, solid-state nuclear magnetic resonance spectroscopy and Fourier transform infrared spectroscopy. The results revealed that the O-H...O hydrogen bonding interaction between the phenolic hydroxyl group of TNP and nitro groups of CL-20, as well as nitro...π, nitro...nitro and O_NO2...π(N)_NO2 interactions, based on the benzene ring and nitro groups, are the main interactions occurring in the cocrystal.

Optimizing the molecular structure and packing style of a crystal by intramolecular cyclization from picrylhydrazone to indazole

Crystal engineering has prompted the development of energetic materials in recent years.

High-density HNIW/TNT cocrystal synthesized using a green chemical method

The main challenge for achieving better energetic materials is to increase their density. In this paper, cocrystals of HNIW (2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane, often referred to as CL-20) with TNT (2,4,6-trinitrotoluene) were synthesized using ethanol in a green chemical method. The cocrystal was formulated as C13H11N15O18 and possesses a higher density (1.934 g cm-3) than published previously (1.846 g cm-3). This high-density cocrystal possesses a new structure, which can be substantiated by the different types of hydrogen bonds. The predominant driving forces that connect HNIW with TNT in the new cocrystal were studied at ambient conditions using single-crystal X-ray diffraction, powder X-ray diffraction, Fourier transform-infrared spectroscopy and Raman spectroscopy. The results reveal that the structure of the new HNIW/TNT cocrystals consists of three one-dimensional hydrogen-bonded chains exploiting the familiar HNIW-TNT multi-component supramolecular structure, in which two hydrogen-bonded chains are between -NO2 (HNIW) and -CH (TNT), and one hydrogen-bonded chain is between -CH (HNIW) and -NO2 (TNT). The changes to the electron binding energy and type of element in the new cocrystal were traced using X-ray photoelectron spectroscopy. Meanwhile, the physicochemical characteristics alter after cocrystallization due to the hydrogen bonding. It was found that the new HNIW/TNT cocrystal is more thermodynamically stable than HNIW. Thermodynamic aspects of new cocrystal decomposition are investigated in order to explain this observation. The detonation velocity of new HNIW/TNT cocrystals is 8631 m s-1, close to that of HNIW, whereas the mechanical sensitivity is lower than HNIW.

Synthesis, Characterization, and Sensitivity of a CL-20/PNCB Spherical Composite for Security

Highly energetic materials have received significant attention, particularly 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20). However, the application of this material was limited due to its high sensitivity. It is well known that the shape, size, and structure of energetic materials (EMs) significantly influence their sensitivity. At present, there are several ways to reduce the sensitivity of CL-20, such as spheroidization, ultrafine processing, and composite technology. However, only one or two of the abovementioned methods have been reported in the literature, and the obtained sensitivity effect was unsatisfactory. Thus, we tried to further reduce the sensitivity of CL-20 by combining the above three methods. The as-prepared composite was precipitated from the interface between two solutions of water and ethyl acetate, and the composite was insensitive compared with other reported CL-20-based EMs. The H₅₀ value for the composite ranged up to 63 cm. This approach opens new prospects for greatly reducing the sensitivity of high Ems.

Theoretical investigation of the effects of the molar ratio and solvent on the formation of the pyrazole-nitroamine cocrystal explosive 3,4-dinitropyrazole (DNP)/2,4,6,8,10,12-hexanitrohexaazaisowurtzitane (CL-20)

The effects of the molar ratio, temperature, and solvent on the formation of the cocrystal explosive DNP/CL-20 were investigated using molecular dynamics (MD) simulation. The cocrystal structure was predicted through Monte Carlo (MC) simulation and using first-principles methods. The results showed that the DNP/CL-20 cocrystal might be more stable in the molar ratio 1:1 near to 318 K, and the most probable cocrystal crystallizes in the triclinic crystal system with the space group P[Formula: see text]. Cocrystallization was more likely to occur in methanol and ethanol at 308 K as a result of solvent effects. The optimized structure and the reduced density gradient (RDG) of the DNP/CL-20 complex confirmed that the main driving forces for cocrystallization were a series of hydrogen bonds and van der Waals forces. Analyses of the trigger bonds, the charges on the nitro groups, the electrostatic surface potential (ESP), and the free space per molecule in the cocrystal lattice were carried out to further explore their influences on the sensitivity of CL-20. The results indicated that the DNP/CL-20 complex tended to be more stable and insensitive than pure CL-20. Moreover, an investigation of the detonation performance of the DNP/CL-20 cocrystal indicated that it possesses high power. Graphical abstract DNP/CL-20 cocrystal models with different molar ratios were investigated at different temperatures using molecular dynamics (MD) simulation methods. Binding energies and mechanical properties were probed to determine the stability and performance of each cocrystal model. Solvated DNP/CL-20 models were established by adding solvent molecules to the cocrystal surface. The binding energies of the models in various solvents were calculated in order to identify the most suitable solvent and temperature for preparing the cocrystal explosive DNP/CL-20.