Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Propagation dynamics of Janus vortex waves

Open Access Open Access

Abstract

We study the propagation dynamics of Janus vortex wave under the action of a focusing lens based upon the formula of focused circular vortex Airy beams. Two dark foci would be generated under the action of a lens, and thus a perfect light hollow bottle could be formed. By controlling corresponding parameters, we could control the focal position and the relative intensity between the two focal intensities. The off-axis optical vortex (OV) would rotate rapidly in two focal regions, but keep still in the lens focus region. The angular displacement of OV in each focusing process is nearly π/2. (Note that the angular displacement for an off-axis OV in single focusing process of Gaussian beam is nearly π.) Two same OVs would repel to each other, but two opposite OVs would attract each other and annihilate at first focus plane in Janus vortex waves.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

The Janus waves were proposed in 2016 [1]. These light waves could generate two distinct foci under the action of a focusing lens, thus have potential great application in many fields. The circular Airy beam (CAB) is one of the most common members in the family of Janus waves [1]. It’s well known that CAB could abruptly autofocus at a real focus in + z direction [2], but in fact, it still has a virtual focus in -z direction. So, under the action of the lens, two new real images could be generated by objects located at “real focus” and “virtual focus”. This property is called dual focus property in this paper. The abruptly autofocusing property of CAB or other similar beams are widely investigated in these years [315]. But there are few studies about the dual focus property of CAB [16,17]. Recently, Papazoglou’s group demonstrate that this property could be applied in strong field science for phase memory preserving [18].

Optical vortex (OV) has attracted intensive interest in last decades. Many interesting phenomena and important applications about OV have been investigated, such as optical rotation, confinement of cold atoms, orbital angular momentum of photons and so on [1925]. The OV imposed in different kinds of special light beams, such as Bessel beams [26], Airy beams [2731], partially coherent beams [3234], perfect light ring [3537], circle Pearcey beams [14], have been widely investigated. These studies have greatly expanded the application of these special light beams; and they are also helpful for us to understand the physical behavior of OV. To our best knowledge, the OV imposed in Janus waves, and its effect on the double foci, has not been investigated yet. It is natural to ask could the beam keep its special dual focus property and how to control this property by different parameters, when OVs are imposed in. On the other hand, it is well-known that the off-axis vortex could rotate in single focus light beam; and a pair of opposite OVs would collide and annihilate exactly at the focus [5,38,39]. So, it is also interesting to investigate what will happen when this pair of OVs are imposed in dual focus light beams. Will they annihilate at the first or the second focus, or at the lens focus?

In this paper, we have investigated the propagation characteristics of Janus vortex waves in detail, based upon the expression of circular Airy vortex beams (CAVB). The propagation dynamics of CAVB through ABCD paraxial optical system are calculated by numerical Hankel transform method; the influences of different parameters on dual focus property are investigated. After that, the rotation of an off-axis OV in Janus vortex wave is studied, some interesting results are found. Finally, we investigate the case when two same or opposite OVs are imposed in CAB, the repel and attraction of two OVs in dual focus light beam are analyzed.

2. Theoretical simulation of propagation of Janus vortex wave

The CAB is demonstrated as a type of Janus wave [1]. The CAB with an OV located at (rk, φk) in polar coordinates could be expressed as [2,5],

$${u_0}({r,\varphi ,z = 0} )= G \cdot Ai\left( {\frac{{{r_0} - r}}{w}} \right)\exp \left( {a\frac{{{r_0} - r}}{w}} \right) \cdot {({r{e^{i\varphi }} - {r_k}{e^{i{\varphi_k}}}} )^n},$$
where r0 is a parameter related with the initial radius, w is a scaling parameter, a is a decaying parameter, n is the topological charge, G is a constant related with the incident power. The propagation of CAVB through any ABCD paraxial optical system could be calculated by Collins formula [40],
$$u({x,y,z} )= \frac{{{e^{ikz}}}}{{i\lambda B}}\int\!\!\!\int {{u_0}({{x_0},{y_0}} )\exp \left\{ {\frac{{ikA}}{{2B}}\left[ {({x_0^2 + y_0^2} )+ \frac{D}{A}({{x^2} + {y^2}} )- \frac{2}{A}({x{x_0} + y{y_0}} )} \right]} \right\}d{x_0}d{y_0}}$$
$$\begin{array}{l} \mbox{ = }\frac{{{e^{ikz}}}}{{i\lambda B}}\exp \left[ {\frac{{ik}}{{2B}}\left( {D - \frac{1}{A}} \right)({{x^2} + {y^2}} )} \right]\\ \mbox{ } \times \int\!\!\!\int {{u_0}({{x_0},{y_0}} )\exp \left\{ {\frac{{ik}}{{2AB}}[{{{({A{x_0} - x} )}^2} + {{({A{y_0} - y} )}^2}} ]} \right\}d{x_0}d{y_0}} \end{array}$$
$$\begin{array}{l} \mbox{ = }\frac{{{e^{ikz}}}}{{i\lambda B}}\exp \left[ {\frac{{ik}}{{2B}}\left( {D - \frac{1}{A}} \right)({{x^2} + {y^2}} )} \right]\\ \mbox{ } \times \int\!\!\!\int {\frac{1}{{{A^2}}}{u_0}\left( {\frac{{{x_t}}}{A},\frac{{{y_t}}}{A}} \right)\exp \left\{ {\frac{{ik}}{{2AB}}[{{{({{x_t} - x} )}^2} + {{({{y_t} - y} )}^2}} ]} \right\}d{x_t}d{y_t}} , \end{array}$$
where xt=Ax0; and (x,y) and (x0,y0) are Cartesian coordinates of the output plane and input plane, respectively. By using the convolution theory of Fourier transform, Eq. (2) could be written as,
$$u(x,y,z) = \mbox{exp}\left[ {\frac{{ikC}}{{2A}}({{x^2} + {y^2}} )} \right] \times F{T^{ - 1}}[{{{\tilde{u}}_0}(A{f_x},A{f_y}) \cdot prop} ],$$
where k is the wave number, fx and fy are spatial frequency, ${\tilde{u}_0} = FT({u_0})$ is the Fourier transform of light fields in the initial plane, FT−1 denotes the inverse Fourier transform; prop is the famous propagating factor, whose formula in ABCD system in Eq. (3) should be,
$$prop = FT\left( {\frac{{{e^{ikz}}}}{{i\lambda B}}\exp \left[ {\frac{{ik}}{{2AB}}({{x^2} + {y^2}} )} \right]} \right) = \frac{{{e^{ikz}}}}{A}\exp [{ - iBA\lambda \pi f_r^2} ],$$
where ${f_r} = \sqrt {f_x^2 + f_y^2} $ is the radial spatial frequency. When A = 1, B = z, we can get the propagating factor in the free space.

Usually, the split-step Fourier transform method can be applied to numerically calculate Eq. (3) [41]. However, this method is quite time-consuming and requires a large computer memory, especially when some lenses are included in the optical system; sometimes the result would be inaccurate. Fortunately, the input light fields in our study is radially symmetric (or could be expanded to the sum of radially symmetric fields), so the 2D Fourier transform in Eq. (3) can be simplified to Hankel transform. Thus, a simple and accurate quasi-discrete Hankel transform method could be applied here [42,43]. However, we should extend this method to ABCD optical system first.

According to the binomial theorem, the initial electric fields u0 could be expanded to the sum of radially symmetric fields ${u_0} = \sum\limits_{l = 0}^n {{P_l}{u_{0l}}} $, where,

$${P_l} = G{({ - 1} )^{n - l}}\left( {\begin{array}{{c}} n\\ l \end{array}} \right){r_k}^{n - l}{e^{i({n - l} ){\varphi _k}}},$$
$${u_{0l}}({r,\varphi } )= Ai\left( {\frac{{{r_0} - r}}{w}} \right)\exp \left( {a\frac{{{r_0} - r}}{w}} \right){r^l}{e^{il\varphi }}.$$
For a radially symmetric light field, u0l=g(r)exp(ilφ), Eq. (3) could be expressed as [43,44],
$$\tilde{g}({{f_r}} )= 2\pi \int\limits_0^\infty {g\left( {\frac{r}{{|A |}}} \right){J_l}({i2\pi {f_r}r} )rdr} ,$$
$${\tilde{u}_{\mbox{0}l}}(A{f_x},A{f_y}) = \tilde{g}({{f_r}} )\exp \left[ {il\left( {\varphi + \frac{{1 - sign(A )}}{2}\pi } \right)} \right],$$
$${u_l}({r,\varphi ,z} )= 2\pi \exp \left[ {il\left( {\varphi + \frac{{1 - sign(A )}}{2}\pi } \right)} \right]\int\limits_0^\infty {\tilde{g}({{f_r}} ){J_l}({i2\pi {f_r}r} )\cdot prop \cdot d{f_r}} .$$
Note that Eqs. (7) and (9) could not be applied when A = 0. When A is 0, the output light beam is just the Fourier transform of the initial light field, and could be easily obtained from Eq. (7). Finally, the light field in the output plane could be calculated by,
$$u({r,\varphi ,z} )= \sum\limits_{l = 0}^n {{P_l}{u_l}({r,\varphi ,z} )} .$$
And the output intensity in this paper is I=|u|2.

The thin lens system is shown in Fig. 1, its optical matrix is,

$$\left[ {\begin{array}{{cc}} A&B\\ C&D \end{array}} \right] = \left[ {\begin{array}{{cc}} {1 - \frac{z}{f}}&{ - \frac{{zs}}{f} + s + z}\\ { - \frac{1}{f}}&{1 - \frac{s}{f}} \end{array}} \right],$$
where, f is the focal length of the lens, z is the distance between the output plane and the lens, s is the distance between the input plane and the lens.

 figure: Fig. 1.

Fig. 1. The CAVB with + 1 OV in the input plane (a) and the lens system (b). The focus of the lens is f, f1 and f2 denote the two foci of focused CAVB.

Download Full Size | PDF

Substituting Eq. (1) and Eq. (11) into Eqs. (7)–(10), we can obtain the propagation dynamics of Janus vortex waves, and analyze its dual focus property. In our simulation, we choose r0=1 mm, w = 0.2 mm, a = 0.1, f = 0.3 m, s = 0 unless stated otherwise; the wavelength is 632.8 nm, and the incident power is 10 mW.

The dual focus property of common CAB under the action of a thin lens is shown in Fig. 2(a). The position where light intensity reaches a peak value and the radius of beam size reaches a minimum value is defined as focus position. It is obvious that, in Fig. 2(a), two foci could be found at z = f1=0.26 m and z = f2=0.35 m, respectively. When an on-axis OV with different topological charge n is imposed in, the CAVB will keep the dual focus property [Figs. 2(d) and 2(g)]. As shown in Figs. 2(d)–2(i), two dark foci could be observed. In each dark focal plane, the intensity increases to a maximum value and the radius of the light ring decreases to a minimum value. Thus, a perfect light hollow bottle could be formed, which might be useful in optical micromanipulation [45,46].

 figure: Fig. 2.

Fig. 2. Propagation dynamics of focused CAVB (the first column), the intensity pattern at the first focus (the second column) and the second focus (the third column). (a)–(c) n = 0; (d)–(f) n = 1; (g)–(i) n = 2.

Download Full Size | PDF

The CAVBs with different n almost focus at two same regions as CAB, if other parameters are not changed. But in fact, as shown in Fig. 3(a), the accurate focus positions change with n. As n increases, the distance between two foci |f2-f1| would be slightly shortened: f1 becomes larger, and f2 becomes smaller. And the focal intensities at two foci also decrease with n. The initial radius r0 can also influence the dual focus property. As shown in Fig. 3(b), as r0 increases, the distance between two dark foci would also be shortened, but unlike the influence of n mentioned before, here, the focal intensity would increase as |f2-f1| decreases.

 figure: Fig. 3.

Fig. 3. Changes of the maximum intensity along z axis with different parameters: (a) with different n, when r0=1 mm, s = 0; (b) with different r0, when n = 1, s = 0; (c) with different s, when n = 1, r0=1 mm.

Download Full Size | PDF

From the imaging theory of a thin lens, it is known that positions of two foci of Janus wave could be controlled by varying the distance between the input plane and the lens s. As s increases, the virtual (or real) object distance would decrease (or increase), leading to the decrease of the real image distance f1 (or f2). The same conclusion could be found in [1] and Fig. 3(c). It is also worth noting that the relative intensity at two foci could be effectively controlled by varying s. Usually, when s = 0, the intensity at first focus is much larger than the second [Figs. 3(a)–3(c)]; but as s increases to 0.3 m, the two focal intensities are nearly the same; when s is 0.6 m, the focal intensity at the second focus exceeds the first one [Fig. 3(c)]. This is of great use for the application of Janus vortex waves in the dual-trap optical tweezes [47,48]. By changing s, we can easily control the trapping force of each trap.

3. Rotation of an off-axis OV in Janus vortex wave

It is demonstrated that an off-axis OV would be forced to move to the center and rotated by some angle in the single-focusing process [5]. It is interesting to investigate the behavior of off-axis OV in the dual-focusing process. Figure 4 shows that the CAVB could keep its dual focus property when an off-axis OV is located at rk=0.5mm: the radial coordinate of OV reaches minimum values in two focus planes [Figs. 4(c) and 4(g)]; and intensity reaches maximum values in two focal planes (Fig. 5). Similar with the case of on-axis OV, the two foci for the beam with off-axis OV would also become closer with a larger topological charge n. In the two focal regions, the positive OV would rotate anticlockwise around the beam axis. Similarly, the negative OV would rotate clockwise. In other positions, the off-axis OV would keep nearly at the same angle position, only the radial position is changed.

 figure: Fig. 4.

Fig. 4. Intensity patterns [(a)–(h)] and phase patterns [(a1)–(h1)] of focused CAVB at different z position. (a) z = 0; (b) z = 0.24 m; (c) first focus plane, z = 0.26 m; (d) z = 0.28 m; (e) lens focus plane, z = 0.3 m; (f) z = 0.32 m; (g) second focus plane, z = 0.35 m; (h) z = 0.4 m. The inset arrows in the third column denote the rotation direction of OV.

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. Changes of maximum intensity along z axis for the focused CAVB with an off-axis OV.

Download Full Size | PDF

The positions of the off-axis OV could be accurately detected from phase patterns at different output planes. The trajectory of the off-axis OV in focused CAVB is shown in Fig. 6, and is compared with the single focus light beam (i.e., focused Gaussian beams). For Gaussian beam (GB), the rotation of OV is greatly accelerated in the focal region of the lens (z = f), but the OV in CAVB almost keeps still in this region. The OV in CAVB rotates rapidly in its two own focal regions (z = f1 and z = f2). Moreover, for GB, the OV could not finish even one full rotation, it slowly rotates to the angle position of π after the focus plane. But for CAVB, in our case, the OV would finish 4 complete rotation at each focal plane. At the first focal region, the OV rotates from 0 to a final position of π/2; at the second focal region, from π/2 to a final position of π. So, the OV for both beams would finally reach a same angle position when the propagation distance is large enough.

 figure: Fig. 6.

Fig. 6. Comparison of angle position of OV in dual focus light beams (i.e., focused CAVB) and in single focus light beams (i.e., focused GB).

Download Full Size | PDF

4. Propagation dynamics of Janus vortex wave with two OVs

For the single-focus beam, it is proved that the two OV could repel or attract each other, when they are very close to each other in the focal region [5,38]. Similar conclusions could also be found in dual-focus light beams. We assume the two OVs are located symmetrically in x-axis, and we only consider n = 1 here for simplicity,

$${u_0}({r,\varphi ,z = 0} )= G \cdot Ai\left( {\frac{{{r_0} - r}}{w}} \right)\exp \left( {a\frac{{{r_0} - r}}{w}} \right)({r{e^{i\varphi }} - {r_k}} )({r{e^{ {\pm} i\varphi }}\mbox{ + }{r_k}} ).$$
Similar with the calculation for the beam with an off-axis OV, by substituting the expansion of Eq. (12) into Eqs. (7)–(10), we can get the intensity patterns and phase patterns at different positions. Here we use (+1, +1) to denote two same OVs with charge + 1; and (+1, −1) for two opposite OVs.

As shown in Fig. 7, when two same OVs (+1, +1) are imposed on, the beam would keep hollow in the whole propagation process; and two OVs would rotate around each other in two focal regions. The two same OVs would move closer in two focal regions, but they cannot completely overlap each other, thus cannot combine into a charge-2 OV under paraxial condition. The same conclusion could also be found in Fig. 8, where rOV denotes the distance between + 1 OV and beam axis. Because of the dual focusing effect, rOV reaches minimum values in two focal regions, and becomes larger in lens focal region and other positions. Comparing with the case of single off-axis OV, the repelling effect of two same OV in focused CAVB is obvious, especially in lens focal region. As shown in Fig. 8, for the single + 1 OV, the distance between OV and beam center is 28µm in lens focal plane z = f; but for the case of two same OVs, this distance would increase to 38µm because of the repelling effect. So, this pair of OVs do not just rotate rigidly as nested in GB [38].

 figure: Fig. 7.

Fig. 7. Intensity profiles [(a)–(f)] and phase patterns [(a1)–(f1)] when two same OVs (+1, +1) are imposed in focused CAVB. (a) z = 0; (b) z = 0.2 m; (c) first focus plane, z = 0.26 m; (d) lens focus plane, z = 0.3 m; (e) second focus plane, z = 0.35 m; (f) z = 0.4 m.

Download Full Size | PDF

 figure: Fig. 8.

Fig. 8. Changes of the distance between OV and beam axis rOV with z in dual focusing process, rk=0.5 mm for both cases.

Download Full Size | PDF

Figure 9 shows the intensity patterns and phase patterns of focused CAVB with two opposite OVs (+1, −1) in different output planes. As we can see in this figure, the two opposite OVs would rotate in different direction, attract each other and annihilate exactly in the first focal plane z = f1 [Fig. 9(c1)]. Because of the disappearance of OV, the focal intensity profile would become Gaussian-like [Fig. 9(c)]. The beam could keep the nonhollow profile in the lens focus plane [Fig. 9(d)] and the second focus plane [Fig. 9(e)]. After the second focal plane, the beam would become hollow and unsymmetrical again [Fig. 9(f)], but the two opposite OVs would never reappear [Fig. 9(f1)]. Because of the collision of two opposite OVs, the beam with (+1, −1) OVs shows similar dual focus property as CAB without OV [ Fig. 10(a)]. However, because of the introduction of OVs, the beam with opposite OVs shows some interesting differences. As shown in Fig. 10(b), when two opposite OVs are imposed in, the axial focal widths would be broadened, and the intensity in lens focal region would be greatly reduced.

 figure: Fig. 9.

Fig. 9. Intensity profiles [(a)–(f)] and phase patterns [(a1)–(f1)] when two opposite OVs (−1, +1) are imposed in focused CAVB. (a) z = 0; (b) z = 0.2 m; (c) first focus plane, z = 0.26 m; (d) lens focus plane, z = 0.3 m; (e) second focus plane z = 0.35 m; (f) z = 0.4 m.

Download Full Size | PDF

 figure: Fig. 10.

Fig. 10. Propagation dynamics of focused CAVB with (+1, −1) OVs in y-z plane (a) and the comparison of intensity distribution along beam axis with focused CAB without OV (b).

Download Full Size | PDF

5. Conclusions

The dual focus property of Janus vortex waves is investigated in this paper based upon the expression of focused CAVB. A convenient and accurate numerical method based upon Hankel transform is presented to simulate the propagation dynamics of light beams through ABCD optical system. By using this method, the propagation characteristics of focused CAVB with an on-axis OV, an off-axis OV and a pair of OVs are investigated in detail.

When an on-axis OV is imposed in Janus vortex waves, two dark foci could be generated, thus a perfect light hollow bottle could be formed. This hollow bottle might be useful in optical micromanipulation. As the topological charge n or the initial radius r0 increases, the two dark foci would become closer. The focal intensities at two foci would decrease with n, but increase with r0. By varying the distance between the input plane and the thin lens s, the focal positions would be changed. And most importantly, the relative intensity between the two focal intensities could be controlled by s. When s is small, the intensity at first focus is larger; but as s increases, the second focal intensity would increase and even exceeds the first focal intensity. These modulation methods might be of great help for the application of Janus vortex wave in dual-trap optical tweezers or other fields.

The off-axis OV would rotate rapidly in two focal regions of Janus vortex waves, but keeps nearly still in the lens focal region. The off-axis OV in Janus vortex wave shows interesting behaviors when compared with the common GB. The off-axis OV in GB could not finish even one complete rotation, it slowly reaches to the π position after the focus. But the OV in CAVB experiences several complete rotations in two focal regions; and in each focusing process, the angular displacement of OV is nearly π/2. Two same OVs would rotate around each other and exhibit repelling effect in focusing process. Two opposite OVs would annihilate exactly at the first focus, and the beam would become nonhollow. After the second focus, the beam would become hollow again, but the OV would never reappear.

We believe that the Janus vortex wave proposed in this paper is helpful to understand the interesting behavior of OV; and could also be used in optical micromanipulation, optics communication and other fields..

Funding

Natural Science Foundation of Zhejiang Province (LY20A040006); National Natural Science Foundation of China (11504274).

Disclosures

The authors declare no conflicts of interest.

References

1. D. G. Papazoglou, V. Y. Fedorov, and S. Tzortzakis, “Janus waves,” Opt. Lett. 41(20), 4656–4659 (2016). [CrossRef]  

2. N. K. Efremidis and D. N. Christodoulides, “Abruptly autofocusing waves,” Opt. Lett. 35(23), 4045–4047 (2010). [CrossRef]  

3. D. G. Papazoglou, N. K. Efremidis, D. N. Christodoulides, and S. Tzortzakis, “Observation of abruptly autofocusing waves,” Opt. Lett. 36(10), 1842–1844 (2011). [CrossRef]  

4. J. A. Davis, D. M. Cottrell, and D. Sand, “Abruptly autofocusing vortex beams,” Opt. Express 20(12), 13302–13310 (2012). [CrossRef]  

5. Y. Jiang, K. Huang, and X. Lu, “Propagation dynamics of abruptly autofocusing Airy beams with optical vortices,” Opt. Express 20(17), 18579–18584 (2012). [CrossRef]  

6. S. Liu, M. Wang, P. Li, P. Zhang, and J. Zhao, “Abrupt polarization transition of vector autofocusing Airy beams,” Opt. Lett. 38(14), 2416–2418 (2013). [CrossRef]  

7. P. Panagiotopoulos, D. G. Papazoglou, A. Couairon, and S. Tzortzakis, “Sharply autofocused ring-Airy beams transforming into non-linear intense light bullets,” Nat. Commun. 4(1), 2622 (2013). [CrossRef]  

8. N. Li, Y. Jiang, K. Huang, and X. Lu, “Abruptly autofocusing property of blocked circular Airy beams,” Opt. Express 22(19), 22847–22853 (2014). [CrossRef]  

9. P. Li, S. Liu, T. Peng, G. Xie, X. Gan, and J. Zhao, “Spiral autofocusing Airy beams carrying power-exponent-phase vortices,” Opt. Express 22(7), 7598–7606 (2014). [CrossRef]  

10. B. Chen, C. Chen, X. Peng, Y. Peng, M. Zhou, and D. Deng, “Propagation of sharply autofocused ring Airy Gaussian vortex beams,” Opt. Express 23(15), 19288–19298 (2015). [CrossRef]  

11. Y. Jiang, X. Zhu, W. Yu, H. Shao, W. Zheng, and X. Lu, “Propagation characteristics of the modified circular Airy beam,” Opt. Express 23(23), 29834–29841 (2015). [CrossRef]  

12. Z. Zhao, C. Xie, D. Ni, Y. Zhang, Y. Li, F. Courvoisier, and M. Hu, “Scaling the abruptly autofocusing beams in the direct-space,” Opt. Express 25(24), 30598–30605 (2017). [CrossRef]  

13. D. Amaya, Ó Martínez-Matos, and P. Vaveliuk, “Abruptly autofocusing beams from phase perturbations having forced symmetry,” Opt. Lett. 44(15), 3733–3736 (2019). [CrossRef]  

14. X. Chen, D. Deng, G. Wang, X. Yang, and H. Liu, “Abruptly autofocused and rotated circular chirp Pearcey Gaussian vortex beams,” Opt. Lett. 44(4), 955–958 (2019). [CrossRef]  

15. N. K. Efremidis, Z. Chen, M. Segev, and D. N. Christodoulides, “Airy beams and accelerating waves: an overview of recent advances,” Optica 6(5), 686–701 (2019). [CrossRef]  

16. M. Hashemi, M. M. Shanei, and C. J. Zapata-Rodríguez, “Flexible design of multifocal metalenses based on autofocused Airy beams,” J. Opt. Soc. Am. B 35(11), 2864–2868 (2018). [CrossRef]  

17. V. Y. Fedorov, D. G. Papazoglou, and S. Tzortzakis, “Transformation of ring-Airy beams during efficient harmonic generation,” Opt. Lett. 44(12), 2974–2977 (2019). [CrossRef]  

18. A. D. Koulouklidis, D. G. Papazoglou, V. Y. Fedorov, and S. Tzortzakis, “Phase Memory Preserving Harmonics from Abruptly Autofocusing Beams,” Phys. Rev. Lett. 119(22), 223901 (2017). [CrossRef]  

19. S. Sato, M. Ishigure, and H. Inaba, “Optical trapping and rotational manipulation of microscopic particles and biological cells using higher-order mode Nd:YAG laser beams,” Electron. Lett. 27(20), 1831–1832 (1991). [CrossRef]  

20. K. Gahagan and G. A. Swartzlander, “Optical vortex trapping of particles,” Opt. Lett. 21(11), 827–829 (1996). [CrossRef]  

21. D. G. Grier, “A revolution in optical manipulation,” Nature 424(6950), 810–816 (2003). [CrossRef]  

22. J. E. Curtis and D. G. Grier, “Structure of Optical Vortices,” Phys. Rev. Lett. 90(13), 133901 (2003). [CrossRef]  

23. M. R. Dennis, K. O’Holleran, and M. J. Padgett, “Chapter 5 Singular Optics: Optical Vortices and Polarization Singularities,” in Progress in Optics, E. Wolf, ed. (Elsevier, 2009).

24. A. M. Yao and M. J. Padgett, “Orbital angular momentum: origins, behavior and applications,” Adv. Opt. Photonics 3(2), 161–204 (2011). [CrossRef]  

25. K. Toyoda, F. Takahashi, S. Takizawa, Y. Tokizane, K. Miyamoto, R. Morita, and T. Omatsu, “Transfer of Light Helicity to Nanostructures,” Phys. Rev. Lett. 110(14), 143603 (2013). [CrossRef]  

26. V. V. Kotlyar, A. A. Kovalev, and V. A. Soifer, “Asymmetric Bessel modes,” Opt. Lett. 39(8), 2395–2398 (2014). [CrossRef]  

27. H. Dai, Y. Liu, D. Luo, and X. Sun, “Propagation dynamics of an optical vortex imposed on an Airy beam,” Opt. Lett. 35(23), 4075–4077 (2010). [CrossRef]  

28. H. Dai, Y. Liu, D. Luo, and X. Sun, “Propagation properties of an optical vortex carried by an Airy beam: experimental implementation,” Opt. Lett. 36(9), 1617–1619 (2011). [CrossRef]  

29. Y. Izdebskaya, T. Lu, D. Neshev, and A. Desyatnikov, “Dynamics of three-Airy beams carrying optical vortices,” Appl. Opt. 53(10), B248–B253 (2014). [CrossRef]  

30. Z. Fang, Y. Chen, Y. Ren, L. Gong, R. Lu, A. Zhang, H. Zhao, and P. Wang, “Interplay between topological phase and self-acceleration in a vortex symmetric Airy beam,” Opt. Express 26(6), 7324–7335 (2018). [CrossRef]  

31. Y. Qian, Y. Shi, W. Jin, F. Hu, and Z. Ren, “Annular arrayed-Airy beams carrying vortex arrays,” Opt. Express 27(13), 18085–18093 (2019). [CrossRef]  

32. Y. Cai and S. Zhu, “Orbital angular moment of a partially coherent beam propagating through an astigmatic ABCD optical system with loss or gain,” Opt. Lett. 39(7), 1968–1971 (2014). [CrossRef]  

33. R. K. Singh, A. M. Sharma, and P. Senthilkumaran, “Vortex array embedded in a partially coherent beam,” Opt. Lett. 40(12), 2751–2754 (2015). [CrossRef]  

34. X. Peng, L. Liu, F. Wang, S. Popov, and Y. Cai, “Twisted Laguerre-Gaussian Schell-model beam and its orbital angular moment,” Opt. Express 26(26), 33956–33969 (2018). [CrossRef]  

35. A. S. Ostrovsky, C. Rickenstorff-Parrao, and V. Arrizón, “Generation of the perfect optical vortex using a liquid-crystal spatial light modulator,” Opt. Lett. 38(4), 534–536 (2013). [CrossRef]  

36. P. Vaity and L. Rusch, “Perfect vortex beam: Fourier transformation of a Bessel beam,” Opt. Lett. 40(4), 597–600 (2015). [CrossRef]  

37. M. Chen, M. Mazilu, Y. Arita, E. M. Wright, and K. Dholakia, “Dynamics of microparticles trapped in a perfect vortex beam,” Opt. Lett. 38(22), 4919–4922 (2013). [CrossRef]  

38. G. Indebetouw, “Optical Vortices and Their Propagation,” J. Mod. Opt. 40(1), 73–87 (1993). [CrossRef]  

39. D. Rozas, C. T. Law, and J. G. A. Swartzlander, “Propagation dynamics of optical vortices,” J. Opt. Soc. Am. B 14(11), 3054–3065 (1997). [CrossRef]  

40. S. A. Collins, “Lens-System Diffraction Integral Written in Terms of Matrix Optics,” J. Opt. Soc. Am. 60(9), 1168–1177 (1970). [CrossRef]  

41. J. D. Schmidt, Numerical Simulation of Optical Wave Propagation with Examples in MATLAB (SPIE, 2010).

42. L. Yu, M. Huang, M. Chen, W. Chen, W. Huang, and Z. Zhu, “Quasi-discrete Hankel transform,” Opt. Lett. 23(6), 409–411 (1998). [CrossRef]  

43. M. Guizar-Sicairos and J. C. Gutiérrez-Vega, “Computation of quasi-discrete Hankel transforms of integer order for propagating optical wave fields,” J. Opt. Soc. Am. A 21(1), 53–58 (2004). [CrossRef]  

44. J. W. Goodman, Introduction to Fourier Optics (Roberts & Company, 2005).

45. P. Zhang, Z. Zhang, J. Prakash, S. Huang, D. Hernandez, M. Salazar, D. N. Christodoulides, and Z. Chen, “Trapping and transporting aerosols with a single optical bottle beam generated by moiré techniques,” Opt. Lett. 36(8), 1491–1493 (2011). [CrossRef]  

46. H. Ye, C. Wan, K. Huang, T. Han, J. Teng, Y. Ping, and C. Qiu, “Creation of vectorial bottle-hollow beam using radially or azimuthally polarized light,” Opt. Lett. 39(3), 630–633 (2014). [CrossRef]  

47. E. Fällman and O. Axner, “Design for fully steerable dual-trap optical tweezers,” Appl. Opt. 36(10), 2107–2113 (1997). [CrossRef]  

48. X. Xu, C. Cheng, Y. Zhang, H. Lei, and B. Li, “Dual focused coherent beams for three-dimensional optical trapping and continuous rotation of metallic nanostructures,” Sci. Rep. 6(1), 29449 (2016). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (10)

Fig. 1.
Fig. 1. The CAVB with + 1 OV in the input plane (a) and the lens system (b). The focus of the lens is f, f1 and f2 denote the two foci of focused CAVB.
Fig. 2.
Fig. 2. Propagation dynamics of focused CAVB (the first column), the intensity pattern at the first focus (the second column) and the second focus (the third column). (a)–(c) n = 0; (d)–(f) n = 1; (g)–(i) n = 2.
Fig. 3.
Fig. 3. Changes of the maximum intensity along z axis with different parameters: (a) with different n, when r0=1 mm, s = 0; (b) with different r0, when n = 1, s = 0; (c) with different s, when n = 1, r0=1 mm.
Fig. 4.
Fig. 4. Intensity patterns [(a)–(h)] and phase patterns [(a1)–(h1)] of focused CAVB at different z position. (a) z = 0; (b) z = 0.24 m; (c) first focus plane, z = 0.26 m; (d) z = 0.28 m; (e) lens focus plane, z = 0.3 m; (f) z = 0.32 m; (g) second focus plane, z = 0.35 m; (h) z = 0.4 m. The inset arrows in the third column denote the rotation direction of OV.
Fig. 5.
Fig. 5. Changes of maximum intensity along z axis for the focused CAVB with an off-axis OV.
Fig. 6.
Fig. 6. Comparison of angle position of OV in dual focus light beams (i.e., focused CAVB) and in single focus light beams (i.e., focused GB).
Fig. 7.
Fig. 7. Intensity profiles [(a)–(f)] and phase patterns [(a1)–(f1)] when two same OVs (+1, +1) are imposed in focused CAVB. (a) z = 0; (b) z = 0.2 m; (c) first focus plane, z = 0.26 m; (d) lens focus plane, z = 0.3 m; (e) second focus plane, z = 0.35 m; (f) z = 0.4 m.
Fig. 8.
Fig. 8. Changes of the distance between OV and beam axis rOV with z in dual focusing process, rk=0.5 mm for both cases.
Fig. 9.
Fig. 9. Intensity profiles [(a)–(f)] and phase patterns [(a1)–(f1)] when two opposite OVs (−1, +1) are imposed in focused CAVB. (a) z = 0; (b) z = 0.2 m; (c) first focus plane, z = 0.26 m; (d) lens focus plane, z = 0.3 m; (e) second focus plane z = 0.35 m; (f) z = 0.4 m.
Fig. 10.
Fig. 10. Propagation dynamics of focused CAVB with (+1, −1) OVs in y-z plane (a) and the comparison of intensity distribution along beam axis with focused CAB without OV (b).

Equations (14)

Equations on this page are rendered with MathJax. Learn more.

u 0 ( r , φ , z = 0 ) = G A i ( r 0 r w ) exp ( a r 0 r w ) ( r e i φ r k e i φ k ) n ,
u ( x , y , z ) = e i k z i λ B u 0 ( x 0 , y 0 ) exp { i k A 2 B [ ( x 0 2 + y 0 2 ) + D A ( x 2 + y 2 ) 2 A ( x x 0 + y y 0 ) ] } d x 0 d y 0
 =  e i k z i λ B exp [ i k 2 B ( D 1 A ) ( x 2 + y 2 ) ]   × u 0 ( x 0 , y 0 ) exp { i k 2 A B [ ( A x 0 x ) 2 + ( A y 0 y ) 2 ] } d x 0 d y 0
 =  e i k z i λ B exp [ i k 2 B ( D 1 A ) ( x 2 + y 2 ) ]   × 1 A 2 u 0 ( x t A , y t A ) exp { i k 2 A B [ ( x t x ) 2 + ( y t y ) 2 ] } d x t d y t ,
u ( x , y , z ) = exp [ i k C 2 A ( x 2 + y 2 ) ] × F T 1 [ u ~ 0 ( A f x , A f y ) p r o p ] ,
p r o p = F T ( e i k z i λ B exp [ i k 2 A B ( x 2 + y 2 ) ] ) = e i k z A exp [ i B A λ π f r 2 ] ,
P l = G ( 1 ) n l ( n l ) r k n l e i ( n l ) φ k ,
u 0 l ( r , φ ) = A i ( r 0 r w ) exp ( a r 0 r w ) r l e i l φ .
g ~ ( f r ) = 2 π 0 g ( r | A | ) J l ( i 2 π f r r ) r d r ,
u ~ 0 l ( A f x , A f y ) = g ~ ( f r ) exp [ i l ( φ + 1 s i g n ( A ) 2 π ) ] ,
u l ( r , φ , z ) = 2 π exp [ i l ( φ + 1 s i g n ( A ) 2 π ) ] 0 g ~ ( f r ) J l ( i 2 π f r r ) p r o p d f r .
u ( r , φ , z ) = l = 0 n P l u l ( r , φ , z ) .
[ A B C D ] = [ 1 z f z s f + s + z 1 f 1 s f ] ,
u 0 ( r , φ , z = 0 ) = G A i ( r 0 r w ) exp ( a r 0 r w ) ( r e i φ r k ) ( r e ± i φ  +  r k ) .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.