# Optical catastrophes of the swallowtail and butterfly beams

###### Abstract

We experimentally realize higher-order catastrophic structures in light fields to access the rich class of caustic swallowtail and butterfly beams. These beams present solutions of paraxial diffraction catastrophe integrals that are determined by potential functions, whose singular mapping manifests as caustic hypersurfaces in control parameter space. We systematically analyze the swallowtail and butterfly beams’ caustics analytically and observe their field distributions experimentally in real and Fourier space. Their spectra can be described by polynomials or expressions with rational exponents capable to form a cusp.

^{†}

^{†}affiliationtext: Institute of Applied Physics and Center for Nonlinear Science (CeNoS), University of Münster, 48149 Münster, Germany

^{†}

^{†}affiliationtext: a.zannotti@uni-muenster.de

Keywords: Caustic, catastrophe theory, nonlinear system, swallowtail, butterfly, paraxial beams.

## 1 Introduction

Caustics are ubiquitous in nature and our daily life, motivating scientists to investigate and utilize this general phenomenon in versatile fields, reaching from astrophysics [1, 2, 3] to models in social sciences [4, 5, 6, 7], and oceanography [8, 9, 10]. Also, they are present in imaging of any kind of imperfect lens as well as in microscopic imaging [11, 12, 13, 14, 15]. For example, caustics play an important role in astrophysics, as the probability density function of particles (light, matter) distributed in a caustic network bears similarities to mass distribution in galaxies [1, 3] and may be related to predictions of the millennium simulation. Thus, understanding how these caustic networks form, their relation to extreme wave effects in the framework of freak and Rogue waves, as they occur, for instance, in oceans, is of high interest for a broad community of scientists in many fields of research [16, 17].

In particular, the occurrence of caustics in light structures presents one of the oldest and most fundamental phenomena in optics [18]. These incoherent ray-optic effects [19, 20] appear as high-intensity lines or surfaces. They are associated with supernumerary arcs close to rainbows [19] or may occur as bright lines of high intensity at the floor of shallow waters. Similar to their formation behind refractive lenses with imperfections, corresponding effects have been observed for numerous kinds of lenses with importance in optics, and surface analytics, in particular gravitational [21, 22] and electro-magnetic [13, 14, 15] lenses. Overall, many fields will benefit from finding approaches to design high-contrast computational caustic mappings [23, 24, 25].

Caustics in light are the manifestation of so-called catastrophes: Singularities in the gradient map of nonlinear potential functions that form as geometrically stable structures when perturbing external control parameters. The most prominent representative of catastrophe [26] that manifests as caustic wave package is the fold catastrophe that was shown to be a solution of the quantum mechanical Schrödinger equation in form of the Airy function [20], and has been transferred and established in optics in 2007 as the paraxial Airy beam. One of the Airy beam’s most exotic properties is the transverse invariant, accelerated movement of its caustic on a parabolic trajectory during propagation, yielding a wealth of applications [27, 28, 29]. The Airy beam thereby represents a spatial light structure besides many new light fields that have been explored as solutions of the paraxial Helmholtz equation in recent years capable to open new avenues in photonics. Among the most famous solutions we find generalizations of Gaussian beams [30, 31], as well as nondiffracting beams [32], like for instance Bessel beams [33, 34], as well as Mathieu and Weber beams [35], that were also designed to be self-accelerated [36]. However, the introduction of the fold catastrophe as paraxial Airy beam acts as a prelude for the renaissance to transfer higher-order catastrophes to optics. Subsequent to the Airy beam, in 2012 the Pearcey beam [37], which represents the diffraction of a cusp catastrophe, emerges as a paraxial beam with intriguing features, like its form-invariance and auto-focusing propagation.

In general, the dimensionality of the investigated catastrophe causes caustics to occur as points, lines, surfaces or hypersurfaces, whose structure manifests as singular mapping of corresponding nonlinear potentials. Catastrophe theory predicts and analyses these potential functions and determines caustics as degenerated critical positions which show abrupt transitions of the observed system [26]. Suggested by Arnol’d, a classification of seven elementary catastrophes determined by their potential functions has been theoretically established, corresponding to seven generic structures of bifurcation geometries [38]. These are, increasing in the order of catastrophe, given by: fold, cusp, swallowtail, butterfly, elliptical umbilic, hyperbolic umbilic, and parabolic umbilic [39, 40, 41]. Although all of these catastrophes have been transferred by theoretical considerations to optical light structures [42, 41, 43], only the creation the lowest-order catastrophes as Airy (fold) and Pearcey (cusp) beams have been experimentally realized in the paraxial approximation and utilized for applications. However, controlling caustics in all its parameters is essential for advanced, e. g. superresolution imaging and precision material processing technologies. Until now, they have not been artificially created beyond their occurrence in natural light fields.

Thus, with this work we present higher-order swallowtail and butterfly catastrophes as paraxial beams and systematically investigate the connection to their appearance and characteristics predicted by catastrophe theory. Therefore, we evaluate theoretically corresponding diffraction integrals and analyze their potential functions in order to find caustics as critical degenerated positions that strongly influence the beams’ properties. Moreover, we demonstrate experimentally that their Fourier components are located on distributions that obey polynomial expressions.

We start with analyzing theoretically the class of caustic beams as well as their Fourier properties. Specifically, we contribute to analytical, numerical and experimental investigations with respect to the swallowtail and butterfly catastrophes by creating optical swallowtail and butterfly beams. Thereby, we describe their caustics with respect to control parameter space that defines the beam’s properties.

## 2 Embedding the swallowtail catastrophe in paraxial beams

We consider the caustic beam depending on one state variable emerging from the canonical diffraction catastrophe integral

(1) |

whose properties are completely determined by the properties of the canonical potential function in the oscillating integrand

(2) |

where . The potential function depends on a state variable and control parameters [42], which are components of the control parameter vector with corresponding scaling factors to make the exponent dimensionless. spans the control parameter space.

Caustics are defined as abrupt transitions of the optical ray system equivalent to bifurcations in nonlinear science by changing the number of crossing rays in each point in space [41]. Depending on the dimensionality of the control parameter space, the caustic arises as a point, line, surface or hypersurface. At these transitions, is stationary in for fixed . The caustics exist at points where the gradient mapping becomes singular, i. e. that the solutions of stationary are degenerated. Therefore, also the Hessian determinant of in has to vanish [42, 41]. That is, we calculate caustic structures corresponding to

(3) |

Note that the well known Airy and Pearcey beams emerge as special cases of since and . In the case of the Pearcey beam, we are able to identify the control parameters with the transverse spatial coordinates to directly obtain the transverse electric field distribution [37]. Instead, the famous 2D Airy beam is constructed by multiplying separable, orthogonal coordinates: [44, 45] and therefore does not present a generic structure predicted by catastrophe theory.

In the following, we develop a construction scheme to create higher-order () caustic beams by identifying two control parameters with spatial transverse coordinates , where and . The remaining control parameters are chosen as being constant and play an important role in controlling further properties of the realized beams.

Swallowtail beams result as solutions , where is one of the control parameters that was chosen to be constant. Thus, three different swallowtail beams are can be constructed, from which each exhibits unique features. However, all share characteristic properties of a swallowtail catastrophe. Furthermore, we create and analyze butterfly beams by calculating . Here, are two control parameters that stay constant. Consequently, six different butterfly beams complete the class of butterfly beams. All of them exhibit characteristic properties determined by the butterfly catastrophe.

We start our investigations on the swallowtail beams demonstrating its of fundamental transverse field distributions, and subsequently connect it with the control parameter space.

The three swallowtail beams of this case of artificially created and tailored caustic beams are shown in their most fundamental appearance in Fig. 1. In each case we set the constant control parameter to zero. The beams are constructed by numerically solving equation using adoptions of methods as described in [46]. These numerical results are shown in A, and were as well realized in experiments, shown in B C, where beneath the real space distributions of B, Fourier spectra C have additionally been recorded. Insets show corresponding calculated and measured spatial phase distributions. Light fields have been created using a phase-only spatial light modulator (SLM) that addresses light fields in real space. The transverse scale is chosen to be , and characteristic scale factors were . Experimental details are given in the SM.

White lines in A are calculated solutions of equation and therefore represent cross sections through the caustic surface in parameter space . We performed a parametrization of the caustic surface in dependence of one parametrization parameter and the remaining constant control parameter. The resulting expressions are given in the SM.

A 2D Fourier transform of Eq. (1) for the swallowtail () and butterfly beams () are performed analytically. For the explicit derivation of the Fourier transforms of these caustic beams we refer to the Supplemental Material (SM).

We pronounce the simplicity of the Fourier spectrum of caustic beams, since these are located on polynomial expressions of degree up to or expressions with rational exponents that are capable of forming cusps. For instance, the Pearcey beam’s Fourier components are located on a parabola with phase distribution [37]. In the SM we show that the beam (C3) exhibits a parabolic , the beam (C2) a cubic , and the beam (C1) even a cuspoid spectral distribution with corresponding phase functions, depending on which cross section in control parameter space (i.e. choice of ) is regarded. Note that by increasing the order of a caustic beam, new distributions of spectral components do not always arise. A good example is the swallowtail beam where const. whose Fourier components were located on a parabola similar to those of a Pearcey beam. Similarly, we show in the SM that butterfly beams have, beneath parabolic and cubic distributions, additionally quartic and cuspoid Fourier component distributions.

Catastrophes form geometrically stable structures where the bifurcation characteristics of the nonlinear potential function suddenly change. The transfer to optics leads to caustics that mark separated areas of different numbers of crossing light rays. Here we clearly see that the beam in A3 exhibits an area of two-beam interference with spatially varying -vectors on the left side of the caustic and performs an abrupt transition at the caustic interface to zero beams on the right side of low intensity.

Numerically calculated swallowtail beams are in very good agreement with the experimental realizations in both, spatial intensity and phase distributions. The predicted caustics indicated as white lines have been parametrized according to equation and fit very well the numerical realizations in scale and form. One of the most striking characteristics of the beams, strongly determining their properties, is the distribution of Fourier coefficients, which is perfectly predicted by the theoretical description stated in the SM.

Note that we intentionally omit any discussion of apodizing exponential or Gaussian functions in order to keep the total intensity of the light fields finite, as was suggested for both, the Airy [44] and Pearcey beam [37]. The fundamental physics of the infinite power beams discussed here resembles that of the apodized ones, and analytic expressions would become extensively large. Nevertheless, for the experimental realization the finite transverse dimensions of the light fields are determined by the size of the SLM.

## 3 Cross sections of the swallowtail caustic

Each caustic beam’s properties are connected to the distribution and kind of their correspondent caustics, may they consist of fold, cusp, or a combination of these catastrophes, which can be found by analyzing the potential functions. In order to investigate and image the caustic surface of the swallowtail catastrophe of control parameter space in different light structures, we arbitrarily chose to be constant and create beams with different values of . Thus, by varying a sliced scan of the surface in control parameter space can be performed whose properties manifest at the light fields.

The resulting swallowtail beams are depicted in Fig. 2 A B with varying control parameter , which performs a transition from positive (A1) to zero (A2) to negative (A3) values. This allows us to trace and illustrate fundamental properties of the swallowtail catastrophe that now manifests in the corresponding light structure: the caustic surface in control parameter space is defined by abrupt transitions of the system between areas of different numbers of crossing beams (in the regime of ray optics) [47]. By varying the control parameter , we have demonstrated the transition of the system from a swallowtail beam that consists of areas with 2 beam interference and areas with 0 beams () over the cuspoid point at to a swallowtail beam that clearly shows areas of 2 as well as 4 beams interference and low intensity areas (0 beams) for , as depicted schematically in C. These system transitions of the swallowtail catastrophe are additionally imaged by the caustic surface in control parameter space in Fig. 2D, which is in perfect agreement with previous discussions [26, 43, 41]. The caustic surface results from the parametrization for constant described in the SM.

## 4 The butterfly caustic

Each beam of the class of caustic beams that depends on one state variable , e.g. Airy, Pearcey or the three swallowtail beams, is striking due to its unique properties. Their static intensity and phase distributions in real and Fourier space are closely connected to catastrophe theory.

Similarities in the properties of these different beams manifest predominantly in their Fourier spectra. Discussions in the SM show that by increasing the order of the caustic beam, new properties will always emerge, since completely new distributions of Fourier components, like for instance higher-order polynomials, will arise.

Therefore, we introduce the six butterfly beams that complete the transfer from catastrophe theory to optics for all catastrophes that belong to the elementary catastrophes and are determined by one state variable . Those butterfly beams, whose Fourier spectra cannot be expressed in terms of lower-order caustic beams, are worth investigating thoroughly, since these promising light fields will show new individual propagation characteristics.

Fig. 3 shows the transverse light field distributions for each of the six butterfly beams in both, real and Fourier space. Numerical calculations of equation (1) are in very good agreement with the experimental realizations. Insets show corresponding calculated or measured spatial phases. Again, the experimentally obtained Fourier spectra are represented by polynomials and expressions with rational exponents, as predicted in the SM.

Similar to the treatment of the swallowtail catastrophes, we state a parametrization for the butterfly catastrophe in the SM. The parametrization was performed for one parameter and the respective two constant control parameters. For this higher-order catastrophe, control parameter space is four-dimensional. We restricted ourselves to image caustic lines (indicated in Fig. 3A as white lines) by choosing corresponding cross sections in parameter space. Areas with different numbers of crossing rays are clearly observable.

For constant the intensity of the spectra always diverges at the origin of Fourier space, which was also observed experimentally where the intensity increases extremely (but finitely). Again, the transverse picture section is chosen to be , and characteristic scale factors were . All constant control parameters are zero.

## 5 Conclusion

To conclude, we extended the diversity of caustic structures in paraxial light by numerically and experimentally embedding swallowtail and butterfly caustics in artificially tailored transverse light structures. In the framework of catastrophe theory we proved that the key properties of the related catastrophes are preserved in the paraxial regime. Thus, we demonstrated that emerging characteristics of higher-order swallowtail and butterfly catastrophes manifest at these optical paraxial beams in terms of abrupt transitions of the number of crossing beams. Dynamics of the potential function of the diffraction integral, which highly influences the beams’ properties, were investigated with respect to control parameter space. Thereby, we showed cross sections through three- and four-dimensional control parameter space and shaped the beams according to corresponding parametrized caustic (hyper-) surfaces.

Fourier spectra of these intriguing beams were calculated analytically and observed experimentally. We demonstrated, that spectral components of these caustics beams are distributed on polynomial expressions of increasing degree or are located on expressions with rational exponents. Some spectra are calculated to exhibit diverging energies at the origin, resulting in highly increasing but finite energies in experimental realizations.

By accessing higher-order catastrophic light structures, our findings allow investigating fundamental characteristics and properties of catastrophes in optics and in particular to utilize the caustics’ high intensities and propagation properties beyond the established caustic Airy and Pearcey beams. The trajectories and focusing effects of caustic light are promising to be discussed for designing tailored waveguides, particle manipulation, material processing, or improved super imaging in microscopic and sub-diffractive applications.

## Acknowledgement

We gratefully acknowledge fruitful discussions with Sir Michael Berry and Mark Dennis.

## References

- [1] Diaferio A 1999 Mon. Not. R. Astron. Soc. 309 610–622
- [2] Nesvadba N, Lehnert M, Eisenhauer F, Genzel R, Seitz S, Davies R, Saglia R, Lutz D, Tacconi L, Bender R and Abuter R 2006 The Astrophysical Journal 650 661–668
- [3] Gifford D, Miller C and Kern N 2013 The Astrophysical Journal 773 116
- [4] Oliva T A, Peters M H and Murthy H S K 1981 Behavioural Science 26 153–162
- [5] Guastello S J 1982 Behav. Sci. 27 259–272
- [6] Sheridan J E and Abelson M A 1983 AMJ 26 418–436
- [7] Hardy L 1996 The Sport Psychologist 10 140–156
- [8] Chao Y Y and Pierson W J 1972 Journal of Geophysical Research 77 4545 – 4554
- [9] LeBlond P H and Mysak L A 1978 Waves in the Ocean (Elsevier Oceanography Series) chap Waves near a caustic
- [10] Mei C C 1989 Applied Dynamics of Ocean Surface Waves (Advanced Series on Ocean Engineering, V. 1) (World Scientific Publishing Co Pte Ltd) chap The Neighborhood of a Straight Caustic
- [11] Kanaya K, Oho E, Adachi K, Yamamoto Y and Doi H 1990 Micron And Microscopica Acta 21 57–68
- [12] Patterson E A and Whelan M P 2008 Nanotechnology 19 105502
- [13] Kennedy S M, Zheng C X, Tang W X, Paganin D M and Jesson D E 2011 Ultramicroscopy 111 356–363
- [14] Petersen T C, Weyland M, Paganin D M, Simula T P, Eastwood S A and Morgan M J 2013 Phys. Rev. Lett. 110 1–5
- [15] Tavabi A H, Migunov V, Dwyer C, Dunin-Borkowski R E and Pozzi G 2015 Ultramicroscopy 157 57–64
- [16] Höhmann R, Kuhl U, Stöckmann H J, Kaplan L and Heller E J 2010 Phys. Rev. Lett. 104 093901
- [17] Mathis A, Froehly L, Toenger S, Dias F, Genty G and Dudley J M 2015 Scientific Reports 5 12822
- [18] Born M and Wolf E 1970 Principles of Optics 4th ed (Pergamon Press)
- [19] Berry M V 1975 J. Phys. A: Math. Gen. 8 566–584
- [20] Berry M V and Balazs N L 1979 Am. J. Phys. 47 264–267
- [21] Blandford R and Narayan R 1986 The Astrophysical Journal 568–582
- [22] Wang Y and Turner E 1997 Mon. Not. R. Astron. Soc. 292 863–870
- [23] Kiser T and Pauly M 2012 Caustic Art (EPFL Technical Report)
- [24] Ball P 2013 New Scientist 217 40–43
- [25] Schwartzburg Y, Testuz R, Tagliasacchi A and Pauly M 2014 ACM Trans. Graph. 33 74
- [26] Arnol’d V I and Wassermann G S 2003 Catastrophe Theory 3rd ed (Springer Verlag)
- [27] Baumgartl J, Mazilu M and Dholakia K 2008 Nature Photon. 2 675–678
- [28] Polynkin P, Kolesik M, Moloney J V, Siviloglou G A and Christodoulides D N 2009 Science 324 229–232
- [29] Diebel F, Bokic B M, Timotijevic D V, Jovic Savic D M and Denz C 2015 Opt. Express 23 24351
- [30] Woerdemann M, Alpmann C and Denz C 2011 Appl. Phys. Lett. 98 111101
- [31] Alpmann C, Schöler C and Denz C 2015 Appl. Phys. Lett. 106 241102
- [32] Rose P, Boguslawski M and Denz C 2012 New J. Phys. 14 033018
- [33] Durnin J, Miceli J J and Eberly J H 1987 Phys. Rev. Lett. 58 1499–1501
- [34] Diebel F, Rose P, Boguslawski M and Denz C 2014 Appl. Phys. Lett. 104 191101
- [35] Diebel F, Rose P, Boguslawski M and Denz C 2016 New J. Phys. 18 053038
- [36] Zhang Y, Liu J, Wen F, Li C, Zhang Z, Zhang Y and Belic M R 2016 arXiv 1605.09525v1
- [37] Ring J D, Lindberg J, Mourka A, Mazilu M, Dholakia K and Dennis M R 2012 Optics Express 20 18955–18966
- [38] Zeeman E C 1979 Structural Stability in Physics (Berlin, Heidelberg: Springer Berlin Heidelberg) chap Catastrophe Theory, pp 12–22
- [39] Arnol’d V I 1975 Russ. Math. Survs. 30 1–75
- [40] Thom R 1975 Trans. by D. Fowler. Reading, Mass.: Benjamin
- [41] Berry M V and Upstill C 1980 Prog. in Optics 28 257–346
- [42] Berry M V, Nye J F and Wright F J 1979 Phil. Trans. R. Soc. Lond. A 291 453–484
- [43] Nye J 2006 Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences 462 2299–2313
- [44] Siviloglou G A and Christodoulides D N 2007 Opt. Lett. 32 979–981
- [45] Siviloglou G A, Broky J, Dogariu A and Christodoulides D N 2007 Phys. Rev. Lett. 99 213901
- [46] Connor J N L and Curtis P R 1982 J. Phys. A: Math. Gen. 15 1179–1190
- [47] Kravtsov Y A and Orlov Y I 1983 Usp. Fiz. Nauk 141 591–627

## Supplementary Material | Optical catastrophes of the swallowtail and butterfly beams

### Fourier spectra of swallowtail and butterfly beams

We perform a 2D Fourier transform of equation (1) in the main article regarding two control parameters . Corresponding Fourier frequencies of with components are with components , i.e. here we introduce . With we denote all remaining control parameters except . We neglected scaling factors for reasons of clarity.

We identify two control parameters with spatial transverse coordinates and treat remaining control parameters as constants. Thus, in case of the swallowtail beam Sw , which is defined by three control parameters , three different transverse field distributions can be created by keeping constant one control parameter. All resulting fields show characteristics of the swallowtail catastrophe.

(1) | |||||

(2) | |||||

(3) |

Thereby, . The Kronecker delta equals if and equals for . Thus, only allows positive (including ) to contribute in the corresponding term. The spectra for the beams, as calculated above, exhibit Fourier components that are each only located in two quadrants: For const., the spectrum is restricted to . For const., Fourier components are only allowed for and , whereas for const., the spectrum is restricted to .

Expanding our approach to butterfly beams Bu and still identifying two control parameters with spatial coordinates, in this case six transverse field distributions can be created.

(4) | |||||

(5) | |||||

(6) | |||||

(7) | |||||

(8) | |||||

(9) |

In general, the number of transverse field distributions for each catastrophe can be calculated by considering standard combinatorial methods: , where is the co-dimension of the corresponding potential function .

### Parametrization of the caustics at cross-sections of the swallowtail and butterfly catastrophe

We performed a parametrization of the caustic surface in dependence of one parameter and the remaining constant control parameters. The resulting expressions are given in table 1 and 2, where the constant control parameters are indicated at the first column. These equations therefore describe cross sections through the three- and four-dimensional control parameter spaces.

### Experimental Details

Fig. 4 shows the experimental setup: An expanded, collimated frequency-doubled Nd:YVO laser beam is split up in two parts, which are linearly polarized. One serves as structure beam, whereas the other represents the reference beam. The structure beam is modulated by a HOLOEYE HEO 1080P reflective LCOS phase-only SLM with full HD resolution, which we have applied to simultaneously modulate both amplitude and phase. Therefore, an appropriate Fourier filter (FF) was used. As the imaging system, we mounted a camera and microscope objective with fixed distances on a -shift moving stage that allows us to scan the light field in the longitudinal direction. Measuring the spatial phase distribution is possible by superimposing a tilted reference beam, which can be enabled on demand. By temporarily installing lens L, the Fourier spectrum becomes accessible.