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 ).

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.

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.

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.

Efficient and accurate description of adsorption in zeolites

Accurate theoretical methods are needed to correctly describe adsorption on solid surfaces or in porous materials. The random phase approximation (RPA) with singles corrections scheme and the second order Møller-Plesset perturbation theory (MP2) are two schemes, which offer high accuracy at affordable computational cost. However, there is little knowledge about their applicability and reliability for different adsorbates and surfaces. Here, we calculate adsorption energies of seven different molecules in zeolite chabazite to show that RPA with singles corrections is superior to MP2, not only in terms of accuracy but also in terms of computer time. Therefore, RPA with singles is a suitable scheme for obtaining highly accurate adsorption energies in porous materials and similar 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.

Quantum Monte Carlo Study of the Water Dimer Binding Energy and Halogen-π Interactions

Halogen-π systems are involved with competition between halogen bonding and π-interaction. Using the diffusion quantum Monte Carlo (DMC) method, we compare the interaction of benzene with diatomic halogens (X₂: Cl₂/Br₂) with the typical hydrogen bonding in the water dimer, taking into account explicit correlations of up to three bodies. The benzene-Cl₂/Br₂ binding energies (13.07 ± 0.42/16.62 ± 0.02 kJ/mol) attributed to both halogen bonding and dispersion are smaller than but comparable to the typical hydrogen bonding in the water dimer binding energy (20.88 ± 0.27 kJ/mol). All of the above values are in good agreement with those from the coupled-cluster with single, double, and noniterative triple excitations (CCSD(T)) results at the complete basis set limit (benzene-Cl₂/Br₂: 12.78/16.17 kJ/mol; water dimer: 21.0 kJ/mol).

Ion-Hydroxyl Interactions: From High-Level Quantum Benchmarks to Transferable Polarizable Force Fields

Ion descriptors in molecular mechanics models are calibrated against reference data on ion-water interactions. It is then typically assumed that these descriptors will also satisfactorily describe interactions of ions with other functional groups, such as those present in biomolecules. However, several studies now demonstrate that this transferability assumption produces, in many different cases, large errors. Here we address this issue in a representative polarizable model and focus on transferability of cationic interactions from water to a series of alcohols. Both water and alcohols use hydroxyls for ion-coordination, and, therefore, this set of molecules constitutes the simplest possible case of transferability. We obtain gas phase reference data systematically from "gold-standard" quantum Monte Carlo and CCSD(T) methods, followed by benchmarked vdW-corrected DFT. We learn that the original polarizable model yields large gas phase water → alcohol transferability errors - the RMS and maximum errors are 2.3 and 5.1 kcal/mol, respectively. These errors are, nevertheless, systematic in that ion-alcohol interactions are overstabilized, and systematic errors typically imply that some essential physics is either missing or misrepresented. A comprehensive analysis shows that when both low- and high-field responses of ligand dipole polarization are described accurately, then transferability improves significantly - the RMS and maximum errors in the gas phase reduce, respectively, to 0.9 and 2.5 kcal/mol. Additionally, predictions of condensed phase transfer free energies also improve. Nevertheless, within the limits of the extrathermodynamic assumptions necessary to separate experimental estimates of salt dissolution into constituent cationic and anionic contributions, we note that the error in the condensed phase is systematic, which we attribute, at least, partially to the parametrization in long-range electrostatics. Overall, this work demonstrates a rational approach to boosting transferability of ionic interactions that will be applicable broadly to improving other polarizable and nonpolarizable models.

The S66 Non-Covalent Interactions Benchmark Reconsidered Using Explicitly Correlated Methods Near the Basis Set Limit

The S66 benchmark for non-covalent interactions has been re-evaluated using explicitly correlated methods with basis sets near the one-particle basis set limit. It is found that post-MP2 ‘high-level corrections’ are treated adequately well using a combination of CCSD(F12*) with (aug-)cc-pVTZ-F12 basis sets on the one hand, and (T) extrapolated from conventional CCSD(T)/heavy-aug-cc-pV{D,T}Z on the other hand. Implications for earlier benchmarks on the larger S66×8 problem set in particular, and for accurate calculations on non-covalent interactions in general, are discussed. At a slight cost in accuracy, (T) can be considerably accelerated by using sano-V{D,T}Z+ basis sets, whereas half-counterpoise CCSD(F12*)(T)/cc-pVDZ-F12 offers the best compromise between accuracy and computational cost.

A Computational Scheme To Evaluate Hamaker Constants of Molecules with Practical Size and Anisotropy

We propose a computational scheme to evaluate Hamaker constants, A, of molecules with practical sizes and anisotropies. Upon the increasing feasibility of diffusion Monte Carlo (DMC) methods to evaluate binding curves for such molecules to extract the constants, we discussed how to treat the averaging over anisotropy and how to correct the bias due to the nonadditivity. We have developed a computational procedure for dealing with the anisotropy and reducing statistical errors and biases in DMC evaluations, based on possible validations on predicted A. We applied the scheme to cyclohexasilane molecule, Si₆H₁₂, used in "printed electronics" fabrications, getting A ≈ 105 ± 2 zJ, being in plausible range supported even by other possible extrapolations. The scheme provided here would open a way to use handy ab initio evaluations to predict wettabilities as in the form of materials informatics over broader molecules.

Bias cancellation in one-determinant fixed-node diffusion Monte Carlo: Insights from fermionic occupation numbers

The accuracy of the fixed-node diffusion Monte Carlo (FNDMC) depends on the node location of the supplied trial state Ψ_{T}. The practical FNDMC approaches available for large systems rely on compact yet effective Ψ_{T}, most often containing an explicitly correlated single Slater determinant (SD). However, SD nodes may be better suited to one system than to another, which may possibly lead to inaccurate FNDMC energy differences. It remains a challenge how to estimate nonequivalence or appropriateness of SDs. Here we use the differences of a measure based on the Euclidean distance between the natural orbital occupation number (NOON) vector of the SD and the exact solution in the NOON vector space, which can be viewed as a measure of SD nonequivalence and as a qualitative measure of the expected degree of nondynamic-correlation-related bias in FNDMC energy differences. This is explored on a set of small noncovalent complexes and covalent bond breaking of Si_{2} vs N_{2}. It turns out that NOON-based measures well reflect the magnitude and sign of the bias present in the data available, thus providing insights into the nature of bias cancellation in SD FNDMC energy differences.

Hexagonal boron nitride and water interaction parameters

The study of hexagonal boron nitride (hBN) in microfluidic and nanofluidic applications at the atomic level requires accurate force field parameters to describe the water-hBN interaction. In this work, we begin with benchmark quality first principles quantum Monte Carlo calculations on the interaction energy between water and hBN, which are used to validate random phase approximation (RPA) calculations. We then proceed with RPA to derive force field parameters, which are used to simulate water contact angle on bulk hBN, attaining a value within the experimental uncertainties. This paper demonstrates that end-to-end multiscale modeling, starting at detailed many-body quantum mechanics and ending with macroscopic properties, with the approximations controlled along the way, is feasible for these systems.

Diffusion Monte Carlo Perspective on the Spin-State Energetics of [Fe(NCH)6](2.)

The energy difference between the high spin and the low spin state of the model compound [Fe(NCH)6](2+) is investigated by means of Diffusion Monte Carlo (DMC), where special attention is dedicated to analyzing the effect of the fix node approximation on the accuracy of the results. For this purpose, we compare several Slater-Jastrow and multireference Slater-Jastrow trial wave functions. We found that a Slater-Jastrow trial wave function constructed with the generalized Kohn-Sham orbitals from hybrid DFT represents the optimal choice. This is understood by observing that hybrid functionals account for the subtle balance between exchange and correlation effects and the respective orbitals accurately describe the ligand-metal hybridization as well as the charge reorganization accompanying the spin transition. Finally the DMC results are compared with those obtained by Hartree-Fock, DFT, CASSCF, and CASPT2. While there is no clear reference value for the high spin-low spin energy difference, DMC and high level CCSD(T) calculations agree within around 0.3 eV.

Stabilizing and Modulating Color by Copigmentation: Insights from Theory and Experiment

Natural anthocyanin pigments/dyes and phenolic copigments/co-dyes form noncovalent complexes, which stabilize and modulate (in particular blue, violet, and red) colors in flowers, berries, and food products derived from them (including wines, jams, purees, and syrups). This noncovalent association and their electronic and optical implications constitute the copigmentation phenomenon. Over the past decade, experimental and theoretical studies have enabled a molecular understanding of copigmentation. This review revisits this phenomenon to provide a comprehensive description of the nature of binding (the dispersion and electrostatic components of π-π stacking, the hydrophobic effect, and possible hydrogen-bonding between pigment and copigment) and of spectral modifications occurring in copigmentation complexes, in which charge transfer plays an important role. Particular attention is paid to applications of copigmentation in food chemistry.

The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus

Sensitive dependence of the electronic structure on the number of layers in few-layer phosphorene raises a question about the true nature of the interlayer interaction in so-called "van der Waals (vdW) solids". We performed quantum Monte Carlo calculations and found that the interlayer interaction in bulk black phosphorus and related few-layer phosphorene is associated with a significant charge redistribution that is incompatible with purely dispersive forces and not captured by density functional theory calculations with different vdW corrected functionals. These findings confirm the necessity of more sophisticated treatment of nonlocal electron correlation in total energy calculations.

Extensions and applications of the A24 data set of accurate interaction energies

The A24 data set (Řezáč and Hobza, J. Chem. Theory Comput. 2013, 9, 2151-2155) is a set of noncovalent complexes large enough to showcase various types of interactions yet small enough to make highly accurate calculations possible. It is intended for the testing of accurate computational methods which are then used as a benchmark in larger model systems. In this work, we improve the best estimate of the interaction energies in the set by updating their CCSD(T)/CBS and CCSDT(Q) components with calculations in larger basis sets. The data set is then used to test a large number of composite CCSD(T) approaches. Special attention is paid to the use of the explicitly correlated MP2-F12 method in these composite calculations. It is shown that an accuracy of 1-2% can be achieved with setups applicable to larger molecules. The effect of frozen natural orbital approximation on the accuracy of composite CCSD(T)/CBS calculations is also quantified. In four trivially saturated complexes where CCSDT(Q)/CBS data are now available, the convergence of the many-body correlation effects is assessed by fixed-node diffusion Monte Carlo (FN-DMC) calculations. A good agreement is achieved between FN-DMC and high-level coupled-cluster which represents an important cross-check of both approaches.

Choosing a density functional for modeling adsorptive hydrogen storage: reference quantum mechanical calculations and a comparison of dispersion-corrected density functionals

Various dispersion-corrected density functionals are compared with high level QM data for several model complexes for adsorptive hydrogen storage.