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.

TurboGenius: Python suite for high-throughput calculations of ab initio quantum Monte Carlo methods.

TurboGenius is an open-source Python package designed to fully control ab initio quantum Monte Carlo (QMC) jobs using a Python script, which allows one to perform high-throughput calculations combined with TurboRVB [Nakano et al. J. Phys. Chem. 152, 204121 (2020)]. This paper provides an overview of the TurboGenius package and showcases several results obtained in a high-throughput mode. For the purpose of performing high-throughput calculations with TurboGenius, we implemented another open-source Python package, TurboWorkflows, that enables one to construct simple workflows using TurboGenius. We demonstrate its effectiveness by performing (1) validations of density functional theory (DFT) and QMC drivers as implemented in the TurboRVB package and (2) benchmarks of Diffusion Monte Carlo (DMC) calculations for several datasets. For (1), we checked inter-package consistencies between TurboRVB and other established quantum chemistry packages. By doing so, we confirmed that DFT energies obtained by PySCF are consistent with those obtained by TurboRVB within the local density approximation (LDA) and that Hartree-Fock (HF) energies obtained by PySCF and Quantum Package are consistent with variational Monte Carlo energies obtained by TurboRVB with the HF wavefunctions. These validation tests constitute a further reliability check of the TurboRVB package. For (2), we benchmarked the atomization energies of the Gaussian-2 set, the binding energies of the S22, A24, and SCAI sets, and the equilibrium lattice parameters of 12 cubic crystals using DMC calculations. We found that, for all compounds analyzed here, the DMC calculations with the LDA nodal surface give satisfactory results, i.e., consistent either with high-level computational or with experimental reference values.

TREXIO: A file format and library for quantum chemistry.

TREXIO is an open-source file format and library developed for the storage and manipulation of data produced by quantum chemistry calculations. It is designed with the goal of providing a reliable and efficient method of storing and exchanging wave function parameters and matrix elements, making it an important tool for researchers in the field of quantum chemistry. In this work, we present an overview of the TREXIO file format and library. The library consists of a front-end implemented in the C programming language and two different back-ends: a text back-end and a binary back-end utilizing the hierarchical data format version 5 library, which enables fast read and write operations. It is compatible with a variety of platforms and has interfaces for Fortran, Python, and OCaml programming languages. In addition, a suite of tools have been developed to facilitate the use of the TREXIO format and library, including converters for popular quantum chemistry codes and utilities for validating and manipulating data stored in TREXIO files. The simplicity, versatility, and ease of use of TREXIO make it a valuable resource for researchers working with quantum chemistry data.

Toward Chemical Accuracy Using the Jastrow Correlated Antisymmetrized Geminal Power Ansatz.

Herein, we report accurate atomization energy calculations for 55 molecules in the Gaussian-2 (G2) set using lattice regularized diffusion Monte Carlo (LRDMC). We compare the Jastrow-Slater determinant ansatz with a more flexible JsAGPs (Jastrow correlated antisymmetrized geminal power with singlet correlation) ansatz. AGPs is built from pairing functions, which explicitly include pairwise correlations among electrons, and hence, this ansatz is expected to be more efficient in recovering the correlation energy. The AGPs wave functions are first optimized at the variational Monte Carlo (VMC) level, which includes both the Jastrow factor and the nodal surface optimization. This is followed by the LRDMC projection of the ansatz. Remarkably, for many molecules, the LRDMC atomization energies obtained using the JsAGPs ansatz reach chemical accuracy (∼1 kcal/mol), and for most other molecules, the atomization energies are accurate within ∼5 kcal/mol. We obtained a mean absolute deviation of 1.6 kcal/mol with JsAGPs and 3.2 kcal/mol with JDFT (Jastrow factor + Slater determinant with DFT orbitals) ansatzes. This work shows the effectiveness of the flexible AGPs ansatz for atomization energy calculations and electronic structure simulations in general.

Space-warp coordinate transformation for efficient ionic force calculations in quantum Monte Carlo.

Ab initio quantum Monte Carlo (QMC) methods are a state-of-the-art computational approach to obtaining highly accurate many-body wave functions. Although QMC methods are widely used in physics and chemistry to compute ground-state energies, calculation of atomic forces is still under technical/algorithmic development. Very recently, force evaluation has started to become of paramount importance for the generation of machine-learning force-field potentials. Nevertheless, there is no consensus regarding whether an efficient algorithm is available for the QMC force evaluation, namely, one that scales well with the number of electrons and the atomic numbers. In this study, we benchmark the accuracy of all-electron variational Monte Carlo (VMC) and lattice-regularized diffusion Monte Carlo (LRDMC) forces for various mono- and heteronuclear dimers (1 ≤ Z ≤ 35, where Z is the atomic number). The VMC and LRDMC forces were calculated with and without the so-called space-warp coordinate transformation (SWCT) and appropriate regularization techniques to remove the infinite variance problem. The LRDMC forces were computed with the Reynolds (RE) and variational-drift (VD) approximations. The potential energy surfaces obtained from the LRDMC energies give equilibrium bond lengths (req) and harmonic frequencies (ω) very close to the experimental values for all dimers, improving the corresponding VMC results. The LRDMC forces with the RE approximation improve the VMC forces, implying that it is worth computing the DMC forces beyond VMC despite the higher computational cost. The LRDMC forces with the VD approximations also show improvement, which unfortunately comes at a much higher computational cost in all-electron calculations. We find that the ratio of computational costs between QMC energy and forces scales as Z∼2.5 without the SWCT. In contrast, the application of the SWCT makes the ratio independent of Z. As such, the accessible QMC system size is not affected by the evaluation of ionic forces but governed by the same scaling as the total energy one.

Magnetism and Charge Order in the Honeycomb Lattice.

Despite being relevant to better understand the properties of honeycomblike systems, as graphene-based compounds, the electron-phonon interaction is commonly disregarded in theoretical approaches. That is, the effects of phonon fields on interacting Dirac electrons is an open issue, in particular when investigating long-range ordering. Thus, here we perform unbiased quantum Monte Carlo simulations to examine the Hubbard-Holstein model (HHM) in the half-filled honeycomb lattice. By performing careful finite-size scaling analysis, we identify semimetal-to-insulator quantum critical points, and determine the behavior of the antiferromagnetic and charge-density wave phase transitions. We have, therefore, established the ground state phase diagram of the HHM for intermediate interaction strength, determining its behavior for different phonon frequencies. Our findings provide quantitative and qualitative descriptions of the model at intermediate coupling strengths, and may shed light on the emergence of many-body properties in honeycomblike systems.

The nature of the chemical bond in the dicarbon molecule.

The molecular dissociation energy has often been explained and discussed in terms of singlet bonds, formed by bounded pairs of valence electrons. In this work, we use a highly correlated resonating valence bond ansatz, providing a consistent paradigm for the chemical bond, where spin fluctuations are shown to play a crucial role. Spin fluctuations are known to be important in magnetic systems and correspond to the zero point motion of the spin waves emerging from a magnetic broken symmetry state. Within our ansatz, a satisfactory description of the carbon dimer is determined by the magnetic interaction of two carbon atoms with antiferromagnetically ordered S = 1 magnetic moments. This is a first step that, thanks to the highly scalable and efficient quantum Monte Carlo techniques, may open the door for understanding challenging complex systems containing atoms with large spins (e.g., transition metals).

General Correlated Geminal Ansatz for Electronic Structure Calculations: Exploiting Pfaffians in Place of Determinants.

We propose here a single Pfaffian correlated variational ansatz that dramatically improves the accuracy with respect to the single determinant one, while remaining at a similar computational cost. A much larger correlation energy is indeed determined by the most general two electron pairing function, including both singlet and triplet channels, combined with a many-body Jastrow factor, including all possible spin-spin, spin-density, and density-density terms. The main technical ingredient to exploit this accuracy is the use of the Pfaffian for antisymmetrizing a highly correlated pairing function, thus recovering the Fermi statistics for electrons with an affordable computational cost. Moreover, the application of the diffusion Monte Carlo, within the fixed node approximation, allows us to obtain very accurate binding energies for the first preliminary calculations reported in this study: C2, N2, and O2 and the benzene molecule. This is promising and remarkable, considering that they represent extremely difficult molecules even for computationally demanding multideterminant approaches, and opens therefore the way for realistic and accurate electronic simulations with an algorithm scaling at most as the fourth power of the number of electrons.

TurboRVB: A many-body toolkit for ab initio electronic simulations by quantum Monte Carlo.

TurboRVB is a computational package for ab initio Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC) and diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions (WFs), capable of describing several materials with very high accuracy, even when standard mean-field approaches [e.g., density functional theory (DFT)] fail. The electronic WF is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal power with spin-singlet pairing or as a Pfaffian, including both singlet and triplet correlations. This WF can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) Ansatz, first proposed by Pauling and Anderson in quantum chemistry [L. Pauling, The Nature of the Chemical Bond (Cornell University Press, 1960)] and condensed matter physics [P.W. Anderson, Mat. Res. Bull 8, 153 (1973)], respectively. The RVB Ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. Therefore, its application to large systems is computationally feasible. The WF is expanded in a localized basis set. Several basis set functions are implemented, such as Gaussian, Slater, and mixed types, with no restriction on the choice of their contraction. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization (Jastrow and antisymmetric parts together) is made possible, thanks to state-of-the-art stochastic algorithms for energy minimization. In the optimization procedure, the first guess can be obtained at the mean-field level by a built-in DFT driver. The code has been efficiently parallelized by using a hybrid MPI-OpenMP protocol, which is also an ideal environment for exploiting the computational power of modern Graphics Processing Unit accelerators.

All-Electron Quantum Monte Carlo with Jastrow Single Determinant Ansatz: Application to the Sodium Dimer.

In this work, we report potential energy surfaces (PESs) of the sodium dimer calculated by variational (VMC) and lattice-regularized diffusion Monte Carlo (LRDMC). The VMC calculation is accurate for determining the equilibrium distance and the qualitative shape of the experimental PES. Remarkably, after the application of the LRDMC projection to this single determinant ansatz, namely, the Jastrow Antisymmetrized Geminal Power (JAGP), chemical accuracy (∼1 kcal/mol) is reached in the binding energy, and the obtained equilibrium internuclear distance and harmonic vibrational frequency are in very good agreement with the experimental ones. This outcome is crucially dependent on the quality of the optimization used to determine the best possible trial function within the chosen ansatz. The strategy adopted in this work is to minimize the variational energy by initializing the trial function with the density functional theory (DFT) single determinant ansatz expanded exactly in the same atomic basis used for the corresponding VMC and LRDMC calculations. This atomic basis is reshaped ad-hoc for QMC calculations. Indeed, we multiply the standard Gaussian-type atomic orbitals by a one-body Jastrow factor, satisfying, in this way, the electron-ion cusp conditions. In order to achieve these important advantages, we have defined a very efficient DFT algorithm in the mentioned basis, by estimating the corresponding matrix elements on a mesh, and by using a much finer mesh grid in the vicinity of nuclei.

Photoinduced η Pairing in the Hubbard Model.

By employing unbiased numerical methods, we show that pulse irradiation can induce unconventional superconductivity even in the Mott insulator of the Hubbard model. The superconductivity found here in the photoexcited state is due to the η-pairing mechanism, characterized by staggered pair-density-wave oscillations in the off-diagonal long-range correlation, and is absent in the ground-state phase diagram; i.e., it is induced neither by a change of the effective interaction of the Hubbard model nor by simple photocarrier doping. Because of the selection rule, we show that the nonlinear optical response is essential to increase the number of η pairs and thus enhance the superconducting correlation in the photoexcited state. Our finding demonstrates that nonequilibrium many-body dynamics is an alternative pathway to access a new exotic quantum state that is absent in the ground-state phase diagram, and also provides an alternative mechanism for enhancing superconductivity.

Angle-resolved photoemission spectroscopy from first-principles quantum Monte Carlo.

Angle-resolved photoemission spectroscopy allows one to visualize in momentum space the probability weight maps of electrons subtracted from molecules deposited on a substrate. The interpretation of these maps usually relies on the plane wave approximation through the Fourier transform of single particle orbitals obtained from density functional theory. Here we propose a first-principle many-body approach based on quantum Monte Carlo (QMC) to directly calculate the quasi-particle wave functions (also known as Dyson orbitals) of molecules in momentum space. The comparison between these correlated QMC images and their single particle counterpart highlights features that arise from many-body effects. We test the QMC approach on the linear C2H2, CO2, and N2 molecules, for which only small amplitude remodulations are visible. Then, we consider the case of the pentacene molecule, focusing on the relationship between the momentum space features and the real space quasi-particle orbital. Eventually, we verify the correlation effects present in the metal CuCl 4 2 - planar complex.

Correlation-Driven Dimerization and Topological Gap Opening in Isotropically Strained Graphene.

The phase diagram of isotropically expanded graphene cannot be correctly predicted by ignoring either electron correlations, or mobile carbons, or the effect of applied stress, as was done so far. We calculate the ground state enthalpy (not just energy) of strained graphene by an accurate off-lattice quantum Monte Carlo correlated ansatz of great variational flexibility. Following undistorted semimetallic graphene at low strain, multideterminant Heitler-London correlations stabilize between ≃8.5% and ≃15% strain an insulating Kekulé-like dimerized (DIM) state. Closer to a crystallized resonating-valence bond than to a Peierls state, the DIM state prevails over the competing antiferromagnetic insulating state favored by density-functional calculations which we conduct in parallel. The DIM stressed graphene insulator, whose gap is predicted to grow in excess of 1 eV before failure near 15% strain, is topological in nature, implying under certain conditions 1D metallic interface states lying in the bulk energy gap.

Phase Diagram of Hydrogen and a Hydrogen-Helium Mixture at Planetary Conditions by Quantum Monte Carlo Simulations.

Understanding planetary interiors is directly linked to our ability of simulating exotic quantum mechanical systems such as hydrogen (H) and hydrogen-helium (H-He) mixtures at high pressures and temperatures. Equation of state (EOS) tables based on density functional theory are commonly used by planetary scientists, although this method allows only for a qualitative description of the phase diagram. Here we report quantum Monte Carlo (QMC) molecular dynamics simulations of pure H and H-He mixture. We calculate the first QMC EOS at 6000 K for a H-He mixture of a protosolar composition, and show the crucial influence of He on the H metallization pressure. Our results can be used to calibrate other EOS calculations and are very timely given the accurate determination of Jupiter's gravitational field from the NASA Juno mission and the effort to determine its structure.

Carbon nanotubes as excitonic insulators.

Fifty years ago Walter Kohn speculated that a zero-gap semiconductor might be unstable against the spontaneous generation of excitons-electron-hole pairs bound together by Coulomb attraction. The reconstructed ground state would then open a gap breaking the symmetry of the underlying lattice, a genuine consequence of electronic correlations. Here we show that this excitonic insulator is realized in zero-gap carbon nanotubes by performing first-principles calculations through many-body perturbation theory as well as quantum Monte Carlo. The excitonic order modulates the charge between the two carbon sublattices opening an experimentally observable gap, which scales as the inverse of the tube radius and weakly depends on the axial magnetic field. Our findings call into question the Luttinger liquid paradigm for nanotubes and provide tests to experimentally discriminate between excitonic and Mott insulators.

Fully Quantum Description of the Zundel Ion: Combining Variational Quantum Monte Carlo with Path Integral Langevin Dynamics.

We introduce a novel approach for a fully quantum description of coupled electron-ion systems from first principles. It combines the variational quantum Monte Carlo solution of the electronic part with the path integral formalism for the quantum nuclear dynamics. On the one hand, the path integral molecular dynamics includes nuclear quantum effects by adding a set of fictitious classical particles (beads) aimed at reproducing nuclear quantum fluctuations via a harmonic kinetic term. On the other hand, variational quantum Monte Carlo can provide Born-Oppenheimer potential energy surfaces with a precision comparable to the most-advanced post-Hartree-Fock approaches, and with a favorable scaling with the system size. In order to cope with the intrinsic noise due to the stochastic nature of quantum Monte Carlo methods, we generalize the path integral molecular dynamics using a Langevin thermostat correlated according to the covariance matrix of quantum Monte Carlo nuclear forces. The variational parameters of the quantum Monte Carlo wave function are evolved during the nuclear dynamics, such that the Born-Oppenheimer potential energy surface is unbiased. Statistical errors on the wave function parameters are reduced by resorting to bead grouping average, which we show to be accurate and well-controlled. Our general algorithm relies on a Trotter breakup between the dynamics driven by ionic forces and the one set by the harmonic interbead couplings. The latter is exactly integrated, even in the presence of the Langevin thermostat, thanks to the mapping onto an Ornstein-Uhlenbeck process. This framework turns out to be also very efficient in the case of noiseless (deterministic) ionic forces. The new implementation is validated on the Zundel ion (H5O2+) by direct comparison with standard path integral Langevin dynamics calculations made with a coupled cluster potential energy surface. Nuclear quantum effects are confirmed to be dominant over thermal effects well beyond room temperature, giving the excess proton an increased mobility by quantum tunneling.

Accelerating ab initio Molecular Dynamics and Probing the Weak Dispersive Forces in Dense Liquid Hydrogen.

We propose an ab initio molecular dynamics method, capable of dramatically reducing the autocorrelation time required for the simulation of classical and quantum particles at finite temperatures. The method is based on an efficient implementation of a first order Langevin dynamics modified by means of a suitable, position dependent acceleration matrix S. Here, we apply this technique to both Lennard-Jones models, to demonstrate the accuracy and speeding-up of the sampling, and within a quantum Monte Carlo based wave function approach, for determining the phase diagram of high-pressure hydrogen with simulations much longer than the autocorrelation time. With the proposed method, we are able to equilibrate in a few hundred steps even close to the liquid-liquid phase transition (LLT). Within our approach, we find that the LLT transition is consistent with recent density functionals predicting a much larger transition pressure when the long range dispersive forces are taken into account.

Toward Accurate Adsorption Energetics on Clay Surfaces.

Clay minerals are ubiquitous in nature, and the manner in which they interact with their surroundings has important industrial and environmental implications. Consequently, a molecular-level understanding of the adsorption of molecules on clay surfaces is crucial. In this regard computer simulations play an important role, yet the accuracy of widely used empirical force fields (FF) and density functional theory (DFT) exchange-correlation functionals is often unclear in adsorption systems dominated by weak interactions. Herein we present results from quantum Monte Carlo (QMC) for water and methanol adsorption on the prototypical clay kaolinite. To the best of our knowledge, this is the first time QMC has been used to investigate adsorption at a complex, natural surface such as a clay. As well as being valuable in their own right, the QMC benchmarks obtained provide reference data against which the performance of cheaper DFT methods can be tested. Indeed using various DFT exchange-correlation functionals yields a very broad range of adsorption energies, and it is unclear a priori which evaluation is better. QMC reveals that in the systems considered here it is essential to account for van der Waals (vdW) dispersion forces since this alters both the absolute and relative adsorption energies of water and methanol. We show, via FF simulations, that incorrect relative energies can lead to significant changes in the interfacial densities of water and methanol solutions at the kaolinite interface. Despite the clear improvements offered by the vdW-corrected and the vdW-inclusive functionals, absolute adsorption energies are often overestimated, suggesting that the treatment of vdW forces in DFT is not yet a solved problem.

Correlation Effects in Scanning Tunneling Microscopy Images of Molecules Revealed by Quantum Monte Carlo.

Scanning tunneling microscopy (STM) and spectroscopy probe the local density of states of single molecules electrically insulated from the substrate. The experimental images, although usually interpreted in terms of single-particle molecular orbitals, are associated with quasiparticle wave functions dressed by the whole electron-electron interaction. Here we propose an ab initio approach based on quantum Monte Carlo to calculate the quasiparticle wave functions of molecules. Through the comparison between Monte Carlo wave functions and their uncorrelated Hartree-Fock counterparts we visualize the electronic correlation embedded in the simulated STM images, highlighting the many-body features that might be observed.

Assessing the orbital selective Mott transition with variational wave functions.

We study the Mott metal-insulator transition in the two-band Hubbard model with different hopping amplitudes t1 and t2 for the two orbitals on the two-dimensional square lattice by using non-magnetic variational wave functions, similarly to what has been considered in the limit of infinite dimensions by dynamical mean-field theory. We work out the phase diagram at half filling (i.e. two electrons per site) as a function of R = t2/t1 and the on-site Coulomb repulsion U, for two values of the Hund's coupling J = 0 and J/U = 0.1. Our results are in good agreement with previous dynamical mean-field theory calculations, demonstrating that the non-magnetic phase diagram is only slightly modified from infinite to two spatial dimensions. Three phases are present: a metallic one, for small values of U, where both orbitals are itinerant; a Mott insulator, for large values of U, where both orbitals are localized because of the Coulomb repulsion; and the so-called orbital-selective Mott insulator (OSMI), for small values of R and intermediate Us, where one orbital is localized while the other one is still itinerant. The effect of the Hund's coupling is two-fold: on one side, it favors the full Mott phase over the OSMI; on the other side, it stabilizes the OSMI at larger values of R.