Time Propagation and Spectroscopy of Fermionic Systems Using a Stochastic Technique.

We present a stochastic method for solving the time-dependent Schrödinger equation, generalizing a ground state full configuration interaction quantum Monte Carlo method. By performing the time integration in the complex plane close to the real-time axis, the numerical effort is kept manageable and the analytic continuation to real frequencies is efficient. This allows us to perform ab initio calculation of electron spectra for strongly correlated systems. The method can be used as a cluster solver for embedding schemes

Orders of magnitude increased accuracy for quantum many-body problems on quantum computers via an exact transcorrelated method

Transcorrelated methods provide an efficient way of partially transferring the description of electronic correlations from the ground-state wave function directly into the underlying Hamiltonian. In particular, Dobrautz et al. [Phys. Rev. B 99, 075119 (2019)2469-995010.1103/PhysRevB.99.075119] have demonstrated that the use of momentum-space representation, combined with a nonunitary similarity transformation, results in a Hubbard Hamiltonian that possesses a significantly more "compact"ground-state wave function, dominated by a single Slater determinant. This compactness/single-reference character greatly facilitates electronic structure calculations. As a consequence, however, the Hamiltonian becomes non-Hermitian, posing problems for quantum algorithms based on the variational principle. We overcome these limitations with the Ansatz-based quantum imaginary-time evolution algorithm and apply the transcorrelated method in the context of digital quantum computing. We demonstrate that this approach enables up to four orders of magnitude more accurate and compact solutions in various instances of the Hubbard model at intermediate interaction strength (U/t=4), enabling the use of shallower quantum circuits for wave-function Ans\ue4tzes. In addition, we propose a more efficient implementation of the quantum imaginary-time evolution algorithm in quantum circuits that is tailored to non-Hermitian problems. To validate our approach, we perform hardware experiments on the ibmq_lima quantum computer. Our work paves the way for the use of exact transcorrelated methods for the simulations of ab initio systems on quantum computers

Spin-Pure Stochastic-CASSCF via GUGA-FCIQMC Applied to Iron-Sulfur Clusters.

Funder: Max-Planck-GesellschaftIn this work, we demonstrate how to efficiently compute the one- and two-body reduced density matrices within the spin-adapted full configuration interaction quantum Monte Carlo (FCIQMC) method, which is based on the graphical unitary group approach (GUGA). This allows us to use GUGA-FCIQMC as a spin-pure configuration interaction (CI) eigensolver within the complete active space self-consistent field (CASSCF) procedure and hence to stochastically treat active spaces far larger than conventional CI solvers while variationally relaxing orbitals for specific spin-pure states. We apply the method to investigate the spin ladder in iron-sulfur dimer and tetramer model systems. We demonstrate the importance of the orbital relaxation by comparing the Heisenberg model magnetic coupling parameters from the CASSCF procedure to those from a CI-only (CASCI) procedure based on restricted open-shell Hartree-Fock orbitals. We show that the orbital relaxation differentially stabilizes the lower-spin states, thus enlarging the coupling parameters with respect to the values predicted by ignoring orbital relaxation effects. Moreover, we find that, while CASCI results are well fit by a simple bilinear Heisenberg Hamiltonian, the CASSCF eigenvalues exhibit deviations that necessitate the inclusion of biquadratic terms in the model Hamiltonian

Ferromagnetic domains in the large- U Hubbard model with a few holes: A full configuration interaction quantum Monte Carlo study

Two-dimensional Hubbard lattices with two or three holes are investigated as a function of U in the large-U limit. In the so-called Nagaoka limit (one-hole system at infinite U), it is known that the Hubbard model exhibits a ferromagnetic ground state. Here, by means of exact full configuration interaction quantum Monte Carlo simulations applied to periodic lattices up to 24 sites, we compute spin-spin correlation functions as a function of increasing U. The correlation functions clearly demonstrate the onset of ferromagnetic domains, centered on individual holes. The overall total spin of the wave functions remains the lowest possible (0 or 12, depending on the number of holes). The ferromagnetic domains appear at interaction strengths comparable to the critical interaction strengths of the Nagaoka transition in finite systems with strictly one hole. The existence of such ferromagnetic domains is the signature of Nagaoka physics in Hubbard systems with a small (but greater than 1) number of holes

Combined unitary and symmetric group approach applied to low-dimensional Heisenberg spin systems

A novel combined unitary and symmetric group approach is used to study the spin-1/2 Heisenberg model and related Fermionic systems in a total spin-adapted representation, using a linearly-parameterised Ansatz for the many-body wave function. We show that a more compact ground-state wave function representation-indicated by a larger leading ground-state coefficient-is obtained when combining the symmetric group S-n, in the form of permutations of the underlying lattice site ordering, with the cumulative spin coupling based on the unitary group U(n). In one-dimensional systems the observed compression of the wave function is reminiscent of block-spin renormalization group approaches, and allows us to study larger lattices (here taken up to 80 sites) with the spin-adapted full configuration interaction quantum Monte Carlo method, which benefits from the sparsity of the Hamiltonian matrix and the corresponding sampled eigenstates that emerge from the reordering. We find that in an optimal lattice ordering the configuration state function with highest weight already captures with high accuracy the spin-spin correlation function of the exact ground-state wave function. This feature is found for more general lattice models, such as the Hubbard model, and ab initio quantum chemical models, exemplified by one-dimensional hydrogen chains. We also provide numerical evidence that the optimal lattice ordering for the unitary group approach is not generally equivalent to the optimal ordering obtained for methods based on matrix-product states, such as the density-matrix renormalization group approach

Ab Initio Transcorrelated Method enabling accurate Quantum Chemistry on near-term Quantum Hardware

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 limited connectivity) hamper the straightforward implementation of most quantum algorithms and call for more noise-resilient solutions. In quantum chemistry, the limited number of available qubits and gate operations is particularly restrictive since, for each molecular orbital, one needs, in general, two qubits. In this study, we propose an explicitly correlated Ansatz based on the transcorrelated (TC) approach, which transfers -- without any approximation -- correlation from the wavefunction directly into the Hamiltonian, thus reducing the number of resources needed to achieve accurate results with noisy, near-term quantum devices. In particular, we show that the exact transcorrelated approach not only allows for more shallow circuits but also improves the convergence towards the so-called basis set limit, providing energies within chemical accuracy to experiment with smaller basis sets and, therefore, 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 just 4 and 6 qubits, respectively. Conventional methods require at least ten times more qubits for the same accuracy

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 cost, 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.Comment: Submitted to J Chem Phy

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