Kicked Tonks–Girardeau Gas Dynamics

• Kicked Tonks–Girardeau gas is a one-dimensional system of impenetrable bosons driven by time-dependent potentials, exhibiting strong correlations via Bose–Fermi mapping.
• It employs periodic kicks or pulse-like drives to explore dynamical localization, fermionization, and the transition from localized to diffusive energy growth.
• Experimental realizations with ultracold atoms validate these effects through observable momentum distributions, coherence measures, and energy saturation signatures.

A kicked Tonks–Girardeau (TG) gas is a paradigm of one-dimensional quantum dynamics in strongly correlated, driven, and out-of-equilibrium systems, realized by subjecting a gas of impenetrable bosons to time-periodic (kicked) or pulse-like driving potentials. The TG limit ($g\rightarrow\infty$) enables a mapping between hard-core bosons and free fermions (Bose–Fermi mapping), juxtaposing strong local correlations (bosonic statistics) with dynamical phenomena such as Anderson localization, effective thermalization, and dynamical fermionization. Kicked TG gases constitute an experimentally tractable setting for exploring many-body localization, delocalization, dynamical phase transitions, high-resolution quantum correlation spectroscopy, and nonthermal excitation, with deep connections to the quantum kicked rotor (QKR) and generalized Floquet systems.

1. Model Hamiltonians and Dynamical Protocols

A kicked TG gas is typically modeled as $N$ bosons of mass $m$ on a ring of length $L=2\pi$ (or in a harmonic trap), with contact interactions in the Lieb–Liniger form,

$$H_{\mathrm{LL}} = \sum_{j=1}^{N}\left[-\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x_j^2} + g\sum_{i<j}\delta(x_i-x_j)\right].$$

In the TG regime ($g\rightarrow\infty$), the system maps to free fermions up to symmetrization: $$\Psi_B(x_1,\ldots,x_N;t) = \prod_{i<j}\operatorname{sgn}(x_i-x_j)\Psi_F(x_1,\ldots,x_N;t).$$ The essential ingredient is a time-dependent driving potential affecting all particles, most notably:

• Periodic or delta-kick drive: $$K \sum_{j=1}^{N} \cos x_j \sum_{n\in\mathbb{Z}} \delta(t-nT)$$
• Quasi-periodic generalization: $$K[1+\varepsilon\cos(\omega_1 t)\cos(\omega_2 t)] \sum_{j=1}^{N} \cos x_j \sum_{n\in\mathbb{Z}} \delta(t-n)$$
• Gaussian pulse drive in a trap: $V_{\mathrm{pulse}}(x,t) = -x S(t) \theta(t)$ with Gaussian envelope $$S(t) = S_A \exp[-(t-t_0)^2/\Delta] \sin[\omega(t-t_0)+\varphi_0]$$

The time evolution is governed stroboscopically via the Floquet operator,

$U_F = e^{-i H_K} e^{-i H_{\mathrm{LL}} T},$

yielding

$\Psi(nT) = (U_F)^n \Psi(0).$

Mapping enables calculation of observables including the one-body density matrix, momentum distributions, and correlation functions (Vuatelet et al., 2022, Li et al., 2024, Yang et al., 17 Jan 2026).

2. Delta-Kick Cooling, Dynamical Fermionization, and Scale-Invariant Evolution

In a harmonic trap with time-dependent frequency $\omega(t)$, the many-body wave function exhibits scale-invariant dynamics: $$\Psi(x_1,\ldots,x_N;t) = b(t)^{-N/2} \exp\left\{ i \frac{\dot{b}}{2b} \sum_j x_j^2 - i\int_0^t \frac{E(0)}{b(t')^2} dt' \right\}\Psi(x_1/b,\ldots,x_N/b;0),$$ where $b(t)$ obeys the Ermakov equation,

$$\ddot{b} + \omega^2(t) b = \omega_0^2 / b^3, \quad b(0)=1,\, \dot{b}(0)=0,$$

and $\omega_0 = \omega(0)$ (2206.13015).

Dynamical fermionization (DF) is observed during expansion—momentum distribution $n(k;t)$ evolves toward the Fermi profile. Delta-kick cooling (DKC) is implemented by a short pulse (harmonic potential of strength $\omega_{\text{kick}}^2$ at $t_k$ for duration $\tau_k$) that erases the accumulated dynamical phase. The protocol rescales the momentum distribution: $n(k;t_k^+) = b(t_k) n\big(b(t_k) k;0\big).$ For $b(t_k)>1$, the initially narrow bosonic momentum profile is restored ("zoom in"); for $b(t_k)<1$, an inverted kick is required. Post-kick, the evolution remains scale-invariant, and $n(k;t)$ is globally rescaled. This enables a "momentum microscope" for quantum correlations, by selecting the scaling parameter $b(t_k)$ (2206.13015).

3. Many-Body Dynamical Localization and Delocalization Transitions

When driven periodically (e.g., kicked rotor protocol; $\varepsilon=0$), the TG gas demonstrates many-body dynamical localization (MBDL): the kinetic energy $E(t) = \langle p^2/2 \rangle$ saturates after initial transient growth,

$\lim_{t\rightarrow\infty} dE/dt \rightarrow 0,$

accompanied by exponential localization in momentum space with localization length $p_{\text{loc}}$ (Vuatelet et al., 2022, Yang et al., 17 Jan 2026).

Under quasi-periodic kicking (e.g., three-frequency QPQKR, $\varepsilon>0$), a dynamical transition akin to the Anderson metal–insulator transition is observed:

• For $K<K_c(\varepsilon)$, the gas remains localized, $E_\infty < \infty$;
• For $K>K_c(\varepsilon)$, energy grows unbounded (diffusive heating), $E(t)\propto t$, corresponding to delocalization;
• At $K\simeq K_c$, sub-diffusive scaling $E(t)\sim t^{2/3}$ is found.

The localization/delocalization threshold is robust against interactions in the TG limit; the phase boundaries for noninteracting and TG systems coincide (Vuatelet et al., 2022). In the localized regime, the steady state is described by an effective Gibbs ensemble: $$\hat{\rho}_{\rm ss} \propto \exp\left[-\frac{H_{\mathrm{TG}} - \mu_{\rm eff} N}{T_{\rm eff}}\right],$$ with $T_{\rm eff}$, $\mu_{\rm eff}$ set by conservation of total particle number and energy.

4. Momentum Distributions, Coherence, and Scaling Phenomena

In the long-time regime under periodic driving, the bosonic momentum distribution $n^B(k)$ shows an exponential low-$k$ decay and a universal $k^{-4}$ Tan-contact tail at large $k$: $$n^B(k) \sim e^{-k/p_{\text{loc}}} \quad (k\to0), \qquad n^B(k) \sim \mathcal{C}_{ss}/k^4.$$ The steady-state contact $\mathcal{C}_{ss}$ links directly to the saturated energy, $\mathcal{C}_{ss}=8N E/L^2$ (Vuatelet et al., 2022, Yang et al., 17 Jan 2026). Coherence length $r_c$ extracted from the one-body density matrix $g^{(1)}_B(r)$ satisfies $g^{(1)}_B(r)\sim \exp(-2r/r_c)$; $r_c$ shrinks as localization strengthens or as initial temperature increases.

At the localization–delocalization transition, the TG system exhibits a breakdown of conventional one-parameter scaling: in the critical regime, two scaling forms emerge for $n_k(t)$, at small/intermediate $k$ and large $k$, with the crossover scale $k_c(t) \sim (N t)^{1/3}$. The fermionic momentum distribution complies with standard scaling, $n_k^F(t) \sim t^{-1/3}f(k t^{-1/3})$ (Airy profile), but bosonic momentum displays distinct exponents and scaling, reflecting strong interactions (Vuatelet et al., 2022).

5. Effects of Finite Temperature

Finite initial temperature $T_0$ modifies, but does not destroy, many-body localization in the kicked TG gas. MBDL persists up to $T_0 \sim \varepsilon_F$ (Fermi energy), though with broader saturated energy and shortened coherence length. Effective thermalization at long times yields emergent temperature $T_{\rm eff}$ and chemical potential $\mu_{\rm eff}$, determined by the late-time energy and density: $$n^F(k) = 1/[e^{(E_k - \mu_{\rm eff}) / T_{\rm eff}} + 1].$$ The scaling of $p_{\text{loc}}/p_F$ with $T_{\rm eff}$ is altered by $T_0$; in the low-$T$ regime (Sommerfeld expansion),

$$p_{\text{loc}}/p_F \simeq \frac{\pi}{2\sqrt{3}}\sqrt{\left(\frac{T_{\rm eff}}{\varepsilon_F}\right)^2 - \left(\frac{T_0}{\varepsilon_F}\right)^2}.$$

No sharp critical $T_c$ is observed for localization under periodic kicks, but in the quasi-periodic regime, the phase boundaries shift only weakly with $T_0$ up to $\sim \varepsilon_F$ (Yang et al., 17 Jan 2026).

6. Pulse-Driven Dynamics, Population Inversion, and Nonthermal States

Applying pulse-like drives (e.g., single or periodic Gaussian pulses) to a harmonically trapped TG gas yields distinct dynamical regimes (Li et al., 2024):

• High-frequency (short) pulses: The system re-phases to its initial ground state after the pulse, with level occupations $n_j(t\gg\mathrm{pulse})\to 1$ for $j=0,\ldots,N-1$. Density and momentum distribution return to those of the TG ground state.
• Low-frequency or long pulses: The system attains a nonthermal, phase-coherent excited state characterized by population inversion: low-lying trap levels ($j< N$) are depopulated while higher levels are populated. The occupation pattern $p_j$ becomes multi-Gaussian, with $N$ distinct excited wavepackets in the trap-level ladder.
• Periodic driving: Observables settle into a Floquet steady state with strictly periodic time-dependence, never thermalizing.

The momentum distribution in all cases retains a single peak around $p=0$ but becomes broader with increasing excitation. The phase-coherent nature of the evolved state is preserved, as evidenced by the persistence of large off-diagonal elements in the single-particle density matrix, even during population inversion (Li et al., 2024).

7. Experimental Realizations and Diagnostics

Realization of the kicked TG gas uses ultracold atoms (e.g., $^{87}$Rb, $^{133}$Cs) in quasi-1D waveguides, with strong transverse confinement ($\omega_\perp$) generating $g\rightarrow\infty$. Axial kicks are implemented via pulsed standing-wave lasers. Temperatures $T_0\sim 50$–$200$ nK correspond to $T_0/\varepsilon_F \sim 0.1$–$1$ in typical experiments.

Experimental observables include:

• Time-of-flight measurements yielding $n^B(k)$ to identify localization, delocalization, and quantum correlations.
• Matter-wave interferometry to extract $g^{(1)}$, providing access to bosonic coherence and correlation lengths.
• Population analysis via spectroscopic means to reconstruct level occupations and confirm multi-Gaussian population inversion in pulse-driven systems.

Key signatures are the saturation (or lack thereof) of kinetic energy under periodic (delocalized under quasi-periodic) drives, scaling of momentum distribution tails, dependence on initial temperature, and direct observation of nonthermal, phase-coherent states following pulse excitation (Vuatelet et al., 2022, 2206.13015, Li et al., 2024, Yang et al., 17 Jan 2026).