Jastrow-Slater Wave Functions Research Articles

Accurate Quantum Monte Carlo Forces for Machine-Learned Force Fields: Ethanol as a Benchmark.

Quantum Monte Carlo (QMC) is a powerful method to calculate accurate energies and forces for molecular systems. In this work, we demonstrate how we can obtain accurate QMC forces for the fluxional ethanol molecule at room temperature by using either multideterminant Jastrow-Slater wave functions in variational Monte Carlo or just a single determinant in diffusion Monte Carlo. The excellent performance of our protocols is assessed against high-level coupled cluster calculations on a diverse set of representative configurations of the system. Finally, we train machine-learning force fields on the QMC forces and compare them to models trained on coupled cluster reference data, showing that a force field based on the diffusion Monte Carlo forces with a single determinant can faithfully reproduce coupled cluster power spectra in molecular dynamics simulations.

Quantification of electron correlation for approximate quantum calculations.

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater-Jastrow wave functions.

Charge density waves in kagome-lattice extended Hubbard models at the van Hove filling

The Hubbard model on the kagome lattice is presently often considered as a minimal model to describe the rich low-temperature behavior of $A{\mathrm{V}}_{3}{\mathrm{Sb}}_{5}$ compounds (with $A=\mathrm{K},\mathrm{Rb},\mathrm{Cs}$), including charge density waves (CDWs), superconductivity, and possibly broken time-reversal symmetry. Here, we investigate, via variational Jastrow-Slater wave functions, the properties of its ground state when both on-site $U$ and nearest-neighbor $V$ Coulomb repulsions are considered at the van Hove filling. Our calculations reveal the presence of different interaction-driven CDWs and, contrary to previous renormalization-group studies, the absence of ferromagnetism and charge- or spin-bond order. No signatures of chiral phases are detected. Remarkably, the CDWs triggered by the nearest-neighbor repulsion possess charge disproportionations that are not compatible with the ones observed in $A{\mathrm{V}}_{3}{\mathrm{Sb}}_{5}$. As an alternative mechanism to stabilize charge-bond order, we consider the electron-phonon interaction, modeled by coupling the hopping amplitudes to quantum phonons, as in the Su-Schrieffer-Heeger model. Our results show the instability towards a trihexagonal distortion with $2\ifmmode\times\else\texttimes\fi{}2$ periodicity, in closer agreement with experimental findings.

Systematic lowering of the scaling of Monte Carlo calculations by partitioning and subsampling

We propose to compute physical properties by Monte Carlo calculations using conditional expectation values. The latter are obtained on top of the usual Monte Carlo sampling by partitioning the physical space in several subspaces or fragments, and subsampling each fragment (i.e., performing side walks) while freezing the environment. No bias is introduced and a zero-variance principle holds in the limit of separability, i.e., when the fragments are independent. In practice, the usual bottleneck of Monte Carlo calculations-the scaling of the statistical fluctuations as a function of the number of particles N-is relieved for extensive observables. We illustrate the method in variational Monte Carlo on the two-dimensional Hubbard model and on metallic hydrogen chains using Jastrow-Slater wave functions. A factor O(N) is gained in numerical efficiency.

Reference Excitation Energies of Increasingly Large Molecules: A QMC Study of Cyanine Dyes.

We revisit here the lowest vertical excitations of cyanine dyes using quantum Monte Carlo and leverage recent developments to systematically improve on previous results. In particular, we employ a protocol for the construction of compact and accurate multideterminant Jastrow-Slater wave functions for multiple states, which we have recently validated on the excited-state properties of several small prototypical molecules. Here, we obtain quantum Monte Carlo excitation energies in excellent agreement with high-level coupled cluster for all the cyanines where the coupled cluster method is applicable. Furthermore, we push our protocol to longer chains, demonstrating that quantum Monte Carlo is a viable methodology to establish reference data at system sizes which are hard to reach with other high-end approaches of similar accuracy. Finally, we determine which ingredients are key to an accurate treatment of these challenging systems and rationalize why a description of the excitation based on only active π orbitals lacks the desired accuracy for the shorter chains.

Extraction of a One-Particle Reduced Density Matrix from a Quantum Monte Carlo Electronic Density: A New Tool for Studying Nondynamic Correlation

In this work, we present a method to build a first order reduced density matrix (1-RDM) of a molecule from variational Quantum Monte Carlo (VMC) computations by means of a given correlated mapping wave function. Such a wave function is modeled on a Generalized Valence Bond plus Complete Active Space Configuration Interaction form and fits at best the density resulting from the Slater-Jastrow wave function of VMC. The accuracy of the method proposed has been proved by comparing the resulting kinetic energy with the corresponding VMC value. This 1-RDM is used to analyze the amount of correlation eventually captured in Kohn-Sham calculations performed in an unrestricted approach (UKS-DFT) and with different energy functionals. We performed test calculations on a selected set of molecules that show a significant multireference character. In this analysis, we compared both local and global indicators of nondynamic and dynamic correlation. Moreover, following the natural orbital decomposition of the 1-RDM, we also compared the effective temperatures of the corresponding Fermi-like distributions. Although there is a general agreement between UKS-DFT and VMC, we found the best match with the functional LC-BLYP.

Tailoring CIPSI Expansions for QMC Calculations of Electronic Excitations: The Case Study of Thiophene.

The perturbatively selected configuration interaction scheme (CIPSI) is particularly effective in constructing determinantal expansions for quantum Monte Carlo (QMC) simulations with Jastrow–Slater wave functions: fast and smooth convergence of ground-state properties and balanced descriptions of ground and excited states of different symmetries have been reported. In particular, accurate excitation energies have been obtained by the pivotal requirement of using CIPSI expansions with similar second-order perturbation corrections for each state, that is, a similar estimated distance to the full configuration interaction limit. Here, we elaborate on the CIPSI selection criterion for excited states of the same symmetry as the ground state, generating expansions from a common orbital set. Using these expansions in QMC as determinantal components of Jastrow–Slater wave functions, we compute the lowest, bright excited state of thiophene, which is challenging due to its significant multireference character. The resulting vertical excitation energies are within 0.05 eV of the best theoretical estimates, already with expansions of only a few thousand determinants. Furthermore, we relax the ground- and excited-state structures following the corresponding root in variational Monte Carlo and obtain bond lengths that are accurate to better than 0.01 Å. Therefore, while the full treatment at the CIPSI level of this system is quite demanding, in QMC, we can compute high-quality excitation energies and excited-state structural parameters building on affordable CIPSI expansions with relatively few, well-chosen determinants.

Ab Initio Dot Structures Beyond the Lewis Picture.

The empirical Lewis picture of the chemical bond dominates the view chemists have of molecules, of their stability and reactivity. Within the mathematical framework of quantum mechanics, all this chemical information is hidden in the many-particle wave function . Thus, to reveal and understand it, there is great interest in enhancing the Lewis model and connecting it to computable quantities. As has previously been shown, the Lewis picture can often be recovered from the probability density with probabilities in agreement with valence bond weights: the structures appear as most likely positions in the all-electron configuration space. Here, we systematically expand this topological probability density analysis to molecules with multiple bonds and lone pairs, employing correlated Slater-Jastrow wave functions. In contrast to earlier studies, non-Lewis structures are obtained that disagree with the prevalent picture and have a potentially better predictive capability. While functional groups are still recovered with these ab initio structures, the boundary between bonds and lone pairs is mostly blurred or non-existent. In order to understand the newly found structures, the Lewis electron pairs are replaced with spin-coupled electron motifs as the fundamental electronic fragment. These electron motifs—which coincide with Lewis’ electron pairs for many single bonds—arise naturally from the generally applicable analysis presented. An attempt is made to rationalize the geometry of the newly-found structures by considering the Coulomb force and the Pauli repulsion.

Quantum Monte Carlo Compton profiles of solid and liquid lithium

We computed the Compton profile of solid and liquid lithium using quantum Monte Carlo (QMC) and compared with recent experimental measurements obtaining good agreement. Importantly, we find it crucial to account for proper core-valence orthogonalization and to address density differences when comparing with experiment. To account for disorder effects, we sampled finite-temperature configurations using molecular dynamics (MD), then performed diffusion Monte Carlo (DMC) simulations on each configuration. We used Slater-Jastrow wavefunctions and grand-canonical twist-averaged boundary conditions. A QMC pseudopotential correction, derived from an all-electron DMC simulation of the perfect crystal was also used. Our calculations provide the first all-electron QMC benchmark for the Compton profile of lithium crystal and pseudopotential-corrected QMC Compton profiles for both the liquid and solid.

Wavepacket golden rule treatment of interparticle Coulombic decay in paired quantum dots.

The interparticle Coulombic decay process in paired quantum dots is studied by electron dynamics calculations. We consider a pair of Coulomb-coupled one-electron charged gallium arsenide quantum dots embedded in a nanowire. The two-electron decay process is approximately described by a single active electron model. Within this model, we employ the time-dependent wavepacket approach to the Fermi golden rule (introduced in the context of vibrational predissociation) to calculate autoionization rates, which are compared to exact rates obtained from fully correlated two-electron dynamics calculations. We found that the approximate decay rates agree well with the exact results in the limit of sufficiently separated quantum dots. Finally, we explore whether the short-range behavior of the new model can be further enhanced by the inclusion of local exchange effects by means of regularization of the Coulomb-potential based on a Jastrow-Slater wavefunction. The proposed method may open a route to study the interparticle Coulombic decay in more intricate systems, e.g., paired metal-nanoparticle-quantum dot systems.

Metallic and insulating stripes and their relation with superconductivity in the doped Hubbard model

The dualism between superconductivity and charge/spin modulations (the so-called stripes) dominates the phase diagram of many strongly-correlated systems. A prominent example is given by the Hubbard model, where these phases compete and possibly coexist in a wide regime of electron dopings for both weak and strong couplings. Here, we investigate this antagonism within a variational approach that is based upon Jastrow-Slater wave functions, including backflow correlations, which can be treated within a quantum Monte Carlo procedure. We focus on clusters having a ladder geometry with MM legs (with MM ranging from 22 to 1010) and a relatively large number of rungs, thus allowing us a detailed analysis in terms of the stripe length. We find that stripe order with periodicity \lambda=8λ=8 in the charge and 2\lambda=162λ=16 in the spin can be stabilized at doping \delta=1/8δ=1/8. Here, there are no sizable superconducting correlations and the ground state has an insulating character. A similar situation, with \lambda=6λ=6, appears at \delta=1/6δ=1/6. Instead, for smaller values of dopings, stripes can be still stabilized, but they are weakly metallic at \delta=1/12δ=1/12 and metallic with strong superconducting correlations at \delta=1/10δ=1/10, as well as for intermediate (incommensurate) dopings. Remarkably, we observe that spin modulation plays a major role in stripe formation, since it is crucial to obtain a stable striped state upon optimization. The relevance of our calculations for previous density-matrix renormalization group results and for the two-dimensional case is also discussed.

Excited States with Selected Configuration Interaction-Quantum Monte Carlo: Chemically Accurate Excitation Energies and Geometries.

We employ quantum Monte Carlo to obtain chemically accurate vertical and adiabatic excitation energies, and equilibrium excited-state structures for the small, yet challenging, formaldehyde and thioformaldehyde molecules. A key ingredient is a robust protocol to obtain balanced ground- and excited-state Jastrow–Slater wave functions at a given geometry, and to maintain such a balanced description as we relax the structure in the excited state. We use determinantal components generated via a selected configuration interaction scheme which targets the same second-order perturbation energy correction for all states of interest at different geometries, and fully optimize all variational parameters in the resultant Jastrow–Slater wave functions. Importantly, the excitation energies as well as the structural parameters in the ground and excited states are converged with very compact wave functions comprising few thousand determinants in a minimally augmented double-ζ basis set. These results are obtained already at the variational Monte Carlo level, the more accurate diffusion Monte Carlo method yielding only a small improvement in the adiabatic excitation energies. We find that matching Jastrow–Slater wave functions with similar variances can yield excitation energies compatible with our best estimates; however, the variance-matching procedure requires somewhat larger determinantal expansions to achieve the same accuracy, and it is less straightforward to adapt during structural optimization in the excited state.

Two-dimensional mixture of dipolar fermions: Equation of state and magnetic phases

We study a two-component mixture of fermionic dipoles in two dimensions at zero temperature, interacting via a purely repulsive $1/r^3$ potential. This model can be realized with ultracold atoms or molecules, when their dipole moments are aligned in the confinement direction orthogonal to the plane. We characterize the unpolarized mixture by means of the Diffusion Monte Carlo technique. Computing the equation of state, we identify the regime of validity for a mean-field theory based on a low-density expansion and compare our results with the hard-disk model of repulsive fermions. At high density, we address the possibility of itinerant ferromagnetism, namely whether the ground state can be fully polarized in the fluid phase. Within the fixed-node approximation, we show that the accuracy of Jastrow-Slater trial wave functions, even with the typical two-body backflow correction, is not sufficient to resolve the relevant energy differences. By making use of the iterative-backflow improved trial wave functions, we observe no signature of a fully-polarized ground state up to the freezing density.

Orbital-dependent backflow wave functions for real-space quantum Monte Carlo

We present and motivate an efficient way to include orbital dependent many--body correlations in trial wave function of real--space Quantum Monte Carlo methods for use in electronic structure calculations. We apply our new orbital--dependent backflow wave function to calculate ground state energies of the first row atoms using variational and diffusion Monte Carlo methods. The systematic overall gain of correlation energy with respect to single determinant Jastrow-Slater wave functions is competitive with the best single determinant trial wave functions currently available. The computational cost per Monte Carlo step is comparable to that of simple backflow calculations.

Metal-insulator transitions, superconductivity, and magnetism in the two-band Hubbard model

We explore the ground-state properties of the two-band Hubbard model with degenerate electronic bands, parametrized by nearest-neighbor hopping $t$, intra- and inter-orbital on-site Coulomb repulsions $U$ and $U^\prime$, and Hund coupling $J$, focusing on the case with $J>0$. Using Jastrow-Slater wave functions, we consider both states with and without magnetic/orbital order. Electron pairing can also be included in the wave function, in order to detect the occurrence of superconductivity for generic electron densities $n$. When no magnetic/orbital order is considered, the Mott transition is continuous for $n=1$ (quarter filling); instead, at $n=2$ (half filling), it is first order for small values of $J/U$, while it turns out to be continuous when the ratio $J/U$ is increased. A significant triplet pairing is present in a broad region around $n=2$. By contrast, singlet superconductivity (with $d$-wave symmetry) is detected only for small values of the Hund coupling and very close to half filling. When including magnetic and orbital order, the Mott insulator acquires antiferromagnetic order for $n=2$; instead, for $n=1$ the insulator has ferromagnetic and antiferro-orbital orders. In the latter case, a metallic phase is present for small values of $U/t$ and the metal-insulator transition becomes first order. In the region with $1<n<2$, we observe that ferromagnetism (with no orbital order) is particularly robust for large values of the Coulomb repulsion and that triplet superconductivity is strongly suppressed by the presence of antiferromagnetism. The case with $J=0$, which has an enlarged SU(4) symmetry due to the interplay between spin and orbital degrees of freedom, is also analyzed.

Perturbatively Selected Configuration-Interaction Wave Functions for Efficient Geometry Optimization in Quantum Monte Carlo.

We investigate the performance of a class of compact and systematically improvable Jastrow–Slater wave functions for the efficient and accurate computation of structural properties, where the determinantal component is expanded with a perturbatively selected configuration interaction scheme (CIPSI). We concurrently optimize the molecular ground-state geometry and full wave function—Jastrow factor, orbitals, and configuration interaction coefficients—in variational Monte Carlo (VMC) for the prototypical case of 1,3-trans-butadiene, a small yet theoretically challenging π-conjugated system. We find that the CIPSI selection outperforms the conventional scheme of correlating orbitals within active spaces chosen by chemical intuition: it gives significantly better variational and diffusion Monte Carlo energies for all but the smallest expansions, and much smoother convergence of the geometry with the number of determinants. In particular, the optimal bond lengths and bond-length alternation of butadiene are converged to better than 1 mÅ with just a few thousand determinants, to values very close to the corresponding CCSD(T) results. The combination of CIPSI expansion and VMC optimization represents an affordable tool for the determination of accurate ground-state geometries in quantum Monte Carlo.

Fixed-Node Diffusion Quantum Monte Carlo Method on Dissociation Energies and Their Trends for R-X Bonds (R = Me, Et, i-Pr, t-Bu).

Achieving both bond dissociation energies (BDEs) and their trends for the R-X bonds with R = Me, Et, i-Pr, and t-Bu reliably is nontrivial. Density functional theory (DFT) methods with traditional exchange-correlation functionals usually have large error on both the BDEs and their trends. The M06-2X functional gives rise to reliable BDEs, but the relative BDEs are determined not as accurately. More demanding approaches such as some double-hybrid functionals, for example, G4 and CCSD(T), are generally required to achieve the BDEs and their trends reliably. The fixed-node diffusion quantum Monte Carlo method (FN-DMC) is employed to calculated BDEs of these R-X bonds with X = H, CH3, OCH3, OH, and F in this work. The single Slater-Jastrow wave function is adopted as trial wave function, and pseudopotentials (PPs) developed for quantum Monte Carlo calculations are chosen. Error of these PPs is modest in wave function methods, while it is more pronounced in DFT calculations. Our results show that accuracy of BDEs with FN-DMC is similar to that of M06-2X and G4, and trends in BDEs are calculated more reliably than M06-2X. Both BDEs and trends in BDEs of these bonds are reproduced reasonably with FN-DMC. FN-DMC using PPs can thus be applied to BDEs and their trends of similar chemical bonds in larger molecules reliably and provide valuable information on properties of these molecules.

Optimizing the Energy with Quantum Monte Carlo: A Lower Numerical Scaling for Jastrow-Slater Expansions.

We present an improved formalism for quantum Monte Carlo calculations of energy derivatives and properties (e.g., the interatomic forces), with a multideterminant Jastrow-Slater function. As a function of the number Ne of Slater determinants, the numerical scaling of O(Ne) per derivative we have recently reported is here lowered to O(Ne) for the entire set of derivatives. As a function of the number of electrons N, the scaling to optimize the wave function and the geometry of a molecular system is lowered to O(N3) + O(NNe), the same as computing the energy alone in the sampling process. The scaling is demonstrated on linear polyenes up to C60H62 and the efficiency of the method is illustrated with the structural optimization of butadiene and octatetraene with Jastrow-Slater wave functions comprising as many as 200 000 determinants and 60 000 parameters.

Simple formalism for efficient derivatives and multi-determinant expansions in quantum Monte Carlo.

We present a simple and general formalism to compute efficiently the derivatives of a multi-determinant Jastrow-Slater wave function, the local energy, the interatomic forces, and similar quantities needed in quantum Monte Carlo. Through a straightforward manipulation of matrices evaluated on the occupied and virtual orbitals, we obtain an efficiency equivalent to algorithmic differentiation in the computation of the interatomic forces and the optimization of the orbital parameters. Furthermore, for a large multi-determinant expansion, the significant computational gain afforded by a recently introduced table method is here extended to the local value of any one-body operator and to its derivatives, in both all-electron and pseudopotential calculations.

Spontaneous symmetry breaking in correlated wave functions

We show that Jastrow-Slater wave functions, in which a density-density Jastrow factor is applied onto an uncorrelated fermionic state, may possess long-range order even when all symmetries are preserved in the wave function. This fact is mainly related to the presence of a sufficiently strong Jastrow term (also including the case of full Gutzwiller projection, suitable for describing spin models). Selected examples are reported, including the spawning of N\'eel order and dimerization in spin systems, and the stabilization of charge and orbital order in itinerant electronic systems.