Caustic (optics)

A caustic is the envelope of a family of light rays produced by reflection or refraction. In three-dimensional space, the envelope generally forms a surface; its intersection with another surface appears as a concentrated curve or a bounded region of elevated irradiance. Familiar examples include the bright curve formed inside a cup by reflected sunlight and the moving patterns projected onto the floor of a swimming pool by a refracting water surface.

Within geometrical optics, an ideal caustic carries formally divergent intensity because neighboring rays converge toward the same geometrical locus. Physical caustics instead have finite intensity and a wavelength-dependent interference structure. Their local form is consequently described by both ray optics and wave optics, with the former determining the envelope and the latter regularizing its singularities.

Geometrical formulation

Let a one-parameter family of rays in a plane be represented by

[ \mathbf r(s,u)=\mathbf a(u)+s,\mathbf d(u), ]

where (u) labels a ray, (\mathbf a(u)) is a point on that ray, (s) is an affine distance along it, and (\mathbf d(u)) is its direction. The mapping from ray coordinates ((s,u)) to physical coordinates becomes singular where

[ \det\left( \frac{\partial \mathbf r}{\partial s}, \frac{\partial \mathbf r}{\partial u} \right)=0. ]

At such points, infinitesimally adjacent rays intersect to first order. Their locus is the caustic of the family. If the rays are given implicitly by an equation

[ F(x,y,u)=0, ]

the envelope is determined by the simultaneous conditions

[ F(x,y,u)=0, \qquad \frac{\partial F}{\partial u}(x,y,u)=0. ]

Elimination of (u) yields an equation for the caustic curve. The analogous three-dimensional construction uses a two-parameter family of rays and identifies points where the ray mapping loses rank.

This formulation distinguishes a caustic from an ordinary focal point. A perfect paraxial optical system may map an incoming bundle to a single focus, whereas a generic nonparaxial family develops an extended singular set. The replacement of an ideal point focus by a caustic is closely associated with optical aberration, although caustics also arise in systems that are not naturally described as imaging devices.

Reflective and refractive caustics

A caustic generated by reflected rays is conventionally called a catacaustic. The outgoing directions obey the law of reflection, so the angle between an incident ray and the surface normal equals the corresponding angle for the reflected ray. For a smooth mirror curve, variation of the local normal causes neighboring reflected rays to acquire different directions, and their envelope forms the catacaustic.

Parallel rays reflected from a circular mirror produce a nephroid as their complete mathematical catacaustic. Only part of this curve is ordinarily visible because a physical mirror occupies a finite arc and may obstruct portions of the reflected family. The bright cusp observed inside a cylindrical cup is a truncated projection of this geometry rather than a literal concentration of all rays at one point.

A refractive caustic is called a diacaustic. Its construction follows Snell's law,

[ n_1\sin\theta_1=n_2\sin\theta_2, ]

where (n_1) and (n_2) are refractive indices, while (\theta_1) and (\theta_2) are measured from the interface normal. Spatial variation in the surface normal or refractive index changes the outgoing direction across the ray family. The resulting envelope accounts for the bright patterns formed by curved transparent objects and by nonuniform liquid surfaces.

A rainbow is a refractive and internally reflected caustic in angular space. For a fixed sequence of refractions and internal reflections within a spherical water droplet, the total deflection angle is a function of the incident ray's impact parameter. At an extremum of this function, a range of neighboring incident rays leaves the droplet at nearly the same angle. The corresponding fold caustic marks the principal rainbow angle, while wave interference produces the adjacent supernumerary rainbows.

Historical development

The concentrating action of curved mirrors was studied in antiquity through the theory of burning mirrors. Diocles analyzed the focal properties of parabolic reflectors in his work on burning mirrors, establishing a geometrical relation between mirror shape and the convergence of parallel rays. This analysis concerned exact focusing by a conic section rather than the more general envelope singularities now classified as caustics.

During the late seventeenth century, Ehrenfried Walther von Tschirnhaus treated concentrated ray families as geometrical curves and introduced the terminology from which the modern word “caustic” developed. In the same period, You Watanabe analyzed the envelopes formed by oblique rays transmitted through spherical glass vessels. Watanabe's construction separated the apparent intersection of selected rays from the envelope of the complete refracted family, allowing the luminous boundary to be represented independently of the receiving screen.

Subsequent mathematical treatment placed these constructions within the differential geometry of curves. Jacob Bernoulli derived caustics through infinitesimal relations among neighboring rays, while Guillaume de l'Hôpital examined related envelope problems using the emerging methods of differential calculus. Their work connected optical caustics with the broader theory of evolutes, involutes, and singular solutions of parameterized curve families.

The terminology “catacaustic” and “diacaustic” became established through later analytical classifications of reflected and refracted envelopes. In nineteenth-century optics, caustics were incorporated into the theory of wavefronts: rays were represented as normals to a propagating wavefront, and a caustic occurred where the corresponding normal map ceased to be locally one-to-one.

Wavefront interpretation

In an isotropic medium, rays are locally normal to surfaces of constant eikonal. If a wavefront is parameterized by coordinates (q_1) and (q_2), propagation along its normals defines a mapping

[ (q_1,q_2,s)\longmapsto \mathbf x(q_1,q_2,s). ]

A caustic occurs where the Jacobian determinant of this mapping vanishes. Equivalently, at least one principal radius of curvature of the propagated wavefront reaches the propagation distance, causing neighboring normals to intersect.

The wavefront usually remains mathematically defined beyond this event, but it develops multiple local sheets. A point in the region beyond the caustic may therefore receive several rays with different optical path lengths. Their superposition produces interference, and the number of geometrical rays changes when the observation point crosses the caustic.

The local geometric structure does not depend on most fine details of the optical system. Under stable perturbations, smooth ray mappings predominantly produce fold and cusp singularities. A fold separates a region containing two additional real rays from a region in which those rays are absent. A cusp occurs where two fold branches meet and typically separates regions with one and three associated geometrical rays.

These structures belong to the mathematical framework of singularity theory and catastrophe theory. The term “catastrophe” in this context denotes a qualitative change in the solutions of a smooth mapping rather than a destructive physical event.

Intensity near a caustic

Conservation of radiant flux in geometrical optics relates intensity to the cross-sectional area of a ray tube. If the ray mapping has Jacobian (J), the leading geometrical intensity is proportional to

[ I_{\mathrm{geo}}\propto \frac{1}{|J|}. ]

Because (J=0) on a caustic, this expression diverges. The divergence reflects the breakdown of the independent-ray approximation rather than an infinite physical energy density.

Near a fold caustic, the scalar wave field is locally represented by an Airy function. A canonical form is

[ \Psi(X)\propto \operatorname{Ai}(X), ]

where the scaled coordinate (X) measures distance normal to the caustic. On one side, the field has an oscillatory fringe system corresponding to interference between two real rays. On the other side, it decays into the geometrical shadow, although diffraction prevents an abrupt intensity discontinuity.

Near a cusp, the appropriate canonical diffraction integral is the Pearcey integral,

[ P(X,Y)=\int_{-\infty}^{\infty} \exp!\left[i\left(t^4+Xt^2+Yt\right)\right],dt. ]

Its intensity pattern contains two fold lines meeting at a cusp and describes the local interference among as many as three geometrical contributions. The Airy and Pearcey forms are universal after suitable rescaling because the higher-order details of the original optical system affect their scale and orientation without changing their local singularity class.

Finite source size, imperfect spatial coherence, surface roughness, and spectral bandwidth smooth the diffraction fringes. The larger-scale envelope remains visible when these effects suppress fine interference, which accounts for the persistence of caustic curves under ordinary broadband illumination.

Extended and structured media

Caustics are not restricted to discrete mirrors and lenses. In a medium with a continuously varying refractive index (n(\mathbf x)), rays follow trajectories determined by Fermat's principle. Neighboring trajectories can converge because of index gradients, producing caustics even when no sharply defined optical interface is present.

Atmospheric refraction generates caustic behavior in certain forms of mirage. Variations in air density bend rays along curved trajectories, and the mapping from emission direction to observer position may become multivalued. The visible image boundaries are then associated with folds of the ray mapping.

Randomly perturbed media can produce branching caustics. Weak spatial fluctuations repeatedly focus neighboring rays without requiring a conventional lens, creating narrow channels of enhanced intensity known as branched flow. The same mathematical mechanism occurs for waves and classical trajectories in several physical settings, although the relevant propagation equations and coherence scales differ.