ORIGINAL RESEARCH article
Superoscillations in Finite-Energy Airy Beams
- Department of Physics, Guangdong University of Petrochemical Technology, Maoming, China
Superoscillations naturally arise in optical fields with dense packing of nodal points of amplitude. Airy wave packets are highly oscillatory and rich of phase singularities. We study to the best of our knowledge, for the first time, the superoscillatory behavior in a band-limited Airy beam whose spectrum is sharply truncated. Our results show that not as expected, the superoscillations occur outside of the Airy-like region, but in regions above a defining line where the beam stops being Airy-like. The degree of superoscillation can be very high there.
The phenomenon that a band-limited signal can oscillate faster, although locally somewhere in its domain, than its largest frequency is called superoscillation. The discovery of this phenomenon dates back to the studies in quantum mechanics by Aharonov et al. . This topic stimulates recent interest in optics [2-4] because of its close connection to, in particular, super-resolution imaging. The latter is usually achieved in the scheme of recovering evanescent waves. More interestingly, the desired large local frequency, or local wave number in optical case, naturally occurs near phase singularities, where it diverges. This connection leads to a recent method of achieving superoscillation using optical vortexed waves .
In fact, the superoscillatory behavior is ubiquitous in optics. Dennis et al.  showed that for certain speckle patterns, the area of the superoscillating region can be of 1/3 in fraction. Makris et al.  constructed diffraction-less beams that can transport subwavelength features into the far field. These studies stimulate us to study the superoscillations embedded in optical Airy beams. Our motivation comes from the facts that Airy beams are oscillatory and rich in optical singularities, hinting that super large local wave numbers may occur even though the beam is inherently paraxial. First discovered by Berry and Balazs  as a set of solutions to the Schrodinger equation, Airy wave packets raise great interest [8, 9] in the context of beam optics due to their peculiar features as diffraction-less propagation, self-acceleration, and self-healing. However, a complete Airy beam is of unlimited spectrum. To facilitate our current task, after briefly reviewing the theory of and the criterion for superoscillation in Section 2, we introduce in Section 3 the method of sharply truncating the Airy beam spectrum; in Section 4, we give the method to compute the local wave number through its relation to optical currents, and then study the superoscillatory area in both a longisection and a transverse section; finally, we give a conclusion of the study.
2 Superoscillations in Light Fields
Suppose we have a general optical wave field, denoted by ϕ(r), propagating in free space, then the vector r denotes the position in three-dimensional space. We write ϕ(r) in terms of its amplitude ρ(r) and phase φ(r):
If ϕ is a plane wave, that is, the eigenfunction of momentum, associated with ϕ, then there is a unique and global wavenumber k, with corresponding spatial frequency v = k/2π. However, for general waves that have complicated distribution of φ(r), no such global wave number can be designated. What we have instead is the local wave vector k, defined as the gradient of the phase function:
For most complex fields, even φ itself cannot be explicitly given. For such cases, one can use a relationship between the phase gradient and the optical current  J,
where the symbol
3 Optical Incomplete Airy Beams
In this article, we will formulate the propagation problem in the 2 D scenario, that is, one transverse plus one longitudinal dimension. We will assign the x coordinate to be the transverse position and ξ coordinate to be the longitudinal one. A full 3 D formulation can be obtained by incorporating another transverse dimension (the y coordinate), which is formally identical to x. Let ϕ(x,ξ) be the amplitude of a light wave, satisfying the normalized paraxial wave
One way to construct the Airy beam is to conceive that the wave is propagated from an initial, say ξ = 0 plane, with its distribution ϕ(x,0) being the Airy function Ai(x):
This equation indicates that the spatial spectrum of the Airy function is given by the unimodular cubic phase factor exp(is3/3). To obtain the beam profile at any longitudinal position Ai(x,ξ), we model the paraxial propagation through the method of angular spectrum. The wave number along ξ, denoted as kξ, can be computed under the Fresnel approximation:
Thus, we have
The second line of Eq. (8) shows that the beam is propagation invariant up to translation, as long as we set the parameters a and b in the integration to be
When exp(iξ) is decreased, we finally have the complex envelope of the Airy beam.
From Eq. (10), we can readily see that the point of peak intensity of the wave packet moves from x = 0 in the initial plane to the point x = ξ2/4 for a particular ξ. This phenomenon is the so-called self-acceleration of Airy beams. The path of the peak intensity in the x−ξ plane is evidently a parabola.
So far what we obtain from Eq. (10) is an Airy beam with unlimited spectrum. However, for the purpose of tracking the superoscillatory behavior, we are in need of a band-limited light field with a sharp spectrum limit. There are different ways of truncating the beam, leading to different modulations of the spectrum. The exponential truncation introduced in  will introduce a long-tailed Gaussian filter to the spectrum. A sharp limiting spectrum can be achieved by imposing a rectangular window, as was done in . We choose the maximum spectrum to be S and use (11), and now the Airy beam changes to the incomplete one:
For the sake of clarity, we omit the integration limits of Eq. (11) and rewrite it into a compact form:
We use the subscript S to denote the maximum transverse wave number. Note that due to the truncation, the second term in Eq. (12) is generally no longer a real function of x and ξ, as is the complete Airy function.
4 Superoscillations in Incomplete Airy Beams
4.1 Superoscillations in a Longisection
Now our task is to compute the phase gradient of the incomplete Airy beam given by (12). As the second term of (12) is complex, the overall phase φ(x,ξ) of the beam is not explicit, and we invoke formula 4:
We have numerically solved Eqs. (11) and (15) for an incomplete Airy beam truncated at S = 10, with x being in the range of 2π[−12.8,12.8] and the propagation range
FIGURE 1. Density plot of the normalized beam intensity of the incomplete Airy beam, Ais=10(x,ξ). The white dashed line is the defining line x = Sξ−S2, separating Airy-like and non–Airy-like regions. The beam retains Airy features to the right of the line. The white solid line is the main lobe parabola, x = ξ2/4.
According to , regions satisfying x > Sξ−S2, that is, fields above the dashed white line, lost Airy-like features.
Both the total and the x− component phase gradients are computed and shown in Figure 2. Please note that for the total gradient
FIGURE 2. Density plots of the local phase gradients: a,
4.2 Superoscillations in Transverse Cross Sections
Now, we turn to the superoscillations in cross sections of a three-dimensional Airy beam. Without loss of generality, we concentrate on the x–y symmetric beam, that is, the beam is just the simple product of two incomplete one-dimensional Airy beams, one along each of the two transverse coordinates. The two-factor beams have the same truncations and other parameters. We denote the three-dimensional beam as Ai3, and it is given in the simple separable form:
We denote the overall phase function of Ai3 by Ѱ(x,y) in a transverse cross section specified by ξ = ξ0. Simple calculations give the following relation:
and we can use the computational results in the longisection to construct the local phase gradients in the transverse plane.
We used the same parameter set as in the 2 D case, for both the x and y directions. As shown in Figure 3, the main lobe position in each sub-figures is the intersection point of the two white solid line, which is computed using the parabola equation for a specified ξ0 position. The dashed black lines are the defining lines separating different regions. Regions to the right of and above the horizontal black lines are Airy-like. As shown in Figure 1, the dashed and solid lines get tangent to each other at a particular value of ξ.
FIGURE 3. Shaded contour plots of the transverse phase gradient
It is obvious that as the beam propagates, the Airy-like region (squares circumscribed by the white and black lines) gets increasingly smaller, while the superoscillating region gets increasingly larger. For small distance, ξ = 0.176, part of the lower left corner enclosed by the parabola is Airy-like, and part is not, while superoscillations occur only on discrete lines. For moderate distance, the beam region becomes totally non–Airy-like, but continuous area emerges as superoscillatory. For all the three cases, the Airy-like regions are always non-superoscillatory, agreeing with cases in the longisection.
In conclusion, we studied the superoscillatory behavior in an incomplete Airy beam in both a longisection and transverse cross sections. Sharply truncating the beam spectrum leads to the emergence of non–Airy-like regions in the beam area which gradually enlarge during propagation. However, in such regions, the fields can oscillate much faster than those in both the largest frequency and the Airy-like region. No finite area in the Airy-like region is found to support superoscillations.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author/s.
YH conceived the research and derived equations about the incomplete Airy beams and the phase gradients, and XC and LN wrote the computational codes, computed the numerical results, and made the figure plots; the three authors wrote the manuscript.
The research presented in this article is financially supported by the Natural Science Foundation of Guangdong Province, China (Grant No. 2016A030307026).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors thank Tianle Li and Weiling Zhu in our department for their helpful discussions.
Keywords: airy beam, superoscillation, propagation, paraxial wave equation, diffraction-free
Citation: He Y, Chen X and Niu L (2021) Superoscillations in Finite-Energy Airy Beams. Front. Phys. 9:650160. doi: 10.3389/fphy.2021.650160
Received: 06 January 2021; Accepted: 25 January 2021;
Published: 21 April 2021.
Edited by:Hao Chen, Shanghai Institute of Microsystem and Information Technology (CAS), China
Reviewed by:Guoquan Zhou, Zhejiang Agriculture and Forestry University, China
Rui-Pin Chen, Zhejiang Sci-Tech University, China
Lingfei Xu, Shanghai Electronical Mechanisim Engineering Institude, China
Copyright © 2021 He, Chen and Niu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yingjun He, email@example.com