skip to main content


Title: Quantum simulation of electronic structure with a transcorrelated Hamiltonian: improved accuracy with a smaller footprint on the quantum computer
Quantum simulations of electronic structure with a transformed Hamiltonian that includes some electron correlation effects are demonstrated. The transcorrelated Hamiltonian used in this work is efficiently constructed classically, at polynomial cost, by an approximate similarity transformation with an explicitly correlated two-body unitary operator. This Hamiltonian is Hermitian, includes no more than two-particle interactions, and is free of electron–electron singularities. We investigate the effect of such a transformed Hamiltonian on the accuracy and computational cost of quantum simulations by focusing on a widely used solver for the Schrödinger equation, namely the variational quantum eigensolver method, based on the unitary coupled cluster with singles and doubles (q-UCCSD) Ansatz. Nevertheless, the formalism presented here translates straightforwardly to other quantum algorithms for chemistry. Our results demonstrate that a transcorrelated Hamiltonian, paired with extremely compact bases, produces explicitly correlated energies comparable to those from much larger bases. For the chemical species studied here, explicitly correlated energies based on an underlying 6-31G basis had cc-pVTZ quality. The use of the very compact transcorrelated Hamiltonian reduces the number of CNOT gates required to achieve cc-pVTZ quality by up to two orders of magnitude, and the number of qubits by a factor of three.  more » « less
Award ID(s):
1800348 1550456
NSF-PAR ID:
10328082
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
22
Issue:
42
ISSN:
1463-9076
Page Range / eLocation ID:
24270 to 24281
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The global minima of urea and thiourea were characterized along with other low‐lying stationary points. Each structure was optimized with the CCSD(T) method and triple‐ζcorrelation consistent basis sets followed by harmonic vibrational frequency computations. Relative energies evaluated near the complete basis set limit with both canonical and explicitly correlated CCSD(T) techniques reveal several subtle but important details about both systems. These computations resolve a discrepancy by demonstrating that the electronic energy of the C2vsecond‐order saddle point of urea lies at least 1.5 kcal mol−1above the C2global minimum regardless of whether the structures were optimized with MP2, CCSD, or CCSD(T). Additionally, urea effectively has one minimum instead of two because the electronic barrier for inversion at one amino group in the Cslocal minimum vanishes at the CCSD(T) CBS limit. Characterization of both systems with the same ab initio methods and large basis sets conclusively establishes that the electronic barriers to inversion at one or both NH2groups in thiourea are appreciably smaller than in urea. CCSDT(Q)/cc‐pVTZ computations show higher‐order electron correlation effects have little impact on the relative energies and are consistently offset by core correlation effects of opposite sign and comparable magnitude.

     
    more » « less
  2. We present molecular dynamics (MD), polarizability driven MD (α-DMD), and pump–probe simulations of Raman spectra of the protonated nitrogen dimer N4H+, and some of its isotopologues, using the explicitly correlated coupled-cluster singles and doubles with perturbative triples [CCSD(T)]-F12b/aug-cc-pVTZ based potential energy surface in permutationally invariant polynomials (PIPs) of Yu et al. [J. Phys. Chem. A 119, 11623 (2015)] and a corresponding PIP-derived CCSD(T)/aug-cc-pVTZ-tr (N:spd, H:sp) polarizability tensor surface (PTS), the latter reported here for the first time. To represent the PTS in terms of a PIP basis, we utilize a recently described formulation for computing the polarizability using a many-body expansion in the orders of dipole–dipole interactions while generating a training set using a novel approach based on linear regression for potential energy distributions. The MD/α-DMD simulations reveal (i) a strong Raman activity at 260 and 2400 cm−1, corresponding to the symmetric N–N⋯H bend and symmetric N–N stretch modes, respectively; (ii) a very broad spectral region in the 500–2000 cm−1 range, assignable to the parallel N⋯H+⋯N proton transfer overtone; and (iii) the presence of a Fermi-like resonance in the Raman spectrum near 2400 cm−1 between the Σg+ N–N stretch fundamental and the Πu overtone corresponding to perpendicular N⋯H+⋯N proton transfer.

     
    more » « less
  3. ABSTRACT

    This paper presents rate coefficients for transitions between rotational levels of the A-type and E-type nuclear spin modifications of methanol induced by collisions with molecular hydrogen. These rate coefficients are required for an accurate determination of methanol abundance in the interstellar medium, where local thermodynamic equilibrium conditions generally do not apply. Time-independent close-coupling quantum scattering calculations have been employed to calculate cross-sections and rate coefficients for the (de-)excitation of methanol in collisions with para- and ortho-H2. These calculations utilized a potential energy surface (PES) for the interaction of methanol with H2 recently computed by the explicitly correlated CCSD(T)-F12a coupled-cluster method that employed a correlation-consistent aug-cc-pVTZ basis. Rate coefficients for temperatures ranging from 3 to 250 K were calculated for all transitions among the first 76 rotational levels of both A-type and E-type methanol, whose energies are less than or equal to 170 K. These rate coefficients are compared with those by Rabli and Flower who carried out coupled-state calculations using a PES computed by second-order many-body perturbation theory. Simple radiative transfer calculations using the present set of rate coefficients are also reported and compared with such calculations using the rate coefficients previously computed by Rabli and Flower.

     
    more » « less
  4. Abstract

    Quantum simulation of chemical systems is one of the most promising near-term applications of quantum computers. The variational quantum eigensolver, a leading algorithm for molecular simulations on quantum hardware, has a serious limitation in that it typically relies on a pre-selected wavefunction ansatz that results in approximate wavefunctions and energies. Here we present an arbitrarily accurate variational algorithm that, instead of fixing an ansatz upfront, grows it systematically one operator at a time in a way dictated by the molecule being simulated. This generates an ansatz with a small number of parameters, leading to shallow-depth circuits. We present numerical simulations, including for a prototypical strongly correlated molecule, which show that our algorithm performs much better than a unitary coupled cluster approach, in terms of both circuit depth and chemical accuracy. Our results highlight the potential of our adaptive algorithm for exact simulations with present-day and near-term quantum hardware.

     
    more » « less
  5. We explore how to build quantum circuits that compute the lowest energy state corresponding to a given Hamiltonian within a symmetry subspace by explicitly encoding it into the circuit. We create an explicit unitary and a variationally trained unitary that maps any vector output by ansatz A(α→) from a defined subspace to a vector in the symmetry space. The parameters are trained varitionally to minimize the energy, thus keeping the output within the labelled symmetry value. The method was tested for a spin XXZ Hamiltonian using rotation and reflection symmetry and H2 Hamiltonian within Sz=0 subspace using S2 symmetry. We have found the variationally trained unitary gives good results with very low depth circuits and can thus be used to prepare symmetry states within near term quantum computers. 
    more » « less