skip to main content


This content will become publicly available on June 12, 2024

Title: Avoiding symmetry roadblocks and minimizing the measurement overhead of adaptive variational quantum eigensolvers
Quantum simulation of strongly correlated systems is potentially the most feasible useful application of near-term quantum computers. Minimizing quantum computational resources is crucial to achieving this goal. A promising class of algorithms for this purpose consists of variational quantum eigensolvers (VQEs). Among these, problem-tailored versions such as ADAPT-VQE that build variational ansätze step by step from a predefined operator pool perform particularly well in terms of circuit depths and variational parameter counts. However, this improved performance comes at the expense of an additional measurement overhead compared to standard VQEs. Here, we show that this overhead can be reduced to an amount that grows only linearly with the numbernof qubits, instead of quartically as in the original ADAPT-VQE. We do this by proving that operator pools of size2n2can represent any state in Hilbert space if chosen appropriately. We prove that this is the minimal size of such complete pools, discuss their algebraic properties, and present necessary and sufficient conditions for their completeness that allow us to find such pools efficiently. We further show that, if the simulated problem possesses symmetries, then complete pools can fail to yield convergent results, unless the pool is chosen to obey certain symmetry rules. We demonstrate the performance of such symmetry-adapted complete pools by using them in classical simulations of ADAPT-VQE for several strongly correlated molecules. Our findings are relevant for any VQE that uses an ansatz based on Pauli strings.

 
more » « less
Award ID(s):
1839136
NSF-PAR ID:
10475599
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Quantum Journal
Date Published:
Journal Name:
Quantum
Volume:
7
ISSN:
2521-327X
Page Range / eLocation ID:
1040
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Many quantum algorithms are developed to evaluate eigenvalues for Hermitian matrices. However, few practical approach exists for the eigenanalysis of non-Hermintian ones, such as arising from modern power systems. The main difficulty lies in the fact that, as the eigenvector matrix of a general matrix can be non-unitary, solving a general eigenvalue problem is inherently incompatible with existing unitary-gate-based quantum methods. To fill this gap, this paper introduces a Variational Quantum Universal Eigensolver (VQUE), which is deployable on noisy intermediate scale quantum computers. Our new contributions include: (1) The first universal variational quantum algorithm capable of evaluating the eigenvalues of non-Hermitian matrices—Inspired by Schur’s triangularization theory, VQUE unitarizes the eigenvalue problem to a procedure of searching unitary transformation matrices via quantum devices; (2) A Quantum Process Snapshot technique is devised to make VQUE maintain the potential quantum advantage inherited from the original variational quantum eigensolver—With additional$$O(log_{2}{N})$$O(log2N)quantum gates, this method efficiently identifies whether a unitary operator is triangular with respect to a given basis; (3) Successful deployment and validation of VQUE on a real noisy quantum computer, which demonstrates the algorithm’s feasibility. We also undertake a comprehensive parametric study to validate VQUE’s scalability, generality, and performance in realistic applications.

     
    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. Abstract

    The road to computing on quantum devices has been accelerated by the promises that come from using Shor’s algorithm to reduce the complexity of prime factorization. However, this promise hast not yet been realized due to noisy qubits and lack of robust error correction schemes. Here we explore a promising, alternative method for prime factorization that uses well-established techniques from variational imaginary time evolution. We create a Hamiltonian whose ground state encodes the solution to the problem and use variational techniques to evolve a state iteratively towards these prime factors. We show that the number of circuits evaluated in each iteration scales as$$O(n^{5}d)$$O(n5d), wherenis the bit-length of the number to be factorized anddis the depth of the circuit. We use a single layer of entangling gates to factorize 36 numbers represented using 7, 8, and 9-qubit Hamiltonians. We also verify the method’s performance by implementing it on the IBMQ Lima hardware to factorize 55, 65, 77 and 91 which are greater than the largest number (21) to have been factorized on IBMQ hardware.

     
    more » « less
  4. Abstract

    One of the cornerstone effects in spintronics is spin pumping by dynamical magnetization that is steadily precessing (around, for example, thez-axis) with frequencyω0due to absorption of low-power microwaves of frequencyω0under the resonance conditions and in the absence of any applied bias voltage. The two-decades-old ‘standard model’ of this effect, based on the scattering theory of adiabatic quantum pumping, predicts that componentISzof spin current vector(ISx(t),ISy(t),ISz)ω0is time-independent whileISx(t)andISy(t)oscillate harmonically in time with a single frequencyω0whereas pumped charge current is zeroI0in the same adiabaticω0limit. Here we employ more general approaches than the ‘standard model’, namely the time-dependent nonequilibrium Green’s function (NEGF) and the Floquet NEGF, to predict unforeseen features of spin pumping: namely precessing localized magnetic moments within a ferromagnetic metal (FM) or antiferromagnetic metal (AFM), whose conduction electrons are exposed to spin–orbit coupling (SOC) of either intrinsic or proximity origin, will pump both spinISα(t)and chargeI(t) currents. All four of these functions harmonically oscillate in time at both even and odd integer multiplesNω0of the driving frequencyω0. The cutoff order of such high harmonics increases with SOC strength, reachingNmax11in the one-dimensional FM or AFM models chosen for demonstration. A higher cutoffNmax25can be achieved in realistic two-dimensional (2D) FM models defined on a honeycomb lattice, and we provide a prescription of how to realize them using 2D magnets and their heterostructures.

     
    more » « less
  5. Abstract

    The proximity of many strongly correlated superconductors to density-wave or nematic order has led to an extensive search for fingerprints of pairing mediated by dynamical quantum-critical (QC) fluctuations of the corresponding order parameter. Here we study anisotropics-wave superconductivity induced by anisotropic QC dynamical nematic fluctuations. We solve the non-linear gap equation for the pairing gap$$\Delta (\theta ,{\omega }_{m})$$Δ(θ,ωm)and show that its angular dependence strongly varies below$${T}_{{\rm{c}}}$$Tc. We show that this variation is a signature of QC pairing and comes about because there are multiples-wave pairing instabilities with closely spaced transition temperatures$${T}_{{\rm{c}},n}$$Tc,n. Taken alone, each instability would produce a gap$$\Delta (\theta ,{\omega }_{m})$$Δ(θ,ωm)that changes sign$$8n$$8ntimes along the Fermi surface. We show that the equilibrium gap$$\Delta (\theta ,{\omega }_{m})$$Δ(θ,ωm)is a superposition of multiple components that are nonlinearly induced below the actual$${T}_{{\rm{c}}}={T}_{{\rm{c}},0}$$Tc=Tc,0, and get resonantly enhanced at$$T={T}_{{\rm{c}},n}\ <\ {T}_{{\rm{c}}}$$T=Tc,n<Tc. This gives rise to strong temperature variation of the angular dependence of$$\Delta (\theta ,{\omega }_{m})$$Δ(θ,ωm). This variation progressively disappears away from a QC point.

     
    more » « less