Framework for constructing generic Jastrow correlation factors

Toward Real Chemical Accuracy on Current Quantum Hardware Through the Transcorrelated Method

Quantum computing is emerging as a new computational paradigm with the potential to transform several research fields including quantum chemistry. However, current hardware limitations (including limited coherence times, gate infidelities, and connectivity) hamper the implementation of most quantum algorithms and call for more noise-resilient solutions. We propose an explicitly correlated Ansatz based on the transcorrelated (TC) approach to target these major roadblocks directly. This method transfers, without any approximation, correlations from the wave function directly into the Hamiltonian, thus reducing the resources needed to achieve accurate results with noisy quantum devices. We show that the TC approach allows for shallower circuits and improves the convergence toward the complete basis set limit, providing energies within chemical accuracy to experiment with smaller basis sets and, thus, fewer qubits. We demonstrate our method by computing bond lengths, dissociation energies, and vibrational frequencies close to experimental results for the hydrogen dimer and lithium hydride using two and four qubits, respectively. To demonstrate our approach's current and near-term potential, we perform hardware experiments, where our results confirm that the TC method paves the way toward accurate quantum chemistry calculations already on today's quantum hardware.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Transcorrelated selected configuration interaction in a bi-orthonormal basis and with a cheap three-body correlation factor

In this work, we develop a mathematical framework for a selected configuration interaction (SCI) algorithm within a bi-orthogonal basis for transcorrelated (TC) calculations. The bi-orthogonal basis used here serves as the equivalent of the standard Hartree-Fock (HF) orbitals. However, within the context of TC, it leads to distinct orbitals for the left and right vectors. Our findings indicate that the use of such a bi-orthogonal basis allows for a proper definition of the frozen core approximation. In contrast, the use of HF orbitals results in bad error cancellations for ionization potentials and atomization energies (AE). Compared to HF orbitals, the optimized bi-orthogonal basis significantly reduces the positive part of the second-order energy (PT2), thereby facilitating the use of standard extrapolation techniques of hermitian SCI. While we did not observe a significant improvement in the convergence of the SCI algorithm, this is largely due to the use in this work of a simple three-body correlation factor introduced in a recent study. This correlation factor, which depends only on atomic parameters, eliminates the need for re-optimization of the correlation factor for molecular systems, making its use straightforward and user-friendly. Despite the simplicity of this correlation factor, we were able to achieve accurate results on the AE of a series of 14 molecules on a triple-zeta basis. We also successfully broke a double bond until the full dissociation limit while maintaining the size consistency property. This work thus demonstrates the potential of the BiO-TC-SCI approach in handling complex molecular systems.

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.

Density Matrix Renormalization Group for Transcorrelated Hamiltonians: Ground and Excited States in Molecules

We present the theory of a density matrix renormalization group (DMRG) algorithm which can solve for both the ground and excited states of non-Hermitian transcorrelated Hamiltonians and show applications in molecular systems. Transcorrelation (TC) accelerates the basis set convergence rate by including known physics (such as, but not limited to, the electron-electron cusp) in the Jastrow factor used for the similarity transformation. It also improves the accuracy of approximate methods such as coupled cluster singles and doubles (CCSD) as shown by recent studies. However, the non-Hermiticity of the TC Hamiltonians poses challenges for variational methods like DMRG. Imaginary-time evolution on the matrix product state (MPS) in the DMRG framework has been proposed to circumvent this problem, but this is currently limited to treating the ground state and has lower efficiency than the time-independent DMRG (TI-DMRG) due to the need to eliminate Trotter errors. In this work, we show that with minimal changes to the existing TI-DMRG algorithm, namely, replacing the original Davidson solver with the general Davidson solver to solve the non-Hermitian effective Hamiltonians at each site for a few low-lying right eigenstates, and following the rest of the original DMRG recipe, one can find the ground and excited states with improved efficiency compared to the original DMRG when extrapolating to the infinite bond dimension limit in the same basis set. An accelerated basis set convergence rate is also observed, as expected, within the TC framework.

A hybrid approach to excited-state-specific variational Monte Carlo and doubly excited states

We extend our hybrid linear-method/accelerated-descent variational Monte Carlo optimization approach to excited states and investigate its efficacy in double excitations. In addition to showing a superior statistical efficiency when compared to the linear method, our tests on small molecules show good energetic agreement with benchmark methods. We also demonstrate the ability to treat double excitations in systems that are too large for a full treatment by using selected configuration interaction methods via an application to 4-aminobenzonitrile. Finally, we investigate the stability of state-specific variance optimization against collapse to other states' variance minima and find that symmetry, Ansatz quality, and sample size all have roles to play in achieving stability.

Tail-regression estimator for heavy-tailed distributions of known tail indices and its application to continuum quantum Monte Carlo data

Standard statistical analysis is unable to provide reliable confidence intervals on expectation values of probability distributions that do not satisfy the conditions of the central limit theorem. We present a regression-based estimator of an arbitrary moment of a probability distribution with power-law heavy tails that exploits knowledge of the exponents of its asymptotic decay to bypass this issue entirely. Our method is applied to synthetic data and to energy and atomic force data from variational and diffusion quantum Monte Carlo calculations, whose distributions have known asymptotic forms [J. R. Trail, Phys. Rev. E 77, 016703 (2008)PLEEE81539-375510.1103/PhysRevE.77.016703; A. Badinski et al., J. Phys.: Condens. Matter 22, 074202 (2010)JCOMEL0953-898410.1088/0953-8984/22/7/074202]. We obtain convergent, accurate confidence intervals on the variance of the local energy of an electron gas and on the Hellmann-Feynman force on an atom in the all-electron carbon dimer. In each of these cases the uncertainty on our estimator is 45% and 60 times smaller, respectively, than the nominal (ill-defined) standard error.

Complementary first and second derivative methods for ansatz optimization in variational Monte Carlo

This Perspective contrasts first and second derivative methods in variational Monte Carlo and presents a hybrid optimization approach that combines their advantages.

Phase diagram of a symmetric electron-hole bilayer system: a variational Monte Carlo study

We study the phase diagram of a symmetric electron-hole bilayer system at absolute zero temperature and in zero magnetic field within the quantum Monte Carlo approach. In particular, we conduct variational Monte Carlo simulations for various phases, i.e. the paramagnetic fluid phase, the ferromagnetic fluid phase, the anti-ferromagnetic Wigner crystal phase, the ferromagnetic Wigner crystal phase and the excitonic phase, to estimate the ground-state energy at different values of in-layer density and inter-layer spacing. Slater-Jastrow style trial wave functions, with single-particle orbitals appropriate for different phases, are used to construct the phase diagram in the (r _s , d) plane by finding the relative stability of trial wave functions. At very small layer separations, we find that the fluid phases are stable, with the paramagnetic fluid phase being particularly stable at [Formula: see text] and the ferromagnetic fluid phase being particularly stable at [Formula: see text]. As the layer spacing increases, we first find that there is a phase transition from the ferromagnetic fluid phase to the ferromagnetic Wigner crystal phase when d reaches 0.4 a.u. at r _s   =  20, and before there is a return to the ferromagnetic fluid phase when d approaches 1 a.u. However, for r _s   <  20 and [Formula: see text] a.u., the excitonic phase is found to be stable. We do not find that the anti-ferromagnetic Wigner crystal is stable over the considered range of r _s and d. We also find that as r _s increases, the critical layer separations for Wigner crystallization increase.

Evidence from Quantum Monte Carlo Simulations of Large-Gap Superfluidity and BCS-BEC Crossover in Double Electron-Hole Layers

We report quantum Monte Carlo evidence of the existence of large gap superfluidity in electron-hole double layers over wide density ranges. The superfluid parameters evolve from normal state to BEC with decreasing density, with the BCS state restricted to a tiny range of densities due to the strong screening of Coulomb interactions, which causes the gap to rapidly become large near the onset of superfluidity. The superfluid properties exhibit similarities to ultracold fermions and iron-based superconductors, suggesting an underlying universal behavior of BCS-BEC crossovers in pairing systems.

Generic expansion of the Jastrow correlation factor in polynomials satisfying symmetry and cusp conditions

Jastrow correlation factors play an important role in quantum Monte Carlo calculations. Together with an orbital based antisymmetric function, they allow the construction of highly accurate correlation wave functions. In this paper, a generic expansion of the Jastrow correlation function in terms of polynomials that satisfy both the electron exchange symmetry constraint and the cusp conditions is presented. In particular, an expansion of the three-body electron-electron-nucleus contribution in terms of cuspless homogeneous symmetric polynomials is proposed. The polynomials can be expressed in fairly arbitrary scaling function allowing a generic implementation of the Jastrow factor. It is demonstrated with a few examples that the new Jastrow factor achieves 85%-90% of the total correlation energy in a variational quantum Monte Carlo calculation and more than 90% of the diffusion Monte Carlo correlation energy.

Three-particle complexes in two-dimensional semiconductors

We evaluate binding energies of trions X±, excitons bound by a donor or acceptor charge X^{D(A)}, and overcharged acceptors or donors in two-dimensional atomic crystals by mapping the three-body problem in two dimensions onto one particle in a three-dimensional potential treatable by a purposely developed boundary-matching-matrix method. We find that in monolayers of transition metal dichalcogenides the dissociation energy of X^{±} is typically much larger than that of localized exciton complexes, so that trions are more resilient to heating, despite the fact that their recombination line in optics is less redshifted from the exciton line than the line of X^{D(A)}.

Compression algorithm for multideterminant wave functions

A compression algorithm is introduced for multideterminant wave functions which can greatly reduce the number of determinants that need to be evaluated in quantum Monte Carlo calculations. We have devised an algorithm with three levels of compression, the least costly of which yields excellent results in polynomial time. We demonstrate the usefulness of the compression algorithm for evaluating multideterminant wave functions in quantum Monte Carlo calculations, whose computational cost is reduced by factors of between about 2 and over 25 for the examples studied. We have found evidence of sublinear scaling of quantum Monte Carlo calculations with the number of determinants when the compression algorithm is used.

Molecular Properties by Quantum Monte Carlo: An Investigation on the Role of the Wave Function Ansatz and the Basis Set in the Water Molecule

Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely, the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudopotential, and the basis set for QMC calculations. We also introduce a new method for the computation of forces with finite variance on open systems and a new strategy for the definition of the atomic orbitals involved in the Jastrow-Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.