Transcorrelated method for electronic systems coupled with variational Monte Carlo calculation

Non-iterative Triples for Transcorrelated Coupled Cluster Theory

We present an implementation of a perturbative triples correction for the coupled cluster ansatz including single and double excitations based on the transcorrelated Hamiltonian. Transcorrelation introduces explicit electron correlation in the electronic Hamiltonian through similarity transformation with a correlation factor. Due to this transformation, the transcorrelated Hamiltonian includes up to three-body couplings and becomes non-Hermitian. Since the conventional coupled cluster equations are solved by projection, it is well suited to harbor non-Hermitian Hamiltonians. The arising three-body operator, however, creates a huge memory bottleneck and increases the runtime scaling of the coupled cluster equations. As it has been shown that the three-body operator can be approximated, by expressing the Hamiltonian in the normal-ordered form, we investigate this approximation for the perturbative triples correction. Results are compared with a code-generation based transcorrelated coupled cluster implementation up to quadruple excitations.

Striking the right balance of encoding electron correlation in the Hamiltonian and the wavefunction ansatz

Multi-configurational electronic structure theory delivers the most versatile approximations to many-electron wavefunctions, flexible enough to deal with all sorts of transformations, ranging from electronic excitations, to open-shell molecules and chemical reactions. Multi-configurational models are therefore essential to establish universally applicable, predictive ab initio methods for chemistry. Here, we present a discussion of explicit correlation approaches which address the nagging problem of dealing with static and dynamic electron correlation in multi-configurational active-space approaches. We review the latest developments and then point to their key obstacles. Our discussion is supported by new data obtained with tensor network methods. We argue in favor of simple electron-only correlator expressions that may allow one to define transcorrelated models in which the correlator does not bear a dependence on molecular structure.

Orbital optimisation in xTC transcorrelated methods

We present a combination of the bi-orthogonal orbital optimisation framework with the recently introduced xTC version of transcorrelation. This allows us to implement non-iterative perturbation based methods on top of the transcorrelated Hamiltonian. Additionally, the orbital optimisation influences results of other truncated methods, such as the distinguishable cluster with singles and doubles. The accuracy of these methods in comparison to standard xTC methods is demonstrated, and the advantages and disadvantages of the orbital optimisation are discussed.

Efficient Exploitation of Numerical Quadrature with Distance-Dependent Integral Screening in Explicitly Correlated F12 Theory: Linear Scaling Evaluation of the Most Expensive RI-MP2-F12 Term

We present a linear scaling atomic orbital based algorithm for the computation of the most expensive exchange-type RI-MP2-F12 term by employing numerical quadrature in combination with CABS-RI to avoid six-center-three-electron integrals. Furthermore, a robust distance-dependent integral screening scheme, based on integral partition bounds [Thompson, T. H.; Ochsenfeld, C. J. Chem. Phys. 2019, 150, 044101], is used to drastically reduce the number of the required three-center-one-electron integrals substantially. The accuracy of our numerical quadrature/CABS-RI approach and the corresponding integral screening is thoroughly assessed for interaction and isomerization energies across a variety of numerical integration grids. Our method outperforms the standard density fitting/CABS-RI approach with errors below 1 μEh even for small grid sizes and moderate screening thresholds. The choice of the grid size and screening threshold allows us to tailor our ansatz to a desired accuracy and computational efficiency. We showcase the approach's effectiveness for the chemically relevant system valinomycin, employing a triple-ζ F12 basis set combination (C54H90N6O18, 5757 AO basis functions, 10,266 CABS basis functions, 735,783 grid points). In this context, our ansatz achieves higher accuracy combined with a 135× speedup compared to the classical density fitting based variant, requiring notably less computation time than the corresponding RI-MP2 calculation. Additionally, we demonstrate near-linear scaling through calculations on linear alkanes. We achieved an 817-fold acceleration for C80H162 and an extrapolated 28,765-fold acceleration for C200H402, resulting in a substantially reduced computational time for the latter─from 229 days to just 11.5 min. Our ansatz may also be adapted to the remaining MP2-F12 terms, which will be the subject of future work.

Studies on the Transcorrelated Method

We investigate the possibility of using a transcorrelated (TC) Hamiltonian to describe electron correlation. A method to obtain TC wavefunctions was developed based on the mathematical framework of the bi-variational principle. This involves the construction of an effective TC Hamiltonian matrix, which can be solved in a self-consistent manner. This was optimized using a method we call second-order-moment minimization and demonstrate that it is possible to obtain highly accurate energies for some closed-shell atoms and helium-like ions. The effects of certain correlator terms on the description of electron-electron and electron-nuclear cusps were also examined graphically, and some TC wavefunctions were compared against near-exact Hylleraas wavefunctions.

Biorthonormal Orbital Optimization with a Cheap Core-Electron-Free Three-Body Correlation Factor for Quantum Monte Carlo and Transcorrelation

We introduce a novel three-body correlation factor that is designed to vanish in the core region around each nucleus and approach a universal two-body correlation factor for valence electrons. The transcorrelated Hamiltonian is used to optimize the orbitals of a single Slater determinant within a biorthonormal framework. The Slater-Jastrow wave function is optimized on a set of atomic and molecular systems containing both second-row elements and 3d transition metal elements. The optimization of the correlation factor and the orbitals, along with an increase in the basis set, results in a systematic lowering of the variational Monte Carlo energy for all systems tested. Importantly, the optimal parameters of the correlation factor obtained for atomic systems can be transferred to molecules. Additionally, the present correlation factor is computationally efficient and uses a mixed analytical-numerical integration scheme that reduces the costly numerical integration from R6 to R3.

Optimizing Jastrow factors for the transcorrelated method

We investigate the optimization of flexible tailored real-space Jastrow factors for use in the transcorrelated (TC) method in combination with highly accurate quantum chemistry methods, such as initiator full configuration interaction quantum Monte Carlo (FCIQMC). Jastrow factors obtained by minimizing the variance of the TC reference energy are found to yield better, more consistent results than those obtained by minimizing the variational energy. We compute all-electron atomization energies for the challenging first-row molecules C2, CN, N2, and O2 and find that the TC method yields chemically accurate results using only the cc-pVTZ basis set, roughly matching the accuracy of non-TC calculations with the much larger cc-pV5Z basis set. We also investigate an approximation in which pure three-body excitations are neglected from the TC-FCIQMC dynamics, saving storage and computational costs, and show that it affects relative energies negligibly. Our results demonstrate that the combination of tailored real-space Jastrow factors with the multi-configurational TC-FCIQMC method provides a route to obtaining chemical accuracy using modest basis sets, obviating the need for basis-set extrapolation and composite techniques.

Optimization of Large Determinant Expansions in Quantum Monte Carlo

We present a new method for the optimization of large configuration interaction (CI) expansions in the quantum Monte Carlo (QMC) framework. The central idea here is to replace the nonorthogonal variational optimization of CI coefficients performed in usual QMC calculations by an orthogonal non-Hermitian optimization thanks to the so-called transcorrelated (TC) framework, the two methods yielding the same results in the limit of a complete basis set. By rewriting the TC equations as an effective self-consistent Hermitian problem, our approach requires the sampling of a single quantity per Slater determinant, leading to minimal memory requirements in the QMC code. Using analytical quantities obtained from both the TC framework and the usual CI-type calculations, we also propose improved estimators which reduce the statistical fluctuations of the sampled quantities by more than an order of magnitude. We demonstrate the efficiency of this method on wave functions containing 10⁵-10⁶ Slater determinants, using effective core potentials or all-electron calculations. In all the cases, a sub-milli-Hartree convergence is reached within only two or three iterations of optimization.

Transcorrelated coupled cluster methods

Transcorrelated coupled cluster and distinguishable cluster methods are presented. The Hamiltonian is similarity transformed with a Jastrow factor in the first quantization, which results in up to three-body integrals. The coupled cluster with singles and doubles equations on this transformed Hamiltonian are formulated and implemented. It is demonstrated that the resulting methods have a superior basis set convergence and accuracy to the corresponding conventional and explicitly correlated methods. Additionally, approximations for three-body integrals are suggested and tested.

A new form of transcorrelated Hamiltonian inspired by range-separated DFT

The present work introduces a new form of explicitly correlated factor in the context of the transcorrelated methods. The new correlation factor is obtained from the r₁₂ ≈ 0 mathematical analysis of the transcorrelated Hamiltonian, and its analytical form is obtained such that the leading order in 1/r₁₂ of the scalar part of the effective two-electron potential reproduces the long-range interaction of the range-separated density functional theory. The resulting correlation factor exactly imposes the cusp and is tuned by a unique parameter μ, which controls both the depth of the coulomb hole and its typical range in r₁₂. The transcorrelated Hamiltonian obtained with such a new correlation factor has a straightforward analytical expression depending on the same parameter μ, and its physical contents continuously change by varying μ: One can change from a non-divergent repulsive Hamiltonian at large μ to a purely attractive one at small μ. We investigate the convergence of the ground state eigenvalues and right eigenvectors of such a new transcorrelated Hamiltonian as a function of the basis set and as a function of μ on a series of two-electron systems. We found that the convergence toward the complete basis set is much faster for quite a wide range of values of μ. We also propose a specific value of μ, which essentially reproduces the results obtained with the frozen Gaussian geminal introduced by Ten-no [Chem. Phys. Lett. 330, 169 (2000)].

σ-SCF: A direct energy-targeting method to mean-field excited states

The mean-field solutions of electronic excited states are much less accessible than ground state (e.g., Hartree-Fock) solutions. Energy-based optimization methods for excited states, like Δ-SCF (self-consistent field), tend to fall into the lowest solution consistent with a given symmetry-a problem known as "variational collapse." In this work, we combine the ideas of direct energy-targeting and variance-based optimization in order to describe excited states at the mean-field level. The resulting method, σ-SCF, has several advantages. First, it allows one to target any desired excited state by specifying a single parameter: a guess of the energy of that state. It can therefore, in principle, find all excited states. Second, it avoids variational collapse by using a variance-based, unconstrained local minimization. As a consequence, all states-ground or excited-are treated on an equal footing. Third, it provides an alternate approach to locate Δ-SCF solutions that are otherwise hardly accessible by the usual non-aufbau configuration initial guess. We present results for this new method for small atoms (He, Be) and molecules (H₂, HF). We find that σ-SCF is very effective at locating excited states, including individual, high energy excitations within a dense manifold of excited states. Like all single determinant methods, σ-SCF shows prominent spin-symmetry breaking for open shell states and our results suggest that this method could be further improved with spin projection.

Correlated Band Structure of a Transition Metal Oxide ZnO Obtained from a Many-Body Wave Function Theory

Obtaining accurate band structures of correlated solids has been one of the most important and challenging problems in first-principles electronic structure calculation. There have been promising recent active developments of wave function theory for condensed matter, but its application to band-structure calculation remains computationally expensive. In this Letter, we report the first application of the biorthogonal transcorrelated (BITC) method: self-consistent, free from adjustable parameters, and systematically improvable many-body wave function theory, to solid-state calculations with d electrons: wurtzite ZnO. We find that the BITC band structure better reproduces the experimental values of the gaps between the bands with different characters than several other conventional methods. This study paves the way for reliable first-principles calculations of the properties of strongly correlated materials.

Real-space finite-difference approach for multi-body systems: path-integral renormalization group method and direct energy minimization method

The path-integral renormalization group and direct energy minimization method of practical first-principles electronic structure calculations for multi-body systems within the framework of the real-space finite-difference scheme are introduced. These two methods can handle higher dimensional systems with consideration of the correlation effect. Furthermore, they can be easily extended to the multicomponent quantum systems which contain more than two kinds of quantum particles. The key to the present methods is employing linear combinations of nonorthogonal Slater determinants (SDs) as multi-body wavefunctions. As one of the noticeable results, the same accuracy as the variational Monte Carlo method is achieved with a few SDs. This enables us to study the entire ground state consisting of electrons and nuclei without the need to use the Born-Oppenheimer approximation. Recent activities on methodological developments aiming towards practical calculations such as the implementation of auxiliary field for Coulombic interaction, the treatment of the kinetic operator in imaginary-time evolutions, the time-saving double-grid technique for bare-Coulomb atomic potentials and the optimization scheme for minimizing the total-energy functional are also introduced. As test examples, the total energy of the hydrogen molecule, the atomic configuration of the methylene and the electronic structures of two-dimensional quantum dots are calculated, and the accuracy, availability and possibility of the present methods are demonstrated.

Quantum Monte Carlo study of first-row atoms using transcorrelated variational Monte Carlo trial functions

The effect of using the transcorrelated variational Monte Carlo (TC-VMC) approach to construct a trial function for fixed node diffusion Monte Carlo (DMC) energy calculations has been investigated for the first-row atoms, Li to Ne. The computed energies are compared with fixed node DMC energies obtained using trial functions constructed from Hartree-Fock and density functional levels of theory. Despite major VMC energy improvement with TC-VMC trial functions, no improvement in DMC energy was observed using these trial functions for the first-row atoms studied. The implications of these results on the nodes of the trial wave functions are discussed.

Correlation energy extrapolation by intrinsic scaling. V. Electronic energy, atomization energy, and enthalpy of formation of water

The method of correlation energy extrapolation by intrinsic scaling, recently introduced to obtain accurate molecular electronic energies, is used to calculate the total nonrelativistic electronic ground state energy of the water molecule. Accurate approximations to the full configuration interaction energies are determined for Dunning's [J. Chem. Phys. 90, 1007 (1989)] correlation-consistent double-, triple- and quadruple-zeta basis sets and then extrapolated to the complete basis set limit. The approach yields the total nonrelativistic energy -76.4390+/-0.0004 hartree, which compares very well with the value of -76.4389 hartree derived from experiment. The energy of atomization is recovered within 0.1 mh. The enthalpy of formation, which is obtained in conjunction with our previous calculation of the dissociation energy of the oxygen molecule, is recovered within 0.05 mh.

A practical treatment for the three-body interactions in the transcorrelated variational Monte Carlo method: Application to atoms from lithium to neon

We suggest a practical solution to dealing with the three-body interactions in the transcorrelated variational Monte Carlo method (TC-VMC). In the TC-VMC method, which was suggested in our previous paper [N. Umezawa and S. Tsuneyuki, J. Chem. Phys. 119, 10015 (2003)], the Jastrow-Slater-type wave function is efficiently optimized through a self-consistent procedure by minimizing the variance of the local energy. The three-body terms in the transcorrelated self-consistent-field equation, which have been simply ignored in our previous works, are efficiently calculated by the Monte Carlo numerical integration. We found that our treatment for the three-body interactions is successful for atoms from Li to Ne.

Correlation energy extrapolation by intrinsic scaling. IV. Accurate binding energies of the homonuclear diatomic molecules carbon, nitrogen, oxygen, and fluorine

The method of extrapolation by intrinsic scaling, recently introduced to obtain correlation energies, is generalized to multiconfigurational reference functions and used to calculate the binding energies of the diatomic molecules C2, N2, O2, and F2. First, accurate approximations to the full configuration interaction energies of the individual molecules and their constituent atoms are determined, employing Dunning's correlation consistent double-, triple- and quadruple zeta basis sets. Then, these energies are extrapolated to their full basis set limits. Chemical accuracy is attained for the binding energies of all molecules.

A ground state potential energy surface for H2 using Monte Carlo methods

Using variational Monte Carlo and a simple explicitly correlated wave function we have computed the Born-Oppenheimer energy of the H2 ground state (X 1Sigmag+) at 24 internuclear distances. We have also calculated the diagonal correction to the Born-Oppenheimer approximation and the lowest-order relativistic corrections at each distance using variational Monte Carlo techniques. The nonadiabatic values are evaluated from numerical derivatives of the wave function with respect to the nuclear coordinates. With this potential energy surface we have computed several of the lowest vibrational-rotational energies for this system. Our results are in good agreement with the best values found in the literature.

Excited electronic state calculations by the transcorrelated variational Monte Carlo method: Application to a helium atom

We have implemented the excited electronic state calculations for a helium atom by the transcorrelated variational Monte Carlo (TC-VMC) method. In this method, Jastrow-Slater-type wave function is efficiently optimized not only for the Jastrow factor but also for the Slater determinant. Since the formalism for the TC-VMC method is based on the variance minimization, excited states as well as the ground state calculations are feasible. It is found that both the first and the second excitation energies given by TC-VMC are much closer to the experimental data than those given by the variational Monte Carlo method with using the Hartree-Fock orbitals. The successful results in the TC-VMC method are considered to be due to the nodal optimization of the wave functions.

Optimized Jastrow–Slater wave functions for ground and excited states: Application to the lowest states of ethene

A quantum Monte Carlo method is presented for determining multideterminantal Jastrow-Slater wave functions for which the energy is stationary with respect to the simultaneous optimization of orbitals and configuration interaction coefficients. The approach is within the framework of the so-called energy fluctuation potential method which minimizes the energy in an iterative fashion based on Monte Carlo sampling and a fitting of the local energy fluctuations. The optimization of the orbitals is combined with the optimization of the configuration interaction coefficients through the use of additional single excitations to a set of external orbitals. A new set of orbitals is then obtained from the natural orbitals of this enlarged configuration interaction expansion. For excited states, the approach is extended to treat the average of several states within the same irreducible representation of the pointgroup of the molecule. The relationship of our optimization method with the stochastic reconfiguration technique by Sorella et al. is examined. Finally, the performance of our approach is illustrated with the lowest states of ethene, in particular with the difficult case of the 1(1)B(1u) state.

Explicitly correlated second order perturbation theory: Introduction of a rational generator and numerical quadratures

A rational generator, which fulfills the cusp conditions for singlet and triplet electron pairs, is proposed and applied to explicitly correlated second order Møller-Plesset perturbation theory calculations. It is shown that the generator in conjunction with frozen geminals improves the convergence of correlation energy without introducing any variational parameters in explicitly correlated functions. A new scheme for three-electron integrals based on numerical quadratures is also illustrated. The method is tested for the convergence of reaction enthalpies with various basis sets.