Next Article in Journal
Rational Design of Non-Noble Metal Single-Atom Catalysts in Lithium–Sulfur Batteries through First Principles Calculations
Previous Article in Journal
The Enhanced Performance of Oxide Thin-Film Transistors Fabricated by a Two-Step Deposition Pressure Process
Previous Article in Special Issue
Chiral Mechanical Effect of the Tightly Focused Chiral Vector Vortex Fields Interacting with Particles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optical Force Effects of Rayleigh Particles by Cylindrical Vector Beams

1
School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025, China
2
Moji-Nano Technology Co., Ltd., Yantai 264006, China
*
Authors to whom correspondence should be addressed.
Nanomaterials 2024, 14(8), 691; https://doi.org/10.3390/nano14080691
Submission received: 25 March 2024 / Revised: 11 April 2024 / Accepted: 15 April 2024 / Published: 17 April 2024
(This article belongs to the Special Issue Advances in Optical Nanomanipulation)

Abstract

:
High-order cylindrical vector beams possess flexible spatial polarization and exhibit new effects and phenomena that can expand the functionality and enhance the capability of optical systems. However, building a general analytical model for highly focused beams with different polarization orders remains a challenge. Here, we elaborately develop the vector theory of high-order cylindrical vector beams in a high numerical aperture focusing system and achieve the vectorial diffraction integrals for describing the tight focusing field with the space-variant distribution of polarization orders within the framework of Richards–Wolf diffraction theory. The analytical formulae include the exact three Cartesian components of electric and magnetic distributions in the tightly focused region. Additionally, utilizing the analytical formulae, we can achieve the gradient force, scattering force, and curl-spin force exerted on Rayleigh particles trapped by high-order cylindrical vector beams. These results are crucial for improving the design and engineering of the tightly focused field by modulating the polarization orders of high-order cylindrical vector beams, particularly for applications such as optical tweezers and optical manipulation. This theoretical analysis also extends to the calculation of complicated optical vortex vector fields and the design of diffractive optical elements with high diffraction efficiency and resolution.

1. Introduction

In recent years, there has been an extensive exploration of cylindrical vector (CV) beams due to their spatially variant polarization orders [1,2,3,4,5,6]. Researchers have developed numerous systems for generating CV beams with exotic properties [7,8,9,10,11,12] and have applied them in various applications, such as optical manipulation [13,14,15,16,17,18,19], the spin-orbit hall effect [20,21,22], multiplex communication [23], and super-resolution imaging [24]. These applications have sparked growing interest in the high numerical aperture (NA) focusing on the characteristics of CV beams and their properties in optical capture and manipulation. CV beams are vector-beam solutions of Maxwell’s equations that exhibit axial symmetry in both amplitude and phase. Furthermore, the additional degree of freedom in high-order CV beam allows for engineering the characteristics of the focal field to tailor the distribution of the optical force, thereby addressing specific requirements of optical tweezer systems. In 2004, Zhan demonstrated that metallic Rayleigh particles with strong scattering and absorption forces could be trapped using tightly focused radial polarization beams [25]. In 2017, Li et al. integrated high-order CV beams with variable polarization topological charge to analyze the spin momentum and torque exerted on the trapped particles [26]. Moreover, forces partially associated with the spin of the internal energy flow can induce particle rotation and influence their motion [27,28]. The aforementioned examples underscore the significance of studying vector diffraction theory near the focus of high-order CV beams.
Thus, the different analytical models of vector formalism in tightly focused fields have been proposed further to explore the focusing characteristics of highly focused CV beams. To achieve high diffraction efficiency and high resolution of diffractive optical elements, computational approaches involving various vector models have been utilized in the vicinity of the focus [29,30]. Rui et al. presented a novel strategy to manipulate the captured metallic nanoparticles by using an engineered azimuthally polarized optical field and gave the complete electric field expression to predict the motion behavior of the trapped nanoparticles [31]. However, it does not include the case of higher-order polarization. Additionally, Chen et al. proved that an optical pulling force could pull a particle toward the source via a backward scattering force [32], but the complete analytical models for the fast calculation were not provided. In summary, there is still a lack of an analytical model based on rigorous diffraction theory to further investigate the multi-dimensional force aroused by the CV beams, including arbitrary polarization orders.
In this paper, we elaborately calculate the analytical formulae of electric and magnetic field distributions based on the Richards–Wolf vector diffraction integral theory [33]. In addition, the gradient force, scattering force, and curl-spin force induced by the high-order CV beams are calculated in detail. Further, we find that high-order polarization vector beams can produce pulling force along the axis, which helps understand the propagation behavior of scattering force in tightly focused high-order beams. The universality of the expressions in arbitrary polarization order CV beams simplifies multi-parameter calculation in the tightly focused region.

2. Theoretical Model

2.1. Calculation of Cylindrical Vector Beam with Arbitrary Polarization Distribution

As depicted in Figure 1, we assume that an exemplar paraxial CV beam passes through an aperture with radius R and is focused by a Fourier lens. A Rayleigh particle with radius a is captured near the focal region, where a is significantly smaller than the wavelength of the beam.
The paraxial cylindrical vector beam propagating in the direction of the z axis has a transverse electric field in the x–y plane perpendicular to the direction of propagation, which can be expressed as follows:
E i = cos ( m ϕ + ϕ 0 ) e x + sin ( m ϕ + ϕ 0 ) e y ,
where m is the polarization order, ϕ is the azimuthal angle, and the initial phase ϕ0 is constant. The ex and ey are the unit vectors in the x and y directions.
Any polarization order of the beam in a transverse field can be decomposed into two basic polarizations: radial and azimuthal. Due to the polarization singularity at the center, different CV beams exhibit a doughnut-shaped intensity distribution [34]. The weight of the radial and azimuth polarization components varies along the azimuthal angle. Figure 2 shows the polarization distribution of the CV beams with m = 0, 1, −1, 2, and ϕ0 = 0. It demonstrates that the polarization structure is determined by the polarization orders, which dictate the spatial distribution and weight of radial and azimuthal polarization components. Thus, the different types of CV beams can be simulated by adjusting the polarization order m. In the transmission field,
e p = cos ϕ sin ϕ 0 , e s = sin ϕ cos ϕ 0 , e r = cos ϕ cos θ sin ϕ cos θ sin θ ,
where ep, es, and er are the unit vectors as shown in Figure 1, and θ is the deflection angle. Then, the transmission field can be expressed as follows:
E t ( θ , ϕ ) = t p ( E i e p ) e r + t s ( E i e s ) e s ,
by substituting from Equations (1) and (2) into Equation (3), we can obtain the following:
E t = sin δ sin ϕ + cos δ cos θ cos ϕ sin δ cos ϕ + cos δ cos θ cos ϕ cos δ sin θ .
Here, δ = (m − 1) ϕ + ϕ0. For a high NA focusing system, the electric field vector in the focal area can be obtained by the diffraction integral of a. Here, a represents the ‘strength factors’ of the electric field, and the relation between a and the electric field in image space is as follows [33]:
a = f l cos θ E t ,
where f is the focal length, and l is the amplitude factor that depends on the power. On substituting Equation (4) into Equation (5), we obtain the following:
a x = f l cos θ ( sin δ sin ϕ + cos δ cos θ cos ϕ ) , a y = f l cos θ ( sin δ cos ϕ + cos δ cos θ cos ϕ ) , a z = f l cos θ cos δ sin θ .
The ‘strength factors’ b of the magnetic field is related by the formula:
b = κ ( s × a ) .
Here, κ = 1/iωμ0. ω is the angular frequency of the optical field, and μ0 is the permeability in vacuum. s is the unit vector along the ray in the image space. Its three components can be expressed as s = (sinθcosφ, sinθsinφ, cosθ). Substituting Equation (6) into Equation (7), we obtain the Cartesian components of the strength vectors b:
b x = κ f l cos θ ( cos δ sin ϕ sin δ cos θ cos ϕ ) , b y = κ f l cos θ ( cos δ cos ϕ sin δ cos θ cos ϕ ) , b z = κ f l cos θ ( sin δ sin θ ) .
According to the Richards–Wolf vector diffraction integral theory [33], the electric and magnetic fields for an aberration elimination system at any point near the focal point can be written as follows:
E x , y , z = 1 i λ Ω a s x , s y s z exp ( i k s x x + s y y + s z z ) d s x d s y , H ( x , y , z ) = 1 i λ Ω b ( s x , s y ) s z exp ( i k s x x + s y y + s z z ) d s x d s y ,
where k = 2π/λ is the wave number, and λ is the wavelength in the vacuum.
To describe the electric and magnetic fields in detail, the quantity dsxdsy/sz, which enters the basic diffraction integral, represents the element dΩ of the solid angle and is given by dsxdsy/sz = dΩ = sinθdθdφ. Substituting from Equation (6) into Equation (9), we obtain the following expressions for the Cartesian components Ex, Ey, and Ez:
E x = 0 2 π 0 α u ( sin δ sin ϕ + cos δ cos θ cos ϕ ) d ϕ d θ , E y = 0 2 π 0 α u ( sin δ cos ϕ + cos δ cos θ sin ϕ ) d ϕ d θ ,   E z = 0 2 π 0 α u ( cos δ sin θ ) d ϕ d θ ,
where α is the maximum allowed incident angle determined by the numerical aperture of the objective lens, and u = −iπfl/(πλ)sinθcos1/2θexp{ik[rsinθcos(ϕ − φ) + zcosθ]}. In the same way, on substituting from Equation (8) into Equation (9), we obtain the following expressions of magnetic fields for the Cartesian components Hx, Hy, and Hz:
H x = κ 0 2 π 0 α u ( cos δ sin ϕ sin δ cos θ cos ϕ ) d ϕ d θ , H y = κ 0 2 π 0 α u ( cos δ cos ϕ sin δ cos θ cos ϕ ) d ϕ d θ , H z = κ 0 2 π 0 α u ( sin δ sin θ ) d ϕ d θ .
The integration concerning ϕ can immediately be carried out with the help of the following formulae, which are valid for any integral value of n [35]:
0 2 π exp ( i n ϕ ) exp [ i k r sin θ cos ( ϕ φ ) ] d ϕ = 2 π i n J n ( k r sin θ ) exp ( i n φ ) .
Here, Jn is the Bessel function of the first kind with order n.
From Equations (10)–(12), we finally obtain the following expressions for the components of the electric field vectors in the focused field:
E x = i m 1 e i ϕ 0 ( I 0 + I 1 ) i 1 m e i ϕ 0 ( I 2 + I 3 ) , E y = i m e i ϕ 0 ( I 1 I 0 ) i m e i ϕ 0 ( I 3 I 2 ) , E z = i m e i ϕ 0 I 4 i m e i ϕ 0 I 5 ,
H x = κ [ i m e i ϕ 0 ( I 0 I 1 ) + i m e i ϕ 0 ( I 2 I 3 ) ] , H y = κ [ i m 1 e i ϕ 0 ( I 0 + I 1 ) + i 1 m e i ϕ 0 ( I 2 I 3 ) ] , H z = κ [ i m 1 e i ϕ 0 I 4 + i 1 m e i ϕ 0 I 5 ] ,
where
I 0 = 0 α u 1 ( 1 + cos θ ) J m ( k r sin θ ) d θ , I 1 = e ( m 2 ) i φ 0 α u 1 ( 1 cos θ ) J m 2 ( k r sin θ ) d θ , I 2 = e ( 2 m ) i φ 0 α u 1 ( 1 cos θ ) J 2 m ( k r sin θ ) d θ , I 3 = e i m φ 0 α u 1 ( 1 + cos θ ) J m ( k r sin θ ) d θ , I 4 = 2 e ( m 1 ) i φ 0 α u 1 sin θ J m 1 ( k r sin θ ) d θ , I 5 = 2 e ( 1 m ) i φ 0 α u 1 sin θ J 1 m ( k r sin θ ) d θ ,
with u1 = π fl/(2λ) sinθcos1/2θexp(ikzcosθ). Hence, we can calculate the distribution of the electromagnetic field of a high-order CV beam with arbitrary polarization orders in the tight focusing system. Therefore, the distribution of optical forces can be obtained. We will discuss these characteristics in the following sections.

2.2. Theory of Optical Force on a Rayleigh Particle

It is assumed that a Rayleigh particle with a radius smaller than the wavelength is near the focusing region, and its motion is controlled by the force of light [36]. Furthermore, in non-magnetoelectric dipolar particles, there are no magnetic dipoles [37]. The particle can be considered simply as a dipole in the calculation [38,39]. The expression of the time-averaged optical force can be written as follows [40,41,42]:
F = 1 4 Re { α 1 } | E | 2 + σ c S + σ c × L S .
The first part of Equation (16) is the gradient force, which depends on the strength of the electric field, and it can be written as follows:
F g = 1 4 Re α 1 E 2 .
Here, α1 is the polarizability given by the following:
α 1 = α 0 1 i α 0 k 3 / ( 6 π ε 0 ) ,
with α0 = 4πε0a3[(ε − 1)/(ε + 2)], and ε is the relative permittivity of the nanoparticle, and a is the radius of the nanoparticle.
The second part in Equation (16) is the scattering force, and it is proportional to the time-averaged Poynting vector:
F s = σ c S ,
where S = 1 2 Re { E × H * } is the time-averaged Poynting vector. σ = kIm{α1}/ε0 is the total cross-section of the particle. c is the speed of light in vacuum.
The third part of the right hand in Equation (16) is the curl-spin force, which is associated with the spin density of the light field. The expression is as follows:
F c = σ c × L S ,
where L S = ε 0 4 ω i { E × E * } is the time-averaged spin density of a transverse electromagnetic field [40]. ω refers to the angular frequency of the time-harmonic. For low numerical aperture illumination, this term is usually ignored, and for plane wave illumination with linear polarization, this term is zero.
In summary, light forces on Rayleigh particles can be described as the sum of three terms: the gradient force, the radiation pressure force, and the curl-spin force. We will discuss these forces in detail by employing the derived formulae and the variation in these forces in vector beams with different polarization orders; then, the trapping behavior of Rayleigh particles in tightly focused systems can be fully verified theoretically.

3. Results and Discussion

3.1. Electric and Magnetic Fields of High-Order Vector Beams

Now, we investigate the electromagnetic intensity distribution of the light field in tightly focused systems. We use the derived integral of the electric field to simulate the electric field intensity distribution, as shown in Figure 3. It shows the Cartesian coordinate components and total electric field intensity distribution of CV beams with polarization orders of 0, 1, −1, and 2 in the focusing field. The simulations are for the objective with NA = 0.9, and the initial phase ϕ0 = 0. The objective filling factor was assumed to be one.
As shown in Figure 3, the electric field presents the different distributions of the different polarization orders. Ex, Ey, and Ez represent the x-component, the y-component, the z-component, respectively. Etotal = (Ex2 + Ey2 + Ez2)1/2 represents the total of the intensity patterns. In the case of the polarization order m = 0 (shown in the top row), the electric field intensity distribution is the same as linear polarization [43]. When the polarization order m = 1 in the second row, the electric field intensity distribution is the same as radial polarization. Furthermore, it is clear to see that several petals are arranged symmetrically on the ring when m = −1 in the third row. Theoretically, the petal-like patterns with the number of 2(|m| + 1) will appear for m < 0. And when m > 1, the total intensity with the arrangement of hot spots with the number of 2(m − 1) in the bottom row. This kind of intensity structure has a potential application value in multiple-particle trapping.
Figure 4 shows Hx, Hy, and Hz that represent the x-component, the y-component, and the z-component, respectively. Htotal = (Hx2 + Hy2 + Hz2)1/2 represents the total magnetic field of the high-order vector beams. Like the electric field intensity distribution, it is observed that several hot spots are arranged symmetrically for m < 0 or m > 1 on the ring. By comparing Figure 3 and Figure 4, it is demonstrated that the electric and magnetic components are perpendicular to each other. Differently, when m = 1, the distribution of electric field intensity shows the solid core, while the distribution of magnetic field intensity has a hollow center.
Figure 5 depicts the three-dimensional intensity distribution of the electric and magnetic fields under tight focusing with different polarization orders. The red part represents the electric field, while the blue part represents the magnetic field. From Figure 3, Figure 4 and Figure 5, it is evident that we can obtain the intensity distribution of the electric and magnetic fields by employing the expressions derived above.
Figure 6 illustrates the intensity contour of the electric and magnetic distributions with different polarization orders when the initial phase ϕ0 is π/4 and π/2, respectively. The arrows represent the polarization distributions. In the top row, where the polarization order m = 0 is depicted, the electric field intensity distribution represents linear polarization. When the polarization order m = 1 and the initial phase ϕ0 = π/2, as shown in the right hand of the second row, the electric and magnetic vectors exhibit the azimuthal and radial distributions, respectively. For m = 0, −1, and 2, the intensity contour varies with the initial phase change.

3.2. Optical Forces on Rayleigh Particle in Tight Focusing of CV Beams

In this section, numerical results are presented to demonstrate the properties of optical forces produced by highly focused partially cylindrical vector (CV) beams exerted on a Rayleigh particle. When a particle is immersed near the focus of light, a compact bright focal spot is expected to capture dielectric particles with a refractive index higher than the ambient index, while a dark core focal spot is expected to capture dielectric particles with a refractive index lower than the ambient index [44]. For numerical calculations, we assume an input power P = 100 mW, a wavelength λ = 840 nm, and an objective lens with a numerical aperture NA = 0.9. The captured gold nanoparticles are assumed to have a radius of 30 nm. The relative permittivity of gold (Au) at a wavelength of 840 nm is ε = −25.096 + 2.045i [45].
Figure 7 illustrates the gradient force distributions by tightly focused CV beams in the transverse and longitudinal planes, respectively. Here, the forces on the particles are linear [46]. The transverse force in the figure represents the vector superposition of the force in the x and y directions, defined as (Fx2 + Fy2)1/2. On the optical axis, due to the symmetry of light, the gradient force is axially symmetric in the focal plane, and from the color bar of the different images in Figure 7, it is observed that the maximum transverse gradient force decreases with the increase in m [47]. Since the intensity region of the electric field enlarges with the increase in m, the average amplitude decreases with the increase in m at a given input power. In Figure 7, arrows denote the direction and magnitude of the transverse gradient force and longitudinal gradient force. The gradient force results in force balance, allowing the Rayleigh particle to be bound.
Figure 8 illustrates the longitudinal scattering force in the focal plane for different orders m of cylindrical vector (CV) beams. Clearly, the scattering force is primarily generated by the energy flow along the propagation direction. In this study, we do not consider the transverse components of the scattering force in the focal plane. Particularly, it is observed that when m = 1, the annular illumination with radial polarization provides a higher gradient force as well as lower scattering forces. Along the optical axis, the scattering force is nearly zero, eliminating a major potential cause of trap destabilization. Furthermore, when m = 2, a force opposite to the light propagation direction is formed near the optical axis at the focus position, causing the light not to push the particles towards the light propagation direction. In this case, the effect of light on the particles resembles that of ‘tractor beams’ [48]. The use of light fields with distinct intensity minima where particles can be retained offers an elegant solution. Essentially, any doughnut beam with a phase or polarization singularity and zero light intensity along its axis satisfies the above criterion, facilitating the formation of a stronger pulling force.
Figure 9 depicts the longitudinal curl-spin force for different orders m of cylindrical vector beams (CVBs). The curl-spin force arises from the gradient in the spatial distribution of beam polarization, resulting from the non-uniform spatial distribution of spin density in the focal region. When the incident light is linearly polarized light, the curl-spin force is zero [49]. However, in a tightly focused system, the curl-spin force is influenced by the polarization of the incident light. As shown in Figure 9, when the polarization of the CVBs changes, the direction of the curl-spin force also changes accordingly. Specifically, when m = −1, 1, 2, the direction of the curl-spin force is opposite to that without polarization, both producing forces toward the center of the light field.
It is evident from Figure 7, Figure 8 and Figure 9 that the transverse gradient force significantly outweighs other forces. Consequently, the positions and number of optical trapping sites can be determined by the transverse gradient force. The gradient force is influenced by the intensity distribution, indicating that by altering m and thus modifying the distribution of the light field, multiple trapping sites could be achieved.

4. Conclusions

In conclusion, we have explored the generation, focusing, and optical force of CV beams. Utilizing the Richards–Wolf vector diffraction integral theory, we derived the integral expressions for the electric field intensity and the magnetic field intensity. Through numerical simulations, we investigated the intensity distribution of the light field and the magnetic field for different orders, providing insights into the 3D intensity distribution under tight focusing. Our findings reveal that by adjusting the order, the spatial distribution and intensity of radial and azimuthal polarization components can be altered, leading to the generation of diverse distributions of electric field components. Furthermore, based on the distribution of the light field and the magnetic field, we characterized the gradient force and explored the scattering force, considering the effects of polarization. We observed that vector beams with arbitrary polarization distributions can generate stable optical trap force potentials, and polarization can influence the direction and magnitude of the force. The non-uniform spatial polarization state of high-order vector beams holds a significant application potential in various fields such as optical control, micro-machining, remote sensing, and space optical communication.

Author Contributions

Investigation, methodology, writing—original draft, formal analysis, and software, Y.Z. and L.Z. (Liqiang Zhou); funding acquisition, project administration, resources, supervision, and validation, L.Z. (Linwei Zhu) and Q.S.; conceptualization, methodology, software, and data curation, X.J. and L.Z. (Linwei Zhu). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (62174073, 61675093), the Program of Science and Technology Development of Yantai (2020XDRH095), and the Project of Taishan industry leading talents (tscx202211051).

Data Availability Statement

All data that support the findings of this study are included within the article.

Conflicts of Interest

Author Qiang Shi was employed by the company Moji-Nano Technology Co., Ltd., Yantai, China. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships.

References

  1. Zhan, Q. Cylindrical vector beams: From mathematical concepts to applications. Adv. Opt. Photonics 2009, 1, 1–57. [Google Scholar] [CrossRef]
  2. Cheng, W.; Haus, J.W.; Zhan, Q. Propagation of vector vortex beams through a turbulent atmosphere. Opt. Express 2009, 17, 17829–17836. [Google Scholar] [CrossRef] [PubMed]
  3. Wang, X.; Li, Y.; Chen, J.; Guo, C.; Ding, J.; Wang, H. A new type of vector fields with hybrid states of polarization. Opt. Express 2010, 18, 10786–10795. [Google Scholar] [CrossRef] [PubMed]
  4. Ren, Z.; Cheng, Z.; Wang, X.; Ding, J.; Wang, H. Polarization interferometric prism: A versatile tool for generation of vector fields, measurement of topological charges, and implementation of a spin–orbit controlled-Not gate. Appl. Phys. Lett. 2021, 118, 011105. [Google Scholar] [CrossRef]
  5. Wang, Q.; Tu, C.; Li, Y.; Wang, H. Polarization singularities: Progress, fundamental physics, and prospects. APL Photonics 2021, 6, 040901. [Google Scholar] [CrossRef]
  6. Zhu, L.; Cao, Y.; Chen, Q.; Ouyang, X.; Xu, Y.; Hu, Z.; Qiu, J.; Li, X. Near-perfect fidelity polarization-encoded multilayer optical data storage based on aligned gold nanorods. Opto-Electron. Adv. 2021, 4, 210002. [Google Scholar] [CrossRef]
  7. Han, W.; Yang, Y.; Cheng, W.; Zhan, Q. Vectorial optical field generator for the creation of arbitrarily complex fields. Opt. Express 2013, 21, 20692–20706. [Google Scholar] [CrossRef] [PubMed]
  8. Wang, X.; Ding, J.; Ni, W.; Guo, C.; Wang, H. Generation of arbitrary vector beams with a spatial light modulator and a common path interferometric arrangement. Opt. Lett. 2007, 32, 3549–3551. [Google Scholar] [CrossRef] [PubMed]
  9. Rui, G.; Zhan, Q. Trapping of resonant metallic nanoparticles with engineered vectorial optical field. Nanophotonics 2014, 3, 351–361. [Google Scholar] [CrossRef]
  10. Wang, X.; Chen, J.; Li, Y.; Ding, J.; Guo, C.; Wang, H. Optical orbital angular momentum from the curl of polarization. Phys. Rev. Lett. 2010, 105, 253602. [Google Scholar] [CrossRef]
  11. Zhang, X.; Chen, R.; Wang, A. Focusing properties of cylindrical vector vortex beams. Opt. Commun. 2018, 414, 10–15. [Google Scholar] [CrossRef]
  12. Hnatovsky, C.; Shvedov, V.; Krolikowski, W.; Rode, A. Revealing local field structure of focused ultrashort pulses. Phys. Rev. Lett. 2011, 106, 123901. [Google Scholar] [CrossRef]
  13. Li, M.; Yan, S.; Yao, B.; Lei, M.; Yang, Y.; Min, J.; Dan, D. Intrinsic optical torque of cylindrical vector beams on Rayleigh absorptive spherical particles. J. Opt. Soc. Am. A 2014, 31, 1710–1715. [Google Scholar] [CrossRef] [PubMed]
  14. Gong, L.; Gu, B.; Rui, G.; Cui, Y.; Zhu, Z.; Zhan, Q. Optical forces of focused femtosecond laser pulses on nonlinear optical Rayleigh particles. Photonics Res. 2018, 6, 138–143. [Google Scholar] [CrossRef]
  15. Lehmuskero, A.; Johansson, P.; Rubinsztein-Dunlop, H.; Tong, L.; Kall, M. Laser trapping of colloidal metal nanoparticles. ACS Nano 2015, 9, 3453–3469. [Google Scholar] [CrossRef] [PubMed]
  16. Gao, D.; Ding, W.; Nieto-Vesperinas, M.; Ding, X.; Rahman, M.; Zhang, T.; Lim, C.; Qiu, C.W. Optical manipulation from the microscale to the nanoscale: Fundamentals, advances and prospects. Light Sci. Appl. 2017, 6, e17039. [Google Scholar] [CrossRef] [PubMed]
  17. Grier, D.G. A revolution in optical manipulation. Nature 2003, 424, 810–816. [Google Scholar] [CrossRef]
  18. Ashkin, A.; Dziedzic, J.M.; Bjorkholm, J.E.; Chu, S. Observation of a single-beam gradient force optical trap for dielectric particles. Opt. Lett. 1986, 11, 288–290. [Google Scholar] [CrossRef] [PubMed]
  19. Pan, Y.; Ren, Z.; Kong, L.; Tu, C.; Li, Y.; Wang, H. Theoretical analysis based on mirror symmetry for tightly focused vector optical fields. Opt. Express 2020, 28, 23416–23432. [Google Scholar] [CrossRef] [PubMed]
  20. Bliokh, K.Y.; Rodríguez-Fortuño, F.J.; Nori, F.; Zayats, A.V. Spin–orbit interactions of light. Nat. Photonics 2015, 9, 796–808. [Google Scholar] [CrossRef]
  21. Fu, Y.; Wang, Y.; Zhang, Y.; He, Y.; Min, C.; Yuan, X. Detecting cylindrical vector beams with an on-chip plasmonic spin-Hall metalens. Opt. Express 2022, 30, 10758–10769. [Google Scholar] [CrossRef] [PubMed]
  22. Li, H.; Ma, C.; Wang, J.; Tang, M.; Li, X. Spin-orbit Hall effect in the tight focusing of a radially polarized vortex beam. Opt. Express 2021, 29, 39419–39427. [Google Scholar] [CrossRef] [PubMed]
  23. Chen, S.; Xie, Z.; Ye, H.; Wang, X.; Guo, Z.; He, Y.; Li, Y.; Yuan, X.; Fan, D. Cylindrical vector beam multiplexer/demultiplexer using off-axis polarization control. Light Sci. Appl. 2021, 10, 222. [Google Scholar] [CrossRef] [PubMed]
  24. Bautista, G.; Huttunen, M.J.; Makitalo, J.; Kontio, J.M.; Simonen, J.; Kauranen, M. Second-harmonic generation imaging of metal nano-objects with cylindrical vector beams. Nano Lett. 2012, 12, 3207–3212. [Google Scholar] [CrossRef] [PubMed]
  25. Zhan, Q. Trapping metallic Rayleigh particles with radial polarization. Opt. Express 2004, 12, 3377–3382. [Google Scholar] [CrossRef] [PubMed]
  26. Li, M.; Yan, S.; Liang, Y.; Zhang, P.; Yao, B. Transverse spinning of particles in highly focused vector vortex beams. Phys. Rev. A 2017, 95, 053802. [Google Scholar] [CrossRef]
  27. Bekshaev, A.Y.; Angelsky, O.V.; Hanson, S.G.; Zenkova, C.Y. Scattering of inhomogeneous circularly polarized optical field and mechanical manifestation of the internal energy flows. Phys. Rev. A 2012, 86, 023847. [Google Scholar] [CrossRef]
  28. Angelsky, O.V.; Bekshaev, A.Y.; Maksimyak, P.P.; Maksimyak, A.P.; Mokhun, I.I.; Hanson, S.G.; Zenkova, C.Y.; Tyurin, A.V. Circular motion of particles suspended in a Gaussian beam with circular polarization validates the spin part of the internal energy flow. Opt. Express 2012, 20, 11351–11356. [Google Scholar] [CrossRef] [PubMed]
  29. Jabbour, T.G.; Kuebler, S.M. Vector diffraction analysis of high numerical aperture focused beams modified by two- and three-zone annular multi-phase plates. Opt. Express 2006, 14, 1033–1043. [Google Scholar] [CrossRef] [PubMed]
  30. Canales, V.F.; Oti, J.E.; Cagigal, M.P. Three-dimensional control of the focal light intensity distribution by analytically designed phase masks. Opt. Commun. 2005, 247, 11–18. [Google Scholar] [CrossRef]
  31. Rui, G.; Wang, X.; Gu, B.; Zhan, Q.; Cui, Y. Manipulation metallic nanoparticle at resonant wavelength using engineered azimuthally polarized optical field. Opt. Express 2016, 24, 7212–7223. [Google Scholar] [CrossRef] [PubMed]
  32. Chen, J.; Ng, J.; Lin, Z.; Chan, C.T. Optical pulling force. Nat. Photonics 2011, 5, 531–534. [Google Scholar] [CrossRef]
  33. Richards, B.; Wolf, E. Electromagnetic diffraction in optical systems II. Structure of the image field in an aplanatic system. Proc. R. Soc. Ser. A 1959, 253, 358–379. [Google Scholar] [CrossRef]
  34. Hnatovsky, C.; Shvedov, V.G.; Shostka, N.; Rode, A.V.; Krolikowski, W. Polarization-dependent ablation of silicon using tightly focused femtosecond laser vortex pulses. Opt. Lett. 2012, 37, 226–228. [Google Scholar] [CrossRef] [PubMed]
  35. Youngworth, K.; Brown, T. Focusing of High Numerical Aperture Cylindrical-Vector Beams. Opt. Express 2000, 7, 77–87. [Google Scholar] [CrossRef] [PubMed]
  36. Gordon, J.P. Radiation Forces and Momenta in Dielectric Media. Phys. Rev. A 1973, 8, 14–21. [Google Scholar] [CrossRef]
  37. Nieto-Vesperinas, M.; Xu, X. Reactive helicity and reactive power in nanoscale optics: Evanescent waves. Kerker conditions. Optical theorems and reactive dichroism. Phys. Rev. Res. 2021, 3, 043080. [Google Scholar] [CrossRef]
  38. Li, M.; Yan, S.; Yao, B.; Liang, Y.; Han, G.; Zhang, P. Optical trapping force and torque on spheroidal Rayleigh particles with arbitrary spatial orientations. J. Opt. Soc. Am. A 2016, 33, 1341–1347. [Google Scholar] [CrossRef] [PubMed]
  39. Li, M.; Yan, S.; Yao, B.; Lei, M.; Yang, Y.; Min, J.; Dan, D. Trapping of Rayleigh spheroidal particles by highly focused radially polarized beams. J. Opt. Soc. Am. B 2015, 32, 468–472. [Google Scholar] [CrossRef]
  40. Albaladejo, S.; Marques, M.I.; Laroche, M.; Saenz, J.J. Scattering forces from the curl of the spin angular momentum of a light field. Phys. Rev. Lett. 2009, 102, 113602. [Google Scholar] [CrossRef]
  41. Zhang, Y.; Shen, J.; Min, C.; Jin, Y.; Jiang, Y.; Liu, J.; Zhu, S.; Sheng, Y.; Zayats, A.V.; Yuan, X. Nonlinearity-Induced Multiplexed Optical Trapping and Manipulation with Femtosecond Vector Beams. Nano Lett. 2018, 18, 5538–5543. [Google Scholar] [CrossRef]
  42. Marago, O.M.; Jones, P.H.; Gucciardi, P.G.; Volpe, G.; Ferrari, A.C. Optical trapping and manipulation of nanostructures. Nat. Nanotechnol. 2013, 8, 807–819. [Google Scholar] [CrossRef] [PubMed]
  43. Jeong, T.M.; Bulanov, S.; Weber, S.; Korn, G. Analysis on the longitudinal field strength formed by tightly-focused radially-polarized femtosecond petawatt laser pulse. Opt. Express 2018, 26, 33091–33107. [Google Scholar] [CrossRef] [PubMed]
  44. Shu, J.; Chen, Z.; Pu, J. Radiation forces on a Rayleigh particle by highly focused partially coherent and radially polarized vortex beams. J. Opt. Soc. Am. A Opt. Image Sci. Vis. 2013, 30, 916–922. [Google Scholar] [CrossRef] [PubMed]
  45. Johnson, P.B.; Christy, R.W. Optical Constants of the Noble Metals. Phys. Rev. B 1972, 6, 4370–4379. [Google Scholar] [CrossRef]
  46. du Preez-Wilkinson, N.; Stilgoe, A.B.; Alzaidi, T.; Rubinsztein-Dunlop, H.; Nieminen, T.A. Forces due to pulsed beams in optical tweezers: Linear effects. Opt. Express 2015, 23, 7190–7208. [Google Scholar] [CrossRef] [PubMed]
  47. Zhang, L.; Qiu, X.; Zeng, L.; Chen, L. Multiple trapping using a focused hybrid vector beam. Chin. Phys. B 2019, 28, 094202. [Google Scholar] [CrossRef]
  48. Liu, H.; Panmai, M.; Peng, Y.; Lan, S. Optical pulling and pushing forces exerted on silicon nanospheres with strong coherent interaction between electric and magnetic resonances. Opt. Express 2017, 25, 12357–12371. [Google Scholar] [CrossRef] [PubMed]
  49. Liberal, I.; Ederra, I.; Gonzalo, R.; Ziolkowski, R.W. Near-field electromagnetic trapping through curl-spin forces. Phys. Rev. A 2013, 87, 063807. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of focusing system geometry.
Figure 1. Schematic diagram of focusing system geometry.
Nanomaterials 14 00691 g001
Figure 2. Polarization distributions of incident vector beam with different orders: (a) m = 0, ϕ0 = 0, (b) m = 1, ϕ0 = 0, (c) m = −1, ϕ0 = 0, and (d) m = 2, ϕ0 = 0. Arrows indicate the direction of polarization at the arrow location.
Figure 2. Polarization distributions of incident vector beam with different orders: (a) m = 0, ϕ0 = 0, (b) m = 1, ϕ0 = 0, (c) m = −1, ϕ0 = 0, and (d) m = 2, ϕ0 = 0. Arrows indicate the direction of polarization at the arrow location.
Nanomaterials 14 00691 g002
Figure 3. The intensity distribution on the focal plane of vector beams Ex, Ey, Ez and Etotal with different orders: (ad) m = 0, (eh) m = 1, (il) m = −1, (mp) m = 2.
Figure 3. The intensity distribution on the focal plane of vector beams Ex, Ey, Ez and Etotal with different orders: (ad) m = 0, (eh) m = 1, (il) m = −1, (mp) m = 2.
Nanomaterials 14 00691 g003
Figure 4. The intensity distribution on the focal plane of vector beams Hx, Hy, Hz and Htotal with different orders: (ad) m = 0, (eh) m = 1, (il) m = −1, (mp) m = 2.
Figure 4. The intensity distribution on the focal plane of vector beams Hx, Hy, Hz and Htotal with different orders: (ad) m = 0, (eh) m = 1, (il) m = −1, (mp) m = 2.
Nanomaterials 14 00691 g004
Figure 5. The iso-surface intensity distribution of the electric and magnetic fields in the tightly focused region with E = 0.5Emax and H = 0.5Hmax.
Figure 5. The iso-surface intensity distribution of the electric and magnetic fields in the tightly focused region with E = 0.5Emax and H = 0.5Hmax.
Nanomaterials 14 00691 g005
Figure 6. The intensity contour of the electric and magnetic fields in the focal plane for the vector beams: (a,b) m = 0, ϕ0 = π/4, (c,d) m = 0, ϕ0 = π/2, (e,f) m = 1, ϕ0 = π/4, (g,h) m = 1, ϕ0 = π/2, (i,j) m = −1, ϕ0 = π/4, (k,l) m = −1, ϕ0 =π/2, (m,n) m = 2, ϕ0 = π/4, (o,p) m = 2, ϕ0 = π/2. The short lines indicate the polarization distributions.
Figure 6. The intensity contour of the electric and magnetic fields in the focal plane for the vector beams: (a,b) m = 0, ϕ0 = π/4, (c,d) m = 0, ϕ0 = π/2, (e,f) m = 1, ϕ0 = π/4, (g,h) m = 1, ϕ0 = π/2, (i,j) m = −1, ϕ0 = π/4, (k,l) m = −1, ϕ0 =π/2, (m,n) m = 2, ϕ0 = π/4, (o,p) m = 2, ϕ0 = π/2. The short lines indicate the polarization distributions.
Nanomaterials 14 00691 g006
Figure 7. The gradient force distributions produced by highly focused CV beams with m equal to 0, 1, −1, and 2. (ad) represent the transverse gradient force distributions. (eh) represent the longitudinal gradient force distributions. Arrows denote the direction and magnitude of the gradient force.
Figure 7. The gradient force distributions produced by highly focused CV beams with m equal to 0, 1, −1, and 2. (ad) represent the transverse gradient force distributions. (eh) represent the longitudinal gradient force distributions. Arrows denote the direction and magnitude of the gradient force.
Nanomaterials 14 00691 g007
Figure 8. The longitudinal scattering force in the focal plane for different orders m of CV beam: (a) m = 0, (b) m = 1, (c) m = −1, (d) m = 2. Arrows denote the direction and magnitude of the scattering force.
Figure 8. The longitudinal scattering force in the focal plane for different orders m of CV beam: (a) m = 0, (b) m = 1, (c) m = −1, (d) m = 2. Arrows denote the direction and magnitude of the scattering force.
Nanomaterials 14 00691 g008
Figure 9. The longitudinal curl-spin force for different orders m of CV beam: (a) m = 0, (b) m = 1, (c) m = −1, (d) m = 2. Arrows denote the direction and magnitude of the curl-spin force.
Figure 9. The longitudinal curl-spin force for different orders m of CV beam: (a) m = 0, (b) m = 1, (c) m = −1, (d) m = 2. Arrows denote the direction and magnitude of the curl-spin force.
Nanomaterials 14 00691 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, Y.; Zhou, L.; Jiang, X.; Zhu, L.; Shi, Q. Optical Force Effects of Rayleigh Particles by Cylindrical Vector Beams. Nanomaterials 2024, 14, 691. https://doi.org/10.3390/nano14080691

AMA Style

Zhao Y, Zhou L, Jiang X, Zhu L, Shi Q. Optical Force Effects of Rayleigh Particles by Cylindrical Vector Beams. Nanomaterials. 2024; 14(8):691. https://doi.org/10.3390/nano14080691

Chicago/Turabian Style

Zhao, Yuting, Liqiang Zhou, Xiaotong Jiang, Linwei Zhu, and Qiang Shi. 2024. "Optical Force Effects of Rayleigh Particles by Cylindrical Vector Beams" Nanomaterials 14, no. 8: 691. https://doi.org/10.3390/nano14080691

APA Style

Zhao, Y., Zhou, L., Jiang, X., Zhu, L., & Shi, Q. (2024). Optical Force Effects of Rayleigh Particles by Cylindrical Vector Beams. Nanomaterials, 14(8), 691. https://doi.org/10.3390/nano14080691

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