Fixed-node Diffusion Monte Carlo Method Research Articles

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

Sc2C(OH)2 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.3eV, suggesting the validity of both methods. Such a bandgap makes Sc2C(OH)2 a 2D semiconductor suitable for optoelectronic applications. The absorbance spectra and the first exciton binding energy (0.63eV), 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.

CrI3 revisited with a many-body ab initio theoretical approach

${\mathrm{CrI}}_{3}$ has recently been shown to exhibit low-dimensional, long-range magnetic ordering from few layers to single layers of ${\mathrm{CrI}}_{3}$. The properties of ${\mathrm{CrI}}_{3}$ bulk and few-layered systems are uniquely defined by a combination of short-range intralayer and long-range interlayer interactions, including strong correlations, exchange, and spin-orbit coupling. Unfortunately, both the long-range van der Waals interactions, which are driven by dynamic, many-body electronic correlations, and the competing strong intralayer correlations, present a formidable challenge for the local or semilocal mean-field approximations employed in workhorse electronic structure approaches like density-functional theory. In this paper we employ a sophisticated many-body approach that can simultaneously describe long- and short-range correlations. We establish that the fixed-node diffusion Monte Carlo (FNDMC) method reproduces the experimental interlayer separation distance of bulk ${\mathrm{CrI}}_{3}$ for the high-temperature monoclinic phase with a reliable prediction of the interlayer binding energy. We subsequently employed the FNDMC results to benchmark the accuracy of several density-functional theory exchange-correlation approximations.

Fundamental gap of fluorographene by many-body GW and fixed-node diffusion Monte Carlo methods.

Fluorographene (FG) is a promising graphene-derived material with a large bandgap. Currently existing predictions of its fundamental gap (Δf) and optical gap (Δopt) significantly vary when compared with experiment. We provide here an ultimate benchmark of Δf for FG by many-body GW and fixed-node diffusion Monte Carlo (FNDMC) methods. Both approaches independently arrive at Δf ≈ 7.1 ± 0.1 eV. In addition, the Bethe-Salpeter equation enabled us to determine the first exciton binding energy, Eb = 1.92 eV. We also point to the possible misinterpretation problem of the results obtained for gaps of solids by FNDMC with single-reference trial wave functions of Bloch orbitals. We argue why instead of Δopt, in the thermodynamic limit, such an approach results in energy differences that rather correspond to Δf, and we also outline conditions when this case actually applies.

Fractional Charge by Fixed-Node Diffusion MonteCarlo Method.

Fixed-node diffusion MonteCarlo (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.

Quantum Monte Carlo simulation for the many-body decomposition of the interaction energy and electron correlation of small superalkali lithium clusters.

Using the fixed-node diffusion Monte Carlo (FN-DMC) method, we calculate the total energy of small lithium clusters Lin (n = 2-6) to obtain the many-body decomposition of the interaction energy of 2- up to 6-body interactions. The obtained many-body decomposition of the interaction energy shows an alternating series with even and odd terms of attractive and repulsive contributions, respectively. The two-body attractive interactions guarantee the stability of the Li2, Li3, and Li4 clusters. For larger clusters Li5 and Li6, the 4-body attractive interactions are required for their stabilization once the strength of the 3-body repulsive interactions overwhelms that of the 2-body attractive ones. With increasing the cluster size, the additive and nonadditive contributions to the interaction energy increase linearly in magnitude but with different slopes for the two-dimensional (2D) planar and three-dimensional (3D) cagelike clusters. The significant increment in nonadditive effects from the 4-atom to the 5-atom cluster has driven the structural transition from 2D to 3D. Combining the FN-DMC calculations with the Hartree-Fock many-body decomposition of the interaction energy, we extract the correlation effects, showing that an odd-even competition pattern in the many-body repulsive and attractive interactions is crucial for the stabilization of the clusters.

Density functional orbitals in quantum Monte Carlo: The importance of accurate densities.

There has been significant recent attention surrounding the accuracy of electronic densities produced by modern parameterized density functional approximations (DFAs). Here, we investigate the impact of using orbitals from density functional calculations in fixed-node Diffusion Monte Carlo (DMC) methods, which is common practice in the calculation of large systems. We find that the accuracy of the density is a strong indicator of the quality of the many-body nodal surface produced by a determinant of the corresponding Kohn-Sham orbitals. Functionals which produce the most accurate electronic densities also produce the lowest variational DMC energies, while functionals that produce poor densities lead to significantly higher energies. This result simplifies the process of choosing orbitals for DMC calculations of large systems and suggests that prioritizing accurate densities in the future development of DFAs would also contribute to the continued improvement of DMC.

Influence of pseudopotentials on excitation energies from selected configuration interaction and diffusion Monte Carlo

Due to their diverse nature, the faithful description of excited states within electronic structure theory methods remains one of the grand challenges of modern theoretical chemistry. Quantum Monte Carlo (QMC) methods have been applied very successfully to ground state properties but still remain generally less effective than other non-stochastic methods for electronically excited states. Nonetheless, we have recently reported accurate excitation energies for small organic molecules at the fixed-node diffusion Monte Carlo (FN-DMC) within a Jastrow-free QMC protocol relying on a deterministic and systematic construction of nodal surfaces using the selected configuration interaction (sCI) algorithm known as CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively). Albeit highly accurate, these all-electron calculations are computationally expensive due to the presence of core electrons. One very popular approach to remove these chemically-inert electrons from the QMC simulation is to introduce pseudopotentials (also known as effective core potentials). Taking the water molecule as an example, we investigate the influence of Burkatzki-Filippi-Dolg (BFD) pseudopotentials and their associated basis sets on vertical excitation energies obtained with sCI and FN-DMC methods. Although these pseudopotentials are known to be relatively safe for ground state properties, we evidence that special care may be required if one strives for highly accurate vertical transition energies. Indeed, comparing all-electron and valence-only calculations, we show that using pseudopotentials with the associated basis sets can induce differences of the order of 0.05 eV on the excitation energies. Fortunately, a reasonable estimate of this shift can be estimated at the sCI level.

Composite Boson Description of a Low-Density Gas of Excitons

Ground state properties of a fermionic Coulomb gas are calculated using the fixed-node diffusion Monte Carlo method. The validity of the composite boson description is tested for different densities. We extract the exciton-exciton $s$-wave scattering length by solving the four-body problem in a harmonic trap and mapping the energy to that of two trapped bosons. The equation of state is consistent with the Bogoliubov theory for composite bosons interacting with the obtained $s$-wave scattering length. The perturbative expansion at low density has contributions physically coming from (a) exciton binding energy, (b) mean-field Gross-Pitaevskii interaction between excitons, (c) quantum depletion of the excitonic condensate (Lee-Huang-Yang terms for composite bosons). In addition, for low densities we find a good agreement with the Bogoliubov bosonic theory for the condensate fraction of excitons. The equation of state in the opposite limit of large density is found to be well described by the perturbative theory including (a) mixture of two ideal Fermi gases (b) exchange energy. We find that for low densities both energetic and coherent properties are correctly described by the picture of composite bosons (excitons).

Fixed-node diffusion Monte Carlo description of nitrogen defects in zinc oxide

Using the fixed-node diffusion Monte Carlo (FN-DMC) method, we evaluate the formation energies and charge transition levels of substitutional nitrogen defects in the wide-band-gap semiconductor zinc oxide. The use of a direct-solution, many-body approach inherently secures a good description of electron-electron interactions, achieving high accuracy without adjustable parameters. According to FN-DMC nitrogen is a deep acceptor with a charge transition level 1.0(3) eV above the valence-band maximum when 72-atom supercells are used. This result falls on the lower end of typically reported hybrid density functional results for the same size supercells, which range from 1.0 to 1.8 eV. Further, residual finite-size effects due to charged defect image interactions in the 72-atom supercells are estimated by supercell extrapolation within hybrid density functional theory. When the finite-size correction is included, we obtain a deep acceptor at 1.6(3) eV. This result is in good agreement with recent experimental measurements. We also analyze the local compressibility of charge according to FN-DMC and common density functionals and find that the use of hybrid functionals obtains compressibilities in better agreement with the many-body theory. Our work illustrates the application of the FN-DMC method to a challenging point defect problem, demonstrating that uncertainties and approximations can be well controlled.

Composite-boson approach to molecular Bose-Einstein condensates in mixtures of ultracold Fermi gases

We show that an ansatz based on independent composite bosons [Phys. Rep. 463, 215 (2008)] accurately describes the condensate fraction of molecular Bose-Einstein condensates in ultracold Fermi gases. The entanglement between the fermionic constituents of a single Feshbach molecule then governs the many-particle statistics of the condensate, from the limit of strong interaction to close to unitarity. This result strengthens the role of entanglement as the indispensable driver of composite-boson behavior. The condensate fraction of fermion pairs at zero temperature that we compute matches excellently previous results obtained by means of fixed-node diffusion Monte Carlo methods and the Bogoliubov depletion approximation. This paves the way towards the exploration of the BEC-BCS crossover physics in mixtures of cold Fermi gases with an arbitrary number of fermion pairs as well as the implementation of Hong-Ou-Mandel-like interference experiments proposed within coboson theory.

One-dimensional multicomponent Fermi gas in a trap: quantum Monte Carlo study

A one-dimensional world is very unusual as there is an interplay between quantum statistics and geometry, and a strong short-range repulsion between atoms mimics Fermi exclusion principle, fermionizing the system. Instead, a system with a large number of components with a single atom in each, on the opposite acquires many bosonic properties. We study the ground-state properties of a multicomponent repulsive Fermi gas trapped in a harmonic trap by a fixed-node diffusion Monte Carlo method. The interaction between all components is considered to be the same. We investigate how the energetic properties (energy, contact) and correlation functions (density profile and momentum distribution) evolve as the number of components is changed. It is shown that the system fermionizes in the limit of strong interactions. Analytical expressions are derived in the limit of weak interactions within the local density approximation for an arbitrary number of components and for one plus one particle using an exact solution.

How large are nonadiabatic effects in atomic and diatomic systems?

With recent developments in simulating nonadiabatic systems to high accuracy, it has become possible to determine how much energy is attributed to nuclear quantum effects beyond zero-point energy. In this work, we calculate the non-relativistic ground-state energies of atomic and molecular systems without the Born-Oppenheimer approximation. For this purpose, we utilize the fixed-node diffusion Monte Carlo method, in which the nodes depend on both the electronic and ionic positions. We report ground-state energies for all systems studied, ionization energies for the first-row atoms and atomization energies for the first-row hydrides. We find the ionization energies of the atoms to be nearly independent of the Born-Oppenheimer approximation, within the accuracy of our results. The atomization energies of molecular systems, however, show small effects of the nonadiabatic coupling between electrons and nuclei.

Practical Schemes for Accurate Forces in Quantum Monte Carlo.

While the computation of interatomic forces has become a well-established practice within variational Monte Carlo (VMC), the use of the more accurate Fixed-Node Diffusion Monte Carlo (DMC) method is still largely limited to the computation of total energies on structures obtained at a lower level of theory. Algorithms to compute exact DMC forces have been proposed in the past, and one such scheme is also put forward in this work, but remain rather impractical due to their high computational cost. As a practical route to DMC forces, we therefore revisit here an approximate method, originally developed in the context of correlated sampling and named here the Variational Drift-Diffusion (VD) approach. We thoroughly investigate its accuracy by checking the consistency between the approximate VD force and the derivative of the DMC potential energy surface for the SiH and C2 molecules and employ a wide range of wave functions optimized in VMC to assess its robustness against the choice of trial function. We find that, for all but the poorest wave function, the discrepancy between force and energy is very small over all interatomic distances, affecting the equilibrium bond length obtained with the VD forces by less than 0.004 au. Furthermore, when the VMC forces are approximate due to the use of a partially optimized wave function, the DMC forces have smaller errors and always lead to an equilibrium distance in better agreement with the experimental value. We also show that the cost of computing the VD forces is only slightly larger than the cost of calculating the DMC energy. Therefore, the VD approximation represents a robust and efficient approach to compute accurate DMC forces, superior to the VMC counterparts.

Dynamic structure function of a cold Fermi gas at unitarity

We present a theoretical study of the dynamic structure function of a resonantly interacting two-component Fermi gas at zero temperature. Our approach is based on dynamic many-body theory able to describe excitations in strongly correlated Fermi systems. The fixed-node diffusion Monte Carlo method is used to produce the ground-state correlation functions which are used as an input for the excitation theory. Our approach reproduces recent Bragg scattering data in both the density and the spin channel. In the BCS regime, the response is close to that of the ideal Fermi gas. On the BEC side, the Bose peak associated with the formation of dimers dominates the density channel of the dynamic response. When the fraction of dimers is large our theory departs from the experimental data, mainly in the spin channel.

Quantum Monte Carlo Methods Describe Noncovalent Interactions with Subchemical Accuracy.

An accurate description of noncovalent interaction energies is one of the most challenging tasks in computational chemistry. To date, nonempirical CCSD(T)/CBS has been used as a benchmark reference. However, its practical use is limited due to the rapid growth of its computational cost with the system complexity. Here, we show that the fixed-node diffusion Monte Carlo (FN-DMC) method with a more favorable scaling is capable of reaching the CCSD(T)/CBS within subchemical accuracy (<0.1 kcal/mol) on a testing set of six small noncovalent complexes including the water dimer. In benzene/water, benzene/methane, and the T-shape benzene dimer, FN-DMC provides interaction energies that agree within 0.25 kcal/mol with the best available CCSD(T)/CBS estimates. The demonstrated predictive power of FN-DMC therefore provides new opportunities for studies of the vast and important class of medium/large noncovalent complexes.

Study of dipole moments of LiSr and KRb molecules by quantum Monte Carlo methods

Heteronuclear dimers are of significant interest to experiments seeking to exploit ultracold polar molecules in a number of novel ways, including precision measurement, quantum computing and quantum simulation. We calculate highly accurate Born–Oppenheimer total energies and electric dipole moments as a function of internuclear separation for two such dimers, LiSr and KRb. We apply fully correlated, high-accuracy quantum Monte Carlo methods for evaluating these molecular properties in a many-body framework. We use small-core effective potentials combined with multi-reference Slater–Jastrow trial wave functions to provide accurate nodes for the fixed-node diffusion Monte Carlo method. For reference and comparison, we calculate the same properties with Hartree–Fock and with restricted Configuration Interaction methods, and carefully assess the impact of the recovered many-body correlations on the calculated quantities. For LiSr, we find a highly non-linear dipole moment curve, which may make this molecule’s dipole moment tunable through vibrational state control.

Quantum Monte Carlo study of a resonant Bose-Fermi mixture.

We study a resonant Bose-Fermi mixture at zero temperature by using the fixed-node diffusion Monte Carlo method. We explore the system from weak to strong boson-fermion interaction, for different concentrations of the bosons relative to the fermion component. We focus on the case where the boson density n(B) is smaller than the fermion density n(F), for which a first-order quantum phase transition is found from a state with condensed bosons immersed in a Fermi sea, to a Fermi-Fermi mixture of composite fermions and unpaired fermions. We obtain the equation of state and the phase diagram, and we find that the region of phase separation shrinks to zero for vanishing n(B).

Ground state structure of BNB using fixed-node diffusion Monte Carlo

BNB is a challenging example of artifactual symmetry-breaking effects due to its susceptibility to a pseudo second-order Jahn–Teller interaction, which results in a structure with unequal BN bondlengths. The fixed-node diffusion Monte Carlo method is employed to calculate the potential curves along the symmetric and asymmetric stretching coordinates. With a multi-determinant wavefunction, the symmetric and asymmetric structures are found degenerate within statistical errors, with the asymmetric configuration lower in energy. The energy difference between the two structures is smaller than the multi-reference coupled cluster result obtained with four determinants, and supports previous conclusions that BNB has a floppy quasi-symmetric ground-state.