Kramers-Henneberger-like Transformation

• Kramers-Henneberger-like transformation is a time-dependent method that shifts light–matter coupling into a displaced binding potential via classical quiver motion.
• The method reformulates the dipole-coupled Hamiltonian into a KH frame, facilitating analysis of stabilization, momentum confinement, and nonadiabatic effects under laser fields.
• It extends to various regimes including short pulses, arbitrary frequencies, curved manifolds, and operator-valued displacements, impacting high-harmonic generation and scattering studies.

The Kramers-Henneberger-like transformation denotes a class of time-dependent unitary or canonical frame changes that follow the laser-driven quiver motion of a charged particle and transfer the explicit light-matter coupling into a time-dependent displacement of the binding potential. In its standard form, the transformation maps a dipole-coupled Hamiltonian to a KH-frame Hamiltonian of the form $H_{\mathrm{KH}}(t)=p^2/(2m)+V(r+\alpha(t))$, where $\alpha(t)$ is the classical quiver displacement. In current usage, “KH-like” also covers extensions to short pulses and arbitrary carrier frequencies, particles constrained to curved manifolds, and operator-valued displacements in quantum-optical settings (Medišauskas et al., 2017).

1. Canonical structure and gauge-equivalent formulations

The standard KH construction begins from a nonrelativistic Hamiltonian in either the length gauge or the velocity gauge, within the dipole approximation. The defining ingredient is the classical displacement $\alpha(t)$, determined by the free-particle equation of motion $m\ddot{\alpha}(t)=-eE(t)$, or equivalently by $\alpha(t)=-\int^t A(\tau)\,d\tau$ in the velocity-gauge representation. The KH transformation is then implemented as a time-dependent translation generated by momentum, often supplemented by a Galilean boost and a scalar phase. In a representative formulation for a linearly polarized monochromatic field, one may write

$$U_{\mathrm{KH}}(t)=\exp\!\left[-\frac{i}{\hbar}m\dot{\alpha}(t)\cdot r\right]\exp\!\left[-\frac{i}{\hbar}\alpha(t)\cdot p\right],$$

which yields, up to a purely time-dependent $c$-number, the transformed Hamiltonian

$H_{\mathrm{KH}}(t)=\frac{p^2}{2m}+V(r+\alpha(t)).$

In one dimension and atomic units, the same content appears as $x=x_e-(E_0/\omega^2)\cos(\omega t)$ and $p=p_e+(E_0/\omega)\sin(\omega t)$, giving $\alpha(t)$0 with $\alpha(t)$1 (Floriani et al., 2024).

The transformation is gauge-equivalent rather than gauge-specific. In the velocity gauge, an additional scalar phase is needed to remove the $\alpha(t)$2 and $\alpha(t)$3 terms and recover the same shifted-potential form. One explicit length-to-KH unitary used in reduced-dimensionality simulations is

$\alpha(t)$4

after which the KH Hamiltonian contains no explicit kinetic ponderomotive offset; the ponderomotive contribution is absorbed into a global phase. This equivalence underlies the recurrent identification of the KH frame with the oscillating, accelerated, or quiver frame (Aynul et al., 2024).

2. Cycle-averaged KH Hamiltonians and validity regimes

For a periodic or nearly periodic drive, the shifted potential $\alpha(t)$5 can be decomposed into Fourier harmonics. The standard high-frequency approximation retains only the zeroth harmonic,

$\alpha(t)$6

leading to the time-independent KH Hamiltonian

$\alpha(t)$7

Its eigenstates define the KH basis and its eigenvalues define KH quasibound levels in the time-averaged picture. In a one-dimensional short-range model studied via the time-dependent Schrödinger equation and Wigner functions, the authors explicitly checked the validity conditions $\alpha(t)$8 and $\alpha(t)$9, using $\alpha(t)$0 a.u., $\alpha(t)$1 a.u., and $\alpha(t)$2 a.u., and concluded that the KH approximation was appropriate for those parameters (Aynul et al., 2024).

This approximation is controlled by more than frequency alone. The same reduced-dimensionality work emphasizes that the vector potential $\alpha(t)$3 and the displacement $\alpha(t)$4 should vanish at the pulse boundaries so that the pulse carries zero net momentum and zero net displacement. Short-pulse formulations likewise require the envelope to vary slowly enough relative to the carrier if one wishes to justify cycle averaging. Where these conditions fail, the time-averaged Hamiltonian ceases to be a complete description and the nonzero Fourier components must be retained.

The most systematic mathematical assessment in the supplied literature is based on Bogolyubov’s averaging theorem. Applied to the KH Hamiltonian, it yields explicit estimates for the difference between full and averaged trajectories as functions of the laser parameters and the region of phase space. A modified Hamiltonian version of the theorem is shown to be better suited to KH systems than the scalar form. The same analysis also makes clear that control deteriorates near the core, where the curvature of the potential is large, so the theorem does not provide rigorous validation of the KH approximation at or very near the minimum of the KH potential, even though it is much more informative away from that region (Floriani et al., 2021).

3. Stabilization, phase space, and momentum confinement

Within the time-averaged KH potential, coherent superpositions of KH eigenstates can generate slow internal motion superimposed on the rapid laser cycle. In the Wigner-phase-space analysis of a short-range KH atom, an equally weighted coherent superposition of the two KH bound states produces beating at the frequency $\alpha(t)$5, and the corresponding Wigner flow performs a cyclic, clockwise motion between the central region and the vicinities of the two KH wells while remaining confined to small $\alpha(t)$6. In the fully time-dependent dynamics, the same cyclic motion survives but is delayed relative to the time-averaged description because turn-on transients distort the wave packet before stabilization sets in (Aynul et al., 2024).

A notable conclusion of that study is that stabilization is associated primarily with momentum confinement rather than with strict spatial trapping. In the KH atom considered there, the Wigner function remains within a small-$\alpha(t)$7 corridor bounded by classical equienergy curves of the time-averaged KH Hamiltonian, even while it spans a broad range in position. Tail-shaped leakage of the Wigner function toward higher momentum, accompanied by alternating-sign interference fringes, is interpreted as a signature of ionization. The same work further reports that the most stable propagation strategy for the full dynamics is initial preparation in the KH ground state, for which the Wigner function stays localized near $\alpha(t)$8 and ionization tails are minimal (Aynul et al., 2024).

The sensitivity of this picture to the shape and symmetry of the underlying potential is illustrated by the Morse-Soft-Coulomb model. In that system, the standard soft-Coulomb potential develops a symmetric double-well KH effective potential and a broad stabilization window, whereas introducing a repulsive Morse branch breaks left-right symmetry and yields only a single effective minimum in the KH potential. The reported consequence is a strong suppression of the stabilization window, together with reduced trapping regions and modified phase-space transport structures (Forlevesi et al., 17 Jun 2026).

4. Beyond naive averaging: short pulses, arbitrary frequencies, and nonperturbative structure

A major extension of the KH paradigm replaces pure cycle averaging by a two-timescale Floquet construction. In that formulation, the slow envelope time $\alpha(t)$9 and fast carrier time $m\ddot{\alpha}(t)=-eE(t)$0 are separated explicitly, and the KH-frame Hamiltonian is expanded only in the fast time. In a plane-wave basis, the Fourier blocks acquire a Toeplitz structure in both momentum and Floquet indices, producing a Block Toeplitz matrix with Toeplitz Blocks. The resulting matrix-vector products are implemented by two-dimensional FFTs on a circulant embedding, giving exact finite-size convolution and $m\ddot{\alpha}(t)=-eE(t)$1 scaling. This permits large Floquet expansions while still exploiting the cycle-averaged KH basis for interpretation (Medišauskas et al., 2017).

The same work argues that pure cycle averaging fails for short pulses and arbitrary frequencies because envelope-induced nonadiabaticity and strong inter-channel couplings are then important. Its two-timescale KH-Floquet formalism was developed precisely to treat that regime. The formalism recovers high-frequency adiabatic stabilization when appropriate, but also describes envelope-driven nonadiabatic excitations, Landau-Zener-Stückelberg interference, and low-frequency virtual multi-channel mixing.

A deeper nonperturbative critique appears in the periodic-orbit approach to “KH scars.” There, the traditional identification of KH states with minima of the cycle-averaged KH potential is challenged directly. The paper shows that the existence of a double well in the averaged KH potential is neither necessary nor sufficient for the existence of a KH state in the full dynamics. Instead, the relevant organizing object is an asymmetric period-$m\ddot{\alpha}(t)=-eE(t)$2 classical periodic orbit localized near one quiver radius. The KH state is effective when that orbit is elliptic, equivalently when the monodromy matrix satisfies $m\ddot{\alpha}(t)=-eE(t)$3 or Greene’s residue satisfies $m\ddot{\alpha}(t)=-eE(t)$4. In the interpolation between the averaged and full Hamiltonians, the fixed point associated with a KH-potential minimum is reported to be destroyed through a saddle-node collision at $m\ddot{\alpha}(t)=-eE(t)$5, well before the full time dependence is restored (Floriani et al., 2024).

5. Generalized KH-like frames in geometric, quantum, and inhomogeneous settings

On a two-dimensional Riemannian manifold embedded in three-dimensional space, the KH transformation ceases to be a global translation. The generalized unitary becomes space- and time-dependent,

$m\ddot{\alpha}(t)=-eE(t)$6

with a simplified generator

$m\ddot{\alpha}(t)=-eE(t)$7

where $m\ddot{\alpha}(t)=-eE(t)$8 is the tangent quiver displacement. The transformed Schrödinger-like equation contains not only the shifted geometric potential but also new time-averaged coefficients of differential operators and operator-valued perturbation terms. At leading order in the high-frequency regime, the effective KH Hamiltonian becomes $m\ddot{\alpha}(t)=-eE(t)$9, where both $\alpha(t)=-\int^t A(\tau)\,d\tau$0 and $\alpha(t)=-\int^t A(\tau)\,d\tau$1 depend explicitly on the geometry through $\alpha(t)=-\int^t A(\tau)\,d\tau$2, $\alpha(t)=-\int^t A(\tau)\,d\tau$3, curvature, and $\alpha(t)=-\int^t A(\tau)\,d\tau$4 (Bendin et al., 2024).

A different generalization quantizes the displacement itself. In the quantum KH construction, the trap center is an operator $\alpha(t)=-\int^t A(\tau)\,d\tau$5 rather than a prescribed classical function. After a time-ordered translation, a squeezing stage, and a momentum-translation stage, the final equation reads

$\alpha(t)=-\int^t A(\tau)\,d\tau$6

Unlike the classical case, the term $\alpha(t)=-\int^t A(\tau)\,d\tau$7 cannot be absorbed as a trivial scalar phase; it generates squeezing of the shaker mode. Because $\alpha(t)=-\int^t A(\tau)\,d\tau$8 and $\alpha(t)=-\int^t A(\tau)\,d\tau$9 do not commute at different times, the renormalized acceleration $$U_{\mathrm{KH}}(t)=\exp\!\left[-\frac{i}{\hbar}m\dot{\alpha}(t)\cdot r\right]\exp\!\left[-\frac{i}{\hbar}\alpha(t)\cdot p\right],$$0 acquires genuine QED corrections controlled by the small parameter $$U_{\mathrm{KH}}(t)=\exp\!\left[-\frac{i}{\hbar}m\dot{\alpha}(t)\cdot r\right]\exp\!\left[-\frac{i}{\hbar}\alpha(t)\cdot p\right],$$1 (Argüello-Luengo et al., 17 Jul 2025).

Spatially inhomogeneous fields generate yet another KH-like variant. In focused beams, the cycle-averaged ponderomotive energy is no longer spatially uniform, and the effective KH-frame Hamiltonian acquires a Stark-like correction $$U_{\mathrm{KH}}(t)=\exp\!\left[-\frac{i}{\hbar}m\dot{\alpha}(t)\cdot r\right]\exp\!\left[-\frac{i}{\hbar}\alpha(t)\cdot p\right],$$2. In the quasistationary approximation used for focused hydrogen, the governing equation becomes

$$U_{\mathrm{KH}}(t)=\exp\!\left[-\frac{i}{\hbar}m\dot{\alpha}(t)\cdot r\right]\exp\!\left[-\frac{i}{\hbar}\alpha(t)\cdot p\right],$$3

The reported result is that even a very small projection of the ponderomotive force along the polarization axis lifts the near-degeneracy of the dichotomic KH ground state and converts the two-lobe density into a single-lobe distribution (Wei, 2018). In a related experimental interpretation of neutral-atom acceleration, the force on the dressed atom includes not only $$U_{\mathrm{KH}}(t)=\exp\!\left[-\frac{i}{\hbar}m\dot{\alpha}(t)\cdot r\right]\exp\!\left[-\frac{i}{\hbar}\alpha(t)\cdot p\right],$$4 but also the KH binding-energy term $$U_{\mathrm{KH}}(t)=\exp\!\left[-\frac{i}{\hbar}m\dot{\alpha}(t)\cdot r\right]\exp\!\left[-\frac{i}{\hbar}\alpha(t)\cdot p\right],$$5, and including that term brought calculated maximum velocities into close agreement with measured values over the full range of pulse durations considered (Wei et al., 2016).

6. Applications, observables, and experimental diagnostics

In high-harmonic generation, the accelerated KH frame shifts the emphasis from the usual three-step picture to the response of a bound density in an oscillating potential. Starting from the velocity gauge and applying the accelerated KH unitary, the Hamiltonian becomes

$$U_{\mathrm{KH}}(t)=\exp\!\left[-\frac{i}{\hbar}m\dot{\alpha}(t)\cdot r\right]\exp\!\left[-\frac{i}{\hbar}\alpha(t)\cdot p\right],$$6

Within the strong-field approximation developed in that frame, the leading-order dipole acceleration is expressed through the static bound-state density and the oscillating force field, and for a one-color linearly polarized drive it reproduces the odd-harmonic selection rule for isotropic targets. The same formulation also shows that, at leading order, the accelerated-KH SFA does not generate the usual long plateau and sharp cutoff associated with the length- or velocity-gauge three-step model (Madsen, 2021).

KH states have also been proposed as experimentally affirmable objects rather than merely formal dressed states. A bichromatic pump-probe strategy uses a high-frequency pump to prepare the KH state and a second probe to ionize it. In the reported VUV single-photon scheme, the double-slit structure of the KH atom is mapped directly into the photoelectron momentum distribution. In the IR tunneling scheme, streaking in the anisotropic KH potential produces a characteristic momentum drift. The same work proposes a third route based on a non-Abelian geometric phase: adiabatic loops in the laser-parameter space can generate a complete spin-flipping transition in the degenerate KH subspace (He et al., 2019).

The KH-like transformation also appears in scattering and reaction theory outside atomic stabilization. In hot quantum plasmas it is used as the semiclassical counterpart of the Block-Nordsieck transformation, converting laser-dressed electron-ion scattering into a stationary problem with a KH-shifted effective potential. There the amplitude of the laser beam, or equivalently the free-electron oscillation amplitude, is found to play an important role in the entanglement fidelity of elastic electron-ion collisions (RoozehdarMogaddam et al., 2020). In deuteron-triton fusion, a KH cycle-averaged treatment of the relative coordinate predicts that the peak of the total fusion cross section shifts from $$U_{\mathrm{KH}}(t)=\exp\!\left[-\frac{i}{\hbar}m\dot{\alpha}(t)\cdot r\right]\exp\!\left[-\frac{i}{\hbar}\alpha(t)\cdot p\right],$$7 keV to $$U_{\mathrm{KH}}(t)=\exp\!\left[-\frac{i}{\hbar}m\dot{\alpha}(t)\cdot r\right]\exp\!\left[-\frac{i}{\hbar}\alpha(t)\cdot p\right],$$8 keV for $$U_{\mathrm{KH}}(t)=\exp\!\left[-\frac{i}{\hbar}m\dot{\alpha}(t)\cdot r\right]\exp\!\left[-\frac{i}{\hbar}\alpha(t)\cdot p\right],$$9, while the angular differential peaks move from zero inclination angle toward $c$0 as $c$1 increases, and the corresponding astrophysical $c$2-factors are enhanced by several times in amplitude (Wu et al., 2021).

Taken together, these formulations show that the KH-like transformation is not a single approximation but a family of frame constructions. What remains invariant across that family is the central idea: the field is traded for a displaced or generalized potential, and the resulting dynamics is organized by the geometry of that dressed frame, whether the relevant structure is a cycle-averaged double well, a periodic orbit, a geometric Laplace-Beltrami operator, a ponderomotive symmetry-breaking term, or an operator-valued quantum displacement.