Abstract
We present an algorithm for manipulating and controlling 3-D field patterns, with energy confined to the narrow vicinity of predefined 3-D trajectories in free-space, which are of arbitrary curvature and torsion. This is done by setting the aperture field’s phase to form smooth caustic surfaces that include the desired trajectory. The aperture amplitude distribution is constructed to manipulate both the on-axis intensity profile and the off-axis beam-width, and is updated iteratively. Once the aperture distribution is calculated, the radiation from a finite sampled aperture is computed numerically using a Fast Fourier Transform-based scheme. This allows for both verification of the design and examination of its sensitivity to parameters of realistic discrete implementation. The algorithm is demonstrated for the cases of an Airy beam of a planar trajectory, as well as for helical and conical-helical trajectory beams.
© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
Enhanced intensity field distributions that follow curved trajectories in free-space have been a subject of intense research in recent years. These field intensity profiles, often termed accelerating beams [1–10], have gained attention due to their weakly-diffractive curved and self-healing nature [11–13]. With their energy confined to the narrow vicinity of predefined trajectories, such field patterns have been shown useful for various applications. These include particle manipulation [14–18], super-resolution imaging [19,20], plasmon excitation [21,22], light-sheet microscopy [23], optical coherence tomography [24], and plasma-channel generation [25], to name but a few. The ray construction is nondispersive, and hence the theory can be formulated a priori also in the short-pulse time domain [26,27].
Several methods have been proposed, aiming at tailoring the on-axis intensity profile of curved beams. In [17], analytical solutions of Maxwell’s equations with losses, propagating in absorbing media, were presented. These wave objects maintain their peak intensity along the beam-axis for large distances. Arbitrarily tailoring of the intensity distribution along the beam-axis of two-dimensional (2-D) surface-plasmon beams was implemented in [28]. The generation of such field patterns requires the design of source field distributions, often defined over wide apertures. For example, strong intensity profiles can be obtained at focal points of a bundle of geometrical optics (GO) rays (Bessel-type beams) [29,30] or along a smooth caustic of GO rays [6,7,31–33]. For these cases, aperture field distributions (phase and amplitude) that give rise to the desired pattern can be calculated analytically.
Several methods have been proposed for shaping the accelerating beam trajectory and intensity profile. In [34], non-paraxial accelerating beams over arbitrary convex trajectories were designed. Some control over the intensity profiles along the caustic was obtained by changing the aperture field amplitude. In [35], independent control of both the trajectory and the maximum amplitude along it was achieved for 2-D caustic beams (CBs) using catastrophe theory. A practical algorithm for designing 2-D caustic beams, that meander along trajectories, consisting of both convex and concave segments, was introduced in [32]. The implementation was based on the generation of a caustic at the desired trajectory. This was done by “back-tracing” rays from the predefined beam trajectory to the source aperture, to determine the source distribution’s phase. The aperture amplitude distribution was set to form a uniform intensity of the produced CB along its trajectory. These field design approaches are limited to particular 2-D cases with limited control over the beams’ parameters. Yet, the principles used in these works inspire the design of more general field patterns that are of intensity confined to the vicinity of arbitrary three-dimensional (3-D) trajectories. Specifically, we are interested in the design of aperture distributions that produce patterns with 3-D trajectories of various curvature and torsion parameters. Moreover, along these trajectories, we would like to control the intensity profile and the beam-width.
In this work, we present a 3-D extension of the algorithm in [32], for the design of beam field patterns, along desired beam-axis trajectories with arbitrary curvature and torsion (as defined in section 2). In the proposed scheme, confining the field to the vicinity of the predefined beam-axis trajectory is done by designing the aperture field’s phase, to give rise to a smooth caustic surface (section 3.1) that formed around the beam-axis "skeleton". Following the 2-D methodology (and using the definitions in section 3.2), we design the aperture amplitude distribution, with the aim of manipulating the on-axis intensity profile (section 3.3). In the 3-D case, another controllable feature of the pattern is the off-axis beam-width. It can be manipulated via the corresponding aperture distribution’s effective off-axis width (section 3.4). The proposed technique is validated via numerical experimentation; Once the aperture distribution is calculated, the radiation from a finite sampled aperture is computed numerically using a Fast Fourier Transform (FFT)-based scheme (section 4). This technique also enables the examination of the design’s sensitivity to parameters of realistic discrete implementation, in an experimental set-up, for example by using a finite emitters array. The algorithm is demonstrated for the cases of Airy beam planar trajectory, as well as for helical and conical-helical trajectories, in sections 5.1, 5.2, and 5.3, respectively.
2. Problem definition
Given a trajectory in a 3-D domain
with $\sigma \in [\sigma _{\textrm {min}},\sigma _{\textrm {max}}]$ denoting the arc length along the beam-axis trajectory, we aim at designing a source (aperture field) in the $z=0$ plane that produces a 3-D beam-field, with a beam-axis that follows this curved trajectory, in free-space. The trajectory can be intricate, with arbitrary curvature and torsion. In this work, we assume that the characteristics of the trajectory are large on a wavelength scale and limit the discussion to trajectories with tangent vectors of positive $z$-component. The beam-field’s desired transverse width is denoted $W_{\textrm {b}}(\sigma )$ (see Fig. 1) and $I_{\textrm {b}}(\sigma )$ denotes the desired intensity profile along the beam-axis. Our design approach relies on the observation that the radiated field’s intensity can be confined to the vicinity of a given trajectory ${\mathbf {r}}_{\textrm {b}}(\sigma )$, if the aperture phase, asymptotically, gives rise to a caustic surface that contains ${\mathbf {r}}_{\textrm {b}}(\sigma )$. In other words, the approach produces a 3-D CB along the trajectory. To that end, the designed aperture field, denoted $u_0(x',y')$, is sought in the form of a ray-field where $k=\omega /c$ is the wavenumber, with $c$ being the free-space wave velocity. Also in (2), $(x',y')$ are coordinates in the $z=0$ plane, and $A$ and $S$ are identified as the aperture field’s amplitude and phase (Eikonal) distributions, respectively. In what follows, we present a method for designing the two functions, $A$ and $S$, given the desired beam-axis trajectory, ${\mathbf {r}}_{\textrm {b}}(\sigma )$, intensity, $I_{\textrm {b}}(\sigma )$, and beam-width, $W_{\textrm {b}}(\sigma )$.3. Aperture field design
The design of an aperture that gives rise to a CB of desired properties consists of three parts: first, an aperture phase distribution is designed, in accordance with the desired CB’s beam-axis trajectory. It is worth noting that, for planar trajectories, this part reduces to application of the 2-D technique in [32], for computing the phase along the trajectory’s projection onto the $z=0$ plane, followed by an extension of the phase in the normal direction. However, for treatment of arbitrary 3-D trajectories, both steps of forming the phase distribution must be generalized. Then, the aperture amplitude distribution is designed to control (i) the CB’s on-axis intensity profile and (ii) the CB’s width. While the phase distribution is computed only once per given trajectory, the aperture amplitude distribution can be further adjusted, iteratively.
3.1 Beam-axis manipulations
The proposed method attempts to form a caustic on a surface containing the beam-axis ${\mathbf {r}}_{\textrm {b}}(\sigma )$. This is done by setting the aperture phase to give rise, at each point on the aperture, to a ray that participates in forming the caustic surface. First, we construct the phase that forms the beam-axis, which can be though of as the caustic "skeleton". Refereing to Fig. 1, from each point ${\mathbf {r}}_{\textrm {b}}(\sigma )$ along the beam-axis, we follow the tangent ray back to its emanating point on the aperture, denoted ${\mathbf {r}}'_{\textrm {a}}(\sigma )$. The ray angles formed between the ray and the corresponding axes, $\theta _{x,y}$, are obtained via the unit vector tangent to the parametric trajectory in (1). i.e.,
Following this, we define the local curve coordinates, $(\sigma ',n')$, on the aperture plane, where $\sigma '$ is the arc length along the curve ${\mathbf {r}}'_{\textrm {a}}(\sigma ')$ and $n'$ denotes the distance of a point $(x',y')$ from ${\mathbf {r}}'_{\textrm {a}}(\sigma ')$ (in the normal direction), with corresponding unit vectors, ${\hat {\boldsymbol \sigma }}'$ and ${\hat {\mathbf {n}}}'$ (see Fig. 2). Here and henceforth, we use hat over bold fonts to denote unit vectors. Thus, the parameter $\sigma \in [\sigma _{\textrm {min}},\sigma _{\textrm {max}}]$ is mapped to $\sigma '\in [0,\sigma '_{\textrm {max}}]$ via (4). This mapping is denoted by $\sigma '(\sigma )$. From each point $(x',y')$ on the curve ${\mathbf {r}}'_{\textrm {a}}(\sigma ')$, the ray that forms the caustic over the beam-axis trajectory emanates in the direction of the unit vector
Using these definitions, we first evaluate the phase, $S(\sigma ',n')$, along the curve ${\mathbf {r}}'_{\textrm {a}}(\sigma ')$, i.e., $S(\sigma ',0)$. By using the angles $[\theta _x(\sigma '),\theta _y(\sigma ')]$, we evaluate the phase gradient along the curve ${\mathbf {r}}'_{\textrm {a}}(\sigma ')$ via [36]
Next, the phase is extended to other points on the aperture. The extension is not unique and it dictates the caustic surface’s shape. First we consider the phase extension for the special case of planar beam-axis where $y(\sigma )=0$ in (1). In this case, $\cos \theta _y(\sigma ')=0$ and therefore, using (6), the normal phase derivative $\partial S /\partial n' =\partial S /\partial y= 0$ for all points along ${\mathbf {r}}'_{\textrm {a}}(\sigma ')$. Any extension of the phase $S(\sigma ',n')$ is subject to this boundary condition. The simplest form of extension is the linear one in which $\partial S /\partial n'=0$ for all points over the aperture. This choice will reproduce the caustic curve for rays in any constant $y$ plane, forming a $\partial /\partial _y=0$ caustic surface (within the beam width domain), see Fig. 3 (though our method of phase construction from tangent rays guarantees the formation of a caustic surface, this can be easily verified by evaluating the ray Jacobian using (3) and setting its determinant to zero). Other continuous extension, subject to the boundary condition, are possible but yield intricate caustic surfaces (a special interesting limit case is the conventional 3-D Airy beam where the caustic takes the form of a cusped double-layered caustic - see [9]).
The phase extension for the general (non-planar) beam-axis is performed in a manner similar to that of the planar beam extension: since the integration in (7) reconstructs the phase along the curve ${\mathbf {r}}'_{\textrm {a}}(\sigma ')$, it sets the directional derivative along the curve to $dS/d\sigma '={\mathbf \nabla }' S \cdot {\hat {\boldsymbol \sigma }}'$. In order to satisfy (6) for points over ${\mathbf {r}}'_{\textrm {a}}(\sigma ')$, the directional phase derivative normal to the curve should satisfy the boundary condition
3.2 Local beam coordinates
In order to manipulate and control the on-axis intensity and the beam-width, we apply the local beam coordinates that are associated with the beam-axis in (1). Denoting $\sigma$ the arc length along the beam-axis trajectory, the beam local coordinates are defined by the unit vectors ${\hat{\textbf t}}(\sigma ), {\hat{\textbf n}}(\sigma )$, and ${\hat{\textbf n}}_{\textrm {b}}={\hat{\textbf t}}\times {\hat{\textbf n}}$, being the tangent, normal, and bi-normal to ${\mathbf {r}}_{\textrm {b}}$ directions, respectively at a point ${\mathbf {r}}_{\textrm {b}}(\sigma )$ on ${\mathbf {r}}_{\textrm {b}}$. They are related by the Frenet-Serrer equations
By using these definitions, we can discuss the controlling of the beam intensity profile along each of the three local coordinates, $\sigma$, $n$ and $n_{\textrm {b}}$. Since ${\mathbf {r}}_{\textrm {b}}$ is a trajectory over a caustic surface, we identify the mechanism of enhanced intensity along the local coordinate $n$ as the caustic local structure. This implies that, asymptotically (and paraxially), the intensity profile in the direction normal to the caustic surface is, as in any simple caustic, of the form of the (squared) Airy function profile (not to be confused with the Airy beam trajectory). Therefore, in the proposed method, control over the width of the beam intensity in the local ${\hat {\mathbf {n}}}$ direction is not possible. However, the beam-width in the transverse $n_{\textrm {b}}$ coordinate and the intensity profile along $\sigma$ are controllable. In what follows, we propose a simple technique for separately controlling each of the two. To that end, in our implementation, we design the aperture field amplitude in (2) to be of the general form
We term $w_{\textrm {t}}$ and $w_{\textrm {l}}$ the transverse filter and longitudinal window, respectively. In the next sections we discuss the properties of these filter and window and their usage for controlling features of the beam intensity.3.3 On-axis beam intensity profile manipulations
The longitudinal (amplitude-equalizing) window $w_{\textrm {l}}(\sigma ')$ controls the field intensity profile along the beam-axis. The desired intensity profile, $I_{\textrm {b}}(\sigma )$, can be any continuous function of $\sigma$, pending on the application. Finer adjustments of the on-axis intensity profile can be carried out iteratively, after setting the transverse filter according to section 3.4. In this procedure, we start by setting
In each ($n^{\textrm {th}}$) step, the aperture field is propagated into the $z>0$ half space (by using the method described, later, in section 4). The field’s intensity, denoted $I_{\textrm {l}}^{(n)}(\sigma )$, is sampled along the beam-axis. We then use the mapping $\sigma '(\sigma )$ and adjust the amplitude-equalizing window via multiplication by the (real) function3.4 Off-axis beam-width manipulations
The tapering of the beam intensity along the coordinate $n_{\textrm {b}}$ is controlled by the transverse filter, $w_{\textrm {t}}(n',\sigma ')$ in (12) - a real function of finite support centered about $n'=0$. There is a trade-off between the desire for a sharp decay of the field intensity away from the beam-axis and the diffraction caused by abrupt changes in the aperture field. Therefore, the window profile, which determines the intensity profile parameters (width, decay rate, diffraction artifacts, etc.), should be selected carefully. Good candidates for the transverse filter, $w_{\textrm {t}}$, include the Gaussian windows, as well as other gracefully decaying window functions.
In the examples presented in section 5, we use a scaled Tukey window, in order to minimize the diffraction. This window is given by
4. Validation
In order to verify the proposed CB synthesis approach, the fields radiated by the aperture distribution are computed in the half-space $z>0$. Ideally, this should be done by applying the Kirchhoff-Huygens integral to the aperture distribution at $z=0$, such that
The computation of (22) for all points on the observer grid is accelerated by using a 2-D FFT -based algorithm, similar to that in [39]. The 2-D FFTs of $u_{0}[m',n']$ and $\partial _{z'}G[m,n,z]$ are first computed, and are then multiplied in the spectral domain, and the result is transformed back to spatial domain by using a 2-D inverse FFT to obtain $u[m,n,z]$. Note that, for an aperture that is modeled as an array of uniformly spaced point sources and measurements on a regular grid, this fast scheme for computing the field is exact (that is, in contrast to the $z$-propagation of sampled plane-wave spectra of the aperture fields).
In this work, we model the aperture as an array of such point sources, weighted by the sampled values of the aperture distribution. We use a typical discretization length of $\Delta x = \Delta y = \lambda /10$. For more accurate modeling of a particular realistic implementation, one could replace $\partial _{n'}G({\mathbf {r}}-{\mathbf {r}}')$ by a function describing the field produced by an element of the array. This field representation enables the study of the proposed synthesis approach’s sensitivity to implementation specifications dictated by realistic system limitations, e.g., array length, element spacing, and elemental sources’ field [4,11]. Finally, it should be noted that the need for fast radiation integral computation is even more pronounced when repeated field evaluation is required. That is the case when an iterative procedure for adjusting the beam’s (on-axis) intensity and (off-axis) width in sections are used, as in sections 3.3 and 3.4.
5. Examples
In this section, the proposed method is demonstrated for several representative examples. The approach is first validated for the simple case of a planar (beam-axis) trajectory. This case is geometrically very similar to the 2-D case, in terms of the design of aperture phase distribution. Therefore, it enables us to conveniently study the sensitivity of the design to the aperture’s discretization. Then, we present results for more complex cases of helical and helical-conical CBs. These cases challenge the proposed scheme and allow us to showcase its advantages.
5.1 Planar beam: analytical example
For the first example, we chose an Airy-type planar (torsion free) CB, for which the beam-axis is defined by the parametric trajectory
This beam-axis describes a planar trajectory on the $x=y$ plane, where $z_{\textrm {b}}(t)=\alpha \sqrt {t}$ is identified as the well-known Airy beam trajectory. In this example, we synthesize a constant-width beam with $W_{\textrm {b}}=15\lambda$ for $\alpha =5$. By using (3)–(6) we evaluate for ${\mathbf {r}}'_{\textrm {a}}(\sigma ')=\frac {-1}{\sqrt {2}}[\sigma ',\sigma ']$:In this example, the beam-axis in (23) is defined for $\sigma \in [\sigma _{\textrm {min}},\sigma _{\textrm {max}}]$. In order to avoid diffraction effects at the end of the beam-axis, ${\mathbf {r}}_{\textrm {b}}(\sigma _{\textrm {max}})$, and to obtain the desired beam-width there, we extend the beam-axis to the interval $\sigma \in [\sigma _{\textrm {min}},1.25\sigma _{max}]$ and calculate the phase and amplitude of the sources accordingly. This extension is used in the following two examples as well.
Following the procedure in section 3.4, the initial CB aperture distribution is obtained by setting $W_{\textrm {t}}=W_{\textrm {b}}/e_{\textrm {a}}=30\lambda$ with $r=1$ in (15). By using a grid spacing of $\lambda /10$, we obtain the aperture field distribution in Fig. 5. The corresponding CB-field intensity is evaluated according to section 4. Figure 6(a) depicts the field intensity over the $x=y$ plane, demonstrating that the on-axis intensity profile follows the desired trajectory (marked by the black dashed curve). Figure 6(b) plots the projection onto the $z=0$ plane, of points in the half-space $z>0$ of field intensity greater than half the maximum at their respective $z$-plane. The solid black line indicates the desired beam trajectory, whereas the dashed black lines to each side of the trajectory indicate the desired -3dB width. This point of view shows clearly the control of the beam-width in the ${\hat{\textbf n}}_{\textrm {b}}$ local direction.
This example is also a convenient test case for studying the sensitivity of the synthesized aperture distribution to constraints imposed by a realistic implementation. Specifically, the spacing between array elements that aim to mimic the continuous source distribution in the aperture may be larger than of the order of a wavelength. Such is the case, for example, of many existing spatial light modulators (SLMs) [1][4][11][40]. These devices for controlling the aperture phase are designed to operate at 300nm–1000nm wavelength, but often use a typical pixel edge length (also termed pixel pitch) that is, at best, currently, as small as $0.72$ $\mu \textrm {m}$ [41]. Antenna arrays are typically designed with a $\lambda /2$ spacing between elements. To evaluate the effect of aperture (under-) sampling, we examine the fields produced by aperture distributions sampled at various rates. For these fields, we calculate the standard deviations of the -3dB beam-width in the local normal axis. These are plotted in Fig. 7 as a function of the grid spacing (Note that, to maintain same grid spacing for both the aperture and the measurements, we use a fixed grid spacing of $\lambda /20$ (smaller than the numerical rule-of-thumb spacing of $\lambda /10$) for both aperture and measurement grids all cases, and for the coarser sampling cases, sparse filling of the aperture is used as needed). It can be seen that the ability to maintain the desired beam-width drops sharply for grid spacings greater than $\lambda /2$. This agrees with the well known $\lambda /2$-spacing requirement for antenna pattern design. Fig. 8 shows, for various cases of grid spacings, the points of maximal intensity for each $z$ plane, alongside the desired beam trajectory. It can be seen that even when exceeding the $\lambda /2$-spacing limit, e. g., for $\Delta =\lambda$, the field pattern’s maxima continue to follow the desired trajectory. It is only for spacings greater than $2\lambda$ that the pattern begins failing to follow the desired trajectory. In other words, in the $\sigma '$ direction, some under-sampling is tolerable, as the trajectory is dictated by the geometrical nature of the ray-field. That is not the case for sampling in the $n'$ direction. Aiming for the design of arbitrary trajectories, which are not necessary within a plane, the $\lambda /2$ restriction will hold for the entire aperture.
5.2 Helical beam
With the basic features of the proposed technique demonstrated, we now turn to examine its adequacy for the construction of CBs with more complex beam-axis trajectory (i.e., a trajectory with arbitrary curvature and torsion) and intensity profile. As our second example, we choose a CB that follows a right-handed helical beam-axis that is described by
where $R$ denotes the helix radius and the parameter $P$ sets the pitch to be $2\pi P$. The arc length of this trajectory is $\sigma = (t-t_{\textrm {min}}) \sqrt {R^2+P^2}$. In this example, the desired beam-width profile was $W_{\textrm {b}}(t)=A[1+sin\left (2\pi \frac {t-t_{\textrm {min}}}{t_{\textrm {max}}-t_{\textrm {min}}}\right )/5]$ and the on-axis intensity profile was $I_{\textrm {b}}(t)=1$.By setting $R=P=20\lambda$ and following the procedure in section 3, we numerically calculate $S(x',y')$. Next, we calculate the amplitude $A(x',y')$ by setting $r=0.8$ in (12), with the transverse filter $w_{\textrm {t}}(n',\sigma ')$ in (15) and the longitudinal window $w_{\textrm {l}}(\sigma ')$ in (13). Finally, we calculate the propagating fields in the $z>0$ half-space via (22). At every $\sigma$, the beam-width is defined as the Euclidean distance between the -3dB of the peak intensity points in the bi-normal local direction (see Fig. 9). Using the initial aperture design, we evaluate the intensity. In Fig. 10, the resulting iso-surfaces of -3dB of the peak intensity is plotted. The beam-width and the on-axis intensity profile are shown (red dot markers) in Fig. 10(c,d), respectively. Clearly the resulting intensity profile and beam-width deviate substantially for the desired ones (black, solid). In order to improve the result, we update the aperture field according to the desired on-axis profile and beam-width, following the procedures in 3.3 and 3.4. The results for the updated aperture fields are shown in Fig. 10(b) (for the -3dB iso-surface) and in Fig. 10(c,d) (red, dotted lines). Clearly, in both beam-width and intensity profiles, the updated aperture distribution meets the design requirement better.
In order to quantify the improved performance of the updated design in the on- and off-axis profiles, we use a criterion similar to the one introduced in [32], i.e., the deviation from the desired profile is assessed via an index, defined as the standard deviations of the obtained parameter minus the desired one. The computed off-axis and on-axis index are $2.87\lambda$ and $1.05$, respectively for the initial CB, and $0.39\lambda$ and $0.03$, for the updated CB. These numbers and Fig. 10, indicate the significant control capabilities over the desired features of the methods.
5.3 Conical-helical beam
For the last example, we chose a CB with a conical-helical beam-axis (a helix of gradually decreasing radius). This example sheds light over the aperture field structure obtained for complex ${\mathbf {r}}'_{\textrm {a}}$ curves. The desired beam-axis is given by
The curve ${\mathbf {r}}'_{\textrm {a}}(\sigma ')$ on the aperture, corresponding to the beam-axis in (27), forms closed loops (see Fig. 11(a)). Here, we extend the approach proposed in [32] for the areas close to the intersection points. For these areas, two sets of rays need to be synthesized, each forming different sections of the caustic surface (see Fig. 11(b)). Therefore, the aperture field is a sum of ray fields of the form in (2). In the iterative process of adjusting the on-axis intensity profile or the off-axis width, we apply the procedure to each of the terms separately, according to the mapping discussed after (3). In this example, we choose $W_{\textrm {t}}(\sigma ')=20\lambda$ and $I(\sigma )=1$. The resulting -3dB iso-surface of the CB intensity after applying the on- axis and off-axis procedure is plotted in Fig. 12.
6. Conclusions
In this paper, we introduced a practical algorithm for the construction of aperture field distributions that correspond to 3D-CBs. The procedure involves the construction of the aperture field phase and amplitude, and application of an iterative algorithm to manipulate and control the intensity profile along the beam-axis as well as the beam-width profile. Numerical examples demonstrate the algorithm’s capability to synthesize aperture field distributions that give rise to beams with arbitrary curvature and torsion. Along these trajectories, control of the on-axis intensity profile and beam-width in the bi-normal direction was demonstrated. The representation of the aperture and measured fields as samples on uniform Cartesian grids allowed for the convenient acceleration of the computation of exact field integrals. This was used for both validation of the approach and as a part of the iterative improvement of the design. A preliminary study shows that the sensitivity, in terms of beam-width deviations from the design specifications, to the aperture sampling rate, is similar to that acceptable in antenna array design. However, if only the trajectory formed by the peak intensity is of interest, the restrictions on the aperture sampling can be relaxed. This suggests that an experimental validation set-up should include the lowest pixel-pitch SLM available to date [41]. Future work may involve the design of optimization schemes to improve the control of the beam parameters manipulated in this work, i.e., the beam-width in the $n_{\textrm {b}}$ local direction and the intensity profile as well as for gaining control over the beam-width in the $n$ local direction.
Funding
Israel Science Foundation (531/19).
Acknowledgments
This research was supported by the ISRAEL SCIENCE FOUNDATION (grant No. 531/19).
Disclosures
The authors declare no conflicts of interest.
References
1. G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Observation of accelerating airy beams,” Phys. Rev. Lett. 99(21), 213901 (2007). [CrossRef]
2. G. A. Siviloglou and D. N. Christodoulides, “Accelerating finite energy airy beams,” Opt. Lett. 32(8), 979–981 (2007). [CrossRef]
3. I. M. Besieris and A. M. Shaarawi, “A note on an accelerating finite energy airy beam,” Opt. Lett. 32(16), 2447–2449 (2007). [CrossRef]
4. G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Ballistic dynamics of airy beams,” Opt. Lett. 33(3), 207–209 (2008). [CrossRef]
5. M. Bandres, “Accelerating beams,” Opt. Lett. 34(24), 3791–3793 (2009). [CrossRef]
6. L. Froehly, F. Courvoisier, A. Mathis, M. Jacquot, L. Furfaro, R. Giust, P. A. Lacourt, and J. M. Dudley, “Arbitrary accelerating micron-scale caustic beams in two and three dimensions,” Opt. Express 19(17), 16455–16465 (2011). [CrossRef]
7. Y. Kaganovsky and E. Heyman, “Wave analysis of airy beams,” Opt. Express 18(8), 8440–8452 (2010). [CrossRef]
8. Y. Hu, G. A. Siviloglou, P. Zhang, N. K. Efremidis, D. N. Christodoulides, and Z. Chen, “Self-accelerating airy beams: Generation, control, and applications,” in Nonlinear Photonics and Novel Optical Phenomena, Z. Chen and R. Morandotti, eds. (Springer New York, New York, NY, 2012), pp. 1–46.
9. Y. Kaganovsky and E. Heyman, “Nonparaxial wave analysis of three-dimensional airy beams,” J. Opt. Soc. Am. A 29(5), 671–688 (2012). [CrossRef]
10. B. K. Singh, R. Remez, Y. Tsur, and A. Arie, “Measurement of acceleration and orbital angular momentum of airy beam and airy-vortex beam by astigmatic transformation,” Opt. Lett. 40(22), 5411–5414 (2015). [CrossRef]
11. J. Broky, G. A. Siviloglou, A. Dogariu, and D. N. Christodoulides, “Self-healing properties of optical airy beams,” Opt. Express 16(17), 12880–12891 (2008). [CrossRef]
12. P. Zhang, Y. Hu, T. Li, D. Cannan, X. Yin, R. Morandotti, Z. Chen, and X. Zhang, “Nonparaxial mathieu and weber accelerating beams,” Phys. Rev. Lett. 109(19), 193901 (2012). [CrossRef]
13. A. Aiello, G. S. Agarwal, M. Paúr, B. Stoklasa, Z. Hradil, J. Řeháček, P. de la Hoz, G. Leuchs, and L. L. Sánchez-Soto, “Unraveling beam self-healing,” Opt. Express 25(16), 19147–19157 (2017). [CrossRef]
14. J. Baumgartl, M. Mazilu, and K. Dholakia, “Optically mediated particle clearing using airy wavepackets,” Nat. Photonics 2(11), 675–678 (2008). [CrossRef]
15. D. N. Christodoulides, “Optical trapping: Riding along an airy beam,” Nat. Photonics 2(11), 652–653 (2008). [CrossRef]
16. Z. Zheng, B.-F. Zhang, H. Chen, J. Ding, and H.-T. Wang, “Optical trapping with focused airy beams,” Appl. Opt. 50(1), 43–49 (2011). [CrossRef]
17. R. Schley, I. Kaminer, E. Greenfield, R. Bekenstein, Y. Lumer, and M. Segev, “Loss-proof self-accelerating beams and their use in non-paraxial manipulation of particles’ trajectories,” Nat. Commun. 5(1), 5189 (2014). [CrossRef]
18. W. Liu, X. Yang, and J. Gao, “Optical transportation and accumulation of microparticles by self-accelerating cusp beams,” Phys. Rev. A 99(4), 043839 (2019). [CrossRef]
19. S. Jia and X. Zhuang, “Super-resolution imaging with airy beams,” Opt. Photonics News 25, 35 (2014).
20. S. Jia, J. C. Vaughan, and X. Zhuang, “Isotropic 3d super-resolution imaging with a self-bending point spread function,” Nat. Photonics 8(4), 302–306 (2014). [CrossRef]
21. D. Kuang, Y. Cao, T. Lepine, and W. Mi, “Curved surface plasmon polariton excitation with shaped beam by fifth-power phase mask,” IEEE Photonics J. 7(6), 1–5 (2015). [CrossRef]
22. A. E. Minovich, A. E. Klein, D. N. Neshev, T. Pertsch, Y. S. Kivshar, and D. N. Christodoulides, “Airy plasmons: non-diffracting optical surface waves,” Laser Photonics Rev. 8(2), 221–232 (2014). [CrossRef]
23. T. Vettenburg, H. I. C. Dalgarno, J. Nylk, C. Coll-Llado, D. E. K. Ferrier, T. Cizmar, F. J. Gunn-Moore, and K. Dholakia, “Light-sheet microscopy using an Airy beam,” Nat. Methods 11(5), 541–544 (2014). [CrossRef]
24. M. Zhang, Z. Ren, and P. Yu, “Improve depth of field of optical coherence tomography using finite energy airy beam,” Opt. Lett. 44(12), 3158–3161 (2019). [CrossRef]
25. P. Polynkin, M. Kolesik, J. V. Moloney, G. A. Siviloglou, and D. N. Christodoulides, “Curved plasma channel generation using ultraintense airy beams,” Science 324(5924), 229–232 (2009). [CrossRef]
26. Y. Kaganovsky and E. Heyman, “Airy pulsed beams,” J. Opt. Soc. Am. A 28(6), 1243–1255 (2011). [CrossRef]
27. P. Saari, “Laterally accelerating airy pulses,” Opt. Express 16(14), 10303–10308 (2008). [CrossRef]
28. I. Epstein, R. Remez, Y. Tsur, and A. Arie, “Generation of intensity-controlled two-dimensional shape-preserving beams in plasmonic lossy media,” Optica 3(1), 15–19 (2016). [CrossRef]
29. I. D. Chremmos, Z. Chen, D. N. Christodoulides, and N. K. Efremidis, “Bessel-like optical beams with arbitrary trajectories,” Opt. Lett. 37(23), 5003–5005 (2012). [CrossRef]
30. J. Zhao, P. Zhang, D. Deng, J. Liu, Y. Gao, I. D. Chremmos, N. K. Efremidis, D. N. Christodoulides, and Z. Chen, “Observation of self-accelerating bessel-like optical beams along arbitrary trajectories,” Opt. Lett. 38(4), 498–500 (2013). [CrossRef]
31. S. Vo, K. Fuerschbach, K. P. Thompson, M. A. Alonso, and J. P. Rolland, “Airy beams: a geometric optics perspective,” J. Opt. Soc. Am. A 27(12), 2574–2582 (2010). [CrossRef]
32. T. Melamed and A. Shlivinski, “Practical algorithm for custom-made caustic beams,” Opt. Lett. 42(13), 2499–2502 (2017). [CrossRef]
33. Y. Wen, Z. Liu, S. Lin, Y. Chen, Y. Zhang, and S. Yu, “Construction, characteristics, and constraints of accelerating beams based on caustic design,” Opt. Express 26(25), 32728–32738 (2018). [CrossRef]
34. R.-S. Penciu, V. Paltoglou, and N. K. Efremidis, “Closed-form expressions for nonparaxial accelerating beams with pre-engineered trajectories,” Opt. Lett. 40(7), 1444–1447 (2015). [CrossRef]
35. M. Goutsoulas and N. Efremidis, “Precise amplitude, trajectory, and beam-width control of accelerating and abruptly autofocusing beams,” Phys. Rev. A 97(6), 063831 (2018). [CrossRef]
36. M. Born and E. Wolf, Principles of Optics (Pergamon Press, New York, 1964).
37. V. Babič and V. Buldyrev, Short-Wavelength Diffraction Theory: Asymptotic Methods (Spring–Verlag, Berlin, 1991).
38. S. Karimkashi and A. A. Kishk, “Focused microstrip array antenna using a dolph-chebyshev near-field design,” IEEE Trans. Antennas Propag. 57(12), 3813–3820 (2009). [CrossRef]
39. A. F. Peterson, S. L. Ray, and R. Mittra, Computational methods for electromagnetics, vol. 24, Chapter 4.12 (IEEE press, New York, 1998).
40. Z. Pi, Y. Hu, Z. Chen, and J. Xu, “Large-scale sharply bending paraxial beams,” APL Photonics 4(5), 056101 (2019). [CrossRef]
41. P. K. Shrestha, Y. T. Chun, and D. Chu, “A high-resolution optically addressed spatial light modulator based on zno nanoparticles,” Light: Sci. Appl. 4(3), e259 (2015). [CrossRef]