Next Article in Journal
The Micro Topology and Statistical Analysis of the Forces of Walking and Failure of an ITAP in a Femur
Next Article in Special Issue
Micro-Hole Generation by High-Energy Pulsed Bessel Beams in Different Transparent Materials
Previous Article in Journal
A Microfluidic Device for Automated High Throughput Detection of Ice Nucleation of Snomax®
Previous Article in Special Issue
Bessel Beam: Significance and Applications—A Progressive Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generation of Complex Transverse Energy Flow Distributions with Autofocusing Optical Vortex Beams

by
Svetlana N. Khonina
1,2,*,
Alexey P. Porfirev
1,2,
Andrey V. Ustinov
1 and
Muhammad Ali Butt
2,*
1
IPSI RAS—Branch of the FSRC “Crystallography and Photonics” RAS, Molodogvardeyskaya 151, 443001 Samara, Russia
2
Department of Technical Cybernetics, Samara National Research University, MoskovskoyeShosse 34, 443086 Samara, Russia
*
Authors to whom correspondence should be addressed.
Micromachines 2021, 12(3), 297; https://doi.org/10.3390/mi12030297
Submission received: 17 February 2021 / Revised: 3 March 2021 / Accepted: 8 March 2021 / Published: 12 March 2021
(This article belongs to the Special Issue Non-diffractive Beams for the State of the Art Applications)

Abstract

:
Optical vortex (OV) beams are widely used for the generation of light fields with transverse energy flow inducing orbital motion of the nano- and microparticles in the transverse plane. Here, we present some new modifications of OV beams with autofocusing properties for shaping complex transverse energy flow distributions varying in space. The angular component of the complex amplitude of these beams is defined by the superpositions of OV beams with different topological charges. The proposed approach provides a convenient method to control the three-dimensional structure of the generated autofocusing OV beams. The control of the transverse distribution of an autofocusing beam provides a wide variety of generated fields with both rotating and periodic properties, which can be used in the field of laser manipulation and laser material processing. Thus, the obtained numerical results predict different types of motion of the trapped particles for the designed OV autofocusing beams. The experimental results agree with modeling results and demonstrate the principal possibility to shape such laser beams using spatial light modulators.

1. Introduction

Structured optical beams with given 3D distribution and characteristics are in demand in various applications, information on which can be found in research and review publications [1,2,3,4,5,6,7,8,9,10]. An interesting type of beams with controlled 3D characteristics is those with autofocusing properties. Such beams can autofocus without lenses or nonlinearity in the optical medium. The main characteristics of autofocusing, such as the curvature of the caustic trajectory, focusing sharpness, and focusing distance, are determined by the type of dependence of the input function on the radius. At present, the various types of autofocusing beams are: circular Airy beams [11,12,13,14,15], circular Pearcey beams [16,17,18], aberration beams [19], as well as a mirror [20,21] and generalized [22,23,24,25] Airy beams.
Among the main applications of autofocusing beams are optical trapping and manipulation [26,27,28] and laser material processing [29,30]. In these applications, not only is the character of autofocusing (sharpness and focus position) important, but also the transverse field distribution. For example, the control of the transverse structure of the intensity distribution allows fabricating different nano-microstructures with the desired shape using pulsed laser radiation [31]. The trapped particles move in the direction of the intensity and phase gradient [32]. One of the methods to vary this distribution is to supplement an axisymmetric autofocusing beam with amplitude and phase modulation [13,14,15,18,19,20,23,33,34,35,36], among which the most interesting is the vortex phase since it is preserved during propagation. Inserting a vortex phase singularity into a laser beam can lead to various polarization and spin-orbit transformations [37,38,39,40].
An important characteristic of the generated field for optical trapping problems is the transverse energy flow density (TEFD), which is proportional to the product of the field intensity and phase gradient [41,42]. Inserting a vortex singularity leads to variations in the phase gradient and provides an opportunity to control the TEFD [32].
In this paper, we consider fields with a complex amplitude in the form of the product of the circular Airy function and a combination of optical vortices. The first factor ensures the fact of autofocusing, and optical vortices change the transverse structure of the field. This approach is convenient for experimental implementation since separate optical elements can be used to control the autofocusing properties and the transverse beam structure.
The superposition of optical vortices allows the formation of rotating (spiral) autofocusing beams during propagation, as well as varying the TEFD in the focal area. Autofocusing beams rotating during propagation have attracted the attention of researchers [23,35,36] due to the transverse beam structure varying depending on the distance, which expands the possibilities of using such beams in various applications. Note that in these works, to obtain the rotation effect, either the spatial combination of different beams was used, for example, by sectors [35], rings [36], or inserting a nonlinear vortex phase into a beam [23]. The second approach is more convenient for experimental implementation; however, the nonlinear vortex phase provides the formation of one type of transverse distribution in the form of a spiral [43].
In this paper, we consider the control of the transverse distribution of an autofocusing beam due to the coaxial superposition of several optical vortices, which can provide a wide variety of generated fields with both rotating and periodic properties [44,45]. The results of modeling the distribution of intensity and TEFD and the experimental results of measuring the intensity are in good agreement.

2. Theoretical Analysis

Let us consider vortex fields with autofocusing properties in the following form:
f ( r , φ )   =   A ( r ) Φ ( φ )   =   A ( r ) l = 1 L exp ( i m l φ )
where A ( r ) is a real function depending only on the radius (for example, circular Airy function, circular Pearcey function, or radial sinus function) and exp ( i m l φ ) is the optical vortex with the topological charge m l .
The main characteristics of autofocusing, such as the curvature of the caustic trajectory and the focusing distance, are determined by the function A ( r ) in Equation (1). However, in applications such as microparticles trapping and laser structuring, the transverse field distribution in the focal region is also important. As a rule, variations in the transverse distribution are achieved by amplitude and phase modulation of an axisymmetric beam A ( r ) [13,14,15,18,19,20,23,33,34,35,36].
The representation in Equation (1) is convenient for experimental implementation since it is possible to use separate optical elements to control the autofocusing properties provided by the function A ( r ) and the transverse structure of the beam due to the superposition of optical vortices provided by the function Φ ( φ ) . Note that variations in Φ ( φ ) can be performed dynamically using the spatial light modulator (SLM).
Factorization of the field into radial and angular components makes it possible to separately perform a theoretical analysis of the autofocusing properties of the field and its transverse structure. Since autofocusing properties that depend on the radial part of the field A ( r ) have been investigated in previous works listed, for example, in the Introduction, in this section, we focus on the analysis of the transverse properties associated with the angular part Φ ( φ ) .
Let us consider the paraxial propagation of the beam (1) in free space at a distance z using the Kirchhoff–Fresnel integral:
E ( ρ , θ , z )   =   i k 2 π z 0 2 π 0 R f ( r , φ ) exp { i k 2 z [ ρ 2   +   r 2     2 ρ r cos ( φ     θ ) ] } r d r d φ
where k = 2π/λ is the wavenumber for laser radiation with a wavelength λ.
If we substitute Equation (1) into Equation (2), then the integration d φ can be performed analytically. For one l-term, we get:
E l ( ρ , θ , z )   =   ( i ) m l + 1 k z exp ( i m l θ ) exp ( i k 2 z ρ 2 ) 0 R A ( r ) J m l ( k ρ r z ) exp ( i k 2 z r 2 ) r d r
The full field will be:
E ( ρ , θ , z )   =   k z exp ( i k 2 z ρ 2 ) l = 1 L [ ( i ) m l + 1 exp ( i m l θ ) 0 R A ( r ) J m l ( k ρ r z ) exp ( i k 2 z r 2 ) r d r ]
Let’s rewrite Equation (4) as follows:
E ( ρ , θ , z )   =   B ( ρ , z ) l = 1 L C l ( ρ , z ) exp ( i m l θ )
where B ( ρ , z )   =   k z exp ( i k 2 z ρ 2 ) ,
C l ( ρ , z )   =   ( i ) m l + 1 0 R A ( r ) J m l ( k ρ r z ) exp ( i k 2 z r 2 ) r d r
It is seen from Equation (5) that the field defined by Equation (1), when propagating, conserves the same set of vortex functions that were at the input plane (at z = 0), but the weights of these functions change depending on the radius of the observation point ρ and the propagation distance z.
For a detailed analysis of transverse distribution, let us consider the field defined by Equation (5) in some fixed plane z   =   z 0 :
E ( ρ , θ , z 0 )   =   B ( ρ , z 0 ) l = 1 L C l ( ρ , z 0 ) exp ( i m l θ )
The TEFD in the plane z   =   z 0 is proportional to the product of the field intensity and phase gradient [32,42,43]:
F ( ρ , θ , z 0 ) | E ( ρ , θ , z 0 ) | 2 ψ ( ρ , θ , z 0 )
where ψ ( ρ , θ , z 0 )   =   arg [ E ( ρ , θ , z 0 ) ] is the phase of the field in the plane z   =   z 0 .
The field intensity defined from Equation (7) has the following form:
| E ( ρ , θ , z 0 ) | 2   =   ( k z 0 ) 2 { [ l = 1 L D l ( ρ ) cos [ m l θ + ψ l ( ρ ) ] ] 2 + [ l = 1 L D l ( ρ ) sin [ m l θ + ψ l ( ρ ) ] ] 2 }
where D l ( ρ )   =   | C l ( ρ , z 0 ) | and ψ l ( ρ )   =   arg [ C l ( ρ , z 0 ) ] .
Since we are interested in the radial and angular energy flow directions, we consider the phase gradient in polar projections:
ψ ( ρ , θ , z 0 )   =   ψ ( ρ , θ , z 0 ) ρ e ρ   +   ψ ( ρ , θ , z 0 ) θ e θ
The TEFD in the radial and angular directions is of particular importance for applications such as optical trapping and manipulation. The radial energy flow, provided by the autofocusing properties, directs the trapped particles to the optical axis [46], and the angular energy flow associated with the vortex structure of the phase leads to the rotation of the particles [47].
The field distribution in the autofocusing plane z   =   z f o c near the optical axis is of particular interest since most of the energy is concentrated there. In this case, the angular component TEFD corresponding to the phase derivative concerning the angle in Equation (10) becomes significant. Let us consider it separately. In general, this expression is rather cumbersome:
ψ ( ρ , θ , z 0 ) θ   =   1 | E ( ρ , θ , z 0 ) | 2   × ×   { [ l = 1 L D l ( ρ ) m l cos [ m l θ   +   ψ l ( ρ ) ] ] [ l = 1 L D l ( ρ ) cos [ m l θ   +   ψ l ( ρ ) ] ]   + +   [ l = 1 L D l ( ρ ) m l sin [ m l θ   +   ψ l ( ρ ) ] ] [ l = 1 L D l ( ρ ) sin [ m l θ   +   ψ l ( ρ ) ] ] }
Special cases are discussed in detail in the next section.

3. Results of Modeling and Experiment

We consider the circular Airy functions as functions depending on the radius:
A ( r )   =   Ai ( r 0     r w ) circ ( r R )
where Ai(x) is the Airy function [48], circ(r/R) is the circle function with unit amplitude, radius R, r 0 is the radial displacement parameter, and w is the normalizing parameter.
Instead of function in Equation (12), one can use other circular beams with autofocusing properties, for example, circular Pearcey beams [17,18], chirped beams [22,25,49], and generalized Airy beams [24]. This mainly affects autofocusing characteristics such as focal trajectory curvature, sharpness and focusing distance.
The following parameters were used in the simulation: wavelength λ = 0.532 μm, w = 5 mm, r 0 = 2 mm, and R = 1 mm.
The experimental parameters coincide with the mentioned simulation parameters. The used experimental setup is shown in Figure 1. The initial linearly polarized laser beam (λ = 532 nm) was extended and spatially filtered by a system composed of a pinhole (PH) (aperture size of 40 µm) and lens L1 (focal length of 350 mm). The collimated laser beam was directed onto a display of a reflective SLM (HOLOEYE, PLUTO VIS with a 1920 × 1080 pixel resolution), which was used to realize the phase masks of the designed elements for the generation of autofocusing optical vortex (OV) beams. Then, the laser beam modulated by the SLM was spatially filtered with a 4f imaging optical system consisting of lenses L2 and L3 (focal lengths of 350 and 150 mm, respectively) and a circular aperture D. This system allowed blocking a portion of laser radiation that was not modulated by SLM due to its pixilated structure. A video camera (CAM) mounted on the optical rail was used to record the generated intensity distributions at different distances from the plane z = 0. The longitudinal intensity distributions were reconstructed from arrays of cross-sections of the recorded transverse intensity distributions along the vertical axis.
Figure 2 shows the results of propagation modeling for the field defined by Equation (1) with the radial function defined by Equation (13), which does not have any angular dependence. The longitudinal distribution pattern shows that autofocusing occurs at the distance z f o c   =   200   mm (marked with the dashed line). Figure 2 also shows the transverse patterns of the field amplitude at different distances, as well as the corresponding distributions of the TEFD in these planes.
In the considered case, the TEFD only had the radial component. At distances z   <   z f o c , the flow was directed to the optical axis (Figure 2d,e), and after focusing, i.e., at distances z   >   z f o c , the flow was directed away from the optical axis (Figure 2f).

3.1. Circular Airy Beams with Vortex Superposition

In this section, we consider the field in Equation (1) with the radial function A ( r ) as the circular Airy function (12) and the angular function as a vortex superposition Φ ( φ ) = l = 1 L exp ( i m l φ ) . A particular case of a beam with a single vortex phase is Φ ( φ ) = exp ( i m 1 φ ) . In this case, the angular component of the TEFD defined by Equation (11) has the following form:
F θ ( ρ , θ , z 0 ) | E ( ρ , θ , z 0 ) | 2 ψ ( ρ , θ , z 0 ) θ   =   | C 1 ( ρ , z 0 ) | 2 m 1
Thus, the speed of rotation of the trapped particle will be proportional to the order m 1 of the vortex singularity present in the beam.
Figure 3 shows the results of propagation modeling for the field defined by Equation (1) in the presence of a single vortex phase of the order m 1 . In this case, the TEFD, in addition to the radial component, had an angular component (see Equation (13)). Note that in the autofocusing plane z f o c   =   200   mm , the radial component was close to zero, so the main energy flow was directed along the ring (Figure 3e). There were both radial and angular components before and after the focal plane, so the energy flow was spiraled (Figure 3f).
Next, we considered the field containing two vortex terms:
f ( r , φ )   =   A ( r ) [ exp ( i m 1 φ )   +   exp ( i m 2 φ ) ]
In this case, the field intensity in a certain plane z   =   z 0 was:
| E ( ρ , θ , z 0 ) | 2   =   ( k z 0 ) 2 { | C 1 ( ρ , z 0 ) | 2   +   | C 2 ( ρ , z 0 ) | 2   + +   2 | C 1 ( ρ , z 0 ) | | C 2 ( ρ , z 0 ) | cos [ ( m 1     m 2 ) θ   +   arg [ C 1 ( ρ , z 0 ) ]     arg [ C 1 ( ρ , z 0 ) ] ] }
As follows from Equation (15), the angular structure of the field depend only on the difference in the orders of optical vortices ( m 1     m 2 ) , whereas variations in the radial functions C l ( ρ , z 0 ) lead only to a scale change and rotation of this structure. A similar result was obtained in [44,45].
The phase is calculated based on the expression:
tan [ ψ ( ρ , θ , z 0 ) ]   =   D 1 ( ρ ) sin [ m 1 θ   +   ψ 1 ( ρ ) ]   +   D 2 ( ρ ) sin [ m 2 θ   +   ψ 2 ( ρ ) ] D 1 ( ρ ) cos [ m 1 θ   +   ψ 1 ( ρ ) ]   +   D 2 ( ρ ) cos [ m 2 θ   +   ψ 2 ( ρ ) ]
Since D 1 ( ρ )     D 2 ( ρ ) , the form of the formula for the phase gradient defined by Equation (11) and intensity defined by Equation (9) was not simplified compared with the general case L > 2.
Table 1 shows the results of simulations and experiments for the field defined by Equation (14) with different orders of vortex terms m 1 ,   m 2 .
The asymmetry concerning sign change m 2 is easily explained based on Equation (15). When using a superposition of two vortex fields, the symmetry of the generated field is determined by the difference | m 1     m 2 | . In the examples considered, m 1   =   2 ,   m 2   =   1 , the input field amplitude (first column of Table 1) should be equal to zero at φ   =   180 ° , and for m 1   =   2 ,   m 2   =   1 , at φ   =   60 ° ,   180 ° ,   300 ° . The transverse structure was preserved during propagation with an accuracy of scale and rotation, which is fully confirmed by both the simulation and experiment.
The corresponding distributions of the TEFD show that particles trapped in the region of maximum intensity will rotate spirally as a whole structure, decreasing on a scale when moving to the focal plane and expanding after it. The number of maxima and the symmetry of the transverse pattern are entirely determined by the function Φ ( φ ) , which can be dynamically varied using SLM.
The use of more than three vortex terms in Equation (1) will complicate the formed patterns and increase the number of degrees of freedom in variations of the transverse field pattern [50]; however, two terms are sufficient to form a wide variety of structures [51].
The proposed approach based on the representation of the angular part of the field defined by Equation (1) in the form of a superposition of optical vortices provides a simple and convenient method to control the 3D structure of autofocusing beams.

3.2. Azimuthally Modulated Circular Vortex Airy Beams

In [14], azimuthally modulated circular Airy beams were investigated, which can be considered a special case of the field defined by Equation (14), with m 2   =   m 1   =   q (q is an arbitrary integer):
f ( r , φ )   =   A ( r ) cos ( q φ )
In this case, d ψ / d θ   =   0 (except for the lines of phase jumps), since the phase is piecewise constant (has values 0 and π).
Figure 4 shows the simulation results for the field defined by Equation (17) with azimuthal modulation of the order of q = 2. The TEFD only has a radial component: before focusing ( z   <   z f o c ), the energy flow is directed toward the optical axis (Figure 4d), and after focusing ( z   >   z f o c ), the energy flow is directed from the optical axis (Figure 4f). The difference compared with classical circular autofocusing beams is that the bright region is separated by zero-intensity stripes in angular sectors, and the motion of trapped particles is limited not only in the radial but also in the angular direction. In the autofocusing plane ( z   =   z f o c ) the energy flow has the opposite direction on the adjacent rings (Figure 4e), which corresponds to its absence in the transverse direction. Note that these rings were not observed in the field amplitude picture, since they arose due to the phase gradient.
Let us consider the beam defined by Equation (17) with an additional vortex phase singularity:
f ( r , φ )   =   A ( r ) cos ( q φ ) exp ( i m φ )
The field defined by Equation (18) can be represented as the sum of two vortex terms:
f ( r , φ )   =   A ( r ) [ exp ( i p 1 φ )   +   exp ( i p 2 φ ) ]
where p 1   =   q   +   m , p 2   =   q     m .
Equation (19) coincides with Equation (14) up to notation. However, the situation considered in this section differs from that considered in the previous section. We initially had a beam structure divided into angular sectors (azimuthally modulated autofocusing beam) and varied it due to the additionally inserted vortex phase.
Equation (19) shows that for m = ±q, one of the terms becomes a constant, which provide a nonzero value of the intensity on the optical axis in the autofocusing region. This is another difference from the situation discussed in the previous section.
Table 2 shows the results of the simulation and experiment for the field defined by Equation (18) with azimuthal modulation of the order of q = 2 and with different orders of additional vortex phase singularity m.
It follows from Equation (18) that at q = 2, regardless of m, the amplitude of the input field should be equal to zero φ   =   45 ° ,   135 ° ,   225 ° ,   315 ° . The m value only affects the phase structure. Both dependencies are fully confirmed by the pictures in the first column.
We argue that the quantity q determines the order of symmetry of the transverse intensity distribution, which is equal to 2q, while the detailed picture of the distribution, including the value of the field at the optical axis in the focal region, depends on the quantity m. When m = ±q, the central spot was at the optical axis.
Detailed results of the TEFD investigation in the focal region for the considered beams are shown in Table 3. The effects of the additional vortex phase singularity are most significant in the focal plane but are also important in the region after autofocusing. This can be explained by the interference interaction of different parts of the beam, which is not present before the focal region. In the areas before and after the focus, the spiral nature of the energy flow is visible. In the plane of focusing, one can observe a more complex structure associated with the influence of an additional vortex, including the presence of several regions with different directions of the angular motion of the energy flow.
As seen from the results in Table 2 and Table 3, inserting an additional vortex phase singularity into the beam defined by Equation (17) significantly changes the field distribution in the autofocusing region.
Comparison of the results of Section 3.1 and Section 3.2 (Table 1 and Table 3) shows that using the proposed approach, different scenarios for controlling the transverse intensity structure and TEFD of autofocusing beams are possible. In the first case, a set of bright spots can be formed in a cross-section, which rotates as a single structure when the longitudinal distance changes. In the second case, a certain symmetry of the transverse structure is fixed, and the inserting of an additional vortex phase allows it to be changed locally, ensuring the presence or absence of a central light spot. The latter property is useful in the formation of optical bottles based on autofocusing beams.

4. Conclusions

In this paper, we considered the possibility of controlling the transverse intensity structure and TEFD of autofocusing beams by inserting a combination of optical vortices into the beam. This approach is convenient for experimental implementation since a separate optical element can be used to form an autofocusing beam, while variations in the transverse distribution can be performed dynamically using SLM.
Based on the proposed approach, we showed, numerically and experimentally, the possibility of forming autofocusing beams, the transverse intensity pattern of which rotates as a whole structure during beam propagation. The corresponding distributions of the TEFD show that particles trapped in the light spots rotate spirally as a single structure, decreasing on a scale when moving to the focal plane and expanding after it. The structure and symmetry of the transverse pattern are determined by the angular function, which is a superposition of optical vortices.
Another option is considered when an azimuthally modulated autofocusing beam with an additional vortex phase is used. In this case, a certain symmetry of the transverse structure is fixed, and the inserting of the vortex phase allows it to be changed locally, ensuring the presence or absence of a central light spot, which can be useful in the formation of optical bottles. The effect of the additional vortex phase is most significant near the focusing region. Complex distributions of the TEFD are formed, including those with several areas with different directions of angular energy flow.
This variety of configurations is a useful property in the implementation of the trapping/confinement of particles, especially when using dynamic optical elements, in particular SLM. For example, a sufficiently large annular region is convenient for initial trapping, after which the element switches to creating sector regions that narrow the particle localization area. The proposed approach provides a simple and convenient way to control the 3D structure of autofocusing beams.

Author Contributions

Conceptualization, S.N.K.; formal analysis, A.V.U. and S.N.K.; investigation, S.N.K. and A.P.P.; resources, M.A.B.; writing—original draft preparation, S.N.K. and A.P.P.; writing—review and editing, M.A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Foundation for Basic Research, grant numbers 20-07-00505, 18-58-14001.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work was partly funded by the Russian Foundation for Basic Research under grants No. 20-07-00505 (simulation results) and 18-58-14001 (experimental results), the Ministry of Science and Higher Education within the government project of FSRC Crystallography and Photonics RAS, under agreement 007-GZ/Ch3363/26 (theoretical part).

Conflicts of Interest

The authors declare that there are no conflict of interest related to this article.

References

  1. Grier, D.G. A revolution in optical manipulation. Nature 2003, 424, 810–816. [Google Scholar] [CrossRef]
  2. Zhan, Q. Cylindrical vector beams: From mathematical concepts to applications. Adv. Opt. Photonics 2009, 1, 1–57. [Google Scholar] [CrossRef]
  3. Fahrbach, F.O.; Simon, P.; Rohrbach, A. Microscopy with self-reconstructing beams. Nat. Photonics 2010, 4, 780–785. [Google Scholar] [CrossRef]
  4. Padgett, M.; Bowman, R. Tweezers with a twist. Nat. Photonics 2011, 5, 343–348. [Google Scholar] [CrossRef]
  5. Andrews, D.L. Structured Light and Its Applications: An Introduction to Phase-Structured Beams and Nanoscale Optical Forces; Academic Press: San Diego, CA, USA, 2011; p. 400. [Google Scholar]
  6. Wördemann, M. Structured Light Fields: Applications in Optical Trapping, Manipulation, and Organisation; Springer Science and Business Media: Berlin/Heidelberg, Germany, 2012; p. 136. [Google Scholar]
  7. Litchinitser, N.M. Structured light meets structured matter. Science 2012, 337, 1054–1055. [Google Scholar] [CrossRef]
  8. Rubinsztein-Dunlop, H.; Forbes, A.; Berry, M.; Dennis, M.; Andrews, D.L.; Mansuripur, M.; Denz, C.; Alpmann, C.; Banzer, P.; Bauer, T.; et al. Roadmap on structured light. J. Opt. 2017, 19, 013001. [Google Scholar] [CrossRef]
  9. Khonina, S.N.; Kazanskiy, N.L.; Karpeev, S.V.; Butt, M.A. Bessel beam: Significance and applications—A progressive review. Micromachines 2020, 11, 997. [Google Scholar] [CrossRef]
  10. Kotlyar, V.V.; Kovalev, A.A.; Nalimov, A.G.; Porfirev, A.P. Evolution of an optical vortex with an initial fractional topological charge. Phys. Rev. A 2020, 102, 023516. [Google Scholar] [CrossRef]
  11. Efremidis, N.K.; Christodoulides, D.N. Abruptly autofocusing waves. Opt. Lett. 2010, 35, 4045–4047. [Google Scholar] [CrossRef]
  12. Papazoglou, D.G.; Efremidis, N.K.; Christodoulides, D.N.; Tzortzakis, S. Observation of abruptly autofocusing waves. Opt. Lett. 2011, 36, 1842–1844. [Google Scholar] [CrossRef]
  13. Davis, J.A.; Cottrell, D.M.; Sand, D. Abruptly autofocusing vortex beam. Opt. Express 2012, 20, 13302–13310. [Google Scholar] [CrossRef]
  14. Porfirev, A.P.; Khonina, S.N. Generation of the azimuthally modulated circular superlinear Airy beams. J. Opt. Soc. Am. B 2017, 34, 2544–2549. [Google Scholar] [CrossRef]
  15. Jiang, Y.; Zhao, S.; Yu, W.; Zhu, X. Abruptly autofocusing property of circular Airy vortex beams with different initial launch angles. J. Opt. Soc. Am. A 2018, 35, 890–894. [Google Scholar] [CrossRef]
  16. Ring, J.; Lindberg, J.; Mourka, A.; Mazilu, M.; Dholakia, K.; Dennis, M. Auto-focusing and self-healing of Pearcey beams. Opt. Express 2012, 20, 18955–18966. [Google Scholar] [CrossRef] [PubMed]
  17. Chen, X.; Deng, D.; Zhuang, J.; Yang, X.; Liu, H.; Wang, G. Nonparaxial propagation of abruptly autofocusing circular Pearcey Gaussian beams. Appl. Opt. 2018, 57, 8418–8423. [Google Scholar] [CrossRef]
  18. Sun, C.; Deng, D.; Wang, G.; Yang, X.; Hong, W. Abruptly autofocusing properties of radially polarized circle Pearcey vortex beams. Opt. Commun. 2020, 457, 124690. [Google Scholar] [CrossRef]
  19. Khonina, S.N.; Ustinov, A.V.; Porfirev, A.P. Aberration laser beams with autofocusing properties. Appl. Opt. 2018, 57, 1410–1416. [Google Scholar] [CrossRef] [PubMed]
  20. Khonina, S.N. Specular and vortical Airy beams. Opt. Commun. 2011, 284, 4263–4271. [Google Scholar] [CrossRef]
  21. Vaveliuk, P.; Lencina, A.; Rodrigo, J.A.; Matos, O.M. Symmetric Airy beams. Opt. Lett. 2014, 39, 2370–2373. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Chremmos, I.; Efremidis, N.K.; Christodoulides, D.N. Pre-engineered abruptly autofocusing beams. Opt. Lett. 2011, 36, 1890–1892. [Google Scholar] [CrossRef]
  23. Li, P.; Liu, S.; Peng, T.; Xie, G.; Gan, X.; Zhao, J. Spiral autofocusing Airy beams carrying power-exponent phase vortices. Opt. Express 2014, 22, 7598–7606. [Google Scholar] [CrossRef] [PubMed]
  24. Khonina, S.N.; Ustinov, A.V. Fractional Airy beams. J. Opt. Soc. Am. A 2017, 34, 1991–1999. [Google Scholar] [CrossRef] [PubMed]
  25. Khonina, S.N.; Porfirev, A.P.; Ustinov, A.V. Sudden autofocusing of superlinear chirp beams. J. Opt. 2018, 20, 025605. [Google Scholar] [CrossRef]
  26. Ddholakia, K.; Čižmár, T. Shaping the future of manipulation. Nat. Photonics 2011, 5, 335–342. [Google Scholar] [CrossRef]
  27. Jiang, Y.F.; Huang, K.K.; Lu, X.H. Radiation force of abruptly autofocusing Airy beams on a Rayleigh particle. Opt. Express 2013, 21, 24413–24421. [Google Scholar] [CrossRef]
  28. Lu, W.; Sun, X.; Chen, H.; Liu, S.; Lin, Z. Abruptly autofocusing property and optical manipulation of circular Airy beams. Phys. Rev. A 2019, 99, 013817. [Google Scholar] [CrossRef]
  29. Panagiotopoulos, P.; Papazoglou, D.G.; Couairon, A.; Tzortzakis, S. Abruptly autofocusing beams enable advanced multiscale photo-polymerization. Nat. Commun. 2013, 4, 2622. [Google Scholar] [CrossRef] [Green Version]
  30. Panagiotopoulos, P.; Couairon, A.; Kolesik, M.; Papazoglou, D.G.; Moloney, J.V.; Tzortzakis, S. Nonlinear plasma-assisted collapse of ring-Airy wave packets. Phys. Rev. A 2016, 93, 033808. [Google Scholar] [CrossRef] [Green Version]
  31. Liu, D.; Wang, Y.; Zhai, Z.; Fang, Z.; Tao, Q.; Perrie, W.; Edwarson, S.P.; Dearden, G. Dynamic laser beam shaping for material processing using hybrid holograms. Opt. Laser Technol. 2018, 102, 68–73. [Google Scholar] [CrossRef]
  32. Ustinov, A.V.; Khonina, S.N.; Khorin, P.A.; Porfirev, A.P. Control of the intensity distribution along the light spiral generated by a generalized spiral phase plate. J. Opt. Soc. Am. B 2021, 38, 420–427. [Google Scholar] [CrossRef]
  33. Zhang, Y.; Li, P.; Liu, S.; Han, L.; Cheng, H.; Zhao, J. Manipulating spin-dependent splitting of vector abruptly autofocusing beam by encoding cosine-azimuthal variant phases. Opt. Express 2016, 24, 28409–28418. [Google Scholar] [CrossRef]
  34. Yan, X.; Guo, L.; Cheng, M.; Chai, S. Free-space propagation of autofocusing Airy vortex beams with controllable intensity gradients. Chin. Opt. Lett. 2019, 17, 040101. [Google Scholar] [CrossRef]
  35. Brimis, A.; Makris, K.G.; Papazoglou, D.G. Tornado waves. Opt. Lett. 2020, 45, 280–283. [Google Scholar] [CrossRef]
  36. Wu, Y.; Xu, C.; Lin, Z.; Qiu, H.; Fu, X.; Chen, K.; Deng, D. Abruptly autofocusing polycyclic tornado ring Airy beam. New J. Phys. 2020, 22, 093045. [Google Scholar] [CrossRef]
  37. Zhu, J.; Chen, Y.; Zhang, Y.; Cai, X.; Yu, S. Spin and orbital angular momentum and their conversion in cylindrical vector vortices. Opt. Lett. 2014, 39, 4435–4438. [Google Scholar] [CrossRef] [PubMed]
  38. Bliokh, K.Y.; Rodriguez-Fortuno, F.; Nori, F.; Zayats, A.V. Spin-orbit interactions of light. Nat. Photonics 2015, 9, 796–808. [Google Scholar] [CrossRef]
  39. Porfirev, A.P.; Ustinov, A.V.; Khonina, S.N. Polarization conversion when focusing cylindrically polarized vortex beams. Sci. Rep. 2016, 6, 1–9. [Google Scholar] [CrossRef] [Green Version]
  40. Shi, P.; Du, L.; Yuan, X. Structured spin angular momentum in highly focused cylindrical vector vortex beams for optical manipulation. Opt. Express 2018, 26, 23449–23459. [Google Scholar] [CrossRef]
  41. Bekshaev, A.Y.; Soskin, M.S. Transverse energy flows in vectorial fields of paraxial beams with singularities. Opt. Commun. 2007, 271, 332–348. [Google Scholar] [CrossRef]
  42. Bekshaev, A.Y. Internal energy flows and instantaneous field of a monochromatic paraxial light beam. Appl. Opt. 2012, 51, C13–C16. [Google Scholar] [CrossRef] [PubMed]
  43. Khonina, S.N.; Ustinov, A.V.; Logachev, V.I.; Porfirev, A.P. Properties of vortex light fields generated by generalized spiral phase plates. Phys. Rev. A 2020, 101, 043829. [Google Scholar] [CrossRef]
  44. Kotlyar, V.V.; Soifer, V.A.; Khonina, S.N. An algorithm for the generation of laser beams with longitudinal periodicity: Rotating images. J. Mod. Opt. 1997, 44, 1409–1416. [Google Scholar] [CrossRef]
  45. Khonina, S.N.; Kotlyar, V.V.; Soifer, V.A.; Lautanen, J.; Honkanen, M.; Turunen, J. Generating a couple of rotating nondiffracting beams using a binary-phase DOE. Optik 1999, 110, 137–144. [Google Scholar]
  46. Zhang, P.; Prakash, J.; Zhang, Z.; Mills, M.S.; Efremidis, N.K.; Christodoulides, D.N.; Chen, Z.G. Trapping and guiding microparticles with morphing autofocusing Airy beams. Opt. Lett. 2011, 36, 2883–2885. [Google Scholar] [CrossRef] [PubMed]
  47. Chen, M.; Huang, S.; Liu, X.; Chen, Y.; Shao, W. Optical trapping and rotating of micro-particles using the circular Airy vortex beams. Appl. Phys. B 2019, 125, 184. [Google Scholar] [CrossRef]
  48. Vallee, O. Airy Functions and Applications in Physics; Vallee, O., Soares., M., Eds.; Imperial College Press: London, UK, 2004; p. 194. ISBN 978-1-86094-478-9. [Google Scholar]
  49. Degtyarev, S.A.; Volotovsky, S.G.; Khonina, S.N. Sublinearly chirped metalenses for forming abruptly autofocusing cylindrically polarized beams. J. Opt. Soc. Am. B 2018, 35, 1963–1969. [Google Scholar] [CrossRef]
  50. Almazov, A.A.; Khonina, S.N.; Kotlyar, V.V. Using phase diffraction optical elements to shape and select laser beams consisting of a superposition of an arbitrary number of angular harmonics. J. Opt. Technol. 2005, 72, 391–399. [Google Scholar] [CrossRef]
  51. Khonina, S.N.; Kotlyar, V.V.; Soifer, V.A.; Jefimovs, K.; Turunen, J. Generation and selection of laser beams represented by a superposition of two angular harmonics. J. Mod. Opt. 2004, 51, 761–773. [Google Scholar] [CrossRef]
Figure 1. Experimental setup for the investigation of the generation and propagation of designed autofocusing optical vortex (OV) beams: PH, pinhole (aperture size of 40 µm); L1, L2, and L3 are lenses (f1 = 350, f2 = 350, and f3 = 150 mm, respectively); SLM, spatial light modulator (HOLOEYE, PLUTO VIS with a 1920 × 1080 pixel resolution); D, circular aperture; and CAM, video camera.
Figure 1. Experimental setup for the investigation of the generation and propagation of designed autofocusing optical vortex (OV) beams: PH, pinhole (aperture size of 40 µm); L1, L2, and L3 are lenses (f1 = 350, f2 = 350, and f3 = 150 mm, respectively); SLM, spatial light modulator (HOLOEYE, PLUTO VIS with a 1920 × 1080 pixel resolution); D, circular aperture; and CAM, video camera.
Micromachines 12 00297 g001
Figure 2. Simulation results for the field defined by Equation (1), which does not have any angular dependence: (a) amplitude and (b) phase of the input field ( x , y     [ 1   mm ,   1   mm ] ); (c) longitudinal distribution of the amplitude ( z     [ 20   mm ,   500   mm ] , y     [ 1   mm ,   1   mm ] ), and transverse patterns ( x , y     [ 0.25   mm ,   0.25   mm ] ) of the field amplitude at different distances (top line): (d) z = 150 mm, (e) z f o c   =   200   mm , and (f) z = 250 mm, as well as the corresponding distributions of the transverse energy flow density (TEFD) amplitude in these planes (bottom line). The red color corresponds to the radial component | F ρ ( ρ , θ ) | , the blue color corresponds to the angular component | F θ ( ρ , θ ) | , and the arrows show the direction of flow.
Figure 2. Simulation results for the field defined by Equation (1), which does not have any angular dependence: (a) amplitude and (b) phase of the input field ( x , y     [ 1   mm ,   1   mm ] ); (c) longitudinal distribution of the amplitude ( z     [ 20   mm ,   500   mm ] , y     [ 1   mm ,   1   mm ] ), and transverse patterns ( x , y     [ 0.25   mm ,   0.25   mm ] ) of the field amplitude at different distances (top line): (d) z = 150 mm, (e) z f o c   =   200   mm , and (f) z = 250 mm, as well as the corresponding distributions of the transverse energy flow density (TEFD) amplitude in these planes (bottom line). The red color corresponds to the radial component | F ρ ( ρ , θ ) | , the blue color corresponds to the angular component | F θ ( ρ , θ ) | , and the arrows show the direction of flow.
Micromachines 12 00297 g002
Figure 3. Simulation results for the field defined by Equation (1) in the presence of a single vortex phase of order m 1   =   1 (the rest of the description is as in Figure 2).
Figure 3. Simulation results for the field defined by Equation (1) in the presence of a single vortex phase of order m 1   =   1 (the rest of the description is as in Figure 2).
Micromachines 12 00297 g003
Figure 4. Simulation results for the field defined by Equation (17) with azimuthal modulation of the order of q = 2 (the rest of the description is as in Figure 2).
Figure 4. Simulation results for the field defined by Equation (17) with azimuthal modulation of the order of q = 2 (the rest of the description is as in Figure 2).
Micromachines 12 00297 g004
Table 1. Results of simulations and experiments for the field defined by Equation (14) with different orders of vortex terms m1, m2.
Table 1. Results of simulations and experiments for the field defined by Equation (14) with different orders of vortex terms m1, m2.
Input Amplitude and Phase
(2 mm × 2 mm)
Longitudinal Intensity Distribution
(2 mm × 300 mm)
Transverse Distribution
z = 150 mmz = 200 mmz = 250 mm
m1 = 2, m2 = 1
Micromachines 12 00297 i001
Micromachines 12 00297 i002
Simulation
Micromachines 12 00297 i003
Experiment
Micromachines 12 00297 i004
Simulation: intensity (1.6 mm × 1.6 mm)
Micromachines 12 00297 i005
Simulation: TEFD (1 mm × 1 mm)
Micromachines 12 00297 i006
Experiment: intensity (1.6 mm × 1.6 mm)
Micromachines 12 00297 i007
m1 = 2, m2 = −1
Micromachines 12 00297 i008
Micromachines 12 00297 i009
Simulation
Micromachines 12 00297 i010
Experiment
Micromachines 12 00297 i011
Simulation: intensity (1.6 mm × 1.6 mm)
Micromachines 12 00297 i012
Simulation: TEFD (1 mm × 1 mm)
Micromachines 12 00297 i013
Experiment: intensity (1.6 mm × 1.6 mm)
Micromachines 12 00297 i014
Table 2. Results of modeling and experiment for the field defined by Equation (18) with azimuthal modulation of order q = 2 with different orders of additional vortex phase singularity m.
Table 2. Results of modeling and experiment for the field defined by Equation (18) with azimuthal modulation of order q = 2 with different orders of additional vortex phase singularity m.
Input Amplitude and Phase
(2 mm × 2 mm)
Longitudinal Intensity Distribution
(1.6 mm × 300 mm)
Transverse Intensity Distribution
(1.6 mm × 1.6 mm)
z = 100 mmz = 200 mmz = 300 mm
q = 2, m = 1 SimulationSimulation
Micromachines 12 00297 i015 Micromachines 12 00297 i016 Micromachines 12 00297 i017
ExperimentExperiment
Micromachines 12 00297 i018 Micromachines 12 00297 i019
q = 2, m = 2 SimulationSimulation
Micromachines 12 00297 i020 Micromachines 12 00297 i021 Micromachines 12 00297 i022
ExperimentExperiment
Micromachines 12 00297 i023 Micromachines 12 00297 i024
Table 3. Results of TEFD simulation for the field defined by Equation (18) with azimuthal modulation of order q = 2 and with different orders of additional vortex phase singularity m.
Table 3. Results of TEFD simulation for the field defined by Equation (18) with azimuthal modulation of order q = 2 and with different orders of additional vortex phase singularity m.
Input Amplitude and Phase
(2 mm × 2 mm)
Transverse Distribution (1 mm × 1 mm)
Value z = 150 mm z = 175 mm z = 200 mm z = 225 mm z = 250 mm
q = 2, m = 1
Micromachines 12 00297 i025
| E ( ρ , θ ) | Micromachines 12 00297 i026
| F ( ρ , θ ) | Micromachines 12 00297 i027
q = 2, m = 2
Micromachines 12 00297 i028
| E ( ρ , θ ) | Micromachines 12 00297 i029
| F ( ρ , θ ) | Micromachines 12 00297 i030
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khonina, S.N.; Porfirev, A.P.; Ustinov, A.V.; Butt, M.A. Generation of Complex Transverse Energy Flow Distributions with Autofocusing Optical Vortex Beams. Micromachines 2021, 12, 297. https://doi.org/10.3390/mi12030297

AMA Style

Khonina SN, Porfirev AP, Ustinov AV, Butt MA. Generation of Complex Transverse Energy Flow Distributions with Autofocusing Optical Vortex Beams. Micromachines. 2021; 12(3):297. https://doi.org/10.3390/mi12030297

Chicago/Turabian Style

Khonina, Svetlana N., Alexey P. Porfirev, Andrey V. Ustinov, and Muhammad Ali Butt. 2021. "Generation of Complex Transverse Energy Flow Distributions with Autofocusing Optical Vortex Beams" Micromachines 12, no. 3: 297. https://doi.org/10.3390/mi12030297

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop