Floquet Prethermalization

• Floquet prethermalization is the phenomenon in driven many-body systems where effective static Hamiltonians yield long-lived non-equilibrium states with exponentially suppressed heating.
• Its framework relies on high-frequency expansions of the time-evolution operator, capturing emergent dynamical features like Rabi oscillations, time-crystalline order, and protected symmetries.
• These insights underpin advances in Floquet quantum engineering, stabilizing sensors and quantum simulators while enabling exotic non-equilibrium phase realizations.

Floquet prethermalization is a phenomenon in periodically driven many-body systems where, under suitable conditions, the system avoids immediate thermalization and instead occupies long-lived, non-equilibrium quasi-steady states. These prethermal states are governed by emergent static effective Hamiltonians derived from high-frequency expansions, exhibit exponentially suppressed heating rates, and may host unique dynamical features such as coherent oscillations, time-crystalline order, or protected symmetries. Ultimately, at longer times or stronger drive, the system crosses over to conventional heating and approaches the infinite-temperature ensemble. Floquet prethermalization underpins current approaches to Floquet engineering, stabilization of quantum sensors, and realization of exotic non-equilibrium phases in both quantum and classical materials.

1. Fundamental Theory and High-Frequency Expansion

Floquet prethermalization is anchored in the Floquet–Magnus (or van Vleck) expansion of the time-evolution operator for a system under a time-periodic Hamiltonian $H(t) = H(t + T)$, where $T = 2\pi/\omega$ is the period and $\omega$ is the drive frequency. The evolution over one cycle is

$$U(T) = \mathcal{T} \exp\left[-i \int_0^T H(t') dt'\right] = \exp(-i H_{\rm eff} T),$$

where $H_{\rm eff}$ is the effective Floquet Hamiltonian. For $\omega$ much larger than the typical interaction and bandwidth scales ($\omega \gg J, U$), one can expand $H_{\rm eff}$ in powers of $1/\omega$:

$$H_{\rm eff} = H^{(0)} + \frac{1}{\omega} [H^{(-1)}, H^{(+1)}] + O\left(\frac{1}{\omega^2}\right),$$

where $H^{(m)}$ are Fourier components of $H(t)$ over $[0, T]$. This expansion truncates optimally at an order $n^* \sim \omega / \Lambda$ ($\Lambda$ maximum local energy), yielding an error suppressed as $\sim e^{-c\omega/\Lambda}$ in the remainder and ensuring exponentially small heating rates for times $t \lesssim \tau^* \sim e^{c \omega/\Lambda}$ (Okamoto et al., 2020, Ho et al., 2022, Rubio-Abadal et al., 2020, Beckmann, 2018).

2. Criteria for Prethermal Regimes and Lifetime Scaling

The principal condition for robust Floquet prethermalization is a hierarchy of energy scales: the drive frequency must greatly exceed local interaction strengths and bandwidths:

$\omega \gg \max\{J, U\}$

for quantum lattice models, or the equivalent for classical spins and long-range interactions. In this regime, direct photon absorption is off-resonant from any many-body transition, and heating is governed by high-order processes in $1/\omega$, giving a prethermal plateau with lifetime

$\tau_{\rm pre} \sim \exp(c \omega / J).$

In resonant scenarios where $\omega \sim U$ or matches other transitions, prethermal states can persist for small drive amplitudes ($A_0 < A_{\rm th}$), provided the wavefunction remains confined to a low-dimensional subspace (Okamoto et al., 2020, Herrmann et al., 2017).

For classical systems, the effective Hamiltonian from the Floquet–Magnus expansion accurately governs ensemble-averaged dynamics over exponentially long times, with chaos limiting single-trajectory predictions but preserving Gibbs equilibrium of coarse-grained observables (Ye et al., 2021, Mori, 2018).

A generalized summary table:

Regime	Conditions	Lifetime Scaling
Off-resonant	$\omega \gg J, U$	$\tau_{\rm pre} \sim e^{c \omega / J}$
Resonant, small amplitude	$\omega \approx U$, $A_0 < A_{\rm th}$	Long-lived Rabi oscillations, lifetime set by $A_0$
Resonant, strong drive	$\omega \approx U$, $A_0 \gg A_{\rm th}$	Rapid heating, $\tau_{\rm th}$ short

3. Dynamical Signatures: Rabi Oscillations, Prethermal Plateaux, and Emergent Symmetries

In both quantum and classical systems, Floquet prethermal states manifest as quasi-stationary plateaux in observables (e.g., double occupancy in the Hubbard model, magnetization in spin chains, momentum distributions in rotor systems). At resonance and for small amplitudes, coherent Rabi oscillations dominate the dynamics:

• In two-site Hubbard models, the field couples ground and excited states with a Rabi frequency $\Omega_R = c_1 E_0$, leading to persistent oscillations in double occupation and kinetic energy (Okamoto et al., 2020).
• For $n$-photon resonances, $\Omega_R \propto E_0^n$.
• In full clusters, observables oscillate at $\Omega_R$ about a quasi-stationary average until heating events broaden the active state manifold.

Emergent symmetries—such as discrete time-translation symmetry, $\mathbb{Z}_N$ order from strong pulsed drives, or space-time symmetries—protect additional prethermal phases, including time crystals. In the prethermal regime, spontaneous symmetry breaking of these emergent symmetries can yield robust subharmonic responses or even higher-order time crystals with lifetimes scaling exponentially in drive (Mizuta et al., 2019, Ye et al., 2021, Na et al., 2024).

4. Transition to Thermalization and Role of Drive Strength

Prethermalization persists until the drive amplitude $A_0$ (or field strength $E_0$) exceeds a threshold $A_{\rm th}$, at which multi-photon absorption involving high-order processes becomes significant, dramatically broadening the state's spread over the eigenbasis and initiating heating toward the infinite-temperature ensemble. For the resonantly driven Hubbard cluster, $A_{\rm th} \sim 0.06$–$0.1$ in units $J_0 = 1$ at $\omega \sim 4.5$–$6$; beyond this, Rabi oscillations are quickly damped, and thermal observables saturate ($d \to 0.25$, $E_{\rm kin} \to 0$) (Okamoto et al., 2020).

In classical systems, while ensemble-averaged observables remain well-described by the prethermal Hamiltonian, true chaotic trajectories diverge after short times, and the breakdown for local observables is set by the Lyapunov exponent (Mori, 2018).

5. Analytical Frameworks: Perturbation Theory and Noise Spectroscopy

Time-dependent perturbation theory clarifies the transition from prethermalization to heating:

• First-order terms in the drive yield narrow-band photon absorption and coherent transitions (one-photon Rabi processes).
• Second-order and higher-order terms produce multi-photon resonances, inter-band heating, and eventual approach to infinite temperature (Okamoto et al., 2020, Haldar et al., 2021).
• In disordered, nonlinear Floquet systems (e.g., kicked Gross-Pitaevskii rotors), prethermalization requires narrow, well-separated Floquet bands ($K + W < \pi/2$), as inter-band heating channels are suppressed (Haldar et al., 2021).

Noise spectroscopy via Laplace inversion has provided powerful experimental diagnostics: continuous detection of decay channels, identification of multiple prethermal relaxation rates, and quantification of lifetime extensions in spin ensembles (Harkins et al., 2024). These spectral features correspond to eigenvalues of the dynamical matrix and can be directly mapped to underlying relaxation mechanisms—e.g., dipolar flip-flops, electronic bath coupling, or higher-order Magnus heating.

6. Applications: Floquet Quantum Engineering and Sensing

Floquet prethermalization underpins a broad range of quantum control and metrological applications:

• Stabilization of continuous, reinitialization-free quantum sensors (e.g., $^{13}$C spin-based AC magnetometers, achieving $T_2' \sim 800$ s), where quantum control protocols exploit off-resonance, small-angle pulse trains to suppress heating and enable long-time coherent interrogation (Harkins et al., 2024).
• Quantum simulators and digital quantum computers have accessed and benchmarked Floquet prethermal plateaux in gauge theory, cluster qubit devices, and high-fidelity error-mitigated circuits (Hayata et al., 2024, Morningstar et al., 2021).
• Engineered Floquet protocols allow for the realization of prethermal time crystals, emergent symmetry-protected topological phases, and order-by-disorder selection in frustrated magnets (where dynamical fluctuations favor discrete ordering from degenerate ground states) (Jin et al., 2024, Mizuta et al., 2019).

7. Experimental Evidence and Future Directions

Robust experimental evidence supports Floquet prethermalization:

• Nuclear magnetic resonance techniques in dipolar spin chains and bulk hyperpolarized solids have observed minute-scale plateaux in transverse magnetization, with exponential scaling of lifetimes (Beatrez et al., 2021, Harkins et al., 2024, Peng et al., 2019).
• Ultracold atom experiments in periodically driven Bose-Hubbard systems reported exponential suppression of heating rates and direct validation of statistical resonance theories, even for systems with unbounded local spectra (Rubio-Abadal et al., 2020, Torre et al., 2020).
• Superconducting qubit processors have simulated gauge theories and complex many-body Floquet evolution, experimentally accessing prethermal regimes for hundreds of qubits (Hayata et al., 2024).

Ongoing research explores extensions to quasi-periodic drives (where smoothness in the drive profile crucially determines heating suppression), adiabatic parameter sweeps with dynamical symmetry protection, and the interplay of Floquet engineering with disorder and many-body localization (Ho et al., 2022, He et al., 2022). A frontier direction is the systematic exploitation of dynamical space-time symmetries in prethermal regime design, enabling tailored control of micromotion and the realization of time-domain topological phases (Na et al., 2024).

---

Floquet prethermalization thus offers a rigorous, experimentally validated framework for suppressing heating in periodically driven many-body systems, sustaining non-trivial dynamical states over exponential time windows, and realizing new nonequilibrium phases with tailored Hamiltonians, symmetries, and application-specific functionalities.