The crystal growth of urea was analyzed
with all-atom molecular
dynamics (MD) simulation for the (001) and (110) faces in contact
with aqueous solutions. The local environment of a crystallizing molecule
was treated in terms of the numbers of crystalline neighbors and the
orientation relative to the crystal surface, and the molecular-level
inhomogeneity of a growing surface was addressed by decomposing the
overall rate of growth into a sum of the contributions conditioned
by the local structure and orientation mode. The contrast of the growth
mechanism between the (001) and (110) faces was then evidenced by
the local contributions, and the roles of the outer layers of the
crystal toward the liquid region were pointed out for (001). The effect
of the additive species in the liquid on the crystal growth of urea
was investigated with biuret, N,N-dimethylformamide (DMF), and acetone. The growth was observed to
be suppressed more strongly in the order of biuret > DMF > acetone,
and it was found that the ordering of suppression by the additive
is common irrespective of the local environment of a crystallizing
urea. This finding implies that the additive’s effect on the
crystal growth can be predicted by treating the flat surface, which
is a convenient system for detailed analyses at atomic resolution.
The correspondence to the free energy of adsorption of the additive
was then examined for the additive-induced modulation of the growth
rate. It was seen that the adsorption free energy correlates to the
extent of modulation of the growth rate, and the interaction components
that govern the adsorption propensity were identified.