pith. sign in

arxiv: 1907.10551 · v1 · pith:35JAHE2Inew · submitted 2019-07-24 · ⚛️ physics.optics

Morphological properties of 2D symmetric Airy beams extracted from the stationary wave approximation

Pith reviewed 2026-05-24 16:43 UTC · model grok-4.3

classification ⚛️ physics.optics
keywords Airy beamsstationary phase approximationphase singularitiesoptical vorticesnonparaxial regimecausticsstructured lightmorphological properties
0
0 comments X

The pith

Stationary phase approximation locates and classifies phase singularities in symmetric Airy beams.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines the structure of symmetric Airy beams using geometrical optics to find caustics, then applies the stationary phase approximation to the wave field to identify interference conditions at critical points. This approach works for both paraxial cases where roots are up to order 3 and nonparaxial where they reach order 6, including longitudinal field components. By doing so, it provides ways to count interfering waves and classify singularities such as optical vortices and dislocations, with confirmation from numerical simulations. A reader would care because it offers a general algorithm for analyzing any structured light field.

Core claim

The stationary wave approximation, built on critical points from caustics, yields conditions to identify the number of waves interfering constructively or destructively at key positions in 2D symmetric Airy beams. This enables distinguishing and classifying phase singularities like optical vortices and dislocations in both paraxial and nonparaxial regimes, as verified by numerical simulations.

What carries the argument

Stationary phase approximation applied to the light field at critical points identified by geometrical optics caustics.

Load-bearing premise

The stationary phase approximation remains accurate for locating and classifying morphological features even when the nonparaxial regime produces roots up to order 6 and when a longitudinal electric-field component is included.

What would settle it

Numerical integration of the nonparaxial Airy beam field including the longitudinal component showing phase singularity types or locations that differ from stationary phase predictions at the critical points.

Figures

Figures reproduced from arXiv: 1907.10551 by F. Camas-Aquino, P. A. Quinto-Su, R. J\'auregui.

Figure 2
Figure 2. Figure 2: FIG. 2: Intensity pattern and caustic (white dashed [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 1
Figure 1. Figure 1: FIG. 1: Intensity pattern and caustic (black dashed [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: (a,c) Caustic (black lines) and rays, (b,d) [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Results of applying the algorithm described in [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Results of applying the algorithm described in [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Intensity pattern and caustic (white dashed [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Results of applying the algorithm described in [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
read the original abstract

We explore the morphological properties of symmetric Airy beams in the paraxial and nonparaxial regimes. We consider a 2D electromagnetic realization with a single transverse component of the electric field, and in the nonparaxial regime, the longitudinal component along the optic axis. The general structure of these beams is analyzed with the combination of several approaches: geometrical optics through the use of caustics, the asymptotic wave properties of the light field using the stationary wave approximation and numerical integration. The geometrical optics approach involves locating the critical points that are later used in the stationary phase approximation. In the paraxial regime the highest order of the roots is 3, while in the nonparaxial regime, the order can be of up to 6. The technique yields conditions to identify interesting features on the beam, like the number of waves interfering constructively/destructively at the critical positions. The results are confirmed by the numerical simulations. In this way it is possible to distinguish and classify phase singularities like optical vortices and dislocations. The developed algorithm could be used to study any structured light field.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The manuscript analyzes morphological properties of 2D symmetric Airy beams in paraxial and nonparaxial regimes via geometrical optics (caustics to locate critical points), stationary-phase approximation for wave asymptotics, and numerical integration. Paraxial roots reach order 3 and nonparaxial up to order 6; the method supplies explicit conditions for constructive/destructive interference at those points and classifies phase singularities (vortices, dislocations). Results are stated to be confirmed by numerics, with the algorithm positioned as applicable to general structured light fields.

Significance. If the stationary-phase treatment remains quantitatively accurate for the higher-order degeneracies, the work supplies a concrete procedure for extracting interference counts and singularity types from beam profiles. The explicit numerical cross-checks and the inclusion of the longitudinal E_z component constitute positive elements that strengthen the central claim.

major comments (1)
  1. [nonparaxial regime analysis] Nonparaxial regime (roots of order up to 6): the stationary-phase conditions for locating interference features and classifying singularities are derived from the elementary stationary-phase formula. Standard stationary-phase asymptotics apply to isolated simple or low-order stationary points (Airy/Pearcey); an order-6 degeneracy requires a distinct higher-order canonical integral whose leading term differs. The manuscript does not indicate that uniform asymptotic expansions or higher-order corrections were employed when the longitudinal E_z component is restored. This directly affects the reported interference counts and singularity classifications, which are load-bearing for the central claim.
minor comments (1)
  1. The abstract states that numerical integration confirms the stationary-phase predictions but does not specify the grid resolution, error metrics, or whether the comparison was performed before or after parameter tuning.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback on our manuscript. We address the single major comment below and outline the revisions we will make.

read point-by-point responses
  1. Referee: Nonparaxial regime (roots of order up to 6): the stationary-phase conditions for locating interference features and classifying singularities are derived from the elementary stationary-phase formula. Standard stationary-phase asymptotics apply to isolated simple or low-order stationary points (Airy/Pearcey); an order-6 degeneracy requires a distinct higher-order canonical integral whose leading term differs. The manuscript does not indicate that uniform asymptotic expansions or higher-order corrections were employed when the longitudinal E_z component is restored. This directly affects the reported interference counts and singularity classifications, which are load-bearing for the central claim.

    Authors: We agree that the elementary stationary-phase formula for isolated simple stationary points does not directly furnish the leading asymptotic term for an order-6 degeneracy, which would instead require the appropriate higher-order canonical integral. In the manuscript the stationary-phase condition is used only to locate the critical points (vanishing first derivative of the total phase) and to obtain the relative phase values that distinguish constructive from destructive interference. The reported interference counts and the classification of phase singularities (vortices versus dislocations) are extracted from these phase relations and are then cross-validated by direct numerical integration of the full vector field, including the restored longitudinal component E_z. We acknowledge that the text does not explicitly discuss the limitations of the elementary formula for high-order points. In the revised manuscript we will insert a short clarifying paragraph in the methods section that states the scope of the stationary-phase approximation employed and notes that amplitude scaling at order-6 points is not claimed from the asymptotic expression but is instead confirmed numerically. This constitutes a partial revision that leaves the central claims unchanged while improving transparency. revision: partial

Circularity Check

0 steps flagged

No circularity: standard asymptotic methods plus independent numerical checks form self-contained chain

full rationale

The derivation locates critical points via caustics (geometrical optics), applies stationary-phase formulas to obtain interference and singularity conditions, and validates against separate numerical integration. No equations reduce a claimed prediction or classification to a fitted parameter or self-defined quantity from the same data; the stationary-phase step is the standard asymptotic expansion applied to externally located roots. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior author work appear in the provided text. The approach is externally falsifiable via the numerical benchmarks and uses only conventional methods.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard mathematical approximations in wave optics without introducing new fitted constants, postulated particles, or ad-hoc entities.

axioms (2)
  • domain assumption Stationary phase approximation accurately captures asymptotic interference behavior of the electromagnetic field near caustics
    Invoked to extract morphological properties from critical points located by geometrical optics
  • domain assumption Numerical integration of the wave equation provides an independent verification of the stationary-phase predictions
    Used to confirm results in both paraxial and nonparaxial regimes

pith-pipeline@v0.9.0 · 5735 in / 1434 out tokens · 21887 ms · 2026-05-24T16:43:36.874717+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

  1. [1]

    This is physically equivalent to select the optical aperture

    The paraxial condition can be mathematically implemented by the proper selection of the integration region. This is physically equivalent to select the optical aperture. An alternative is to choose an angular spectrum that is not negligible only on the appropriate wavevector region. The bidimensional finite energy Airy beam [19] is usu- ally defined by the ...

  2. [2]

    (5) Notice the presence of the parametera1 in the expressions of the overall phase

    The overall phase in k space for a finite energy Airy beam in the paraxial regime is Φpx F (ks1; s1, ζ) = kss3 + k3 s1 /3− a2 1ks1 + ks1 s1− k2 s1 ζ/2, (4) while out of this regime Φnpx F (ks1; s1, s3) = k3 s1 /3− a2 1ks1 + ks1 s1 + √ k2s− k2s1 s3. (5) Notice the presence of the parametera1 in the expressions of the overall phase. Similarly, the bidimensio...

  3. [3]

    Singular Optics,

    M. S. Soskin and M. V. Vasnetsov, “Singular Optics,” Progress in Optics 42 , 219276 (2001)

  4. [4]

    Dynamics of microparticles trapped in a per- fect vortex beam,

    M. Chen, M. Mazilu, Y. Arita, E. M. Wright, and K. Dholakia, “Dynamics of microparticles trapped in a per- fect vortex beam,” Opt. Letts. 38, 49194922 (2013)

  5. [5]

    Quantized Rotation of Atoms from Photons with Orbital Angular Momentum,

    M. F. Andersen, C. Ryu, P. Clad´ e, V. Natarajan, A. Vaziri, K. Helmerson, and W. D. Phillips, “Quantized Rotation of Atoms from Photons with Orbital Angular Momentum,” Phys. Rev. Letts. 97, 170406 (2006)

  6. [6]

    Breaking Abbe’s diffraction resolution limit in fluorescence microscopy with stimulated emission depletion beams of various shapes,

    T. A. Klar, E. Engel, and S. W. Hell, “Breaking Abbe’s diffraction resolution limit in fluorescence microscopy with stimulated emission depletion beams of various shapes,” Phys. Rev. E 64, 066613 (2001)

  7. [7]

    Evolution of quan- tum superoscillations and optical superresolution with- out evanescent waves,

    M. V. Berry and S. Popescu, “Evolution of quan- tum superoscillations and optical superresolution with- out evanescent waves,” J. Phys. A: Math. Gen. 39, 69656977 (2006)

  8. [8]

    Deep microstructure topography characterization with optical vortex interferometer,

    J. Masajada, M. Leniec, E. Jankowska, H. Thienpont, H. Ottevaere, V. Gomez, “Deep microstructure topography characterization with optical vortex interferometer,”Opt. Express 16, 1917919191 (2008)

  9. [9]

    Rotational effects of twisted light on atoms beyond the paraxial approximation,

    R. J´ auregui, “Rotational effects of twisted light on atoms beyond the paraxial approximation,” Phys. Rev. A 70, 033415 (2004)

  10. [10]

    Control of atomic transition rates via laser- light shaping,

    R. J´ auregui, “Control of atomic transition rates via laser- light shaping,” Phys. Rev. A 91, 043842 (2015)

  11. [11]

    Transfer of optical orbital angular momentum to a bound electron,

    C. T. Schmiegelow, J. Schulz, H. Kaufmann, T. Ruster, U. G. Poschinger, and F. Schmidt-Kaler, “Transfer of optical orbital angular momentum to a bound electron,” Nat. Commun. 7, 12998 (2016)

  12. [12]

    Electromagnetic Field Enhancement and Spec- trum Shaping through Plasmonically Integrated Optical Vortices,

    W. Ahn, S. V. Boriskina, Y. Hong, and B. M. Rein- hard, “Electromagnetic Field Enhancement and Spec- trum Shaping through Plasmonically Integrated Optical Vortices,” Nano Lett. 12 219227 (2012)

  13. [13]

    Achieving an Accurate Surface Profile of a Photonic Crystal for Near-Unity Solar Absorption in a Super Thin-Film Architecture,

    P. Kuang, S. Eyderman, M. Hsieh, A. Post, S. John, and S. Lin, “Achieving an Accurate Surface Profile of a Photonic Crystal for Near-Unity Solar Absorption in a Super Thin-Film Architecture,” ACS Nano 10 61166124 (2016)

  14. [14]

    Plasmonic materials for energy: From physics to applications,

    S. V. Boriskina, H. Ghasemi and G. Chen, “Plasmonic materials for energy: From physics to applications,” Mater. Today 16, 375386 (2013)

  15. [15]

    Super-resolution optical telescopes with local light diffraction shrinkage,

    C. Wang, D. Tang, Y. Wang, Z. Zhao, J. Wang, M. Pu, Y. Zhang, W. Yan, P. Gao, and X. Luo, “Super-resolution optical telescopes with local light diffraction shrinkage,” Sci. Rep. 5, 18485 (2015)

  16. [16]

    Symmetric Airy beams,

    P. Vaveliuk, A. Lencina, J. A. Rodrigo, and O. Mar- tinez Matos, “Symmetric Airy beams,” Opt. Lett. 39, 23702373 (2014)

  17. [17]

    On the general prop- erties of symmetric incomplete Airy beams,

    R. J´ auregui and P. A. Quinto-Su, “On the general prop- erties of symmetric incomplete Airy beams,” J. Opt. Soc. Am. A 31, 24842488 (2014)

  18. [18]

    Optical stacking of microparticles in a pyramidal structure created with a symmetric cubic phase,

    P. A. Quinto-Su and R. J´ auregui, “Optical stacking of microparticles in a pyramidal structure created with a symmetric cubic phase,” Opt. Express 22, 1228312288 (2014)

  19. [19]

    Caustics, catastrophes, and symmetries in curved beams,

    P. Vaveliuk, A. Lencina, J. A. Rodrigo, and O. Mar- tinez Matos, “Caustics, catastrophes, and symmetries in curved beams,” Phys. Rev. A 92, 033850 (2015)

  20. [20]

    Nonspreading wave pack- ets,

    M. V. Berry and N. L. Balazs, “Nonspreading wave pack- ets,” Am. J. Phys. 47, 264267 (1979)

  21. [21]

    Acceler- ating finite energy Airy beams,

    G. A. Siviloglou and D. N. Christodoulides, “Acceler- ating finite energy Airy beams,” Opt. Lett. 32, 979981 (2007)

  22. [22]

    Observation of Accelerating Airy Beams,

    G. A. Siviloglou, J. Broky, A. Dogariu and D. N. Christodoulides, “Observation of Accelerating Airy Beams,” Phys. Rev. Lett. 99, 213901 (2007)

  23. [23]

    Catastrophe Optics: Mor- phologies of caustics and their diffraction patterns,

    M. V. Berry and C. Upstill, “Catastrophe Optics: Mor- phologies of caustics and their diffraction patterns,” Progress in Optics XVIII , 257323 (1980)

  24. [24]

    V. I. Arnold, Singularities of Caustics and Wave Fronts , Springer, Mathematics and Its Applications (Soviet Se- ries) Volume 62 (1990)

  25. [25]

    J. F. Nye, Natural Focusing and Fine Structure of Light: Caustics and Wave Dislocations , Taylor & Francis (1999)

  26. [26]

    V. I. Arnold, Catastrophe Theory , Springer, 3rd ed. (1992). 10

  27. [27]

    The structure of an electromagnetic field in the neighbourhood of a cusp of a caustic,

    T. Pearcey, “The structure of an electromagnetic field in the neighbourhood of a cusp of a caustic,” Phil. Mag. 37, 311317 (1946)

  28. [28]

    Zur Darstellung von Wellenfeldern

    W. Braunbek, “Zur Darstellung von Wellenfeldern”, Z. Naturforsch A 6, 12-15 (1951)

  29. [29]

    Topology of optical vortex lines formed by the inter- ference of three, four, and five plane waves

    K. O’Holleran, M. J. Padgett and M. R. Den- nis,“Topology of optical vortex lines formed by the inter- ference of three, four, and five plane waves”,Opt. Express 14, 3039 (2006)

  30. [30]

    Self-healing properties of optical Airy beams,

    J. Broky, G. A. Siviloglou, A. Dogariu, and D. N. Christodoulides, “Self-healing properties of optical Airy beams,” Opt. Express 16, 1288012891 (2008)

  31. [31]

    Airy beams: a geometric optics per- spective,

    S. Vo, K. Fuerschbach, K. P. Thompson, M. A. Alonso, and J. P. Rolland, “Airy beams: a geometric optics per- spective,” J. Opt. Soc. Am. A 27, 25742582 (2010)

  32. [32]

    Auto-focusing and self-healing of Pearcey beams,

    J. D. Ring, J. Lindberg, A. Mourka, M. Mazilu, K. Dho- lakia, and M. R. Dennis, “Auto-focusing and self-healing of Pearcey beams,” Opt. Express 20, 1895518966 (2012)

  33. [33]

    V. P. Maslov, Theoria vozmushchenii i asimptoticheskie metody (Theory of Perturbations and Asymptotic Meth- ods) (Izd. MGU, 1965) [French trans.: V. P. Maslov,The- orie des perturbations et methode asymptotique (Duod, 1972)]

  34. [34]

    Asymptotic theory of wave-propagation

    R. M. Lewis, “Asymptotic theory of wave-propagation”, Arch. Ration. Mech. Anal. 20, 191 (1965)

  35. [35]

    Considerations on the electromagnetic flow in Airy beams based on the Gouy phase,

    C. J. Zapata-Rodr´ ıguez, D. Pastor, and J. J. Miret, “Considerations on the electromagnetic flow in Airy beams based on the Gouy phase, ” Opt. Express 20, 2355323559 (2012)

  36. [36]

    Evolution from a Fraunhofer to a Pearcey diffraction pattern,

    J. F. Nye, “Evolution from a Fraunhofer to a Pearcey diffraction pattern,”J. Opt. A: Pure Appl. Opt. 5, 495502 (2003)

  37. [37]

    Optics of nonuniformly moving media,

    U. Leonhardt and P. Piwnicki, “Optics of nonuniformly moving media,” Phys. Rev. A 60, 4301 (1999)

  38. [38]

    GoosH¨ anchen and Im- bertFedorov shifts for bounded wavepackets of light,

    M. Ornigotti and A. Aiello, “ GoosH¨ anchen and Im- bertFedorov shifts for bounded wavepackets of light,” J. Opt. 15, 014004 (2013)