hardware algorithms simulation

Quantum computer-based simulation of Stark many-body localization in a 1D Fermi-Hubbard model

Curator's Take

AI Commentary

This article demonstrates that a modest‑size superconducting processor can already capture the hallmark crossover between thermalizing and Stark many‑body localized dynamics in a 12‑qubit Fermi‑Hubbard chain, marking one of the first experimental quantum simulations of interacting fermions under a strong linear field. By integrating spin‑resolved Jordan‑Wigner encoding, SWAP networks and a tensor‑network‑driven circuit optimizer, the authors slash two‑qubit gates and depth by roughly 88 %, showing how sophisticated compilation can stretch noisy hardware to tackle non‑trivial many‑body physics. The work not only validates quantum‑hardware predictions against exact numerics but also paves the way for near‑term studies of disorder‑free localization phenomena that are otherwise costly on classical simulators.

— Mark Eatherly

Summary

Many-body localization (MBL) is a dynamical phenomenon that describes the non-ergodicity of isolated quantum many-body systems. In contrast to thermalization, this phenomenon leads to a long-lived memory of initial states of local systems and slow growth of entanglement. In this work, we study Stark MBL in a 12-qubit correlated fermionic system described by the one-dimensional Fermi-Hubbard model using Hamiltonian simulation on an IBM superconducting qubit quantum computer. To enable such a computation on current-day noisy hardware, we combine a series of compilation steps, including the use of the spin-resolved Jordan-Wigner transformation, employing SWAP networks, and integrating a tensor-network-based quantum circuit optimization routine on top of a standard circuit optimization pipeline. As a result, there is approximately an 88$\%$ and 87$\%$ reduction in two-qubit gate count and circuit depth, respectively. Through such simulations of the real-time dynamics using Trotterized quantum circuits, we exhibit a crossover from thermalizing dynamics of the system at a weak tilt of the field to a strongly localized behavior at large tilt with short evolution times. We also benchmark our obtained results with respect to those from exact simulations.