Prediction of Ethanol-Mediated Growth Morphology of Ammonium Dinitramide/Pyrazine-1,4-Dioxide Cocrystal at Different Temperatures

The crystal morphology of high energetic materials plays a crucial role in aspects of their safety performance such as impact sensitivity. In order to reveal the crystal morphology of ammonium dinitramide/pyrazine-1,4-dioxide (ADN/PDO) cocrystal at different temperatures, the modified attachment energy model (MAE) was used at 298, 303, 308, and 313 K to predict the morphology of the ADN/PDO cocrystal under vacuum and ethanol. The results showed that under vacuum conditions, five growth planes of the ADN/PDO cocrystal were given, which were (1 0 0), (0 1 1), (1 1 0), (1 1 −1), and (2 0 −2). Among them, the ratios of the (1 0 0) and (0 1 1) planes were 40.744% and 26.208%, respectively. In the (0 1 1) crystal plane, the value of S was 1.513. The (0 1 1) crystal plane was more conducive to the adsorption of ethanol molecules. The order of binding energy between the ADN/PDO cocrystal and ethanol solvent was (0 1 1) > (1 1 −1) > (2 0 −2) > (1 1 0) > (1 0 0). The radial distribution function analysis revealed that there were hydrogen bonds between the ethanol and the ADN cations, van der Waals interactions with the ADN anions. As the temperature increased, the aspect ratio of the ADN/PDO cocrystal was reduced, making the crystal more spherical, which helped to further reduce the sensitivity of this explosive.


Introduction
High-energy compounds are a class of compounds containing explosive groups that can independently undergo chemical reactions and output energy. They are widely used in exploration, aviation, and other fields. Ammonium dinitramide (ADN) [1] is an ionic energetic compound, which is in the α crystal form at normal temperature and pressure, with a crystal density of 1.812 g/cm 3 [2], an enthalpy of formation of −1.22 mJ/kg [3], and an oxygen balance of 25.8% [4]. Compared with ammonium perchlorate (AP), ADN has a higher enthalpy of formation [5], which can provide higher specific impulse when used in propellants. Since ADN is free of chlorine, its combustion products are clean, and its characteristic signal, such as primary smoke, etc., is low when it burns. However, the application of ADN is limited by the high hygroscopicity at room temperature [6,7].
Cocrystallization of energetic materials can effectively improve the physical and chemical properties [8,9]. Previous studies [10,11] have shown that cocrystallization is an effective way to improve the hygroscopicity of ADN. The moisture absorption rate experiment with an AD/18C6 (18-crown-6) cocrystal [12] showed that the hygroscopicity of the ADN/18C6 cocrystal was significantly lower than that of the ADN. The modified attachment energy was adopted by Xie [13] to explore the influence of ethanol on the morphology of an ADN/18C6 cocrystal at different temperatures. The results revealed morphology of an ADN/18C6 cocrystal at different temperatures. The results reve that the morphology of the ADN/18C6 cocrystal was close to spherical at 293 K, and cocrystal of the ADN/18C6 had lower hygroscopicity than the ADN at this temperat Pyrazine-1,4-dioxide (PDO) was selected by Michael K [14] to construct an ADN/P cocrystal with a molar ratio of 2:1. The hygroscopicity test showed that the hygroscop of the ADN decreased significantly after the cocrystallization of the ADN and PDO. thermal stability of the ADN/PDO cocrystal was higher than the ADN. Meanwhile energetic properties were improved over the ADN.
As an important indicator of crystals, their morphology has been widely consid [15][16][17]. The morphology of an energetic material can significantly affect its sensiti [18], which is affected by saturation, solvent, temperature, etc. [19,20]. In this work "explosive/solvent" double-layer model was constructed, and four temperatures w adopted to perform the molecular dynamics, which was carried out at 298, 303, 308, 313 K. The morphology of the ADN/PDO cocrystal was predicted with the help of a m ified attachment energy model. The binding energy, radial distribution function, and fusion coefficient were used to reveal the growth morphology of the ADN/PDO cocry in an ethanol environment.

Force Field Verification
As an ionic compound, ADN has strong electrostatic interactions between anions cations. In order to accurately calculate the interaction energy between the ADN/PDO ethanol, the DMol 3 program was selected to calculate the Mulliken charge of ADN/PDO unitcell and ethanol, which was adopted as the atomic charge for subsequ calculations. The Perdew-Burke-Ernzerhof (PBE) correlation was used to calculate exchange-correlation energy. The all-electron method was implemented to treat the electrons. The double numerical plus polarization numerical basis set (DNP) was adap The Mulliken charge is shown in Figure 1. In classical molecular dynamics, the force field parameters directly determine the curacy of calculation results. In order to accurately describe the interaction in ADN/PDO cocrystal, COMPASSIII, PCFF, CVFF, Universal, and Dreiding force fields the Mulliken charge were adopted to describe the properties of the atom in the w In classical molecular dynamics, the force field parameters directly determine the accuracy of calculation results. In order to accurately describe the interaction in the ADN/PDO cocrystal, COMPASSIII, PCFF, CVFF, Universal, and Dreiding force fields and the Mulliken charge were adopted to describe the properties of the atom in the whole simulation model. The optimized unit cell parameters of the ADN/PDO cocrystal are given in Table 1. The relative error (RE) of the unit cell parameters in Table 1  by CVFF was smaller than the results under other force fields, and it was closer to the experimental value. Meanwhile, the calculation of Yusop [21] revealed that the CVFF force field was suitable for the molecular dynamics simulation of ethanol. Therefore, the CVFF force field was adopted to perform the molecular dynamics simulation.

Morphology of the ADN/PDO in Vacuum
The morphology of the ADN/PDO in vacuum was calculated based on the AE model, as shown in Figure 2. Under vacuum conditions, the important growth planes and area ratios of the ADN/PDO cocrystal are shown in Table 2. It can be seen from Figure 2 that the vacuum morphology of the ADN/PDO cocrystal was approximately hexagonal prism, and the aspect ratio of the crystal was 2.569. Its crystal morphology was mainly composed of five growth crystal planes (1 0 0), (0 1 1), (1 1 0), (1 1 −1), and (2 0 −2). Among them, the (1 0 0) surface had the largest exposed area, accounting for 40.744%, and had the greatest morphological importance. According to Formula (3), it can be seen that the absolute value of the attachment energy was proportional to the relative growth rate of the corresponding crystal plane. The order of the relative growth rate of each crystal plane of the ADN/PDO was (2 0 −2) > (1 1 −1) > (1 1 0) > (0 1 1) > (1 0 0).
Molecules 2023, 28, x FOR PEER REVIEW 3 of 14 simulation model. The optimized unit cell parameters of the ADN/PDO cocrystal are given in Table 1. The relative error (RE) of the unit cell parameters in Table 1 revealed that the RE optimized by CVFF was smaller than the results under other force fields, and it was closer to the experimental value. Meanwhile, the calculation of Yusop [21] revealed that the CVFF force field was suitable for the molecular dynamics simulation of ethanol. Therefore, the CVFF force field was adopted to perform the molecular dynamics simulation.

Morphology of the ADN/PDO in Vacuum
The morphology of the ADN/PDO in vacuum was calculated based on the AE model, as shown in Figure 2. Under vacuum conditions, the important growth planes and area ratios of the ADN/PDO cocrystal are shown in Table 2. It can be seen from Figure 2 that the vacuum morphology of the ADN/PDO cocrystal was approximately hexagonal prism, and the aspect ratio of the crystal was 2.569. Its crystal morphology was mainly composed of five growth crystal planes (1 0 0), (0 1 1), (1 1 0), (1 1 −1), and (2 0 −2). Among them, the (1 0 0) surface had the largest exposed area, accounting for 40.744%, and had the greatest morphological importance. According to Formula (3), it can be seen that the absolute value of the attachment energy was proportional to the relative growth rate of the corresponding crystal plane. The order of the relative growth rate of each crystal plane of the ADN/PDO was (2 0 −2) > (1 1 −1) > (1 1 0) > (0 1 1) > (1 0 0).   In the modified attachment energy (MAE) model, S was used to characterize the surface roughness of the (h k l) crystal plane [22]. The larger the value of S, the rougher the surface of the corresponding crystal plane, and the more conducive to the adsorption of solvents [18,23]. It can be seen that the values of S corresponding to (1 0 0), (0 1 1), (1 1 0), (1 1 −1), and (2 0 −2) crystal planes were 1.223, 1.513, 1.359, 1.389, and 1.836, respectively. The packing patterns of five crystal planes of the ADN/PDO and the corresponding Connolly surfaces are shown in Figure 3. Based on Figure 3, more chemical groups were exposed in the (0 1 1) and (2 0 −2) crystal planes compared with the other three crystal planes. It can be seen from Figure 3b,e that the (0 1 1) crystal plane was mainly exposed to the PDO and the cations of the ADN, and the (2 0 −2) crystal plane was exposed to both the cations and anions of the ADN and PDO. This results in the (2 0 −2) crystal plane were more bumpy and rougher than the (0 1 1) crystal plane. The (2 0 −2) plane was the roughest and had more adsorption sites. It was conducive to absorbing the molecule of solvent. This had a certain hindering effect on the growth of the crystal plane.  the roughest and had more adsorption sites. It was conducive to absorbing the molecule of solvent. This had a certain hindering effect on the growth of the crystal plane.   When a crystal grows in a solution, the solute molecules first diffuse to the crystal surface through diffusion, then overcome the desorption energy barrier of the solvent, and finally adsorb to the crystal surface to complete the crystal growth [23,24]. Therefore, the crystal morphology can be significantly affected by the adsorption of solvents on the crystal surface. The greater the binding energy between solvent and crystal plane, the stronger the binding effect, the greater the desorption energy barrier that the solute needs to overcome, and the stronger the inhibitory effect of the solvent on the growth of crystal plane. The formula for calculating the binding energy is

Morphology of the AND/PDO in Solvent
where E bind , E int are the binding energy and interaction energy between the crystal plane and the solvent, respectively. E tot is the total energy of the mixed system of crystal plane and solvent, E cry is the energy of the crystal plane in the model, E sol is the energy of the solvent in the model. The calculation results are shown in Table 3. Table 3. Energy details of the layered model for each crystal face of ADN/PDO-EtOH at different temperatures. According to Table 3, it can be seen that the binding energy between the ADN/PDO cocrystal and the ethanol solvent were positive on all important growth crystal planes, indicating that the interaction between the solute and the solvent was dominated by attraction. The variations of the binding energy between the ADN/PDO and the ethanol on each important crystal face of the ADN/PDO cocrystal are shown in Figure 4, which were obtained under different temperatures. The order of binding energy between five crystal planes of the ADN/PDO cocrystal and ethanol was (0 1 1) > (1 1 −1) > (2 0 −2) > (1 1 0) > (1 0 0). It can be revealed that the binding effect between the ethanol and the ADN/PDO was the strongest on the (0 1 1) crystal plane. Meanwhile, ethanol had the strongest inhibitory effect on the growth of this crystal plane. It can be seen from Figure 4 that the binding energy between the ADN/PDO crystal plane and ethanol first increased and then decreased as the temperature increased. Its value was always positive,

Radial Distribution Function
The radial distribution function (RDF) is often adopted to analyze non-bonding interactions between atomic pairs. It represents the probability of finding atom B at a distance r from atom A. In general, non-bonding interactions between atomic pairs include hydrogen bonding (less than 3.1 Å), van der Waals forces (3.1-5.0 Å), and electrostatic forces (greater than 5.0 Å).
Taking the (0 1 1) crystal plane as an example, the RDF results are shown in Figure 5. The black, red, and blue curves corresponded to the RDF between the oxygen atom in ethanol (EtOH-O) and the nitrogen atom connected to two nitro groups in the ADN anion (ADN-N), the RDF between the EtOH-O and the oxygen atom in the PDO (PDO-O), and the RDF between the EtOH-O and the nitrogen atom in the ADN cations (ADN-N1), respectively. The coordinates of the local maximum value of each curve are given in the Figure 5. Based on Figure 5a, the RDF of the EtOH-O and ADN-N1 had the first peak at r = 2.85 Å at 298 K, indicating that there was a hydrogen bond between the ethanol and ADN cations. The peak of the black curve appeared at the position of r = 4.91 Å, indicating that there was a van der Waals interaction between the ethanol and the ADN anions. There were two peaks in the radial distribution function between the EtOH-O and PDO-O, which appeared at r = 2.67and 4.67 Å respectively. This revealed the presence of hydrogen bonds and van der Waals interactions between the ADN and the PDO.

Radial Distribution Function
The radial distribution function (RDF) is often adopted to analyze non-bonding interactions between atomic pairs. It represents the probability of finding atom B at a distance r from atom A. In general, non-bonding interactions between atomic pairs include hydrogen bonding (less than 3.1 Å), van der Waals forces (3.1-5.0 Å), and electrostatic forces (greater than 5.0 Å).
Taking the (0 1 1) crystal plane as an example, the RDF results are shown in Figure 5. Similarly, when the temperatures were 303, 308, and 313 K, there were hydrogen bonds and van der Waals interactions between the ethanol and the PDO. There were hydrogen bonds and van der Waals interactions between the ethanol and the cation and anion of the ADN, respectively. As the temperature increased, the number of local maximum values of the three RDFs remained unchanged, while the abscissa of the local maximum point moved to the right as a whole. This revealed that the increase in temperature led to the change of the hydrogen bond and van der Waals interaction between the ethanol and the ADN, PDO, which first increased and then decreased. It was consistent with the binding energy results in Figure 4.

Diffusion Coefficients
Generally, the diffusion of a solvent is greatly affected by temperature. The diffusion of a solvent on the crystal surface may cause it to adsorb on the crystal surface. Absorption of solute molecules to the crystal surface is blocked by this process. The growth of the corresponding crystal planes is also inhibited. To reveal the effect of temperature on the diffusion of ethanol molecules on the AND/PDO crystal plane, the calculated trajectories were processed, and the corresponding mean square displacement curves were calculated, which is given in Figure 6. Similarly, when the temperatures were 303, 308, and 313 K, there were hydrogen bonds and van der Waals interactions between the ethanol and the PDO. There were hydrogen bonds and van der Waals interactions between the ethanol and the cation and anion of the ADN, respectively. As the temperature increased, the number of local maximum values of the three RDFs remained unchanged, while the abscissa of the local maximum point moved to the right as a whole. This revealed that the increase in temperature led to the change of the hydrogen bond and van der Waals interaction between the ethanol and the ADN, PDO, which first increased and then decreased. It was consistent with the binding energy results in Figure 4.

Diffusion Coefficients
Generally, the diffusion of a solvent is greatly affected by temperature. The diffusion of a solvent on the crystal surface may cause it to adsorb on the crystal surface. Absorption of solute molecules to the crystal surface is blocked by this process. The growth of the corresponding crystal planes is also inhibited. To reveal the effect of temperature on the diffusion of ethanol molecules on the AND/PDO crystal plane, the calculated trajectories were processed, and the corresponding mean square displacement curves were calculated, which is given in Figure 6. According to Einstein's law of diffusion, the formula for calculating the diffusion coefficient is where D is the diffusion coefficient. Therefore, the diffusion coefficient is one-sixth of the slope of the mean square displacement curve. The diffusion coefficients of the ethanol on each crystal surface at different temperatures are given in Table 4. Taking the (0 1 1) surface as an example, when the temperature was 298 K, the diffusion coefficient of the ethanol in the system was 0.52 × 10 −8 m 2 ·s −1 . When the temperature rose to 313 K, the diffusion coefficient of the system was 0.64 × 10 −8 m 2 ·s −1 . As the According to Einstein's law of diffusion, the formula for calculating the diffusion coefficient is where D is the diffusion coefficient. Therefore, the diffusion coefficient is one-sixth of the slope of the mean square displacement curve. The diffusion coefficients of the ethanol on each crystal surface at different temperatures are given in Table 4. Taking the (0 1 1) surface as an example, when the temperature was 298 K, the diffusion coefficient of the ethanol in the system was 0.52 × 10 −8 m 2 ·s −1 . When the temperature rose to 313 K, the diffusion coefficient of the system was 0.64 × 10 −8 m 2 ·s −1 . As the temperature increased, the diffusion coefficient of the ethanol under the same crystal plane increased gradually. At the same temperature, the largest diffusion coefficient of the ethanol in the layered model was obtained under the (1 1 0) crystal plane. However, the final morphology of the crystals in the solvent was determined by the diffusion of the solute and solvent and the competitive adsorption of the solute and solvent on the crystal surface.

Morphology Analysis
According to the modified AE model, the corrected attachment energy and aspect ratio between the crystal plane of the ADN/PDO cocrystal and ethanol are shown in Table 5, which was calculated under different temperatures. The calculated morphology of the ADN/PDO cocrystal in ethanol at different temperatures and the corresponding experimental crystal morphology of the ADN/PDO cocrystal [25] are shown in Figure 7. Based on Table 5, it can be seen that the (0 1 1) crystal surface of the ADN/PDO cocrystal had a greater interaction with the ethanol than the (1 0 0) plane at the four temperatures (298, 303, 308, and 313 K). This indicated that the ethanol had a strong inhibitory effect on the growth of the ADN/PDO (0 1 1) crystal surface.  According to the results of S on each crystal plane of the ADN/PDO cocrystal, the (2 0 −2) crystal plane was the roughest, which was the most favorable for the adsorption of solvents. However, the binding energy between the solvent and the (2 0 −2) crystal plane was weak, and the diffusion coefficient of ethanol on the (2 0 −2) plane was small. Under the comprehensive influence of various factors, the relative growth rate of the (2 0 −2) plane was large, which made the (2 0 −2) crystal surface disappear eventually during the growth process. Similarly, the roughness of the (0 1 1) surface and the binding energy between explosives and ethanol were large. This gave the (0 1 1) surface the largest area ratio in the ethanol solvent.
When the temperature was lower than 308 K, the morphology of the ADN/PDO cocrystal in the ethanol solvent was a quadrangular prism. The three crystal planes (1 1 0), (1 1 −1), and (2 0 −2) eventually disappeared during the growth process owing to the high relative growth rate. When the temperature increased from 298 to 308K, the (0 1 1) surface area ratio of the ADN/PDO cocrystal was greater than 80%, and the aspect ratios of the ADN/PDO cocrystal were 3.324, 2.653, and 4.090, respectively. When the temperature was 313 K, the morphology of the ADN/PDO cocrystal was irregular prism, and the aspect ratio was 2.353. The ADN/PDO cocrystal morphology calculated by the modified AE model was the closest to the experimental one at 313 K. The presented modeling and simulations can also be applied for the prediction of the drug crystal topologies in a pharmaceutical application employing drug delivery carriers [26][27][28].

Modified Attachment Energy Model
The Attachment Energy Model (AE model) was developed by Hartman and Bennema [25,29] on the basis of the Periodic Bond Chain (PBC) theory. In the AE model, the relative growth rate of each crystal plane of the crystal is proportional to the absolute value of the crystal plane attachment energy.
where Rhkl is the relative growth rate, and Eatt is the attachment energy of the crystal plane, kcal·mol −1 . The attachment energy (Eatt) is defined as the energy released by a wafer with a thickness of dhkl attached to the (h k l) crystal plane. The formula is where Eatt, Elatt, Eslice are the attachment energy, lattice energy of crystal, and the energy of growth slice, respectively. According to the results of S on each crystal plane of the ADN/PDO cocrystal, the (2 0 −2) crystal plane was the roughest, which was the most favorable for the adsorption of solvents. However, the binding energy between the solvent and the (2 0 −2) crystal plane was weak, and the diffusion coefficient of ethanol on the (2 0 −2) plane was small. Under the comprehensive influence of various factors, the relative growth rate of the (2 0 −2) plane was large, which made the (2 0 −2) crystal surface disappear eventually during the growth process. Similarly, the roughness of the (0 1 1) surface and the binding energy between explosives and ethanol were large. This gave the (0 1 1) surface the largest area ratio in the ethanol solvent.
When the temperature was lower than 308 K, the morphology of the ADN/PDO cocrystal in the ethanol solvent was a quadrangular prism. The three crystal planes (1 1 0), (1 1 −1), and (2 0 −2) eventually disappeared during the growth process owing to the high relative growth rate. When the temperature increased from 298 to 308K, the (0 1 1) surface area ratio of the ADN/PDO cocrystal was greater than 80%, and the aspect ratios of the ADN/PDO cocrystal were 3.324, 2.653, and 4.090, respectively. When the temperature was 313 K, the morphology of the ADN/PDO cocrystal was irregular prism, and the aspect ratio was 2.353. The ADN/PDO cocrystal morphology calculated by the modified AE model was the closest to the experimental one at 313 K. The presented modeling and simulations can also be applied for the prediction of the drug crystal topologies in a pharmaceutical application employing drug delivery carriers [26][27][28].

Modified Attachment Energy Model
The Attachment Energy Model (AE model) was developed by Hartman and Bennema [25,29] on the basis of the Periodic Bond Chain (PBC) theory. In the AE model, the relative growth rate of each crystal plane of the crystal is proportional to the absolute value of the crystal plane attachment energy.
where R hkl is the relative growth rate, and E att is the attachment energy of the crystal plane, kcal·mol −1 . The attachment energy (E att ) is defined as the energy released by a wafer with a thickness of d hkl attached to the (h k l) crystal plane. The formula is where E att , E latt , E slice are the attachment energy, lattice energy of crystal, and the energy of growth slice, respectively.
By calculating the attachment energy of the crystal face, the habit of the crystal can be predicted with the help of the AE model. However, the external environment of crystal growth is not considered by the AE model. It is difficult to accurately reveal the actual growth process of crystals in solution [30,31]. Therefore, the AE model should be corrected, and the modified formula for calculating the attachment energy is where S is used to describe the surface characteristics and is defined as S = A acc A hkl . Among them, A acc is the accessible area of solvent on the crystal plane unit (h k l), A hkl is the cross-sectional area of the crystal plane unit (h k l). E s represents the effect of the solvent on the crystal plane growth, which is defined as where A box and E int are the cross-sectional area of the simulation box and the interaction energy between the crystal layer and the solvent, respectively. E int is defined as where E tot , E cry , E sol are the total energy of the mixed system of crystal plane and solvent, the energy of crystal plane in the calculation model, and the energy of the solvent in this model, respectively. Meanwhile, under solvent conditions, the relative growth rate of the crystal is proportional to the absolute value of the corrected attachment energy.

Computational Methods
The initial unit cell for the ADN/PDO cocrystal was obtained from the Cambridge Crystallographic Data Centre (CCDC) [14]. The initial unit cell parameters of the ADN/PDO cocrystal were a = 11.592 Å, b = 8.188 Å, c = 7.227 Å, α = γ = 90 • , β = 101.236 • , and the space group was P21/c, belonging to the monoclinic crystal system. The molecular structure of the ADN and PDO and the unit cell structure of the ADN/PDO cocrystal are shown in Figure 8. In Figure 8, the gray, white, red, and blue balls correspond to carbon, hydrogen, oxygen, and nitrogen atoms, respectively.  The flowchart of the calculation model is given in Figure 9. At first, the stable structure of the ADN/PDO unit cell was obtained after geometry optimization. The Smart optimization algorithm was used, which is the built-in algorithm in Materials Studio 2020.
The Ewald summation method was selected as the electrostatic force summation method, and the precision was set to 0.001 kcal·mol −1 . The atom-based method was selected as the van der Waals force summation method, and the cut-off radius was set to 12.5 Å. The morphology of the AND/PDO in vacuum was predicted with the help of the AE model. The important growth planes under vacuum conditions were determined. Then, the important crystal planes of the ADN/PDO were cut to expand it into a 4 × 4 × 4 supercell structure. A solvent layer containing 400 ethanol molecules was constructed using the Amorphous Cell module. The Build Layers function was used to construct an "explosive/solvent" double-layer model with a 50 Å vacuum layer above the solvent layer.  The model was optimized using the Forcite module of Materials Studio 2020 lowed by molecular dynamics calculations at specified temperatures (298, 303, 308 313 K). When molecular dynamics calculations were performed, the ensemble an simulation time were set as NVT and 500 ps, respectively. A total of 500,000 calcul steps were performed, and the timestep was set as 1 fs. The Andersen method adopted to control the temperature through the whole simulation process. The initia locity of the particle was randomly assigned according to a Gaussian distribution. Fin the homemade script was used to calculate the interaction energy between the ADN/ The model was optimized using the Forcite module of Materials Studio 2020, followed by molecular dynamics calculations at specified temperatures (298, 303, 308, and 313 K). When molecular dynamics calculations were performed, the ensemble and the simulation time were set as NVT and 500 ps, respectively. A total of 500,000 calculation steps were performed, and the timestep was set as 1 fs. The Andersen method was adopted to control the temperature through the whole simulation process. The initial velocity of the particle was randomly assigned according to a Gaussian distribution. Finally, the homemade script was used to calculate the interaction energy between the ADN/PDO explosive and ethanol, the radial distribution function, the mean square displacement, and the mass density distribution based on the stable part of the calculated trajectory.

Conclusions
In summary, the modified attachment energy (MAE) model was adopted to investigate the morphology of the ADN/PDO cocrystal in vacuum and ethanol at different temperatures (298, 303, 308, and 313 K). The CVFF forcefield and Mulliken charges were used to perform the molecular dynamics simulation. The main conclusions of this work are as follows: (1) The growth morphology of the ADN/PDO is hexagonal prism in vacuum, and the five main crystal surfaces are (1 0 0), (0 1 1), (1 1 0), (1 1 −1), and (2 0 −2). Among them, the (1 0 0) surface has the largest exposed area, accounting for 40.744%. (2) The binding energy between the ADN/PDO cocrystal and ethanol solvent is positive on all important growth planes. The order of binding energy is (0 1 1) > (1 1 −1) > (2 0 −2) > (1 1 0) > (1 0 0). The binding effect between the ethanol and ADN/PDO is strongest on the (0 1 1) crystal plane. (3) The radial distribution function analysis of the (0 1 1) crystal plane showed that there are hydrogen bonds between the ethanol and ADN cations, van der Waals interactions with the ADN anions, and hydrogen bonds and van der Waals interactions with the PDO at the same time. In the (0 1 1) crystal plane, the value of S is 1.513, which indicates that this surface has a large roughness. This is more conducive to the adsorption of ethanol molecules. (4) As the temperature increases, the diffusion coefficient of the ethanol under the same crystal plane increases gradually. Meanwhile, the morphology analysis indicated that increasing the temperature is beneficial to reducing the aspect ratio of the crystal. This is conducive to the reduction of explosive sensitivity.