Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2510.26300

Fermionic dynamics on a trapped-ion quantum computer beyond exact classical simulation

Alam, Faisal, Bosse, Jan Lukas, Čepaitė, Ieva, Chapman, Adrian, Clinton, Laura, Crichigno, Marcos, Crosson, Elizabeth, Cubitt, Toby, Derby, Charles, Dowinton, Oliver, Eassa, Norhan, Faehrmann, Paul K., Flammia, Steve, Flynn, Brian, Gambetta, Filippo Maria, García-Patrón, Raúl, Hunter-Gordon, Max, Jones, Glenn, Khedkar, Abhishek, Klassen, Joel, Kreshchuk, Michael, McMullan, Edward Harry, Mineh, Lana, Montanaro, Ashley, Mora, Caterina, Morton, John J. L., Nocera, Alberto, Patel, Dhrumil, Rolph, Pete, Santos, Raul A., Seddon, James R., Sheridan, Evan, Somogyi, Wilfrid, Svensson, Marika, Vaishnav, Niam, Wang, Sabrina Yue, Wright, Gethin, Chertkov, Eli, Dreyer, Henrik, Foss-Feig, Michael

Published Oct 30, 2025• 0 citations
arXivHTMLPDF

Abstract

Simulation of the time-dynamics of fermionic many-body systems has long been predicted to be one of the key applications of quantum computers. Such simulations -- for which classical methods are often inaccurate -- are critical to advancing our knowledge and understanding of quantum chemistry and materials, underpinning a wide range of fields, from biochemistry to clean-energy technologies and chemical synthesis. However, the performance of all previous digital quantum simulations of fermions has been matched by classical methods, and it has thus far remained unclear whether near-term, intermediate-scale quantum hardware could offer any computational advantage in this area. Here, we implement an efficient quantum simulation algorithm on Quantinuum's System Model H2 trapped-ion quantum computer for the time dynamics of a 56-qubit system that is too complex for exact classical simulation. We focus on the periodic spinful 2D Fermi-Hubbard model and present evidence of spin-charge separation, where the elementary electron's charge and spin decouple. In the limited cases where ground truth is available through exact classical simulation, we find that it agrees with the results we obtain from the quantum device. Employing long-range Wilson operators to study deconfinement of the effective gauge field between spinons and the effective potential between charge carriers, we find behaviour that differs from predictions made by classical tensor network methods. Our results herald the use of quantum computing for simulating strongly correlated electronic systems beyond the capacity of classical computing.

DOI
10.48550/arxiv.2510.26300

Related papers

Programmable digital quantum simulation of 2D Fermi-Hubbard dynamics using 72 superconducting qubits

Faisal Alam, Jan Lukas Bosse, Ieva Čepaitė +41 more

Simulating the time-dynamics of quantum many-body systems was the original use of quantum computers proposed by Feynman, motivated by the critical role of quantum interactions between electrons in the properties of materials and molecules. Accurately simulating such systems remains one of the most promising applications of general-purpose digital quantum computers, in which all the parameters of the model can be programmed and any desired physical quantity output. However, performing such simulations on today's quantum computers at a scale beyond the reach of classical methods requires advances in the efficiency of simulation algorithms and error mitigation techniques. Here we demonstrate programmable digital quantum simulation of the dynamics of the 2D Fermi-Hubbard model -- one of the best-known simplified models of electrons in crystalline solids -- at a scale beyond exact classical state-vector simulation. We implement simulations of this model on lattice sizes up to ${6\times 6}$ using 72 qubits on Google's Willow quantum processor, across a range of physical parameters, including different on-site electron-electron interaction strengths and magnetic flux values, and study phenomena including formation of magnetic polarons (charge carriers surrounded by local magnetic polarisation), dynamical symmetry-breaking in stripe-ordered states, attraction of charge carriers on an entangled background state known as a valence bond solid, and the approach to equilibrium through thermalisation. We validate our results against exact calculations in parameter regimes where these are feasible, and compare them to approximate classical simulations performed using tensor network and operator propagation methods. Our results demonstrate that meaningful programmable digital quantum simulation of many-body interacting electron models is now feasible on state-of-the-art quantum hardware.

arXiv:2510.26845•Oct 30, 2025•0 citations

Hybrid Oscillator-Qubit Quantum Processors: Simulating Fermions, Bosons, and Gauge Fields

Eleanor Crane, Kevin C. Smith, Teague Tomesh +9 more

Key to solving many societally relevant problems in materials science, quantum chemistry, and nuclear physics is understanding their fundamental microscopic quantum behavior. All of these problems can be phrased as fermionic (e.g., electrons, quarks) and bosonic particles (e.g., lattice vibrations, photons) interacting with each other. However, solving these systems of strongly interacting fermions and bosons is out of reach for classical computers and also non-error-corrected qubit-only quantum computers. In this paper, we introduce a hybrid qubit and bosonic-oscillator quantum-simulation framework, which is so efficient that it could simulate some of these strongly correlated problems even without error correction. Specifically, we leverage the available universal bosonic gate set with O ( 1 ) operations per gate as a function of the bosonic cutoff. This enables our framework to provide exact decompositions of higher-order fermion-boson interactions, including notoriously challenging terms such as ( 2 + 1 ) D gauge theories. We show how to compile dynamics, ancilla-free partial error detection, nonlocal observable measurements, and ground-state energy estimation using both qubit-boson variational quantum eigensolvers and quantum signal processing. We introduce a “post-qubit quantum advantage” experiment which is out of reach for both near-term qubit-only quantum computers and classical computers. We justify this proposal by estimating bosonic gate fidelities through high–boson-occupation Lindblad simulation and by comparing to qubit-only algorithms. For the latter, we develop a new qubit-based algorithm for simulating a beam splitter with a boson number cutoff N max , which scales as O ˜ ( N max 1 + o ( 1 ) ) , asymptotically outperforming previously known qubit methods. However, to reach the same circuit fidelity as qubit-boson hardware with 99.9% per-gate fidelity, a qubit quantum computer using the best-known method would still need fidelities of at least 99.9999%, likely requiring substantial error correction. Our framework can be used in superconducting hardware, which we focus on, as well as trapped-ion and neutral-atom hardware. This work establishes digital hybrid qubit-boson quantum simulation as a promising and advantageous method for solving the many challenging, strongly correlated, fermion-boson models of relevance to science and technology.

arXiv:2409.03747•Physical Review X•Jul 23, 2026•9 citations

Fast, accurate, high-resolution simulation of large-scale Fermi-Hubbard models on a digital quantum processor

Gavin S. Hartnett, Khadijeh Sona Najafi, Aleksei Khindanov +5 more

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 $v_c$, $v_s$, 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 $\sim 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\gtrsim 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.

quant-pharXiv:2605.04025•May 6, 2026•0 citations

Characterising many-body systems via quantum dynamics simulations

Lucia Vilchez-Estevez

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.

•Jan 1, 2025•0 citations

A universal and efficient hybrid digital-analog fermionic quantum simulator

Hao-Tian Wei, Kaden R. A. Hazzard

We present a universal framework to harness fermionic ultracold atom platforms for quantum simulation, showing how variational algorithms on existing hardware can simulate many-body systems well beyond the hardware's native Hamiltonian. Our analysis provides evidence that one can quantum simulate the ground-state properties of a broad class of gapless target Hamiltonians of local observables in a quantum evolution time that grows polynomially with the inverse relative error, $T\sim O(\mathrm{poly}(1/ε))$ up to logarithmic corrections, offering an exponential speedup over na{ï}ve classical algorithms such as exact diagonalization. We provide numerical evidence and theoretical argument that this holds for energy density, density-density, and spin-spin correlations in three qualitatively distinct models -- the repulsive Hubbard model; a Hubbard model augmented with nearest-neighbor attractive interactions, which introduces the phenomenon of pairing; and the Hofstadter-Hubbard model, which introduces a gauge field and fractional quantum Hall physics. This work demonstrates quantum simulation using current fermionic platforms far beyond the models natively implemented in the hardware.

cond-mat.quant-gasCross-listedarXiv:2606.05517•Jun 5, 2026•0 citations

Experimental Quantum Simulation of Chemical Dynamics

Tomas Navickas, Ryan J. MacDonell, Christophe H. Valahu +10 more

Accurate simulation of dynamic processes in molecules and reactions is among the most challenging problems in quantum chemistry. Quantum computers promise efficient chemical simulation, but the existing quantum algorithms require many logical qubits and gates, placing practical applications beyond existing technology. Here, we carry out the first quantum simulations of chemical dynamics by employing a more hardware-efficient encoding scheme that uses both qubits and bosonic degrees of freedom. Our trapped-ion device accurately simulates the dynamics of nonadiabatic chemical processes, which are among the most difficult problems in computational chemistry because they involve strong coupling between electronic and nuclear motions. We demonstrate the programmability and versatility of our approach by simulating the dynamics of three different molecules, as well as open-system dynamics in the condensed phase, all with the same quantum resources. Our approach requires orders of magnitude fewer resources than equivalent qubit-only quantum simulations, demonstrating the potential of using hybrid encoding schemes to accelerate quantum simulations of complex chemical processes, which could have applications in fields ranging from energy conversion and storage to biology and drug design.

arXiv:2409.04044•Journal of the American Chemical Society•May 14, 2025•25 citations
Simulation of the time-dynamics of fermionic many-body systems has long been predicted to be one of the key applications of quantum computers. Such simulations -- for which classical methods are often inaccurate -- are critical to advancing our knowledge and understanding of quantum chemistry and materials, underpinning a wide range of fields, from biochemistry to clean-energy technologies and chemical synthesis. However, the performance of all previous digital quantum simulations of fermions has been matched by classical methods, and it has thus far remained unclear whether near-term, intermediate-scale quantum hardware could offer any computational advantage in this area. Here, we implement an efficient quantum simulation algorithm on Quantinuum's System Model H2 trapped-ion quantum computer for the time dynamics of a 56-qubit system that is too complex for exact classical simulation. We focus on the periodic spinful 2D Fermi-Hubbard model and present evidence of spin-charge separation, where the elementary electron's charge and spin decouple. In the limited cases where ground truth is available through exact classical simulation, we find that it agrees with the results we obtain from the quantum device. Employing long-range Wilson operators to study deconfinement of the effective gauge field between spinons and the effective potential between charge carriers, we find behaviour that differs from predictions made by classical tensor network methods. Our results herald the use of quantum computing for simulating strongly correlated electronic systems beyond the capacity of classical computing.
Simulating the time-dynamics of quantum many-body systems was the original use of quantum computers proposed by Feynman, motivated by the critical role of quantum interactions between electrons in the properties of materials and molecules. Accurately simulating such systems remains one of the most promising applications of general-purpose digital quantum computers, in which all the parameters of the model can be programmed and any desired physical quantity output. However, performing such simulations on today's quantum computers at a scale beyond the reach of classical methods requires advances in the efficiency of simulation algorithms and error mitigation techniques. Here we demonstrate programmable digital quantum simulation of the dynamics of the 2D Fermi-Hubbard model -- one of the best-known simplified models of electrons in crystalline solids -- at a scale beyond exact classical state-vector simulation. We implement simulations of this model on lattice sizes up to 6×6{6\times 6}6×6 using 72 qubits on Google's Willow quantum processor, across a range of physical parameters, including different on-site electron-electron interaction strengths and magnetic flux values, and study phenomena including formation of magnetic polarons (charge carriers surrounded by local magnetic polarisation), dynamical symmetry-breaking in stripe-ordered states, attraction of charge carriers on an entangled background state known as a valence bond solid, and the approach to equilibrium through thermalisation. We validate our results against exact calculations in parameter regimes where these are feasible, and compare them to approximate classical simulations performed using tensor network and operator propagation methods. Our results demonstrate that meaningful programmable digital quantum simulation of many-body interacting electron models is now feasible on state-of-the-art quantum hardware.
Key to solving many societally relevant problems in materials science, quantum chemistry, and nuclear physics is understanding their fundamental microscopic quantum behavior. All of these problems can be phrased as fermionic (e.g., electrons, quarks) and bosonic particles (e.g., lattice vibrations, photons) interacting with each other. However, solving these systems of strongly interacting fermions and bosons is out of reach for classical computers and also non-error-corrected qubit-only quantum computers. In this paper, we introduce a hybrid qubit and bosonic-oscillator quantum-simulation framework, which is so efficient that it could simulate some of these strongly correlated problems even without error correction. Specifically, we leverage the available universal bosonic gate set with O ( 1 ) operations per gate as a function of the bosonic cutoff. This enables our framework to provide exact decompositions of higher-order fermion-boson interactions, including notoriously challenging terms such as ( 2 + 1 ) D gauge theories. We show how to compile dynamics, ancilla-free partial error detection, nonlocal observable measurements, and ground-state energy estimation using both qubit-boson variational quantum eigensolvers and quantum signal processing. We introduce a “post-qubit quantum advantage” experiment which is out of reach for both near-term qubit-only quantum computers and classical computers. We justify this proposal by estimating bosonic gate fidelities through high–boson-occupation Lindblad simulation and by comparing to qubit-only algorithms. For the latter, we develop a new qubit-based algorithm for simulating a beam splitter with a boson number cutoff N max , which scales as O ˜ ( N max 1 + o ( 1 ) ) , asymptotically outperforming previously known qubit methods. However, to reach the same circuit fidelity as qubit-boson hardware with 99.9% per-gate fidelity, a qubit quantum computer using the best-known method would still need fidelities of at least 99.9999%, likely requiring substantial error correction. Our framework can be used in superconducting hardware, which we focus on, as well as trapped-ion and neutral-atom hardware. This work establishes digital hybrid qubit-boson quantum simulation as a promising and advantageous method for solving the many challenging, strongly correlated, fermion-boson models of relevance to science and technology.
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=31L=31L=31-site (62-qubit) Néel initial state, we directly observe spin-charge separation to t=9t=9t=9 in natural units and quantitatively extract velocities vcv_cvc​, vsv_svs​, which match classical simulations across model parameters, and analytics in appropriate regimes. We then extend experiments to L=60L=60L=60 fermions (120 qubits) and evolution times to t=6t=6t=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χ= 4096χ=4096 expands the range of evolution times within which agreement has RMSE ∼1%\sim 1\%∼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≳5t\gtrsim 5t≳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χ= 4096χ=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.
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.
We present a universal framework to harness fermionic ultracold atom platforms for quantum simulation, showing how variational algorithms on existing hardware can simulate many-body systems well beyond the hardware's native Hamiltonian. Our analysis provides evidence that one can quantum simulate the ground-state properties of a broad class of gapless target Hamiltonians of local observables in a quantum evolution time that grows polynomially with the inverse relative error, T∼O(poly(1/ε))T\sim O(\mathrm{poly}(1/ε))T∼O(poly(1/ε)) up to logarithmic corrections, offering an exponential speedup over na{ï}ve classical algorithms such as exact diagonalization. We provide numerical evidence and theoretical argument that this holds for energy density, density-density, and spin-spin correlations in three qualitatively distinct models -- the repulsive Hubbard model; a Hubbard model augmented with nearest-neighbor attractive interactions, which introduces the phenomenon of pairing; and the Hofstadter-Hubbard model, which introduces a gauge field and fractional quantum Hall physics. This work demonstrates quantum simulation using current fermionic platforms far beyond the models natively implemented in the hardware.
Accurate simulation of dynamic processes in molecules and reactions is among the most challenging problems in quantum chemistry. Quantum computers promise efficient chemical simulation, but the existing quantum algorithms require many logical qubits and gates, placing practical applications beyond existing technology. Here, we carry out the first quantum simulations of chemical dynamics by employing a more hardware-efficient encoding scheme that uses both qubits and bosonic degrees of freedom. Our trapped-ion device accurately simulates the dynamics of nonadiabatic chemical processes, which are among the most difficult problems in computational chemistry because they involve strong coupling between electronic and nuclear motions. We demonstrate the programmability and versatility of our approach by simulating the dynamics of three different molecules, as well as open-system dynamics in the condensed phase, all with the same quantum resources. Our approach requires orders of magnitude fewer resources than equivalent qubit-only quantum simulations, demonstrating the potential of using hybrid encoding schemes to accelerate quantum simulations of complex chemical processes, which could have applications in fields ranging from energy conversion and storage to biology and drug design.