Molecular mechanism for the cross-nucleation between polymorphs

We use molecular simulations to study the early stages of crystallization in a supercooled liquid of Lennard-Jones particles. We observe the onset of concomitant polymorphism and demonstrate that this phenomenon results from the cross-nucleation of a metastable polymorph on the stable polymorph. We also show that cross-nucleation is selective as it only takes place between polymorphs of almost equivalent free energy. Our simulations provide detailed insights into the molecular mechanism underlying concomitant polymorphism and cross-nucleation between polymorphs.     

Kinetics of cross-nucleation between polymorphs

We report the first kinetic measurement of cross-nucleation between polymorphs, a newly discovered phenomenon important to the theory and control of crystallization. d-Mannitol crystallized from its melt first as the least stable delta polymorph and then as the second least stable alpha polymorph, with alpha nucleating on delta. The kinetics of cross-nucleation was determined from the frequencies of alpha nuclei appearing on delta spherulites, the distances between alpha and delta nuclei, and the growth rate of the delta spherulite. The presence of poly(vinylpyrrolidone), a noncrystallizing, melt-miscible additive, increased the rate of cross-nucleation.     

Polymorph selectivity under nanoscopic confinement

Polymorph selectivity has been achieved during crystallization of anthranilic acid (AA) and 5-methyl-2-[(2-nitrophyenyl)amino]-3-thiophenecarbonitrile (ROY), both considered benchmarks of polymorphic behavior, within nanoporous glass beads and polymer monoliths. Whereas polymorph III of AA crystallizes from the melt on nonporous glass beads or within larger pores, the metastable polymorph II crystallizes in pores with diameters <23 nm, with the selectivity toward this form increasing with decreasing pore size. Of the six ROY polymorphs characterized by single-crystal X-ray diffraction, the yellow form (Y) crystallizes during evaporation of pyridine solutions imbibed by the 30-nm cylindrical pores of porous polycyclohexylethylene (p-PCHE) monoliths. Although both R and ON grow from the melt on the external surfaces of PCHE, only the red form (R) crystallizes in the pores. Amorphous ROY also forms in p-PCHE pores during evaporation from pyridine solutions, subsequently crystallizing to the R nanocrystals upon heating. Although heterogeneous nucleation on the pore walls may play a role, these observations suggest that nucleation and polymorph selectivity is governed by critical size constraints imposed by the ultrasmall pores. The ability to achieve polymorph selectivity in both glass and polymer matrices suggests wide-ranging compatibility with various organic crystalline solids, promising a new approach to controlling polymorphism and searching for unknown polymorphs.     

Crystalline polymorph selection and discovery with polymer heteronuclei

The discovery and selective production of crystalline polymorphs, an outstanding problem in solid-state chemistry, is of great importance industrially in, for example, the manufacture of pharmaceuticals and pigments. Despite considerable efforts, no reliable method exists to produce all of the stable polymorphs of a given compound. Herein, we report methodology to control the phenomenon of crystal polymorphism through the use of diverse libraries of polymer heteronuclei including both commercially available polymers and combinatorially synthesized cross-linked polymers. This new approach for exploring polymorph space offers the advantage of high throughput crystallization to discover multiple polymorphs combined with the ability to selectively produce a given form from a single solvent and temperature condition by simply varying the nature of the polymer substrate. This technique is successfully demonstrated on the pharmaceuticals acetaminophen, sulfamethoxazole, and carbamazepine and on the pharmaceutical intermediate 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecarbonitrile (ROY). High throughput screening, accomplished by optical microscopy and Raman spectroscopy, identified the selective production of the two stable polymorphs of acetaminophen and all six stable forms of ROY. Furthermore, one new form of carbamazepine and two new forms of sulfamethoxazole were discovered; in these cases, single crystals were obtained enabling the structural characterization of two new tetramorphic systems.     

Insights into the molecular mechanism underlying polymorph selection

We use molecular simulations to study polymorph selection during the crystallization of charge-stabilized colloidal suspension. By modifying the conditions of crystallization, we invert the stability of two polymorphs and induce the formation of crystallites whose structure is predominantly that of the stable polymorph. However, our simulations reveal that kinetics play a major role not only during the nucleation step but also in the growth mechanism. The growth of postcritical crystallites of the stable polymorph proceeds through a complex mechanism involving the cross-nucleation of a third metastable polymorph followed by the conversion of this third polymorph into the stable structure.     

New polymorphs of ROY and new record for coexisting polymorphs of solved structures

With six polymorphs coexisting at room temperature, 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecarbonitrile (ROY) is the top system in the current Cambridge Structural Database (Feb. 2005) for the number of polymorphs of solved crystal structures. Here we report two new ROY polymorphs, Y04 and YT04, and the crystal structure of YT04. Y04 is a metastable polymorph that tends to crystallize first from a melt at room temperature, and YT04 is a product of solid-state transformation of Y04. Despite its late discovery, YT04 is the densest among the polymorphs at 25 degrees C and likely the second most stable at 0 K. The conformation of ROY in YT04 is similar to those in the other two yellow polymorphs (Y and YN) but significantly different from those in the orange and red colored polymorphs (ON, OP, ORP, and R). Having escaped years of solution crystallization in several laboratories, Y04 and YT04 exemplify polymorphs that are likely missed by solvent-based screening and discovered through alternative routes.     

How many more polymorphs of ROY remain undiscovered

With 12 crystal forms, 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecabonitrile (a.k.a. ROY) holds the current record for the largest number of fully characterized organic crystal polymorphs. Four of these polymorph structures have been reported since 2019, raising the question of how many more ROY polymorphs await future discovery. Employing crystal structure prediction and accurate energy rankings derived from conformational energy-corrected density functional theory, this study presents the first crystal energy landscape for ROY that agrees well with experiment. The lattice energies suggest that the seven most stable ROY polymorphs (and nine of the twelve lowest-energy forms) on the Z′ = 1 landscape have already been discovered experimentally. Discovering any new polymorphs at ambient pressure will likely require specialized crystallization techniques capable of trapping metastable forms. At pressures above 10 GPa, however, a new crystal form is predicted to become enthalpically more stable than all known polymorphs, suggesting that further high-pressure experiments on ROY may be warranted. This work highlights the value of high-accuracy crystal structure prediction for solid-form screening and demonstrates how pragmatic conformational energy corrections can overcome the limitations of conventional density functionals for conformational polymorphs.

Crystal structure prediction suggests that the low-energy polymorphs of ROY have already been found, but a new high-pressure form is predicted.     

Molecular simulation of cross-nucleation between polymorphs

Recent experiments report that an early nucleating crystalline structure (or polymorph) may nucleate another polymorph. We use molecular dynamics simulations to model this phenomenon known as cross-nucleation. We study the onset of crystallization in a liquid of Lennard-Jones particles cooled at a temperature 22% below the melting temperature. We show that growth proceeds through the successive cross-nucleation of the metastable hexagonal close-packed (hcp) polymorph on the stable face-centered cubic (fcc) polymorph and of the stable fcc polymorph on the metastable hcp polymorph. This finding is in agreement with the experimental results which demonstrated that the cross-nucleation of a stable polymorph on a metastable polymorph is just as likely as the cross-nucleation of a metastable polymorph on a stable polymorph. We then extend our findings established in the case of the homogeneous crystal nucleation to a situation of practical interest, i.e., when a seed of the stable polymorph is used. By studying the crystal growth from the (111) plane of a perfect fcc crystal, we show that, again, growth proceeds through the cross-nucleation of the hcp and fcc structures.     

Effect of self-seed crystal structure on growth of polymorphs in poly(butylene 2,6-naphthalate): A cross-nucleation study

Cross-nucleation between different crystal polymorphs is a particular, self-seed assisted type of heterogeneous nucleation, where a fast-growing polymorph nucleates at a pre-existing crystal surface of another polymorph. Here, we present a study on cross-nucleation between different crystalline phases of poly(butylene 2,6-naphthalate) (PBN), employing hotstage polarized-light optical microscopy and temperature-resolved wide-angle X-ray scattering as analysis tools. PBN forms α-crystals at relatively low temperature and β′-crystals at rather high temperature, with cross-nucleation experiments designed such to first obtain few α- or β′-seed crystals (mother phase) which then are transferred to higher or lower temperature, respectively, to monitor the continuation of the crystallization process and possible growth of the other polymorph. In case of cooling β′-crystals to lower temperature where typically α-crystals form in the non-seeded isotropic melt, β′-crystals nucleate growth of α-crystals, following many examples of cross-nucleation in the literature. In contrast, if low-temperature-generated α-crystals are heated to a temperature where β′-crystals form in a non-seeded melt, the cross-nucleation efficacy is reduced as, beside growth of cross-nucleated β′-crystals, also growth of the mother phase is observed. This unexpected result demonstrates the importance of the structure of the nucleating substrate and the interplay between kinetic and thermodynamic aspects of crystal growth.     

Nucleation of one polymorph by another

Nucleation of one polymorph by another has been observed directly in the melt crystallization of d-mannitol and d-sorbitol. The new polymorph nucleated on an existing one and grew to dominate the end product. This effect is relevant to controlling polymorphism in the manufacture of specialty chemicals and to developing theories of polymorphic nucleation and concomitant polymorphs, for which current efforts emphasize homogeneous nucleation of alternative polymorphs.     

Direct evidence for distinct colour origins in ROY polymorphs

ROY is one of the most well-studied families of crystal structures owing to it being the most polymorphic organic material on record. The various red, orange, and yellow colours of its crystal structures are widely-believed to originate from molecular conformation, though the orange needle (ON) polymorph is thought to be an exception. We report high-pressure, single-crystal X-ray measurements which provide direct experimental evidence that the colour origin in ON is intermolecular, revealing that the molecule undergoes minimal deformation but still exhibits a pronounced, reversible, pale orange → dark red colour change between ambient pressure and 4.18 GPa. Our experimental data are rationalised with band structures, calculated using an accurate hybrid DFT approach, where we are able to account for the variation in colour for five polymorphs of ROY. We highlight the outlier behaviour of ON which shows marked π⋯π stacking interactions that are directly modified through application of pressure. Band structure calculations confirm these intermolecular interactions as the origin of the colour change.

Alternative colour origins in ROY polymorphs are conclusively determined for the first time, using high-pressure diffraction and hybrid DFT.     

Size-dependent growth of polymorphs in nanopores and Ostwald's step rule of stages

More than 100 years after Ostwald postulated his step rule of stages, predictive understanding as to early crystallization stages of polymorphic materials is still premature. We studied crystallization of the polymorphic pharmaceutical acetaminophen in nanoporous glasses as a model for early stages of bulk crystallization since the surface energy significantly contributes to the total Gibbs free energy of nanosized crystals in both cases. Systematic studies of transitions between different polymorphs inside nanoporous glasses show that the thermodynamic stability of the polymorphs depends on the crystal size. Accordingly, the transient occurrence of different polymorphs during crystal growth in bulk systems can be related to surface energy contributions to the total Gibbs free energy of the developing crystals. In nanosized early-stage crystals with high surface-to-volume ratios other polymorphs may be stable than in large crystals with low surface-to-volume ratios. Improved control of the crystallization of polymorphic materials by imposing well-defined confinement is a promising strategy to tailor release of polymorphic drugs and to optimize optical, electronic, magnetic and ferroelectric properties of polymorphic materials.     

Gel-induced selective crystallization of polymorphs

Although nanoporous materials have been explored for controlling crystallization of polymorphs in recent years, polymorphism in confined environments is still poorly understood, particularly from a kinetic perspective, and the role of the local structure of the substrate has largely been neglected. Herein, we report the use of a novel material, polymer microgels with tunable microstructure, for controlling polymorph crystallization from solution and for investigating systematically the effects of nanoconfinement and interfacial interactions on polymorphic outcomes. We show that the polymer microgels can improve polymorph selectivity significantly. The polymorphic outcomes correlate strongly with the gel-induced nucleation kinetics and are very sensitive to both the polymer microstructure and the chemical composition. Further mechanistic investigations suggest that the nucleation-templating effect and the spatial confinement imposed by the polymer network may be central to achieving polymorph selectivity. We demonstrate polymer microgels as promising materials for controlling crystal polymorphism. Moreover, our results help advance the fundamental understanding of polymorph crystallization at complex interfaces, particularly in confined environments.     

Luminescence Color‐Tuning through Polymorph Doping: Preparation of a White‐Emitting Solid from a Single Gold(I)–Isocyanide Complex by Simple Precipitation

We report the luminescent color tuning of a new complex, 2‐benzothiophenyl(4‐methoxyphenyl isocyanide)gold(I) ( 1 ), by using a new “polymorph doping” approach. The slow crystallization of the complex 1 afforded three different pure polymorphic crystals with blue, green, and orange emission under UV‐light irradiation. The crystal structures of pure polymorphs of 1 were investigated in detail by means of single‐crystal X‐ray analyses. Theoretical calculations based on the single‐crystal structures provided qualitative explanation of the difference in the excited energy‐levels of the three polymorphs of 1 . In sharp contrast, the rapid precipitation of 1 , with the optimized conditions reproducibly afforded homogeneous powder materials showing solid‐state white‐emission with Commission Internationale de l’Éclairage (CIE) 1931 chromaticity coordinates of (0.33, 0.35), which is similar to pure white. New “polymorphic doping” has been revealed to be critical to this white emission through spectroscopic and X‐ray diffraction analyses. The coexistence of the multiple polymorphs of 1 within the homogeneous powder materials and the ideal mixing of multiple luminescent colors gave its white emission accompanied with energy transfer from the predominant green‐emitting polymorph to the minor orange‐emitting polymorph.     

Study of dynamics and crystallization kinetics of 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecarbonitrile at ambient and elevated pressure

The organic liquid ROY, i.e., 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecarbonitrile, has been a subject of detailed study in the last few years. One interest in ROY lies in its polymorph-dependent fast crystal growth mode below and above the glass transition temperature. This growth mode is not diffusion controlled, and the possibility that it is enabled by secondary relaxation had been suggested. However, a previous study by dielectric relaxation spectroscopy had not been able to find any resolved secondary relaxation. The present paper reports new dielectric measurements of ROY in the liquid and glassy states at ambient pressure and elevated pressure, which were performed to provide more insight into the molecular dynamics as well as the crystallization tendency of ROY. In the search of secondary relaxation, a special glassy state of ROY was prepared by applying high pressure to the liquid state, from which secondary relaxation was possibly resolved. Thus, the role of secondary relaxation in crystallization of ROY remains to be clarified. Notwithstanding, the secondary relaxation present is not necessarily the sole enabler of crystallization. In an effort to search for possible cause of crystallization other than secondary relaxation, we also performed crystallization kinetics studies of ROY at different T and P combinations while keeping the structural relaxation time constant. The results show that crystallization of ROY speeds up with pressure, opposite to the trend found in the crystallization of ibuprofen studied up to 1 GPa. The dielectric relaxation and thermodynamic properties of ROY with phenolphthalein dimethylether (PDE) are similar in many respects, but PDE does not crystallize. Taking all the above into account, besides the secondary relaxation, the specific chemical structure, molecular interactions and packing of the molecules are additional factors that could affect the kinetics of crystallization found in ROY.     

Inclusion complexes of the natural product gossypol. The formation of different gossypol polymorphs on decomposition of channel type inclusion complexes

The existence of seven gossypol polymorphs has been established. Two of them are obtained by direct crystallization from solution. The remaining five polymorphs are the products of desolvation of channel type complexes (tubulates). Each isostructural group of the complexes on decomposition gives one polymorph. Gossypol thus possesses specific peculiarities in terms of the decomposition of its tubulates, and also the absence of thermotropic polymorphic transitions.     

Measuring free-energy difference between crystal polymorphs through eutectic melting

We describe a method to measure the free-energy difference, DeltaG, between crystal polymorphs from their calorimetric data of eutectic melting with a common additive. The use of different additives yields DeltaG as a function of temperature. The method is suitable for crystals that chemically decompose or physically transform before melting. It applies to not only true polymorphs but also pairs of racemate and conglomerate of resolvable enantiomers. We illustrate the method with the polymorphs of glycine, d-mannitol, and tazofelone and report a new value (123 degrees C) for the enantiotropic transition temperature of alpha and gamma glycine. We show how different additives (including a liquid additive, water) can be used for different compounds. The DeltaG data thus obtained are important for structure-stability studies and controlling crystallization in polymorphic systems.     

Crystalline polymorphism and molecular structure of sodium pravastatin

In this work different crystallization processes of sodium pravastatin are explored and a new polymorph is obtained. The analytical results of powder X-ray diffraction (PXRD) and thermal analysis for this new polymorph indicate that it is different from the polymorphs previously reported. This new crystal form shows different physical-chemical properties than the previous forms, such as crystallographic structure, thermal behavior, and melting point, 181.5 degrees C. Besides, all crystallization processes previously reported use an aprotic solvent as antisolvent. However, we propose a new crystallization process for sodium pravastatin that uses only protic solvents, overcoming industrial scaling and environmental problems. Variable-temperature PXRD experiments show a transformation between different crystal forms in the range of 80-120 degrees C. Solid-state 13C NMR, reported in this work for the first time, and Fourier transform infrared (FT-IR) studies of some polymorphs show some differences in intermolecular interactions, especially with carboxylate and hydroxyl groups. Quantum mechanical calculations of the pravastatin molecule are also presented for the first time, obtaining a molecular structure similar to the experimental structure existing within the crystal lattice of the tert-octylamonium salt of pravastatin.     

Glycine exists mainly as monomers, not dimers, in supersaturated aqueous solutions: implications for understanding its crystallization and polymorphism

Glycine, the simplest amino acid, is described as existing as hydrogen-bonded cyclic dimers in supersaturated aqueous solutions and, as a result, crystallizing in a centrosymmetric polymorph (polymorph alpha) for which the dimer can be viewed as the building unit, in favor of other polymorphs of polar structures. In exhibiting this relation between polymorphic selectivity and self-association in solution, glycine is thought to illustrate a general principle. We measured the freezing-point depression of glycine-water up to 30% supersaturation and found that glycine exists mainly as monomers, not dimers, and that the dimer stability constant K D is smaller than 0.1 kg of H 2O/mol if the observed osmotic abnormality is attributed to dimerization. We also revisited a report cited as evidence for glycine dimerization: the slowdown of glycine diffusion with solution age. Pulsed gradient spin-echo NMR spectroscopy was used in place of the previous method of Gouy interferometry to avoid perturbations to sloution structure caused by the interdiffusion between two solutions of different concentrations. No aging effect was observed on glycine diffusion up to 24% supersaturation after five days. The solute size calculated from diffusivity, viscosity, and the Stokes-Einstein relation showed no increase with concentration or solution age. We conclude that glycine exists in supersaturated aqueous solutions mainly as monomers, not dimers, and remains so upon aging. This result does not invalidate the theories of how pH and additives affect glycine's polymorphic preference, because they begin with the assumption that alpha glycine is the preferred polymorph under normal conditions, but requires a new explanation for that assumption itself.     

Role of liquid polymorphism during the crystallization of silicon

Using molecular simulation, we establish the pivotal role played by liquid polymorphs during the crystallization of silicon. When undercooled at a temperature 20% below the melting point, a silicon melt is under the form of the highly coordinated, high-density liquid (HDL) polymorph. We find that crystallization starts with the formation, within the HDL liquid, of a nanosized droplet of the least stable liquid polymorph, known as the almost tetracoordinated low-density liquid (LDL) polymorph. We then show that the crystalline embryo forms within the LDL droplet, close to the interface with the surrounding HDL liquid, thereby following a pathway associated with a much lower free energy barrier than the direct formation of the crystalline embryo from the HDL liquid would have required. This implies that, for substances exhibiting liquid polymorphs, theories, like the classical nucleation theory, and empirical rules, like Ostwald's rule, should be modified to account for the role of liquid polymorphs in the nucleation process.