Wave equations provide a natural testbed for near-term quantum simulation of partial differential equations, but hardware demonstrations have remained limited in spatial dimension, equation class, system size, and physically meaningful output. We develop and benchmark structure-preserving, Fourier-based quantum circuits for the one- and two dimensional acoustic wave equations and Dirac dynamics with variable mass on the Quantinuum H2-2 trapped-ion processor. The experiments include one-dimensional grids with up to \(1024\) points and \(32\times32\) two-dimensional grids, corresponding to an encoded state-space dimension of up to \(4096\). Rather than reconstructing the full fields, we estimate subdomain kinetic energies directly from measurement samples. Across all tested acoustic and Dirac dynamics problems, the H2-2 results track the classical kinetic-energy dynamics with mean absolute errors between \(5.9\times10^{-3}\) and \(2.4\times10^{-2}\). At fixed retained bandwidth, the compiled gate counts grow approximately quadratically with the number of grid qubits; the acoustic circuit sizes are essentially independent of evolution time, whereas the cost also grows with the number of product-formula steps. These results provide hardware-level evidence that accurate observable dynamics can remain resolvable for structured wave problems with thousands of encoded degrees of freedom on a present-day trapped-ion processor.
Wave equations provide a natural testbed for near-term quantum simulation of partial differential equations, but hardware demonstrations have remained limited in spatial dimension, equation class, system size, and physically meaningful output. We develop and benchmark structure-preserving, Fourier-based quantum circuits for the one- and two dimensional acoustic wave equations and Dirac dynamics with variable mass on the Quantinuum H2-2 trapped-ion processor. The experiments include one-dimensional grids with up to (1024) points and (32\times32) two-dimensional grids, corresponding to an encoded state-space dimension of up to (4096). Rather than reconstructing the full fields, we estimate subdomain kinetic energies directly from measurement samples. Across all tested acoustic and Dirac dynamics problems, the H2-2 results track the classical kinetic-energy dynamics with mean absolute errors between (5.9\times10^{-3}) and (2.4\times10^{-2}). At fixed retained bandwidth, the compiled gate counts grow approximately quadratically with the number of grid qubits; the acoustic circuit sizes are essentially independent of evolution time, whereas the cost also grows with the number of product-formula steps. These results provide hardware-level evidence that accurate observable dynamics can remain resolvable for structured wave problems with thousands of encoded degrees of freedom on a present-day trapped-ion processor.
We investigate one-dimensional quantum wave-packet dynamics using a grid-based encoding of the wave function onto qubit registers. Time evolution is implemented using a split operator approach, with the kinetic energy operator realized via the Quantum Fourier Transform (QFT) and arbitrary discretized potentials represented by commuting Pauli-Z gates. The algorithm is benchmarked against classical grid-based simulations and validated using a noiseless Qiskit Statevector simulator before being implemented on current quantum hardware (IBM Quantum and IonQ) for qubit counts ranging from two to five. Good agreement with benchmark calculations is obtained for two- and three-qubit implementations, while deviations increase at larger qubit counts due to hardware noise and increased circuit depth. Among the devices examined, the algorithm maintains closer agreement with the benchmark on the trapped-ion IonQ Forte processor as the number of qubits increases.
Explicit quantum circuits for the acoustic and elastic wave equations, obtained by Hamiltonian simulation of their first-order forms, conserve energy and therefore describe closed domains: an outgoing wave reflects off the edge of the computational region or, with periodic boundaries, wraps around and re-enters from the opposite side. We open the domain with the simplest physical absorber, a graded sponge layer in which every wave field is damped at a rate that grows towards the outer edge, and implement the resulting dissipative dynamics on a quantum simulator by Schrödingerisation, which turns a non-unitary linear evolution into a family of unitary ones by means of a single auxiliary variable. We derive the whole chain from the partial differential equation to the circuit, in one, two and three dimensions. An independent finite-difference solver confirms that damping all wave fields absorbs far better than the common velocity-only damping mask. For the sponge the quantum analysis is sharp: recovery is exact on the whole non- negative half of the auxiliary axis at every time; that axis can be truncated periodically at a length set by the largest damping rate and the simulated time; and the success probability is a fixed constant, about 0.4, times the fraction of the energy not yet absorbed, with only one or two auxiliary qubits measured. On the circuit side, the damping step, which generic synthesis compiles into a number of controlled-NOT gates exponential in the number of qubits, factorises exactly into polynomially many multi-controlled phase gates. Gate-level simulations on up to 25 qubits agree with the classical solutions of the damped wave equations.
We design and implement quantum circuits for the simulation of the one-dimensional wave equation on the Quantinuum H1-1 quantum computer. The circuit depth of our approach scales as O ( n 2 ) for n qubits representing the solution on 2 n grid points, and leads to infidelities of O ( 2 − 4 n t 2 ) for simulation time t assuming smooth initial conditions. By varying the qubit count we study the interplay between the algorithmic and physical gate errors to identify the optimal working point of minimum total error. Our approach to simulating the wave equation can be used with appropriate state preparation algorithms across different quantum processors and serve as an application-oriented benchmark. Published by the American Physical Society 2024
The quantum lattice Boltzmann scheme of Succi and Dellar advances a four-component Dirac spinor on a lattice by a fixed sequence of local, exactly norm-preserving operations: a basis rotation, a collision, a streaming shift, and the inverse rotation. While this correspondence has been noted qualitatively, it has remained unclear whether it can be realized exactly and unitarily as a gate-model quantum circuit. This unitarity is a structural property of the scheme, which should, in principle, enable an operator-by-operator port onto quantum hardware. We give a gate-level construction of every operation of the three-dimensional Dirac QLB scheme: the fixed rotation gates, the collision gate, the streaming shift as a controlled increment on a position register, the position-dependent potential as a phase oracle, and periodic and reflecting boundary conditions as unitary circuits. Together, these constructions provide the first complete demonstration that the scheme is exactly representable in the gate model without any approximation beyond the classical method itself. We then compose them into single-axis, two- and three-dimensional time steps. On a state-vector emulator the resulting circuits reproduce the classical QLB solver to machine precision, thus verifying our quantum implementation. Our aim is representability and correctness rather than computational speedup; questions of state preparation, asymptotic cost, and hardware noise are identified as open avenues for future work. Appendices restore physical units, derive the scheme's Zitterbewegung in closed form and compare it with a trapped-ion simulation, and execute a reduced two-site instance on a Rigetti superconducting processor, where choosing the qubits from the live calibration lowers the state-preparation error. All operators, circuits, and figures are reproducible from the open-source quantumKineticMethods library.
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].
Measurement-based quantum simulation (MBQS)---a recently proposed architecture for simulating lattice gauge theories---implements Hamiltonian dynamics by consuming a model-specific entangled resource state with adaptive mid-circuit measurements, rather than by a gate-based circuit. The local constraints in lattice gauge theories are mirrored by the higher-form symmetries of the resource state. Here we report, to our knowledge, the first experimental realization of MBQS of real-time dynamics in the (2+1)-dimensional Z2 gauge theory using the Quantinuum System Model H2 trapped-ion processor. We observe coherent evolution of gauge-invariant observables on 2×2 and 3×3 spatial lattices, consuming virtual three-dimensional cluster states of 200 and 288 resource-state qubits that are generated from instantaneous blocks of 48 and 54 qubits within the 56-qubit register by measurement, reset, and re-entanglement. The measurement record that drives the evolution simultaneously provides one-form-symmetry syndromes at no additional cost, enabling postselection that strongly suppresses observed Gauss-law violations and improves aggregate agreement with ideal Trotterized dynamics. Our results demonstrate that MBQS is a viable, symmetry-aware architecture for simulating lattice field theories on present-day hardware.