Kicked Harmonic Oscillator Dynamics

• Kicked harmonic oscillator is a system where a harmonic potential is impulsively driven at discrete times, producing varied regimes such as bounded motion and ballistic energy growth.
• It employs diverse kicking methods—including momentum transfer, frequency jumps, and measurement-induced operations—to induce squeezing and parametric resonance.
• Analytical tools like Floquet maps and transfer matrices illuminate transitions in energy growth and coherence, offering insights for both theoretical and experimental studies.

A kicked harmonic oscillator is a harmonic degree of freedom subjected to impulsive perturbations applied at isolated times. In the literature represented here, the “kick” may be a direct momentum transfer, a discontinuous jump of the oscillator frequency or other control parameter, a periodic $\delta$-kicked drive, a non-Hermitian $\mathcal{PT}$-symmetric pulse, or a stroboscopic measurement-induced operation on an ancillary system. Across these realizations, the central problems are the stroboscopic map generated by the kicks, the fate of coherence and thermal mixtures, the onset of squeezing and parametric resonance, and the transition between bounded motion, ballistic energy growth, diffusive behavior, and damped oscillation [(Steuernagel, 2011); (Dodonov et al., 2019); (Prado et al., 2016); (Li, 24 Jan 2025); (Montenegro et al., 2021)].

1. Model classes and canonical formulations

The term covers several non-equivalent but structurally related models. In a direct-kick formulation, a mirror of mass $M$ is initially trapped in a harmonic potential,

$$\hat H_{\rm trap}=\frac{\hat p^2}{2M}+\frac12 K_0 \hat x^2,$$

then released into free motion or a weak trap,

$$\hat H_{\rm free}=\frac{\hat p^2}{2M}+\frac12 k \hat x^2,$$

and subsequently driven into a momentum-superposition state by a single-photon interferometric kick (Steuernagel, 2011). In a parametric-kick formulation, the canonical variables $(q,p)$ remain continuous, while the control parameter in the potential changes abruptly, for example by $\omega_1\to \omega_2$ in

$H(q,p;\omega)=\frac{p^2}{2}+\frac12 \omega^2 q^2,$

or by a train of $\delta$-kicks in $\omega^2(t)$ [(Andresas et al., 2013); (Dodonov et al., 2019)].

Periodic kicking naturally leads to Floquet or Floquet-like descriptions. For the $\mathcal{PT}$0-kicked frequency oscillator one writes

$\mathcal{PT}$1

where $\mathcal{PT}$2 is the base frequency, $\mathcal{PT}$3 the kick strength, and $\mathcal{PT}$4 the kick period (Dodonov et al., 2019). In a non-Hermitian extension, the time-dependent Hamiltonian takes the form

$\mathcal{PT}$5

with $\mathcal{PT}$6 the kick strength, $\mathcal{PT}$7 the gain/loss parameter, and $\mathcal{PT}$8 determining resonant versus non-resonant behavior (Li, 24 Jan 2025).

A distinct class replaces Hamiltonian impulses by measurement backaction. In the spin-kicked construction, the background oscillator Hamiltonian remains static,

$\mathcal{PT}$9

while the “kick” is produced by stroboscopic measurements on an ancillary spin after entangling evolution under

$M$0

The resulting map is Floquet-like but not the canonical quantization of a time-dependent kicked Hamiltonian (Montenegro et al., 2021).

2. Impulsive state preparation and parametric action jumps

In the released-mirror problem, the kick is not a change of $M$1 but a coherent momentum superposition generated by a balanced interferometer. The two photon paths impart opposite momenta $M$2, and after tracing out the photon the mirror wave function is multiplied by

$M$3

with controllable interferometric phase $M$4 (Steuernagel, 2011). The initial mirror state is a Boltzmann mixture over trap eigenstates,

$M$5

where $M$6 and $M$7 (Steuernagel, 2011). This formulation makes explicit that thermal occupation broadens the ensemble before the kick, whereas the post-kick dynamics remains unitary.

The parametric-kick formulation emphasizes adiabatic invariants. For the one-dimensional harmonic oscillator the action is

$M$8

Under an instantaneous jump $M$9, the phase-space point $$\hat H_{\rm trap}=\frac{\hat p^2}{2M}+\frac12 K_0 \hat x^2,$$0 is unchanged through the kick, but the energy and action are re-evaluated with the post-kick Hamiltonian (Andresas et al., 2013). Averaging over the initial microcanonical phase yields

$$\hat H_{\rm trap}=\frac{\hat p^2}{2M}+\frac12 K_0 \hat x^2,$$1

with equality only for the trivial case $$\hat H_{\rm trap}=\frac{\hat p^2}{2M}+\frac12 K_0 \hat x^2,$$2 (Andresas et al., 2013). In the same treatment, the Gibbs entropy is written as

$$\hat H_{\rm trap}=\frac{\hat p^2}{2M}+\frac12 K_0 \hat x^2,$$3

so the increase of the mean action implies non-decrease of the ensemble Gibbs entropy under the idealized instantaneous kick (Andresas et al., 2013).

These two formulations clarify a common ambiguity: a “kick” may refer either to an abrupt state-dependent impulse, as in momentum-superposition generation, or to an instantaneous parameter jump that leaves $$\hat H_{\rm trap}=\frac{\hat p^2}{2M}+\frac12 K_0 \hat x^2,$$4 continuous while altering the Hamiltonian itself.

3. Periodic kicking, Floquet maps, and resonance structure

For periodic $$\hat H_{\rm trap}=\frac{\hat p^2}{2M}+\frac12 K_0 \hat x^2,$$5-kicked frequency modulation, the exact dynamics can be reduced to a unimodular transfer matrix. Writing the classical auxiliary solution $$\hat H_{\rm trap}=\frac{\hat p^2}{2M}+\frac12 K_0 \hat x^2,$$6 piecewise between kicks, one obtains a one-period matrix $$\hat H_{\rm trap}=\frac{\hat p^2}{2M}+\frac12 K_0 \hat x^2,$$7 with

$$\hat H_{\rm trap}=\frac{\hat p^2}{2M}+\frac12 K_0 \hat x^2,$$8

and

$$\hat H_{\rm trap}=\frac{\hat p^2}{2M}+\frac12 K_0 \hat x^2,$$9

where $$\hat H_{\rm free}=\frac{\hat p^2}{2M}+\frac12 k \hat x^2,$$0 are Chebyshev polynomials of the second kind (Dodonov et al., 2019). This exact representation organizes the dynamics into three regimes. When $$\hat H_{\rm free}=\frac{\hat p^2}{2M}+\frac12 k \hat x^2,$$1, the energy oscillates but stays finite. When $$\hat H_{\rm free}=\frac{\hat p^2}{2M}+\frac12 k \hat x^2,$$2, the mean energy grows exponentially. At the critical boundary $$\hat H_{\rm free}=\frac{\hat p^2}{2M}+\frac12 k \hat x^2,$$3, parametric resonance gives the quadratic law

$$\hat H_{\rm free}=\frac{\hat p^2}{2M}+\frac12 k \hat x^2,$$4

so that $$\hat H_{\rm free}=\frac{\hat p^2}{2M}+\frac12 k \hat x^2,$$5 (Dodonov et al., 2019).

A two-jump problem provides a complementary exact result. If the frequency changes from $$\hat H_{\rm free}=\frac{\hat p^2}{2M}+\frac12 k \hat x^2,$$6 to $$\hat H_{\rm free}=\frac{\hat p^2}{2M}+\frac12 k \hat x^2,$$7 at $$\hat H_{\rm free}=\frac{\hat p^2}{2M}+\frac12 k \hat x^2,$$8 and back to $$\hat H_{\rm free}=\frac{\hat p^2}{2M}+\frac12 k \hat x^2,$$9 at $(q,p)$0, the oscillator is in a squeezed state at any time $(q,p)$1 when starting from the fundamental state. During $(q,p)$2 the squeezing parameter is

$(q,p)$3

which oscillates with period $(q,p)$4; after the second jump the amplitude freezes, $(q,p)$5, while only the squeeze phase continues to rotate (Tibaduiza et al., 2020).

In the periodically kicked oscillator coupled to a heat bath, the stroboscopic description remains central. The closed system exhibits ballistic energy growth at exact resonance and linear-in-$(q,p)$6 growth off resonance, while dissipation eventually arrests the growth and produces a quasi-stationary cyclic evolution that can be analyzed through Wigner functions at long times (Prado et al., 2016). In the $(q,p)$7-symmetric kicked harmonic oscillator, the Floquet operator is explicitly factored as $(q,p)$8 (Li, 24 Jan 2025). There the resonance condition is not a small perturbation of the Hermitian case: irrational $(q,p)$9 leads to directed current of momentum and ballistic diffusion of energy, whereas integer $\omega_1\to \omega_2$0 yields damped oscillations of momentum and energy with identical frequencies (Li, 24 Jan 2025).

4. Coherence, visibility, and dissipation

The released-mirror problem addresses a specific controversy about “decoherence without dissipation.” The experimentally relevant visibility at the center is defined by comparing the diagonal density matrix element at $\omega_1\to \omega_2$1 for the two interferometric phases $\omega_1\to \omega_2$2 and $\omega_1\to \omega_2$3,

$\omega_1\to \omega_2$4

Although $\omega_1\to \omega_2$5 and later decays because different thermal components accumulate different dynamical phases, the density matrix is still evolving unitarily, so no true decoherence occurs in the sense of irreversible loss of off-diagonal elements (Steuernagel, 2011). If the mirror is released into a weak trap, then at each half-period

$\omega_1\to \omega_2$6

the probability distribution $\omega_1\to \omega_2$7 and the corresponding visibility revive exactly to their initial value, demonstrating perfect recoherence (Steuernagel, 2011).

The same work distinguishes thermal washing-out from genuine loss of coherence. The initial thermal state can obscure interference fringes, but this obscuration does not imply destruction of the underlying superposition. Indeed, the analysis shows that higher temperature can aid pattern formation when higher-$\omega_1\to \omega_2$8 components provide larger spatial width and faster expansion, even though the dynamics remains coherent (Steuernagel, 2011). The associated “entrainment scheme,” in which a single-photon interferometric kick and rapid readout are phase-locked, is described as relatively insensitive to temperature provided the packet is wide enough to carry the full fringe pattern at the moment of kicking (Steuernagel, 2011).

When a heat bath is present, dissipation and decoherence enter explicitly through the Caldeira-Leggett master equation,

$\omega_1\to \omega_2$9

with damping rate $H(q,p;\omega)=\frac{p^2}{2}+\frac12 \omega^2 q^2,$0 and diffusion coefficient $H(q,p;\omega)=\frac{p^2}{2}+\frac12 \omega^2 q^2,$1 (Prado et al., 2016). In this open-system KHO, phase-space methods make it possible to compute high resolution Wigner functions at long times, and the system approaches a quasi-stationary cyclic evolution whose thermodynamic properties can be studied through work, heat, and entropy production per cycle (Prado et al., 2016). A linear “Fourier-law” relation for the heat current is reported in dimensionless variables, and bath-induced decoherence suppresses interference effects that would otherwise sustain ballistic growth or dynamical localization (Prado et al., 2016).

5. Squeezing, tomograms, and exact state reconstruction

Frequency-kicked oscillators provide an exact route to squeezed states. In the two-jump problem, the evolved state for $H(q,p;\omega)=\frac{p^2}{2}+\frac12 \omega^2 q^2,$2 is a squeezed vacuum of the original $H(q,p;\omega)=\frac{p^2}{2}+\frac12 \omega^2 q^2,$3 oscillator, and for $H(q,p;\omega)=\frac{p^2}{2}+\frac12 \omega^2 q^2,$4 the final state has overlap only with even Fock states. The vacuum-persistence probability is

$H(q,p;\omega)=\frac{p^2}{2}+\frac12 \omega^2 q^2,$5

and the transition amplitudes satisfy $H(q,p;\omega)=\frac{p^2}{2}+\frac12 \omega^2 q^2,$6 (Tibaduiza et al., 2020). The physical interpretation given there is that the sudden jump $H(q,p;\omega)=\frac{p^2}{2}+\frac12 \omega^2 q^2,$7 mixes creation and annihilation operators and produces squeezing, while the second jump freezes the squeezing amplitude and leaves only a phase rotation (Tibaduiza et al., 2020).

In the damped Caldirola-Kanai model with a single $H(q,p;\omega)=\frac{p^2}{2}+\frac12 \omega^2 q^2,$8-kick of the frequency,

$H(q,p;\omega)=\frac{p^2}{2}+\frac12 \omega^2 q^2,$9

the kick acts as an instantaneous quadratic phase and hence as a squeezing “kick” (Chernega et al., 2018). The observable content is conveniently described in the tomographic representation,

$\delta$0

which gives the probability distribution of the rotated-scaled quadrature $\delta$1 (Chernega et al., 2018). For weak damping, $\delta$2, suitable $\delta$3 can produce genuine squeezing, meaning $\delta$4 at some times. For strong damping, $\delta$5, a single $\delta$6-kick cannot induce transient squeezing, and the free-particle limit $\delta$7 likewise exhibits no squeezing under a single kick (Chernega et al., 2018).

The periodic Kronig-Penney excitation extends this picture to repeated kicks. Starting from the vacuum, the exact wave function remains a squeezed vacuum state; starting from a coherent state, it becomes a squeezed coherent state, because the $\delta$8-kicks induce squeezing but no displacement (Dodonov et al., 2019). This exact solvability is significant because it ties together transfer-matrix stability, squeezing parameters, and mean-energy growth within one analytic construction.

6. Non-Hermitian, measurement-driven, and related extensions

The $\delta$9-symmetric quantum kicked harmonic oscillator introduces a complex kick potential and thereby a non-Hermitian mechanism for transport and energy growth. In the non-resonant case, where $\omega^2(t)$0 is irrational, the dynamics can be mapped onto an effective tight-binding Hamiltonian in the momentum basis with Hermitian and anti-Hermitian nearest-neighbor hopping terms. The result is a directed momentum current,

$\omega^2(t)$1

together with ballistic energy growth,

$\omega^2(t)$2

(Li, 24 Jan 2025). Under resonant conditions, by contrast, both momentum and energy oscillate as damped cosine functions with identical frequencies, and for the representative parameter set $\omega^2(t)$3, $\omega^2(t)$4 the fitted oscillation frequency is reported as $\omega^2(t)$5 (Li, 24 Jan 2025). The same study states that there is no sharp $\omega^2(t)$6 separating diffusion from oscillation; rather, the relevant boundary is set by rational versus irrational $\omega^2(t)$7 (Li, 24 Jan 2025).

Measurement-induced kicking yields a different nonclassical extension. After each free-evolution interval $\omega^2(t)$8, a projective spin measurement generates Kraus operators

$\omega^2(t)$9

and the ensemble map is

$\mathcal{PT}$00

(Montenegro et al., 2021). This system has “no classical analogous” because the kicks arise from invasive quantum measurements rather than any classical impulsive force (Montenegro et al., 2021). Nevertheless, it exhibits Floquet-like structures: for rational $\mathcal{PT}$01 the Husimi function forms crystalline phase-space patterns, such as a square lattice for $\mathcal{PT}$02 and a 10-fold pattern for $\mathcal{PT}$03, while irrational $\mathcal{PT}$04 produces quasicrystalline and frustrated structures (Montenegro et al., 2021). The same model shows an energy-growth resonance defined by

$\mathcal{PT}$05

with $\mathcal{PT}$06, and an ensemble Loschmidt echo that decays approximately as

$\mathcal{PT}$07

for the parameters of one numerical example (Montenegro et al., 2021).

A broader, but not Hamiltonian, periodically-kicked-oscillator literature supplies a useful contrast. The dissipative limit-cycle model on $\mathcal{PT}$08 studied by Lin and Young has kicks applied to a shear flow rather than to a harmonic Hamiltonian. Its stroboscopic map displays the geometric mechanism “shear + stretch + fold,” horseshoes, positive Lyapunov exponents, and SRB measures on suitable parameter sets (Lin et al., 2010). This suggests that many phenomena often associated with kicked harmonic oscillators—resonances, strange attractors, and sharp transitions in stroboscopic dynamics—belong to a larger theory of impulsively driven oscillatory systems, even though the strictly harmonic and Hamiltonian cases retain their own exact invariants, squeezing structures, and coherence questions (Lin et al., 2010).