This work could help to better understand the solvent effects on crystal habits and aspect ratio changes at the molecular level, which provide some guidance for solvent selection in industrial crystallization processes. With the catechol crystal habits acquired using both experimental and simulation methods in isopropanol, methyl acetate and ethyl acetate, solvent effects on crystal morphology were explored based on the modified attachment energy model. Firstly, morphologically dominant crystal faces were obtained with the predicted crystal habit in vacuum. Then, modified attachment energies were calculated by the molecular dynamics simulation to modify the crystal shapes in a real solvent environment, and the simulation results were in agreement with the experimental ones. Meanwhile, the surface properties such as roughness and the diffusion coefficient were introduced to analyze the solvent adsorption behaviors and the radial distribution function curves were generated to distinguish diverse types of interactions like hydrogen bonds and van der Waals forces. Results show that the catechol crystal habits were affected by the combination of the attachment energy, surface structures and molecular interaction types. Moreover, the changing aspect ratios of catechol crystals are closely related to the existence of hydrogen bonds which contribute to growth inhibition on specific faces.
Molecular dynamics simulations of solvent effects on the crystal morphology of lithium carbonate
