From crystal structure prediction to polymorph prediction: interpreting the crystal energy landscape.

Many organic molecules are emerging as having many crystalline forms, including polymorphs and solvates, as more techniques are being used to generate and characterise the organic solid state. The fundamental scientific and industrial interest in controlling crystallisation is inspiring the development of computational methods of predicting which crystal structures are thermodynamically feasible. Sometimes, computing this crystal energy landscape will reveal that a molecule has one way of packing with itself that is sufficiently more favourable than any other so only this crystal structure will be observed. More frequently, there will be many energy minima that are energetically feasible, showing approximately equi-energetic compromises between the various intermolecular interactions allowed by the conformational flexibility. Such cases generally lead to multiple solid forms. At the moment, we usually calculate the lattice energy landscape, an approximation to the real crystal energy landscape at 0 K. Despite its limitations, many studies show that this is a valuable complement to solid form screening, which helps in discovering new structures as well as rationalising the solid forms that are found in experimental searches. The range of factors that can determine which of the thermodynamically feasible crystal structures are observed polymorphs, shows the many further challenges in developing crystal energy landscapes as a tool for control of the organic solid state.

[1]  Dario Braga,et al.  Crystal Engineering: From Molecules and Crystals to Materials , 1999 .

[2]  K. Shankland,et al.  Control and prediction of packing motifs: a rare occurrence of carbamazepine in a catemeric configuration , 2006 .

[3]  A. Gavezzotti,et al.  Polymorphic Forms of Organic Crystals at Room Conditions: Thermodynamic and Structural Implications , 1995 .

[4]  J. Steed,et al.  Comment on “On the presence of multiple molecules in the crystal asymmetric unit (Z′ > 1)” by Gautam R. Desiraju, CrystEngComm, 2007, 9, 91 , 2007 .

[5]  Lian Yu.  Survival of the fittest polymorph: how fast nucleater can lose to fast grower , 2007 .

[6]  G. Desiraju,et al.  Tautomeric polymorphism in omeprazole. , 2007, Chemical communications.

[7]  Harold A Scheraga,et al.  Exercises in prognostication: crystal structures and protein folding. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[8]  G. Desiraju,et al.  On the polymorphism of aspirin: crystalline aspirin as intergrowths of two "polymorphic" domains. , 2007, Angewandte Chemie.

[9]  M. Zaworotko,et al.  From molecules to crystal engineering: supramolecular isomerism and polymorphism in network solids. , 2001, Chemical reviews.

[10]  G. P. Stahly.  Diversity in Single- and Multiple-Component Crystals. The Search for and Prevalence of Polymorphs and Cocrystals , 2007 .

[11]  T. C. Lewis,et al.  The polymorphism of progesterone: stabilization of a 'disappearing' polymorph by co-crystallization. , 2007, Journal of pharmaceutical sciences.

[12]  F. Allen.  The Cambridge Structural Database: a quarter of a million crystal structures and rising. , 2002, Acta crystallographica. Section B, Structural science.

[13]  Sarah L. Price,et al.  Toward Crystal Structure Prediction for Conformationally Flexible Molecules: The Headaches Illustrated by Aspirin , 2004 .

[14]  S. Price,et al.  A Systematic Experimental and Theoretical Study of the Crystalline State of Six Chloronitrobenzenes , 2008 .

[15]  T. C. Lewis,et al.  Investigating Unused Hydrogen Bond Acceptors Using Known and Hypothetical Crystal Polymorphism , 2005 .

[16]  W. Jiskoot,et al.  Selective Production of Polymorphs and Pseudomorphs Using Supercritical Fluid Crystallization from Aqueous Solutions , 2007 .

[17]  P Verwer,et al.  A test of crystal structure prediction of small organic molecules. , 2000, Acta crystallographica. Section B, Structural science.

[18]  Sarah L. Price,et al.  Energy minimization of crystal structures containing flexible molecules , 2006 .

[19]  P. Karamertzanis,et al.  Molecular conformations and relative stabilities can be as demanding of the electronic structure method as intermolecular calculations. , 2006, The journal of physical chemistry. A.

[20]  Alan R. Kennedy,et al.  Targeted crystallisation of novel carbamazepine solvates based on a retrospective Random Forest classification , 2008 .

[21]  C. Bugg,et al.  Relationship between the mutagenic and base-stacking properties of halogenated uracil derivatives. The crystal structures of 5-chloro- and 5-bromouracil. , 1975, Biochimica et biophysica acta.

[22]  K. Shankland,et al.  Single-crystal X-ray and neutron powder diffraction investigation of the phase transition in tetrachlorobenzene. , 2006, Acta crystallographica. Section B, Structural science.

[23]  G. Day,et al.  Importance of Molecular Shape for the Overall Stability of Hydrogen Bond Motifs in the Crystal Structures of Various Carbamazepine-Type Drug Molecules , 2007 .

[24]  A. Gavezzotti.  A solid-state chemist's view of the crystal polymorphism of organic compounds. , 2007, Journal of pharmaceutical sciences.

[25]  A. Matzger,et al.  Crystalline polymorph selection and discovery with polymer heteronuclei. , 2005, Journal of the American Chemical Society.

[26]  Peddy Vishweshwar,et al.  The predictably elusive form II of aspirin. , 2005, Journal of the American Chemical Society.

[27]  J. Bernstein,et al.  Graph-set analysis of hydrogen-bond patterns in organic crystals. , 1990, Acta crystallographica. Section B, Structural science.

[28]  G. Day,et al.  Atomistic calculations of phonon frequencies and thermodynamic quantities for crystals of rigid organic molecules , 2003 .

[29]  Sarah L Price,et al.  Challenges of crystal structure prediction of diastereomeric salt pairs. , 2005, The journal of physical chemistry. B.

[30]  William Jones,et al.  Beyond the isotropic atom model in crystal structure prediction of rigid molecules: atomic multipoles versus point charges , 2005 .

[31]  A. Myerson,et al.  Supersaturation and Polarization Dependence of Polymorph Control in the Nonphotochemical Laser-Induced Nucleation (NPLIN) of Aqueous Glycine Solutions , 2006 .

[32]  K. Harris,et al.  Alteration of polymorphic selectivity through different crystallization mechanisms occurring in the same crystallization solution. , 2007, The journal of physical chemistry. B.

[33]  S. Price,et al.  A new polymorph of 5-fluorouracil found following computational crystal structure predictions. , 2005, Journal of the American Chemical Society.

[34]  A. Stone,et al.  Atom–atom potentials from ab initio calculations , 2007 .

[35]  A. Gavezzotti,et al.  Molecular Aggregation: Structure Analysis and Molecular Simulation of Crystals and Liquids , 2007 .

[36]  J. Anwar,et al.  Computer Simulation of Crystallization from Solution , 1998 .

[37]  J. Bernstein,et al.  Groth's original concomitant polymorphs revisited , 2002 .

[38]  A. Gavezzotti,et al.  Are Crystal Structures Predictable , 1994 .

[39]  A. Gavezzotti,et al.  Ten years of experience in polymorph prediction: what next? , 2002 .

[40]  Liang Chen,et al.  Grid Service Orchestration Using the Business Process Execution Language (BPEL) , 2005, Journal of Grid Computing.

[41]  G. Desiraju,et al.  On the polymorphism of aspirin. , 2007, Angewandte Chemie.

[42]  Marc-Antoine Perrin,et al.  Energy ranking of molecular crystals using density functional theory calculations and an empirical van der waals correction. , 2005, The journal of physical chemistry. B.

[43]  Anthony J Stone,et al.  Distributed Multipole Analysis:  Stability for Large Basis Sets. , 2005, Journal of chemical theory and computation.

[44]  Graeme M. Day,et al.  Elastic Constant Calculations for Molecular Organic Crystals , 2001 .

[45]  P. W. Cains,et al.  Toward the computational design of diastereomeric resolving agents: an experimental and computational study of 1-phenylethylammonium-2-phenylacetate derivatives. , 2007, The journal of physical chemistry. B.

[46]  Jacco van de Streek,et al.  Searching the Cambridge Structural Database for the `best' representative of each unique polymorph , 2006 .

[47]  Gautam R. Desiraju,et al.  On the presence of multiple molecules in the crystal asymmetric unit ( Z ′ > 1) , 2007 .

[48]  Roger J. Davey,et al.  Concerning the Relationship between Structural and Growth Synthons in Crystal Nucleation: Solution and Crystal Chemistry of Carboxylic Acids As Revealed through IR Spectroscopy , 2006 .

[49]  R. Davey,et al.  A Whole Output Strategy for Polymorph Screening: Combining Crystal Structure Prediction, Graph Set Analysis, and Targeted Crystallization Experiments in the Case of Diflunisal , 2003 .

[50]  Sarah L Price,et al.  Toward the Prediction of Organic Hydrate Crystal Structures. , 2007, Journal of chemical theory and computation.

[51]  K. Shankland,et al.  Crystallization and crystal energy landscape of hydrochlorothiazide , 2007 .

[52]  Gordon J. T. Tiddy,et al.  Crystal engineering – nucleation, the key step , 2002 .

[53]  S. Price,et al.  An automated parallel crystallisation search for predicted crystal structures and packing motifs of carbamazepine. , 2006, Journal of pharmaceutical sciences.

[54]  K. Harris,et al.  A multi-technique approach for probing the evolution of structural properties during crystallization of organic materials from solution. , 2007, Faraday discussions.

[55]  Reginald B. H. Tan,et al.  Stable polymorphs: difficult to make and difficult to predict , 2007 .

[56]  M. Spackman,et al.  Solvent inclusion in the structural voids of form II carbamazepine: single-crystal X-ray diffraction, NMR spectroscopy and Hirshfeld surface analysis , 2007 .

[57]  K. Shankland,et al.  An automated platform for parallel crystallization of small organic molecules , 2006 .

[58]  Christopher S Towler,et al.  Impact of molecular speciation on crystal nucleation in polymorphic systems: the conundrum of gamma glycine and molecular 'self poisoning'. , 2004, Journal of the American Chemical Society.

[59]  T. C. Lewis,et al.  Computational prediction and X-ray determination of the crystal structures of 3-oxauracil and 5-hydroxyuracil—an informal blind test , 2005 .

[60]  D. Allan,et al.  An exploration of the polymorphism of piracetam using high pressure , 2005 .

[61]  Enrico Drioli,et al.  Selective Glycine Polymorph Crystallization by Using Microporous Membranes , 2007 .

[62]  Alfred Y Lee,et al.  Factors Affecting the Polymorphic Outcome of Glycine Crystals Constrained on Patterned Substrates , 2006 .

[63]  Michele Parrinello,et al.  Freezing of a Lennard-Jones fluid: from nucleation to spinodal regime. , 2006, Physical review letters.

[64]  Joel Bernstein,et al.  Polymorphism in Molecular Crystals , 2002 .

[65]  Michael J. Cima,et al.  Elucidation of crystal form diversity of the HIV protease inhibitor ritonavir by high-throughput crystallization , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[66]  Frank J. J. Leusen,et al.  A study of different approaches to the electrostatic interaction in force field methods for organic crystals , 2003 .

[67]  B. Kariuki,et al.  Characterization of complicated new polymorphs of chlorothalonil by X-ray diffraction and computer crystal structure prediction. , 2004, Journal of the American Chemical Society.

[68]  C. Pantelides,et al.  Ab initio crystal structure prediction. II. Flexible molecules , 2007 .

[69]  C. A. Mitchell,et al.  Selective nucleation and discovery of organic polymorphs through epitaxy with single crystal substrates. , 2001, Journal of the American Chemical Society.

[70]  Sarah L. Price,et al.  Quantifying intermolecular interactions and their use in computational crystal structure prediction , 2004 .

[71]  T. Threlfall.  Crystallisation of Polymorphs: Thermodynamic Insight into the Role of Solvent , 2000 .

[72]  G. Day,et al.  Dynamics in crystals of rigid organic molecules: contrasting the phonon frequencies calculated by molecular dynamics with harmonic lattice dynamics for imidazole and 5-azauracil , 2004 .

[73]  J. Gale,et al.  An approach to developing a force field for molecular simulation of martensitic phase transitions between phases with subtle differences in energy and structure. , 2004, Journal of the American Chemical Society.

[74]  K. Shankland,et al.  Search for a predicted hydrogen bonding motif--a multidisciplinary investigation into the polymorphism of 3-azabicyclo[3.3.1]nonane-2,4-dione. , 2007, Journal of the American Chemical Society.

[75]  C. Catlow,et al.  Kinetic insights into the role of the solvent in the polymorphism of 5-fluorouracil from molecular dynamics simulations. , 2006, The journal of physical chemistry. B.

[76]  I. Weissbuch,et al.  Solvent effect on crystal polymorphism: why addition of methanol or ethanol to aqueous solutions induces the precipitation of the least stable beta form of glycine. , 2005, Angewandte Chemie.

[77]  Gautam R. Desiraju,et al.  Crystal engineering : the design of organic solids , 1989 .

[78]  T. C. Lewis,et al.  A third blind test of crystal structure prediction. , 2005, Acta crystallographica. Section B, Structural science.

[79]  Chick C. Wilson,et al.  Comparing entire crystal structures: structural genetic fingerprinting , 2007 .

[80]  Michele Parrinello,et al.  Exploring polymorphism : The case of benzene , 2005 .

[81]  Bouke P. van Eijck,et al.  Ab initio crystal structure predictions for flexible hydrogen-bonded molecules. Part III. Effect of lattice vibrations , 2001, J. Comput. Chem..

[82]  K. Shankland,et al.  A catemer-to-dimer structural transformation in cyheptamide , 2008 .

[83]  K. Shankland,et al.  Discovery of three polymorphs of 7-fluoroisatin reveals challenges in using computational crystal structure prediction as a complement to experimental screening , 2008 .

[84]  Jukka Rantanen,et al.  Crystallization of glycine with ultrasound. , 2006, International journal of pharmaceutics.

[85]  G. Day,et al.  Solvent inclusion in form II carbamazepine. , 2007, Chemical communications.

[86]  M. Alderton,et al.  Distributed multipole analysis , 2006 .

[87]  G. Desiraju,et al.  Searching for a polymorph: second crystal form of 6-amino-2-phenylsulfonylimino-1,2-dihydropyridine. , 2003, Angewandte Chemie.

[88]  Shinji Saito,et al.  Molecular dynamics simulation of the ice nucleation and growth process leading to water freezing , 2002, Nature.

[89]  Sarah L Price,et al.  A nonempirical anisotropic atom-atom model potential for chlorobenzene crystals. , 2003, Journal of the American Chemical Society.

[90]  A. Gavezzotti,et al.  A Statistical Study of Density and Packing Variations among Crystalline Isomers , 2000 .

[91]  J. Bernstein,et al.  Serendipity and four polymorphic structures of benzidine, C12H12N2. , 2006, Journal of the American Chemical Society.

[92]  John B. O. Mitchell,et al.  Can we predict lattice energy from molecular structure? , 2003, Acta Crystallographica Section B Structural Science.

[93]  W. David,et al.  Effect of high pressure on the crystal structures of polymorphs of glycine , 2005 .

[94]  Victor E Bazterra,et al.  A Distributed Computing Method for Crystal Structure Prediction of Flexible Molecules:  An Application to N-(2-Dimethyl-4,5-dinitrophenyl) Acetamide. , 2007, Journal of chemical theory and computation.

[95]  S. Chemburkar,et al.  Dealing with the Impact of Ritonavir Polymorphs on the Late Stages of Bulk Drug Process Development , 2000 .

[96]  J. Dunitz.  Are crystal structures predictable? , 2003, Chemical communications.

[97]  T. Threlfall,et al.  Structural and thermodynamic explanations of Ostwald's rule , 2003 .

[98]  M. Schmidt,et al.  Determination of the structure of the violet pigment C22H12Cl2N6O4 from a non-indexed X-ray powder diagram. , 2005, Acta Crystallographica Section B Structural Science.

[99]  Robin Taylor,et al.  Mercury: visualization and analysis of crystal structures , 2006 .

[100]  C. Hunter,et al.  Exciton coupling in porphyrin dimers , 1989 .

[101]  Michael B. Hursthouse,et al.  A versatile procedure for the identification, description and quantification of structural similarity in molecular crystals , 2005 .

[102]  Ian Brown,et al.  Providing an effective data infrastructure for the simulation of complex materials , 2006 .

[103]  A. Gavezzotti,et al.  Molecular recognition in organic crystals: directed intermolecular bonds or nonlocalized bonding? , 2005, Angewandte Chemie.

[104]  Sarah L Price,et al.  Crystal structure prediction of small organic molecules: a second blind test. , 2002, Acta crystallographica. Section B, Structural science.

[105]  Christer B. Aakeröy,et al.  Crystal engineering : Strategies and architectures , 1997 .

[106]  G. Day,et al.  A study of the known and hypothetical crystal structures of pyridine: why are there four molecules in the asymmetric unit cell? , 2002 .

[107]  D. Tocher,et al.  The Discovery of New Crystal Forms of 5-Fluorocytosine Consistent with the Results of Computational Crystal Structure Prediction , 2006 .

[108]  Jeffrey S. Tan,et al.  Capillary Precipitation of a Highly Polymorphic Organic Compound , 2003 .

[109]  S. Price,et al.  Validation of a search technique for crystal structure prediction of flexible molecules by application to piracetam. , 2005, Acta crystallographica. Section B, Structural science.