skip to main content


Title: Nearly tight Trotterization of interacting electrons
We consider simulating quantum systems on digital quantum computers. We show that the performance of quantum simulation can be improved by simultaneously exploiting commutativity of the target Hamiltonian, sparsity of interactions, and prior knowledge of the initial state. We achieve this using Trotterization for a class of interacting electrons that encompasses various physical systems, including the plane-wave-basis electronic structure and the Fermi-Hubbard model. We estimate the simulation error by taking the transition amplitude of nested commutators of the Hamiltonian terms within the η -electron manifold. We develop multiple techniques for bounding the transition amplitude and expectation of general fermionic operators, which may be of independent interest. We show that it suffices to use ( n 5 / 3 η 2 / 3 + n 4 / 3 η 2 / 3 ) n o ( 1 ) gates to simulate electronic structure in the plane-wave basis with n spin orbitals and η electrons, improving the best previous result in second quantization up to a negligible factor while outperforming the first-quantized simulation when n = η 2 − o ( 1 ) . We also obtain an improvement for simulating the Fermi-Hubbard model. We construct concrete examples for which our bounds are almost saturated, giving a nearly tight Trotterization of interacting electrons.  more » « less
Award ID(s):
1839204
NSF-PAR ID:
10298206
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Quantum
Volume:
5
ISSN:
2521-327X
Page Range / eLocation ID:
495
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The Hubbard model is an essential tool for understanding many-body physics in condensed matter systems. Artificial lattices of dopants in silicon are a promising method for the analog quantum simulation of extended Fermi-Hubbard Hamiltonians in the strong interaction regime. However, complex atom-based device fabrication requirements have meant emulating a tunable two-dimensional Fermi-Hubbard Hamiltonian in silicon has not been achieved. Here, we fabricate 3 × 3 arrays of single/few-dopant quantum dots with finite disorder and demonstrate tuning of the electron ensemble using gates and probe the many-body states using quantum transport measurements. By controlling the lattice constants, we tune the hopping amplitude and long-range interactions and observe the finite-size analogue of a transition from metallic to Mott insulating behavior. We simulate thermally activated hopping and Hubbard band formation using increased temperatures. As atomically precise fabrication continues to improve, these results enable a new class of engineered artificial lattices to simulate interactive fermionic models. 
    more » « less
  2. Abstract

    The quantum simulation of quantum chemistry is a promising application of quantum computers. However, forNmolecular orbitals, the$${\mathcal{O}}({N}^{4})$$O(N4)gate complexity of performing Hamiltonian and unitary Coupled Cluster Trotter steps makes simulation based on such primitives challenging. We substantially reduce the gate complexity of such primitives through a two-step low-rank factorization of the Hamiltonian and cluster operator, accompanied by truncation of small terms. Using truncations that incur errors below chemical accuracy allow one to perform Trotter steps of the arbitrary basis electronic structure Hamiltonian with$${\mathcal{O}}({N}^{3})$$O(N3)gate complexity in small simulations, which reduces to$${\mathcal{O}}({N}^{2})$$O(N2)gate complexity in the asymptotic regime; and unitary Coupled Cluster Trotter steps with$${\mathcal{O}}({N}^{3})$$O(N3)gate complexity as a function of increasing basis size for a given molecule. In the case of the Hamiltonian Trotter step, these circuits have$${\mathcal{O}}({N}^{2})$$O(N2)depth on a linearly connected array, an improvement over the$${\mathcal{O}}({N}^{3})$$O(N3)scaling assuming no truncation. As a practical example, we show that a chemically accurate Hamiltonian Trotter step for a 50 qubit molecular simulation can be carried out in the molecular orbital basis with as few as 4000 layers of parallel nearest-neighbor two-qubit gates, consisting of fewer than 105non-Clifford rotations. We also apply our algorithm to iron–sulfur clusters relevant for elucidating the mode of action of metalloenzymes.

     
    more » « less
  3. The Hartree–Fock–Bogoliubov (HFB) theory is the starting point for treating superconducting systems. However, the computational cost for solving large scale HFB equations can be much larger than that of the Hartree–Fock equations, particularly when the Hamiltonian matrix is sparse, and the number of electrons N is relatively small compared to the matrix size N b . We first provide a concise and relatively self-contained review of the HFB theory for general finite sized quantum systems, with special focus on the treatment of spin symmetries from a linear algebra perspective. We then demonstrate that the pole expansion and selected inversion (PEXSI) method can be particularly well suited for solving large scale HFB equations. For a Hubbard-type Hamiltonian, the cost of PEXSI is at most 𝒪( N b 2 ) for both gapped and gapless systems, which can be significantly faster than the standard cubic scaling diagonalization methods. We show that PEXSI can solve a two-dimensional Hubbard-Hofstadter model with N b up to 2.88 × 10 6 , and the wall clock time is less than 100 s using 17 280 CPU cores. This enables the simulation of physical systems under experimentally realizable magnetic fields, which cannot be otherwise simulated with smaller systems. 
    more » « less
  4. Understanding H 2 binding and activation is important in the context of designing transition metal catalysts for many processes, including hydrogenation and the interconversion of H 2 with protons and electrons. This work reports the first thermodynamic and kinetic H 2 binding studies for an isostructural series of first-row metal complexes: NiML, where M = Al ( 1 ), Ga ( 2 ), and In ( 3 ), and L = [N( o -(NCH 2 P i Pr 2 )C 6 H 4 ) 3 ] 3− . Thermodynamic free energies (Δ G °) and free energies of activation (Δ G ‡ ) for binding equilibria were obtained via variable-temperature 31 P NMR studies and lineshape analysis. The supporting metal exerts a large influence on the thermodynamic favorability of both H 2 and N 2 binding to Ni, with Δ G ° values for H 2 binding found to span nearly the entire range of previous reports. The non-classical H 2 adduct, (η 2 -H 2 )NiInL ( 3 -H 2 ), was structurally characterized by single-crystal neutron diffraction—the first such study for a Ni(η 2 -H 2 ) complex or any d 10 M(η 2 -H 2 ) complex. UV-Vis studies and TD-DFT calculations identified specific electronic structure perturbations of the supporting metal which poise NiML complexes for small-molecule binding. ETS-NOCV calculations indicate that H 2 binding primarily occurs via H–H σ-donation to the Ni 4p z -based LUMO, which is proposed to become energetically accessible as the Ni(0)→M( iii ) dative interaction increases for the larger M( iii ) ions. Linear free-energy relationships are discussed, with the activation barrier for H 2 binding (Δ G ‡ ) found to decrease proportionally for more thermodynamically favorable equilibria. The Δ G ° values for H 2 and N 2 binding to NiML complexes were also found to be more exergonic for the larger M( iii ) ions. 
    more » « less
  5. Because Fermi liquids are inherently non-interacting states of matter, all electronic levels below the chemical potential are doubly occupied. Consequently, the simplest way of breaking the Fermi-liquid theory is to engineer a model in which some of those states are singly occupied, keeping time-reversal invariance intact. We show that breaking an overlooked1 local-in-momentum space ℤ2 symmetry of a Fermi liquid does precisely this. As a result, although the Mott transition from a Fermi liquid is correctly believed to arise without breaking any continuous symmetry, a discrete symmetry is broken. This symmetry breaking serves as an organizing principle for Mott physics whether it arises from the tractable Hatsugai–Kohmoto model or the intractable Hubbard model. Through a renormalization-group analysis, we establish that both are controlled by the same fixed point. An experimental manifestation of this fixed point is the onset of particle–hole asymmetry, a widely observed2,3,4,5,6,7,8,9,10 phenomenon in strongly correlated systems. Theoretically, the singly occupied region of the spectrum gives rise to a surface of zeros of the single-particle Green function, denoted as the Luttinger surface. Using K-homology, we show that the Bott topological invariant guarantees the stability of this surface to local perturbations. Our proof demonstrates that the strongly coupled fixed point only corresponds to those Luttinger surfaces with co-dimension p + 1 with odd p. We conclude that both Hubbard and Hatsugai–Kohmoto models lie in the same high-temperature universality class and are controlled by a quartic fixed point with broken ℤ2 symmetry. 
    more » « less