Optical forces induced behavior of a particle in a non-diffracting vortex beam

An interaction between a light field with complex field spatial distribution and a micro-particle leads to forces that drag the particle in space and may confine it in a stable position or a trajectory. The particle behavior is determined by its size with respect to the characteristic length of the spatially periodic or symmetric light field distribution. We study theoretically and experimentally the behavior of a microparticle near the center of an optical vortex beam in a plane perpendicular to the beam propagation. We show that such particle may be stably trapped either in a dark spot on the vortex beam axis, or in one of two points placed off the optical axis. It may also circulate along a trajectory having its radius smaller or equal to the radius of the first bright vortex ring. © 2012 Optical Society of America OCIS codes: (140.3300) Laser beam shaping; (170.4520) Optical confinement and manipulation; (350.4855) Optical tweezers or optical manipulation. References and links 1. A. Ashkin, J. M. Dziedzic, J. E. Bjorkholm, and S. Chu, “Observation of a single-beam gradient force optical trap for dielectric particles,” Opt. Lett. 11, 288–290 (1986). 2. J. E. Curtis, B. A. Koss, and D. G. Grier, “Dynamic holographic optical tweezers,” Opt. Commun. 207, 169–175 (2002). 3. D. Grier and Y. Roichman, “Holographic optical trapping,” Appl. Opt. 45, 880–887 (2006). 4. D. V. Petrov, “Raman spectroscopy of optically trapped particles,” J. Opt. A: Pure Appl. Opt. 9, S139–S156 (2007). 5. M. M. Wang, E. Tu, D. E. Raymond, J. M. Yang, H. Zhang, N. Hagen, B. Dees, E. M. Mercer, A. H. Forster, I. Kariv, P. J. Marchand, and W. F. Butler, “Microfluidic sorting of mammalian cells by optical force switching,” Nature Biotechnol. 23, 83–87 (2005). 6. J. Guck, S. Schinkinger, B. Lincoln, F. Wottawah, S. Ebert, M. Romeyke, D. Lenz, H. M. Erickson, R. Ananthakrishnan, D. Mitchell, J. Käs, S. Ulvick, and C. Bilby, “Optical deformability as an inherent cell marker for testing malignant transformation and metastatic competence,” Biophys. J. 88, 3689–3698 (2005). 7. W. J. Greenleaf, M. T. Woodside, and S. M. Block, “High-resolution, single-molecule measurements of biomolecular motion,” Annu. Rev. Biophys. Biomol. Struct. 36, 171–190 (2007). 8. J. R. Moffitt, Y. R. Chemla, S. B. Smith, and C. Bustamante, “Recent advances in optical tweezers,” Annual Review of Biochemistry 77, 205–228 (2008). 9. P. Zemánek, A. Jonáš, and M. Liška, “Simplified description of optical forces acting on a nanoparticle in the Gaussian standing wave,” J. Opt. Soc. Am. A 19, 1025–1034 (2002). 10. M. P. MacDonald, G. C. Spalding, and K. Dholakia, “Microfluidic sorting in an optical lattice,” Nature 426, 421–424 (2003). 11. T. Čižmár, M. Šiler, M. Šerý, P. Zemánek, V. Garcés-Chávez, and K. Dholakia, “Optical sorting and detection of sub-micron objects in a motional standing wave,” Phys. Rev. B 74, 035105 (2006). 12. I. Ricárdez-Vargas, P. Rodrı́guez-Montero, R. Ramos-Garcı́a, and K. Volke-Sepúlveda, “Modulated optical sieve for sorting of polydisperse microparticles,” Appl. Phys. Lett. 88, 121116 (2006). 13. P. Jákl, T. Čižmár, M. Šerý, and P. Zemánek, “Static optical sorting in a laser interference field,” Appl. Phys. Lett. 92, 161110 (2008). #174096 $15.00 USD Received 10 Aug 2012; revised 30 Sep 2012; accepted 1 Oct 2012; published 9 Oct 2012 (C) 2012 OSA 22 October 2012 / Vol. 20, No. 22 / OPTICS EXPRESS 24304 14. K. Dholakia, M. P. MacDonald, P. Zemánek, and T. Čižmár, “Cellular and colloidal separation using optical forces,” Methods Cell Biology 82, 467–495 (2007). 15. L. Allen, M. Beijersbergen, R. Spreeuw, and J. Woerdman, “Orbital angular-momentum of light and the transformation of Laguerre-Gaussian laser modes,” Phys. Rev. A 45, 8185–8189 (1992). 16. R. Beth, “Mechanical detection and measurement of the angular momentum of light,” Phys. Rev. 50, 115–125 (1936). 17. E. Santamato, B. Daino, M. Romagnoli, M. Settembre, and Y. Shen, “Collective rotation of molecules driven by the angular-momentum of light in a nematic film,” Phys. Rev. Lett. 57, 2423–2426 (1986). 18. M. Friese, T. Nieminen, N. Heckenberg, and H. Rubinsztein-Dunlop, “Optical alignment and spinning of lasertrapped microscopic particles,” Nature 394, 348–350 (1998). 19. M. Padgett and R. Bowman, “Tweezers with a twist,” Nat. Photonics 5, 343–348 (2011). 20. M. Babiker, W. Power, and L. Allen, “Light-induced torque on moving atoms,” Phys. Rev. Lett. 73, 1239–1242 (1994). 21. W. Power, L. Allen, M. Babiker, and V. Lembessis, “Atomic motion in light beams possessing orbital angular momentum,” Phys. Rev. A 52, 479–488 (1995). 22. J. Arlt, K. Dholakia, J. Soneson, and E. M. Wright, “Optical dipole traps and atomic waveguides based on Bessel light beams,” Phys. Rev. A 63, 063602 (2001). 23. M. Babiker, C. Bennett, D. Adrews, and L. Romero, “Orbital angular momentum exchange in the interaction of twisted light with molecules,” Phys. Rev. Lett 89, 143601 (2002). 24. D. S. Bradshaw and D. L. Andrews, “Interactions between spherical nanoparticles optically trapped in LaguerreGaussian modes,” Opt. Lett. 30, 3039–3041 (2005). 25. A. O. Santillán, K. Volke-Sepúlveda, and A. Flores-Pérez, “Wave fields with a periodic orbital angular momentum gradient along a single axis: a chain of vortices,” New. J. Phys. 11, 043004 (2009). 26. K. Volke-Sepúlveda and R. Jauregui, “All-optical 3D atomic loops generated with Bessel light fields,” J. Phys. B 42, 085303 (2009). 27. I. Ricardez-Vargas and K. Volke-Sepúlveda, “Experimental generation and dynamical reconfiguration of different circular optical lattices for applications in atom trapping,” J. Opt. Soc. Am. B 27, 948–955 (2010). 28. V. Garcés-Chávez, K. Volke-Sepúlveda, S. Chávez-Cerda, W. Sibbett, and K. Dholakia, “Transfer of orbital angular momentum to an optically trapped low-index particle,” Phys. Rev. A 66, 063402 (2002). 29. K. Volke-Sepúlveda, S. Chávez-Cerda, V. Garcés-Chávez, and K. Dholakia, “Three-dimensional optical forces and transfer of orbital angular momentum from multiringed light beams to spherical microparticles,” J. Opt. Soc. Am. B 21, 1749–1757 (2004). 30. A. Ohta and Y. Kawata, “Analyses of radiation force and torque on a spherical particle near a substrate illuminated by a focused Laguerre-Gaussian beam,” Opt. Commun. 274, 269–273 (2007). 31. S. H. Simpson and S. Hanna, “Rotation of absorbing spheres in Laguerre-Gaussian beams,” J. Opt. Soc. Am. 26, 173–183 (2009). 32. J. Ng, Z. Lin, and C. Chan, “Theory of optical trapping by an optical vortex beam,” Phys. Rev. Lett. 104, 103601 (2010). 33. D. Zhang, X. Yuan, S. Tjin, and S. Krishnan, “Rigorous time domain simulation of momentum transfer between light and microscopic particles in optical trapping,” Opt. Express 12, 2220–2230 (2004). 34. M. Dienerowitz, M. Mazilu, P. J. Reece, T. F. Krauss, and K. Dholakia, “Optical vortex trap for resonant confinement of metal nanoparticles,” Opt. Express 16, 4991–4999 (2008). 35. Y. Zhao, G. Milne, J. S. Edgar, G. D. M. Jeffries, D. McGloin, and D. T. Chiu, “Quantitative force mapping of an optical vortex trap,” Appl. Phys. Lett. 92, 161111 (2008). 36. H. He, M. Friese, N. Heckenberg, and H. Rubinsztein-Dunlop, “Direct observation of transfer of angularmomentum to absorptive particles from a laser-beam with a phase singularity,” Phys. Rev. Lett. 75, 826–829 (1995). 37. M. Friese, J. Enger, H. Rubinsztein-Dunlop, and N. Heckenberg, “Optical angular-momentum transfer to trapped absorbing particles,” Phys. Rev. A 54, 1593–1596 (1996). 38. N. B. Simpson, K. Dholakia, L. Allen, and M. J. Padgett, “Mechanical equivalence of spin and orbital angular momentum of light: an optical spanner,” Opt. Lett. 22, 52–54 (1997). 39. K. Volke-Sepúlveda, V. Garcés-Chávez, S. Chávez-Cerda, J. Arlt, and K. Dholakia, “Orbital angular momentum of a high-order Bessel light beam,” J. Opt. B: Quantum Semiclass. Opt. 4, S82–S89 (2002). 40. A. O’Neil, I. MacVicar, L. Allen, and M. Padgett, “Intrinsic and extrinsic nature of the orbital angular momentum of a light beam,” Phys. Rev. Lett. 88, 053601 (2002). 41. D. G. Grier, “A revolution in optical manipulation,” Nature 424, 810–816 (2003). 42. V. Garcés-Chávez, D. McGloin, M. J. Padgett, W. Dultz, H. Schmitzer, and K. Dholakia, “Observation of the transfer of the local angular momentum density of a multiringed light beam to an optically trapped particle,” Phys. Rev. Lett. 91, 093602 (2003). 43. J. E. Curtis and D. G. Grier, “Structure of optical vortices,” Phys. Rev. Lett. 90, 133901 (2003). 44. K. Ladavac and D. Grier, “Microoptomechanical pumps assembled and driven by holographic optical vortex arrays,” Opt. Express 12, 1144–1149 (2004). #174096 $15.00 USD Received 10 Aug 2012; revised 30 Sep 2012; accepted 1 Oct 2012; published 9 Oct 2012 (C) 2012 OSA 22 October 2012 / Vol. 20, No. 22 / OPTICS EXPRESS 24305 45. S. Parkin, G. Knöner, T. A. Nieminen, N. R. Heckenberg, and H. Rubinsztein-Dunlop, “Measurement of the total optical angular momentum transfer in optical tweezers,” Opt. Express 14, 6963–6970 (2006). 46. S.-H. Lee and D. G. Grier, “Giant colloidal diffusivity on corrugated optical vortices,” Phys. Rev. Lett. 96, 190601 (2006). 47. C. Schmitz, K. Uhring, J. Spatz, and J. Curtis, “Tuning the orbital angular momentum in optical vortex beams,” Opt. Express 14, 6604–6612 (2006). 48. M. Reichert and H. Stark, “Circling particles and drafting in optical vortices,” J. Phys.: Condens. Matter 16, S4085–S4094 (2004). 49. K. Ladavac and D. G. Grier, “Colloidal hydrodynamic coupling in concentric optical vortices,” Europhys. Lett. 70, 548–554 (2005). 50. C.-S. Guo, Y.-N. Yu, and Z. Hong, “Optical sorting using an array of optical vortices with fractional topological charge,” Opt. Commun. 283, 1889–1893 (2010). 51. V. R. Daria, M. A. Go, and H.-A. Bachor, “Simultaneous transfer of linear and orbital angular momentum to multiple low-index particles,” J. Opt. 13, 044004, (2011). 52. T. Čižmár, D. I. C. Dalgarno, P. C. Ashok, F. J. Gunn-Moore, 


Introduction
An optical force arises from the interaction between light and a microparticle as the result of light scattering by the microparticle.The optical force is studied and employed within the framework of optical micromanipulation techniques and has been utilized, for example, in the following instruments: optical tweezers [1], holographic optical tweezers [2,3], Raman tweez-ers [4], optical cell sorters [5], optical stretchers [6], or optical pikotenzometers [7,8].The above mentioned applications mainly use the transfer of the linear momentum of light from laser beams to an object.In the case of more complex spatial light distributions the particle behavior is strongly determined by its size with respect to the characteristic field pattern.This phenomena is sometimes called the size effect and has been observed in various field geometries.The particle behavior in the standing wave represents the simplest case [9].Here the particle is pushed with its center to the intensity maximum or minimum depending on its size.However, particles of particular sizes are not pushed at all and the overall optical force is negligible.Such strong dependence of the optical force on the particle size has been employed in various methods of passive optical sorting of microparticles in one-and two-dimensional optical lattices [10][11][12][13][14].
Except the linear momentum, light can posses also spin and orbital angular momentum [15].Spin angular momentum is associated with the polarization of light [16] and its change, for example due to the light transmittance through a birefringent microobject, results in a torque rotating the microobject around its axis [17,18].It has been pointed out that there exists the orbital angular momentum which is associated with the spatial field distribution in optical vortex beams [15,19].The transfer of the angular momentum from a vortex laser beam to an object and its subsequent behavior has been mainly investigated theoretically at the level of atoms, nanoparticles and molecules [20][21][22][23][24][25][26][27] or particles much larger that the light wavelength within the scope of the ray optics [28,29].More complex theoretical approaches, such as generalized Lorenz-Mie theory [30][31][32], finite-difference time-domain approach [33], finite element method [34], or multidipole approximation of the particle [35], are used less frequently but due to their wider range of applicability enable to express optical forces and torques acting upon a dielectric or metal particle of sizes comparable to the laser wavelength and illuminated by beams of structured spatial field distribution (e.g.optical vortex beams).The experimental verification of the transfer of angular momentum upon a particle has been demonstrated by rotation of an absorbing particle in the dark center of the vortex beam [36,37] or orbiting of microparticles around the beam axis in the high-intensity ring of the vortex beam [19,[38][39][40][41][42][43][44][45][46][47][48][49][50][51][52].Even though in these cases the sizes of such particles were comparable to the laser beam wavelength, no quantitative comparison between the theory and the experiment has been performed.Outside the optical domain the transfer of angular momentum from acoustic waves upon a macroscopic object led to its rotation, too [53].
Majority of examples referenced above focused either on experimental or theoretical aspects of angular momentum transfer.In this paper, we combine both theoretical and experimental approaches, we present a parametric theoretical study of the behavior of a single particle in the high-order Bessel beam and we compare the theoretical predictions with experimental observations.As the principle novelty, we identified three different regimes of particle behavior: the particle can orbit along the high intensity ring of the vortex beam, it can be stably localized at the dark spot placed on the vortex beam axis or at one of two points placed off the optical axis.

Beam description and calculation of the optical forces
In this paper we assume an ideal high-order Bessel beam generated behind an axicon that is illuminated with a plane wave [54][55][56][57][58][59][60].In the case of a plane wave linearly polarized along the x-axis the vector electric field of such Bessel beam (BB) can be described as [61,62]: where and e x,y,z are the base vectors along x, y, z Cartesian coordinate axes, ρ = x 2 + y 2 is radial distance from beam axis, φ is the azimuthal angle, α 0 is a semiapex angle of a cone along which the BB is formed, k and k r = k sin(α 0 ) are the wave vector and its projection into radial direction, respectively, and m is the topological charge of the optical vortex beam.The zeroorder BB corresponds to m = 0.For practical reasons we define the radius ρ m of the vortex beam with |m| > 0 using the radius of the first maximum of the intensity in the radial direction, obtained from dJ m (k r ρ m )/dρ m = 0.The vortex beam radius ρ m can be related to the radius ρ 0 of the core of the zero-order BB which is defined as the radius of the first intensity minimum in the radial direction and obtained from J 0 (k r ρ 0 ) = 0.For the topological charges up to 5 we obtain ρ 1 = 0.7656 ρ 0 , ρ 2 = 1.27 ρ 0 , ρ 3 = 1.747 ρ 0 , ρ 4 = 2.2112 ρ 0 , ρ 5 = 2.6678 ρ 0 , where [63] ρ In the following parametric studies we assume the zero-order BB core radii in the range ρ 0 = 0.3 − 1.5 μm.
The power carried by the innermost high intensity core of the BB for m = 0, ±1, ±2, ... can be expressed in the paraxial case [63] as where ω 0 is the light angular frequency, μ 0 is the vacuum permeability, σ m is the first off-axis root of the Bessel function of the m-th order (e.g.σ 0 = 2.4048, σ 1 = 3.8317, σ 2 = 5.1356, etc.).The minus sign in Eq. ( 4) is compensated by opposite signs of terms J m−1 (σ m ) and J m+1 (σ m ).Therefore, the power carried by the central ring of the BB of topological charge m = 1 is 1.53 × bigger than the power carried by the zero-order BB core.In the case of higher topological charges we obtain P 2 /P 0 = 1.95,P 3 /P 0 = 2.32 etc.
In coincidence with our experimental work, we determine the value of the electric field intensity E 0 (α) from the power P 0,core carried by the central core of the zero-order BB.This approach provides a direct link between the parameters of the idealized BB used for the calculations and the experimental realizations of such beams because P 0,core and ρ 0 can be measured experimentally.The typical experimental value corresponds to P 0,core = 5 mW that we use to obtain the value of E 0 for all the calculations presented in this paper independently on their core radius ρ 0 and topological charges m ≥ 1.This approach is based on the way how we generated experimentally the vortex beams of different topological charges using the spatial light modulator.First we generated zero-order BB and then we switched to high-order BBs.
Using Eq. ( 1) and following the same approach as in Ref. [62] we present analytical formulas for optical forces acting upon a particle much smaller than the trapping wavelength (i.e.induced dipole or Rayleigh particle): where α and α denotes the real and imaginary part of the particle polarizability [62], respectively, ε 0 is the permittivity of vacuum, ε 1 is the relative permittivity of the surrounding medium and we skipped k r ρ in Bessel functions.Even though these equations are of limited validity they provide insight into the behavior of tiny particles.Terms related to α give rise to gradient force that due to a 3 dependence represents the leading force acting upon tiny particles.Terms related to α correspond to scattering force.It depends on the particle radius as a 6 and, therefore, it becomes more pronounced for larger particles (within the validity of this Rayleigh approximation).Terms with P reflects the non-paraxiality of the beam and their influence increases for larger α 0 (i.e. for narrower BB core ρ 0 ).Nice example of the competition between the gradient and scattering force can be shown for F φ .Neglecting all terms with P we end up with only scattering azimuthal force F φ = ε 0 ε 1 α E 2 0 mJ m /(2r) that evidently causes particle orbiting.However, due to a 6 dependence this force is weak for smaller particles and can be overcome by gradient force (the first term in Eq. (6).The sine term in gradient force leads to two azimuthal equilibrium positions placed off the beam axis, as will be discussed below.
Optical forces acting upon a spherical dielectric particle of any size we calculated using the generalized Lorenz-Mie theory [63][64][65][66][67][68].To speed up the numerical calculations presented in the following section we utilized the results of Taylor [69] and we expressed the scattered field coefficients A ln and B ln analytically for the vortex BB [60].It has shortened the computation time in the Matlab environment by two orders of magnitude.

Numerical results
We considered polystyrene particles of refractive index n 2 = 1.59 and radii a in the range 10 nm -1.5 μm surrounded by water (refractive index n 1 = 1.334).Besides water we also considered air as the surrounding medium (n 1 = 1) because optical manipulation in air offers much lower friction and probably represents upcoming direction of further development of this technique.The radial and axial optical forces were calculated at a single axial position of the particle considering BB beam core radii ρ 0 in the range 0.3 − 1.5 μm and topological charges m = 1, 2, 3, 4, 5, 10. Figure 1 shows the force acting upon a polystyrene particle of radius 1 μm placed into the BB having the following core radii ρ 0 = 1, 1.25, 5 μm.The radial component of the force is depicted in Fig. 1(a) as a function of the radial coordinates ρ along parallel (solid) and perpendicular (dashed) direction to the beam polarization.It can be seen that the particle center is stably trapped at the bright part of the vortex fringe (ρ/ρ 1 = 1) only for the widest core radius (ρ 0 = 5 μm), i.e. when a ρ 1 .Moreover, the red curve in Fig. 1(b) shows that the azimuthal force is almost constant and, thus, the particle will orbit along the vortex ring (regime R1).In the case of narrower beam core radius (ρ 0 = 1.25 μm), i.e. the particle radius is still smaller but comparable to ρ 1 , the particle is trapped off-axis at the radial distance ρ > 0, however, this distance is smaller than the vortex beam radius and the azimuthal component of the optical force changes its sign.Therefore, the particle does not orbit along the vortex ring but is stably trapped at certain azimuthal position (regime R2).If the beam core radius is even narrower (ρ 0 = 1 μm), the particle overlaps the first bright vortex fringe of radius ρ 1 and it is trapped with its center at the vortex beam axis, i.e. ρ = 0 (regime R3).In this case the azimuthal force is equal to zero and thus no curve associated to this regime is plotted in Fig. 1(b).Furthermore, Fig. 1(a) demonstrates, that as the beam core radius decreases (i.e.α 0 increases), the non-paraxial vectorial properties of the optical vortex beam become apparent and the radial force component differs along various azimuthal directions (parallel and perpendicular to the beam polarization).
All three different regimes of particle's behavior are summarized in Table 1 and links to subsequent Figs. 1, 2, 3 are established.The results of calculations in the xy plane are visualized in Fig. 2 for the polystyrene particle of radius a = 1 μm and six different radii of the BB beam core ρ 0 .The size of the particle is compared to the radius of the first bright vortex beam fringe and thus the direct link between the particle size and characteristic dimension of the vortex beam can be established for each regime.The particle is initially trapped in the innermost bright fringe while it orbits there (regime R1, see Fig. 2(a)), similarly as in Fig. 1.As the core radius decreases the particle stops orbiting but still being located off the vortex beam axis (Fig. 2(b), regime R2).Two stable positions exist here for the particle, see green curve in Fig. 1 as well.For even thinner core the particle locates itself with its center on the vortex beam axis, i.e. into the dark, light-free region (Fig. 2(c), regime R3).While reducing the core radius even further, i.e. when the particle radius is getting much larger than the radius of the first bright vortex fringe of radius ρ 1 , particle's behavior described by these three regimes repeats(see Figs.2(d-f)).
Figure 3 uses the particle and the BB core radii as the free parameters to distinguish all three different regimes R1, R2, R3 of the particle behavior.The force calculations were performed for beam core values ρ 0 varied in steps 10, 25 or 50 nm, however, near the borders between  The optical force acting upon the polystyrene particle of radius 1 μm in the radial direction.Solid curves show the force along the beam polarization (i.e.x axis) dashed curves show the force in direction perpendicular to the beam polarization (y axis).Red (regime R1), green (regime R2), and blue (regime R3) curves show the forces for vortex beam of corresponding core radii ρ 0 = 5 μm, ρ 0 = 1.25 μm and ρ 0 = 1 μm, respectively.Note that the red curve is multiplied by factor 10. (b) The azimuthal optical force acting upon the same particle placed in the stable radial distance r denoted in Fig. 1(a), i.e. r = ρ 1 = 0.7656ρ 0 = 3.828 μm for beam core radius ρ 0 = 5 μm (red curve) and r = 468 nm for ρ 0 = 1.25 μm (green curve).
all regimes the steps were decreased to 5 nm to ensure the proper border placement.We can see that very small particles (a ≤ 50 nm) are always trapped off-axis in regime R2 which is in agreement with Eq. ( 6).This is caused by the fact that such particles (basically elementary dipoles) are strongly influenced by the azimuthal optical gradients along the vortex innermost rings that arises from the non-paraxial terms with P ⊥ in the beam description.One can also see that this region gets narrower as the core radius increases (i.e.α 0 decreases) and the azimuthal inhomogeneity in optical intensity becomes less pronounced.Once the optical force associated to the optical angular momentum of the vortex ring (it is related to azimuthal scattering force in the case of Rayleigh particles) is strong enough for larger particles and such particle starts to circulate along the innermost fringe (regime R1).At certain particle radii a (∼ 0.82ρ 0 ) the particle either stops in R2 or jumps directly onto the vortex axis into R3.Even bigger particle (a 1.7ρ 0 ) jumps again off axis either into the orbiting regime R3 or into the stable off-axis position R2.Behavior of particles surrounded by air follows similar trend as those immersed in water, however, the boundaries between the regimes are less smooth due to the to larger contrast in refractive indices of air and particle giving rise to stronger morphological dependent resonances.
Figure 4 compares the off-axial equilibrium positions r of the particle in regimes R1 and R2 relatively to the radius ρ 1 of the inner-most bright intensity vortex fringe.Particles much smaller then ρ 1 are confined at radial positions very close to the first high intensity ring of the vortex beam, i.e. r/ρ 1 1, because they are under the dominant influence of the intensity gradient near the bright region of the vortex fringe.As the radius of the particle increases and approaches the beam core radius, the particle starts to overlap larger volume of the vortex beam in contrast to only narrow part near the bright vortex fringe, and thus its radial position gets closer to the beam axis.Particles of radius overlapping several vortex fringes are deviated less from the on-axial position (r/ρ 1 0.5 in water and air).

Regime R1: orbiting particles
We have calculated and analyzed the particle trajectories for the orbiting regime R1 in order to find the average angular velocity ω of the particle movement.We have tracked the particle motion for several orbits and used the Euler formula to solve the equation of motion in the Fig. 2. Optical forces and trajectories of a particle of radius 1 μm placed in an optical vortex beam of topological charge m = 1 having various BB core radii (a: ρ 0 = 5 μm, b: ρ 0 = 1.25 μm, c: ρ 0 = 1 μm, d: ρ 0 = 0.58 μm, e: ρ 0 = 0.5 μm, f: ρ 0 = 0.45 μm).The background pseudo-color plot shows the electric field intensity |E| 2 normalized relative to the maximal intensity in (f).The magenta curves denote the deterministic trajectories of a particle (i.e.without considering the Brownian motion) starting at different locations and following the particle motion towards an equilibrium point or a stable orbit.The blue circle depicts the particle edge, its center is shown by the blue dot.The black and cyan contour represents the zero forces in the radial and azimuthal directions, respectively.If the magenta trajectories follow the black curve (zero radial force), see (a) and (d), this black curve forms a set of equilibrium positions and the particle orbits along the black curve (the particle center is drawn just in one selected position).If black and cyan curves intersect there exist equilibrium positions of the particle off the vortex axis, see (b,d).One of such possible stable positions of the particle center is denoted by the full blue dot.In other cases, see (c,f), the particle is trapped with its center on the vortex beam axis.over-damped case using v = F/γ where F is the optical force acting upon the particle, v is the particle velocity and γ = 6πηa is Stokes drag coefficient, η is the kinematic viscosity of the medium.In most cases the particle orbited with constant speed and we fitted these trajectories by the sinusoidal function x(t) = A sin(ωt) + B in order to obtain the angular velocity ω.The results are shown in Fig. 5.However, in several cases, mostly near the border between orbiting R1 and off-axis stable R2, the particle velocity along the orbit varied significantly and we used the period of the particle orbit T to determine the average angular velocity ω = 2π/T .Figure 5 shows that tiny particles that just started their orbiting move with the lowest angular velocity.ω further increases significantly (note the logarithmic scale of the Fig. 5) for the small radii of the beam core due to the higher optical intensity in the first bright fringe of these beams.However, this high angular velocity sustains even for larger particles that switched back to the rotation after being located on the vortex axis.Figure 5(b) stresses that due to the lower friction in air the particles orbit here with angular velocities two orders of magnitude higher comparing to motion in water under the same conditions.

Regime R2: off-axis lateral confinement
The particle placed in regime R2 may be located in one of two locations placed symmetrically with respect to the vortex axis.The angle Φ between radius vector of the first position and the x axis (corresponding to the direction of the beam polarization) is depicted in Fig. 6.Smaller particles are trapped at positions that are perpendicular to the beam polarization, however as the angular force increases with particle radius, Φ increases until the particle starts its orbiting in regime R1.For particles of radius larger than the core of the beam abrupt changes of Φ appear for water.In contrast particles of all sized surrounded by air are confined mainly at positions with Φ 90 • .

Axial force
In all regimes mentioned above the particle is propelled along the beam propagation axis by the axial force F z .Figure 7 shows the numerical results corresponding to the axial force F z acting upon the particle when it is laterally localized (regimes R2 and R3).In the case of orbiting regime R1 the average axial force is shown.The results show that the particle is axially more strongly propelled if it is settled in regimes R3 or R2.For larger particles the force is stronger and less dependent on the particle positions because the particle overlaps several vortex beam fringes.regime (R1) is preferable for higher topological charges while the off-axis stable regime (R2) vanishes there.Since with increasing m and fixed ρ 0 the radius of the first high-intensity ring increases, we have observed no trapping at the vortex core for m = 10 for the selected particle radius.

Higher topological charges
The normalized radial distance ρ/ρ m of the particle center from the vortex axis is depicted in Fig. 8(b) for the off-axis regimes R1 and R2, where ρ m is the radius of the inner-most bright vortex fringe.The particle orbits or is radially fixed mainly at a distance lower than the corresponding ρ m .This distance slowly decreases under the conditions when the particle tends to move to the axial position (R3).However, for vortex beams of higher topological charge having small core radius (ρ 0 < 0.75 μm) the particle may orbit along a trajectory having radius even larger than ρ m .This is caused by the fact that the distance between the outer intensity fringes decreases for higher topological charges and the particle covers more light intensity while being settled between the first and the second inner-most intensity fringe.
Figure 8(c) shows the average angular velocity ω of the particle.One can see, that ω decreases for increasing topological charge and constant ρ 0 .This result is slightly against the well accepted idea that optical vortex beams of higher topological charge produces faster angular motion due to larger orbital angular momentum.However, since we keep the power in the BB core constant (P 0,core = 5 mW), optical intensity in the ring of the high-order BB decreases as m raises due to the increased ring radius and, therefore, the density of the orbital angular momentum decreases.For example, if we kept the radius of the innermost bright fringe constant (e.g.ρ 1 = ρ 2 = 1.15 μm) for different topological charges, we obtain ρ 0 = 1.5 μm and ρ 0 = 0.9 μm corresponding to m = 1 and m = 2, respectively.Comparing with Fig. 8(c) one can see that for these values of the beam cores ρ 0 , the particle orbits with much higher angular velocity for topological charge m = 2.This trend gets even more pronounced for higher vortex topological charges.

Experimental beam generation and measurement procedures
The BB of the zero-order we formed behind a lens (i.e. in the Fourier space) using a ring-like blazed phase diffraction grating imposed on the spatial light modulator (SLM, Holoeye LC-2500R).The width of its core was determined by the semiapex angle α 0 which was controlled by the radius of the diffraction grating in the form of a ring.The BB of topological charge m was generated by adding an azimuthally linearly increasing phase from 0 to 2mπ to the previous grating.Since the SLM enabled real-time modification of such diffraction grating, it gave us the freedom to generate non-diffracting beam of different orders and also widths.Figure 9 introduces how this key optical element was placed in the experimental setup.The incoming laser beam (Verdi V10, Coherent, λ vac = 532 nm) was spatially filtered using the achromatic doublet L 1 ( f 1 = 19 mm, Thorlabs AC127-019-A) and the pinhole of diameter 10 μm (Thorlabs P10C).The beam was collimated by the achromatic doublet L 2 (Thorlabs AC254-200-A) of focal length f 2 = 200 mm so that the beam fully overlapped the SLM chip.The grating imposed on the SLM diffracts the incident beam and only the first diffraction order was used for subsequent experiments.The zeroth diffraction order was blocked in the focal plane of the lens L 3 (Thorlabs AC508-750-A-ML) of focal length 750 mm while the first diffraction order was collimated by lens L 4 (Thorlabs AC254-300-A, f 4 = 300 mm) and entered the aspheric lens used as the objective (Thorlabs C240TME-A, f O = 8 mm).The lens L 4 and the objective served as a telescope that demagnified the width of the BB approximately 40× in the sample space.The sample space was filled with deionized water with dispersed polystyrene microspheres of radius 1 μm (DukeScientific 4K-02) or 2 μm (DukeScientific 4K-04) and illuminated by the BB.The microspheres behavior was observed with the planachromat microscope objective (Olympus PLCN 60×, NA 0.8) and recorded with the fast CCD camera (Basler piA640-210gm).These images were post-processed and the microspheres positions were obtained in nanometer precision (using the calibrated graticule giving 107 nm/pixel).
During the single experimental procedure, the same microsphere was illuminated by the high-order BB and the microsphere positions were recorded.The beam center was found using the zero-order BB that confined the microsphere at the beam center.Both, the high-order and the zero-order BBs corresponded to the same radius of the diffraction ring imposed on the SLM.The radial positions of the particle in high-order and zero-order BBs were recorded over 20 successive measurements.The above described experimental procedure was repeated for BBs of topological charges m = 1 − 5 and for different radii of the diffraction ring imposed on the SLM corresponding to different radii ρ 0 of the BB cores in the sample plane.
Further on, the lateral profiles of the beam optical intensity were recorded for each of the above mentioned combinations when a region with no microspheres was illuminated and the notch filter was removed from the imaging path.The measured profiles were fitted with the expected theoretical ones (see Eq. Fig. 10.Examples of particle (radius 2 μm) motion when the BB is changed from the zero-order to high-order with m = 1 − 5, ρ 0 = 870 nm in all cases.Top row: Theoretical predictions of particle trajectories following the conventions of Fig. 2. Bottom row: Experimental observations.Background plot shows the measured intensity profile of the vortex beam.Yellow spots denote the particle trapped in the zero-order BB core and magenta curves show the motion of the particle when its is illuminated by the high-order BB.The beams are polarized along horizontal axis.
for all studied topological charges.Figure 10 compares the theoretically expected behavior to the experimental observations.The figure shows examples of particles trajectories in the BB corresponding to ρ 0 = 870 nm when the topological charge was changed from m = 0 at the beginning to m = 1 − 5.The yellow spots in the bottom row indicate places where the particle is trapped in the zero-order BB and the magenta curves denote the trajectories of the particle to its new equilibrium position when the high-order BB is established.The particle remains near the vortex beam axis for the topological charges m = 1, 2, 3 while it moves to the bright ring for topological charges m = 4, 5.This behavior coincides with the theoretical calculations that predict that the particle should stay on the beam axis (i.e. in the regime R3) for topological charges m = 1 − 3 while it should orbit around the bright fringe in regime R1 for topological charges m = 4, 5.We have not observed particle orbiting over the whole circle even though we used the in-situ aberrations correction method [70] to approach the vortex beam profile expected theoretically as close as possible.Comparison of the BB intensity profiles shown at the top and bottom rows in Fig. 10 indicates that we were not able to suppress the aberrations completely and the intensity variation along the vortex fringe survived.Such variation probably caused that we have not been able to observe the orbiting of the single trapped particle (regime R1) but we only observed particle motion towards the vortex fringe followed by its motion less then half a circle along the fringe.The aberrations became more pronounced for small diameter vortices of low topological charges 1 or 2 while their influence decreased for vortices of charges 3 to 5. The main contribution to imperfect aberration corrections comes from the poor surface flatness of the SLM chip [71] Since we were able to distinguish experimentally whether the particle settled into the vortex dark core (R3) or off the beam axis, we were able to determine the particle distance from the vortex beam axis for different core radii and topological charges.Figure 11 shows the average distance of the particle center from the vortex axis relative to the radius of the innermost vortex fringe.The left and right columns show the experimental data for the polystyrene particles having radius 1 and 2 μm, respectively.The red curves show the theoretical prediction of the particle distance from the vortex axis, similarly to Fig. 8(b).
Due to the imperfect compensation of the beam aberrations discussed above, the coincidence between the experimental results and the theoretical predictions is worse for m = 1 and gets better for m = 2 and a = 1 μm.For other parameters the coincidence is slightly worse for lower ρ 0 where the radius of the innermost vortex beam ring is in the range of 1 to 3 wavelengths (in water) because small disturbances in the beam propagation cause intensity or phase inhomogeneities that are not considered in the theoretical description of the beam and cause differences between the predicted and the observed particle behavior.Very good coincidence is obtained for the topological charges m = 3, 4, 5 where the actual radius of the vortex ring is much higher than in the previous cases.

Conclusion
The influence of the particle size on optical forces localizing such particle in the spatially structured beam has been already demonstrated for the zero-order BB, standing waves and 2D optical lattices.Our new results demonstrate that optical forces acting upon a particle illuminated by the high-order BB are strongly influenced by particle size with respect to the radius of the innermost bright vortex fringe.Using the generalized Lorenz-Mie theory we have identified three different regimes of particles behavior in the lateral plane that have been summarized in Table 1.The particle either orbits close to the high intensity ring of the vortex beam (regime R1), or it is trapped off the beam axis in one of two stable azimuthal positions (regime R2), or settles with its center in the zero intensity center of the beam (regime R3).We have shown how these regimes depend on the radius of the beam core, radius of the spherical particle, and topological charge of the beam.The smallest particles follow the R2 regime, larger particles switch to R1 and orbits in the vortex beam with certain angular frequency.In both cases the particles' distance from the beam center is smaller than the radius of the high intensity fringe.With increasing size of the particle with respect to the radius of the innermost vortex fringe the regimes switch between each other and the particles positions, angular frequency or axial optical force are influenced.The regime R1 becomes dominant for higher topological charges under the range of investigated parameters.Experimental investigations have been performed for 24 beam widths and 5 topological charges and the results coincide with the theoretical radial positions of the particles better for wider beams and higher topological charges (m > 2).The worse coincidence was observed for lower topological charges and narrower cores due to imperfect compensation of the beam aberrations.

Table 1 .
Different regimes of particle's behavior in non-diffracting vortex beam c,f) blue in the dark beam center

Fig. 1 .
Fig. 1. (a)The optical force acting upon the polystyrene particle of radius 1 μm in the radial direction.Solid curves show the force along the beam polarization (i.e.x axis) dashed curves show the force in direction perpendicular to the beam polarization (y axis).Red (regime R1), green (regime R2), and blue (regime R3) curves show the forces for vortex beam of corresponding core radii ρ 0 = 5 μm, ρ 0 = 1.25 μm and ρ 0 = 1 μm, respectively.Note that the red curve is multiplied by factor 10. (b) The azimuthal optical force acting upon the same particle placed in the stable radial distance r denoted in Fig.1(a), i.e. r = ρ 1 = 0.7656ρ 0 = 3.828 μm for beam core radius ρ 0 = 5 μm (red curve) and r = 468 nm for ρ 0 = 1.25 μm (green curve).

#(Fig. 3 .Fig. 4 .Fig. 5 .(Fig. 6 .Fig. 7 .
Fig.3.Phase map summarizing three different regimes of behavior of polystyrene particles as a function of the particle radii a and the BB core radii ρ 0 in the optical vortex beam of topological charge m = 1.The particle orbits along the circular trajectories (red areas, R1, it corresponds to cases in Fig.2(a,d)), settles in one of two off-axis positions (green areas, R2, similar to cases in Fig.2(b,e)) or in the dark center of the vortex beam (blue areas, R3, it corresponds to cases in Fig.2(c,f)).The particle is surrounded either by water (a) or air (b).

Figure 8 Fig. 8 .
Figure 8  demonstrates the behavior of a polystyrene particle (a = 1 μm) placed into the vortex beams of m = 1, 2, 3, 5, 10 and corresponding BB core radii ρ 0 in the range 0.3 -1.5 μm.The phase curves identifying all three regimes are depicted in Fig.8(a).Obviously, the orbiting
(1)) to find the beam parameters.This procedure was repeated

Fig. 11 .
Fig.11.Measured stable radial position r of the particle in the vortex beam relative to the radius of the innermost vortex ring ρ m of different corresponding ρ 0 .Blue (positive topological charge m > 0) and green (negative topological charge m < 0) points correspond to the measured data, the error-bars indicate 95 % confidence level of the average.The red curves denote the theoretical prediction for the measured beam parameters.The BBs of various corresponding beam core radius ρ 0 and topological charges m = 1 − 5 (rows from the top to the bottom) were generated by the setup shown in Fig.9.The left and right column corresponds to polystyrene particle of radii 1 μm and 2 μm, respectively.