Anharmonic and Quantum Fluctuations in Molecular Crystals: A First-Principles Study of the Stability of Paracetamol

Abstract

Molecular crystals often exist in multiple competing polymorphs, showing significantly different physico-chemical properties. Computational crystal structure prediction is key to interpret and guide the search for the most stable or useful form: A real challenge due to the combinatorial search space, and the complex interplay of subtle effects that work together to determine the relative stability of different structures. Here we take a comprehensive approach based on different flavors of thermodynamic integration in order to estimate all contributions to the free energies of these systems with density-functional theory, including the oft-neglected anharmonic contributions and nuclear quantum effects. We take the two main stable forms of paracetamol as a paradigmatic example. We find that anharmonic contributions, different descriptions of van der Waals interactions, and nuclear quantum effects all matter to quantitatively determine the stability of different phases. Our analysis highlights the many challenges inherent in the development of a quantitative and predictive framework to model molecular crystals. However, it also indicates which of the components of the free energy can benefit from a cancellation of errors that can redeem the predictive power of approximate models, and suggests simple steps that could be taken to improve the reliability of ab initio crystal structure prediction.

References in corpus (13)

• Separable Dual Space Gaussian Pseudo-potentials
• i-PI: A Python interface for ab initio path integral molecular dynamics simulations
• Colored-noise thermostats à la carte
• Efficient first-principles calculation of the quantum kinetic energy and momentum distribution of nuclei
• Interatomic Methods for the Dispersion Energy Derived from the Adiabatic Connection Fluctuation-Dissipation Theorem
• Efficient methods and practical guidelines for simulating isotope effects
• Unraveling quantum mechanical effects in water using isotopic fractionation
• Effect of quantum nuclear motion on hydrogen bonding
• Inverse Temperature Dependence of Nuclear Quantum Effects in DNA Base Pairs
• Accurate molecular dynamics and nuclear quantum effects at low cost by multiple steps in real and imaginary time: using density functional theory to accelerate wavefunction methods
• Stability of Complex Biomolecular Structures: Vander Waals, Hydrogen Bond Cooperativity, and Nuclear Quantum Effects
• Anharmonic nuclear motion and the relative stability of hexagonal and cubic ice
• Impact of vibrational entropy on the stability of unsolvated peptide helices with increasing length

Cited by in corpus (31)

• Nuclear quantum effects enter the mainstream
• i-PI 2.0: A Universal Force Engine for Advanced Molecular Simulations
• Reliable and Practical Computational Prediction of Molecular Crystal Polymorphs
• Computing the Absolute Gibbs Free Energy in Atomistic Simulations: Applications to Defects in Solids
• A complete description of thermodynamic stabilities of molecular crystals
• i-PI 3.0: a flexible and efficient framework for advanced atomistic simulations
• Anharmonic lattice dynamics via the special displacement method
• All-Electron, Real-Space Perturbation Theory for Homogeneous Electric Fields: Theory, Implementation, and Application within DFT
• Hydrogen Diffusion and Trapping in α-Iron: The Role of Quantum and Anharmonic Fluctuations
• Ab initio structural optimization at finite temperatures based on anharmonic phonon theory: Application to the structural phase transitions of BaTiO${}_3$
• How accurate are simulations and experiments for the lattice energies of molecular crystals?
• Progress and Challenges in Ab Initio Simulations of Quantum Nuclei in Weakly Bonded Systems
• Nuclear Quantum Effects in liquid water at near classical computational cost using the adaptive Quantum Thermal Bath
• Anharmonic Effects in the Low-Frequency Vibrational Modes of Aspirin and Paracetamol Crystals
• Finite-temperature materials modeling from the quantum nuclei to the hot electrons regime
• Phonon-induced localization of excitons in molecular crystals from first principles
• The importance of nuclear quantum effects for NMR crystallography
• Quantum vibronic effects on the electronic properties of solid and molecular carbon
• Hybrid Classical/Machine-Learning Force Fields for the Accurate Description of Molecular Condensed-Phase Systems
• Adsorption and Diffusion of H Atoms on beta-PtO2 Surface: The Role of Nuclear Quantum Effects
• A Quantum Computing Implementation of Nuclear-Electronic Orbital (NEO) Theory: Towards an Exact pre-Born-Oppenheimer Formulation of Molecular Quantum Systems
• Accurate and efficient machine learning interatomic potentials for finite temperature modeling of molecular crystals
• A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals
• Anharmonic Vibrational States of Double-Well Potentials in the Solid State from DFT Calculations
• Aqueous Solution Chemistry In Silico and the Role of Data Driven Approaches
• Beyond Gaussian fluctuations of quantum anharmonic nuclei. The case of rotational degrees of freedom
• Acid Base Chemistry of Short Hydrogen Bonds: A Tale of Schrödinger's Cat in Glutamine-Derived Crystals
• Molecular Energy Learning Using Alternative Blackbox Matrix-Matrix Multiplication Algorithm for Exact Gaussian Process
• Accurate Machine Learning Interatomic Potentials for Polyacene Molecular Crystals: Application to Single Molecule Host-Guest Systems
• Harmonic-to-anharmonic thermodynamic integration made simple using REG TI
• A Novel NPT Thermodynamic Integration Scheme to Derive Rigorous Gibbs Free Energies for Crystalline Solids