Complex Caustics

• Complex caustics are singular structures defined by envelopes of rays and wave fronts, arising through multivalued evolution and catastrophe theory.
• They manifest in optics, fluid dynamics, quantum systems, and cosmology, with classifications such as folds, cusps, swallowtails, and umbilics.
• Analytical methods using complex scalar fields and geometric invariants enable the regularization of singularities and practical exploration of diverse physical phenomena.

A complex caustic is not a single object with one disciplinary definition, but a family of related singular structures associated with envelopes of rays, wave fronts, particle trajectories, or Lagrangian maps whose geometry, topology, or dynamics become nontrivial through multivalued evolution, catastrophe-theoretic classification, or explicit complexification. In current literature, the expression spans caustic singularities of pressureless perfect fluid and shift-symmetric k-essence completed by a canonical complex scalar field, caustics by reflection analyzed in complex projective geometry, many-body quantum caustics with Airy scaling, turbulent-aerosol and active-matter caustics, and the caustic skeleton of the cosmic web (Babichev et al., 2017, Uskova, 2 Jun 2026, Roy et al., 2024, Meibohm et al., 2022, Feldbrugge et al., 2017).

1. Terminological scope and core definitions

In geometrical optics and billiards, a caustic is an envelope. For an oval $\gamma \subset \mathbb{R}^2$ and a point $O$ inside, the $n$-th caustic by reflection, $\Gamma_n$, is the envelope of the family of rays from $O$ that have undergone $n$ reflections in $\gamma$ (Uskova, 2 Jun 2026). In the symplectic/contact formulation of optics, a caustic surface is the set of critical values of the projection from a Lagrangian or Legendrian submanifold to physical space,

$$\mathcal{C}=\left\{\mathbf{q}\mid \mathrm{rank}(d(\pi|_L))<2\right\},$$

or, for a parametrized family $\mathbf{q}(s,t)$, the locus where $$\|\partial_s\mathbf{q}\times \partial_t\mathbf{q}\|=0$$ (Shang et al., 3 Dec 2025).

In fluid, particle, and cosmological settings, caustics appear when a map from phase space or Lagrangian coordinates to configuration space becomes singular. For turbulent aerosols, caustics are folds in the spatial distribution of particles where the mapping from phase space to configuration space becomes singular, leading to multivalued particle velocities at a point; for the cosmic web, they are regions where density formally diverges due to shell-crossing and phase-space folding (Meibohm et al., 2022, Feldbrugge et al., 2017). In Bloch-sphere catastrophes, caustic surfaces appear when a real or complex field is mapped to its order parameter manifold, and the mapping becomes locally non-invertible (Eastwood et al., 2014).

The qualifier “complex” is used in several technically distinct senses. In elliptical billiards, the problem is complexified by treating the boundary as a complex algebraic curve in $O$0 and extending the reflection law with the complexified Euclidean metric (Uskova, 2 Jun 2026). In caustic-free completion of pressureless perfect fluid and k-essence, the relevant completion is a canonical complex scalar field, and the collapse time becomes complex (Babichev et al., 2017). In order-parameter mappings, the source object may itself be a complex field (Eastwood et al., 2014). This suggests that “complex caustic” is best understood as a cross-disciplinary label for caustics whose analysis requires complex fields, complex geometry, or singular multivalued structure.

2. Catastrophe-theoretic classification and invariant structure

A large part of the modern theory organizes caustics by catastrophe type. In three-dimensional optics, the stable singularities are those of codimension $O$1: fold $O$2, cusp $O$3, swallowtail $O$4, hyperbolic umbilic $O$5, and elliptic umbilic $O$6, with generating functions such as $O$7 for the fold and $O$8 for the cusp (Shang et al., 3 Dec 2025). In the billiard setting, an ordinary cusp is a singularity locally equivalent, holomorphically, to $O$9 (Uskova, 2 Jun 2026).

The cosmic-web literature adopts a parallel but dynamically specialized classification. In three dimensions the principal caustic types include $n$0 fold, $n$1 cusp, $n$2 swallowtail, $n$3 butterfly, and the $n$4-family $n$5, $n$6, and $n$7 (Feldbrugge et al., 2017). For Hamiltonian fluids, fold formation occurs when

$n$8

where $n$9 is an eigenvalue of the deformation tensor. The cusp condition adds

$\Gamma_n$0

and the swallowtail condition adds

$\Gamma_n$1

A more general shell-crossing criterion is

$\Gamma_n$2

which emphasizes that eigenvectors, not only eigenvalues, determine the spatial structure of the caustic skeleton (Feldbrugge et al., 2017).

In elliptical billiards, global algebro-geometric invariants play an analogous role. The proof that the first caustic $\Gamma_n$3 has exactly four ordinary cusps uses the genus $\Gamma_n$4, degree $\Gamma_n$5, dual degree $\Gamma_n$6, numbers of nodes $\Gamma_n$7, and cusps $\Gamma_n$8, related by Plücker formulas such as

$\Gamma_n$9

For $O$0, the computed data are $O$1, $O$2, and $O$3, and the only compatible solution gives exactly $O$4 ordinary cusps, with $O$5, $O$6, and $O$7 (Uskova, 2 Jun 2026).

3. Complexification, regularization, and singularity avoidance

One of the sharpest uses of complex structure is the caustic-free completion of pressureless perfect fluid and shift-symmetric k-essence by a canonical complex scalar field. The unifying two-scalar action is

$O$8

For $O$9 it reproduces pressureless perfect fluid, while for $n$0 it reduces to $n$1 with $n$2. After defining a complex scalar $n$3, the action becomes that of a canonical complex scalar field,

$n$4

The central mechanism is that the would-be collapse time is promoted to a complex value, so the singularity is not developed in real time (Babichev et al., 2017).

For the “perfect caustic” example of dust collapse, the pressureless-perfect-fluid velocity is

$n$5

with singularity at $n$6. In the complex completion,

$n$7

which yields

$n$8

Because the denominator never vanishes for real $n$9, $\gamma$0 and its derivatives stay finite. In the non-relativistic limit, the same completion reduces to the Schrödinger equation,

$\gamma$1

and the Madelung formulation introduces the quantum-pressure term

$\gamma$2

which smooths out would-be singularities (Babichev et al., 2017).

A different complex-analytic mechanism appears in the Helmholtz equation. In the Fock-Schwinger proper-time formulation,

$\gamma$3

and cusp caustics are tied to poles of the einbein action,

$\gamma$4

The residue $\gamma$5 vanishes on ghost sources, and cusp caustics lie at the intersection of a ghost source with a smooth caustic. Higher-order caustics are associated, by a proposed map, with essential singularities of the einbein action. The singularities originate from degenerations of a Dirichlet problem as the einbein is varied, and large classes of local perturbations do not move the poles or change their residues (Guralnik et al., 2019).

4. Reflection, refraction, and complex geometry of optical caustics

In elliptical billiards, the first caustic by reflection is studied in the complex projective plane. Reflection is defined as the unique nontrivial involutive complex isometry with respect to the complex quadratic form $\gamma$6 fixing a given non-isotropic line; isotropic lines require a limiting definition (Uskova, 2 Jun 2026). The complexification removes constraints of the real domain and permits generic-position arguments and the use of Plücker formulas. For an ellipse and a non-focal light source, the first reflected caustic $\gamma$7 is a rational curve of degree $\gamma$8 with exactly four ordinary cusps (Uskova, 2 Jun 2026).

The same general area contains a rigorous symplectic/contact description of optical caustics. The optical phase space is modeled as $\gamma$9 with canonical symplectic form

$$\mathcal{C}=\left\{\mathbf{q}\mid \mathrm{rank}(d(\pi|_L))<2\right\},$$0

while the extended phase space carries the contact form

$$\mathcal{C}=\left\{\mathbf{q}\mid \mathrm{rank}(d(\pi|_L))<2\right\},$$1

Light rays in three-dimensional Euclidean space correspond to Reeb orbits in a five-dimensional contact manifold, and wave fronts are Legendrian submanifolds satisfying $$\mathcal{C}=\left\{\mathbf{q}\mid \mathrm{rank}(d(\pi|_L))<2\right\},$$2. This framework classifies stable caustic surfaces and establishes a correspondence with Seidel aberration theory (Shang et al., 3 Dec 2025).

For caustics of wave fronts reflected by a surface, the reflected direction is

$$\mathcal{C}=\left\{\mathbf{q}\mid \mathrm{rank}(d(\pi|_L))<2\right\},$$3

The caustic is described as the focal set of a reflected wave front or, equivalently, the focal surface of a “virtual deformation” of the mirror with modified fundamental forms

$$\mathcal{C}=\left\{\mathbf{q}\mid \mathrm{rank}(d(\pi|_L))<2\right\},$$4

For flat incident wave fronts, the caustic is parametrized by

$$\mathcal{C}=\left\{\mathbf{q}\mid \mathrm{rank}(d(\pi|_L))<2\right\},$$5

where $$\mathcal{C}=\left\{\mathbf{q}\mid \mathrm{rank}(d(\pi|_L))<2\right\},$$6 are the roots of

$$\mathcal{C}=\left\{\mathbf{q}\mid \mathrm{rank}(d(\pi|_L))<2\right\},$$7

For spherical wave fronts, the same quadratic is shifted by $$\mathcal{C}=\left\{\mathbf{q}\mid \mathrm{rank}(d(\pi|_L))<2\right\},$$8 (Yampolsky et al., 25 Jan 2025).

For refraction, the complete caustic by refraction of a circle is the evolute of a Cartesian oval. With radiant point $$\mathcal{C}=\left\{\mathbf{q}\mid \mathrm{rank}(d(\pi|_L))<2\right\},$$9, circle center $\mathbf{q}(s,t)$0, radius $\mathbf{q}(s,t)$1, and an auxiliary point $\mathbf{q}(s,t)$2, the associated family of Cartesian ovals is

$\mathbf{q}(s,t)$3

As the circle tends to a line, the Cartesian oval collapses into an ellipse, so the caustic by refraction of a line is the evolute of an ellipse (Rydell, 2024).

5. Quantum, magnonic, and structured-light caustics

In many-body quantum dynamics, quantum caustics are extended regions of enhanced probability density where the classical ray description would diverge but wave interference smooths the singularity. In the transverse-field Ising model after a local quench, the outward-propagating excitation creates a light-cone-like caustic pattern dressed by Airy fringes. In the weak-coupling limit, the phase is

$\mathbf{q}(s,t)$4

and the continuum approximation near the light-cone edge gives

$\mathbf{q}(s,t)$5

The separation of the first two maxima scales as

$\mathbf{q}(s,t)$6

The exponent $\mathbf{q}(s,t)$7 is universal in the paramagnetic phase, starts varying at the quantum phase transition, and remains robust under weak integrability breaking and under open-boundary edge effects at the level of the overall power-law scaling (Roy et al., 2024).

In magnonics, caustics arise from anisotropic dispersion rather than from curved interfaces alone. In homogeneous in-plane magnetized thin films, the caustic condition is the vanishing of the curvature of the isofrequency contour,

$\mathbf{q}(s,t)$8

The near-field diffraction model computes the spin-wave response as

$\mathbf{q}(s,t)$9

In an extended $$\|\partial_s\mathbf{q}\times \partial_t\mathbf{q}\|=0$$0 thick yttrium iron garnet film, nano-constricted rf waveguides generate non-reciprocal caustic-like spin-wave beams, which were directly observed by spatially resolved micro-focused Brillouin light spectroscopy; the measurements agree with both micromagnetic simulation and the near-field diffraction model (Wagle et al., 2024).

Structured-light work extends caustics from canonical Airy or Pearcey profiles to engineered three-dimensional trajectories and morphing transverse patterns. A compensation phase implemented with 3D-printed metasurfaces permits caustic fields whose in-plane patterns are preserved or morphed from one structure to another during propagation. The field is designed in Fourier space, with focal-curve constraints

$$\|\partial_s\mathbf{q}\times \partial_t\mathbf{q}\|=0$$1

phase

$$\|\partial_s\mathbf{q}\times \partial_t\mathbf{q}\|=0$$2

and compensation phase

$$\|\partial_s\mathbf{q}\times \partial_t\mathbf{q}\|=0$$3

Large-scale fabrication is achieved by two-photon polymerization lithography (Zhou et al., 2023).

A related development in SU(2) structured beams shows that caustic-linked wave packets can visualize geometric phase directly, without interferometers or beam truncation. The beam is synthesized as

$$\|\partial_s\mathbf{q}\times \partial_t\mathbf{q}\|=0$$4

with Pancharatnam-Berry phase

$$\|\partial_s\mathbf{q}\times \partial_t\mathbf{q}\|=0$$5

total SU(2) Pancharatnam-Berry phase

$$\|\partial_s\mathbf{q}\times \partial_t\mathbf{q}\|=0$$6

and periodicity

$$\|\partial_s\mathbf{q}\times \partial_t\mathbf{q}\|=0$$7

The method visualizes both geometric phase and Gouy phase through the evolution of caustic-linked wave packets (Li et al., 17 Nov 2025).

6. Particles, cosmic structure, astrophysical events, and technological uses

For turbulent aerosols, caustics are dynamically important because they enhance collision rates and accelerate coagulation. In three dimensions, with particle velocity-gradient tensor $$\|\partial_s\mathbf{q}\times \partial_t\mathbf{q}\|=0$$8 and fluid velocity-gradient tensor $$\|\partial_s\mathbf{q}\times \partial_t\mathbf{q}\|=0$$9, the evolution law is

$O$00

and caustics form when $O$01 (Meibohm et al., 2022). In the weak-inertia persistent limit, threshold crossing in the $O$02-$O$03 plane is governed by the Vieillefosse line

$O$04

and the most likely route to caustic formation is a unique optimal fluctuation that propagates along its positive branch. The action-minimizing threshold configuration is pure strain,

$O$05

and the optimal trajectory exactly traces the positive branch of the Vieillefosse line (Meibohm et al., 2022).

The small- and large-inertia regimes differ qualitatively. For large particle inertia, caustic formation is analogous to Kramers escape in a white-noise model. For small inertia, it proceeds through an optimal fluctuation characterized by violent strain exceeding a threshold and low vorticity. In two dimensions, the fixed point ceases to exist when the Okubo-Weiss parameter exceeds

$O$06

and the optimal strain burst can be written as

$O$07

(Meibohm et al., 2020).

Active matter provides a zero-inertia analog. Self-propelled particles in vortical flow obey centroid and internal-polarization dynamics

$O$08

$O$09

Caustics arise without inertia because the active term $O$10 permits streamline crossing. Numerical studies in two-dimensional Navier-Stokes turbulence show intense caustics in straining regions, peaking at intermediate levels of self-propulsion (Chajwa et al., 2023).

In cosmology, the caustic skeleton formalism describes walls, filaments, and nodes as singularities of the Lagrangian map $O$11 with displacement field $O$12 and deformation tensor $O$13. The density is

$O$14

so caustics occur where $O$15 for some $O$16 (Feldbrugge et al., 2017). Applied to the Local Universe, the theory predicts two topologically distinct classes of filaments, $O$17 swallowtail and $O$18 umbilic. In the Manticore-Local reconstructions, the Pisces-Perseus Supercluster is revealed to be distinctly $O$19-dominated, while the extended Stickman structure around the Coma Cluster shows increasing relevance of $O$20 filaments only toward smaller scales (Read et al., 24 Apr 2026).

Astrophysical observations also use caustic language in lensing and cluster dynamics. OGLE-2016-BLG-1003 was the first resolved caustic-crossing binary-source event discovered by second-generation microlensing surveys. Its light curve showed two nested pairs of caustic-crossing features, and the decisive MOA coverage of two closely spaced caustic entrances uniquely supported the binary-source/binary-lens interpretation over the single-source/triple-lens alternative (1705.01531). In galaxy clusters, the caustic method uses the trumpet-shaped edge in projected radius–line-of-sight-velocity space to estimate the escape-velocity profile and mass profile. In mock redshift surveys with $O$21 galaxies within $O$22, it recovers $O$23–$O$24 of real substructures, while only $O$25–$O$26 of identified substructures correspond to real substructures of the central cluster (Yu et al., 2015).

Technological work increasingly treats caustics as engineered propagation resources. In wireless systems, representative caustic beams are classified by mathematical origin into Helmholtz eigenmode solutions, catastrophe-theory-based caustics, and generalized-Snell’s-law designs, and they are attributed three propagation properties: self-bending, self-healing, and near-field non-diffracting behavior (Liu et al., 10 Jun 2026). In dynamic display technology, ultrasonically modulated liquid surfaces driven by a $O$27 phased-array transducer at $O$28 generate dynamic caustic patterns. The acoustic pressure field is optimized through

$O$29

with numerical and Digital Twin losses

$O$30

and time averaging

$O$31

The reported system can generate continuous animations and complex caustic patterns at high frequencies, although with lower contrast and resolution than solid-surface methods (Nagakura et al., 22 May 2025).