Analytical nuclear gradients for the range-separated many-body dispersion model of noncovalent interactions

Abstract

Accurate treatment of the long-range electron correlation energy, including van der Waals (vdW) or dispersion interactions, is essential for describing the structure, dynamics, and function of a wide variety of systems. Among the most accurate models for including dispersion into density functional theory (DFT) is the range-separated many-body dispersion (MBD) method [A. Ambrossetti et al., J. Chem. Phys. 140, 18A508 (2014)], in which the correlation energy is modeled at short-range by a semi-local density functional and at long-range by a model system of coupled quantum harmonic oscillators. In this work, we develop analytical gradients of the MBD energy with respect to nuclear coordinates, including all implicit coordinate dependencies arising from the partitioning of the charge density into Hirshfeld effective volumes. To demonstrate the efficiency and accuracy of these MBD gradients for geometry optimizations of systems with intermolecular and intramolecular interactions, we optimized conformers of the benzene dimer and isolated small peptides with aromatic side-chains. We find excellent agreement with the wavefunction theory reference geometries of these systems (at a fraction of the computational cost) and find that MBD consistently outperforms the popular TS and D3(BJ) dispersion corrections. To demonstrate the performance of the MBD model on a larger system with supramolecular interactions, we optimized the C${}_{60}$@C${}_{60}$H${}_{28}$ buckyball catcher host-guest complex. Finally, we find that neglecting the implicit nuclear coordinate dependence arising from the charge density partitioning, as has been done in prior numerical treatments, leads to an unacceptable error in the MBD forces, with relative errors of ~20% (on average) that can extend well beyond 100%.

18 pages including 5 figures; plus 9 pages supplementary information including 5 additional figures

References in corpus (7)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• A Higher-Accuracy van der Waals Density Functional
• Nonlocal van der Waals density functional: The simpler the better
• Van der Waals bonding in layered compounds from advanced first-principles calculations
• Perspective: Advances and challenges in treating van der Waals dispersion forces in density functional theory
• Van der Waals interactions in DFT made easy by Wannier functions
• Analytical nuclear gradients for the range-separated many-body dispersion model of noncovalent interactions

Cited by in corpus (14)

• Advanced capabilities for materials modelling with Quantum ESPRESSO
• Reliable and Practical Computational Prediction of Molecular Crystal Polymorphs
• Non-covalent interactions across organic and biological subsets of chemical space: Physics-based potentials parametrized from machine learning
• Isotope Effects in Liquid Water via Deep Potential Molecular Dynamics
• Equation of state of fluid methane from first principles with machine learning potentials
• Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics I: Theory, Algorithm, and Performance
• Communication: Charge-Population Based Dispersion Interactions for Molecules and Materials
• Analytical nuclear gradients for the range-separated many-body dispersion model of noncovalent interactions
• libMBD: A general-purpose package for scalable quantum many-body dispersion calculations
• Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics II: Extensions to the Isobaric-Isoenthalpic and Isobaric-Isothermal Ensembles
• Thermal Expansion in Dispersion-Bound Molecular Crystals
• Diffusion Monte Carlo Study on Relative Stabilities of Boron Nitride Polymorphs
• Reproducibility of Potential Energy Surfaces of Organic/Metal Interfaces on the Example of PTCDA on Ag(111)
• Smooth Particle Mesh Ewald-integrated stochastic Lanczos Many-body Dispersion algorithm