Airy-Gaussian Geometry: Hybrid Optics & Stochastics

• Airy-Gaussian geometry is a hybrid paradigm blending non-diffracting, self-accelerating Airy beams with Gaussian profiles, exhibiting locally Brownian fluctuations.
• The analytic framework uses paraxial wave and Schrödinger-type equations to describe reversible beam transformations, caustics, and envelope dynamics.
• Practical applications include engineered wavefront shaping, robust terahertz communication, and nonlinear vortex beam dynamics in integrated photonics.

Airy-Gaussian Geometry denotes both a class of hybrid wave phenomena—blending the quintessential features of Airy and Gaussian beams or processes—and the analytic and geometric structures that arise in such hybridizations. The term naturally encompasses the optical context, where beams interpolate or transition between Gaussian and Airy modalities, as well as the stochastic/probabilistic context, where spatial fluctuations governed by Airy processes exhibit locally Gaussian behavior. These dual aspects, unified by the interplay of non-diffraction, self-bending, and probabilistic Brownian properties, underpin a rich landscape spanning mathematical physics, nonlinear optics, statistical mechanics, and information theory.

1. Fundamental Definitions and Prototypical Structures

"Airy–Gaussian geometry" comprises continuous field or process configurations that locally interpolate or mediate between Gaussian (Brownian) and Airy-type behavior.

• Airy Beams and Processes: Airy beams are non-diffracting, self-accelerating solutions to the paraxial wave equation, with spatial intensity profiles proportional to the Airy function $\mathrm{Ai}(\xi)$ and characteristic parabolic trajectories. Airy processes (e.g., Airy$_1$, Airy$_2$) arise in Last-Passage Percolation (LPP) and KPZ-related models as scaling limits of surface fluctuation fields (Pimentel, 2017).
• Gaussian Beams and Processes: Gaussian beams represent the diffraction-limited, minimum-uncertainty solution to the paraxial equation, with symmetric intensity profiles and no acceleration or caustic singularity. In growth models, Gaussian processes (Brownian motion) govern local increments under appropriate scaling limits.

A canonical Airy-Gaussian beam in optics takes the tensor-product or envelope-modulated form: $$E(x,y,0) = e^{-[g_x x^2 + g_y y^2]}\, \mathrm{Ai}[\eta(x+i y)] \qquad (g_x, g_y > 0)$$ Hybrid beams such as the circular Airy-Gaussian vortex (CAGV) and asymmetric Cauchy–Riemann beams generalize this template by introducing vortex phase factors and asymmetric Gaussian envelopes (Hu et al., 2022, Korneev et al., 2024).

In stochastic geometry, the “Airy–Gaussian geometry” refers to the rigorous result that small increments of the Airy process or the Airy sheet are locally Gaussian (Brownian) under the natural scaling (Pimentel, 2017): $A(x+h) - A(x)\ \approx\ \sqrt{2h}\, N(0,1)$ where $A(\cdot)$ denotes the Airy process.

2. Governing Equations and Analytic Frameworks

Optical Airy–Gaussian Beams

In paraxial optics, Airy–Gaussian hybrid beams satisfy the paraxial wave equation: $i\,\partial_z E + \frac{1}{2k}\nabla_\perp^2 E = 0$ subject to initial conditions that blend the Airy and Gaussian forms. For instance, an asymmetric Airy–Gaussian beam propagates according to the operator method: $$E(x, y, z) = \frac{\exp\left[ -\left(\frac{g_x x^2}{w_x(z)} + \frac{g_y y^2}{w_y(z)} \right) \right]}{\sqrt{w_x(z) w_y(z)}}\, \exp[ \cdots ]\, \mathrm{Ai}(\zeta(x, y, z))$$ where the $w_q(z)$, $\Delta(z)$, and $\zeta(x, y, z)$ encode the envelope’s breathing and caustic deformation (Korneev et al., 2024).

In engineered waveguides with a linear gradient-index profile mimicking a gravitational field, the transition between Gaussian and Airy geometries is governed by a Schrödinger-type equation: $$i \partial_z \psi = -\frac{1}{2 n_0 k_0} \partial_x^2 \psi - k_0 [a x + (b - n_0)] \psi$$ This admits Airy-function solutions in the linear region and enacts a reversible Gaussian $\rightarrow$ Airy $\rightarrow$ Gaussian transformation on the beam cross-section (Wang et al., 2018).

Stochastic Airy Processes

In spatial stochastic models, the Airy process $A^b(u)$ arises as the scaling limit of the rescaled LPP weights: $$A_n^b(u) = 2^{-1/3} n^{-1/3} \left\{L_b([u n^{2/3}]) - L_b([0])\right\}$$ with $L_b$ denoting boundary-perturbed last-passage weights. The local Gaussianity (Brownian scaling) and continuity of $A^b$ are established via coupling techniques and tightness in Skorohod space (Pimentel, 2017).

3. Geometric Properties and Parameter Dependencies

Airy–Gaussian fields are characterized by their unique beam and fluctuation geometries, controlled by envelope, truncation, and scaling parameters.

• Trajectory and Caustics: Airy–Gaussian beams exhibit main lobe trajectories described by generalized parabolas, modified by asymmetric envelope widths ($g_x \neq g_y$) or vortex phase parameters (for CAGV). The caustic surface in circular Airy–Gaussian vortex beams (CAGVBs) follows:

$$r_\mathrm{caustic}(z) = r_0 - \frac{z^2}{4 \omega^2}$$

Abrupt autofocusing occurs at a finite $z_f$ (Hu et al., 2022).

• Width and Asymmetry: The Gaussian parameters $g_x$, $g_y$ (or their equivalents $w$, $\omega$) govern the transverse profile squeezing and rate of envelope “breathing.” For $g_x \neq g_y$ the focal spot is stretched into an asymmetric caustic cross-section with drift toward the lower $g_q$ direction (Korneev et al., 2024).
• Vortex Charge and Nonlinearity: In nonlinear media, the CAGVB’s focal ring and distance scale with fractional Laplacian exponent $\alpha$, distribution parameter $b$, topological charge $m$, and input power $P_0$ as:

$$R_{\nabla} \propto \alpha^{0.3} b^{0.5} (m+1)^{0.4} P_0^{-0.2},\quad z_\text{focus} \propto \alpha^{-0.8} b^{0.6} (m+1)^{0.1} P_0^{-0.5}$$

(He et al., 2020).

• Stochastic Increments: For Airy processes $A(x)$ or sheets $A(x, y)$, local increments are asymptotically Gaussian:

$$A(x+h) - A(x) \approx \sqrt{2h} N(0, 1),\quad A(x+h, y+k) - A(x, y) \approx \sqrt{2}[B_1(h) + B_2(k)]$$

(Pimentel, 2017).

4. Optical Realizations, Transitions, and Applications

Engineered Transitions

• Reversible Wavefront Shaping: By constructing a gradient-index waveguide whose effective index varies linearly with $x$, one realizes a linear potential in the paraxial equation, enabling a Gaussian input beam to break up into an accelerating Airy lobe and then revert to the Gaussian form. Experiment and full-wave simulations confirm the exact theoretical Airy–Gaussian cycle; the wavefront is shaped reversibly over a $z \lesssim 100$ μm region (Wang et al., 2018).
• Vector and Vortex Beams: Circular Airy–Gaussian vortex beams (CAGV, CAGVB) encode both OAM (vortex) and accelerating caustic behavior, with the possibility of abrupt autofocusing or autodefocusing, tunable via launch angle and envelope parameters. Vector modes—nonseparable in polarization and spatial degrees—display spatially evolving polarization and may exhibit phase rotation linked to OAM difference (Hu et al., 2022, He et al., 2020).
• Terahertz Communication: Airy beams and Airy–Gaussian hybrids have been proposed for robust quasi-line-of-sight Terahertz (THz) wireless links in data centers, outperforming Gaussian beams in propagation tolerance, self-bending around obstructions, and side-lobe-mediated self-healing. The cascaded geometrical-wave-based channel model (CGWCM) captures these unique propagation effects (Zhao et al., 29 Apr 2025).

Parameter Effects in Nonlinear Regimes

• Fractional Diffraction: In media governed by the fractional nonlinear Schrödinger equation (FNSE), the abruptness, radius, and distance to autofocusing are controlled by the Lévy index $\alpha$, envelope width, vortex charge, and input intensity, resulting in rich autofocusing and autodefocusing dynamics not present in standard (integer exponent) cases (He et al., 2020).
• Off-Axis Dynamics: Off-axis multiple CAGVBs with positive vortex pairs exhibit spiral convergence (autofocusing), shape deformation, and eventual outward divergence due to the interplay of nonlinearity, diffraction, and vortex interactions (He et al., 2020).

5. Mathematical Characterization and Local Tangent Geometry

The most rigorous formulation of Airy–Gaussian geometry is through the modulus-of-continuity and local limit theorems for rescaled spatial processes of Airy type.

• Continuity and Local Gaussianity: For any subsequential limit of last-passage percolation with a suitable profile, the Airy process $A^b(\cdot)$ satisfies a Brownian modulus-of-continuity and, after rescaling, converges locally to standard Brownian motion:

$$|A^b(x)-A^b(y)| = O\left(|x-y|^{1/2} \sqrt{\ln \frac{1}{|x-y|}}\right) \text{ a.s.}$$

and

$$\epsilon^{-1/2} (A^b(\epsilon x) - A^b(0)) \Rightarrow B(x)$$

for standard Brownian $B$, as $\epsilon \to 0$ (Pimentel, 2017).

• Airy Sheet: The two-parameter Airy sheet $A(x, y)$ exhibits increments that, post-rescaling, become sums of independent Brownian motions in each coordinate direction. Explicitly:

$$A(x+h, y+k) - A(x, y) \approx \sqrt{2}[B_1(h) + B_2(k)]$$

(Pimentel, 2017).

This implies that at infinitesimal scales, all Airy fluctuation fields in the KPZ universality class possess a “tangent” geometry equivalent to a Gaussian field. The explicit correspondence enables detailed analysis of fluctuations, coalescence behavior, and regularity properties.

6. Comparison Table: Airy, Gaussian, and Airy–Gaussian Features

Feature	Gaussian	Airy	Airy–Gaussian (Hybrid)
Diffraction/Spreading	Yes; $w(z) \propto z$	Non-diffracting main lobe	Intermediate/parameter-tunable
Trajectory	Straight	Parabolic, self-accelerating	Parabola/ellipse, possibly asymmetric
Side Lobes	None	Decaying oscillatory tail	Present, tunable by truncation
Self-healing	No	Yes (side-lobe mediated)	Yes/Partial (depends on truncation)
Autofocusing	No	Yes (caustic collapse)	Tunable abruptness/focal location
Local Stochastic Fluct.	Gaussian	Locally Gaussian (Brownian limit)	Brownian at small scales
Vortex/Polarization	Optional (LG beams)	Yes (vortex Airy), OAM possible	Yes, full vector mode families possible

7. Outlook and Ongoing Research

The Airy–Gaussian geometry framework is central to contemporary explorations in beam shaping, nonparaxial wave control, and stochastic growth. Emerging directions include:

• Engineering novel hybrid beams for integrated photonics and high-precision particle manipulation using curvature-controlled Airy–Gaussian profiles (Wang et al., 2018, Zhao et al., 29 Apr 2025).
• Detailed study of vortex and vectorial Airy–Gaussian beams in nonlinear and fractional-diffraction regimes, with focus on tunable autofocusing, autodefocusing, and topological effects (Hu et al., 2022, He et al., 2020).
• Further mathematical elucidation of Airy–type processes, sheets, and higher-dimensional analogues, leveraging their precise local Gaussian structure for universality classification and fluctuation estimates in disordered systems (Pimentel, 2017).

The interplay of analytic tractability, experimental realizability, and probabilistic universality continues to make Airy-Gaussian geometry a point of convergence for diverse theoretical and applied fields.