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

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.

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 Micromachines 2021, 12, 297 2 of 12 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.

Theoretical Analysis
Let us consider vortex fields with autofocusing properties in the following form: 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(im 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: 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: The full field will be: Let's rewrite Equation (4) as follows: 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 : The TEFD in the plane z = z 0 is proportional to the product of the field intensity and phase gradient [32,42,43]: 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: 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: 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 Micromachines 2021, 12, 297 4 of 12 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 oc 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: Special cases are discussed in detail in the next section.

Results of Modeling and Experiment
We consider the circular Airy functions as functions depending on the radius: 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 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.      (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 oc = 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.
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 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 200 mm foc z = (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 foc z z < , the flow was directed to the optical axis (Figure 2d,e), and after focusing, i.e., at distances foc z z > , the flow was directed away from the optical axis (Figure 2f).

Circular Airy Beams with Vortex Superposition
In this section, we consider the field in Equation (1)   In the considered case, the TEFD only had the radial component. At distances z < z f oc , the flow was directed to the optical axis ( Figure 2d,e), and after focusing, i.e., at distances z > z f oc , the flow was directed away from the optical axis ( Figure 2f).

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 In this case, the angular component of the TEFD defined by Equation (11) has the following form: 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 oc = 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). 1 dition to the radial component, had an angular component (see Equation (13)). Note that in the autofocusing plane 200 mm foc z = , 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).  (1) in the presence of a single vortex phase of order 1 1 m = (the rest of the description is as in Figure 2).
Next, we considered the field containing two vortex terms: In this case, the field intensity in a certain plane 0 z z = was: As follows from Equation (15) [44,45].
The phase is calculated based on the expression: ρ , 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.  (1) in the presence of a single vortex phase of order m 1 = 1 (the rest of the description is as in Figure 2).
Next, we considered the field containing two vortex terms: In this case, the field intensity in a certain plane z = z 0 was: 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: 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 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 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 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 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 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 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.  (14)  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. 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. 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. 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 corresponding distributions of the TEFD show that particles trapped in the re-gion 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. gion 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.

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 2 1 m m q = − = (q is an arbitrary integer): In this case, 0 d d ψ θ = (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 azi- 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.

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ϕ) (17) 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 oc ), the energy flow is directed toward the optical axis (Figure 4d), and after focusing (z > z f oc ), 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 oc ) 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.
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.

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 2 1 m m q = − = (q is an arbitrary integer): In this case, 0 d d ψ θ = (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)  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.  (17) with azimuthal modulation of the order of q = 2 (the rest of the description is as in Figure 2).
Let us consider the beam defined by Equation (17) with an additional vortex phase singularity:  (17) with azimuthal modulation of the order of q = 2 (the rest of the description is as in Figure 2). Let us consider the beam defined by Equation (17) with an additional vortex phase singularity: The field defined by Equation (18) can be represented as the sum of two vortex terms: 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. Table 2. Results of modeling and experiment for the field defined by Equation (18)  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 Tables 2 and 3, inserting an additional vortex phase singularity into the beam defined by Equation (17) significantly changes the field distribution in the autofocusing region. Table 2. Results of modeling and experiment for the field defined by Equation (18)  Experiment Experiment 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 Tables 2 and 3, inserting an additional vortex phase singularity into the beam defined by Equation (17) significantly changes the field distribution in the autofocusing region. Table 2. Results of modeling and experiment for the field defined by Equation (18)  Experiment Experiment 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 Tables 2 and 3, inserting an additional vortex phase singularity into the beam defined by Equation (17) significantly changes the field distribution in the autofocusing region. Table 2. Results of modeling and experiment for the field defined by Equation (18)

Experiment Experiment
Experiment Experiment 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 Tables 2 and 3, inserting an additional vortex phase singularity into the beam defined by Equation (17) significantly changes the field distribution in the autofocusing region. Table 2. Results of modeling and experiment for the field defined by Equation (18)  Experiment Experiment 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 Tables 2 and 3, inserting an additional vortex phase singularity into the beam defined by Equation (17) significantly changes the field distribution in the autofocusing region.  Table 3. Results of TEFD simulation for the field defined by Equation (18) Table 3. Results of TEFD simulation for the field defined by Equation (18) Table 3. Results of TEFD simulation for the field defined by Equation (18) Table 3. Results of TEFD simulation for the field defined by Equation (18) Table 3. Results of TEFD simulation for the field defined by Equation (18)  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 Tables 2 and 3, inserting an additional vortex phase singularity into the beam defined by Equation (17) significantly changes the field distribution in the autofocusing region. Table 3. Results of TEFD simulation for the field defined by Equation (18) Table 3. Results of TEFD simulation for the field defined by Equation (18) Comparison of the results of Sections 2.1 and 2.2 (Tables 1 and 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.

Conclusions
|E(ρ, θ)| Table 3. Results of TEFD simulation for the field defined by Equation (18) Comparison of the results of Sections 2.1 and 2.2 (Tables 1 and 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.

Conclusions
|F(ρ, θ)| Table 3. Results of TEFD simulation for the field defined by Equation (18) Comparison of the results of Sections 2.1 and 2.2 (Tables 1 and 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. Table 3. Results of TEFD simulation for the field defined by Equation (18)

Conclusions
Comparison of the results of Sections 2.1 and 2.2 (Tables 1 and 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. Table 3. Results of TEFD simulation for the field defined by Equation (18)

|E(ρ, θ)|
Comparison of the results of Sections 2.1 and 2.2 (Tables 1 and 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.

Conclusions
|F(ρ, θ)| Table 3. Results of TEFD simulation for the field defined by Equation (18) (Tables 1 and 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.

Conclusions
Comparison of the results of Sections 3.1 and 3.2 (Tables 1 and 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.

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.