Noncovalent Interactions by Quantum Monte Carlo

Beyond Single-Reference Fixed-Node Approximation in Ab Initio Diffusion Monte Carlo Using Antisymmetrized Geminal Power Applied to Systems with Hundreds of Electrons

Diffusion Monte Carlo (DMC) is an exact technique to project out the ground state (GS) of a Hamiltonian. Since the GS is always bosonic, in Fermionic systems, the projection needs to be carried out while imposing antisymmetric constraints, which is a nondeterministic polynomial hard problem. In practice, therefore, the application of DMC on electronic structure problems is made by employing the fixed-node (FN) approximation, consisting of performing DMC with the constraint of having a fixed, predefined nodal surface. How do we get the nodal surface? The typical approach, applied in systems having up to hundreds or even thousands of electrons, is to obtain the nodal surface from a preliminary mean-field approach (typically, a density functional theory calculation) used to obtain a single Slater determinant. This is known as single reference. In this paper, we propose a new approach, applicable to systems as large as the C60 fullerene, which improves the nodes by going beyond the single reference. In practice, we employ an implicitly multireference ansatz (antisymmetrized geminal power wave function constraint with molecular orbitals), initialized on the preliminary mean-field approach, which is relaxed by optimizing a few parameters of the wave function determining the nodal surface by minimizing the FN-DMC energy. We highlight the improvements of the proposed approach over the standard single-reference method on several examples and, where feasible, the computational gain over the standard multireference ansatz, which makes the methods applicable to large systems. We also show that physical properties relying on relative energies, such as binding energies, are affordable and reliable within the proposed scheme.

Diffusion Monte Carlo method for barrier heights of multiple proton exchanges and complexation energies in small water, ammonia, and hydrogen fluoride clusters

Proton exchange reactions are of key importance in many processes in water. However, it is nontrivial to achieve reliable barrier heights for multiple proton exchanges and complexation energies in hydrogen-bonded systems theoretically. Performance of the fixed-node diffusion quantum Monte Carlo (FN-DMC) with the single-Slater-Jastrow trial wavefunction on total energies, barrier heights of multiple proton exchanges, and complexation energies of small water, ammonia, and hydrogen fluoride clusters is investigated in this study. Effects of basis sets and those of locality approximation (LA), T-move approximation (T-move), and determinant localization approximation (DLA) schemes in dealing with the nonlocal part of pseudopotentials on FN-DMC results are evaluated. According to our results, diffuse basis functions are important in achieving reliable barrier heights and complexation energies with FN-DMC, although the cardinal number of the basis set is more important than diffuse basis functions on total energies of these systems. Our results also show that the time step bias with DLA and LA is smaller than T-move; however, the time step bias of DMC energies with respect to time steps using the T-move is roughly linear up to 0.06 a.u., while this is not the case with LA and DLA. Barrier heights and complexation energies with FN-DMC using these three schemes are always within chemical accuracy. Taking into account the fact that T-move and DLA are typically more stable than LA, FN-DMC calculations with the T-move or DLA scheme and basis sets containing diffuse basis functions are suggested for barrier heights of multiple proton exchanges and complexation energies of hydrogen-bonded clusters.

Toward accurate modeling of structure and energetics of bulk hexagonal boron nitride

Materials that exhibit both strong covalent and weak van der Waals interactions pose a considerable challenge to many computational methods, such as DFT. This makes assessing the accuracy of calculated properties, such as exfoliation energies in layered materials like hexagonal boron nitride (h-BN) problematic, when experimental data are not available. In this paper, we investigate the accuracy of equilibrium lattice constants and exfoliation energy calculation for various DFT-based computational approaches in bulk h-BN. We contrast these results with available experiments and reference fixed-node diffusion quantum Monte Carlo (QMC) results. From our reference QMC calculation, we obtained an exfoliation energy of - 33 ± 2 meV/atom (-0.38 ± 0.02 J/m2 ).

Molecules in Environments: Toward Systematic Quantum Embedding of Electrons and Drude Oscillators

We develop a quantum embedding method that enables accurate and efficient treatment of interactions between molecules and an environment, while explicitly including many-body correlations. The molecule is composed of classical nuclei and quantum electrons, whereas the environment is modeled via charged quantum harmonic oscillators. We construct a general Hamiltonian and introduce a variational Ansatz for the correlated ground state of the fully interacting molecule-environment system. This wave function is optimized via the variational Monte Carlo method and the ground state energy is subsequently estimated through the diffusion Monte Carlo method. The proposed scheme allows an explicit many-body treatment of electrostatic, polarization, and dispersion interactions between the molecule and the environment. We study solvation energies and excitation energies of benzene derivatives, obtaining excellent agreement with explicit ab initio calculations and experiments.

Dynamic and Nondynamic Electron Correlation Energy Decomposition Based on the Node of the Hartree-Fock Slater Determinant

Distinguishing between dynamic and nondynamic electron correlation energy is a fundamental concept in quantum chemistry. It can be challenging to make a clear distinction between the two types of correlation energy or to determine their actual contributions in specific cases using wave function theory. This is because both single-reference and multireference methods cover both types of correlation energy to some extent. Fixed-node diffusion quantum Monte Carlo (FNDMC) accurately covers dynamic correlations, but it is limited in overall accuracy by the node of the trial wave function. We introduce a methodology for partitioning an exact electron correlation energy into its dynamic and nondynamic components. This is accomplished by restricting a ground-state solution from sharing its node with a spin-restricted Hartree-Fock Slater determinant. The FNDMC method is used as a tool to conveniently project out a lowest-energy state obeying such a boundary condition. The proposed approach provides an unambiguous and useful procedure for separating electron correlation energy, as demonstrated on multiple systems, including the He atom, bond breaking of H2, the parametric H2-H2 system, the Be-Ne atomic series with low- and high-spin states for C, N, and O atoms, and small molecules such as BH, HF, and CO at both equilibrium and elongated configurations, respectively.

Regularized and Opposite Spin-Scaled Functionals from Møller-Plesset Adiabatic Connection─Higher Accuracy at Lower Cost

Noncovalent interactions (NCIs) play a crucial role in biology, chemistry, material science, and everything in between. To improve pure quantum-chemical simulations of NCIs, we propose a methodology for constructing approximate correlation energies by combining an interpolation along the Møller-Plesset adiabatic connection (MP AC) with a regularization and spin-scaling strategy applied to MP2 correlation energies. This combination yields cosκos-SPL2, which exhibits superior accuracy for NCIs compared to any of the individual strategies. With the N4 formal scaling, cosκos-SPL2 is competitive or often outperforms more expensive dispersion-corrected double hybrids for NCIs. The accuracy of cosκos-SPL2 particularly shines for anionic halogen bonded complexes, where it surpasses standard dispersion-corrected DFT by a factor of 3 to 5.

Quantifying the Impact of Halogenation on Intermolecular Interactions and Binding Modes of Aromatic Molecules

Halogenation of aromatic molecules is frequently used to modulate intermolecular interactions with ramifications for optoelectronic and mechanical properties. In this work, we accurately quantify and understand the nature of intermolecular interactions in perhalogenated benzene (PHB) clusters. Using benchmark binding energies from the fixed-node diffusion Monte Carlo (FN-DMC) method, we show that generalized Kohn-Sham semicanonical projected random phase approximation (GKS-spRPA) plus approximate exchange kernel (AKX) provides reliable interaction energies with mean absolute error (MAE) of 0.23 kcal/mol. Using the GKS-spRPA+AXK method, we quantify the interaction energies of several binding modes of PHB clusters ((C₆X₆)_n; X = F, Cl, Br, I; n = 2, 3). For a given binding mode, the interaction energies increase 3-4 times from X = F to X = I; the X-X binding modes have energies in the range of 2-4 kcal/mol, while the π-π binding mode has interaction energies in the range of 4-12 kcal/mol. SAPT-DFT-based energy decomposition analysis is then used to show that the equilibrium geometries are dictated primarily by the dispersion and exchange interactions. Finally, we test the accuracy of several dispersion-corrected density functional approximations and show that only the r2SCAN-D4 method has a low MAE and correct long-range behavior, which makes it suitable for large-scale simulations and for developing structure-function relationships of halogenated aromatic systems.

Accuracy of Noncovalent Interactions Involving d-Elements by the 1-Determinant Fixed-Node Diffusion Monte Carlo Method with Effective Core Potentials

A critical assessment of effective core potential (ECP)-based single-determinant (SD) fixed-node diffusion quantum Monte Carlo (FNDMC) accuracy in prototypical noncovalent closed-shell systems involving d-elements is presented. Careful analysis of biases and elimination of possible bias sources leads to two findings of practical importance for SD FNDMC in these systems. First, in some systems (HCu:HCu, HCu:CuH), SD FNDMC reveals large biases of interaction energy differences (significantly exceeding the target 2% relative error) vs a reliable coupled-cluster CCSD(T)/CBS (complete basis set) reference. Second, the leading error of SD FNDMC with ECPs was attributed to a higher nuclear charge Z of d-group (pseudo) atoms, when compared to sp elements, in line with a previously reported finding that aggregate SD FNDMC bias tends to increase in systems with higher electronic densities. Therefore, SD FNDMC should only be used with caution in systems with a large Z.

Benchmarking fundamental gap of Sc2C(OH)2 MXene by many-body methods

Sc₂C(OH)₂ is a prototypical non-magnetic member of MXenes, a promising transition-metal-based 2D material family, with a direct bandgap. We provide here a benchmark of its fundamental gap Δ obtained from many-body GW and fixed-node diffusion Monte Carlo methods. Both approaches independently arrive at a similar value of Δ ∼ 1.3 eV, suggesting the validity of both methods. Such a bandgap makes Sc₂C(OH)₂ a 2D semiconductor suitable for optoelectronic applications. The absorbance spectra and the first exciton binding energy (0.63 eV), based on the Bethe-Salpeter equation, are presented as well. The reported results may serve to delineate experimental uncertainties and enable selection of reasonable approximations such as density functional theory functionals, for use in modeling of related MXenes.

Barrier heights, reaction energies and bond dissociation energies for RH + HO2 reactions with coupled-cluster theory, density functional theory and diffusion quantum Monte Carlo methods

Hydrogen abstraction reactions by the HO₂ radical from hydrocarbon molecules are an important class of reactions in the autoignition of hydrocarbon fuels. Performance of DLPNO-CC and DFT methods using three hybrids and four double hybrids as well as FN-DMC with the single-Slater-Jastrow trial wavefunction on barrier heights and reaction energies of RH + HO₂ reactions as well as bond dissociation energies of the involved X-H molecules is evaluated by comparison with the highly accurate CCSD(T)-F12b/CBS results in this study. Our results show that the DLPNO-CCSD(T)-F12 method can achieve highly accurate barrier heights, reaction energies and X-H bond energies for RH + HO₂ reactions at a relatively low computational cost, and it is applicable to the H-abstraction reactions of larger molecules. Among all DFAs, MN15 and the employed double hybrids can achieve accurate barrier heights and reaction energies with MADs of less than or around 2 kJ mol^-1, but their error on X-H bond energies is more pronounced. Only DSD-BLYP and DSD-PBEB95 can provide X-H bond energies with MADs less than 4 kJ mol^-1. Considering dispersion correction in DFT calculations does not improve these barrier heights and reaction energies. The error of FN-DMC on barrier heights and reaction energies is slightly larger than that of MN15 and those of double hybrids, but it can achieve results within chemical accuracy for these reactions and the X-H bond energies.

Reliability of Computing van der Waals Bond Lengths of Some Rare Gas Diatomics

When the bond lengths of 11 molecules containing van der Waals bonds are optimized by 572 methods and 20 basis sets, it is found that the best mean absolute deviations (MADs) of density-functional theory (DFT) methods are 0.005 Å (shown by APFD/6-311++G**), 0.007 Å (B2PLYPD3(Full)/aug-cc-pVQZ), and 0.010 Å (revDSDPBEP86/aug-cc-pVQZ), while the best MADs of ab initio methods are 0.008 Å (BD(T)/aug-cc-pVTZ) and 0.016 Å (MP4/aug-cc-pVQZ). Moreover, the best MADs calculated by 54 selected methods in combination with 60 other basis sets (such as 6-311++G, 6-31++G(3d'f,3p'd), and UGBS1V++) are not better. Therefore, these bond lengths can be calculated with extremely high accuracy by some special methods and basis sets, and CCSD(T) is also not as good as expected because its best MAD is only 0.023 Å (CCSD(T)/aug-cc-pVQZ).

Probing and Manipulating Noncovalent Interactions in Functional Polymeric Systems

Noncovalent interactions, which usually feature tunable strength, reversibility, and environmental adaptability, have been recognized as driving forces in a variety of biological and chemical processes, contributing to the recognition between molecules, the formation of molecule clusters, and the establishment of complex structures of macromolecules. The marriage of noncovalent interactions and conventional covalent polymers offers the systems novel mechanical, physicochemical, and biological properties, which are highly dependent on the binding mechanisms of the noncovalent interactions that can be illuminated via quantification. This review systematically discusses the nanomechanical characterization of typical noncovalent interactions in polymeric systems, mainly through direct force measurements at microscopic, nanoscopic, and molecular levels, which provide quantitative information (e.g., ranges, strengths, and dynamics) on the binding behaviors. The fundamental understandings of intermolecular and interfacial interactions are then correlated to the macroscopic performances of a series of noncovalently bonded polymers, whose functions (e.g., stimuli-responsiveness, self-healing capacity, universal adhesiveness) can be customized through the manipulation of the noncovalent interactions, providing insights into the rational design of advanced materials with applications in biomedical, energy, environmental, and other engineering fields.

Recognition in the Domain of Molecular Chirality: From Noncovalent Interactions to Separation of Enantiomers

It is not a coincidence that both chirality and noncovalent interactions are ubiquitous in nature and synthetic molecular systems. Noncovalent interactivity between chiral molecules underlies enantioselective recognition as a fundamental phenomenon regulating life and human activities. Thus, noncovalent interactions represent the narrative thread of a fascinating story which goes across several disciplines of medical, chemical, physical, biological, and other natural sciences. This review has been conceived with the awareness that a modern attitude toward molecular chirality and its consequences needs to be founded on multidisciplinary approaches to disclose the molecular basis of essential enantioselective phenomena in the domain of chemical, physical, and life sciences. With the primary aim of discussing this topic in an integrated way, a comprehensive pool of rational and systematic multidisciplinary information is provided, which concerns the fundamentals of chirality, a description of noncovalent interactions, and their implications in enantioselective processes occurring in different contexts. A specific focus is devoted to enantioselection in chromatography and electromigration techniques because of their unique feature as "multistep" processes. A second motivation for writing this review is to make a clear statement about the state of the art, the tools we have at our disposal, and what is still missing to fully understand the mechanisms underlying enantioselective recognition.

Exploring Accurate Potential Energy Surfaces via Integrating Variational Quantum Eigensolver with Machine Learning

The potential energy surface (PES) is crucial for interpreting a variety of chemical reaction processes. However, predicting accurate PESs with high-level electronic structure methods is a challenging task due to the high computational cost. As an appealing application of quantum computing, we show in this work that variational quantum algorithms can be integrated with machine learning (ML) techniques as a promising scheme for exploring accurate PESs. Different from using a ML model to represent the potential energy, we encode the molecular geometry information into a deep neural network (DNN) to represent parameters of the variational quantum eigensolver (VQE), leaving the PES to the wave function ansatz. Once the DNN model is trained, the variational optimization procedure that hinders the application of the VQE to complex systems is avoided, and thus the evaluation of PESs is significantly accelerated. Numerical results demonstrate that a simple DNN model is able to reproduce accurate PESs for small molecules.

Coupled Cluster Benchmarking of Large Noncovalent Complexes in L7 and S12L as Well as the C60 Dimer, DNA-Ellipticine, and HIV-Indinavir

In this work, we report the benchmark binding energies of the seven complexes within the L7 data set, six host-guest complexes from the S12L data set, a C₆₀ dimer, the DNA-ellipticine intercalation complex, and the largest system of the study, the HIV-indinavir system, which contained 343 atoms or 139 heavy atoms. The high-quality values reported were obtained via a focal point method that relies on the canonical form of second-order Møller-Plesset theory and the domain-based local pair natural orbital scheme for the coupled cluster with single double and perturbative triple excitations [DLPNO-CCSD(T)] extrapolated to the complete basis set (CBS) limit. The results in this work not only corroborate but also improve upon some previous benchmark values for large noncovalent complexes albeit at a relatively steep cost. Although local CCSD(T) and the largely successful fixed-node diffusion Monte Carlo (FN-DMC) have been shown to generally agree for small- to medium-size systems, a discrepancy in their reported binding energy values arises for large complexes, where the magnitude of the disagreement is a definite cause for concern. For example, the largest deviation in the L7 data set was 2.8 kcal/mol (∼10%) on the low end in C3GC. Such a deviation only grows worse in the S12L set, which showed a difference of up to 10.4 kcal/mol (∼25%) by a conservative estimation in buckycatcher-C₆₀. The DNA-ellipticine complex also generated a disagreement of 4.4 kcal/mol (∼10%) between both state-of-the-art methods. The disagreement between local CCSD(T) and FN-DMC in large noncovalent complexes shows that it is urgently needed to have the canonical CCSD(T), the Monte Carlo CCSD(T), or the full configuration interaction quantum Monte Carlo approaches available to large systems on the hundred-atom scale to solve this dilemma. In addition, the performances of cheaper popular computational methods were assessed for the studied complexes with respect to DLPNO-CCSD(T)/CBS. r²SCAN-3c, B97M-V, and PBE0+D4 work well in large noncovalent complexes in this work, and GFN2-xTB performs well in π-π stacking complexes. B97M-V is the most reliable computationally efficient approach to predicting noncovalent interactions for large complexes, being the only one to have binding errors within the so-called 1 kcal/mol "chemical accuracy". The benchmark interaction energies of these host-guest complexes, molecular materials, and biological systems with electronic and medicinal implications provide crucial reference data for the improvement of current and future lower-cost methods.

Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface

Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. We discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.

Diffusion quantum Monte Carlo method on diradicals using single- and multi-determinant-Jastrow trial wavefunctions and different orbitals

In this work, the diffusion quantum Monte Carlo (DMC) method is employed to calculate the energies of singlet and triplet states for a series of organic diradicals and diatomic diradicals with π2 configuration. Single-determinant-Jastrow (SDJ) trial wavefunctions for triplet states, two-determinant-Jastrow (2DJ) trial wavefunctions for the singlet states, and multi-determinant-Jastrow (MDJ) trial wavefunctions are employed in DMC calculations using restricted open-shell B3LYP (ROB3LYP) orbitals, complete-active-space self-consistent field (CASSCF) orbitals, state-average CASSCF orbitals, or frozen-CASSCF orbitals. Our results show that DMC energies using either SDJ/2DJ or MDJ with ROB3LYP orbitals are close to or lower than those with the other orbitals for organic diradicals, while they are not very sensitive to the employed orbitals for diatomic diradicals. Furthermore, using MDJ can reduce DMC energies to some extent for most of the investigated organic diradicals and some diatomic diradicals. The importance of MDJ on DMC energies can be estimated based on the percentage of main determinants in the CASCI wavefunction. On the other hand, singlet–triplet gaps can be calculated reasonably with DMC using MDJ with a mean absolute error of less than 2 kcal/mol with all these orbitals. CASCI wavefunctions using density functional theory orbitals are preferred in constructing MDJ trial wavefunctions in practical DMC calculations since it is easier to obtain such wavefunctions than CASSCF methods.

Importance of Van der Waals Interactions in Hydrogen Adsorption on a Silicon-carbide Nanotube Revisited with vdW-DFT and Quantum Monte Carlo

Density functional theory (DFT) is a valuable tool for calculating adsorption energies toward designing materials for hydrogen storage. However, dispersion forces being absent from the local/semi-local theory, it remains unclear as to how the consideration of van der Waals (vdW) interactions affects such calculations. For the first time, we applied diffusion Monte Carlo (DMC) to evaluate the adsorption characteristics of a hydrogen molecule on a (5,5) armchair silicon-carbide nanotube (H₂-SiCNT). Within the DFT framework, we benchmarked various exchange-correlation functionals, including those recently developed for treating dispersion or vdW interactions. We found that the vdW-corrected DFT methods agree well with DMC, whereas the local (semilocal) functional significantly over (under)-binds. Furthermore, we fully optimized the H₂-SiCNT geometry within the DFT framework and investigated the correlation between the structure and charge density. The vdW contribution to the adsorption was found to be non-negligible at ∼1 kcal/mol per hydrogen molecule, which amounts to 9-29% of the ideal adsorption energy required for hydrogen storage applications.

Monte Carlo MP2-F12 for Noncovalent Interactions: The C60 Dimer

A scalable stochastic algorithm is presented that can evaluate explicitly correlated (F12) second-order many-body perturbation (MP2) energies of weak, noncovalent, intermolecular interactions. It first transforms the formulas of the MP2 and F12 energy differences into a short sum of high-dimensional integrals of Green's functions in real space and imaginary time. These integrals are then evaluated by the Monte Carlo method augmented by parallel execution, redundant-walker convergence acceleration, direct-sampling autocorrelation elimination, and control-variate error reduction. By sharing electron-pair walkers across the supermolecule and its subsystems spanned by the joint basis set, the statistical uncertainty is reduced by one to 2 orders of magnitude in the MP2 binding energy corrected for the basis-set incompleteness and superposition errors. The method predicts the MP2-F12/aug-cc-pVDZ binding energy of 19.1 ± 4.0 kcal mol^-1 for the C₆₀ dimer at the center distance of 9.748 Å.

Coupled cluster benchmarks of large noncovalent complexes: The L7 dataset as well as DNA-ellipticine and buckycatcher-fullerene

In this work, benchmark binding energies for dispersion-bound complexes in the L7 dataset, the DNA-ellipticine intercalation complex, and the buckycatcher-C₆₀ complex with 120 heavy atoms using a focal-point method based on the canonical form of second-order Møller-Plesset theory (MP2) and the domain based local pair natural orbital scheme for the coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set (CBS) limit are reported. This work allows for increased confidence given the agreement with respect to values recently obtained using the local natural orbital CCSD(T) for L7 and the canonical CCSD(T)/CBS result for the coronene dimer (C2C2PD). Therefore, these results can be considered pushing the CCSD(T)/CBS binding benchmark to the hundred-atom scale. The disagreements between the two state-of-the-art methods, CCSD(T) and fixed-node diffusion Monte Carlo, are substantial with at least 2.0 (∼10%), 1.9 (∼5%), and 10.3 kcal/mol (∼25%) differences for C2C2PD in L7, DNA-ellipticine, and buckycatcher-C₆₀, respectively. Such sizable discrepancy above "chemical accuracy" for large noncovalent complexes indicates how challenging it is to obtain benchmark binding interactions for systems beyond small molecules, although the three up-to-date density functionals, PBE0+D4, ωB97M-V, and B97M-V, agree better with CCSD(T) for these large systems. In addition to reporting these values, different basis sets and various CBS extrapolation parameters for Hartree-Fock and MP2 correlation energies were tested for the first time in large noncovalent complexes with the goal of providing some indications toward optimal cost effective routes to approach the CBS limit without substantial loss in quality.

Noncovalent Interactions by the Quantum Monte Carlo Method: Strong Influence of Isotropic Jastrow Cutoff Radii

We present a paradigmatic example of a strong effect of Jastrow cutoff radii setup on the accuracy of noncovalent interaction energy differences within one-determinant Slater-Jastrow fixed-node diffusion Monte Carlo (1FNDMC) simulations using isotropic Jastrow terms and effective-core potentials. Analysis of total energies, absolute and relative errors, and local energy variance of energy differences vs the reference results suggests a simple procedure to marginalize the related biases. The presented data showcase improvements in dispersion-bounded systems within such a 1FNDMC method.

Interactions between large molecules pose a puzzle for reference quantum mechanical methods

Quantum-mechanical methods are used for understanding molecular interactions throughout the natural sciences. Quantum diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] are state-of-the-art trusted wavefunction methods that have been 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. We show that CCSD(T) and DMC interaction energies are not consistent for a set of polarizable supramolecules. Whilst there is agreement for some of the complexes, in a few key systems disagreements of up to 8 kcal mol^-1 remain. These findings thus indicate that more caution is required when aiming at reproducible non-covalent interactions between extended molecules.

Toward quantum Monte Carlo forces on heavier ions: Scaling properties

Quantum Monte Carlo (QMC) forces have been studied extensively in recent decades because of their importance with spectroscopic observables and geometry optimization. Here, we benchmark the accuracy and computational cost of QMC forces. The zero-variance zero-bias (ZVZB) force estimator is used in standard variational and diffusion Monte Carlo simulations with mean-field based trial wavefunctions and atomic pseudopotentials. Statistical force uncertainties are obtained with a recently developed regression technique for heavy tailed QMC data [P. Lopez Rios and G. J. Conduit, Phys. Rev. E 99, 063312 (2019)]. By considering selected atoms and dimers with elements ranging from H to Zn (1 ≤ Z_eff ≤ 20), we assess the accuracy and the computational cost of ZVZB forces as the effective pseudopotential valence charge, Z_eff, increases. We find that the costs of QMC energies and forces approximately follow simple power laws in Z_eff. The force uncertainty grows more rapidly, leading to a best case cost scaling relationship of approximately Z_eff ^6.5(3) for diffusion Monte Carlo. We find that the accessible system size at fixed computational cost scales as Z_eff ^-2, insensitive to model assumptions or the use of the "space warp" variance-reduction technique. Our results predict the practical cost of obtaining forces for a range of materials, such as transition metal oxides where QMC forces have yet to be applied, and underscore the importance of further developing force variance-reduction techniques, particularly for atoms with high Z_eff.

Noncovalent Interactions from Models for the Møller-Plesset Adiabatic Connection

Given the omnipresence of noncovalent interactions (NCIs), their accurate simulations are of crucial importance across various scientific disciplines. Here we construct accurate models for the description of NCIs by an interpolation along the Møller-Plesset adiabatic connection (MP AC). Our interpolation approximates the correlation energy, by recovering MP2 at small coupling strengths and the correct large-coupling strength expansion of the MP AC, recently shown to be a functional of the Hartree-Fock density. Our models are size consistent for fragments with nondegenerate ground states, have the same cost as double hybrids, and require no dispersion corrections to capture NCIs accurately. These interpolations greatly reduce large MP2 errors for typical π-stacking complexes (e.g., benzene-pyridine dimers) and for the L7 data set. They are also competitive with state-of-the-art dispersion enhanced functionals and can even significantly outperform them for a variety of data sets, such as CT7 and L7.

Fragmentation of natural orbital occupation numbers-based diagnostic of differential multireference character in complexes with hydrogen bonds

We explore the possible route to approximate natural orbital occupation numbers-based diagnostic of differential multireference character of noncovalent energy differences by techniques based on many-body expansion. It turns out that two-body fragmentation of monomers may lead to a reasonable approximation of such a diagnostic in hydrogen-bonded complexes. The results are useful, for example, for assessment of the expected bias cancellation in energy differences of larger systems obtained by single-reference methods.

Accurate Reduced-Cost CCSD(T) Energies: Parallel Implementation, Benchmarks, and Large-Scale Applications

The accurate and systematically improvable frozen natural orbital (FNO) and natural auxiliary function (NAF) cost-reducing approaches are combined with our recent coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] implementations. Both of the closed- and open-shell FNO-CCSD(T) codes benefit from OpenMP parallelism, completely or partially integral-direct density-fitting algorithms, checkpointing, and hand-optimized, memory- and operation count effective implementations exploiting all permutational symmetries. The closed-shell CCSD(T) code requires negligible disk I/O and network bandwidth, is MPI/OpenMP parallel, and exhibits outstanding peak performance utilization of 50-70% up to hundreds of cores. Conservative FNO and NAF truncation thresholds benchmarked for challenging reaction, atomization, and ionization energies of both closed- and open-shell species are shown to maintain 1 kJ/mol accuracy against canonical CCSD(T) for systems of 31-43 atoms even with large basis sets. The cost reduction of up to an order of magnitude achieved extends the reach of FNO-CCSD(T) to systems of 50-75 atoms (up to 2124 atomic orbitals) with triple- and quadruple-ζ basis sets, which is unprecedented without local approximations. Consequently, a considerably larger portion of the chemical compound space can now be covered by the practically "gold standard" quality FNO-CCSD(T) method using affordable resources and about a week of wall time. Large-scale applications are presented for organocatalytic and transition-metal reactions as well as noncovalent interactions. Possible applications for benchmarking local CCSD(T) methods, as well as for the accuracy assessment or parametrization of less complete models, for example, density functional approximations or machine learning potentials, are also outlined.

Diffusion Monte Carlo method on small boron clusters using single- and multi- determinant-Jastrow trial wavefunctions

Multireference character in some small boron clusters could be significant, and a previous all-electron fixed-node diffusion quantum Monte Carlo (FN-DMC) calculation with the single-determinant-Jastrow (SDJ) trial wavefunction shows that the atomization energy (AE) of B₄ ⁺ is overestimated by about 1.4 eV compared with the coupled cluster method with single, doubles, and perturbative triples [CCSD(T)] results. All-electron FN-DMC calculations and those with the pseudopotential (PP) using SDJ and multi-determinant-Jastrow (MDJ) trial wavefunctions with B3LYP orbitals as well as CC calculations at different levels are carried out on B_n ^Q (n = 1-5, Q = -1, 0, 1) clusters. The obtained FN-DMC energies indicate that the node error of the employed SDJ trial wavefunction in all-electron calculations is different from that with the PP for some clusters. The error of AEs and dissociation energies (DEs) from all-electron FN-DMC calculations is larger than that with the PP when the SDJ trial wavefunction is employed, while errors of CC methods do not depend on whether the PP is used. AEs and DEs of the boron clusters are improved significantly when MDJ trial wavefunctions are used in both all-electron calculations and those with the PP, and their error is similar to that of CCSD(T) compared with CCSDT(Q) results. On the other hand, reasonable adiabatic electron detachment energies (ADEs) and ionization potentials (AIPs) are achieved with FN-DMC using SDJ trial wavefunctions and MDJ is less effective on ADEs and AIPs. Furthermore, the relative energy between two structures of B₉ ^- is predicted reliably with FN-DMC using the SDJ trial wavefunction and the effect of MDJ is negligible, while density functional theory results using different exchange-correlation functionals differ significantly.

On the role of steric and exchange–correlation effects in halogenated complexes

Density functional theory formalisms of energy partitioning schemes are utilized to find out what energetic components govern interactions in halogenated complexes.

Perspectives on van der Waals Density Functionals: The Case of TiS2

The van der Waals interaction is of foundational importance for a wide variety of physical systems. In particular, van der Waals forces lie at the heart of potential device technologies that may be realized from the functional organization of layered two-dimensional (2D) nanomaterials. For intermediate to large-scale applications modeling, van der Waals density functionals have become the de facto choice for first-principles calculations. In particular, the vdW-DF family of functionals have provided a systematic approach to this theoretically challenging problem. While much progress has been made, there remains room for improvement in the microscopic description of vdW forces from these density functionals. In this work, we compute benchmark results for the binding energy and the electronic density response to binding in TiS₂ via accurate diffusion quantum Monte Carlo calculations. We compare these benchmark data to results obtained from local, semilocal, and van der Waals functionals. In particular, we gauge the quality of the original vdW-DF/vdW-DF2 functionals, as well as updated variants such as vdW-DF-C09, vdW-DF-optB88, vdW-DF-optB86b, and vdW-DF2-B86R. We find a close relationship between the accuracy of predicted interlayer separation distances and binding energies for TiS₂, with the vdW-DF-optB88 functional performing very well in terms of both quantities. In general, the more recently developed functionals are systematic improvements over older ones. However, when considering the response of the electron density to binding, we find that local-density approximation (LDA) and PBEsol generally outperform the vdW-DF functionals in describing the interlayer charge accumulation with vdW-DF-C09 variants performing the best overall.

MAP: An MP2 Accuracy Predictor for Weak Interactions from Adiabatic Connection Theory

Second-order Møller-Plesset perturbation theory (MP2) approximates the exact Hartree-Fock (HF) adiabatic connection (AC) curve by a straight line. Thus, by using the deviation of the exact curve from the linear behavior, we construct an indicator for the accuracy of MP2. We then use an interpolation along the HF AC to transform the exact form of our indicator into a highly practical MP2 accuracy predictor (MAP) that comes at a negligible additional computational cost. We show that this indicator is already applicable to systems that dissociate into fragments with a nondegenerate ground state, and we illustrate its usefulness by applying it to the S22 and S66 datasets.

On the Possibility of Helium Adsorption in Nitrogen Doped Graphitic Materials

The potassium salt of polyheptazine imide (K-PHI) is a promising photocatalyst for various chemical reactions. From powder X-ray diffraction data an idealized structural model of K-PHI has been derived. Using atomic coordinates of this model we defined an energetically optimized K-PHI structure, in which the K ions are present in the pore and between the PHI-planes. The distance between the anion framework and K+ resembles a frustrated Lewis pair-like structure, which we denote as frustrated Coulomb pair that results in an interesting adsorption environment for otherwise non-adsorbing, non-polar gas molecules. We demonstrate that even helium (He) gas molecules, which are known to have the lowest boiling point and the lowest intermolecular interactions, can be adsorbed in this polarized environment with an adsorption energy of  - 4.6 kJ mol-1 per He atom. The interaction between He atoms and K-PHI is partially originating from charge transfer, as disclosed by our energy decomposition analysis based on absolutely localized molecular orbitals. Due to very small charge transfer interactions, He gas adsorption saturates at 8 at%, which however can be subject to further improvement by cation variation.

Modified Density Functional Dispersion Correction for Inorganic Layered MFX Compounds (M = Ca, Sr, Ba, Pb and X = Cl, Br, I)

MFX (M = Ca, Ba, Sr, Pb and X = Cl, Br, I) compounds have received considerable attention due to their technological application as X-ray detectors, pressure sensors, and optical data storage materials, when doped with rare-earth ions. MFX compounds belong to the class of layered materials with a tetragonal Matlockite crystal structure, characterized by weakly stacked double-halide layers along the crystallographic c-axis. These layers predominantly determine phase transitions, elastic, and mechanical properties. However, the correct description of the lattice parameter c is a challenge for most standard DFT functionals, which tend to overestimate the lattice parameter c. Because of the weak interactions between the halide layers, dispersion-corrected functionals seem to be a better choice. We investigated 11 different inorganic layered MFX compounds for which experimental data are available, with standard and dispersion-corrected functionals to assess their performance in reproducing the lattice parameter c, structural, and vibrational properties of the MFX compounds. Our results revealed that these functionals do not describe the weak interactions between the halide layers in a balanced way. Therefore, we modified Grimme's popular DFT-D2 dispersion correction scheme in two different ways by (i) replacing the dispersion coefficients and van der Waals radii with those of noble gas atoms or (ii) increasing the van der Waals radii of the MFX atoms up to 40%. Comparison with the available experimental data revealed that the latter approach applied to the PBE (Perdew-Burke-Ernzerhof)-D2 functional with 30% increased van der Waals radii, which we coined PBE-D2* (S_rvdw 1.30) is best suited to fine-tune the description of the weak interlayer interactions in MFX compounds, thus significantly improving the description of their structural, vibrational, and mechanical properties. Work is in progress applying this new, computationally inexpensive scheme to other inorganic layered compounds and periodic systems with weakly stacked layers.

Fractional Charge by Fixed-Node Diffusion Monte Carlo Method

Fixed-node diffusion Monte Carlo (FNDMC) method is a stochastic quantum many-body approach that has a great potential in electronic structure theory. We examine how FNDMC total energy E(N) satisfies exact constraints, linearity and derivative discontinuity, versus fractional electron number N, if combined with mean-field trial wave functions that miss such features. H and Cl atoms with fractional charge reveal that FNDMC method is well able to restore the piecewise linearity of E(N). The method uses ensemble and projector ingredients to achieve the correct charge localization. A water-solvated Cl^{-} complex illustrates superior performance of FNDMC method for charged noncovalent systems.

Singlet-triplet gaps in diradicals obtained with diffusion quantum Monte Carlo using a Slater-Jastrow trial wavefunction with a minimum number of determinants

Diradicals are essential species in a wide range of chemical processes, whereas the computational study of their electronic structure often remains a challenge due to near-degeneracy of the frontier molecular orbitals. The fixed-node diffusion quantum Monte Carlo (FN-DMC) method is employed to calculate adiabatic energy gaps of some typical diradicals with the Slater-Jastrow trial wavefunction. The antisymmetrized part of the trial wavefunction is taken to be a linear combination of a minimum number of determinants using RB3LYP orbitals from the closed-shell singlet state or ROB3LYP orbitals from the triplet state. Our results show that using the two-determinant-Jastrow trial wavefunction is necessary to achieve reliable energy differences between closed-shell singlet states. The energy of the triplet state with MS = 1 is calculated to be lower than that with MS = 0 with FN-DMC even using trial wavefunctions with spin-pure states as their antisymmetrized parts and this difference is reduced with better orbitals. This indicates that the fixed-node error is smaller for the triplet state with MS = 1. Adiabatic energy gaps obtained from the present FN-DMC calculations are in reasonable agreement with available experimental values. Compared with results of the high level EOM-SF-CC method, energy gaps of FN-DMC with RB3LYP orbitals are slightly better than those using ROB3LYP orbitals and results of EOM-SF-CCSD. The present FN-DMC calculations using the simplest ansatz for the trial wavefunction can achieve reasonable results for these diradicals and they can readily be applied to large diradicals.