**3. Radiation force induced by CVB**

#### **3.1. Optical force on Rayleigh particles**

In the Rayleigh regime, particles can be considered as infinitesimal induced dipoles which interact with incident beam. Here we assume that the particle is a microsphere. RPF will be decomposed into scattering force and gradient force.

The oscillating dipole, which is induced by time-harmonic fields, can be considered as an antenna. The antenna will radiate energy. The difference between energy removed from incident beam and energy radiated by the antenna accounts for the change of momentum flux, and hence rusults in a scattering force. The scattering force can be expressed as

$$\vec{F}\_{scat} = \hat{\mathbf{e}}\_z \mathbf{C}\_{\rho\nu} \mathbf{n}\_1^2 \mathbf{e}\_0 \left| E \right|^2 \tag{4}$$

with

$$C\_{pr} = C\_{scat} = \frac{8}{3} \pi (ka)^4 a^2 \left(\frac{m^2 - 1}{m^2 + 2}\right)^2 \tag{5}$$

and

$$m = n\_2 \, / \, n\_1 \tag{6}$$

where *n1* is the refractive index of surrounding media, and *n2* is the refractive index of the particle. *a* is the radius of microsphere. *ε<sup>0</sup>* is the dielectric constant in the vacumm. Note that the scattering force always points in the direction of incident beam.

When a particle is illuminated by a non-uniform electric field, it will experience a gradient force.

$$\vec{F}\_{gud} = \pi n\_{\text{i}}^2 \varepsilon\_0 a^3 \left(\frac{m^2 - 1}{m^2 + 2}\right) \nabla \left|E\right|^2\tag{7}$$

For a time-harmonic field, the gradient force can also be expressed in terms of the intensity *I* of incident beam:

$$\vec{F}\_{gud} = \frac{2\pi n\_i a^3}{c} \left(\frac{m^2 - 1}{m^2 + 2}\right) \nabla I \tag{8}$$

where *c* is the speed of light in the vacumm. It is obvious that the gradient force depends on the gradient of the intensity. By sbustituting Eqs. (2) and (3) into Eqs. (4) and (7), we can obtain the scattering and gradient force of vector Bessel-Gaussian beams exerted on a microsphere. For radially polarized Bessel-Gaussian beam, they can be expressed as

$$\vec{F}\_{scat} = \hat{\mathbf{e}}\_z C\_{pr} n\_1^2 \varepsilon\_0 E\_0^2 \frac{\rho^2}{\mathcal{W}\_0^2} \left( e^{-\frac{\rho^2}{\mathcal{W}\_0^2}} \right)^2 \tag{9}$$

$$\vec{F}\_{grad} = \pi n\_1^2 \varepsilon\_0 a^3 \left(\frac{m^2 - 1}{m^2 + 2}\right) \left[\frac{2E\_0^2 \rho}{\nu\_0^2} \left(e^{-\frac{\rho^2}{w\_0^2}}\right)^2 - \frac{4E\_0^2 \rho^3}{\nu\_0^4} \left(e^{-\frac{\rho^2}{w\_0^2}}\right)^2\right] \hat{e}\_\rho \tag{10}$$

From Eq. (10), we can find that the gradient force has only ρ component. This is because *|E|2* is only dependent on ρ. Here we give only the force for radially polarized beam incidence, and that for azimuthally polarized beam incidence can be derived in the same way.

#### **3.2. Radiation force exerted on Mie particles**

with

and

force.

of incident beam:

<sup>2</sup> <sup>2</sup>

*<sup>k</sup> <sup>m</sup> CC a* (5)

2 1 *mnn* = / (6)

2

*m*

4 2

where *n1* is the refractive index of surrounding media, and *n2* is the refractive index of the particle. *a* is the radius of microsphere. *ε<sup>0</sup>* is the dielectric constant in the vacumm. Note that

When a particle is illuminated by a non-uniform electric field, it will experience a gradient

<sup>2</sup> <sup>2</sup> 2 3

*m a m*

è ø

For a time-harmonic field, the gradient force can also be expressed in terms of the intensity *I*

3 2

æ ö - <sup>=</sup> ç ÷

*n m c m*

2 2 1

where *c* is the speed of light in the vacumm. It is obvious that the gradient force depends on the gradient of the intensity. By sbustituting Eqs. (2) and (3) into Eqs. (4) and (7), we can obtain the scattering and gradient force of vector Bessel-Gaussian beams exerted on a microsphere.

è ø

2

r

*w*

2 2 2 3

2 4 r

0 0 <sup>1</sup> <sup>ˆ</sup> <sup>2</sup>

è ø èø èø ë û

2 2 100 2 0

r

*m E E a ee mw w*


= ç ÷

e

r *<sup>w</sup>*

2 3 0 0 1 0 22 4

r *w w*

+

2

Ñ

*e nF E e* (9)

 r

r

2 2 2 2 0 0 2 2

 r

r

æ ö -

ç ÷ è ø

*F n e* (10)

1

p

For radially polarized Bessel-Gaussian beam, they can be expressed as

ˆ

*scat z pr C*


p e

*grad*

<sup>æ</sup> - <sup>Ñ</sup> +

1 2

*F E* (7)

*<sup>F</sup> <sup>a</sup> <sup>I</sup>* (8)

1 0 2

<sup>ö</sup> <sup>=</sup> ç ÷

p ( ) æ ö - = = ç ÷ è ø <sup>+</sup> *pr scat <sup>a</sup>*

the scattering force always points in the direction of incident beam.

202 Advances in Optical Fiber Technology: Fundamental Optical Phenomena and Applications

p e

*grad n*

r

r *grad* 8 1 3 2

> Many practical particles manipulated with optical tweezers, such as bioloical cells, are Mie particles, whose size is in the order of the wavelength of trapping beam. To calculate the radiation force exerted on such particles, a rigorous electromagnetic theory based on the Maxwell equations must be considered. Generalized Lorenz-Mie Thoery (GLMT) developed by Gouesbet et al. can solve the interaction between homogeneous spheres and focused beams with any shape, and has been entended to solve the scattering of shaped beam by multilayered spheres, homogeneous and multilayered cylinders, and homogeneous and multilayered spheroids. GLMT has been applied to the rigorous calculation of radiation pressure and optical torque. In GLMT, the incident beam is described by a set of beam shape coefficients(BSCs), which can be evaluated by integral localized approximation (ILA) [34].

> This section is devoted to the GLMT for radiation force exerted on a sphere illuminated by a vetor Bessel-Gaussian beam. The general theory for radiation force based on electromagnetic scattering theory is followed by BSCs for CVB. To clarify the physical interpretation of various features of RPF that are implicit in the GLMT, Debye Series Expansion (DSE) is introduced.

#### *3.2.1. Generization Lorenz-Mie theory*

Consider a sphere with radius *a* and refractive index *m1* illuminated by a CVB of wavelength λ in the surrounding media. The center of the sphere is located at OP, origin of the Cartesian coordinate system *OP-xyz*. The beam center is at *OG*, origin of coordinate system *OG-uvw*, with *u* axis parallel to *x* and similarly for the others. The coordinates of *OG* in the system *OP-xyz* is (*x0*, *y0*, *z0*). The refractive index of surrounding media is *m2*. The other parameters are defined in Fig. 2.

When a sphere is illuminated by focused beam, the RPF is proportional to the net momentum removed from the incident beam, and can be expressed in terms of the surface integration of Maxwell stress tensor

$$\mathbf{x} < \mathbf{F} > \dots < \widehat{\oint\_{S} \mathbf{n} \cdot \tilde{T} dS} > \tag{11}$$

where < > represents a time average, *n* ^ the outward normal unit vector, and *S* a surface enclosing the particle. The Maxwell stress tensor *T* ↔ is given by

$$<\vec{T}> = \frac{1}{4\pi} \left( \mathcal{s} \mathbf{E} \mathbf{E} + \mathbf{H} \mathbf{H} - \frac{1}{2} \left( \mathcal{s} \boldsymbol{E}^2 + \boldsymbol{H}^2 \right) \vec{I} \right) \tag{12}$$

**Figure 2.** Coordinate systems in GLMT. *OP-xyz* is attached to the sphere, and *OG-uvw* to the incident beam.

Where the electromagnetic fields *E* and *H* are the total fields, namely the sum of the incident and scattered fields, given by

$$\mathbf{E} = \mathbf{E}\_i + \mathbf{E}\_s,\\\mathbf{H} = \mathbf{H}\_i + \mathbf{H}\_s \tag{13}$$

*Ei* and *H<sup>i</sup>* are the incident electromagnetic wave, and can be expaned as:

$$E\_r^i = \frac{E\_0}{k\_0^2 r^2} \sum\_{n=1}^n \sum\_{m=-n}^{+\pi} c\_n^{m\nu} \mathbf{g}\_{n, \Im M}^n n(n+1) \wp\_n(k\_0 r) P\_n^{|m|}(\cos \theta) e^{(im\psi)} \tag{14}$$

$$E\_{\theta}^{i} = \frac{E\_{0}}{k\_{0}r} \sum\_{n=1}^{n} \sum\_{m=-n}^{+\pi} c\_{n}^{\prime\prime\prime} \left[ \mathbf{g}\_{n,\Pi t}^{\prime\prime} \boldsymbol{\nu}\_{n}^{\prime}(k\_{0}r) \boldsymbol{\tau}\_{n}^{|m|}(\cos\theta) + m \mathbf{g}\_{n,\Pi t}^{\prime\prime} \boldsymbol{\nu}\_{n}(k\_{0}r) \boldsymbol{\pi}\_{n}^{|m|}(\cos\theta) \right] e^{(im\boldsymbol{\rho})} \tag{15}$$

$$E\_{\rho}^{i} = i \frac{E\_{0}}{k\_{0}r} \sum\_{n=1}^{n} \sum\_{m=-n}^{+\pi} c\_{n}^{\rho n} \left[ \text{mg}\_{n, \text{TM}}^{n} \wp\_{n}^{\prime}(k\_{0}r) \pi\_{n}^{|m|}(\cos \theta) + \text{g}\_{n, \text{TE}}^{\pi} \wp\_{n}(k\_{0}r) \pi\_{n}^{|m|}(\cos \theta) \right] e^{(im\phi)} \tag{16}$$

$$H\_r^i = \frac{H\_0}{k\_0^2 r^2} \sum\_{n=1}^n \sum\_{m=-n}^{+n} c\_n^{mr} \mathbf{g}\_{n, \text{TE}}^n n(n+1) \boldsymbol{\nu}\_n(k\_0 r) P\_n^{|m|}(\cos \theta) e^{(m\boldsymbol{\nu})} \tag{17}$$

$$H\_{\vartheta}^{i} = -\frac{H\_{0}}{k\_{0}r} \sum\_{v=1}^{n} \sum\_{m=-a}^{+\pi} \mathcal{c}\_{\pi}^{mv} \left[ \mathsf{m} \mathsf{g}\_{n, \text{TM}}^{m} \psi\_{n}(k\_{0}r) \pi\_{n}^{|m|}(\cos \theta) - \mathsf{g}\_{n, \text{KL}}^{m} \psi\_{n}^{\prime}(k\_{0}r) \pi\_{n}^{|m|}(\cos \theta) \right] e^{(im\phi)} \tag{18}$$

$$H\_{\varphi}^{\prime} = -\frac{iH\_0}{k\_0 r} \sum\_{n=1}^{n} \sum\_{m=-n}^{+n} \mathcal{c}\_{\pi}^{\prime m} \left[ \mathbf{g}\_{\pi, \text{Dir}}^{\prime \prime} \boldsymbol{\nu}\_{\pi}(k\_0 r) \boldsymbol{\pi}\_{\pi}^{|m|}(\cos \theta) - m \mathbf{g}\_{\pi, \text{IP}}^{\prime \prime} \boldsymbol{\nu}\_{\pi}^{\prime}(k\_0 r) \boldsymbol{\pi}\_{\pi}^{|m|}(\cos \theta) \right] \mathbf{e}^{(m\varphi)} \tag{19}$$

with

Where the electromagnetic fields *E* and *H* are the total fields, namely the sum of the incident

*v*

0 || ( )

0 | | | | ( ) , 0 , 0

0 | | | | ( ) , 0 , 0

0 || ( )

y

 q

 q

= + å å*<sup>n</sup> <sup>i</sup> pw m <sup>m</sup> im*

= + é ù ¢ å å ë û

= + é ù ¢ å å ë û

*i pw m m m m im n n TM n n n TE n n*

= + å å*<sup>n</sup> <sup>i</sup> pw m <sup>m</sup> im*

*i pw m m m m im n n TM n n n TE n n*

y

( 1) ( ) (cos )

( ) (cos ) ( ) (cos )

*<sup>E</sup> c g k r mg k r e k r* (15)

( ) (cos ) ( ) (cos )

*<sup>E</sup> E i c mg k r g kr e k r* (16)

( 1) ( ) (cos )

 yt

 yp

are the incident electromagnetic wave, and can be expaned as:

2 2 , 0

*x*

0 00 (, ,) *Oxyz <sup>G</sup>*

**Figure 2.** Coordinate systems in GLMT. *OP-xyz* is attached to the sphere, and *OG-uvw* to the incident beam.

*w*

204 Advances in Optical Fiber Technology: Fundamental Optical Phenomena and Applications

*r n n TM n n*

yt

yp

2 2 , 0

*r n n TE n n*

0 1

0 1

*H*

¥ + = =-

*nmn*

*E*

*u*

0 1

0 1

*E*

¥ + = =-

*nmn*

¥ + = =-

*nmn*

*n*

*n*

q

j ¥ + = =-

*nmn*

=+ = + , **E E EH H H** *is i s* (13)

*y*

*P*

*z*

q

f

*OP*

*r*

j

> q

> > q

j

 q

*<sup>H</sup> c g nn krP e k r* (17)

j

j

 q

*<sup>E</sup> c g nn krP e k r* (14)

and scattered fields, given by

*Ei* and *H<sup>i</sup>*

$$\left(\pi\_{\pi}^{\pi^{\pi}}(\cos \theta) = \frac{dP\_{\pi}^{\pi^{\pi}}(\cos \theta)}{d\theta}\right) \tag{20}$$

$$\pi\_{\pi}^{\pi^{\pi}}(\cos \theta) = m \frac{P\_{\pi}^{\pi^{\pi}}(\cos \theta)}{\sin \theta} \tag{21}$$

$$c\_{\
u}^{\rho\nu} = (-i)^{\nu+l} \frac{2n+l}{n(n+l)} \tag{22}$$

where *Pn <sup>m</sup>*(cos*θ*) represents the associated Legendre polynomials of degree *n* and order *m, ψ*(· ) is the spherical Ricatti-Bessel functions of first kind, and the prime indicates the derivative of the function with respect to its argument. *gn*,*TM <sup>m</sup>* and *gn*,*TE <sup>m</sup>* are so-called BSCs and will be discussed in next subsection.

Similarly, the scattered fields *Es* and *Hs* have the expression :

$$E\_r^s = -\frac{E\_0}{k\_0^2 r^2} \sum\_{n=1}^{n} \sum\_{m=-n}^{+\pi} c\_n^{m\*} A\_n^m n(n+1) \xi\_n(k\_0 r) P\_n^{|n|}(\cos \theta) e^{i(m\psi)} \tag{23}$$

$$E\_{\theta}^{\prime} = -\frac{E\_0}{k\_0 r} \sum\_{n=1}^{n} \sum\_{m=-n}^{+n} c\_n^{\prime m} \left[ A\_n^m \xi\_n^{\prime}(k\_0 r) \pi\_n^{|m|}(\cos \theta) + m B\_n^m \xi\_n(k\_0 r) \pi\_n^{|m|}(\cos \theta) \right] e^{(im\varphi)} \tag{24}$$

$$E\_{\varphi}^{\prime} = -i \frac{E\_0}{k\_0 r} \sum\_{n=1}^{n} \sum\_{m=-n}^{+n} c\_{\pi}^{\prime m} \left[ m A\_{\pi}^{\prime n} \tilde{\xi}\_{n}^{\prime}(k\_0 r) \pi\_{n}^{|m|}(\cos \theta) + B\_{\pi}^{\prime n} \tilde{\xi}\_{n}(k\_0 r) \pi\_{n}^{|m|}(\cos \theta) \right] e^{(m\rho)} \tag{25}$$

$$H'\_r = -\frac{H\_0}{k\_0^2 r^2} \sum\_{n=1}^n \sum\_{m=-n}^{+n} c\_n^{\prime\prime\prime} B\_n^m n(n+1) \xi\_n(k\_0 r) P\_n^{|m|}(\cos \theta) e^{i(m\psi)} \tag{26}$$

$$H'\_{\vartheta} = \frac{H\_0}{k\_0 r} \sum\_{n=1}^{n} \sum\_{m=-n}^{+n} c\_n^{(m)} \left[ m A\_n'' \xi\_n(k\_0 r) \pi\_n^{|m|}(\cos \theta) - B\_n'' \xi\_n'(k\_0 r) \pi\_n^{|m|}(\cos \theta) \right] e^{(im\phi)} \tag{27}$$

$$H'\_{\varphi} = \frac{iH\_0}{k\_0 r} \sum\_{n=1}^{n} \sum\_{m=-n}^{+n} c\_n^{m\*} \left[ A\_n^m \xi\_n(k\_0 r) \pi\_n^{|m|}(\cos \theta) - m B\_n^m \xi\_{nn}^{'}(k\_0 r) \pi\_n^{|m|}(\cos \theta) \right] e^{(mp)} \tag{28}$$

Where *ξn*(*k*0*r*) is Ricatti-Hankel functions, and scattering coefficients *An <sup>m</sup>* and *Bn <sup>m</sup>* can be expressed by traditional Mie scattering coefficients *an*, *bn* and BSCs *gn*,*TM <sup>m</sup>* , *gn*,*TE <sup>m</sup>* :

$$A\_n^m = a\_n \mathbf{g}\_{n,TM}^m, \qquad B\_n^m = b\_n \mathbf{g}\_{n,\mathcal{IE}}^m \tag{29}$$

with

$$a\_{\pi} = \frac{-m\_{\text{i}}\boldsymbol{\nu}\_{\text{n}}^{\prime}(\mathbf{x})\boldsymbol{\nu}\_{\text{n}}^{\prime}(\mathbf{y}) + m\_{2}\boldsymbol{\nu}^{\prime}(\mathbf{x})\boldsymbol{\nu}\_{\text{n}}^{\prime}(\mathbf{y})}{-m\_{1}\boldsymbol{\xi}\_{\text{n}}^{\prime(\text{i})^{\prime}}(\mathbf{x})\boldsymbol{\nu}\_{\text{n}}^{\prime}(\mathbf{y}) + m\_{2}\boldsymbol{\xi}\_{\text{n}}^{\prime(\text{i})}(\mathbf{x})\boldsymbol{\nu}\_{\text{n}}^{\prime}(\mathbf{y})} \tag{30}$$

$$b\_{\boldsymbol{n}} = \frac{-m\_{2}\boldsymbol{\upmu}\_{\boldsymbol{n}}^{\boldsymbol{\upprime}}(\mathbf{x})\boldsymbol{\upmu}\_{\boldsymbol{n}}(\mathbf{y}) + m\_{1}\boldsymbol{\upmu}(\mathbf{x})\boldsymbol{\upmu}\_{\boldsymbol{n}}^{\boldsymbol{\upprime}}(\mathbf{y})}{-m\_{2}\boldsymbol{\upxi}\_{\boldsymbol{n}}^{\boldsymbol{\upprime}(\mathbf{1})}(\mathbf{x})\boldsymbol{\upmu}\_{\boldsymbol{n}}(\mathbf{y}) + m\_{1}\boldsymbol{\upxi}\_{\boldsymbol{n}}^{\boldsymbol{\upprime}(\mathbf{1})}(\mathbf{x})\boldsymbol{\upmu}\_{\boldsymbol{n}}^{\boldsymbol{\upprime}}(\mathbf{y})} \tag{31}$$

$$\mathbf{x} = m\_2 k\_0 a, \qquad \mathbf{y} = m\_l k\_0 a \tag{32}$$

Substituting Eqs. (14) - (19) and (23) - (28) into Eqs. (11) - (12), and after some algebra, we can get the formula for RPFs which can be characterized by radiation pressure cross section (RPCS):

$$\mathbf{F(r)} = \frac{\mathcal{Q}\_2 I\_0}{c} \left[ \hat{\mathbf{e}}\_\mathbf{r} C\_{pr,x}(\mathbf{r}) + \hat{\mathbf{e}}\_\mathbf{y} C\_{pr,y}(\mathbf{r}) + \hat{\mathbf{e}}\_\mathbf{z} C\_{pr,z}(\mathbf{r}) \right] \tag{33}$$

where RPCS *Cpr*,*<sup>i</sup>* (*i* = *x*, *y*, *z*) has a longitudinal cross section *Cpr*,*<sup>z</sup>*

$$\begin{split} C\_{pr,z} &= \frac{\lambda^2}{\pi} \sum\_{n=1}^{n} \text{Re} \left\{ \frac{1}{n+1} (A\_n \mathbf{g}\_{n,TM}^0 \mathbf{g}\_{n+1,TM}^{\circ \bullet} + B\_n \mathbf{g}\_{n,TL}^0 \mathbf{g}\_{n+1,TL}^{\circ \bullet}) + \sum\_{n=1}^{n} \left[ \frac{1}{(n+1)^2} \frac{(n+m+1)!}{(n-m)!} \right. \\ &\times \left( A\_n \mathbf{g}\_{n,TM}^{\circ \bullet} \mathbf{g}\_{n+1,TM}^{\circ \bullet,TM} + A\_n \mathbf{g}\_{n,TM}^{-\circ \bullet} \mathbf{g}\_{n+1,TM}^{\circ \bullet \bullet} + B\_n \mathbf{g}\_{n,TL}^{-\circ \bullet} \mathbf{g}\_{n+1,TL}^{\circ \bullet \bullet} \right) \\ &+ m \frac{2n+1}{n^2 (n+1)^2} \frac{(n+m)!}{(n-m)!} C\_n \left( \mathbf{g}\_{n,TM}^{\circ \bullet} \mathbf{g}\_{n,TE}^{\circ \bullet \bullet} - \mathbf{g}\_{n,TM}^{-\circ \bullet} \mathbf{g}\_{n,TE}^{-\circ \bullet \bullet} \right) \end{split} \tag{34}$$

and two transverse cross section *Cpr*,*<sup>x</sup>* and *Cpr*,*<sup>y</sup>*

$$C\_{\rho r,x} = \text{Re}(C) \qquad C\_{\rho r,y} = \text{Im}(C) \tag{35}$$

where

#### Fiber-Based Cylindrical Vector Beams and Its Applications to Optical Manipulation http://dx.doi.org/10.5772/59151 207

$$\begin{split} C &= \frac{\lambda^2}{2\pi} \sum\_{n=1}^{n} \left\{ -\frac{(2n+2)!}{(n+1)^2} F\_{\pi}^{n+1} + \sum\_{m=1}^{n} \frac{(n+m)!}{(n-m)!} \frac{1}{(n+1)^2} \\ &\times \left[ F\_{\pi}^{n+1} - \frac{n+m+1}{n-m+1} F\_{\pi}^{m} + \frac{2n+1}{n^2} (C\_n \mathfrak{g}\_{n,TM}^{n-1} \mathfrak{g}\_{n,\overline{n}}^{n\*} \\ &- C\_n \mathfrak{g}\_{n,TM}^{-n+1} \mathfrak{g}\_{n+1,\overline{n}}^{-n+1} + C\_n^\* \mathfrak{g}\_{n,\overline{n}}^{n-1} \mathfrak{g}\_{n,TM}^{n\*} - C\_n^\* \mathfrak{g}\_{n,\overline{n}}^{-n} \mathfrak{g}\_{n,\overline{n}}^{-n+1\*}) \right] \end{split} \tag{36}$$

with

0 | | | | ( ) 0 0

 q

<sup>=</sup> é ù - å å ë û

, , = = , *m m mm*

1 2 (1) ' (1) 1 2

2 1 (1) ' (1) 2 1

( )

 y

 y ( )

() ()

*n n n*

*n n nn m x ym x*

*n n n*

*n n nn m x ymx*

Substituting Eqs. (14) - (19) and (23) - (28) into Eqs. (11) - (12), and after some algebra, we can get the formula for RPFs which can be characterized by radiation pressure cross section (RPCS):

,,,

<sup>1</sup> 1 ( 1)! ( ) 1 ( 1) ( )!

<sup>2</sup> ( ) = ++ é ù ˆˆˆ () () () ë û *pr x pr y pr z*

0 0\* 0 0\* , , 1, , 1, 2 1 1

*n m C Re A g g B g g <sup>n</sup> <sup>n</sup> n m*

( )

= =

*Ag g Ag g Bg g Bg g*

*n m m m m m mm m m n n TM n TM n n TM n TM n n TE n TE n n TE n TE*

å å*<sup>n</sup>*

´ + ++

*n n*

*n nm <sup>m</sup> C g n n n m* } \* \*

\* \*\* \* , 1, , 1, , 1, , 1,


> , ,, ) - - - ù <sup>û</sup> *mm mm TM n TE n TM n TE g gg*

+ +

<sup>ì</sup> <sup>é</sup> + + <sup>=</sup> <sup>í</sup> + + <sup>ê</sup> <sup>î</sup> + - <sup>ë</sup> <sup>+</sup>

( ) ( ) ( ) ( )( )

 yy ¢ ¢

¢ ¢

*y <sup>m</sup> <sup>b</sup> m x <sup>y</sup>* (31)

*<sup>I</sup> CCC <sup>c</sup>* **Fr e r e r e r xyz** (33)

, , = = Re( ) Im( ) *C CC C pr x pr y* (35)

*y <sup>m</sup> <sup>a</sup> m x <sup>y</sup>* (30)

 xy

*y x*

() ()

( ) ( ) ( ) ( )( )

 yy

 xy

*y x*

*s pw m m m m im n nn n n nn n*

xt

206 Advances in Optical Fiber Technology: Fundamental Optical Phenomena and Applications

Where *ξn*(*k*0*r*) is Ricatti-Hankel functions, and scattering coefficients *An*

expressed by traditional Mie scattering coefficients *an*, *bn* and BSCs *gn*,*TM*

'

'

+- <sup>=</sup> - +

y y

x

2 0

where RPCS *Cpr*,*<sup>i</sup>* (*i* = *x*, *y*, *z*) has a longitudinal cross section *Cpr*,*<sup>z</sup>*

*pr z n n TM n TM n n TE n TE*

2 2 ,

2 1 ( )! ( ( 1) ( )!

+ +

+ -

and two transverse cross section *Cpr*,*<sup>x</sup>* and *Cpr*,*<sup>y</sup>*

2

¥

l

p

+

where

+- <sup>=</sup> - +

y y

x

*n*

*n*

( ) (cos ) ( ) (cos )

 x

*iH <sup>H</sup> c A kr mB k r <sup>e</sup> k r* (28)

¢

 p

*n n n TM n n n TE g Ba bA g* (29)

2 0 1 0 *x mk a y mk a* = = , (32)

j

*<sup>m</sup>* , *gn*,*TE <sup>m</sup>* :

*<sup>m</sup>* and *Bn*

*<sup>m</sup>* can be

(34)

 q

0 1

¥ +

*n*

= =-

*nmn*

j

with

$$F\_n^{n} = A\_n \mathbf{g}\_{n, \mathcal{IM}}^{n-1} \mathbf{g}\_{n+1, \mathcal{IM}}^{n^\*} + B\_n \mathbf{g}\_{n, \mathcal{IL}}^{n-1} \mathbf{g}\_{n+1, \mathcal{IL}}^{n^\*} + A\_n^{n^\*} \mathbf{g}\_{n+1, \mathcal{IM}}^{-n} \mathbf{g}\_{n, \mathcal{IM}}^{-n+1^\*} + B\_n^{n^\*} \mathbf{g}\_{n+1, \mathcal{IL}}^{-n} \mathbf{g}\_{n, \mathcal{IL}}^{-n+1^\*} \tag{37}$$

$$A\_u = a\_u + a\_{n+1}^\* - 2a\_n a\_{n+1}^\* \tag{38}$$

$$B\_n = b\_n + b\_{n+1}^\* - 2b\_n b\_{n+1}^\* \tag{39}$$

$$C\_n = -i(a\_n + b\_{n+1}^\* - 2a\_n b\_{n+1}^\*) \tag{40}$$

Note that substituting Eq. (13) into Eqs. (11) - (12) shows that the total RPF can be devided into thress parts:

$$<\mathbf{F}> = <\mathbf{F}\_i> + <\mathbf{F}\_{mix}> + <\mathbf{F}\_i> \tag{41}$$

where < *F<sup>i</sup>* > depends only on the incident fields, < *F<sup>s</sup>* > is associated with the scattered fields, and < *Fmix* > involves the interactions of the incident beam with the scattered field. After a great deal of algebra, we can get that < *F<sup>i</sup>* > =0, which can be understood by the momentum conser‐ vation law for monochromatic fields in free space. The RPCS for < *Fmix* > and < *F<sup>s</sup>* > can be directly given using Eqs. (34) - (40) by changing Eqs. (38) - (40) using

$$\begin{aligned} A\_n &= a\_n + a\_{n+1}^\* \\ B\_n &= b\_n + b\_{n+1}^\* \\ C\_n &= -i(a\_n + b\_{n+1}^\*) \end{aligned} \tag{42}$$

for < *Fmix* >, and

$$\begin{aligned} A\_n &= -2a\_n a\_{n+1}^\* \\ B\_n &= -2b\_n b\_{n+1}^\* \\ C\_n &= 2ia\_n b\_{n+1}^\* \end{aligned} \tag{43}$$

### for < *F<sup>s</sup>* >.

#### *3.2.2. Beam shape coefficients for CVB*

This section is devoted the derivation of BSCs for CVB using ILA. In the ILA, the beam shape coefficients *gn*,*TM <sup>m</sup>* and *gn*,*TE <sup>m</sup>* are obtained respectively from the radial component of electric and magnetic field *Er* and *Hr* according to [34]

$$\mathbf{g}'''\_{\pi, \Pi^{\pm}} = \frac{Z\_{\pi}''}{2\pi H\_{\theta 0}} \int\_{0}^{2\pi} \overline{H\_{r}}(r, \theta, \phi) \mathbf{e}^{-im\phi} d\phi \tag{44}$$

$$\mathbf{g}\_{n,\!TM}^{''} = \frac{Z\_{\
u}^{'''}}{2\pi E\_{\mu 0}} \int\_0^{2\pi} \overline{E\_r}(r, \theta, \phi) \mathbf{e}^{-in\phi} d\phi \tag{45}$$

with

$$Z\_{\kappa}^{m} = \begin{cases} \frac{2n(n+1)i}{2n+1} & m = 0\\ \left(\frac{-2i}{2n+1}\right)^{\|\kappa\|-1} & m \neq 0 \end{cases} \tag{46}$$

*Er* ¯ and *Hr* ¯ are respectively the localized fields of *Er* and *Hr*, and they are obtained by changing *kr* to (*n* + 1 / 2) and *θ* to *π* / 2 in their expression. For a radially polarized Bessel-Gaussian beam, the localized radial component of electric field are derived from Eq. (2):

$$\begin{split} \overline{E}'\_{rad} &= E\_0 \overline{\Omega}\_0 \left[ - (\rho\_\pi \cos \phi - \underline{\xi}\_0) \cos \phi + (\rho\_\pi \sin \phi - \eta\_0) \sin \phi \right] \\ &= E\_0 \overline{\Omega}\_0 \left[ \rho\_\pi - \rho\_0 \sin(\phi + \phi\_0) \right] \end{split} \tag{47}$$

with

$$\overline{\Omega}\_0 = -\sqrt{2} \exp\left[ - (\rho\_\kappa^2 + \xi\_0^2 + \eta\_0^2) \right] \exp\left[ 2\,\rho\_\kappa (\xi\_0 \cos\phi + \eta\_0 \sin\phi) \right] \tag{48}$$

$$
\rho\_{\pi} = \frac{kr}{k\nu\_0} = \frac{1}{k\nu\_0} \left( n + \frac{1}{2} \right) \tag{49}
$$

$$
\xi\_0 = \rho\_0 \sin \phi\_0, \quad \eta\_0 = \rho\_0 \cos \phi\_0 \tag{50}
$$

Substituting Eq. (47) into Eq. (45), and considering the formula of Bessel function

$$J\_{\pi}(\mathbf{x}) = \frac{1}{2\pi} \int\_{-\pi}^{\pi} e^{i(x\sin\theta - a\theta)} d\theta = \frac{1}{2\pi} \int\_{0}^{2\pi} e^{i(x\sin\theta - a\theta)} d\theta \tag{51}$$

we can obtain the final expression of BSCs

$$\log\_{n, \text{TM}}^{n, rad} = \frac{1}{2} Z\_n^n \overline{\Omega}\_n e^{im\phi\_0} \left[ 2\rho\_n J\_n(-2i\rho\_n \rho\_0) + i\rho\_0 (J\_{n-1}(-2i\rho\_n \rho\_0) - J\_{n+1}(-2i\rho\_n \rho\_0)) \right] \tag{52}$$

Here we only derive *gn*,*TM <sup>m</sup>*,*rad* , and *gn*,*TE <sup>m</sup>*,*rad* can be derived in the similar way.

#### *3.2.3. Debye series expansion*

for < *F<sup>s</sup>* >.

with

with

coefficients *gn*,*TM*

*3.2.2. Beam shape coefficients for CVB*

*<sup>m</sup>* and *gn*,*TE*

magnetic field *Er* and *Hr* according to [34]

This section is devoted the derivation of BSCs for CVB using ILA. In the ILA, the beam shape

2

2

2 1

<sup>ï</sup> <sup>+</sup> <sup>=</sup> <sup>í</sup>

*n m nn i <sup>m</sup> <sup>n</sup> <sup>Z</sup>*

*m*

the localized radial component of electric field are derived from Eq. (2):

00 0 0

r fx

=W - + *rad n n*

> r x h

r

xr

r r ff

*<sup>r</sup> E E*

*n*

[ ]

*g Er d*

( , , )e <sup>2</sup> p


( , , )e <sup>2</sup> p


*g Hr d*

, <sup>0</sup> <sup>0</sup>

208 Advances in Optical Fiber Technology: Fundamental Optical Phenomena and Applications

*n TE r B Z*

*m m n im*

p

, <sup>0</sup> <sup>0</sup>

p

*n TM r B Z*

*<sup>m</sup>* are obtained respectively from the radial component of electric and

f

f

 f

 f

*<sup>H</sup>* (44)

*<sup>E</sup>* (45)

 f

 f+ ë û *n n* (48)

*kr <sup>n</sup> kw kw* (49)

= = sin , cos (50)

*<sup>E</sup>* (47)

 fh (46)

qf

qf

1 2 ( 1) <sup>0</sup> 2 1

*Er* ¯ and *Hr* ¯ are respectively the localized fields of *Er* and *Hr*, and they are obtained by changing

*kr* to (*n* + 1 / 2) and *θ* to *π* / 2 in their expression. For a radially polarized Bessel-Gaussian beam,

[ ]

 rx

> f

1 1 2

( cos )cos ( sin )sin

 f r fh

0 0 0 0

[ ] <sup>222</sup> <sup>0</sup> 0 0 0 0 W =- - + + 2 exp ( ) exp 2 ( cos sin ) é ù

0 0

00 000 0

 fhr

Substituting Eq. (47) into Eq. (45), and considering the formula of Bessel function

æ ö == + ç ÷ è ø *<sup>n</sup>*

sin( )

= W- - + -

 <sup>ì</sup> <sup>+</sup> <sup>=</sup> <sup>ï</sup>

æ ö - <sup>ï</sup> <sup>¹</sup> ç ÷ <sup>ï</sup> îè ø +

*<sup>i</sup> <sup>m</sup> <sup>n</sup>*

<sup>2</sup> <sup>0</sup>

GLMT is a rigorous electromagnetic theory, and can exactly predict the RPF exerted on a sphere by focused beam. Whereas the solution is complicated combinations of Bessel functions, and the mathematical complexity obscures the physical interpretation of various features of RPF. The DSE, which is a rigorous electromagnetic theory, expresses the Mie scattering coefficients into a series of Fresnel coefficients and gives physical interpretation of different scattering processes. The DSE is an efficient technique to make explicit the physical interpretation of various features of RPF which are implicit in the GLMT. The DSE is firstly presented by Debye in 1908 for the interaction between electromagnetic waves and cylinders. Since then, the DSE for electromagnetic scattering by homogeneous, coated, multilayered spheres, multilayered cylinders at normal incidence, homogeneous, multilayered cylinder at oblique incidence, and spherical gratings are studied. In our previous work, DSE has been employed to the analysis of RPF exerted on a sphere induced by a Gaussian and Bessel beam.

As shown in Fig. 3, when an incoming spherical multipole wave, which is

$$\Psi = \xi\_n^{(l)}(m\_2kr)P\_n^{\
u}(\cos\theta)\begin{Bmatrix}\cos m\phi\\\sin m\phi\end{Bmatrix},\tag{53}$$

encounters the interface of the sphere at *r=a*, portion of it will be transmitted into the sphere, and another portion will reflected back. The transmitted and reflected waves are respectively:

$$\Psi\_1 = T\_{\
u}^{\;\;21} \xi\_{\mu}^{(l)}(m\_l k r) P\_{\
u}^{\;\;n}(\cos \theta) \begin{cases} \cos m\phi \\ \sin m\phi \end{cases} \qquad r \le a \tag{54}$$

$$\Psi\_2' = \left[ \xi\_n^{(1)}(m\_2kr) + R\_n^{212}\xi\_n^{(2)}(m\_2kr) \right] P\_n^m(\cos\theta) \begin{cases} \cos m\phi \\ \sin m\phi \end{cases} \qquad r \ge a \tag{55}$$

**Figure 3.** Debye model of light scattering by a sphere

Applying the boundary conditions, which reqires continuity of the tangential components of filds at the interface, to the incident, transmitted and reflected waves,we can obtain:

$$T\_{\pi}^{21} = \frac{m\_1}{m\_2} \frac{2i}{D\_{\pi}} \tag{56}$$

$$R\_{\pi}^{212} = \frac{\alpha \xi\_{\pi}^{\varepsilon(l)}(\kappa\_{2})\xi\_{\pi}^{\varepsilon(l)}(\kappa\_{1}) - \beta \xi\_{\pi}^{\varepsilon(l)}(\kappa\_{2})\xi\_{\pi}^{\varepsilon(l)}(\kappa\_{1})}{D\_{\pi}} \tag{57}$$

with

$$D\_n = -\alpha \xi\_n^{\varepsilon(2)^\cdot} (\kappa\_2) \xi\_n^{\varepsilon(l)} (\kappa\_1) + \beta \xi\_n^{\varepsilon(2)} (\kappa\_2) \xi\_n^{\varepsilon(l)^\cdot} (\kappa\_1) \tag{58}$$

$$
\kappa\_j = m\_j ka
\tag{59}
$$

$$\alpha = \begin{cases} \text{l}, & \text{for TE} \\ \frac{m\_{\text{l}}}{m\_{\text{l}}}, & \text{for TM} \end{cases}, \qquad \beta = \begin{cases} \frac{m\_{\text{l}}}{m\_{\text{l}}}, & \text{for TE} \\ \text{m}\_{\text{l}} \\ \text{l}, & \text{for TM} \end{cases} \tag{60}$$

Similarly, the consideration of outgoing multipole waves can get

$$T\_s^{12} = \frac{2i}{D\_s} \tag{61}$$

$$R\_{\boldsymbol{n}}^{121} = \frac{\alpha \xi\_{\boldsymbol{n}}^{\varepsilon(2)}(\kappa\_2)\xi\_{\boldsymbol{n}}^{\varepsilon(2)}(\kappa\_1) - \beta \xi\_{\boldsymbol{n}}^{\varepsilon(2)}(\kappa\_2)\xi\_{\boldsymbol{n}}^{\varepsilon(2)}(\kappa\_1)}{D\_{\boldsymbol{n}}} \tag{62}$$

Substituting all Fresnel coefficients into

$$\left(1 - R\_n^{121}\right)\left(1 - R\_n^{212}\right) - T\_n^{21} T\_n^{12} \tag{63}$$

and after much algebra, we get

*p* = 0

*p* = 2

ax k x k bx k x k

ax k x k bx k x k

2

a

*n*

*R*

with

**Figure 3.** Debye model of light scattering by a sphere

<sup>212</sup> *Rn*

210 Advances in Optical Fiber Technology: Fundamental Optical Phenomena and Applications

21 *Tn*

*p* =1

12 *Tn*

*m D* (56)

*<sup>D</sup>* (57)

*j j* = *m ka* (59)

(60)

*m*1

Applying the boundary conditions, which reqires continuity of the tangential components of

*n*

 () () () () ¢ ¢ - <sup>=</sup> *n n nn*

(1) (1) (1) (1) 212 21 2 1

> (2) (1) (2) (1) 21 2 1

> > b

*for TE m for TE <sup>m</sup> <sup>m</sup> for TM <sup>m</sup> for TM*

k

ì ì ï ï = = í í ï ï î î

1 2

, 1,

1, , ,

 () () () () ¢ ¢

1

*Dnn n nn* = - + (58)

*n*

filds at the interface, to the incident, transmitted and reflected waves,we can obtain:

21 1 2 2

*<sup>m</sup> <sup>i</sup> <sup>T</sup>*

*<sup>n</sup>* =

<sup>121</sup> *Rn*

*m*2

$$\begin{split} \begin{Bmatrix} a\_{\pi} \\ b\_{\pi} \end{Bmatrix} &= \frac{1}{2} \left[ 1 - R\_{\pi}^{212} - \frac{T\_{\pi}^{21} T\_{\pi}^{12}}{1 - R\_{\pi}^{121}} \right] \\ &= \frac{1}{2} \left[ 1 - R\_{\pi}^{212} - \sum\_{\rho=1}^{n} T\_{\pi}^{21} \left( R\_{\pi}^{121} \right)^{\rho - 1} T\_{\pi}^{12} \right] \end{split} \tag{64}$$

where the prime indicates the derivative of the function with respect to its argument. *ξ<sup>n</sup>* (1) ( ·) and *ξ<sup>n</sup>* (2) ( ·) are respectively the spherical Ricatti-Hankel functions of first and second kinds. The definition of all Fresnel coefficients and Debye term *p* are given in Fig. 3. For convenience, we note p = - 1 and p = 0 respectively for the diffraction and direct reflection. In our previous work, we have theoretically and numerically proved that when *p* ranges from 1 to , the Eq. (64) is identical to the traditional Mie scattering coefficients. Here we provide the DSE for homogeneous spheres, and the DSE for multilayered spheres can be found in our previous work.

#### **4. Numerical results and discussions**

In this section, the GLMT and DSE will be employed to analyze the RPF exerted on a homo‐ geneous sphere induced by a radially polarized vector Bessel-Gaussian beam. Xu et al. used GLMT to analyze the RPF exerted on a slightly volatile silocone oil of refractive index , which can be levitated in the air by a beam of wavelength . We first use GLMT to analyze the RPF exerted on such oil induced by vector Bessel-Gaussian beam, and DSE will be employed to the study of the contribution of various scattering process to RPF. In our calculation, the radius of the particle is .

We first explore the influence of beam center location on the RPF. In our calculation, we assume the beam center is located on the x axis so that . Fig. 4 gives the transverse RPCS *∞* versus *m*<sup>1</sup> =1.5 for various beam-waist radius *λ* =0.5*µm*. Here we consider *a* =2.5*µm* and *y*<sup>0</sup> = *z*<sup>0</sup> =0, which are respectively larger, equal and smaller than the radius of the particle. We can find that the particle can not be trapped at the beam center *Cpr*,*<sup>x</sup>* for all beams. This results from the fact all beams have null central intensity. It is worth pointing out that a stable trap corresponds to a particle position where the RPF is zero and its slope is positive. All curves have two equilibrium points, which are symmetric with respect to beam axis (*x*0). This is decided by the intensity maxima of beams. So a vector Bessel-Gaussian beam can simultaneously trap more than one particles. We can also find that the interval of equilibrium points increases with the increasing of *w*0. This can be easily explained from the fact that the interval of intensity peaks increases with the increasing of *w*<sup>0</sup> =5*µm*, 2.5*µm*.

**Figure 4.** Transverse cross-section 1*µm* versus *x*<sup>0</sup> =0 with parameter *x*<sup>0</sup> =0. *w*0, *w*0, *Cpr*,*x*, *x*0, *w*0 and *λ* =0.5*µm*

To clarify the physical explanation of some features of RPCS, it is necessary for us to consider the contribution of each mode *p* to RPCS. The contribution of each *p* mode to RPCS can be computed separately by considering a single term in Eq. (64). Now we consider the contribu‐ tion of a single *p* mode to transverse RPCS *m*<sup>1</sup> =1.5. Here we set beam-waist radius *m*<sup>2</sup> =1.

It is shown in Fig. 5 the transverse RPCS *y*<sup>0</sup> =0 versus *a* =2.5*µm* with parameter *pmax* =*∞*. In the calculation, the beam-waist radius is assumed *Cpr*,*<sup>x</sup>*. Comparison of Fig. 5 with Fig. 4 shows that when *w*<sup>0</sup> =5*µm* the results obtained by DSE are identical to those by GLMT. In fact, when

100 *max p* ,<\$%&?>the<\$%&?>results<\$%&?>of<\$%&?>DSE<\$%&?>is<\$%&?>very<\$%&?>close<\$%&?>to<\$%&?>GLMT

%&?>of<\$%&?> 1 *max p* .<\$%&?>Fig.<\$%&?>5<\$%&?>shows<\$%&?>that<\$%&?>when<\$%&?>

<\$%&?>for<\$%&?> *p* 1 <\$%&?>and<\$%&?> *p* 2 .<\$%&?>Generally,<\$%&?>the<\$%&?>RPF<\$%&?>at<\$%&?> <sup>0</sup> *x* 0 is<\$%&?>zero<\$%&?>for<\$%&?>any<\$%&?>mode<\$%&?>*p*<\$%&?>because<\$%&?>of<\$%&?>the<\$%&?>symmetry.<\$% &?>The<\$%&?>magnitude<\$%&?>of<\$%&?> *Cpr x*, <\$%&?>for<\$%&?> *p* 1

*p* 2 .<\$%&?>This<\$%&?>validates<\$%&?>that<\$%&?>the<\$%&?>transverse<\$%&?>RPCS<\$%&?> *Cpr x*,

<\$%&?>is<\$%&?>dominated<\$%&?>by<\$%&?>the<\$%&?>contributions<\$%&?>of<\$%&?>direct<\$%&?>transmission<\$% &?>( *p* 1 ).<\$%&?>Note<\$%&?>that<\$%&?>the<\$%&?>curve<\$%&?>for<\$%&?> *p* 2

<sup>0</sup> *x* 0 .<\$%&?>Near<\$%&?>the<\$%&?>beam<\$%&?>axis<\$%&?>( <sup>0</sup> *x* 0 ),<\$%&?>the<\$%&?>curvesfor<\$%&?> *p* 1 <\$%&?>has<\$%&?>positive<\$%&?>slope,<\$%&?>while<\$%&?>that<\$%&?>for<\$%&?> *p* 2

has<\$%&?>negative<\$%&?>one.<\$%&?>To<\$%&?>explain<\$%&?>such<\$%&?>phenomena,<\$%&?>we<\$%&?>must<\$%

&?>consider<\$%&?>the<\$%&?>integral<\$%&?>effect<\$%&?>of<\$%&?>all<\$%&?>intensity<\$%&?>peaks.

*m* ,<\$%&?>while<\$%&?>the<\$%&?>curve<\$%&?>for<\$%&?> *p* 1

and<\$%&?> <sup>0</sup> *w m* <sup>5</sup>

.

should<\$%&?>be<\$%&?>negilible.<\$%&?>For<\$%&?>example,<\$%&?>if<\$%&?>

*Cpr*,*<sup>x</sup>* is large enough, the difference between two theories should be negilible. For example, if *x*0, the results of DSE is very close to GLMT results. Special attention should be paid to the case of *pmax*. Fig. 5 shows that when *w*<sup>0</sup> =5*µm*, the agreement between the results of GLMT and DSE is already good. This concludes that main contribution of RPF comes from the scattering processes of diffraction (*pmax* →*∞*), specular reflection (*pmax*) and direct transmission ( *pmax* =100). 1 *max p* ,<\$%&?>the<\$%&?>agreement<\$%&?>between<\$%&?>the<\$%&?>results<\$%&?>of<\$%&?>GLMT<\$%&?>and<\$ %&?>DSE<\$%&?>is<\$%&?>already<\$%&?>good.<\$%&?>This<\$%&?>concludes<\$%&?>that<\$%&?>main<\$%&?>contrib ution<\$%&?>of<\$%&?>RPF<\$%&?>comes<\$%&?>from<\$%&?>the<\$%&?>scattering<\$%&?>processes<\$%&?>of<\$%&?>d iffraction<\$%&?>( *p* 1),<\$%&?>specular<\$%&?>reflection<\$%&?>( *p* 0 )<\$%&?>and<\$%&?>direct<\$%&?>transmissi on<\$%&?>( *p* 1 ).

We first explore the influence of beam center location on the RPF. In our calculation, we assume the beam center is located on the x axis so that . Fig. 4 gives the transverse RPCS *∞* versus *m*<sup>1</sup> =1.5 for various beam-waist radius *λ* =0.5*µm*. Here we consider *a* =2.5*µm* and *y*<sup>0</sup> = *z*<sup>0</sup> =0, which are respectively larger, equal and smaller than the radius of the particle. We can find that the particle can not be trapped at the beam center *Cpr*,*<sup>x</sup>* for all beams. This results from the fact all beams have null central intensity. It is worth pointing out that a stable trap corresponds to a particle position where the RPF is zero and its slope is positive. All curves have two equilibrium points, which are symmetric with respect to beam axis (*x*0). This is decided by the intensity maxima of beams. So a vector Bessel-Gaussian beam can simultaneously trap more than one particles. We can also find that the interval of equilibrium points increases with the increasing of *w*0. This can be easily explained from the fact that the interval of intensity peaks increases


x0 (m)

**Figure 4.** Transverse cross-section 1*µm* versus *x*<sup>0</sup> =0 with parameter *x*<sup>0</sup> =0. *w*0, *w*0, *Cpr*,*x*, *x*0, *w*0 and *λ* =0.5*µm*

To clarify the physical explanation of some features of RPCS, it is necessary for us to consider the contribution of each mode *p* to RPCS. The contribution of each *p* mode to RPCS can be computed separately by considering a single term in Eq. (64). Now we consider the contribu‐ tion of a single *p* mode to transverse RPCS *m*<sup>1</sup> =1.5. Here we set beam-waist radius *m*<sup>2</sup> =1.

It is shown in Fig. 5 the transverse RPCS *y*<sup>0</sup> =0 versus *a* =2.5*µm* with parameter *pmax* =*∞*. In the calculation, the beam-waist radius is assumed *Cpr*,*<sup>x</sup>*. Comparison of Fig. 5 with Fig. 4 shows that when *w*<sup>0</sup> =5*µm* the results obtained by DSE are identical to those by GLMT. In fact, when

with the increasing of *w*<sup>0</sup> =5*µm*, 2.5*µm*.





0.0

C

pr,x

0.2

0.4

0.6

0.8

w0

w0

w0

=1.0 m

212 Advances in Optical Fiber Technology: Fundamental Optical Phenomena and Applications

=2.5 m

=5.0 m

Figure 5. Transverse<\$%&?>cross-section<\$%&?>*Cpr x*, versus<\$%&?> <sup>0</sup> *<sup>x</sup>* <\$%&?>with<\$%&?>parameter<\$%&?> *max <sup>p</sup>* .<\$%&?> **Figure 5.** Transverse cross-section *pmax* =1 versus *pmax* =1 with parameter *p* = −1. *p* =0, *p* =1, *Cpr*,*x*, *x*0, *pmax* and *λ* =0.5*µm*.

 0.5*<sup>m</sup>* ,<\$%&?> <sup>1</sup> *<sup>m</sup>* 1.5 ,<\$%&?> <sup>2</sup> *<sup>m</sup>* <sup>1</sup> ,<\$%&?> <sup>0</sup> *<sup>y</sup>* <sup>0</sup> ,<\$%&?> *a m* 2.5To<\$%&?>clarify<\$%&?>the<\$%&?>physical<\$%&?>explanation<\$%&?>of<\$%&?>some<\$%&?>phenomena<\$%&?>of<\$ %&?>RPF,<\$%&?>it<\$%&?>is<\$%&?>necessary<\$%&?>to<\$%&?>consider<\$%&?>the<\$%&?>contribution<\$%&?>of<\$% &?>each<\$%&?>mode<\$%&?>*p*<\$%&?>to<\$%&?>RPCS,<\$%&?>which<\$%&?>can<\$%&?>be<\$%&?>computed<\$%&?>se parately<\$%&?>by<\$%&?>considering<\$%&?>a<\$%&?>single<\$%&?>term<\$%&?>in<\$%&?>Eq.<\$%&?>(64).<\$%&?>No w<\$%&?>we<\$%&?>consider<\$%&?>the<\$%&?>contribution<\$%&?>of<\$%&?>a<\$%&?>single<\$%&?>*p*<\$%&?>mode<\$ %&?>to<\$%&?>transverse<\$%&?>RPCS<\$%&?> *Cpr x*, .<\$%&?>Fig<\$%&?>.6<\$%&?>depicts<\$%&?>the<\$%&?>transverse<\$%&?>RPCS<\$%&?> *Cpr x*, versus<\$%&?> <sup>0</sup> *x* To clarify the physical explanation of some phenomena of RPF, it is necessary to consider the contribution of each mode *p* to RPCS, which can be computed separately by considering a single term in Eq. (64). Now we consider the contribution of a single *p* mode to transverse RPCS *m*<sup>1</sup> =1.5. Fig.6 depicts the transverse RPCS *m*<sup>2</sup> =1 versus *y*<sup>0</sup> =0 for *a* =2.5*µm* and *w*<sup>0</sup> =5*µm*. Gener‐ ally, the RPF at *Cpr*,*<sup>x</sup>* is zero for any mode *p* because of the symmetry. The magnitude of *Cpr*,*<sup>x</sup>* for *x*0 is much greater than that for *p* =1. This validates that the transverse RPCS *p* =2 is dominated by the contributions of direct transmission (*x*<sup>0</sup> =0). Note that the curve for *Cpr*,*<sup>x</sup>* has two equilibrium points at about *p* =1, while the curve for *p* =2 has only one points at *Cpr*,*<sup>x</sup>*. Near the beam axis (*p* =1), the curvesfor *p* =2 has positive slope, while that for *x*<sup>0</sup> = ± 1.3*µm* has negative one. To explain such phenomena, we must consider the integral effect of all intensity peaks.

<sup>0</sup> *x* 1.3

<\$%&?>is<\$%&?>much<\$%&?>greater<\$%&?>than<\$%&?>that<\$%&?>for<\$%&?>

<\$%&?>has<\$%&?>only<\$%&?>one<\$%&?>points<\$%&?>at<\$%&?>

<\$%&?>has<\$%&?>two<\$%&?>equilibrium<\$%&?>points<\$%&?>at<\$%&?>about<\$%&?>

14 we must consider the integral effect of all intensity peaks.

Note that the curve for *p* 2 has two equilibrium points at about <sup>0</sup> 11 *x* 1.3


C

<sup>2</sup> *m* 1 , <sup>0</sup> *y* 0 , *a m* 2.5

and <sup>0</sup> 3 *w m* 5

1

pr,x

14 **Optical Fiber**

 pmax pmax=100 pmax=1

**Fig. 5.** Transverse cross-section *Cpr x*, versus <sup>0</sup> *x* with parameter *max p* .


 0.5

*m* , while the

x0 (m)

*m* , <sup>1</sup> 2 *m* 1.5 ,

4 To clarify the physical explanation of some phenomena of RPF, it is necessary to consider 5 the contribution of each mode *p* to RPCS, which can be computed separately by considering 6 a single term in Eq. (64). Now we consider the contribution of a single *p* mode to transverse RPCS *Cpr x*, . Fig .6 depicts the transverse RPCS *Cpr x*, versus <sup>0</sup> 7 *x* for *p* 1 and *p* 2 . Generally, the RPF at <sup>0</sup> 8 *x* 0 is zero for any mode *p* because of the symmetry. The magnitude of *Cpr x*, 9 for *p* 1 is much greater than that for *p* 2 . This validates that the transverse RPCS *Cpr x*, 10 is dominated by the contributions of direct transmission ( *p* 1 ).

curve for *p* 1 has only one points at <sup>0</sup> *x* 0 . Near the beam axis ( <sup>0</sup> 12 *x* 0 ), the curvesfor

.

**Figure 6.** Transverse cross-section *p* =1 versus *x*<sup>0</sup> =0 corresponding to single mode *p*. *x*<sup>0</sup> =0, *p* =1, *p* =2, *Cpr*,*<sup>x</sup> x*<sup>0</sup> and *λ* =0.5*µm*.

### **5. Conclusions**

15

Rigorous theories including GLMT and DSE for RPF exerted on spheres induced by CVB is derived. The incident beam is described by a set of BSCs which is calculated by integral localized approximation, and the scattering coefficients are expanded using Debye series. For very small particles, namely Rayleigh particles, an approximation model is also given. Such thoery can be easily extended to the RPF exerted on multilayered sphere, and also to the RPF induced by other beams. Debye series is used to isolate the contribution of various scattering process to the RPF. The results are of special importance for the improvement of optical tweezers system.

#### **Acknowledgements**

The authors acknowledge support from the Natural Science Foundation of China (Grant No. 61101068), the National Science Foundation for Distinguished Young Scholars of China (Grant No.61225002), and the Fundamental Research Funds for the Central Universities.

### **Author details**

Renxian Li\* , Lixin Guo, Bing Wei, Chunying Ding and Zhensen Wu

\*Address all correspondence to: rxli@mail.xidian.edu.cn

School of Physics and Optoelectronic Engineering, Xidian University, China
