Theoretical insights into effects of molar ratios on stabilities, mechanical properties and detonation performance of CL-20/RDX cocrystal explosives by molecular dynamics simulation

Insight into Formation, Synchronized Release and Stability of Co-Amorphous Curcumin-Piperine by Integrating Experimental-Modeling Techniques

Intermolecular interactions between drug and co-former are crucial in the formation, release and physical stability of co-amorphous system. However, the interactions remain difficult to investigate with only experimental tools. In this study, intermolecular interactions of co-amorphous curcumin-piperine (i.e., CUR-PIP CM) during formation, dissolution and storage were explored by integrating experimental and modeling techniques. The formed CUR-PIP CM exhibited the strong hydrogen bond interaction between the phenolic OH group of CUR and the CO group of PIP as confirmed by FTIR, ss 13C NMR and molecular dynamics (MD) simulation. In comparison to crystalline CUR, crystalline PIP and their physical mixture, CUR-PIP CM performed significantly increased dissolution accompanied by the synchronized release of CUR and PIP, which arose from the greater interaction energy of H2O-CUR molecules and H2O-PIP molecules than CUR-PIP molecules, breaking the hydrogen bond between CUR and PIP molecules, and then causing a pair-wise solvation of CUR-PIP CM at the molecular level. Furthermore, the stronger intermolecular interaction between CUR and PIP was revealed by higher binding energy of CUR-PIP molecules, which contributed to the excellent physical stability of CUR-PIP CM over amorphous CUR or PIP. The study provides a unique insight into the formation, release and stability of co-amorphous system from MD perspective. Meanwhile, this integrated technique can be used as a practical methodology for the future design of co-amorphous formulations.

Theoretical study on the structure, electronic properties, intermolecular interactions, and detonation performance of DAF:ADNP co-crystal

CONTEXT

Explosives have a wide range of applications in many fields due to their high energy and high density. Recently, a new synthesized co-crystal explosive DAF:ADNP presents high detonation performance and low sensitivity. This work is aimed to understand how the structure and intermolecular interactions affect the performance of the DAF:ADNP co-crystal. The results indicate that the formed π-π interactions and stronger hydrogen bonds in the co-crystal enhance its stability and its impact sensitivity is reduced. The strong intralayer H···N and H···O interactions and interlayer π-π stacking are the main driving force for the formation of the co-crystal. Compared with the pure crystals, the detonation performance of the co-crystal slightly decreases, while its sensitivity reduces.

METHODS

All calculations were used the DFT-PBE-D method with Vanderbilt-type ultrasoft pseudopotentials and plane wave (340.0 eV) in the CASTEP package. Radial distribution function were calculated by NVT-MD simulations for 100 ps with a time step of 1 fs at 298 K. Hirshfeld surfaces were generated by CrystalExplorer 3.0 and reduced density gradient analyses were performed by Multiwfn 3.0.

Energy and Electronic Properties of Nanostructures Based on the CL-20 Framework with the Replacement of the Carbon Atoms by Silicon and Germanium: A Density Functional Theory Study

We consider Si_nCL-20 and Ge_nCL-20 systems with carbon atoms replaced by silicon/germanium atoms and their dimers. The physicochemical properties of the silicon/germanium analogs of the high-energy molecule CL-20 and its dimers were determined and studied using density functional theory with the B3LYP/6-311G(d,p) level of theory. It was found that the structure and geometry of Si_nCL-20/Ge_nCL-20 molecules change dramatically with the appearance of Si-/Ge-atoms. The main difference between silicon- or germanium-substituted Si_nCL-20/Ge_nCL-20 molecules and the pure CL-20 molecule is that the NO₂ functional groups make a significant rotation relative to the starting position in the classical molecule, and the effective diameter of the frame of the systems increases with the addition of Si-/Ge-atoms. Thus, the effective framework diameter of a pure CL-20 molecule is 3.208 Å, while the effective diameter of a fully silicon-substituted Si₆CL-20 molecule is 4.125 Å, and this parameter for a fully germanium-substituted Ge₆CL-20 molecule is 4.357 Å. The addition of silicon/germanium atoms to the system leads to a decrease in the binding energy. In detail, the binding energies for CL-20/Si₆CL-20/Ge₆CL-20 molecules are 4.026, 3.699, 3.426 eV/atom, respectively. However, it has been established that the framework maintains stability, with an increase in the number of substituting silicon or germanium atoms. In addition, we designed homodesmotic reactions for the CL-20 molecule and its substituted derivatives Si₆CL-20/Ge₆CL-20, and then determined the strain energy to find out in which case more energy would be released when the framework breaks. Further, we also studied the electronic properties of systems based on CL-20 molecules. It was found that the addition of germanium or silicon atoms instead of carbon leads to a decrease in the size of the HOMO-LUMO gap. Thus, the HOMO-LUMO gaps of the CL-20/Si₆CL-20/Ge₆CL-20 molecules are 5.693, 5.339, and 5.427 eV, respectively. A similar dependence is also observed for CL-20 dimers. So, in this work, we have described in detail the dependence of the physicochemical parameters of CL-20 molecules and their dimers on the types of atoms upon substitution.

Study on the effect of solvent on cocrystallization of CL-20 and HMX through theoretical calculations and experiments

Cocrystallization is a helpful method for explosives design. However, lack of understanding of the cocrystallization mechanism leads to inefficiency in cocrystal preparation. Therefore, studying the effects of solvent on cocrystal is of great importance for the efficient application of 2,4,6,8,10,12-hexanitrohexaazaisowurtzitane (CL-20). In this paper, the effect of solvent on cocrystallization is investigated by the CL-20/HMX cocrystal/solvent cluster model, the CL-20/HMX/solvent mixture model, the CL-20/HMX cocrystal/solvent interface model combined with quantum chemistry and molecular dynamic methods. The authors find that the hydrogen bond between CL-20 and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) is the strongest and the binding energy of cocrystal and solvent molecules is the weakest in ethyl acetate (EA) solvent, indicating that CL-20 and HMX tend to be combined together and there is less hindrance by solvent molecules. Analysis of the CL-20/HMX/solvent mixture and mass density distribution studies show that the solvent effect has a great influence on the crystal faces and the cocrystallization rate of CL-20 and HMX is the highest in EA solvent. The XRD and SEM characterization results are consistent with the theoretical calculations. The present work on the effects of solvent on CL-20/HMX cocrystals is beneficial for understanding the mechanism of the growth of energetic cocrystal materials. It is helpful in selecting more suitable theoretical and experimental conditions and makes access to excellent cocrystals more efficient.

Cocrystallization is a helpful method for explosives design. Studying the effects of solvent on cocrystal is of great importance for the efficient application of CL-20/HMX cocrystal.

Investigation of the solid-liquid ternary phase diagrams of 2HNIW·HMX cocrystal

The influence of temperature and solvent on the solid–liquid ternary phase diagrams of the 2HNIW·HMX cocrystal has been investigated. Ternary phase diagrams were constructed for the 2HNIW·HMX cocrystal in acetonitrile and ethyl acetate at 15 °C and 25 °C. HMX and HNIW showed inconsistent dissolution behavior and congruent dissolution behavior in acetonitrile and ethyl acetate, respectively. In the HMX–HNIW–acetonitrile system, the 2HNIW·HMX cocrystal has a narrow thermodynamically stable region at both temperatures. The cocrystal exhibits a wider thermodynamically stable region in the HMX–HNIW–ethyl acetate system. The results show that the choice of solvent has a crucial influence on the dissolution behavior of the cocrystal and the size and position of each region in the phase diagram, while the temperature has no apparent effect on the overall appearance of the phase diagram. By properly selecting the ratios, the 2HNIW·HMX cocrystal could be prepared by the isothermal slurry conversion crystallization method.

The ternary phase diagrams of 2HNIW·HMX cocrystal system in ethyl acetate at 15 °C. In the HMX–HNIW–ethyl acetate system, HMX and HNIW showed congruent dissolution behavior, and the cocrystal exhibits a wider thermodynamically stable region.

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.

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.

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.

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.