Search

We describe a matrix product state (MPS) extension for the Fermionic Quantum Emulator (FQE) software library. We discuss the theory behind symmetry adapted matrix product states for approximating many-body wavefunctions of spin-1/2 fermions, and we present an open-source, MPS-enabled implementation of the FQE interface (MPS-FQE). The software uses the open-source pyblock3 and block2 libraries for most elementary tensor operations, and it can largely be used as a drop-in replacement for FQE that allows for more efficient, but approximate, emulation of larger fermionic circuits. Finally, we show several applications relevant to both near-term and fault-tolerant quantum algorithms where approximate emulation of larger systems is expected to be useful: characterization of state preparation strategies for quantum phase estimation, the testing of different variational quantum eigensolver Ans\"atze, the numerical evaluation of Trotter errors, and the simulation of general quantum dynamics problems. In all these examples, approximate emulation with MPS-FQE allows us to treat systems that are significantly larger than those accessible with a full statevector emulator., Comment: 11 pages, 6 figures

The simulation of chemistry is among the most promising applications of quantum computing. However, most prior work exploring algorithms for block-encoding, time-evolving, and sampling in the eigenbasis of electronic structure Hamiltonians has either focused on modeling finite-sized systems, or has required a large number of plane wave basis functions. In this work, we extend methods for quantum simulation with Bloch orbitals constructed from symmetry-adapted atom-centered orbitals so that one can model periodic \textit{ab initio} Hamiltonians using only a modest number of basis functions. We focus on adapting existing algorithms based on combining qubitization with tensor factorizations of the Coulomb operator. Significant modifications of those algorithms are required to obtain an asymptotic speedup leveraging translational (or, more broadly, Abelian) symmetries. We implement block encodings using known tensor factorizations and a new Bloch orbital form of tensor hypercontraction. Finally, we estimate the resources required to deploy our algorithms to classically challenging model materials relevant to the chemistry of Lithium Nickel Oxide battery cathodes within the surface code.

In this article we present an algorithm to efficiently evaluate the exchange matrix in periodic systems when Gaussian basis set with pseudopotentials are used. The usual algorithm for evaluating exchange matrix scales cubically with the system size because one has to perform O(N2) fast Fourier transforms (FFT). Here we introduce an algorithm that retains the cubic scaling but reduces the prefactor significantly by eliminating the need to do FFTs during each exchange build. This is accomplished by representing the products of Gaussian basis function using a linear combination of an auxiliary basis the number of which scales linearly with the size of the system. We store the potential due to these auxiliary functions in memory which allows us to obtain the exchange matrix without the need to do FFT, albeit at the cost of additional memory requirement. Although the basic idea of using auxiliary functions is not new, our algorithm is cheaper due to a combination of three ingredients: (a) we use robust Pseudospectral method that allows us to use a relatively small number of auxiliary basis to obtain high accuracy (b) we use occ-RI exchange which eliminates the need to construct the full exchange matrix and (c) we use the (interpolative separable density fitting) ISDF algorithm to construct these auxiliary basis that are used in the robust pseudospectral method. The resulting algorithm is accurate and we note that the error in the final energy decreases exponentially rapidly with the number of auxiliary functions.

Obtaining the free energy of large molecules from quantum mechanical energy functions is a longstanding challenge. We describe a method that allows us to estimate, at the quantum mechanical level, the harmonic contributions to the thermodynamics of molecular systems of large size, with modest cost. Using this approach, we compute the vibrational thermodynamics of a series of diamond nanocrystals, and show that the error per atom decreases with system size in the limit of large systems. We further show that we can obtain the vibrational contributions to the binding free energies of prototypical protein-ligand complexes where the exact computation is too expensive to be practical. Our work raises the possibility of routine quantum mechanical estimates of thermodynamic quantities in complex systems.

We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [{\it J. Chem. Theory Comput.} \textbf{2019}, 15, 6137-6253] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest's theorem) for all one-particle properties, while remaining energy conserving for time-independent Hamiltonians. We present the time-dependent orbital-optimized coupled cluster doubles method (Keldysh-OCCD) in analogy with the formalism for zero-temperature dynamics, extended to finite temperatures through the time-dependent action on the Keldysh contour. To demonstrate the conservation property and understand the numerical performance of the method, we apply it to several problems of non-equilibrium finite-temperature dynamics: a 1D Hubbard model with a time-dependent Peierls phase, laser driving of molecular H$_2$, driven dynamics in warm-dense silicon, and transport in the single impurity Anderson model., Comment: 18 pages, 13 figures

We describe a coupled cluster framework for coupled systems of electrons and phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard-Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations which generally perform well for weak to moderate coupling. Finally, we report progress towards an implementation for {\it ab initio} calculations on solids, and present some preliminary results on finite-size models of diamond. We also report the implementation of electron-phonon coupling matrix elements from crystalline Gaussian type orbitals (cGTO) within the PySCF program package.

We discuss the theory and implementation of the finite temperature coupled cluster singles and doubles (FT-CCSD) method including the equations necessary for an efficient implementation of response properties. Numerical aspects of the method including the truncation of the orbital space and integration of the amplitude equations are tested on some simple systems, and we provide some guidelines for applying the method in practice. The method is then applied to the 1D Hubbard model, the uniform electron gas at warm, dense conditions, and some simple materials. The performance on model systems at high temperatures is encouraging: for the 1-dimensional Hubbard model FT-CCSD provides a qualitatively accurate description of finite-temperature correlation effects even at $U = 8$, and it allows for the computation of systematically improvable exchange-correlation energies of the warm, dense UEG over a wide range of conditions. We highlight the obstacles that remain in using the method for realistic {\it ab initio} calculations on materials.

We describe the ground- and excited-state electronic structure of bulk MnO and NiO, two prototypical correlated electron materials, using coupled cluster theory with single and double excitations (CCSD). As a corollary, this work also reports the first implementation of unrestricted periodic ab initio equation-of motion CCSD. Starting from a Hartree-Fock reference, we find fundamental gaps of 3.46 eV and 4.83 eV for MnO and NiO respectively for the 16 unit supercell, slightly overestimated compared to experiment, although finite-size scaling suggests that the gap is more severely overestimated in the thermodynamic limit. From the character of the correlated electronic bands we find both MnO and NiO to lie in the intermediate Mott/charge-transfer insulator regime, although NiO appears as a charge transfer insulator when only the fundamental gap is considered. While the lowest quasiparticle excitations are of metal 3d and O 2p character in most of the Brillouin zone, near the {\Gamma} point, the lowest conduction band quasiparticles are of s character. Our study supports the potential of coupled cluster theory to provide high level many-body insights into correlated solids.

We leverage the Keldysh formalism to extend our implementation of finite temperature coupled cluster theory [\textit{J. Chem. Theory Comput.} 2018, \textit{14}, 5690-5700] to thermal systems that have been driven out of equilibrium. The resulting Keldysh coupled cluster theory is discussed in detail. We describe the implementation of the equations necessary to perform Keldysh coupled cluster singles and doubles calculations of finite temperature dynamics, and we apply the method to some simple systems including a Hubbard model with a Peierls phase and an {\it ab initio} model of warm-dense silicon subject to an ultrafact XUV pulse.

In this comment we address the preprint of Jha and Hirata (arXiv:1809.10316 [physics.chem-ph]) which claims "Numerical Evidence Falsifying Finite-Temperature Many-Body Perturbation Theory." We agree that finite difference differentiation of the exact grand potential is the correct way to verify the terms in the perturbation expansion. However, it is our suspicion that theoretical errors, uncontrolled numerical errors, and/or errors in their implementation have led Jha and Hirata to incorrectly conclude that finite temperature perturbation theory is incorrect. We show numerical evidence from finite difference differentiation of the grand potential that supports the correctness of finite-temperature perturbation theory.

We present a time-dependent formulation of coupled cluster theory. This theory allows for direct computation of the free energy of quantum systems at finite temperature by imaginary time integration and is closely related to the thermal cluster cumulant theory of Mukherjee and co-workers. Our derivation highlights the connection to perturbation theory and zero-temperature coupled cluster theory. We show explicitly how the finite-temperature coupled cluster singles and doubles amplitude equations can be derived in analogy with the zero-temperature theory and how response properties can be efficiently computed using a variational Lagrangian. We discuss the implementation for realistic systems and showcase the potential utility of the method with calculations of the exchange correlation energy of the uniform electron gas at warm dense matter conditions.