Tomonori Shirakawa

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.

We perform digital quantum simulations of the noninteracting Su-Schrieffer-Heeger (SSH) model using a parameterized quantum circuit. The circuit comprises two main components: the first prepares the initial state from the product state $|0\rangle^{\otimes L}$, where $L$ is the system size; the second consists of $M$ layers of brick-wall unitaries simulating time evolution. The evolution times, encoded as the rotation angles of quantum gates in the second part, are optimized variationally to minimize the energy. The SSH model exhibits two distinct topological phases, depending on the relative strengths of inter- and intra-cell hopping amplitudes. We investigate the evolution of the energy, entanglement entropy, and mutual information towards topologically trivial and nontrivial ground states. Our results find the follows: (i) When the initial and target ground states belong to the same topological phase, the variational energy decreases exponentially, the entanglement entropy quickly saturates in a system-size-independent manner, and the mutual information remains spatially localized, as the number of layers increases. (ii) When the initial and target ground states belong to different topological phases, the variational energy decreases polynomially, the entanglement entropy initially grows logarithmically before decreasing, and the mutual information spreads ballistically across the entire system, with increasing the number of layers. Furthermore, by calculating the polarization, we identify a topological phase transition occurring at an intermediate circuit layer when the initial and final target states lie in different topological characters. Finally, we experimentally confirm this topological phase transition in an 18-site system using 19 qubits on a trapped-ion quantum computer provided by Quantinuum.