Molecular dynamic simulations on TKX-50/RDX cocrystal

Molecular dynamics simulation of sensitivity of HMX, FOX-7, and TATB crystals

METHODS

This study used molecular dynamics (MD) to simulate three materials (HMX, FOX-7, and TATB) under the NVT ensemble. Six temperatures (100 K, 200 K, 300 K, 400 K, 500 K, and 600 K) were simulated. In addition, the trigger bond lengths, energy bands, and density of states of three materials were obtained at different temperatures and compared with the calculated results at 0 K.

CONTEXT

The results indicate that the trigger bond lengths of the three materials are very close to the experimental values. Overall, the maximum and average bond lengths of the trigger bonds increase with increasing temperature. The band gap value decreases with increasing temperature. The changes in trigger bond length and band gap value are consistent with the experimental fact that sensitivity increases with increasing temperature. And Eg > 1 eV is consistently found within the temperature range of 0-600 K, indicating that all three materials are non-metallic compounds.

Theoretical study on intra-molecule interactions in TKX-50

To fully and deeply understand the weak interactions in the gaseous structure of the TKX-50 molecule, two conformations I and II of the TKX-50 molecule confirmed in a crystal cell were optimized at the B3LYP/6-311g(d,p) level in the gas state, and the single point energy of the optimized structure was calculated at the M06-2X/ma-TZVPP level. Analyzing methods for weak interactions such as the interaction region indicator (IRI), topological basin analysis, and the extended transition state-natural orbitals for chemical valence (ETS-NOCV) theory with the help of Multiwfn code were employed to reveal the corresponding intramolecular weak interactions. The results showed that there were 5 kinds of intramolecular weak interaction in both conformations. They are two types of H bond, two types of intra-ring weak interaction, and one type of O-N bond within the molecular fragment containing the bis-tetrazole ring. The combined effect of all these weak interactions holds the bis-tetrazole ring of TKX-50 retaining an almost coplanar configuration. Meanwhile, the strength of these weak interactions is significantly different in conformation I and conformation II. The most obvious difference is that conformation II has a significant H transfer between intramolecular fragments due to the mirror rotation of almost 180° of cations (NH3OH)+ perpendicular to the N-O bond axis thereof as compared to the reference conformation I. This conformational difference not only makes the weak interaction between the two conformations very different but also forms a quasi-covalent bond in conformation II with much larger bonding energy than other H bonds, thus resulting in conformation II having lower electron energy and more stable geometry. In addition, the order of breaking various H bonds in the combustion decomposition process of TKX-50 is deduced by comparing various H bonds.

Interfacial engineered RDX/TATB energetic co-particles for enhanced safety performance and thermal stability

1,3,5-trinitro-1,3,5-triazinane (RDX) has attracted considerable attentions in energy related fields. However, the safety performance of RDX needs to be improved in terms of various external stimulations. Herein, such issues of...

Comparative investigation on the thermostability, sensitivity, and mechanical performance of RDX/HMX energetic cocrystal and its mixture

Molecular mechanics (MM) and molecular dynamics (MD) simulation method were applied to explore the impact of temperature (220-380 K) on the thermostability, sensitivity, and mechanical performance of RDX (1,3,5-trinitro-1,3,5-triazacyco-hexane)/HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane) energetic cocrystal and mixture models. The mechanical property, the maximum trigger bond length ([Formula: see text]), binding energy, and cohesive energy density (CED) of the pure RDX, β-HMX crystal, the cocrystal, and mixture models were acquired and compared. The results manifest that temperature has an important impact on the binding capacity between the components of the cocrystal and mixture. The binding energies decrease as the temperature rises, and the cocrystal has larger values than those of mixture. For all the models, the [Formula: see text] increases and the CEDs decrease with the rising temperature, implying that the sensitivity of the explosives increases, while the [Formula: see text] values of the cocrystal are smaller than those of HMX and the CED values are between those of RDX and β-HMX, indicating that the sensitivity has been enhanced through co-crystallization. As the temperature increases, the shear modulus (G), bulk modulus (K), and tensile modulus (E) values of all models have an evident downtrend. Simultaneously, G, K, and E values of the cocrystal model are less than those of RDX and β-HMX, while the K/G ratio and Cauchy pressure (C₁₂-C₄₄) are larger, signifying that co-crystallization can weaken the brittleness and enhance the ductility of the pure crystals. Compared with the mixture, the cocrystal has better ductility and stability.

Theoretical calculation into the effect of molar ratio on the structures, stability, mechanical properties and detonation performance of 1,3,5,7-tetranitro-1,3,5,7-tetrazocane/ 1,3,5-trinitro-1,3,5-triazacyco-hexane cocrystal

Molecular dynamics (MD) simulation was conducted to research the effect of molar ratio on the thermal stability, mechanical properties, and detonation performance of HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane)/RDX (1,3,5-trinitro-1,3,5-triazacyco-hexane) cocrystal explosive at ambient condition. The binding energy, mechanical properties, and the detonation parameters of the pure β-HMX, RDX crystal, and the cocrystal models were got and contrasted. The results demonstrate that molar ratio has a great influence on the properties of the cocrystal system. The binding energy of the cocrystals has the maximum values at the 1:1 molar ratio, indicating that the stability of HMX/RDX(1:1) cocrystal is the best and HMX and RDX may prefer to cocrystallizing at 1:1 molar ratio. What's more, the tensile modulus (E) and shear modulus (G) of the HMX/RDX(1:1) cocrystals have the minimum value, while the C₁₂-C₄₄ and K/G have the maximum value, implying that the cocrystal at 1:1 molar ratio has the best mechanical properties. Simultaneously, the E, K, and G of the cocrystals are all smaller than those of β-HMX's and generally larger than those RDX's, while the Cauchy pressure (C₁₂-C₄₄) and K/G ratio were greater, demonstrating that cocrystallizing can improve the brittleness and enhance the ductility. The detonation velocity (D) and detonation pressure (P) decrease with the rising RDX content, while the properties are still superior to the pure RDX crystal; thus, the energy properties of the cocrystal are still excellent. In a word, HMX/RDX cocrystal at 1:1 molar ratio has the best thermal stability, mechanical properties, and the excellent energetic performance.

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.

Comparative studies on structure, sensitivity and mechanical properties of CL-20/DNDAP cocrystal and composite by molecular dynamics simulation

Molecular dynamics simulation was performed on 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20), 2,4-dinitro-2,4-diazapentane (DNDAP), and CL-20/DNDAP cocrystal and composite under COMPASS force field at different temperatures. The binding energy (E_bind), radial distribution function (RDF), trigger bond length, cohesive energy density (CED) and mechanical properties were studied and compared. The results show that the binding energy of the cocrystal is evidently higher than that of the composite at the same temperature. RDF analysis reveals that hydrogen bonds and vdW forces between CL-20 and DNDAP exist in both CL-20/DNDAP cocrystal and composite, and the interactions in the cocrystal are stronger than those in the composite. The maximum trigger bond length decreases in the order ε-CL-20 > CL-20/DNDAP composite > CL-20/DNDAP cocrystal. Moreover, the rigidity and stiffness of the cocrystal and composite decrease compared to that of CL-20, while the ductility and elasticity are better than that of the two pure components. These results demonstrate that CL-20/DNDAP cocrystal might be very promising in explosive applications.

This work proves that cocrystallization is a more effective modification method than mixing and elucidated the underlying mechanism.

Theoretical investigation of the structures and properties of CL-20/DNB cocrystal and associated PBXs by molecular dynamics simulation

In this work, a CL-20/DNB cocrystal explosive model was established and six different kinds of fluoropolymers, i.e., PVDF, PCTFE, F₂₃₁₁, F₂₃₁₂, F₂₃₁₃ and F₂₃₁₄ were added into the (1 0 0), (0 1 0), (0 0 1) crystal orientations to obtain the corresponding polymer bonded explosives (PBXs). The influence of fluoropolymers on PBX properties (energetic property, stability and mechanical properties) was investigated and evaluated using molecular dynamics (MD) methods. The results reveal a decrease in engineering moduli, an increase in Cauchy pressure (i.e., rigidity and stiffness is lessened), and an increase in plastic properties and ductility, thus indicating that the fluoropolymers have a beneficial influence on the mechanical properties of PBXs. Of all the PBXs models tested, the mechanical properties of CL-20/DNB/F₂₃₁₁ were the best. Binding energies show that CL-20/DNB/F₂₃₁₁ has the highest intermolecular interaction energy and best compatibility and stability. Therefore, F₂₃₁₁ is the most suitable fluoropolymer for PBXs. The mechanical properties and binding energies of the three crystal orientations vary in the order (0 1 0) > (0 0 1) > (1 0 0), i.e., the mechanical properties of the (0 1 0) crystal orientation are best, and this is the most stable crystal orientation. Detonation performance results show that the density and detonation parameters of PBXs are lower than those of the pure CL-20 and CL-20/DNB cocrystal explosive. The power and energetic performance of PBXs are thus weakened; however, these PBXs still have excellent detonation performance and are very promising. The results and conclusions provide some helpful guidance and novel instructions for the design and manufacture of PBXs.

Comparative studies on structures, mechanical properties, sensitivity, stabilities and detonation performance of CL-20/TNT cocrystal and composite explosives by molecular dynamics simulation

To investigate and compare the differences of structures and properties of CL-20/TNT cocrystal and composite explosives, the CL-20/TNT cocrystal and composite models were established. Molecular dynamics simulations were performed to investigate the structures, mechanical properties, sensitivity, stabilities and detonation performance of cocrystal and composite models with COMPASS force field in NPT ensemble. The lattice parameters, mechanical properties, binding energies, interaction energy of trigger bond, cohesive energy density and detonation parameters were determined and compared. The results show that, compared with pure CL-20, the rigidity and stiffness of cocrystal and composite models decreased, while plastic properties and ductility increased, so mechanical properties can be effectively improved by adding TNT into CL-20 and the cocrystal model has better mechanical properties. The interaction energy of the trigger bond and the cohesive energy density is in the order of CL-20/TNT cocrystal > CL-20/TNT composite > pure CL-20, i.e., cocrystal model is less sensitive than CL-20 and the composite model, and has the best safety parameters. Binding energies show that the cocrystal model has higher intermolecular interaction energy values than the composite model, thus illustrating the better stability of the cocrystal model. Detonation parameters vary as CL-20 > cocrystal > composite, namely, the energy density and power of cocrystal and composite model are weakened; however, the CL-20/TNT cocrystal explosive still has desirable energy density and detonation performance. This results presented in this paper help offer some helpful guidance to better understand the mechanism of CL-20/TNT cocrystal explosives and provide some theoretical assistance for cocrystal explosive design.