Abstract:The ability to reliably predict the structures and stabilities of a molecular crystal and its polymorphs without any prior experimental information would be an invaluable tool for a number of fields, with specific and immediate applications in the design and formulation of pharmaceuticals. In this case, detailed knowledge of the polymorphic energy landscape for an active pharmaceutical ingredient yields profound insight regarding the existence and likelihood of late-appearing polymorphs. However, the computational prediction of the structures and stabilities of molecular crystal polymorphs is particularly challenging due to the high dimensionality of conformational and crystallographic space accompanied by the need for relative (free) energies to within $\approx$ 1 kJ/mol per molecule. In this work, we combine the most successful crystal structure sampling strategy with the most accurate energy ranking strategy of the latest blind test of organic crystal structure prediction (CSP), organized by the Cambridge Crystallographic Data Centre (CCDC). Our final energy ranking is based on first-principles density functional theory (DFT) calculations that include three key physical contributions: (i) a sophisticated treatment of Pauli exchange-repulsion and electron correlation effects with hybrid functionals, (ii) inclusion of many-body van der Waals dispersion interactions, and (iii) account of vibrational free energies. In doing so, this combined approach has an optimal success rate in producing the crystal structures corresponding to the five blind-test molecules. With this practical approach, we demonstrate the feasibility of obtaining reliable structures and stabilities for molecular crystals of pharmaceutical importance, paving the way towards an enhanced fundamental understanding of polymorphic energy landscapes and routine industrial application of molecular CSP methods.

First Order Methods for Geometric Optimization of Crystal Structures

First Order Methods for Geometric Optimization of Crystals: Experimental Analysis

Scalable Crystal Structure Relaxation Using an Iteration-Free Deep Generative Model with Uncertainty Quantification

Projected gradient descent algorithm for $\textit{ab initio}$ crystal structure relaxation under a fixed unit cell volume

Optimality guarantees for crystal structure prediction

Accelerating crystal structure search through active learning with neural networks for rapid relaxations

Global optimization in the discrete and variable-dimension conformational space: The case of crystal with the strongest atomic cohesion

Latent Conservative Objective Models for Data-Driven Crystal Structure Prediction

Efficient numerical methods for computing the stationary states of phase field crystal models

Mathematical Crystal Chemistry

Gradient Type Optimization Methods for Electronic Structure Calculations

Crystal structure prediction by combining graph network and optimization algorithm

An Improved Crystal Structure Algorithm for Engineering Optimization Problems

Geometric Deep Learning for Molecular Crystal Structure Prediction

A spur to molecular geometry optimization: Gradient-enhanced universal kriging with on-the-fly adaptive ab initio prior mean functions in curvilinear coordinates

Reliable and Practical Computational Prediction of Molecular Crystal Polymorphs

Crystal Structure Prediction via Oblivious Local Search

Efficient ab initio schemes for finding thermodynamically stable and metastable atomic structures: benchmark of cascade genetic algorithms

First-principles stability ranking of molecular crystal polymorphs with the DFT+MBD approach

Reliable crystal structure predictions from first principles

Crystal structure optimization approach to problem solving in mechanical engineering design