We simulate the dynamics of a spin-1/2 chain with nearest neighbor Ising interactions, quenched disorder and periodic driving over 57 qubits on a current quantum computer. A transition between time crystalline behaviour and non-symmetry breaking dynamics is observed via critical fluctuations in the sub-harmonic frequency response of the system. We probe shorter time scales by implementing the Floquet unitary directly, as well as longer time scales by making use of an equivalent Hamiltonian description. In order to extract the signal from the noisy data produced by current quantum computer devices, we develop a strategy for error mitigation and show that the results are robust under variations of the parameters introduced in this scheme. Our findings are in agreement with previous numerical simulations but add larger system size data points to the phase diagram.
We simulate the dynamics of a spin-1/2 chain with nearest neighbor Ising interactions, quenched disorder and periodic driving over 57 qubits on a current quantum computer. A transition between time crystalline behaviour and non-symmetry breaking dynamics is observed via critical fluctuations in the sub-harmonic frequency response of the system. We probe shorter time scales by implementing the Floquet unitary directly, as well as longer time scales by making use of an equivalent Hamiltonian description. In order to extract the signal from the noisy data produced by current quantum computer devices, we develop a strategy for error mitigation and show that the results are robust under variations of the parameters introduced in this scheme. Our findings are in agreement with previous numerical simulations but add larger system size data points to the phase diagram.
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.
Novel dynamical phases that violate ergodicity have been a subject of extensive research in recent years. A periodically driven system is naively expected to lose all memory of its initial state due to thermalization, yet this can be avoided in the presence of many-body localization. A discrete time crystal represents a driven system whose local observables spontaneously break time translation symmetry and retain memory of the initial state indefinitely. Here we report the observation of a discrete time crystal on a chain consisting of 57 superconducting qubits on a state--of--the--art quantum computer. We probe random initial states and compare the cases of vanishing and finite disorder to distinguish many-body localization from pre-thermal dynamics. We further report results on the dynamical phase transition between the discrete time crystal and a thermal regime, which is observed via critical fluctuations in the system's sub-harmonic frequency response and a significant speed-up of spin depolarisation.
A discrete time crystal is an out-of-equilibrium phase of matter characterized by the spontaneous breaking of discrete time-translation symmetry. Using extensive numerical simulations based on matrix-product-state methods, we provide evidence for discrete time-crystalline behavior in strongly disordered spin chains with Heisenberg interactions, including the isotropic point, subject to periodic driving. Starting from a many-body localized regime, we observe that rotations induced by delta kicks produce a pronounced subharmonic response at half the drive frequency in spin observables. We investigate the stability of this response against rotation-angle errors through entanglement entropy, quantum Fisher information, short-range spin correlations, and restricted-control ergotropy. Increasing the rotation error reveals an intermediate dynamical regime separating the time-crystalline and Floquet-localized responses. In this regime, most observables exhibit signatures of weakly correlated dynamics reminiscent of Anderson localization.
We report on the first experimental observation of a non-equilibrium phase of matter, the discrete time crystal (DTC). A DTC breaks time-translational symmetry and displays spatio-temporal quantum order in all of its eigenstates, a feature dubbed “eigenstate order”. We implement Floquet dynamics on a 1D chain of 20 superconducting qubits [2]. Engineered disorders in the two-qubit couplings allow many-body localization (MBL) to occur despite strong external drive, thereby stabilizing the non-equilibrium phase [3]. We carefully validate the phase structure of the DTC by probing the average response of all eigenstates belonging to the Floquet unitary. Using a suitable choice of order parameter, we further identify the location of the MBL-ergodicity crossover via experimentally observed finite-size effects. These results open a direct path to studying quantum phase transitions and critical phenomena on NISQ quantum processors.
The simulation of complex quantum many-body systems is a promising short-term goal of noisy intermediate-scale quantum (NISQ) devices. However, the limited connectivity of native qubits hinders the implementation of quantum algorithms that require long-range interactions. We present the outcomes of a digital quantum simulation where we overcome the limitations of the qubit connectivity in NISQ devices. Utilizing the universality of quantum processor native gates, we demonstrate how to implement couplings among physically disconnected qubits at the cost of increasing the circuit depth. We apply this method to simulate a Floquet-driven quantum spin chain featuring interactions beyond nearest neighbors. Specifically, we benchmark the prethermal stabilization of the discrete Floquet time-crystalline response as the interaction range increases, a phenomenon never observed experimentally. Our quantum simulation addresses one of the significant limitations of superconducting quantum processors, namely, device connectivity. It reveals that nontrivial physics involving couplings beyond nearest neighbors can be extracted after the impact of noise is properly taken into account in the theoretical model and consequently mitigated from the experimental data.
Periodically driven quantum many-body systems can spontaneously break discrete time-translation symmetry, realizing discrete time crystals. To date, both experimental and theoretical efforts have largely focused on the simplest case of spontaneous period-doubling in Z2 discrete time crystals realized with qubits. This owes, in part, to the challenge of stabilizing eigenstate order in higher discrete symmetry (Zn) time crystals, due to the presence of richer domain wall physics. Here, we demonstrate the realization of a Z3 discrete time crystal by implementing a Floquet chiral clock model in a chain of 15 superconducting qutrits. Unlike the conventional Ising setting, our system features a tunable chiral angle that governs domain-wall dynamics, spectral degeneracies, and crucially, the stability of time-crystalline order. Using disordered nearest-neighbor chiral interactions, we observe robust subharmonic period tripling that persists across a wide range of drive strengths and is independent of initial state. Finally, we highlight the special role that chirality plays in our Z3 discrete time crystal -- in its absence, the system's Floquet dynamics exhibit a marked initial state dependence governed by domain wall degeneracies. Our results establish native qudit hardware as a powerful platform to access a broader landscape of non-equilibrium phases.