Molecular dynamics investigation on the morphology of HNIW affected by the growth condition

The crystal morphology of Hexanitrohexaazaisowurtzitane (HNIW) in vacuum and ethyl acetate was investigated by attachment energy (AE) model through molecular dynamics (MD). The simulated solvent-affected outcome agrees well with the experimental crystal habit cultivated from ethyl acetate. The influence of model size on the modified AE was discussed, and a reasonable model size was obtained. The model size can be broadened to calculate the morphology of other compounds. The radial distribution function (RDF) analyses and diffusion coefficient analyses were performed to explore the adsorption and diffusion behaviors of ethyl acetate molecules on HNIW surfaces. The simulation results can provide some guidance for the further crystallization process of HNIW.