Theoretical investigations on the structures and properties of CL-20/TNT cocrystal and its defective models by molecular dynamics simulation

Crystal morphology prediction of CL-20 and 1,4-DNI co-crystal at different temperatures

CONTEXT

The morphologies of hexanitrohexaazaisowurtzitane (CL-20) and 1,4-dinitroimidazole (1,4-DNI) co-crystal under vacuum or solvent at different temperatures were predicted. The CL-20/1,4-DNI co-crystal has six important growth crystal planes: (002), (011), (101), (11‒1), (110), (111). The areas of (002), (101), and (011) planes account for a relatively large proportion, which are important crystal planes that affect the crystal morphology. The crystal habits at different temperatures were simulated. The simulation results showed that the crystal plane attachment energy of CL-20 and 1,4-DNI co-crystal increases with the increase of temperature, indicating that the increase of temperature is conducive to the growth of crystal planes. The aspect ratio decreases with the increase of temperature and the morphology of co-crystal becomes more spherical at a higher temperature. The theoretical predictions are in good agreement with the experiment. The simulation results can provide guidance for the crystallization of CL-20/1,4-DNI to obtain a nearly spherical crystal morphology.

METHODS

The CL-20/1,4-DNI unit cell structure was geometrically optimized by the COMPASS force field. The AE model was used to predict the morphology of CL-20/1,4-DNI under vacuum, resulting in the most morphologically important growth planes. Ethyl acetate was selected as the solvent. The interaction energy between the solvent and the crystal plane, and the attachment energies in solvent at 298 K, 320 K, 340 K, 360 K, and 380 K were predicted. The NVT ensemble is used in the molecular dynamics calculation process. The simulation step is 1 fs and the total simulation time is 500 ps. The Andersen thermostat is selected as the temperature control method. In the potential energy calculation, the atom-based and Ewald methods were selected to calculate the van der Waals force and the electrostatic interaction force, respectively.

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.