Periodically driven interacting quantum many-body systems can exhibit long-lived prethermal dynamics, where local observables retain coherent structure even as entanglement and operator complexity grow. Accessing this regime at the system sizes and times needed to determine physical properties of the prethermal state remains a central challenge: state-of-the-art classical methods become unreliable, while noise in quantum hardware degrades observable expectation values. Here we overcome these limitations for a Floquet Ising magnet realized on a heavy-hex lattice. Using the advanced error mitigation software QESEM on an IBM Heron r3 superconducting quantum processor, we measure magnetization dynamics with percent-level precision and resolve long-lived subharmonic prethermal oscillations in systems of up to 74 qubits. These experiments reach regimes for which leading tensor-network simulations fail to converge, while sparse Pauli-path simulations remain strongly truncation dependent despite extensive computations on advanced GPUs and the Fugaku supercomputer. Leveraging this quantum-accessible regime, we extend finite-size scaling to larger systems and find an unexpectedly slow decrease of the oscillation amplitude with system size, providing strong evidence that this oscillatory response persists in the thermodynamic limit of heavy-hex ladders. A hierarchy of mitigation and validation tests, including unbiased error mitigation, agreement between independent mitigation estimators, noise-model validation on the superconducting hardware, and cross-platform corroboration at selected Floquet cycles on Quantinuum System Model H2 and Quantinuum Helios trapped-ion hardware, supports the reliability of these findings. Our work establishes error-mitigated quantum processors as quantitative scientific instruments for discovering new physics in non-equilibrium quantum matter.
Comments
50 pages, 38 figures, 5 tables, includes appendices
Periodically driven interacting quantum many-body systems can exhibit long-lived prethermal dynamics, where local observables retain coherent structure even as entanglement and operator complexity grow. Accessing this regime at the system sizes and times needed to determine physical properties of the prethermal state remains a central challenge: state-of-the-art classical methods become unreliable, while noise in quantum hardware degrades observable expectation values. Here we overcome these limitations for a Floquet Ising magnet realized on a heavy-hex lattice. Using the advanced error mitigation software QESEM on an IBM Heron r3 superconducting quantum processor, we measure magnetization dynamics with percent-level precision and resolve long-lived subharmonic prethermal oscillations in systems of up to 74 qubits. These experiments reach regimes for which leading tensor-network simulations fail to converge, while sparse Pauli-path simulations remain strongly truncation dependent despite extensive computations on advanced GPUs and the Fugaku supercomputer. Leveraging this quantum-accessible regime, we extend finite-size scaling to larger systems and find an unexpectedly slow decrease of the oscillation amplitude with system size, providing strong evidence that this oscillatory response persists in the thermodynamic limit of heavy-hex ladders. A hierarchy of mitigation and validation tests, including unbiased error mitigation, agreement between independent mitigation estimators, noise-model validation on the superconducting hardware, and cross-platform corroboration at selected Floquet cycles on Quantinuum System Model H2 and Quantinuum Helios trapped-ion hardware, supports the reliability of these findings. Our work establishes error-mitigated quantum processors as quantitative scientific instruments for discovering new physics in non-equilibrium quantum matter.
The simulation of non-equilibrium quantum many-body systems on noisy intermediate-scale quantum (NISQ) processors represents a central challenge and opportunity in modern condensed matter physics. Periodically driven (Floquet) systems can host exotic phases of matter with no equilibrium counterpart, most notably discrete time crystals (DTCs) characterized by spontaneous breaking of discrete time-translation symmetry and robust subharmonic responses. However, the realization of DTCs faces fundamental obstacles: driven many-body systems generically absorb energy and thermalize to an infinite-temperature state, destroying time-crystalline order. Two primary mechanisms have been identified to stabilize DTCs—many-body localization (MBL) through strong disorder, and Floquet prethermalization wherein high-frequency driving suppresses heating over exponentially long timescales. This work investigates the digital quantum simulation of prethermal DTCs on current superconducting quantum processors, with a focus on understanding and mitigating the impact of hardware noise. Recent experiments on 133-qubit heavy-hexagonal devices have demonstrated the realization of clean two-dimensional DTCs and incommensurately modulated DTCs (IM-DTCs) using kicked Ising models, observing period-doubling oscillations persisting for up to 100 Floquet cycles . Crucially, these simulations reveal that quantum noise, often perceived as detrimental, can in certain regimes stabilize novel prethermal DTC phases—a phenomenon identified as noise-induced DTCs supported by symmetry charge pumping or random gate fluctuations . The heavy-hexagonal lattice geometry, compared to square lattices, exhibits significantly slower quantum scrambling, thereby extending the prethermal plateau and enhancing DTC stability. We address two key challenges limiting the scalability of DTC simulations: limited qubit connectivity and environmental noise. By leveraging the universality of native gate sets, we demonstrate how to implement effective couplings among physically disconnected qubits at the cost of increased circuit depth, enabling the study of interaction-range effects on DTC stabilization—a phenomenon previously unobserved experimentally . Furthermore, we benchmark the performance of error mitigation protocols based on depolarising noise models, achieving excellent agreement between error-mitigated hardware data and tensor-network simulations in regimes where classical computations remain tractable . Our results highlight the complementary roles of quantum hardware and classical simulation: quantum processors provide access to regimes where entanglement growth outpaces tensor-network methods, while careful validation against classical benchmarks ensures the reliability of quantum outcomes. This study establishes a concrete pathway for using NISQ devices as experimental platforms for synthesizing and probing novel non-equilibrium quantum phases, advancing both our fundamental understanding of Floquet many-body dynamics and the practical capabilities of quantum simulation. The identification of noise as an active ingredient in phase stabilization suggests new directions for quantum engineering of spatiotemporal orders on near-term hardware.
Many-body localization (MBL) provides a mechanism by which interacting quantum systems evade thermalization, leading to persistent memory of initial conditions and slow entanglement growth. Probing these dynamical signatures in large systems and at long evolution times remains challenging for both classical simulations and current quantum devices. Here we experimentally investigate the ergodic-MBL crossover in quasiperiodic Floquet Ising systems using up to 144 qubits on an IBM Quantum processor. By implementing deep Floquet circuits reaching up to 5000 cycles, we access long-time many-body dynamics beyond the regime explored in previous quantum computing experiments. Measurements of autocorrelation functions reveal a smooth crossover from rapid thermalization at weak quasiperiodic potential strength to persistent correlations in the strong-disorder regime. Notably, in addition to the one-dimensional system, we also observe clear signatures consistent with localization behavior in the two-dimensional system. Furthermore, the quantum Fisher information exhibits logarithmic growth over thousands of Floquet cycles, providing evidence for slow entanglement spreading characteristic of the MBL regime. These results demonstrate that programmable quantum processors can serve as experimental platforms for probing nonergodic quantum many-body dynamics and exploring localization phenomena in regimes beyond the reach of classical simulations.
Floquet many-body phases such as discrete time crystals (DTCs) are typically fragile to imperfections, and stabilizing them on noisy quantum hardware remains a central challenge in nonequilibrium quantum physics. Here, we use IBM Eagle and Heron superconducting processors to implement Floquet dynamics of a kicked Ising model on two-dimensional Kagome lattices, engineered via ancilla-assisted embeddings into the heavy-hex connectivity of the devices. By combining error-mitigated measurements on quantum hardware with matrix-product-state simulations incorporating an ancilla-noise model constructed from experimental device data, we observe long-lived subharmonic magnetization oscillations that are stabilized -- rather than destroyed -- by structured quantum noise. Across different two-dimensional lattice geometries, increasing cases beyond Kagome lattices, and with or without boundary symmetry-charge pumping, ancilla errors effectively act as spatiotemporal disorder that induces stochastic sign flips of the Ising couplings, providing a unified mechanism for robust period-doubling responses. When symmetry-charge pumping is present, intrinsic boundary-localized π modes cooperate with this disorder to yield a boundary-assisted DTC characterized by suppressed scrambling and sharply localized dynamics. In contrast, in implementations without pumping, the noiseless dynamics rapidly thermalize and exhibit no subharmonic order, whereas the same noise process alone generates a DTC-like long-lived subharmonic response over experimentally accessible time windows. These results identify engineered ancilla noise as a practical control knob for inducing, stabilizing, and geometrically tailoring nonequilibrium dynamical order on scalable superconducting quantum processors.
Understanding how isolated quantum many-body systems thermalize remains a central question in modern physics. We study the onset of ergodicity in a two-dimensional disordered Heisenberg Floquet model using digital quantum simulation on IBM's Nighthawk superconducting processor, reaching system sizes of up to 10×10 qubits. We probe ergodicity across different length scales by coarse-graining the system into spatial patches of varying sizes and introducing a measure based on the collision entropy of each patch, enabling a detailed study of when ergodic behavior emerges across scales. The high sampling rate of superconducting quantum processing units, together with an optimal sample estimator, allow us to access patches of sizes up to 3×3. We observe that as the Heisenberg coupling J increases, the noiseless system undergoes a smooth crossover from subergodic to ergodic behavior, with smaller patches approaching their random-matrix-theory values first, thereby revealing a hierarchy across scales. In the region of parameter space where classical tensor-network simulations are reliable, small patches or small values of J, we find excellent agreement with the error-mitigated quantum simulation. Beyond this regime, volume-law entanglement and contraction complexity growth causes the cost of classical methods to rise sharply. Our results open new directions for the use of quantum computers in the study of quantum thermalization.
Time-dependent drives hold the promise of realizing non-equilibrium many-body phenomena that are absent in undriven systems. Yet, drive-induced heating normally destabilizes the systems, which can be parametrically suppressed in the high-frequency regime by using periodic (Floquet) drives. It remains largely unknown to what extent highly controllable quantum simulators can suppress heating in non-periodically driven systems. Using the 78-qubit superconducting quantum processor, Chuang-tzu 2.0, we report the experimental observation of long-lived prethermal phases in many-body systems with tunable heating rates, driven by structured random protocols, characterized by n-multipolar temporal correlations. By measuring both the particle imbalance and subsystem entanglement entropy, we monitor the entire heating process over 1,000 driving cycles and observe the existence of the prethermal plateau. The prethermal lifetime is ‘doubly tunable’: one way by driving frequency, the other by multipolar order; it grows algebraically with the frequency with the universal scaling exponent 2n+1. Using quantum state tomography on different subsystems, we demonstrate a non-uniform spatial entanglement distribution and observe a crossover from area-law to volume-law entanglement scaling. With 78 qubits and 137 couplers in a 2D configuration, the entire far-from-equilibrium heating dynamics are beyond the reach of simulation using tensor-network numerical techniques. Our work highlights superconducting quantum processors as a powerful platform for exploring universal scaling laws and non-equilibrium phases of matter in driven systems in regimes where classical simulation faces formidable challenges.
Time-dependent drives hold promise for realizing non-equilibrium many-body phenomena that are absent in undriven systems1–3. Yet, drive-induced heating normally destabilizes the systems4,5, which can be parametrically suppressed in the high-frequency regime by using periodic (Floquet) drives6,7. It remains largely unknown to what extent highly controllable quantum simulators can suppress heating in non-periodically driven systems. Here, using the 78-qubit superconducting quantum processor, Chuang-tzu 2.0, we report the experimental observation of long-lived prethermal phases in many-body systems with tunable heating rates, driven by structured random protocols, characterized by n-multipolar temporal correlations. By measuring both the particle imbalance and subsystem entanglement entropy, we monitor the entire heating process over 1,000 driving cycles and observe the existence of the prethermal plateau. The prethermal lifetime is ‘doubly tunable’: one way by driving frequency, the other way by multipolar order; it grows algebraically with the frequency with the universal scaling exponent 2n + 1. Using quantum-state tomography on different subsystems, we demonstrate a non-uniform spatial entanglement distribution and observe a crossover from area-law to volume-law entanglement scaling. With 78 qubits and 137 couplers in a two-dimensional configuration, the entire far-from-equilibrium heating dynamics are beyond the reach of simulation using tensor-network numerical techniques. Our work highlights superconducting quantum processors as a powerful platform for exploring universal scaling laws and non-equilibrium phases of matter in driven systems in regimes where classical simulation faces formidable challenges. The existence of a long-lived, prethermal regime in many-body systems with tunable heating rates, driven by structured random protocols, is observed using a 78-qubit superconducting quantum processor.