Aromatic stacking - a key step in nucleation

Influence of Solvent Selection on the Crystallizability and Polymorphic Selectivity Associated with the Formation of the "Disappeared" Form I Polymorph of Ritonavir

The comparative crystallizability and polymorphic selectivity of ritonavir, a novel protease inhibitor for the treatment of acquired immune-deficiency syndrome, as a function of solvent selection are examined through an integrated and self-consistent experimental and computational molecular modeling study. Recrystallization at high supersaturation by rapid cooling at 283.15 K is found to produce the metastable "disappeared" polymorphic form I from acetone, ethyl acetate, acetonitrile, and toluene solutions in contrast to ethanol which produces the stable form II. Concomitant crystallization of the other known solid forms is not found under these conditions. Isothermal crystallization studies using turbidometric detection based upon classical nucleation theory reveal that, for an equal induction time, the required driving force needed to initiate solution nucleation decreases with solubility in the order of ethanol, acetone, acetonitrile, ethyl acetate, and toluene consistent with the expected desolvation behavior predicted from the calculated solute solvation free energies. Molecular dynamics simulations of the molecular and intermolecular chemistry reveal the presence of conformational interplay between intramolecular and intermolecular interactions within the solution phase. These encompass the solvent-dependent formation of intramolecular O-H...O hydrogen bonding between the hydroxyl and carbamate groups coupled with differing conformations of the hydroxyl's shielding phenyl groups. These conformational preferences and their relative interaction propensities, as a function of solvent selection, may play a rate-limiting role in the crystallization behavior by not only inhibiting to different degrees the nucleation process but also restricting the assembly of the optimal intermolecular hydrogen bonding network needed for the formation of the stable form II polymorph.

The role of alkyl chain length in the melt and solution crystallization of paliperidone aliphatic prodrugs

Fatty acid-derivative prodrugs have been utilized extensively to improve the physicochemical, biopharmaceutical and pharmacokinetic properties of active pharmaceutical ingredients. However, to our knowledge, the crystallization behavior of prodrugs modified with different fatty acids has not been explored. In the present work, a series of paliperidone aliphatic prodrugs with alkyl chain lengths ranging from C4 to C16 was investigated with respect to crystal structure, crystal morphology and crystallization kinetics. The paliperidone derivatives exhibited isostructural crystal packing, despite the different alkyl chain lengths, and crystallized with the dominant (100) face in both melt and solution. The rate of crystallization for paliperidone derivatives in the melt increases with alkyl chain length owing to greater molecular mobility. In contrast, the longer chains prolong the nucleation induction time and reduce the crystal growth kinetics in solution. The results show a correlation between difficulty of nucleation in solution and the interfacial energy. This work provides insight into the crystallization behavior of paliperidone aliphatic prodrugs and reveals that the role of alkyl chain length in the crystallization behavior has a strong dependence on the crystallization method.

Novel Para-Aminobenzoic Acid Analogs and Their Potential Therapeutic Applications

A "building block" is a key component that plays a substantial and critical function in the pharmaceutical research and development industry. Given its structural versatility and ability to undergo substitutions at both the amino and carboxyl groups, para-aminobenzoic acid (PABA) is a commonly used building block in pharmaceuticals. Therefore, it is great for the development of a wide range of novel molecules with potential medical applications. Anticancer, anti-Alzheimer's, antibacterial, antiviral, antioxidant, and anti-inflammatory properties have been observed in PABA compounds, suggesting their potential as therapeutic agents in future clinical trials. PABA-based therapeutic chemicals as molecular targets and their usage in biological processes are the primary focus of this review study. PABA's unique features make it a strong candidate for inclusion in a massive chemical database of molecules having drug-like effects. Based on the current literature, further investigation is needed to evaluate the safety and efficacy of PABA derivatives in clinical investigations and better understand the specific mechanism of action revealed by these compounds.

The Influence of Solvent Selection upon the Crystallizability and Nucleation Kinetics of Tolfenamic Acid Form II

The influence of the solution environment on the solution thermodynamics, crystallizability, and nucleation of tolfenamic acid (TFA) in five different solvents (isopropanol, ethanol, methanol, toluene, and acetonitrile) is examined using an integrated workflow encompassing both experimental studies and intermolecular modeling. The solubility of TFA in isopropanol is found to be the highest, consistent with the strongest solvent-solute interactions, and a concomitantly higher than ideal solubility. The crystallizability is found to be highly dependent on the solvent type with the overall order being isopropanol < ethanol < methanol < toluene < acetonitrile with the widest solution metastable zone width in isopropanol (24.49 to 47.41 °C) and the narrowest in acetonitrile (8.23 to 16.17 °C). Nucleation is found to occur via progressive mechanism in all the solvents studied. The calculated nucleation parameters reveal a considerably higher interfacial tension and larger critical nucleus radius in the isopropanol solutions, indicating the higher energy barrier hindering nucleation and hence lowering the nucleation rate. This is supported by diffusion coefficient measurements which are lowest in isopropanol, highlighting the lower molecular diffusion in the bulk of solution compared to the other solutions. The TFA concentration and critical supersaturation at the crystallization onset is found to be directly correlated with TFA/isopropanol solutions having the highest values of solubility and critical supersaturation. Intermolecular modeling of solute-solvent interactions supports the experimental observations of the solubility and crystallizability, highlighting the importance of understanding solvent selection and solution state structure at the molecular level in directing the solubility, solute mass transfer, crystallizability, and nucleation kinetics.

Molecular Mechanism of Organic Crystal Nucleation: A Perspective of Solution Chemistry and Polymorphism

Crystal nucleation determining the formation and assembly pathway of first organic materials is the central science of various scientific disciplines such as chemical, geochemical, biological, and synthetic materials. However, our current understanding of the molecular mechanisms of nucleation remains limited. Over the past decades, the advancements of new experimental and computational techniques have renewed numerous interests in detailed molecular mechanisms of crystal nucleation, especially structure evolution and solution chemistry. These efforts bifurcate into two categories: (modified) classical nucleation theory (CNT) and non-classical nucleation mechanisms. In this review, we briefly introduce the two nucleation mechanisms and summarize current molecular understandings of crystal nucleation that are specifically applied in polymorphic crystallization systems of small organic molecules. Many important aspects of crystal nucleation including molecular association, solvation, aromatic interactions, and hierarchy in intermolecular interactions were examined and discussed for a series of organic molecular systems. The new understandings relating to molecular self-assembly in nucleating systems have suggested more complex multiple nucleation pathways that are associated with the formation and evolution of molecular aggregates in solution.

Formation and stabilization mechanism of mesoscale clusters in solution

To understand the existence of complex meso-sized solute-rich clusters, which challenge the understanding of phases and phase equilibria, the formation and stabilization mechanisms of clusters in solution during nucleation of crystals and the associated physico-chemical rules are studied in detail. An essential part of the mechanism is the formation of long-lived oligomers between solute molecules. By means of density functional theory simulation and nuclear magnetic resonance experiments, this work showed that the oligomers in solution tend to be π-π stacking dimers. Clusters are formed under the combined effect of diffusion and monomer-dimer reaction. The physically meaningful quantities such as the monomer-dimer reaction rate constants and the diffusion coefficients of both species were obtained by reaction-diffusion kinetics and diffusion-ordered spectroscopy results. The evolution of cluster radius as a function of time, and the qualitative spatial distributions of monomer and dimer densities under steady-state were plotted to better understand the formation process and the nature of the clusters.

Fourier Transform Infrared Spectroscopy and Vibrational Circular Dichroism Assisted Elucidation of the Solution-State Supramolecular Speciation in Racemic and Enantiopure Ketoprofen

The molecular structure and solution-state molecular interactions in the popular non-steroidal anti-inflammatory drug, ketoprofen, are extensively studied with the aim of gaining a better understanding of the chemical behavior of its solution state and its connection to its nucleation pathway and crystallization outcome. Using as reference solid-state X-ray structures of enantiomeric and racemic forms of ketoprofen, a set of self-assembly models underpinned by density functional theory calculations has been considered for the analysis of spectroscopic data, infrared (IR) and vibrational circular dichroism (VCD), obtained for solutions of the samples as a function of composition and solvent. From our results it can be concluded that, contrary to the general belief for generic carboxylic acids, there are no cyclic dimeric structures of ketoprofen present in solution, but rather linear arrays made up of two (in high polar or diluted media) or more units (in low polar or low dilution media). This observation is in line with the idea that the weak contacts (other than H-bonding) would hold the key to molecular self-assembly, in agreement with recent studies on other aromatic carboxylic acids.

Functional pH-responsive polymers containing dynamic enaminone linkages for the release of active organic amines

Dynamic covalent bonds have attracted attention for the development of pH-responsive polymers, however, studies using acid-cleavable enaminone linkages as a means of controlled drug release have been limited. Herein, we...

Pre-Nucleation Clusters Predict Crystal Structures in Models of Chiral Molecules

Kinetics can play an important role in the crystallization of molecules and can give rise to polymorphism, the tendency of molecules to form more than one crystal structure. Current computational methods of crystal structure prediction, however, focus almost exclusively on identifying the thermodynamically stable polymorph. Kinetic factors of nucleation and growth are often neglected because the underlying microscopic processes can be complex and accurate rate calculations are numerically cumbersome. In this work, we use molecular dynamics computer simulations to study simple molecular models that reproduce the crystallization behavior of real chiral molecules, including the formation of enantiopure and racemic crystals, as well as polymorphism. A significant fraction of these molecules forms crystals that do not have the lowest free energy. We demonstrate that at high supersaturation crystal formation can be accurately predicted by considering the similarities between oligomeric species in solution and molecular motifs in the crystal structure. For the case of racemic mixtures, we even find that knowledge of crystal free energies is not necessary and kinetic considerations are sufficient to determine if the system will undergo spontaneous chiral separation. Our results suggest conceptually simple ways of improving current crystal structure prediction methods.

Acid-Catalyzed RNA-Oligomerization from 3',5'-cGMP

The assembly of ancient informational polymers from nucleotide precursors is the central challenge of life's origin on our planet. Among the possible solutions, dry polymerization of 3’,5’‐cyclic guanosine monophosphate (3’,5’‐cGMP) has been proposed as a candidate to create oligonucleotides of 15–20 units in length. However, the reported sensitivity of the reaction to the presence of cations raised questions of whether this chemistry could be relevant in a geological context. The experiments in this study show that the presence of cations is not restrictive as long as the reaction is conducted in an acidic environment, in contrast to previous reports that suggested optimal conditions at pH 9.

Unveiling the self-association and desolvation in crystal nucleation

As the first step in the crystallization process, nucleation has been studied by many researchers. In this work, phenacetin (PHEN) was selected as a model compound to investigate the relationship between the solvent and nucleation kinetics. Induction times at different supersaturation in six solvents were measured. FTIR and NMR spectroscopy were employed to explore the solvent-solute interactions and the self-association properties in solution. Density functional theory (DFT) was adopted to evaluate the strength of solute-solvent interactions and the molecular conformations in different solvents. Based on these spectroscopy data, molecular simulation and nucleation kinetic results, a comprehensive understanding of the relationship between molecular structure, crystal structure, solution chemistry and nucleation dynamics is discussed. Both the solute-solvent interaction strength and the supramolecular structure formed by the self-association of solute molecules affect the nucleation rate. The findings reported here shed new light on the molecular mechanism of nucleation in solution.

Solution and Solid-State Photophysical Properties of Positional Isomeric Acrylonitrile Derivatives with Core Pyridine and Phenyl Moieties: Experimental and DFT Studies

The compounds I (Z)-2-(phenyl)-3-(2,4,5-trimethoxyphenyl)acrylonitrile with one side (2,4,5-MeO-), one symmetrical (2Z,2'Z)-2,2'-(1,4-phenylene)bis(3-(2,4,5-trimethoxyphenyl)acrylonitrile), II (both sides with (2,4,5-MeO-), and three positional isomers with pyridine (Z)-2-(pyridin-2- 3, or 4-yl)-3-(2,4,5-trimethoxyphenyl)acrylonitrile, III-V were synthetized and characterized by UV-Vis, fluorescence, IR, H1-NMR, and EI mass spectrometry as well as single crystal X-ray diffraction (SCXRD). The optical properties were strongly influenced by the solvent (hyperchromic and hypochromic shift), which were compared with the solid state. According to the solvatochromism theory, the excited-state (μe) and ground-state (μg) dipole moments were calculated based on the variation of Stokes shift with the solvent's relative permittivity, refractive index, and polarity parameters. SCXRD analyses revealed that the compounds I and II crystallized in the monoclinic system with the space group, P21/n and P21/c, respectively, and with Z = 4 and 2. III, IV, and V crystallized in space groups: orthorhombic, Pbca; triclinic, P-1; and monoclinic, P21 with Z = 1, 2, and 2, respectively. The intermolecular interactions for compounds I-V were investigated using the CCDC Mercury software and their energies were quantified using PIXEL. The density of states (DOS), molecular electrostatic potential surfaces (MEPS), and natural bond orbitals (NBO) of the compounds were determined to evaluate the photophysical properties.

Use of additives to regulate solute aggregation and direct conformational polymorph nucleation of pimelic acid

Understanding the nucleation pathway and achieving regulation to produce the desired crystals are mutually beneficial. The authors previously proposed a nucleation pathway of conformational polymorphs in which solvation and solute self-assembly could affect the result of the conformational rearrangement and further nucleation outcomes. Based on this, herein α,ω-alkanedi-carb-oxy-lic acids (DAn, where n represents the number of carbon atoms in the molecule, n = 2-6, 8-11) were designed as homologous additives to interfere with the self-assembly of pimelic acid (DA7) to further induce the form II compound, which differs from form I only in conformation. Interestingly, longer-chain additives (DA6-11) have a stronger form II-inducing ability than short-chain ones (DA2-4). In addition, an apparent gradient of the degree of interference with solute self-assembly, consistent with form II-inducing ability, was detected by infrared and nuclear magnetic resonance spectroscopy. The calculated molecular electrostatic potential charges also clearly indicate that additive-solute electrostatic interactions gradually increase with increasing carbon chain length of the additives, reaching a maximum value with DA6-11. This novel use of additives demonstrates a direct link between solute aggregation and conformational polymorph nucleation.

The unexpected discovery of the ninth polymorph of tolfenamic acid

A new polymorph of tolfenamic acid, form IX, has been crystallised from a simple cooling crystallisation experiment raising the question as to why this polymorph had never been reported before.

Can molecular flexibility control crystallization? The case of para substituted benzoic acids

Despite the technological importance of crystallization from solutions almost nothing is known about the relationship between the kinetic process of nucleation and the molecular and crystal structures of a crystallizing solute. Nowhere is this more apparent than in our attempts to understand the behavior of increasingly large, flexible molecules developed as active components in the pharmaceutical arena. In our current contribution we develop a general protocol involving a combination of computation (conformation analysis, lattice energy), and experiment (measurement of nucleation rates), and show how significant advances can be made. We present the first systematic study aimed at quantifying the impact of molecular flexibility on nucleation kinetics. The nucleation rates of 4 para substituted benzoic acids are compared, two of which have substituents with flexible chains. In making this comparison, the importance of normalizing data to account for differing solubilities is highlighted. These data have allowed us to go beyond popular qualitative descriptors such ‘crystallizability’ or ‘crystallization propensity’ in favour of more precise nucleation rate data. Overall, this leads to definite conclusions as to the relative importance of solution chemistry, solid-state interactions and conformational flexibility in the crystallization of these molecules and confirms the key role of intermolecular stacking interactions in determining relative nucleation rates. In a more general sense, conclusions are drawn as to conditions under which conformational change may become rate determining during a crystallization process.

Little is known about the relationship between the kinetic process of nucleation and the molecular and crystal structures of a crystallizing solute. Here we compare the behaviour of a series of benzoic acids with a focus on conformational effects.

Transforming Computed Energy Landscapes into Experimental Realities: The Role of Structural Rugosity

We exploit the possible link between structural surface roughness and difficulty of crystallisation. Polymorphs with smooth surfaces may nucleate and crystallise more readily than polymorphs with rough surfaces. The concept is applied to crystal structure prediction landscapes and reveals a promising complementary way of ranking putative crystal structures.

The affinity of copper(ii) ions towards l-amino acids in the solid-state: a simple route towards mixed complexes

Competitive milling was successfully employed to determine the relative affinity of Cu(ii) ions towards selected l-amino acids (Asn, Gln, His, Phe, Pro, and Trp). Described process opens a simple route towards mixed coordination compounds.

Can solvated intermediates inform us about nucleation pathways? The case of β-pABA

Using crystallography to search for nucleation pathways: α and β polymorphs of p-aminobenzoic acid.

Insight into the role of pre-assembly and desolvation in crystal nucleation: a case of p-nitrobenzoic acid

As one of the most important phenomena in crystallization, the crystal nucleation process has always been the focus of research. In this work, influences of pre-assembly species and the desolvation process on the crystal nucleation process were studied. p-Nitrobenzoic acid (PNBA) was taken as a model compound to investigate the relationship between solution chemistry and nucleation kinetics in seven different solvents. One unsolvated form and four solvates of PNBA were obtained and one of the solvates was newly discovered. The nucleation behaviours and nucleation kinetics of PNBA in the seven solvents were studied and analyzed. Density functional theory (DFT) and solvation energy calculation were adopted to evaluate the strength of solute-solvent interactions. Vibrational spectroscopy combined with molecular simulation was applied to reveal the pre-assembly species in the solution. Based on these results, a comprehensive understanding of the relationship between molecular structure, crystal structure, solution chemistry and nucleation dynamics was proposed and discussed. It was found that the structural similarity between solution chemistry and crystal structure, the interaction between specific sites and the overall strength of solvation will jointly affect the nucleation process.

Investigation of an Unusual Crystal Habit of Hydrochlorothiazide Reveals Large Polar Enantiopure Domains and a Possible Crystal Nucleation Mechanism

The observation of an unusual crystal habit in the common diuretic drug hydrochlorothiazide (HCT), and identification of its subtle conformational chirality, has stimulated a detailed investigation of its crystalline forms. Enantiomeric conformers of HCT resolve into an unusual structure of conjoined enantiomorphic twin crystals comprising enantiopure domains of opposite chirality. The purity of the domains and the chiral molecular conformation are confirmed by spatially revolved synchrotron micro-XRD experiments and neutron diffraction, respectively. Macroscopic inversion twin symmetry observed between the crystal wings suggests a pseudoracemic structure that is not a solid solution or a layered crystal structure, but an unusual structural variant of conglomerates and racemic twins. Computed interaction energies for molecular pairs in the racemic and enantiopure polymorphs of HCT, and the observation of large opposing unit-cell dipole moments for the enantiopure domains in these twin crystals, suggest a plausible crystal nucleation mechanism for this unusual crystal habit.

Thermodynamic vs. Kinetic Basis for Polymorph Selection

Ratios of equilibrium solubilities rarely exceed two-fold for polymorph pairs. A model has been developed based on two intrinsic properties of polymorph pairs, namely the ratio of equilibrium solubilities of the individual pairs (C*me/C*st) and the ratio of interfacial energies (γst/γme) and one applied experimental condition, namely the supersaturation identifies which one of a pair of polymorphs nucleates first. A domain diagram has been developed, which identifies the point where the critical free energy of nucleation for the polymorph pair are identical. Essentially, for a system supersaturated with respect to both polymorphs, the model identifies that low supersaturation with respect to the stable polymorph (Sst) leads to an extremely small supersaturation with respect to the metastable polymorph (Sme), radically driving up the critical free energy with respect to the metastable polymorph. Generally, high supersaturations sometimes much higher than the upper limit of the metastable zone, are required to kinetically favour the metastable polymorph.

Comparative Study on the Preparation and Gas Sensing Properties of Reduced Graphene Oxide/SnO2 Binary Nanocomposite for Detection of Acetone in Exhaled Breath

Reduced graphene oxide/tin dioxide (RGO/SnO₂) binary nanocomposite for acetone sensing performance was successfully studied and applied in exhaled breath detection. The influence of structural characteristics was explored by synthesizing the composite (RGO/SnO₂) using the solvothermal method (GS-I) and the hydrothermal method (GS-II) by the chemical route and mechanical mixing, respectively. The nanocomposites characterized by X-ray diffraction (XRD), high-resolution transmission electron microscopy (HRTEM), Fourier transform-infrared spectroscopy (FT-IR), and Brunauer-Emmett-Teller (BET) revealed that GS-I exhibited better surface area, surface energy and showed enhanced gas response than GS-II at an operating temperature of 200 °C. These sensors exhibited comparable response in humid environment as well, suitable for acetone sensing in exhaled breath that clearly distinguishes between healthy and diabetes subjects. The enhanced response at lower concentrations was attributed to the synergistic effect at the RGO/SnO₂ interface. These results indicate that modification in the structural characteristics of RGO/SnO₂ nanocomposite enhances the sensing property. Furthermore, it proved to be a promising material for potential application for point-of care, noninvasive diabetes detection.

Towards understanding π-stacking interactions between non-aromatic rings

The first systematic study of π interactions between non-aromatic rings, based on the authors' own results from an experimental X-ray charge-density analysis assisted by quantum chemical calculations, is presented. The landmark (non-aromatic) examples include quinoid rings, planar radicals and metal-chelate rings. The results can be summarized as: (i) non-aromatic planar polyenic rings can be stacked, (ii) interactions are more pronounced between systems or rings with little or no π-electron delocalization (e.g. quinones) than those involving delocalized systems (e.g. aromatics), and (iii) the main component of the interaction is electrostatic/multipolar between closed-shell rings, whereas (iv) interactions between radicals involve a significant covalent contribution (multicentric bonding). Thus, stacking covers a wide range of interactions and energies, ranging from weak dispersion to unlocalized two-electron multicentric covalent bonding ('pancake bonding'), allowing a face-to-face stacking arrangement in some chemical species (quinone anions). The predominant interaction in a particular stacked system modulates the physical properties and defines a strategy for crystal engineering of functional materials.

Prediction of self-assembly of adenosine analogues in solution: a computational approach validated by isothermal titration calorimetry

The recent discovery of the role of adenosine-analogues as neuroprotectants and cognitive enhancers has sparked interest in these molecules as new therapeutic drugs. Understanding the behavior of these molecules in solution and predicting their ability to self-assemble will accelerate new discoveries. We propose a computational approach based on density functional theory, a polarizable continuum solvation description of the aqueous environment, and an efficient search procedure to probe the potential energy surface, to determine the structure and thermodynamic stability of molecular clusters of adenosine analogues in solution, using caffeine as a model. The method was validated as a tool for the prediction of the impact of small structural variations on self-assembly using paraxanthine. The computational results were supported by isothermal titration calorimetry experiments. The thermodynamic parameters enabled the quantification of the actual percentage of dimer present in solution as a function of concentration. The data suggest that both caffeine and paraxanthine are present at concentrations comparable to the ones found in biological samples.

Polymorphism in p-aminobenzoic acid

We review the polymorphism of p-aminobenzoic acid (pABA), a model drug compound whose crystallisation and polymorphic behaviour has been extensively studied in recent years.

Pre-nucleation aggregation based on solvent microheterogeneity

The microheterogeneous region of aqueous acetonitrile leads to preferred localisation and aggregation of caffeine and theophylline on the interface.

Rationalising crystal nucleation of organic molecules in solution using artificial neural networks

Artificial neural network modelling is used to analyse and predict primary nucleation based on various physicochemical solute and solvent parameters.

Unexpected solvent effects on the UV/Vis absorption spectra of o-cresol in toluene and benzene: in contrast with non-aromatic solvents

Cresol is a prototype molecule in understanding intermolecular interactions in material and biological systems, because it offers different binding sites with various solvents and protonation states under different pH values. It is found that the UV/Vis absorption spectra of o-cresol in aromatic solvents (benzene, toluene) are characterized by a sharp peak, unlike the broad double-peaks in 11 non-aromatic solvents. Both molecular dynamics simulations and electronic structure calculations revealed the formation of intermolecular π-complexation between o-cresol and aromatic solvents. The thermal movements of solvent and solute molecules render the conformations of o-cresol changing between trans and cis isomers. The π-interaction makes the cis configuration a dominant isomer, hence leading to the single keen-edged UV/Vis absorption peak at approximately 283 nm. The free conformation changes between trans and cis in aqueous solution rationalize the broader absorption peaks in the range of 260-280 nm. The pH dependence of the UV/Vis absorption spectra in aqueous solutions is also rationalized by different protonation states of o-cresol. The explicit solvent model with long-ranged interactions is vital to describe the effects of π-complexation and electrostatic interaction on the UV/Vis absorption spectra of o-cresol in toluene and alkaline aqueous (pH > 10.3) solutions, respectively.

Conformational and structural stability of the single molecule and hydrogen bonded clusters of para aminobenzoic acid in the gas and solution phases

The structures of α- and β-para aminobenzoic acid are deconstructed into their hydrogen bonding molecular structural building blocks, where they are analysed usingab initioquantum mechanical calculations of their conformation and cluster stability in solution.

Crystal nucleation and growth in a polymorphic system: Ostwald's rule, p-aminobenzoic acid and nucleation transition states

Water has a magic effect on the nucleation kinetics of an aminobenzoic acid – stacking interactions baffle Ostwald's rule and offer dual pathways.

A basis for the kinetic selection of polymorphs during solution crystallization of organic compounds

Domain diagram for supersaturation needed for a given polymorph pair to select kinetically the metastable or stable forms.