Interactions between Large Molecules: Puzzle for Reference Quantum-Mechanical Methods

Abstract

Quantum-mechanical methods are widely used for understanding molecular interactions throughout biology, chemistry, and materials science. Quantum diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] are two state-of-the-art and trusted wavefunction methods that have been categorically shown to yield accurate interaction energies for small organic molecules. These methods provide valuable reference information for widely-used semi-empirical and machine learning potentials, especially where experimental information is scarce. However, agreement for systems beyond small molecules is a crucial remaining milestone for cementing the benchmark accuracy of these methods. Approaching such well-converged predictive power in larger molecules has motivated major developments in CCSD(T) as well as DMC algorithms in the past years, resulting in orders of magnitude time-to-solution reductions. Here, we show that CCSD(T) and DMC interaction energies are not in consistent agreement for a set of polarizable supramolecules. Whilst agreement is found for some of the complexes, in a few key systems disagreements of up to 8 kcal/mol remained. This leads to differences of up to 6 orders of magnitude in the corresponding binding association constant at room temperature for systems which are well within the accustomed domain of applicability for both methods. These findings thus indicate that more caution is required when aiming at reproducible non-covalent interactions between extended molecules. Our data contradicts the expectation that the most comprehensive and robust wavefunction methods predict identical non-covalent interactions and indicate an unsolved challenge for benchmark approaches.

References in corpus (14)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Continuum variational and diffusion quantum Monte Carlo calculations
• Weak binding between two aromatic rings: feeling the van der Waals attraction by quantum Monte Carlo methods
• Beyond the locality approximation in the standard diffusion Monte Carlo method
• Smooth relativistic Hartree-Fock pseudopotentials for H to Ba and Lu to Hg
• Norm-conserving Hartree-Fock pseudopotentials and their asymptotic behavior
• Variational and Diffusion Quantum Monte Carlo Calculations with the CASINO Code
• Multi-Determinant Wave-functions in Quantum Monte Carlo
• TurboRVB: a many-body toolkit for {\it ab initio} electronic simulations by quantum Monte Carlo
• Shape and Energy Consistent Pseudopotentials for Correlated Electron systems
• The S66 noncovalent interactions benchmark reconsidered using explicitly correlated methods near the basis set limit
• Toward Accurate Adsorption Energetics on Clay Surfaces
• A comparative study using state-of-the-art electronic structure theories on solid hydrogen phases under high pressures
• Assessing the Accuracy of the Jastrow Antisymmetrized Geminal Power in the H4 Model System

Cited by in corpus (30)

• r2SCAN-D4: Dispersion corrected meta-generalized gradient approximation for general chemical applications
• The Modern Mathematics of Deep Learning
• Measurements as a roadblock to near-term practical quantum advantage in chemistry: resource analysis
• Towards the ground state of molecules via diffusion Monte Carlo on neural networks
• S66x8 Noncovalent Interactions Revisited: New Benchmark and Performance of Composite Localized Coupled-Cluster Methods
• Averting the infrared catastrophe in the gold standard of quantum chemistry
• Towards Pair Atomic Density Fitting for Correlation Energies with Benchmark Accuracy
• libMBD: A general-purpose package for scalable quantum many-body dispersion calculations
• Linear-scaling local natural orbital CCSD(T) approach for open-shell systems: algorithm, benchmarks, and large-scale applications
• Hybrid Classical/Machine-Learning Force Fields for the Accurate Description of Molecular Condensed-Phase Systems
• Quantum computing quantum Monte Carlo with hybrid tensor network for electronic structure calculations
• Performance of Localized-Orbital Coupled Cluster Approaches for the Conformational Energies of Longer n-alkane Chains
• Diffusion Monte Carlo Study on Relative Stabilities of Boron Nitride Polymorphs
• Another Angle on Benchmarking Noncovalent Interactions
• On-Premises Superconducting Quantum Computer for Education and Research
• Improved modularity and new features in ipie: Toward even larger AFQMC calculations on CPUs and GPUs at zero and finite temperatures
• Post-CCSD(T) corrections in the S66 noncovalent interactions benchmark
• Systematic discrepancies between reference methods for non-covalent interactions within the S66 dataset
• Understanding Discrepancies of Wavefunction Theories for Large Molecules
• Phaseless auxiliary field quantum Monte Carlo with projector-augmented wave method for solids
• Basis set incompleteness errors in fixed-node diffusion Monte Carlo calculations on non-covalent interactions
• Approaching the basis-set limit of the dRPA correlation energy with explicitly correlated and Projector Augmented-wave methods
• Reproducibility of fixed-node diffusion Monte Carlo across diverse community codes: The case of water-methane dimer
• Self-Refinement of Auxiliary-Field Quantum Monte Carlo via Non-Orthogonal Configuration Interaction
• The S66 Noncovalent Interaction Benchmark Re-examined: Composite Localized Coupled Cluster Approaches
• Efficient Implementation of the Random Phase Approximation with Domain-based Local Pair Natural Orbitals
• Quantum computing quantum Monte Carlo algorithm
• Beyond single-reference fixed-node approximation in ab initio Diffusion Monte Carlo using antisymmetrized geminal power applied to systems with hundreds of electrons
• A quantum Monte Carlo study of systems with effective core potentials and node nonlinearities
• Accurate Helium-Benzene Potential: from CCSD(T) to Gaussian Process Regression