Kicked Transverse Field Ising Model

• Kicked Transverse Field Ising Model is a driven quantum many-body system characterized by periodic transverse kicks in a spin chain, producing chaotic evolution and high entanglement.
• It employs Floquet dynamics with non-integrable interactions and resonance conditions to reveal detailed quantum chaos and entanglement growth mechanisms.
• The model underpins applications in quantum simulators, enabling benchmarks for random state generation and controlled investigations of quantum heating.

The Kicked Transverse Field Ising Model (kTFIM) refers to a class of quantum many-body spin systems in which a transverse field (typically along the $x$ axis) is applied in a time-periodic, “kicked” (delta-function pulsed) manner, often complemented by a continuous or ramped longitudinal field and/or Ising exchange interactions. These models represent an intersection of quantum Floquet dynamics, quantum chaos, non-integrable many-body evolution, and the study of entanglement and statistical properties characteristic of random pure states.

1. Model Definition and Quantum Map Structure

The Hamiltonian of the kTFIM is typically formulated as

$$\mathcal{H}(t) = J_z \sum_{j=1}^L \sigma_j^z \sigma_{j+1}^z + h_z(t) \sum_{j=1}^L \sigma_j^z + \sum_{k=-\infty}^{\infty} \delta\left(k - \frac{t}{\tau}\right) h_x(t) \sum_{j=1}^L \sigma_j^x,$$

where:

• $J_z$ is the Ising interaction (often set to unity),
• $h_z(t)$ is a longitudinal field, commonly ramped in time, e.g., $h_{z0}\sin(\alpha t)$,
• $h_x(t)$ is the (pulsed or kicked) transverse field, e.g., $h_{x0}\cos(\alpha t)$,
• $\delta$ is the Dirac delta, generating kicks at intervals $\tau$,
• $\sigma_j^{x,z}$ are Pauli operators at site $$\mathcal{H}(t) = J_z \sum_{j=1}^L \sigma_j^z \sigma_{j+1}^z + h_z(t) \sum_{j=1}^L \sigma_j^z + \sum_{k=-\infty}^{\infty} \delta\left(k - \frac{t}{\tau}\right) h_x(t) \sum_{j=1}^L \sigma_j^x,$$0.

Between kicks, the system evolves under $$\mathcal{H}(t) = J_z \sum_{j=1}^L \sigma_j^z \sigma_{j+1}^z + h_z(t) \sum_{j=1}^L \sigma_j^z + \sum_{k=-\infty}^{\infty} \delta\left(k - \frac{t}{\tau}\right) h_x(t) \sum_{j=1}^L \sigma_j^x,$$1. At kick times, all spins experience a non-commuting rotation from the transverse field. The stroboscopic evolution is described by the Floquet (quantum map) sequence: $$\mathcal{H}(t) = J_z \sum_{j=1}^L \sigma_j^z \sigma_{j+1}^z + h_z(t) \sum_{j=1}^L \sigma_j^z + \sum_{k=-\infty}^{\infty} \delta\left(k - \frac{t}{\tau}\right) h_x(t) \sum_{j=1}^L \sigma_j^x,$$2 with

$$\mathcal{H}(t) = J_z \sum_{j=1}^L \sigma_j^z \sigma_{j+1}^z + h_z(t) \sum_{j=1}^L \sigma_j^z + \sum_{k=-\infty}^{\infty} \delta\left(k - \frac{t}{\tau}\right) h_x(t) \sum_{j=1}^L \sigma_j^x,$$3

$$\mathcal{H}(t) = J_z \sum_{j=1}^L \sigma_j^z \sigma_{j+1}^z + h_z(t) \sum_{j=1}^L \sigma_j^z + \sum_{k=-\infty}^{\infty} \delta\left(k - \frac{t}{\tau}\right) h_x(t) \sum_{j=1}^L \sigma_j^x,$$4

This checkerboard of unitary operations reflects the interplay between non-commuting dynamical generators and defines a non-integrable, time-crystalline structure even for fixed $$\mathcal{H}(t) = J_z \sum_{j=1}^L \sigma_j^z \sigma_{j+1}^z + h_z(t) \sum_{j=1}^L \sigma_j^z + \sum_{k=-\infty}^{\infty} \delta\left(k - \frac{t}{\tau}\right) h_x(t) \sum_{j=1}^L \sigma_j^x,$$5 and $$\mathcal{H}(t) = J_z \sum_{j=1}^L \sigma_j^z \sigma_{j+1}^z + h_z(t) \sum_{j=1}^L \sigma_j^z + \sum_{k=-\infty}^{\infty} \delta\left(k - \frac{t}{\tau}\right) h_x(t) \sum_{j=1}^L \sigma_j^x,$$6.

The state at stroboscopic times $$\mathcal{H}(t) = J_z \sum_{j=1}^L \sigma_j^z \sigma_{j+1}^z + h_z(t) \sum_{j=1}^L \sigma_j^z + \sum_{k=-\infty}^{\infty} \delta\left(k - \frac{t}{\tau}\right) h_x(t) \sum_{j=1}^L \sigma_j^x,$$7 is recursively,

$$\mathcal{H}(t) = J_z \sum_{j=1}^L \sigma_j^z \sigma_{j+1}^z + h_z(t) \sum_{j=1}^L \sigma_j^z + \sum_{k=-\infty}^{\infty} \delta\left(k - \frac{t}{\tau}\right) h_x(t) \sum_{j=1}^L \sigma_j^x,$$8

2. Dynamical Protocols, Entanglement Growth, and Random State Statistics

The kTFIM supports distinct regimes depending on field protocols:

Quench Protocol (Pulsed Fields):

• $$\mathcal{H}(t) = J_z \sum_{j=1}^L \sigma_j^z \sigma_{j+1}^z + h_z(t) \sum_{j=1}^L \sigma_j^z + \sum_{k=-\infty}^{\infty} \delta\left(k - \frac{t}{\tau}\right) h_x(t) \sum_{j=1}^L \sigma_j^x,$$9 is swept from a maximal value to zero, while $J_z$0 ramps from zero to a peak—typically $J_z$1, $J_z$2, $J_z$3.
• Initial state is usually the product eigenstate of $J_z$4 (all spins polarized along $J_z$5), closely aligned with the maximal transverse field ground state.
• Observables:
  • Nearest-neighbor concurrence $J_z$6: Two-qubit entanglement measure; vanishes as longitudinal field increases.
  • Q-measure: $J_z$7, quantifies single-site purity loss, rising as multipartite entanglement spreads.
  • Block Entropy ($J_z$8): $J_z$9, measures half-chain entanglement; saturates near the random state value, as given by Page's formula $h_z(t)$0 for $h_z(t)$1-dimensional subsystems [$h_z(t)$2].

Time-evolved states rapidly converge to statistical properties mirroring Haar-random pure states: (i) Block entropy approaches Page’s value, (ii) reduced density matrix eigenvalue distributions match Marchenko–Pastur law, and (iii) wavefunction intensities follow exponential statistics.

Steady Protocol (Fixed Fields):

• Both $h_z(t)$3 and $h_z(t)$4 are time-independent.
• Time evolution is purely by repeated application of a Floquet operator $h_z(t)$5.
• Observed entanglement is non-monotonic (“U”-shaped in kick period $h_z(t)$6) and statistical properties can deviate from random-matrix predictions, especially with periodic boundary conditions.

3. Quantum Resonance and Kicking Interval Effects

A remarkable non-ergodic effect emerges at specific $h_z(t)$7: $h_z(t)$8 At $h_z(t)$9, with periodic boundaries,

$h_{z0}\sin(\alpha t)$0

rendering the interaction trivial (up to a phase). Here, entanglement generation is "frozen," even in a system that is non-integrable away from these points. Such resonance results in stroboscopic "pauses" or suppression of scrambling, paralleling quantum resonances in the kicked rotor.

This freezing is robust even when the full time evolution away from resonance is highly entangling and drives the system to nearly maximal multipartite entanglement.

4. Role of Integrability Breaking and Nonintegrable Dynamics

The addition of a longitudinal field $h_{z0}\sin(\alpha t)$1 breaks the integrability of the free-fermion solution (contrasting the exactly solvable static TFIM). This promotes nontrivial entanglement growth and state complexity. The combined effect of non-commuting pulses (transverse field kicks) and a ramping $h_{z0}\sin(\alpha t)$2 creates a regime in which the conventional mapping to quadratic fermions fails, and eigenstate statistics exhibit Wigner-Dyson level statistics and random matrix universality.

Numerically, as multipartite entanglement grows, the concurrence drops to zero, and the system's density matrices yield Marchenko–Pastur eigenvalue distributions, providing evidence for the generation of random-like states [$h_{z0}\sin(\alpha t)$3].

5. Impact of Boundary Conditions and Spectral Statistics

Recent spectral analysis of related kicked Ising models at self-dual points reveals a strong sensitivity to boundary conditions [$h_{z0}\sin(\alpha t)$4]:

• Periodic Boundaries: The trace of the Floquet propagator behaves as a real Gaussian (higher moments enhanced: $h_{z0}\sin(\alpha t)$5), due to extra conjugate symmetries.
• Open Boundaries: Trivializes extra symmetry, so the trace is complex Gaussian ($h_{z0}\sin(\alpha t)$6), the standard COE prediction.

This boundary-driven crossover affects the interpretation of spectral fluctuations, relevant in assessing quantum chaos.

6. Connections to Quantum Chaos, Ergodicity, and Statistical Mechanics

The kTFIM serves as a paradigmatic system for studying the emergence of quantum chaos in Floquet systems:

• Level Statistics: As in higher-dimensional kicked models [$h_{z0}\sin(\alpha t)$7], level spacing distributions closely follow the Wigner surmise (COE predictions) except at integrable or resonant points.
• Dynamical Susceptibility: Time-averaged autocorrelation functions of local observables (e.g., magnetization) can exhibit sharp transitions between ergodic (decaying) and non-ergodic (persisting) dynamics, sometimes decoupled from global spectral transitions (flatness of spectral density).
• Out-of-time-order correlators (OTOCs) and Lyapunov Exponents: In generalizations with all-to-all or $h_{z0}\sin(\alpha t)$8-body interactions, short-time OTOC growth rates can be extracted from classical Lyapunov exponents, confirming a semiclassical-to-quantum chaos correspondence [$h_{z0}\sin(\alpha t)$9].

Observable	Quench (Pulsed)	Steady (Fixed)
Block Entropy, Q-measure	Grows to Page/random state	Non-monotonic, “U”-shaped
State Intensity Statistics	Exponential law (random state)	Possible deviations
Reduced Density Eigenvalues	Marchenko–Pastur law	Deviations (periodic BCs)
Two-site concurrence	Drops to zero (large $h_x(t)$0)	Oscillatory, dips at resonance

7. Applications, Experimental Realizations, and Extensions

• Quantum Simulators: The protocols described—especially quenching and kicking—map directly onto cold atom, trapped ion, and superconducting qubit platforms where time-dependent control of transverse and longitudinal fields is routine.
• Benchmarking and Random State Generation: The kTFIM dynamically generates multipartite random states, critical for benchmarking quantum devices, testing volume-law entanglement scaling, and as state-initialization for quantum information tasks.
• Probing Quantum Phase Transitions: Through tuning pulse sequences and field ramps, the kTFIM provides a testbed for investigating nonequilibrium Kibble-Zurek scalings and symmetry breaking gap detection protocols.
• Control of Quantum Heating: The identification of resonance points and corresponding dynamical freezing suggests mechanisms for mitigating Floquet heating or for engineering “prethermal” plateaus in driven quantum many-body systems.

8. Summary and Outlook

The kicked transverse field Ising model embodies a tightly controlled avenue for exploring the boundary between fully quantum chaotic, highly entangling evolution and the existence of precise, dynamically stabilized or “frozen” regimes. Its study highlights:

• The role of time-dependent protocols (quench versus steady),
• The importance of competing non-commuting terms,
• The emergence of quantum resonance phenomena and their analytic characterization,
• The dynamical generation of nearly-random quantum states exhibiting maximal multipartite entanglement,
• The sensitivity of dynamical and statistical observables—especially spectral moments and entanglement—to boundary conditions and protocol details.

Ongoing research explores further generalizations (e.g., $h_x(t)$1-body interactions, monitored nonunitary evolution), connections to quantum thermodynamic universality, and implications for the stability of quantum simulation and quantum information processing in kicked many-body systems [$h_x(t)$2, $h_x(t)$3, $h_x(t)$4].