CL-20-Based Cocrystal Energetic Materials: Simulation, Preparation and Performance

The cocrystallization of high-energy explosives has attracted great interests since it can alleviate to a certain extent the power-safety contradiction. 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaaza-isowurtzitane (CL-20), one of the most powerful explosives, has attracted much attention for researchers worldwide. However, the disadvantage of CL-20 has increased sensitivity to mechanical stimuli and cocrystallization of CL-20 with other compounds may provide a way to decrease its sensitivity. The intermolecular interaction of five types of CL-20-based cocrystal (CL-20/TNT, CL-20/HMX, CL-20/FOX-7, CL-20/TKX-50 and CL-20/DNB) by using molecular dynamic simulation was reviewed. The preparation methods and thermal decomposition properties of CL-20-based cocrystal are emphatically analyzed. Special emphasis is focused on the improved mechanical performances of CL-20-based cocrystal, which are compared with those of CL-20. The existing problems and challenges for the future work on CL-20-based cocrystal are discussed.


Introduction
Energetic materials (EMs) are widely used in military and civilian applications, such as weaponry, aerospace explorations and fireworks. However, both the power and safety of energetic materials are the most concern in the field of energetic materials application, but there is an essential contradiction between them: the highly energetic materials are often not safe, and, at present, the rareness of pure low-sensitive and highly energetic explosive has been found [1 -3]. For a long time, in order to obtain EMs with lower sensitivity, the modifications of existing explosives have often focused mainly on recrystallizing with solution and coating with polymer [4][5][6]. However, these traditional methods cannot markedly reduce the sensitivities of existing explosives by only modifying morphology or diluting power, due to the unchanging inherent structures of explosive molecules [7,8]. Due to the stringent requirements for both low-sensitivity and high-power simultaneously, the cocrystallization of explosive, a technique by which a multi-component crystal of several neutral explosive molecules forms in a defined ratio through non-covalent interactions (e.g., H-bond, electrostatic interaction, etc.) [9,10], has attracted great interest since it can alleviate, to a certain extent, the power-safety contradiction [11,12]. The new energetic cocrystals can potentially exhibit decreased sensitivity, higher performance through Simulation details: Andersen was set as the temperature control method (298.15 K). COMPASS force-field was assigned. Summation methods for electrostatic and van der Waals were Ewald and Atom-based, respectively. The accuracy for the Ewald method was 1.0 × 10 −4 kcal·mol −1 . Cutoff distance and buffer width for the Atom-based method were 15.5 Å and 2.0 Å, respectively; 1.0 f. of time step was set for MD processes, and the total dynamic time was performed with 100,000 fs. All the MD calculations were carried out with MS 7.0 [31].
Because the binding energies of cocrystals obtained from the MD calculations are of different numbers of molecules, i.e., different supercells and different molecular molar ratios, those of the same number of molecules should be adopted as a standard for assessing the component interaction strengths and stabilities of cocrystals. Therefore, according to our recent investigation [32], an energy correction formula for binding energy was used to standardize (uniform) the differences caused by diverse supercells and different molecular molar ratios as follows: where Eb * denotes the binding energy after corrected, and Ni and N0 are the number of molecules for different supercells and a standard pattern (molar ratio in 1:1), respectively. The binding energy Eb is calculated by the following formula: where Etot, ECL-20 or ECM is single point energy of cocrystal or monomer; n and m are the number of monomers in cocrystal.

CL-20/TNT System
CL-20 is a complex cage compound, which can be regarded as a six-membered ring and two five-membered rings arranged together, and it has six nitro groups. The preparation of CL-20/TNT cocrystal by the solvent method has attracted much attention worldwide due to its high-energy, insensitive and cheap explosives. The density and melting point of the CL-20/TNT cocrystal explosive (1.91 cm, 136 °C) are between CL-20 (2.04 g/cm 3 , 210 °C) and TNT (1.70 g/cm 3 , 81 °C), and the sensitivity is nearly double reduced that of CL-20 [11]. Figure 1 shows the intermolecular hydrogen bond of CL-20/TNT cocrystal. Simulation details: Andersen was set as the temperature control method (298.15 K). COMPASS force-field was assigned. Summation methods for electrostatic and van der Waals were Ewald and Atom-based, respectively. The accuracy for the Ewald method was 1.0 × 10 −4 kcal·mol −1 . Cutoff distance and buffer width for the Atom-based method were 15.5 Å and 2.0 Å, respectively; 1.0 f. of time step was set for MD processes, and the total dynamic time was performed with 100,000 fs. All the MD calculations were carried out with MS 7.0 [31].
Because the binding energies of cocrystals obtained from the MD calculations are of different numbers of molecules, i.e., different supercells and different molecular molar ratios, those of the same number of molecules should be adopted as a standard for assessing the component interaction strengths and stabilities of cocrystals. Therefore, according to our recent investigation [32], an energy correction formula for binding energy was used to standardize (uniform) the differences caused by diverse supercells and different molecular molar ratios as follows: where Eb * denotes the binding energy after corrected, and Ni and N0 are the number of molecules for different supercells and a standard pattern (molar ratio in 1:1), respectively. The binding energy Eb is calculated by the following formula: where Etot, ECL-20 or ECM is single point energy of cocrystal or monomer; n and m are the number of monomers in cocrystal.

CL-20/TNT System
CL-20 is a complex cage compound, which can be regarded as a six-membered ring and two five-membered rings arranged together, and it has six nitro groups. The preparation of CL-20/TNT cocrystal by the solvent method has attracted much attention worldwide due to its high-energy, insensitive and cheap explosives. The density and melting point of the CL-20/TNT cocrystal explosive (1.91 cm, 136 °C) are between CL-20 (2.04 g/cm 3 , 210 °C) and TNT (1.70 g/cm 3 , 81 °C), and the sensitivity is nearly double reduced that of CL-20 [11]. Figure 1 shows the intermolecular hydrogen bond of CL-20/TNT cocrystal. Simulation details: Andersen was set as the temperature control method (298.15 K). COMPASS force-field was assigned. Summation methods for electrostatic and van der Waals were Ewald and Atom-based, respectively. The accuracy for the Ewald method was 1.0 × 10 −4 kcal·mol −1 . Cutoff distance and buffer width for the Atom-based method were 15.5 Å and 2.0 Å, respectively; 1.0 f. of time step was set for MD processes, and the total dynamic time was performed with 100,000 fs. All the MD calculations were carried out with MS 7.0 [31].
Because the binding energies of cocrystals obtained from the MD calculations are of different numbers of molecules, i.e., different supercells and different molecular molar ratios, those of the same number of molecules should be adopted as a standard for assessing the component interaction strengths and stabilities of cocrystals. Therefore, according to our recent investigation [32], an energy correction formula for binding energy was used to standardize (uniform) the differences caused by diverse supercells and different molecular molar ratios as follows: where Eb * denotes the binding energy after corrected, and Ni and N0 are the number of molecules for different supercells and a standard pattern (molar ratio in 1:1), respectively. The binding energy Eb is calculated by the following formula: where Etot, ECL-20 or ECM is single point energy of cocrystal or monomer; n and m are the number of monomers in cocrystal.

CL-20/TNT System
CL-20 is a complex cage compound, which can be regarded as a six-membered ring and two five-membered rings arranged together, and it has six nitro groups. The preparation of CL-20/TNT cocrystal by the solvent method has attracted much attention worldwide due to its high-energy, insensitive and cheap explosives. The density and melting point of the CL-20/TNT cocrystal explosive (1.91 cm, 136 °C) are between CL-20 (2.04 g/cm 3 , 210 °C) and TNT (1.70 g/cm 3 , 81 °C), and the sensitivity is nearly double reduced that of CL-20 [11]. Figure 1 shows the intermolecular hydrogen bond of CL-20/TNT cocrystal.

FOX-7
Molecules 2020, 25 Simulation details: Andersen was set as the temperature control method (298.15 K). COMPASS force-field was assigned. Summation methods for electrostatic and van der Waals were Ewald and Atom-based, respectively. The accuracy for the Ewald method was 1.0 × 10 −4 kcal·mol −1 . Cutoff distance and buffer width for the Atom-based method were 15.5 Å and 2.0 Å, respectively; 1.0 f. of time step was set for MD processes, and the total dynamic time was performed with 100,000 fs. All the MD calculations were carried out with MS 7.0 [31].
Because the binding energies of cocrystals obtained from the MD calculations are of different numbers of molecules, i.e., different supercells and different molecular molar ratios, those of the same number of molecules should be adopted as a standard for assessing the component interaction strengths and stabilities of cocrystals. Therefore, according to our recent investigation [32], an energy correction formula for binding energy was used to standardize (uniform) the differences caused by diverse supercells and different molecular molar ratios as follows: where Eb * denotes the binding energy after corrected, and Ni and N0 are the number of molecules for different supercells and a standard pattern (molar ratio in 1:1), respectively. The binding energy Eb is calculated by the following formula: where Etot, ECL-20 or ECM is single point energy of cocrystal or monomer; n and m are the number of monomers in cocrystal.

CL-20/TNT System
CL-20 is a complex cage compound, which can be regarded as a six-membered ring and two five-membered rings arranged together, and it has six nitro groups. The preparation of CL-20/TNT cocrystal by the solvent method has attracted much attention worldwide due to its high-energy, insensitive and cheap explosives. The density and melting point of the CL-20/TNT cocrystal explosive (1.91 cm, 136 °C) are between CL-20 (2.04 g/cm 3 , 210 °C) and TNT (1.70 g/cm 3 , 81 °C), and the sensitivity is nearly double reduced that of CL-20 [11]. Figure 1 shows the intermolecular hydrogen bond of CL-20/TNT cocrystal.

TKX-50
Molecules 2020, 25 Simulation details: Andersen was set as the temperature control method (298.15 K). COMPASS force-field was assigned. Summation methods for electrostatic and van der Waals were Ewald and Atom-based, respectively. The accuracy for the Ewald method was 1.0 × 10 −4 kcal·mol −1 . Cutoff distance and buffer width for the Atom-based method were 15.5 Å and 2.0 Å, respectively; 1.0 f. of time step was set for MD processes, and the total dynamic time was performed with 100,000 fs. All the MD calculations were carried out with MS 7.0 [31].
Because the binding energies of cocrystals obtained from the MD calculations are of different numbers of molecules, i.e., different supercells and different molecular molar ratios, those of the same number of molecules should be adopted as a standard for assessing the component interaction strengths and stabilities of cocrystals. Therefore, according to our recent investigation [32], an energy correction formula for binding energy was used to standardize (uniform) the differences caused by diverse supercells and different molecular molar ratios as follows: where Eb * denotes the binding energy after corrected, and Ni and N0 are the number of molecules for different supercells and a standard pattern (molar ratio in 1:1), respectively. The binding energy Eb is calculated by the following formula: where Etot, ECL-20 or ECM is single point energy of cocrystal or monomer; n and m are the number of monomers in cocrystal.

CL-20/TNT System
CL-20 is a complex cage compound, which can be regarded as a six-membered ring and two five-membered rings arranged together, and it has six nitro groups. The preparation of CL-20/TNT cocrystal by the solvent method has attracted much attention worldwide due to its high-energy, insensitive and cheap explosives. The density and melting point of the CL-20/TNT cocrystal explosive (1.91 cm, 136 °C) are between CL-20 (2.04 g/cm 3 , 210 °C) and TNT (1.70 g/cm 3 , 81 °C), and the sensitivity is nearly double reduced that of CL-20 [11]. Figure 1 shows the intermolecular hydrogen bond of CL-20/TNT cocrystal. Simulation details: Andersen was set as the temperature control method (298.15 K). COMPASS force-field was assigned. Summation methods for electrostatic and van der Waals were Ewald and Atom-based, respectively. The accuracy for the Ewald method was 1.0 × 10 −4 kcal·mol −1 . Cutoff distance and buffer width for the Atom-based method were 15.5 Å and 2.0 Å, respectively; 1.0 f. of time step was set for MD processes, and the total dynamic time was performed with 100,000 fs. All the MD calculations were carried out with MS 7.0 [31].
Because the binding energies of cocrystals obtained from the MD calculations are of different numbers of molecules, i.e., different supercells and different molecular molar ratios, those of the same number of molecules should be adopted as a standard for assessing the component interaction strengths and stabilities of cocrystals. Therefore, according to our recent investigation [32], an energy correction formula for binding energy was used to standardize (uniform) the differences caused by diverse supercells and different molecular molar ratios as follows: where Eb * denotes the binding energy after corrected, and Ni and N0 are the number of molecules for different supercells and a standard pattern (molar ratio in 1:1), respectively. The binding energy Eb is calculated by the following formula: where Etot, ECL-20 or ECM is single point energy of cocrystal or monomer; n and m are the number of monomers in cocrystal.

CL-20/TNT System
CL-20 is a complex cage compound, which can be regarded as a six-membered ring and two five-membered rings arranged together, and it has six nitro groups. The preparation of CL-20/TNT cocrystal by the solvent method has attracted much attention worldwide due to its high-energy, insensitive and cheap explosives. The density and melting point of the CL-20/TNT cocrystal explosive (1.91 cm, 136 °C) are between CL-20 (2.04 g/cm 3 , 210 °C) and TNT (1.70 g/cm 3 , 81 °C), and the sensitivity is nearly double reduced that of CL-20 [11]. Figure 1 shows the intermolecular hydrogen bond of CL-20/TNT cocrystal.
Simulation details: Andersen was set as the temperature control method (298.15 K). COMPASS force-field was assigned. Summation methods for electrostatic and van der Waals were Ewald and Atom-based, respectively. The accuracy for the Ewald method was 1.0 × 10 −4 kcal·mol −1 . Cutoff distance and buffer width for the Atom-based method were 15.5 Å and 2.0 Å, respectively; 1.0 f. of time step was set for MD processes, and the total dynamic time was performed with 100,000 fs. All the MD calculations were carried out with MS 7.0 [31].
Because the binding energies of cocrystals obtained from the MD calculations are of different numbers of molecules, i.e., different supercells and different molecular molar ratios, those of the same number of molecules should be adopted as a standard for assessing the component interaction strengths and stabilities of cocrystals. Therefore, according to our recent investigation [32], an energy correction formula for binding energy was used to standardize (uniform) the differences caused by diverse supercells and different molecular molar ratios as follows: where E b * denotes the binding energy after corrected, and N i and N 0 are the number of molecules for different supercells and a standard pattern (molar ratio in 1:1), respectively. The binding energy E b is calculated by the following formula: where E tot , E CL-20 or E CM is single point energy of cocrystal or monomer; n and m are the number of monomers in cocrystal.

CL-20/TNT System
CL-20 is a complex cage compound, which can be regarded as a six-membered ring and two five-membered rings arranged together, and it has six nitro groups. The preparation of CL-20/TNT cocrystal by the solvent method has attracted much attention worldwide due to its high-energy, insensitive and cheap explosives. The density and melting point of the CL-20/TNT cocrystal explosive (1.91 cm, 136 • C) are between CL-20 (2.04 g/cm 3 , 210 • C) and TNT (1.70 g/cm 3 , 81 • C), and the sensitivity is nearly double reduced that of CL-20 [11]. Figure 1 shows the intermolecular hydrogen bond of CL-20/TNT cocrystal.  (13) in TNT are 0.231 nm and 0.244 nm, respectively, which are smaller than the sum of their van der Waals radii (0.272 nm [33]), indicating the existence of intermolecular hydrogen bonds. Through the interaction of C-H…O hydrogen bond between the intermolecular, CL-20 and TNT molecules are connected as the form of "zigzag chain" structure, and combined to form the stable structure of CL-20/TNT cocrystal, indicating that the intermolecular hydrogen bond makes the arrangement of CL-20 and TNT molecules in the crystal more regular, and the molecular packing is close with higher crystal density (1.84 g·cm −3 ).
Meanwhile, the cohesive energy density (CED) of CL-20/TNT cocrystal and its CL-20/TNT composite mixture at different temperatures were investigated [11], results show in Table 2 and Figure 2. With the increase of temperature, CED and its component VDW force and electrostatic force decrease gradually, which is negatively related to the law that the sensitivity increases gradually with the increase of temperature. Under certain conditions, CED, which reflects the energy required for phase transition, can be used to determine the relative thermal sensitivity. Deviations are listed in the parentheses.   (13) in TNT are 0.231 nm and 0.244 nm, respectively, which are smaller than the sum of their van der Waals radii (0.272 nm [33]), indicating the existence of intermolecular hydrogen bonds. Through the interaction of C-H . . . O hydrogen bond between the intermolecular, CL-20 and TNT molecules are connected as the form of "zigzag chain" structure, and combined to form the stable structure of CL-20/TNT cocrystal, indicating that the intermolecular hydrogen bond makes the arrangement of CL-20 and TNT molecules in the crystal more regular, and the molecular packing is close with higher crystal density (1.84 g·cm −3 ).

CL-20/HMX System
Meanwhile, the cohesive energy density (CED) of CL-20/TNT cocrystal and its CL-20/TNT composite mixture at different temperatures were investigated [11], results show in Table 2 and Figure 2. With the increase of temperature, CED and its component VDW force and electrostatic force decrease gradually, which is negatively related to the law that the sensitivity increases gradually with the increase of temperature. Under certain conditions, CED, which reflects the energy required for phase transition, can be used to determine the relative thermal sensitivity. Table 2. The cohesive energy density (CED) of CL-20/TNT cocrystal and its CL-20/TNT composite at different temperatures [11], with permission from Han Neng Cai Liao, 2016.  13) in TNT are 0.231 nm and 0.244 nm, respectively, which are smaller than the sum of their van der Waals radii (0.272 nm [33]), indicating the existence of intermolecular hydrogen bonds. Through the interaction of C-H…O hydrogen bond between the intermolecular, CL-20 and TNT molecules are connected as the form of "zigzag chain" structure, and combined to form the stable structure of CL-20/TNT cocrystal, indicating that the intermolecular hydrogen bond makes the arrangement of CL-20 and TNT molecules in the crystal more regular, and the molecular packing is close with higher crystal density (1.84 g·cm −3 ).

T/K VdW Energy/kJ·cm
Meanwhile, the cohesive energy density (CED) of CL-20/TNT cocrystal and its CL-20/TNT composite mixture at different temperatures were investigated [11], results show in Table 2 and Figure 2. With the increase of temperature, CED and its component VDW force and electrostatic force decrease gradually, which is negatively related to the law that the sensitivity increases gradually with the increase of temperature. Under certain conditions, CED, which reflects the energy required for phase transition, can be used to determine the relative thermal sensitivity. Deviations are listed in the parentheses.

CL-20/HMX System
For HMX, if CL-20 and HMX can form the cocrystal in a certain proportion, the sensitivity of CL-20 can be significantly reduced on the basis of its energy reduce little; on the other hand, the cost of explosive can be significantly reduced and its application range can be expanded [32].
The properties of hexanitrohexaazaisowurtzitane (CL-20)/cyclotetramethylenete-tranitramine (HMX) cocrystal were simulated and compared with those of CL-20/HMX mixture with the molar ratio of CL-20 and HMX as 2:1 [2]. Simulation and calculation results show that the cocrystal process of CL-20/HMX can significantly improve the anti-deformation ability and ductility of the system, and the tensile modulus of the cocrystal structure is greater than that of the mixture.  (Table 3), showing that the two components in the mixture system can be stably adsorbed and have good physical compatibility. The cohesive energy density (CED) value of CL-20/HMX cocrystal structure is far greater than that of CL-20/HMX mixture, and with the increase of temperature, the cohesive energy density of the two components decreases and the structural stability becomes worse ( Table 4). The low sensitivity of CL-20/HMX cocrystal system is caused by the existence of hydrogen bond CH-O with relatively short length.
Five growth faces (0 1 1), (1 01), (1 0 -1), (0 0 2), and (1 1 -1) and random face of ε-, γand β-CL-20 were selected to study the binding energies of the cocrystals with FOX-7 in different molar ratios. These growth faces and a random face of FOX-7were selected to study the binding energies of cocrystals. Based on the MD simulation of substituted models, E b * of the ε-, γ-, and β-CL-20 cocrystal explosives with FOX-7 on the different cocrystal faces in the different molar ratios are calculated ( Table 5). As can be seen, except for the random trend, there is a trend that the strongest binding energies

CL-20/TKX-50 System
TKX-50 is an ionic salt structure with two hydroxylamine cations, and it can easily form hydrogen bonds between H of -NH 3 + in TKX-50 and O of -NO 2 in CL-20, which is in good agreement with the results of surface electrostatic potential energy analysis. The surface electrostatic potential energy analysis of CL-20 and TKX-50 is useful for exploring the formation mechanism of CL-20/TKX-50 cocrystal. CL-20/TKX-50 cocrystals with the mole ratio of 1:1, 1:2, 1:3 and 2:3 were prepared by means of the solvent-nonsolvent method [34]. Compared with the raw materials and the mixture, the XRD spectra of CL-20/TKX-50 cocrystals have no obvious change when the feed ratio is 1:1, 1:3 and 2:3, and it can be considered that the eutectic of CL-20/TKX-50 cocrystals has not been successfully prepared. Finally, the optimal mole ratio is 1:2. Figure 3 shows the surface electrostatic potential energy distribution map of CL-20 (a) and TKX-50 (b), and the equilibrium structure of CL-20/TKX-50 cocrystal model (c), a schematic diagram of intermolecular interaction of CL-20 and TKX-50 model (d).

CL-20/DNB System
As a cheap and insensitive explosive, DNB is often used as an alternative explosive for TNT, but its energy is not ideal. If CL-20 and DNB can be combined in the same lattice through noncovalent bond at the molecular level through cocrystal technology to form an explosive crystal with unique structure, it is expected that the sensitivity and cost of CL-20 will be greatly reduced without significant energy reduction, so as to expand the application range of CL-20 [12]. In order to improve the hazardous performance of CL-20, 1,3-dinitrobenze (DNB) was introduced to CL-20 to form CL-20/DNB cocrystal with 1:1 mole ratio. The primitive cell of CL-20/DNB cocrastal was shown in Figure  4. The molecular dynamics simulation of 2,4,6,8,10,12-hexanitrohexaazaisowurtzitane (CL-20)/1,3dinitrobenzene (DNB) cocrystal and CL-20/DNB cocrystal with two polymer binders, hydroxylterminated polybutatiene (HTPB) and polyethylene glycol (PEG) respectively were conducted [3]. Results indicate that the binding energy of CL-20/DNB/PEG is larger than that of CL-20/DNB/HTPB, predicting that the stability and compatibility of the former is better than those of the latter (Table 6).
In comparison with CL-20/DNB cocrystal, an addition of a small amount of binder (HTPB or PEG) decreases the elastic constants (Cij), tensile modulus (E), bulk (K) and shear modulus (G), while

CL-20/DNB System
As a cheap and insensitive explosive, DNB is often used as an alternative explosive for TNT, but its energy is not ideal. If CL-20 and DNB can be combined in the same lattice through noncovalent bond at the molecular level through cocrystal technology to form an explosive crystal with unique structure, it is expected that the sensitivity and cost of CL-20 will be greatly reduced without significant energy reduction, so as to expand the application range of CL-20 [12]. In order to improve the hazardous performance of CL-20, 1,3-dinitrobenze (DNB) was introduced to CL-20 to form CL-20/DNB cocrystal with 1:1 mole ratio. The primitive cell of CL-20/DNB cocrastal was shown in Figure 4.

CL-20/DNB System
As a cheap and insensitive explosive, DNB is often used as an alternative explosive for TNT, but its energy is not ideal. If CL-20 and DNB can be combined in the same lattice through noncovalent bond at the molecular level through cocrystal technology to form an explosive crystal with unique structure, it is expected that the sensitivity and cost of CL-20 will be greatly reduced without significant energy reduction, so as to expand the application range of CL-20 [12]. In order to improve the hazardous performance of CL-20, 1,3-dinitrobenze (DNB) was introduced to CL-20 to form CL-20/DNB cocrystal with 1:1 mole ratio. The primitive cell of CL-20/DNB cocrastal was shown in Figure  4. The molecular dynamics simulation of 2,4,6,8,10,12-hexanitrohexaazaisowurtzitane (CL-20)/1,3dinitrobenzene (DNB) cocrystal and CL-20/DNB cocrystal with two polymer binders, hydroxylterminated polybutatiene (HTPB) and polyethylene glycol (PEG) respectively were conducted [3]. Results indicate that the binding energy of CL-20/DNB/PEG is larger than that of CL-20/DNB/HTPB, predicting that the stability and compatibility of the former is better than those of the latter (Table 6). In comparison with CL-20/DNB cocrystal, an addition of a small amount of binder (HTPB or PEG) decreases the elastic constants (Cij), tensile modulus (E), bulk (K) and shear modulus (G), while Cauchy pressure (C12-C44) and K/G value increase, showing that the stiffness of the polymerbonded explosives (PBXs) system is weaker, and its ductibility is better. The molecular dynamics simulation of 2,4,6,8,10,12-hexanitrohexaazaisowurtzitane (CL-20)/1,3dinitrobenzene (DNB) cocrystal and CL-20/DNB cocrystal with two polymer binders, hydroxylterminated polybutatiene (HTPB) and polyethylene glycol (PEG) respectively were conducted [3]. Results indicate that the binding energy of CL-20/DNB/PEG is larger than that of CL-20/DNB/HTPB, predicting that the stability and compatibility of the former is better than those of the latter (Table 6). In comparison with CL-20/DNB cocrystal, an addition of a small amount of binder (HTPB or PEG) decreases the elastic constants (C ij ), tensile modulus (E), bulk (K) and shear modulus (G), while Cauchy pressure (C 12 -C 44 ) and K/G value increase, showing that the stiffness of the polymerbonded explosives (PBXs) system is weaker, and its ductibility is better. In addition, the initial thermal decomposition pathways, as well as some important products generating mechanism of CL-20/DNB cocrystal at high temperatures (2000, 2500 K and 3000 K), were studied by reactive molecular dynamics simulations using ReaxFF force field [4]. Results show that with the increasing of temperature during the thermal decomposition process, the time to balance and potential energy decrease, while the quantity of products increases. All the CL-20 molecules decompose faster than that of DNB, and as the temperature rises, the decomposition rate of DNB increases significantly. According to the product identification analysis, the main thermal decomposition products are NO 2 , NO, N 2 , H 2 O, HNO 3 , HON, HONO and CO 2 for the cocrystal. The major initial decomposition mechanism is the breaking of N−NO 2 in the CL-20 and C−NO 2 in the DNB, which contributes to the formation of NO 2 . Then, the number of NO 2 increases to the peak rapidly and decreases subsequently ( Figure 5). After the NO 2 →ONO rearrangement, it participates in other reactions and eventually N 2 , NO, HONO, HON, H 2 O occur, and so on. In addition, the simulation results indicate that carbon-containing clusters formed in the later stage of decomposition at 2500 K and 3000 K, which is a common phenomenon during the detonation of rich carbon-containing explosives.  In addition, the initial thermal decomposition pathways, as well as some important products generating mechanism of CL-20/DNB cocrystal at high temperatures (2000, 2500 K and 3000 K), were studied by reactive molecular dynamics simulations using ReaxFF force field [4]. Results show that with the increasing of temperature during the thermal decomposition process, the time to balance and potential energy decrease, while the quantity of products increases. All the CL-20 molecules decompose faster than that of DNB, and as the temperature rises, the decomposition rate of DNB increases significantly. According to the product identification analysis, the main thermal decomposition products are NO2, NO, N2, H2O, HNO3, HON, HONO and CO2 for the cocrystal. The major initial decomposition mechanism is the breaking of N−NO2 in the CL-20 and C−NO2 in the DNB, which contributes to the formation of NO2. Then, the number of NO2 increases to the peak rapidly and decreases subsequently ( Figure 5). After the NO2→ONO rearrangement, it participates in other reactions and eventually N2, NO, HONO, HON, H2O occur, and so on. In addition, the simulation results indicate that carbon-containing clusters formed in the later stage of decomposition at 2500 K and 3000 K, which is a common phenomenon during the detonation of rich carboncontaining explosives. Meanwhile, the effect of CL-20/DNB cocrystallizing and mixing on the sensitivity, binding energy, mechanical properties and thermal decomposition on the molecular level were investigated under the condition of the COMPASS force field. Results indicate that the cocrystallizing and mixing Meanwhile, the effect of CL-20/DNB cocrystallizing and mixing on the sensitivity, binding energy, mechanical properties and thermal decomposition on the molecular level were investigated under the condition of the COMPASS force field. Results indicate that the cocrystallizing and mixing can reduce the sensitivity of CL-20, increase that of DNB, and the cocrystallizing effect is more obvious and stable [35].

Preparation
For CL-20/TNT cocrystal, the CL-20/TNT cocrystal with 1:1 mole ratio was prepared by means of the recrystallization method at room temperature, and with ethyl acetate as solvent [8]. It was found that the morphology of CL-20/TNT cocrystal is different from that of CL-20 and TNT. The crystal of CL-20 is "spindle" shape and TNT is an irregular block crystal. However, CL-20/TNT cocrystal is a prismatic crystal with a smooth and complete surface and uniform size ( Figure 6). The average particle size is 270 µm, indicating that the designing and preparation of cocrystal can effectively change the shape and size of explosive cocrystals. can reduce the sensitivity of CL-20, increase that of DNB, and the cocrystallizing effect is more obvious and stable [35].

Preparation
For CL-20/TNT cocrystal, the CL-20/TNT cocrystal with 1:1 mole ratio was prepared by means of the recrystallization method at room temperature, and with ethyl acetate as solvent [8]. It was found that the morphology of CL-20/TNT cocrystal is different from that of CL-20 and TNT. The crystal of CL-20 is "spindle" shape and TNT is an irregular block crystal. However, CL-20/TNT cocrystal is a prismatic crystal with a smooth and complete surface and uniform size ( Figure 6). The average particle size is 270 μm, indicating that the designing and preparation of cocrystal can effectively change the shape and size of explosive cocrystals. At the same time, ultrafine CL-20/TNT cocrystal explosive was prepared by a spray-drying method [20]. Results show that the prepared samples are not the mix of CL-20 and TNT but rather ultrafine CL-20/TNT cocrystal explosives. The particle sizes of the cocrystal explosives were lower than 1 μm and they can aggregate into many microparticles, which are spherical in shape and 1-10 μm in size (Figure 7). The thermal decomposition process can be divided into two stages. The peak temperatures of exothermic decomposition for the first and second stages are 218.98 °C and 253.15 °C, respectively ( Figure 8). The characteristic drop height of CL-20/TNT cocrystal explosives is 49.3 cm, which increases by 36.2 cm compared with that of raw CL-20 (Table 7).  At the same time, ultrafine CL-20/TNT cocrystal explosive was prepared by a spray-drying method [20]. Results show that the prepared samples are not the mix of CL-20 and TNT but rather ultrafine CL-20/TNT cocrystal explosives. The particle sizes of the cocrystal explosives were lower than 1 µm and they can aggregate into many microparticles, which are spherical in shape and 1-10 µm in size (Figure 7). The thermal decomposition process can be divided into two stages. The peak temperatures of exothermic decomposition for the first and second stages are 218.98 • C and 253.15 • C, respectively ( Figure 8). The characteristic drop height of CL-20/TNT cocrystal explosives is 49.3 cm, which increases by 36.2 cm compared with that of raw CL-20 (Table 7). can reduce the sensitivity of CL-20, increase that of DNB, and the cocrystallizing effect is more obvious and stable [35].

Preparation
For CL-20/TNT cocrystal, the CL-20/TNT cocrystal with 1:1 mole ratio was prepared by means of the recrystallization method at room temperature, and with ethyl acetate as solvent [8]. It was found that the morphology of CL-20/TNT cocrystal is different from that of CL-20 and TNT. The crystal of CL-20 is "spindle" shape and TNT is an irregular block crystal. However, CL-20/TNT cocrystal is a prismatic crystal with a smooth and complete surface and uniform size ( Figure 6). The average particle size is 270 μm, indicating that the designing and preparation of cocrystal can effectively change the shape and size of explosive cocrystals. At the same time, ultrafine CL-20/TNT cocrystal explosive was prepared by a spray-drying method [20]. Results show that the prepared samples are not the mix of CL-20 and TNT but rather ultrafine CL-20/TNT cocrystal explosives. The particle sizes of the cocrystal explosives were lower than 1 μm and they can aggregate into many microparticles, which are spherical in shape and 1-10 μm in size (Figure 7). The thermal decomposition process can be divided into two stages. The peak temperatures of exothermic decomposition for the first and second stages are 218.98 °C and 253.15 °C, respectively ( Figure 8). The characteristic drop height of CL-20/TNT cocrystal explosives is 49.3 cm, which increases by 36.2 cm compared with that of raw CL-20 (Table 7).   The ultrafine CL-20/HMX cocrystal explosive was prepared by an ultra-highly efficient mixing method [6]. The prepared samples were regular block-like ultrafine CL-20/HMX cocrystal explosives with a uniform particle size of less than 1 μm (Figure 9), which appeared new stronger diffraction peaks at 11.558°, 13.264°, 18.601°, 24.474°, 33.785°, 36.269°. The purity of the CL-20/HMX cocrystal explosive was 92.6%.  The ultrafine CL-20/HMX cocrystal explosive was prepared by an ultra-highly efficient mixing method [6]. The prepared samples were regular block-like ultrafine CL-20/HMX cocrystal explosives with a uniform particle size of less than 1 µm (Figure 9), which appeared new stronger diffraction peaks at 11.558 • , 13.264 • , 18.601 • , 24.474 • , 33.785 • , 36.269 • . The purity of the CL-20/HMX cocrystal explosive was 92.6%.
Moreover, the micro/nano CL-20/HMX energetic cocrystal materials were prepared by mechanical ball milling, and characterized using SEM [7]. Results show that after a milling time of 120 min, CL-20/HMX cocrystals prepared under optimal conditions are spherical in shape and the particle size is 80-250 nm. Compared with respective energetic monomers, the prepared micro/nano CL-20/HMX energetic cocrystal materials exhibit unique crystal structure and thermal decomposition properties ( Figure 10).
In addition, nano-sized CL-20/HMX cocrystals with a mean particle size of 81.6 nm were prepared under the conditions of the 0.3 mm diameter of milling balls [9]. The impact sensitivity was conducted to evaluate the impact safety performance of the CL-20/HMX cocrystal ( Table 7). The micro-morphology of the explosives is near-spherical and the particle size reveals a normal distribution (Figure 11). Before and after milling, the element composition and molecular structure of CL-20/HMX did not change in comparison to the raw CL-20 and HMX, while there are new crystal phases in the formed cocrystal. The characteristic drop height (H 50 ) of CL-20/HMX was 32.62 cm with 5 kg hammer, which is lower than that of raw materials and mixture, indicating that the mechanical ball milling method is not simply to mix the two explosives physically, but to form the cocrystal between the explosive molecules through the hydrogen bond. Meanwhile, the crystal technology can improve the impact safety of explosives.  Moreover, the micro/nano CL-20/HMX energetic cocrystal materials were prepared by mechanical ball milling, and characterized using SEM [7]. Results show that after a milling time of 120 min, CL-20/HMX cocrystals prepared under optimal conditions are spherical in shape and the particle size is 80-250 nm. Compared with respective energetic monomers, the prepared micro/nano CL-20/HMX energetic cocrystal materials exhibit unique crystal structure and thermal decomposition properties ( Figure 10). In addition, nano-sized CL-20/HMX cocrystals with a mean particle size of 81.6 nm were prepared under the conditions of the 0.3 mm diameter of milling balls [9]. The impact sensitivity was conducted to evaluate the impact safety performance of the CL-20/HMX cocrystal ( Table 7). The micro-morphology of the explosives is near-spherical and the particle size reveals a normal distribution (Figure 11). Before and after milling, the element composition and molecular structure of CL-20/HMX did not change in comparison to the raw CL-20 and HMX, while there are new crystal phases in the formed cocrystal. The characteristic drop height (H50) of CL-20/HMX was 32.62 cm with 5 kg hammer, which is lower than that of raw materials and mixture, indicating that the mechanical ball milling method is not simply to mix the two explosives physically, but to form the cocrystal between the explosive molecules through the hydrogen bond. Meanwhile, the crystal technology can improve the impact safety of explosives.   In addition, nano-sized CL-20/HMX cocrystals with a mean particle size of 81.6 nm were prepared under the conditions of the 0.3 mm diameter of milling balls [9]. The impact sensitivity was conducted to evaluate the impact safety performance of the CL-20/HMX cocrystal ( Table 7). The micro-morphology of the explosives is near-spherical and the particle size reveals a normal distribution ( Figure 11). Before and after milling, the element composition and molecular structure of CL-20/HMX did not change in comparison to the raw CL-20 and HMX, while there are new crystal phases in the formed cocrystal. The characteristic drop height (H50) of CL-20/HMX was 32.62 cm with 5 kg hammer, which is lower than that of raw materials and mixture, indicating that the mechanical ball milling method is not simply to mix the two explosives physically, but to form the cocrystal between the explosive molecules through the hydrogen bond. Meanwhile, the crystal technology can improve the impact safety of explosives.  Nano-scale CL-20/HMX cocrystal of CL-20 and HMX in 2:1 molar ratio with a mean size below 200 nm was prepared by bead milling [36]. Figure 12 shows the morphological evolution of the crystal particles. It was found that most of the coformer crystal particles are~1 µm (Figure 12a). The mean particle size of the discrete coformers substantially decreased after 10 min of milling (Figure 12b), with no conversion to the cocrystalline material. Plate-like crystal particles with dimensions less than 500 nm started to appear in the specimen being milled for 20 min, as indicated by the arrow in Figure 12c. The plate-like particles were assigned to the 2CL-20·HMX cocrystal as (1) the 2CL-20·HMX cocrystal is known to have a plate-like morphology; (2) the appearance of these particles and the diffraction peaks of the 2CL-20·HMX cocrystal in the XRD pattern occurred at the same time; and (3) more plate-like particles were observed in the specimens upon further milling (Figure 12d and e, respectively). The observation of these relatively large cocrystal particles seems to be contradicting to the intensive collisions between the grinding media and particles and between the particles themselves occurring during the milling process. It is possible that some growth occurs during the drying of the sampled specimens. diffraction peaks of the 2CL-20·HMX cocrystal in the XRD pattern occurred at the same time; and (3) more plate-like particles were observed in the specimens upon further milling (Figure 12d and e,  respectively). The observation of these relatively large cocrystal particles seems to be contradicting to the intensive collisions between the grinding media and particles and between the particles themselves occurring during the milling process. It is possible that some growth occurs during the drying of the sampled specimens. Preparation procedure: 438.0 mg (1 mmol) CL-20 and 472.3 mg (2 mmol) TKX-50 were dried in a beaker at 40 °C for 3 h, 30 mL DMF was added to the beaker, heat to 80 °C in the water bath, and ultrasonic for 30 min until the drug is completely dissolved. Then, 100 mL chloroform was put into the flask. The completely dissolved solution was dropped into the flask at the rate of 0.8 mL/min, and the magnetic stirring speed was kept at 1000 r/min. After dropping, keep the original speed and continue stirring for 1 h. After standing for 3 h and filtering, the filtered samples were dried in a Preparation procedure: 438.0 mg (1 mmol) CL-20 and 472.3 mg (2 mmol) TKX-50 were dried in a beaker at 40 • C for 3 h, 30 mL DMF was added to the beaker, heat to 80 • C in the water bath, and ultrasonic for 30 min until the drug is completely dissolved. Then, 100 mL chloroform was put into the flask. The completely dissolved solution was dropped into the flask at the rate of 0.8 mL/min, and the magnetic stirring speed was kept at 1000 r/min. After dropping, keep the original speed and continue stirring for 1 h. After standing for 3 h and filtering, the filtered samples were dried in a vacuum for 3 h to obtain CL-20/TKX-50 ultrafine cocrystal samples. The performance was compared with that of CL-20/TKX-50 mechanical mixture in the same mole ratio.
To decrease the sensitivity of CL-20, the CL-20/TKX-50 cocrystal explosive was prepared by solvent-nonsolvent method. The surface electrostatic potentials of CL-20 and TKX-50 were analyzed and the possible non-covalent bonding between cocrystal molecules was predicted [34]. Results shown that the prepared CL-20/TKX-50 cocrystal has a flat sheet shape, the formation, disappearance, shift and change of intensity of peaks been proved the formation of a new lattice structure. The crystal morphology of CL-20 and TKX-50 is nearly spherical shape, and the particle size distribution is relatively uniform, and the particle size is 1 µm. However, the morphology of CL-20/TKX-50 cocrystal shows a long and thin lamellar structure, which is quite different from that of CL-20 and TKX-50, and the particle size of CL-20/TKX-50 cocrystal is 10 µm, indicating that the cocrystal process can not only change the morphology of the original crystal, but also prove the formation of a new crystal ( Figure 13). morphology of CL-20 and TKX-50 is nearly spherical shape, and the particle size distribution is relatively uniform, and the particle size is 1 μm. However, the morphology of CL-20/TKX-50 cocrystal shows a long and thin lamellar structure, which is quite different from that of CL-20 and TKX-50, and the particle size of CL-20/TKX-50 cocrystal is 10 μm, indicating that the cocrystal process can not only change the morphology of the original crystal, but also prove the formation of a new crystal ( Figure  13).

Energetic Performance
In order to study the thermal decomposition properties of CL-20/TNT cocryatal explosives, the samples of CL-20/TNT cocryatals were tested by DSC (TA Instruments, New Castle, DE, USA)at the heating rate of 10 °C·min −1 , the curve shown in Figure 14. As is shown there are three stages of the thermal decomposition process of CL-20/TNT cocrystal, the maximum endothermic peak temperature is 143.5 °C, which is higher than that of the melting point of TNT (81 °C [37]). With the increase of temperature, the hydrogen bond between the cocrystal molecules breaks and the molecular structure is destroyed [38]. There are two processes of exothermic decomposition at 222.6 °C and 250.1 °C, respectively. Compared with the maximum exothermic peak value of CL-20 (321.5 °C [33]) and TNT (245 °C [39]), the exothermic peak of CL-20/TNT cocrystal shifts, indicating that the cocrystal changes the thermal decomposition characteristics of the raw materials and endows the cocrystal with unique thermal decomposition behavior.

Energetic Performance
In order to study the thermal decomposition properties of CL-20/TNT cocryatal explosives, the samples of CL-20/TNT cocryatals were tested by DSC (TA Instruments, New Castle, DE, USA) at the heating rate of 10 • C·min −1 , the curve shown in Figure 14. As is shown there are three stages of the thermal decomposition process of CL-20/TNT cocrystal, the maximum endothermic peak temperature is 143.5 • C, which is higher than that of the melting point of TNT (81 • C [37]). With the increase of temperature, the hydrogen bond between the cocrystal molecules breaks and the molecular structure is destroyed [38]. There are two processes of exothermic decomposition at 222.6 • C and 250.1 • C, respectively. Compared with the maximum exothermic peak value of CL-20 (321.5 • C [33]) and TNT (245 • C [39]), the exothermic peak of CL-20/TNT cocrystal shifts, indicating that the cocrystal changes the thermal decomposition characteristics of the raw materials and endows the cocrystal with unique thermal decomposition behavior. morphology of CL-20 and TKX-50 is nearly spherical shape, and the particle size distribution is relatively uniform, and the particle size is 1 μm. However, the morphology of CL-20/TKX-50 cocrystal shows a long and thin lamellar structure, which is quite different from that of CL-20 and TKX-50, and the particle size of CL-20/TKX-50 cocrystal is 10 μm, indicating that the cocrystal process can not only change the morphology of the original crystal, but also prove the formation of a new crystal ( Figure  13).

Energetic Performance
In order to study the thermal decomposition properties of CL-20/TNT cocryatal explosives, the samples of CL-20/TNT cocryatals were tested by DSC (TA Instruments, New Castle, DE, USA)at the heating rate of 10 °C·min −1 , the curve shown in Figure 14. As is shown there are three stages of the thermal decomposition process of CL-20/TNT cocrystal, the maximum endothermic peak temperature is 143.5 °C, which is higher than that of the melting point of TNT (81 °C [37]). With the increase of temperature, the hydrogen bond between the cocrystal molecules breaks and the molecular structure is destroyed [38]. There are two processes of exothermic decomposition at 222.6 °C and 250.1 °C, respectively. Compared with the maximum exothermic peak value of CL-20 (321.5 °C [33]) and TNT (245 °C [39]), the exothermic peak of CL-20/TNT cocrystal shifts, indicating that the cocrystal changes the thermal decomposition characteristics of the raw materials and endows the cocrystal with unique thermal decomposition behavior.  The impact sensitivity value of CL-20/TNT cocrystal is H 50 = 28 cm. Compared with CL-20 (H 50 = 15 cm), the impact sensitivity of CL-20 is significantly reduced by 87%. Low-sensitivity TNT and high-sensitivity CL-20 combine to form a cocrystal at the molecular scale through eutectic technology, which changes the internal composition and crystal structure of explosives compared with traditional methods. In addition, due to the intermolecular hydrogen bond, on the one hand, it increases the stability of the molecular system of the cocrystal, on the other hand, it improves the anti-vibration of the cocrystal molecule to the mechanical external force, so the desensitization effect is obvious. Through the eutectic technology, it can effectively realize the desensitization of high sensitive explosive and improve its safety performance. The crystal morphology, particle size and sensitivity of prepared CL-20/HMX cocrystal were characterized [6]. The thermal decomposition process of cocrystal explosives had only one exothermic decomposition stage with peak temperatures of 248.3 • C. The enthalpy for the exothermic decomposition of the cocrystal (2912.1 J·g −1 ) was remarkably higher than that of the physical mixture of CL-20 and HMX (1327.3 J·g −1 ) ( Figure 15). The friction sensitivity of CL-20/HMX cocrystal explosive was 84%, which was decreased by 16% compared with original CL-20, and the characteristic height of the cocrystal was increased by 28.6 cm and 11.5 cm compared with original CL-20 and HMX, respectively ( Table 8). The compatibility of CL-20/HMX cocrystal with components of solid propellant shown that the prepared CL-20/HMX cocrystal was compatible with NG/BTTN, AP and Al powder, while incompatible with HGAP, N-100. stability of the molecular system of the cocrystal, on the other hand, it improves the anti-vibration of the cocrystal molecule to the mechanical external force, so the desensitization effect is obvious. Through the eutectic technology, it can effectively realize the desensitization of high sensitive explosive and improve its safety performance. The crystal morphology, particle size and sensitivity of prepared CL-20/HMX cocrystal were characterized [6]. The thermal decomposition process of cocrystal explosives had only one exothermic decomposition stage with peak temperatures of 248.3 °C. The enthalpy for the exothermic decomposition of the cocrystal (2912.1 J·g −1 ) was remarkably higher than that of the physical mixture of CL-20 and HMX (1327.3 J·g −1 ) ( Figure 15). The friction sensitivity of CL-20/HMX cocrystal explosive was 84%, which was decreased by 16% compared with original CL-20, and the characteristic height of the cocrystal was increased by 28.6 cm and 11.5 cm compared with original CL-20 and HMX, respectively ( Table 8). The compatibility of CL-20/HMX cocrystal with components of solid propellant shown that the prepared CL-20/HMX cocrystal was compatible with NG/BTTN, AP and Al powder, while incompatible with HGAP, N-100. Figure 15. DSC curves of different samples [6], with permission from Han Neng Cai Liao, 2020. Figure 15. DSC curves of different samples [6], with permission from Han Neng Cai Liao, 2020. At the same time, the DSC curves of CL-20/HMX cocrystal were prepared by means of mechanical ball milling, and compared with the CL-20/HMX mixture [7]. The mechanical sensitivity of the micro/nano CL-20/HMX energetic cocrystal materials was reduced obviously compared to that of raw HMX, while the energy output property was equivalent to that of raw CL-20.
In order to compare the effect of the prepared CL-20/TKX-50 cocrystal on the thermal decomposition of each component, the DSC curves of CL-20, TKX-50, CL-20/TKX-50 mixture and CL-20/TKX-50 cocrystal were investigated [34] (Figure 16). The exothermic decomposition peaks of CL-20 and TKX-50 were 240.2 and 234.8 • C, respectively, and there were two exothermic decomposition peaks in the decomposition process of the CL-20/TKX-50 mixture, indicating that the decomposition process of the CL-20/TKX-50 mixture is obviously a simple superposition of CL-20 and TKX-50. The first decomposition peak appears at 229.8 • C, corresponding to the exothermic decomposition of part of TKX-50; the second peak appears at 242.5 • C, corresponding to the exothermic decomposition of CL-20. For the CL-20/TKX-50 cocrystal decomposition, with the increase of temperature, the first exothermic decomposition peak of CL-20/TKX-50 cocrystal appears at 171.6 • C, which can be attributed to the destruction of hydrogen in the cocrystal structure and the decomposition of a small amount of TKX-50. Subsequently, a large number of cocrystal products begin to decompose, and the second exothermic decomposition peak appeared at 222.8 • C, indicating the formation of hydrogen bond and the existence of a new structure between CL-20/TKX-50 cocrystal molecules. At the same time, the DSC curves of CL-20/HMX cocrystal were prepared by means of mechanical ball milling, and compared with the CL-20/HMX mixture [7]. The mechanical sensitivity of the micro/nano CL-20/HMX energetic cocrystal materials was reduced obviously compared to that of raw HMX, while the energy output property was equivalent to that of raw CL-20.
In order to compare the effect of the prepared CL-20/TKX-50 cocrystal on the thermal decomposition of each component, the DSC curves of CL-20, TKX-50, CL-20/TKX-50 mixture and CL-20/TKX-50 cocrystal were investigated [34] (Figure 16). The exothermic decomposition peaks of CL-20 and TKX-50 were 240.2 and 234.8 °C, respectively, and there were two exothermic decomposition peaks in the decomposition process of the CL-20/TKX-50 mixture, indicating that the decomposition process of the CL-20/TKX-50 mixture is obviously a simple superposition of CL-20 and TKX-50. The first decomposition peak appears at 229.8 °C, corresponding to the exothermic decomposition of part of TKX-50; the second peak appears at 242.5 °C, corresponding to the exothermic decomposition of CL-20. For the CL-20/TKX-50 cocrystal decomposition, with the increase of temperature, the first exothermic decomposition peak of CL-20/TKX-50 cocrystal appears at 171.6 °C, which can be attributed to the destruction of hydrogen in the cocrystal structure and the decomposition of a small amount of TKX-50. Subsequently, a large number of cocrystal products begin to decompose, and the second exothermic decomposition peak appeared at 222.8 °C, indicating the formation of hydrogen bond and the existence of a new structure between CL-20/TKX-50 cocrystal molecules. In addition, the DSC curves of CL-20, DNB, CL-20/DNB cocrystal were shown in Figure 17. There are three (one is endothermic melting and two exothermic decompositions) thermal decomposition stages. The melting temperature of cocrystal is 136.4 °C, 44.7 °C higher than that of raw material DNB (91.7 °C) in the melting stage, indicating that the cocrystal decomposes into liquid DNB at this point. There are two exothermic peaks at 216.9 °C and 242.9 °C in the exothermic decomposition stage, which was attributed to the exothermic decomposition of two single components after cocrystal decomposition. Moreover, the crystal density of CL-20/DNB is higher than that of CL-20/TNT, and the sensitivity of DNB is lower than that of TNT. Thus, the sensitivity of CL-20/DNB is lower than that of CL-20/TNT. In addition, the cost of the DNB is significantly lower than that of TNT. Therefore, the cocrystal may become an excellent explosive with high-energy, insensitive and low-cost prospective ingredients in explosives and propellants. In addition, the DSC curves of CL-20, DNB, CL-20/DNB cocrystal were shown in Figure 17. There are three (one is endothermic melting and two exothermic decompositions) thermal decomposition stages. The melting temperature of cocrystal is 136.4 • C, 44.7 • C higher than that of raw material DNB (91.7 • C) in the melting stage, indicating that the cocrystal decomposes into liquid DNB at this point. There are two exothermic peaks at 216.9 • C and 242.9 • C in the exothermic decomposition stage, which was attributed to the exothermic decomposition of two single components after cocrystal decomposition. Moreover, the crystal density of CL-20/DNB is higher than that of CL-20/TNT, and the sensitivity of DNB is lower than that of TNT. Thus, the sensitivity of CL-20/DNB is lower than that of CL-20/TNT. In addition, the cost of the DNB is significantly lower than that of TNT. Therefore, the cocrystal may become an excellent explosive with high-energy, insensitive and low-cost prospective ingredients in explosives and propellants.  Table 8. The characteristic drop height of CL-20/TKX-50 cocrystal is lower than that of TKX-50, but significantly higher than that of CL-20, β-HMX and CL-20/TKX-50 mixture, indicating that the impact sensitivity of CL-20/TKX-50 cocrystal is much lower. The detonation velocity and detonation pressure of CL-20/TKX-50 cocrystal are slightly lower than that of CL-20, but the detonation performance is obviously improved compared with β-HMX, indicating CL-20/TKX-50 cocrystal has good detonation performance.

Existing Problems and Challenges
As we all know, it is of great significance for the preparation of energetic materials with highenergy and low-sensitivity to use eutectic technology to modify the single high-energy explosives. While the development of an energetic cocrystal is still in the exploratory stage, the main challenges for cocrystal at present are as follows: (1) the optimal preparation method of energetic cocrystal still needs to be simplified and simulated. The solvent evaporation method is mainly used to explore the preparation conditions of cocrystal by experience at present and many experiments with different stoichiometric ratios cannot meet the engineering requirements. (2) The effective cocrystal formation principle and formation mechanism of energetic compounds under the guidance of thermodynamics are not fully revealed, which greatly hinders the controllable construction of energetic cocrystal compounds [40]. The main challenges for us in the future are possibly: (1) the design and formation mechanism of energetic cocrystal can be draw lessons from the other research fields, the cocrystal is formed by self-assembly of hydrogen bonding and π-π stacking interactions; (2) on the basis of the existing mature eutectic technology, a safe, efficient and practical synthesis method suitable for energetic cocrystal should be improved, promoting the theoretical research of cocrystal formation, and simulate and design the formation of energetic cocrystal from the theoretical level [41]. Institute for the thermal decomposition tests. We are especially thankful to Fengqi Zhao from Science and Technology on Combustion and Explosion Laboratory, Xi'an 710065, China and RuiQi Shen from Nanjing University of Science and Technology for providing many helpful suggestions and English statements.