Quantum Simulation of the Static Yukawa Model via Effective Particles on a Trapped-Ion Quantum Simulator
Qingfeng Wang, Christopher G. Yale, Edward C. Tortorici, Brian K. McFarland, Daniel Lobser, Melissa C. Revelle, Andrew J. Landahl, Susan M. Clark, Peter J. Love
quant-phSubmitted Oct 20, 2026• Announced Oct 21, 2026• 0 citations
Trapped-ion, neutral-atom, and circuit QED quantum computers naturally enable coupled qubit-boson simulations. Boson-fermion interactions arise in quantum field theory, where renormalization is required to obtain finite results. We study one of the simplest examples of a coupled fermion-boson system, the static Yukawa model, on the QSCOUT trapped-ion system. We use similarity renormalization to express the same physics in terms of effective particles described by different values of the renormalization parameter $s$. This enables us to map data taken in the window of bosonic occupancy that is experimentally accessible to larger occupancies of distinct effective particles. We use this to study strongly coupled, high-occupancy regimes despite limited readout ranges. Using a trapped-ion qubit-boson device on Sandia National Laboratories' QSCOUT platform, we demonstrate the time evolution of the static Yukawa model.
Trapped-ion, neutral-atom, and circuit QED quantum computers naturally enable coupled qubit-boson simulations. Boson-fermion interactions arise in quantum field theory, where renormalization is required to obtain finite results. We study one of the simplest examples of a coupled fermion-boson system, the static Yukawa model, on the QSCOUT trapped-ion system. We use similarity renormalization to express the same physics in terms of effective particles described by different values of the renormalization parameter s. This enables us to map data taken in the window of bosonic occupancy that is experimentally accessible to larger occupancies of distinct effective particles. We use this to study strongly coupled, high-occupancy regimes despite limited readout ranges. Using a trapped-ion qubit-boson device on Sandia National Laboratories' QSCOUT platform, we demonstrate the time evolution of the static Yukawa model.
Simulating out-of-equilibrium dynamics of quantum field theories in nature is challenging with classical methods, but is a promising application for quantum computers. Unfortunately, simulating interacting bosonic fields involves a high boson-to-qubit encoding overhead. Furthermore, when mapping to qubits, the infinite-dimensional Hilbert space of bosons is necessarily truncated, with truncation errors that grow with energy and time. A qubit-based quantum computer, augmented with an active bosonic register, and with qubit, bosonic, and mixed qubit-boson quantum gates, offers a more powerful platform for simulating bosonic theories. We demonstrate this capability experimentally in a hybrid analog-digital trapped-ion quantum computer, where qubits are encoded in the internal states of the ions, and the bosons in the ions' motional states. Specifically, we simulate nonequilibrium dynamics of a (1+1)-dimensional Yukawa model, a simplified model of interacting nucleons and pions, and measure fermion- and boson-occupation-state probabilities. These dynamics populate high bosonic-field excitations starting from an empty state, and the experimental results capture well such high-occupation states. This simulation approaches the regime where classical methods become challenging, bypasses the need for a large qubit overhead, and removes truncation errors. Our results, therefore, open the way to achieving demonstrable quantum advantage in qubit-boson quantum computing.
Universal digital quantum computers are able to approximate any Hamiltonian evolution as a digital quantum circuit consisting of discrete time steps. Strongly interacting physics models, for example in high-energy physics or materials science, make this a challenging task for near-term devices as they place a high demand on the quantum resources required [1]. One alternative is analog quantum simulation where the native evolution of a quantum system can be tailored to mimic that of the model of interest. While more efficient at specific problems, there are limits on the types of evolution available and on the observables that can be measured. In this conference paper we present first results towards unifying these two concepts in a trapped-ion experiment. Trapped atomic ions are one of the most advanced experimental platforms for both analog quantum simulation and digital quantum computing. Typically, the internal spin degrees of freedom are used qubits or for the simulation of spin models. The harmonic oscillator modes given by the motion of the ions in the trap are only employed transiently to create quantum gates or spin-spin interactions. However, these motional modes constitute a large Hilbert space themselves and present a bosonic degree of freedom, which is costly to model with qubits. Many systems of interest, for example Lattice gauge theories, require both spin and boson degrees of freedom. It was recently shown that they can be encoded more efficiently using the harmonic oscillator modes directly [2].
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 analyze the efficiency of quantum simulations of fermionic and bosonic models in trapped ions. In particular, we study the optimal time of entangling gates and the required number of total elementary gates. Furthermore, we exemplify these estimations in the light of quantum simulations of quantum field theories, condensed-matter physics, and quantum chemistry. Finally, we show that trapped-ion technologies are a suitable platform for implementing quantum simulations involving interacting fermionic and bosonic modes, paving the way for overcoming classical computers in the near future. PACS Codes: 03.67.Ac, 03.67.Lx, 37.10.Ty.
Abstract Elementary quantum mechanics proposes that a closed physical system consistently evolves in a reversible manner. However, control and readout necessitate the coupling of the quantum system to the external environment, subjecting it to relaxation and decoherence. Consequently, system-environment interactions are indispensable for simulating physically significant theories. A broad spectrum of physical systems in condensed-matter and high-energy physics, vibrational spectroscopy, and circuit and cavity QED necessitates the incorporation of bosonic degrees of freedom, such as phonons, photons, and gluons, into optimized fermion algorithms for near-future quantum simulations. In particular, when a quantum system is surrounded by an external environment, its basic physics can usually be simplified to a spin or fermionic system interacting with bosonic modes. Nevertheless, troublesome factors such as the magnitude of the bosonic degrees of freedom typically complicate the direct quantum simulation of these interacting models, necessitating the consideration of a comprehensive plan. This strategy should specifically include a suitable fermion/boson-to-qubit mapping scheme to encode sufficiently large yet manageable bosonic modes, and a method for truncating and/or downfolding the Hamiltonian to the defined subspace for performing an approximate but highly accurate simulation, guided by rigorous error analysis. In this pedagogical tutorial review, we aim to provide such an exhaustive strategy, focusing on encoding and simulating certain bosonic-related model Hamiltonians, inclusive of their static properties and time evolutions. Specifically, we emphasize two aspects: (1) the discussion of recently developed quantum algorithms for these interacting models and the construction of effective Hamiltonians, and (2) a detailed analysis regarding a tightened error bound for truncating the bosonic modes for a class of fermion-boson interacting Hamiltonians.
Gate-based digital quantum simulations offer an exciting new paradigm for studying the many-body physics of strongly correlated systems. In this context, electron-phonon models are challenging for qubit-only quantum simulators, as bosonic degrees of freedom require costly finite-dimensional encodings. Here, we elaborate on an alternative approach based on a digital-analog circuit QED architecture, where fermions are encoded in transmon qubits while bosons are represented directly by microwave resonators. The central building block of this framework is a qubit-resonator Rabi gate that emulates strong electron-phonon coupling and can be implemented through a sequence of resonant Jaynes-Cummings gates interleaved with layers of single-qubit rotations. Using this Rabi gate as the fundamental unitary operation, we construct quantum circuits for the Hubbard-Holstein and Yukawa-Sachdev-Ye-Kitaev models, which describe, respectively, strongly correlated electrons coupled to phonons and phonon-mediated interactions among Majorana fermions. We further demonstrate how nonclassical phonon physics and signatures of quantum chaos in these models can be probed through circuit simulations, and develop measurement and variational protocols tailored to near-term superconducting quantum hardware.