Abstract Level spectroscopy stands as a powerful method for identifying the transition point that delineates distinct quantum phases. Since each quantum phase exhibits a characteristic sequence of excited states, the crossing of energy levels between low-lying excited states offers a reliable mean to estimate the phase transition point. While approaches like the Variational Quantum Eigensolver are useful for approximating ground states of interacting systems using quantum computing, capturing low-energy excitations remains challenging. In our study, we introduce an equivariant quantum circuit that preserves the total spin and the translational symmetry to accurately describe singlet and triplet excited states in the J 1 – J 2 Heisenberg model on a chain, which are crucial for characterizing its transition point. Additionally, we assess the impact of noise on the variational state, showing that conventional mitigation techniques like Zero Noise Extrapolation reliably restore its physical properties.
Abstract Level spectroscopy stands as a powerful method for identifying the transition point that delineates distinct quantum phases. Since each quantum phase exhibits a characteristic sequence of excited states, the crossing of energy levels between low-lying excited states offers a reliable mean to estimate the phase transition point. While approaches like the Variational Quantum Eigensolver are useful for approximating ground states of interacting systems using quantum computing, capturing low-energy excitations remains challenging. In our study, we introduce an equivariant quantum circuit that preserves the total spin and the translational symmetry to accurately describe singlet and triplet excited states in the J 1 – J 2 Heisenberg model on a chain, which are crucial for characterizing its transition point. Additionally, we assess the impact of noise on the variational state, showing that conventional mitigation techniques like Zero Noise Extrapolation reliably restore its physical properties.
Quantum computing allows for the manipulation of highly correlated states whose properties quickly go beyond the capacity of any classical method to calculate. Thus one natural problem which could lend itself to quantum advantage is the study of ground-states of condensed matter models, and the transitions between them. However, current levels of hardware noise can require extensive application of error-mitigation techniques to achieve reliable computations. In this work, we use several IBM devices to explore a finite-size spin model with multiple “phaselike” regions characterized by distinct ground-state configurations. Using preoptimized Variational Quantum Eigensolver (VQE) solutions, we demonstrate that in contrast to calculating the energy, where zero-noise extrapolation is required in order to obtain qualitatively accurate yet still unreliable results, calculations of the energy derivative, two-site spin correlation functions, and the fidelity susceptibility yield accurate behavior across multiple regions, even with minimal or no application of error-mitigation approaches. Taken together, these sets of observables could be used to identify level crossings in a simple, noise-robust manner which is agnostic to the method of ground state preparation. This work shows promising potential for near-term application to identifying quantum phase transitions, including avoided crossings and nonadiabatic conical intersections in electronic structure calculations. Published by the American Physical Society 2024
Resolving Spectral Gaps of the Frustrated Spin-1/2 J1-J2 Chain via Dual Spin-Sector Variational Quantum Eigensolver
Finding excited states and resolving spectral gaps in frustrated quantum spin systems using Variational Quantum Eigensolvers is complicated by state-ordering inversions. In the frustrated spin-1/2 J1-J2 Heisenberg chain, standard single-sector variational quantum deflation collapses onto intermediate triplet excitation, preventing the isolation of excited singlets and the singlet-singlet gap. To resolve these gaps, a Dual Spin-Sector Variational Quantum Eigensolver framework is introduced for near-term quantum architectures. By partitioning the search across magnetization sectors Sz=0 and Sz=1, the lowest triplet is resolved as an unconstrained ground state in Sz=1 with zero deflation penalty. Leveraging exact global spin-rotation symmetry, this state is mapped analytically into Sz=0 without ancilla overhead, unlocking dual deflation via shallow compute-uncompute circuits. Tracking both singlet-triplet and singlet-singlet gaps across coupling ratios J2/J1 from 0 to 1 for chains up to 16 qubits faithfully captures all critical phenomena: the continuous Berezinskii-Kosterlitz-Thouless transition, exact Majumdar-Ghosh dimer degeneracy, level-ordering inversion in the gapless liquid phase, and incommensurate spiral excitations. Circuit depth scaling up to twenty layers provides an overparameterization mechanism overcoming classical local minima in strongly frustrated regimes, reducing relative energy errors below 1% and elevating state fidelities above 97% on 16 qubits. Projective measurement sampling demonstrates resilience against shot noise, providing a practical blueprint for exploring frustrated quantum magnetism on near-term quantum devices.
Variational quantum algorithms are among the most promising approaches for simulating interacting quantum many-body systems on noisy intermediate-scale quantum (NISQ) devices. However, the practical success of variational quantum eigensolvers (VQE) critically depends on the structure of the chosen variational ansatz. In this work, we investigate the ground-state properties of the one-dimensional antiferromagnetic Heisenberg spin-1/2 chain using both generic hardware-efficient ansatz and physics-informed, symmetry-preserving variational circuits. We benchmark variational results against exact diagonalization and noiseless simulations, and subsequently validate the approach on real IQM Garnet quantum hardware. Our results demonstrate that incorporating physical symmetries into the circuit design leads to significantly improved energy estimates, enhanced robustness against hardware noise, and clearer convergence behavior when compared to hardware-efficient ansatz under identical resource constraints. These findings highlight the importance of problem specific ansatz construction for reliable quantum simulations in the NISQ era.
Mapping quantum phase diagrams through classical simulation is notoriously resource-intensive, as even small systems far from the thermodynamic limit demand prohibitive computational effort. The variational quantum eigensolver (VQE) offers a compelling alternative, exploiting approximate ground states to distinguish phases. An appealing proposal, dubbed as Delta-VQE, determines critical points by contrasting variational energies optimized from reference states of distinct phases. Intriguingly, the diagnostic sharpens as circuit depth decreases, highlighting its promise as a resource-conscious probe of quantum criticality. To probe the broader applicability and underlying mechanisms of this approach, we investigate the one-dimensional transverse-field Ising model with a three-spin cluster interaction, a setting in which the Ising transitions are generally situated beyond the self-dual line. We demonstrate that, whenever dual ansätze are employed, Delta-VQE invariably detects the self-dual points rather than the true criticality. In contrast, when ansätze are carefully tailored to embody the competing phases across the boundary, the genuine Ising critical point can be successfully identified with only minor finite-size effects. Our results establish that, while Delta-VQE provides a resource-efficient probe of quantum criticality without requiring precise ground-state preparation, its diagnostic power is fundamentally contingent upon the judicious selection of physically representative ansätze.
Abstract Mapping out phase diagrams of quantum systems using classical simulations can be challenging or intractable due to the computational resources required to simulate even small quantum systems far away from the thermodynamic limit. We investigate using quantum computers and the variational quantum eigensolver (VQE) for this task. In contrast to the task of preparing the exact ground state using VQE, sketching phase diagrams might require less quantum resources and accuracy, because low fidelity approximations to the ground state may be enough to correctly identify different phases. We used classical numerical simulations of low-depth VQE circuits to compute order parameters for four well-studied spin and fermion models which represent a mix of 1D and 2D, and exactly-solvable and classically hard systems. We find that it is possible to predict the location of phase transitions up to reasonable accuracy using states produced by VQE even when their overlap with the true ground state is small. Further, we introduce a model-agnostic predictor of phase transitions based on the speed with which the VQE energy improves with respect to the circuit depth, and find that in some cases this is also able to predict phase transitions.
We propose a momentum-space based variational quantum eigensolver (VQE) framework for simulating quasiparticle excitations in interacting quantum many-body systems on near-term quantum devices. Leveraging translational invariance and other symmetries of the Hamiltonian, we reconstruct the momentum-resolved quasiparticle excitation spectrum through targeted simulation of low-lying excited states using VQE. We construct a translationally symmetric variational ansatz designed to evolve a free-fermion particle-hole excited state with definite momentum q to an excited state of the interacting system at the same momentum, employing a fermionic fast Fourier transform (FFFT) circuit coupled to a Hamiltonian Variational Ansatz (HVA) circuit. Even though the particle number is not explicitly conserved in the variational ansatz, the correct quasiparticle state is reached by energetic optimization. We benchmark the performance of the proposed VQE implementation on the XXZ Hamiltonian, which maps onto the Tomonaga-Luttinger liquid in the fermionic representation. Our numerical results show that VQE can capture the low-lying excitation spectrum of the bosonic quasiparticle/two-spinon dispersion of this model at various interaction strengths. We estimate the renormalized velocity of the quasiparticles by calculating the slope of the dispersion near zero momentum using the VQE-optimized energies at different system sizes, and demonstrate that it closely matches theoretical results obtained from Bethe ansatz. Finally, we highlight extensions of our proposed VQE implementation to simulate quasiparticles in other interacting quantum many-body systems.