Logarithmically Slow Relaxation in Quasiperiodically Driven Random Spin Chains

We simulate the dynamics of a disordered interacting spin chain subject to a quasiperiodic time-dependent drive, corresponding to a stroboscopic Fibonacci sequence of two distinct Hamiltonians. Exploiting the recursive drive structure, we can efficiently simulate exponentially long times. After an initial transient, the system exhibits a long-lived glassy regime characterized by a logarithmically slow growth of entanglement and decay of correlations analogous to the dynamics at the many-body delocalization transition. Ultimately, at long time scales, which diverge exponentially for weak or rapid drives, the system thermalizes to infinite temperature. The slow relaxation enables metastable dynamical phases, exemplified by a time quasicrystal in which spins exhibit persistent oscillations with a distinct quasiperiodic pattern from that of the drive. We show that in contrast with Floquet systems, a high-frequency expansion strictly breaks down above fourth order, and fails to produce an effective static Hamiltonian that would capture the prethermal glassy relaxation.