Solvent selection for explaining the morphology of nitroguanidine crystal by molecular dynamics simulation

In this article, a method was performed to predict the morphology of needle-shaped crystals by analyzing the growth mechanisms for the various crystal faces. As an example, the crystal morphology of a nitroguanidine (NQ) was investigated via molecular dynamics simulations. The modified attachment energy (MEA) model was constructed by introducing surface chemistry terms and the relevant morphology of the habit crystal faces. The results indicate that the growth morphology of NQ in vacuum is dominated by {2 2 0}, {0 4 0}, {1 1 1}, {1 3 1} and {3 1 1} faces. The {2 2 0} and {0 4 0} faces are parallel to the elongation direction of the crystal, while the other faces are at the needle tips direction. The atoms or atomic groups exposed in crystal surface were used to analyze the relationship between structure and morphology. Compared to the surrounding faces, the needle tip faces have a large number of polar atoms or atomic groups. The needle tip faces have a high electronegativity on N, O atoms via molecular electrostatic potential (ESP) analysis. Furthermore, the protic solvent was used to reduce the attachment energy of the tip surfaces for achieving the purpose of inhibiting the growth of needle tips. Gamma-butyrolactone as the selected solvent inhibited effectively the growth of the needle tip faces. The predicted result is serviceable for the formation design.