Understanding low-energy excitations in fermionic systems is crucial for their characterization. They determine the response of the system to external weak perturbations, its dynamical correlation functions, and provide mechanisms for the emergence of exotic phases of matter. In this work, we study the spin excitation spectra of the 1D Fermi-Hubbard model using a digital quantum processor. Introducing a protocol that is naturally suited for simulation on quantum computers, we recover the retarded spin Green's function from the time evolution of simple observables after a specific quantum quench. We exploit the robustness of the protocol to perturbations of the initial state to minimize the quantum resources required for the initial state preparation, and to allocate the majority of them to a Trotterized time-dynamics simulation. This, combined with the intrinsic resilience to noise of the protocol, allows us to accurately reconstruct the spin excitation spectrum for large instances of the 1D Fermi-Hubbard model without making use of expensive error mitigation techniques, using up to 30 qubits of an IBM Heron r2 device.
Understanding low-energy excitations in fermionic systems is crucial for their characterization. They determine the response of the system to external weak perturbations, its dynamical correlation functions, and provide mechanisms for the emergence of exotic phases of matter. In this work, we study the spin excitation spectra of the 1D Fermi-Hubbard model using a digital quantum processor. Introducing a protocol that is naturally suited for simulation on quantum computers, we recover the retarded spin Green's function from the time evolution of simple observables after a specific quantum quench. We exploit the robustness of the protocol to perturbations of the initial state to minimize the quantum resources required for the initial state preparation, and to allocate the majority of them to a Trotterized time-dynamics simulation. This, combined with the intrinsic resilience to noise of the protocol, allows us to accurately reconstruct the spin excitation spectrum for large instances of the 1D Fermi-Hubbard model without making use of expensive error mitigation techniques, using up to 30 qubits of an IBM Heron r2 device.
We report experimental digital quantum simulation of the one-dimensional Fermi-Hubbard model on a superconducting quantum processor at a scale beyond the reach of exact statevector simulation and challenging for state-of-the-art tensor-network methods. We encode this problem using up to 120 qubits through an efficient mapping that reduces circuit complexity, and we improve accuracy through error suppression to simulate dynamical evolution using up to 90 Trotter steps. From a vacancy defect introduced in the middle of an L=31-site (62-qubit) Néel initial state, we directly observe spin-charge separation to t=9 in natural units and quantitatively extract velocities vc, vs, which match classical simulations across model parameters, and analytics in appropriate regimes. We then extend experiments to L=60 fermions (120 qubits) and evolution times to t=6 using 30 Trotter steps. Quantum-processor outputs agree quantitatively with approximate classical simulations performed using a time-dependent variational principle (TDVP) solver and we observe that increasing the TDVP bond dimension through χ=4096 expands the range of evolution times within which agreement has RMSE ∼1% before the approaches diverge. Owing to the large scale of the simulation and the use of overhead-free error-suppression techniques, for simulated evolution times at the limit of quantum/classical agreement (t≳5 in natural hopping units), the wall-clock runtime of the quantum processor is nearly three orders of magnitude faster than TDVP simulations with χ=4096 and implemented using the ITensor and TeNPy packages. These results establish contemporary digital quantum processors as a versatile, quantitatively accurate, and competitive platform for the study of fermionic many-body dynamics in regimes where leading classical methods can become prohibitively expensive.
The elementary excitation spectrum of a many-body quantum system encodes many key properties, including phenomena as diverse as transport, thermalisation and ground state structure. Excitation spectra of strongly correlated systems are typically encoded in dynamical structure factors, which are demanding to measure experimentally and challenging to compute classically. Here we use quench spectroscopy on a superconducting quantum processor to extract excitation spectra of spin chains of L=101 spins. By tailoring the combination of quench protocol and observable, we selectively access distinct excitation sectors across several phases of the spin-1/2 XXZ chain, resolving free magnons, multi-magnon bound states, and two-spinon continua. Notably, we demonstrate that the protocol does not rely on ground state preparation: in the classically challenging gapless regime, we extract spectra directly from the quench dynamics of easily prepared product states, a procedure that is natural and straightforward on quantum hardware. Our work establishes quench spectroscopy as a fast and flexible probe of many-body excitation spectra on digital quantum hardware, introduces a novel quench protocol that does not require costly state preparation routines, and provides a scalable route towards regimes where classical simulation may become intractable.
Strongly correlated quantum systems give rise to many exotic physical phenomena, including high-temperature superconductivity. Simulating these systems on quantum computers may avoid the prohibitively high computational cost incurred in classical approaches. However, systematic errors and decoherence effects presented in current quantum devices make it difficult to achieve this. Here, we simulate the dynamics of the one-dimensional Fermi-Hubbard model using 16 qubits on a digital superconducting quantum processor. We observe separations in the spreading velocities of charge and spin densities in the highly excited regime, a regime that is beyond the conventional quasiparticle picture. To minimize systematic errors, we introduce an accurate gate calibration procedure that is fast enough to capture temporal drifts of the gate parameters. We also employ a sequence of error-mitigation techniques to reduce decoherence effects and residual systematic errors. These procedures allow us to simulate the time evolution of the model faithfully despite having over 600 two-qubit gates in our circuits. Our experiment charts a path to practical quantum simulation of strongly correlated phenomena using available quantum devices.
Simulating quantum many-body systems is essential for understanding complex phenomena in condensed matter, chemistry, and materials science. Quantum computing offers a way to overcome the limitations of classical methods by operating directly in exponentially large Hilbert spaces. In particular, classical methods often fail to capture the dynamics of strongly correlated systems, which makes quantum algorithms a promising alternative for exploring these challenging regimes. This thesis uses quantum dynamics simulations to probe the behaviour of many-body systems. We first introduce a spectroscopy-based protocol to extract excitation spectra from the time evolution of observables and analyse its performance on near-term quantum devices. We show that the protocol requires only short evolution times and is resilient to errors, facilitating its implementation on real hardware. Building on this, we simulate fermionic excitation spectra using current quantum devices, introducing resource optimisations - including an approximate free-fermion initial state, a unitary quench without ancillae, and reduced Trotter steps - that significantly lower circuit depth. These improvements enable simulations at system sizes competitive with the state of the art, yielding results that match classical benchmarks. While most studies focus on unitary dynamics of isolated systems, open quantum dynamics is increasingly relevant for describing non-equilibrium processes across biology, materials science, and quantum control. We extend our scope to systems coupled to an environment and governed by Lindbladian evolution. By studying simulation complexity through physically meaningful quantities - such as correlation length and mixing time - we explore whether open-system dynamics can be as tractable as their unitary counterparts. These findings provide insights into the computational complexity landscape of open-system simulation and highlight regimes where quantum methods offer a clear advantage.
The famous, yet unsolved, Fermi-Hubbard model for strongly-correlated electronic systems is a prominent target for quantum computers. However, accurately representing the Fermi-Hubbard ground state for large instances may be beyond the reach of near-term quantum hardware. Here we show experimentally that an efficient, low-depth variational quantum algorithm with few parameters can reproduce important qualitative features of medium-size instances of the Fermi-Hubbard model. We address 1 × 8 and 2 × 4 instances on 16 qubits on a superconducting quantum processor, substantially larger than previous work based on less scalable compression techniques, and going beyond the family of 1D Fermi-Hubbard instances, which are solvable classically. Consistent with predictions for the ground state, we observe the onset of the metal-insulator transition and Friedel oscillations in 1D, and antiferromagnetic order in both 1D and 2D. We use a variety of error-mitigation techniques, including symmetries of the Fermi-Hubbard model and a recently developed technique tailored to simulating fermionic systems. We also introduce a new variational optimisation algorithm based on iterative Bayesian updates of a local surrogate model.
The simulation of quantum many-body systems, relevant for quantum chemistry and condensed matter physics, is one of the most promising applications of near-term quantum computers before fault-tolerance. However, since the vast majority of quantum computing technologies are built around qubits and discrete gate-based operations, the translation of the physical problem into this framework is a crucial step. This translation will often be device specific, and a suboptimal implementation will be punished by the exponential compounding of errors on real devices. The importance of an efficient mapping is already revealed for models of spinful fermions in two or three dimensions, which naturally arise when the relevant physics relates to electrons. Using the most direct and well-known mapping, the Jordan-Wigner transformation, leads to a nonlocal representation of local degrees of freedom, and necessities efficient decompositions of nonlocal unitary gates into a sequence of hardware accessible local gates. In this paper, we provide a step-by-step recipe for simulating the paradigmatic two-dimensional Fermi-Hubbard model on a quantum computer using only local operations. To provide the ingredients for such a recipe, we briefly review the plethora of different approaches that have emerged recently but focus on the Derby-Klassen compact fermion mapping to make our discussion concrete. We provide a detailed recipe for an end-to-end simulation including embedding on a physical device, preparing initial states such as ground states, simulation of unitary time evolution, and measurement of observables and spectral functions. We explicitly compute the resource requirements for simulating a global quantum quench and conclude by discussing the challenges and future directions for simulating strongly correlated fermionic matter on quantum computers.