Classical Eikonal: Wavefronts & Scattering

• Classical eikonal is a unifying phase framework linking wavefront propagation and scattering processes across optics, kinetic theory, and relativistic settings.
• It reduces complex wave equations to first-order nonlinear forms, enabling a clear characterization of wavefronts, rays, and transport limits.
• The approach employs symmetry analysis and variable transformations to connect microscopic transport dynamics with macroscopic classical observables in scattering.

Classical eikonal denotes a family of closely related structures that arise when wave propagation or scattering is reduced to phase data. In geometrical optics and Hamilton–Jacobi theory, it is the first-order nonlinear equation satisfied by a rapidly varying phase, typically $|\nabla S(\mathbf r)|^2=n^2(\mathbf r)$ or, in relativistic notation, $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$ (Kulyabov et al., 2019, Xiao et al., 2010). In kinetic theory, a hyperbolic WKB limit of transport with BGK relaxation produces a non-quadratic Hamilton–Jacobi equation that is only asymptotically equivalent to the classical quadratic eikonal for small gradients (Bouin et al., 2012). In high-energy scattering, the same term refers to an impact-parameter phase $\delta(s,b)$ entering $\tilde{\mathcal S}(s,b)\approx e^{2i\delta(s,b)}$, or, in a stricter canonical formulation, to a generator $\chi$ of the classical scattering map $A_{\rm out}=e^{\{\chi,\}}A_{\rm in}$ (Cadoni et al., 20 Apr 2026, Kim et al., 7 Nov 2025). The modern literature therefore treats “classical eikonal” not as a single equation, but as a unifying Hamilton–Jacobi and phase-space paradigm linking wavefronts, rays, transport limits, and classical observables in scattering.

1. Wavefront equations and geometric-optics foundations

In optical and PDE settings, the classical eikonal equation is the leading-order WKB reduction of a wave equation. For the scalar Helmholtz equation,

$$\nabla^2 \psi(\mathbf r)+k_0^2 n^2(\mathbf r)\psi(\mathbf r)=0,$$

the short-wavelength ansatz

$\psi(\mathbf r)=A(\mathbf r)e^{ik_0 S(\mathbf r)}$

yields, at order $O(k_0^2)$,

$|\nabla S(\mathbf r)|^2=n^2(\mathbf r),$

while the next order gives the transport equation

$g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$0

(Kulyabov et al., 2019). The level sets $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$1 are wavefronts, and rays are the characteristic curves orthogonal to those fronts.

The relativistic version appears in the classical theory of wavefronts in Maxwell electrodynamics. If a wavefront is described by

$g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$2

then the Lorentz-invariant eikonal equation is

$g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$3

equivalently

$g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$4

which is the Hamilton–Jacobi equation for massless propagation (Xiao et al., 2010).

A parallel mathematical tradition studies the relativistic and unit eikonal equations

$g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$5

with

$g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$6

in Minkowski signature (Yehorchenko, 2022). In this language, the null equation describes lightlike gradients, while the unit equation fixes the Minkowski norm of the gradient to $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$7. The same source notes that these are classical eikonal equations in the sense of geometrical optics and Hamilton–Jacobi theory.

Context	Representative object	Role
Geometrical optics	$g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$8	Wavefronts and rays
Relativistic PDE	$g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$9, $\delta(s,b)$0	Null and unit eikonal equations
Kinetic transport	$\delta(s,b)$1	Velocity-averaged eikonal analogue
Impact-parameter scattering	$\delta(s,b)$2	Classical phase shift
Canonical scattering map	$\delta(s,b)$3	Classical scattering generator

These forms are not interchangeable, but they are structurally linked by the fact that each encodes propagation through a phase function whose derivatives determine observable trajectories or fronts.

2. Solution theory, characteristics, and symmetry structure

For the unit relativistic eikonal equation $\delta(s,b)$4, a general local solution can be organized by the rank of the Hessian $\delta(s,b)$5 (Yehorchenko, 2022). The rank-zero sector consists of linear solutions

$\delta(s,b)$6

For maximal spatial rank, the general solution is given parametrically by

$\delta(s,b)$7

together with the implicit constraints

$\delta(s,b)$8

where $\delta(s,b)$9 are parameters and $\tilde{\mathcal S}(s,b)\approx e^{2i\delta(s,b)}$0 is an arbitrary sufficiently smooth function. Intermediate ranks $\tilde{\mathcal S}(s,b)\approx e^{2i\delta(s,b)}$1 admit analogous parametric forms involving auxiliary functions $\tilde{\mathcal S}(s,b)\approx e^{2i\delta(s,b)}$2. This rank stratification makes the general solution resemble characteristic data encoded by a generating function.

The same paper develops the solution theory by hodograph and contact transformations. A general contact transformation is written as

$\tilde{\mathcal S}(s,b)\approx e^{2i\delta(s,b)}$3

and is used to convert first-order nonlinear eikonal equations into simpler equations in new variables. In two dimensions, the Euclidean eikonal

$\tilde{\mathcal S}(s,b)\approx e^{2i\delta(s,b)}$4

is reduced to a parametric solution through a Legendre-type transformation; the higher-dimensional relativistic case follows the same pattern.

The symmetry theory is equally rich. For the null equation $\tilde{\mathcal S}(s,b)\approx e^{2i\delta(s,b)}$5, the maximal Lie symmetry algebra is infinite dimensional and contains transformations

$\tilde{\mathcal S}(s,b)\approx e^{2i\delta(s,b)}$6

with $\tilde{\mathcal S}(s,b)\approx e^{2i\delta(s,b)}$7 a conformal Killing field and $\tilde{\mathcal S}(s,b)\approx e^{2i\delta(s,b)}$8 arbitrary smooth functions. A distinctive property is that any function of a solution is again a solution. For $\tilde{\mathcal S}(s,b)\approx e^{2i\delta(s,b)}$9, the maximal Lie symmetry algebra is finite dimensional but contains translations, Lorentz transformations, dilations, and special conformal transformations in an extended $\chi$0-dimensional space where the dependent variable $\chi$1 is treated as an additional coordinate. Discrete symmetries include time reflection, space reflections, sign reversal $\chi$2, permutations of spatial coordinates, and local hodograph transformations.

In optical applications, the same characteristic structure becomes computational. The numerical study of Maxwell and Luneburg lenses rewrites the eikonal equation as a characteristic ODE system. In Cartesian coordinates, with $\chi$3,

$\chi$4

and one further has

$\chi$5

This identifies the ray parameter $\chi$6 with the eikonal $\chi$7 up to an additive constant (Kulyabov et al., 2019).

3. Kinetic and transport generalizations

A nontrivial generalization of the classical eikonal arises from kinetic transport with BGK relaxation (Bouin et al., 2012). The starting model is

$\chi$8

with $\chi$9 bounded and symmetric, and $A_{\rm out}=e^{\{\chi,\}}A_{\rm in}$0 a symmetric normalized Maxwellian. Under the hyperbolic scaling

$A_{\rm out}=e^{\{\chi,\}}A_{\rm in}$1

the rescaled equation becomes

$A_{\rm out}=e^{\{\chi,\}}A_{\rm in}$2

The WKB/Hopf–Cole ansatz

$A_{\rm out}=e^{\{\chi,\}}A_{\rm in}$3

leads to a kinetic phase equation. As $A_{\rm out}=e^{\{\chi,\}}A_{\rm in}$4, $A_{\rm out}=e^{\{\chi,\}}A_{\rm in}$5 converges locally uniformly to a limit $A_{\rm out}=e^{\{\chi,\}}A_{\rm in}$6 independent of $A_{\rm out}=e^{\{\chi,\}}A_{\rm in}$7, and $A_{\rm out}=e^{\{\chi,\}}A_{\rm in}$8 is the viscosity solution of the implicit Hamilton–Jacobi equation

$A_{\rm out}=e^{\{\chi,\}}A_{\rm in}$9

By monotonicity, this is equivalent to

$$\nabla^2 \psi(\mathbf r)+k_0^2 n^2(\mathbf r)\psi(\mathbf r)=0,$$0

where the effective Hamiltonian is defined implicitly by

$$\nabla^2 \psi(\mathbf r)+k_0^2 n^2(\mathbf r)\psi(\mathbf r)=0,$$1

The paper explicitly presents this as a kinetic analogue of the classical eikonal equation obtained from the heat equation with small diffusivity,

$$\nabla^2 \psi(\mathbf r)+k_0^2 n^2(\mathbf r)\psi(\mathbf r)=0,$$2

The relation is asymptotic rather than exact: $$\nabla^2 \psi(\mathbf r)+k_0^2 n^2(\mathbf r)\psi(\mathbf r)=0,$$3 Outside the small-gradient regime, the kinetic Hamiltonian is generally non-quadratic. In one dimension with $$\nabla^2 \psi(\mathbf r)+k_0^2 n^2(\mathbf r)\psi(\mathbf r)=0,$$4 and $$\nabla^2 \psi(\mathbf r)+k_0^2 n^2(\mathbf r)\psi(\mathbf r)=0,$$5,

$$\nabla^2 \psi(\mathbf r)+k_0^2 n^2(\mathbf r)\psi(\mathbf r)=0,$$6

while in the two-velocity model

$$\nabla^2 \psi(\mathbf r)+k_0^2 n^2(\mathbf r)\psi(\mathbf r)=0,$$7

Several structural properties distinguish the kinetic Hamiltonian from the classical quadratic one. It is convex, globally Lipschitz, and satisfies

$$\nabla^2 \psi(\mathbf r)+k_0^2 n^2(\mathbf r)\psi(\mathbf r)=0,$$8

This bounded gradient expresses finite propagation speeds inherited from the bounded velocity set $$\nabla^2 \psi(\mathbf r)+k_0^2 n^2(\mathbf r)\psi(\mathbf r)=0,$$9. A plausible implication is that the kinetic eikonal retains the Hamilton–Jacobi form while encoding microscopic transport and scattering information through the effective Hamiltonian rather than through a universal quadratic norm.

4. Impact-parameter phases and the scattering eikonal

In scattering theory, the classical eikonal appears in high-energy, small-angle, large-impact-parameter regimes where multiple soft exchanges exponentiate into an impact-parameter phase (Cadoni et al., 20 Apr 2026, Bellazzini et al., 2022). For gravitational $\psi(\mathbf r)=A(\mathbf r)e^{ik_0 S(\mathbf r)}$0 scattering, one writes

$\psi(\mathbf r)=A(\mathbf r)e^{ik_0 S(\mathbf r)}$1

with $\psi(\mathbf r)=A(\mathbf r)e^{ik_0 S(\mathbf r)}$2 the transverse impact parameter. In the nonlocal-gravity analysis of scalar scattering, the leading eikonal phase is extracted from the Born amplitude by

$\psi(\mathbf r)=A(\mathbf r)e^{ik_0 S(\mathbf r)}$3

Classical momentum transfer and deflection then follow from

$\psi(\mathbf r)=A(\mathbf r)e^{ik_0 S(\mathbf r)}$4

The gravitational case admits a direct geometric interpretation. In massless scattering, the eikonal phase can be viewed as the phase shift of a particle crossing an Aichelburg–Sexl shock wave, with

$\psi(\mathbf r)=A(\mathbf r)e^{ik_0 S(\mathbf r)}$5

In the nonlocal theory with Gaussian form factor $\psi(\mathbf r)=A(\mathbf r)e^{ik_0 S(\mathbf r)}$6, the resulting eikonal regularizes the short-distance behavior and is reinterpreted in terms of generalized Aichelburg–Sexl geometries and smeared linearized Schwarzschild metrics. This same eikonal data is then used to motivate regular black-hole spacetimes with de Sitter cores (Cadoni et al., 20 Apr 2026).

In AdS–Schwarzschild scattering of a highly energetic light probe off a heavy object, the eikonal phase is the classical action of a null geodesic,

$\psi(\mathbf r)=A(\mathbf r)e^{ik_0 S(\mathbf r)}$7

equivalently

$\psi(\mathbf r)=A(\mathbf r)e^{ik_0 S(\mathbf r)}$8

and its derivatives yield classical observables,

$\psi(\mathbf r)=A(\mathbf r)e^{ik_0 S(\mathbf r)}$9

In the flat-space limit, the phase develops a large imaginary part for sufficiently small impact parameter, and the inelastic cross-section equals the classical absorption cross-section of the black hole (Parnachev et al., 2020).

The large-angular-momentum origin of eikonal exponentiation has also been reformulated group-theoretically. In the large-$O(k_0^2)$0 limit, the contraction

$O(k_0^2)$1

maps the partial-wave description to the isometries of the transverse plane; the Wigner $O(k_0^2)$2-matrices become Bessel functions, and the eikonal transform becomes a two-dimensional Fourier transform in impact-parameter space (Bellazzini et al., 2022). That paper further distinguishes resolvable from non-resolvable corrections: pure quantum-gravity corrections $O(k_0^2)$3 are stated to be never resolvable in the transplanckian eikonal regime, whereas gauge or tidal corrections can dominate higher post-Minkowskian terms.

Infrared structure is part of the same classical story. Combining impact-parameter eikonal exponentiation with momentum-space IR exponentiation, one can determine the divergent part of the two-loop eikonal from tree-level and one-loop elastic amplitudes in $O(k_0^2)$4 supergravity and in general relativity. In that framework, $O(k_0^2)$5 controls conservative dynamics, while the IR-divergent part of $O(k_0^2)$6 is directly linked to soft radiation, the zero-frequency limit of the energy spectrum, and radiation-reaction effects (Heissenberg, 2021).

5. Phase, generator, Magnusian, and generalized channel structure

A central terminological issue is that amplitudes literature uses “eikonal” for two inequivalent objects (Kim et al., 7 Nov 2025). The eikonal phase is defined by

$O(k_0^2)$7

and in the classical limit it coincides with the on-shell action,

$O(k_0^2)$8

By contrast, the eikonal generator or Magnusian is the classical limit of the logarithm of the $O(k_0^2)$9-matrix operator,

$|\nabla S(\mathbf r)|^2=n^2(\mathbf r),$0

and it generates the scattering map in phase space,

$|\nabla S(\mathbf r)|^2=n^2(\mathbf r),$1

These two objects are inequivalent in general. The exact relation is

$|\nabla S(\mathbf r)|^2=n^2(\mathbf r),$2

or equivalently

$|\nabla S(\mathbf r)|^2=n^2(\mathbf r),$3

Only in integrable scattering do $|\nabla S(\mathbf r)|^2=n^2(\mathbf r),$4 and $|\nabla S(\mathbf r)|^2=n^2(\mathbf r),$5 coincide up to a Legendre transform; with spin or radiation, that simplification fails.

The same distinction is operationalized in the Magnus-expansion approach to the classical eikonal (Kim et al., 2024). There the classical eikonal is defined as the generator of the canonical transformation from in-states to out-states,

$|\nabla S(\mathbf r)|^2=n^2(\mathbf r),$6

and computed from a Magnus expansion in oriented tree graphs whose edges encode retarded and advanced worldline propagators. The paper develops a Hopf-algebra-based algorithm that computes the tree coefficients up to the $|\nabla S(\mathbf r)|^2=n^2(\mathbf r),$7th order—“over half a million trees”—in less than an hour, applies the method to the 3PM gravitational eikonal, and emphasizes that the Magnus prescription yields a finite eikonal whereas the naïve time-symmetric one is IR-divergent from 3PM on.

The worldline formulation provides a bridge between amplitude and canonical perspectives. In the worldline/QFT comparison, the eikonal method and WQFT are shown to be equivalent descriptions of the classical limit of $|\nabla S(\mathbf r)|^2=n^2(\mathbf r),$8 scattering with massless mediators; reducible and irreducible WQFT diagrams map one-to-one to the exponentiation structure of the impact-parameter amplitude (Ajith et al., 2024). This suggests that the classical eikonal can be treated either as a phase extracted recursively from amplitudes or as a directly computable generator in worldline variables.

Once spin or inelasticity are admitted, the eikonal becomes matrix-valued or channel-valued. For spinning particles, one assumes a matrix eikonal exponentiation

$|\nabla S(\mathbf r)|^2=n^2(\mathbf r),$9

and classical observables are obtained from derivatives and Poisson brackets of the spinning eikonal, including the shifted variables $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$00 and $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$01 (Luna et al., 2023). For inelastic coupled-channel eikonal scattering, the diagonal elastic block becomes

$g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$02

while off-diagonal blocks encode channel transitions. Unitarity then implies

$g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$03

so the purely elastic unit-modulus phase is recovered only when inelastic channels are absent (Aoude et al., 2024).

6. Causality, polarization, and effective geometry

Classical eikonal data often doubles as a causality diagnostic. In higher-derivative effective theories of gravity, helicity-flip processes can survive the eikonal limit and promote the phase to a $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$04 matrix in helicity space (Huber et al., 2020). For scattering of gravitons or photons of frequency $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$05 off a heavy scalar of mass $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$06, the impact-parameter phase matrix is diagonalized, and its eigenvalues are used to compute the classical deflection angle and Shapiro time delay/advance to 2PM order via

$g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$07

Whenever the classical expectation of helicity conservation is violated, the eikonal eigenvalues become non-degenerate, and time advance becomes a generic possibility at small impact parameter. In that analysis, graviton scattering in $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$08 and $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$09 avoids time advance if the couplings satisfy

$g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$10

whereas graviton scattering in $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$11 and photon scattering in $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$12 exhibit unavoidable time advance in the EFT regime.

A related causality analysis appears in Lorentz-violating Maxwell theory (Xiao et al., 2010). There the wavefront method yields a modified eikonal condition $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$13, equivalent to $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$14, and in a simple anisotropic example the wavefront function $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$15 obeys

$g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$16

The same analysis shows that the wavefront velocity equals the group velocity and the energy-flow velocity, and that there always exists a mode with signal velocity $g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0$17; the paper therefore concludes that causality is violated classically in that model.

The causal interpretation of eikonal scattering in gravity can also be phrased positively. In the large-angular-momentum regime, analyticity and unitarity imply non-linear positivity bounds on the eikonal amplitude, with positivity of time delay as the simplest case (Bellazzini et al., 2022). In AdS–Schwarzschild, the onset of a large imaginary part in the eikonal phase at impact parameters below the photon sphere is not presented as superluminality, but as the transition from elastic scattering to absorption, and the resulting inelastic cross-section equals the classical geometric absorption cross-section (Parnachev et al., 2020).

These examples make clear that the classical eikonal is not merely a calculational shortcut. It encodes effective geometry, polarization dependence, propagation speed, absorption, and the distinction between acceptable and pathological EFT corrections. A plausible implication is that eikonal methods continue to serve as a common language for comparing wavefront causality, amplitude analyticity, and classical observables across PDE theory, kinetic limits, and gravitational scattering.