Reliable and Practical Computational Prediction of Molecular Crystal Polymorphs

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.

References in corpus (7)

• Resolution-of-identity approach to Hartree-Fock, hybrid density functionals, RPA, MP2, and \textit{GW} with numeric atom-centered orbital basis functions
• Long-range correlation energy calculated from coupled atomic response functions
• Anharmonic vibrational properties in periodic systems: energy, electron-phonon coupling, and stress
• Quantum Monte Carlo study of the phase diagram of solid molecular hydrogen at extreme pressures
• Anharmonic and Quantum Fluctuations in Molecular Crystals: A First-Principles Study of the Stability of Paracetamol
• Anharmonic nuclear motion and the relative stability of hexagonal and cubic ice
• Analytical nuclear gradients for the range-separated many-body dispersion model of noncovalent interactions

Cited by in corpus (19)

• QM7-X: A comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules
• A complete description of thermodynamic stabilities of molecular crystals
• Isotope Effects in Liquid Water via Deep Potential Molecular Dynamics
• Quantum Mechanics of Proteins in Explicit Water: The Role of Plasmon-Like Solute-Solvent Interactions
• Progress and Challenges in Ab Initio Simulations of Quantum Nuclei in Weakly Bonded Systems
• NENCI-2021 Part I: A Large Benchmark Database of Non-Equilibrium Non-Covalent Interactions Emphasizing Close Intermolecular Contacts
• Anharmonic Effects in the Low-Frequency Vibrational Modes of Aspirin and Paracetamol Crystals
• Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics II: Extensions to the Isobaric-Isoenthalpic and Isobaric-Isothermal Ensembles
• Hybrid Classical/Machine-Learning Force Fields for the Accurate Description of Molecular Condensed-Phase Systems
• A Multimer Embedding Approach for Molecular Crystals up to Harmonic Vibrational Properties
• Random Phase Approximation Applied to Many-Body Noncovalent Systems
• Accurate and efficient machine learning interatomic potentials for finite temperature modeling of molecular crystals
• Understanding Discrepancies of Wavefunction Theories for Large Molecules
• Automated Machine Learning Pipeline: Large Language Models-Assisted Automated Dataset Generation for Training Machine-Learned Interatomic Potentials
• Accurate Machine Learning Interatomic Potentials for Polyacene Molecular Crystals: Application to Single Molecule Host-Guest Systems
• Modeling Non-Covalent Interatomic Interactions on a Photonic Quantum Computer
• Efficient Computation of the Long-Range Exact Exchange using an Extended Screening Function
• Coulomb Interactions between Dipolar Quantum Fluctuations in van der Waals Bound Molecules and Materials
• Harmonic-to-anharmonic thermodynamic integration made simple using REG TI