Kinetics of Polymorphic Phase Transformations of <i>o</i>-Aminobenzoic Acid: Application of a Dispersive Kinetic Model Plus Molecular Dynamics Simulation of Prenucleation Aggregates

The specific rate at which one crystalline phase converts to another can vary as a function of time under isothermal conditions. This behavior gives rise to sigmoidal kinetic transients that are characteristic of nucleation and growth. Such curves are commonly fitted using the Johnson-Mehl-Avrami-Erofe'ev-Kolmogorov (JMAEK) equation. However, due to the ambiguity surrounding the time exponent in the JMAEK model, we present an alternative two-parameter dispersive kinetic model (DKM) and apply it to the study of solution-mediated polymorphic conversions of the prototypical molecule, <i>o</i>-aminobenzoic acid (<i>o</i>-ABA) [Jiang, S.; Jansens, P. J.; ter Horst, J. H. Control over polymorph formation of <i>o</i>-aminobenzoic acid. <i>Cryst. Growth Des.</i> <b>2010</b>, <i>10</i>, 2541-2547]. Using our DKM, we reconstructed a distribution of activation energies, <i>D</i>(<i>E</i>), from each experimental transient. Then, using <i>D</i>(<i>E</i>), a corresponding particle size distribution (PSD) of the critical nuclei formed during phase transformation is predicted. Lastly, molecular dynamics (MD) simulations are performed to study the prenucleation aggregation behavior of <i>o</i>-ABA in solution, under experimentally relevant conditions, to complement the kinetic modeling of the macroscopic phase conversion in the solid state. We observe that <i>o</i>-ABA molecules weakly associate with each other to form a variety of "loose" aggregates. These aggregates are mostly dimers and trimers exhibiting H-bonding and π-π interactions in various configurations that generally do not conform to any of the known crystal packing arrangements of the most common <i>o</i>-ABA polymorphs. Therefore, the observed molecular self-association is more consistent with a nonclassical nucleation pathway whereby monomer densification occurs ahead of cluster formation and, eventually, structural ordering. Our molecular-level simulations in solution complement the original study performed using experimental measurements on bulk crystals, with the DKM serving to bridge the scale gap between the two approaches by providing a window into the nanoscale species (nuclei), ultimately impacting the overall rate of conversion.