Ab initiocrystal structure prediction. II. Flexible molecules

Role of dispersion interactions in the polymorphism and entropic stabilization of the aspirin crystal

Aspirin has been used and studied for over a century but has only recently been shown to have an additional polymorphic form, known as form II. Since the two observed solid forms of aspirin are degenerate in terms of lattice energy, kinetic effects have been suggested to determine the metastability of the less abundant form II. Here, first-principles calculations provide an alternative explanation based on free-energy differences at room temperature. The explicit consideration of many-body van der Waals interactions in the free energy demonstrates that the stability of the most abundant form of aspirin is due to a subtle coupling between collective electronic fluctuations and quantized lattice vibrations. In addition, a systematic analysis of the elastic properties of the two forms of aspirin rules out mechanical instability of form II as making it metastable.

Complex Polymorphic System of Gallic Acid-Five Monohydrates, Three Anhydrates, and over 20 Solvates

We report the structure of the fifth monohydrate of gallic acid and two additional anhydrate polymorphs and evidence of at least 22 other solvates formed, many containing water and another solvent. This unprecedented number of monohydrate polymorphs and diversity of solid forms is consistent with the anhydrate and monohydrate crystal energy landscapes, showing both a wide range of packing motifs and also some structures differing only in proton positions. By aiding the solution of structures from powder X-ray diffraction data and guiding the screening, the computational studies help explain the complex polymorphism of gallic acid. This is industrially relevant, as the three anhydrates are stable at ambient conditions but hydration/dehydration behavior is very dependent on relative humidity and phase purity.

The polymorphs of ROY: application of a systematic crystal structure prediction technique

We investigate the ability of current ab initio crystal structure prediction techniques to identify the polymorphs of 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecarbonitrile, also known as ROY because of the red, orange and yellow colours of its polymorphs. We use a methodology combining the generation of a large number of structures based on a computationally inexpensive model using the CrystalPredictor global search algorithm, and the further minimization of the most promising of these structures using the CrystalOptimizer local minimization algorithm which employs an accurate, yet efficiently constructed, model based on isolated-molecule quantum-mechanical calculations. We demonstrate that this approach successfully predicts the seven experimentally resolved structures of ROY as lattice-energy minima, with five of these structures being within the 12 lowest energy structures predicted. Some of the other low-energy structures identified are likely candidates for the still unresolved polymorphs of this molecule. The relative stability of the predicted structures only partially matches that of the experimentally resolved polymorphs. The worst case is that of polymorph ON, whose relative energy with respect to Y is overestimated by 6.65 kJ mol(-1). This highlights the need for further developments in the accuracy of the energy calculations.

Towards crystal structure prediction of complex organic compounds--a report on the fifth blind test

Following on from the success of the previous crystal structure prediction blind tests (CSP1999, CSP2001, CSP2004 and CSP2007), a fifth such collaborative project (CSP2010) was organized at the Cambridge Crystallographic Data Centre. A range of methodologies was used by the participating groups in order to evaluate the ability of the current computational methods to predict the crystal structures of the six organic molecules chosen as targets for this blind test. The first four targets, two rigid molecules, one semi-flexible molecule and a 1:1 salt, matched the criteria for the targets from CSP2007, while the last two targets belonged to two new challenging categories - a larger, much more flexible molecule and a hydrate with more than one polymorph. Each group submitted three predictions for each target it attempted. There was at least one successful prediction for each target, and two groups were able to successfully predict the structure of the large flexible molecule as their first place submission. The results show that while not as many groups successfully predicted the structures of the three smallest molecules as in CSP2007, there is now evidence that methodologies such as dispersion-corrected density functional theory (DFT-D) are able to reliably do so. The results also highlight the many challenges posed by more complex systems and show that there are still issues to be overcome.

Which, if any, hydrates will crystallise? Predicting hydrate formation of two dihydroxybenzoic acids

A study of two dihydroxybenzoic acid isomers shows that computational methods can be used to predict hydrate formation, the compound:water ratio and hydrate crystal structures. The calculations also help identify a novel hydrate found in the solid form screening that validates this study.

Successful prediction of a model pharmaceutical in the fifth blind test of crystal structure prediction

The range of target structures in the fifth international blind test of crystal structure prediction was extended to include a highly flexible molecule, (benzyl-(4-(4-methyl-5-(p-tolylsulfonyl)-1,3-thiazol-2-yl)phenyl)carbamate, as a challenge representative of modern pharmaceuticals. Two of the groups participating in the blind test independently predicted the correct structure. The methods they used are described and contrasted, and the implications of the capability to tackle molecules of this complexity are discussed.

Modelling organic crystal structures using distributed multipole and polarizability-based model intermolecular potentials

Crystal structure prediction for organic molecules requires both the fast assessment of thousands to millions of crystal structures and the greatest possible accuracy in their relative energies. We describe a crystal lattice simulation program, DMACRYS, emphasizing the features that make it suitable for use in crystal structure prediction for pharmaceutical molecules using accurate anisotropic atom-atom model intermolecular potentials based on the theory of intermolecular forces. DMACRYS can optimize the lattice energy of a crystal, calculate the second derivative properties, and reduce the symmetry of the spacegroup to move away from a transition state. The calculated terahertz frequency k = 0 rigid-body lattice modes and elastic tensor can be used to estimate free energies. The program uses a distributed multipole electrostatic model (Q, t = 00,...,44s) for the electrostatic fields, and can use anisotropic atom-atom repulsion models, damped isotropic dispersion up to R(-10), as well as a range of empirically fitted isotropic exp-6 atom-atom models with different definitions of atomic types. A new feature is that an accurate model for the induction energy contribution to the lattice energy has been implemented that uses atomic anisotropic dipole polarizability models (alpha, t = (10,10)...(11c,11s)) to evaluate the changes in the molecular charge density induced by the electrostatic field within the crystal. It is demonstrated, using the four polymorphs of the pharmaceutical carbamazepine C(15)H(12)N(2)O, that whilst reproducing crystal structures is relatively easy, calculating the polymorphic energy differences to the accuracy of a few kJ mol(-1) required for applications is very demanding of assumptions made in the modelling. Thus DMACRYS enables the comparison of both known and hypothetical crystal structures as an aid to the development of pharmaceuticals and other speciality organic materials, and provides a tool to develop the modelling of the intermolecular forces involved in molecular recognition processes.

Modeling the interplay of inter- and intramolecular hydrogen bonding in conformational polymorphs

The predicted stability differences of the conformational polymorphs of oxalyl dihydrazide and ortho-acetamidobenzamide are unrealistically large when the modeling of intermolecular energies is solely based on the isolated-molecule charge density, neglecting charge density polarization. Ab initio calculated crystal electron densities showed qualitative differences depending on the spatial arrangement of molecules in the lattice with the greatest variations observed for polymorphs that differ in the extent of inter- and intramolecular hydrogen bonding. We show that accounting for induction dramatically alters the calculated stability order of the polymorphs and reduces their predicted stability differences to be in better agreement with experiment. Given the challenges in modeling conformational polymorphs with marked differences in hydrogen bonding geometries, we performed an extensive periodic density functional study with a range of exchange-correlation functionals using both atomic and plane wave basis sets. Although such electronic structure methods model the electrostatic and polarization contributions well, the underestimation of dispersion interactions by current exchange-correlation functionals limits their applicability. The use of an empirical dispersion-corrected density functional method consistently reduces the structural deviations between the experimental and energy minimized crystal structures and achieves plausible stability differences. Thus, we have established which types of models may give worthwhile relative energies for crystal structures and other condensed phases of flexible molecules with intra- and intermolecular hydrogen bonding capabilities, advancing the possibility of simulation studies on polymorphic pharmaceuticals.