Integral localized approximation description of ordinary Bessel beams and application to optical trapping forces

doi:10.1364/BOE.2.001893

« Show journal navigation

Integral localized approximation description of ordinary Bessel beams and application to optical trapping forces

Leonardo A. Ambrosio and Hugo E. Hernández-Figueroa »View Author Affiliations

Biomedical Optics Express, Vol. 2, Issue 7, pp. 1893-1906 (2011)
http://dx.doi.org/10.1364/BOE.2.001893

View Full Text Article

Acrobat PDF (1481 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
Introduction
The GLMT and the ILA applied to ordinary Bessel beams
Numerical results for the BSC’s of an ordinary Bessel beam
Application to optical force calculations
Conclusions
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (30)
Cited By
Figures (6)
Metrics

Abstract

Ordinary Bessel beams are described in terms of the generalized Lorenz-Mie theory (GLMT) by adopting, for what is to our knowledge the first time in the literature, the integral localized approximation for computing their beam shape coefficients (BSCs) in the expansion of the electromagnetic fields. Numerical results reveal that the beam shape coefficients calculated in this way can adequately describe a zero-order Bessel beam with insignificant difference when compared to other relative time-consuming methods involving numerical integration over the spherical coordinates of the GLMT coordinate system, or quadratures. We show that this fast and efficient new numerical description of zero-order Bessel beams can be used with advantage, for example, in the analysis of optical forces in optical trapping systems for arbitrary optical regimes.

1. Introduction

The generalized Lorenz-Mie theory (GLMT) is an extension of the Lorenz-Mie theory [1

G. Mie, “Beiträge zur Optik Trüber Medien, Speziell Kolloidaler Metallösungen,” Ann. Phys. 330(3), 377–452 (1908). [CrossRef]

] for describing the electromagnetic field components of an arbitrary laser beam in terms of spherical harmonic functions [2

G. Gouesbet and G. Gréhan, “Sur la généralisation de la théorie de Lorenz-Mie,” J. Opt. (Paris) 13, 97–103 (1982).

B. Maheu, G. Gouesbet, and G. Gréhan, “A concise presentation of the generalized Lorenz-Mie theory for arbitrary incident profile,” J. Opt. (Paris) 19, 59–67 (1988).

], the coefficients of which being called the beam shape coefficients (BSCs), responsible for correctly modeling the intensity profile of the beam [4

G. Gouesbet, G. Gréhan, and B. Maheu, “Expressions to compute the coefficients gnm in the generalized Lorenz-Mie theory using finite series,” J. Opt. (Paris) 19, 35–48 (1988).

]. Their numerical evaluation using quadratures [2

G. Gouesbet and G. Gréhan, “Sur la généralisation de la théorie de Lorenz-Mie,” J. Opt. (Paris) 13, 97–103 (1982).

] or finite series [4

G. Gouesbet, G. Gréhan, and B. Maheu, “Expressions to compute the coefficients gnm in the generalized Lorenz-Mie theory using finite series,” J. Opt. (Paris) 19, 35–48 (1988).

], however, can be a pretty time-consuming, lengthy or awkward task, first because of numerical integrations over the spherical coordinates of the adopted coordinate system, and second because of the inexistence of a single expression, in the latter case. Thus, efficient techniques such as the localized interpretation and, as its natural consequence, the integral localized approximation (ILA), have been developed for calculating on and off-axis BSCs with insignificant difference and much faster than the above mentioned numerical integrations [5

G. Gouesbet, G. Gréhan, and B. Maheu, “On the generalized Lorenz-Mie theory: first attempt to design a localized approximation to the computation of the coefficients gnm,” J. Opt. (Paris) 20, 31–43 (1989).

–7

K. F. Ren, G. Gouesbet, and G. Gréhan, “Integral localized approximation in generalized lorenz-mie theory,” Appl. Opt. 37(19), 4218–4225 (1998). [CrossRef] [PubMed]

]. Anyway, regardless of the scheme adopted for evaluating the BSCs, since the development of the GLMT plenty of applications have been benefited by this theoretical methodology (for a resume of the updated subject see, e.g., Refs. [8

J. A. Lock and G. Gouesbet, “Generalized Lorenz-Mie theory and applications,” J. Quant. Spectrosc. Radiat. Transf. 110(11), 800–807 (2009). [CrossRef]

–10

G. Gouesbet, J. A. Lock, and G. Gréhan, “Generalized Lorenz-Mie theories and description of electromagnetic arbitrary shaped beams: localized approximations and localized beam models,” J. Quant. Spectrosc. Radiat. Transf. 112(1), 1–27 (2011). [CrossRef]

]. and references therein).

The localized approximation developed by Gouesbet et al. is based on the localization principle of van de Hulst [11

H. C. van de Hulst, Light Scattering by Small Particles (Wiley, New York, 1957).

] and was rigorously justified for on- and off-axis focused Gaussian beams and laser sheets [12

G. Gouesbet and J. A. Lock, “Rigorous justification of the localized approximation to the beam-shape coefficients in generalized Lorenz-Mie theory. I. On-axis beams,” J. Opt. Soc. Am. A 11(9), 2503–2515 (1994). [CrossRef]

,13

G. Gouesbet and J. A. Lock, “Rigorous justification of the localized approximation to the beam-shape coefficients in generalized Lorenz-Mie theory. II. Off-axis beams,” J. Opt. Soc. Am. A 11(9), 2516–2525 (1994). [CrossRef]

]. Its improved version, viz., the integral localized approximation, has been used during the years for biomedical optics research such as in optical force and torque calculations exerted on both homogeneous and multilayered spherical dielectric particles in optical trapping systems [14

K. F. Ren, G. Gréhan, and G. Gouesbet, “Radiation pressure forces exerted on a particle arbitrarily located in a Gaussian beam by using the generalized Lorenz-Mie theory, and associated resonance effects,” Opt. Commun. 108(4-6), 343–354 (1994). [CrossRef]

,15

H. Polaert, G. Gréhan, and G. Gouesbet, “Forces and torques exerted on a multilayered spherical particle by a focused Gaussian beam,” Opt. Commun. 155(1-3), 169–179 (1998). [CrossRef]

]. Even optical forces exerted on hypothetical negative refractive index metamaterials by Gaussian beams have been analyzed with the ILA in the context of the GLMT [16

L. A. Ambrosio and H. E. Hernández-Figueroa, “Fundamentals of negative refractive index optical trapping: forces and radiation pressures exerted by focused Gaussian beams using the generalized Lorenz-Mie theory,” Biomed. Opt. Express 1(5), 1284–1301 (2010). [CrossRef] [PubMed]

In optical trapping experiments, however, there are actually other laser beams of relevant importance, due to their increasing interest because of their properties such as multiple trapping and angular momentum transfer [17

J. Arlt, V. Garces-Chavez, W. Sibbett, and K. Dholakia, “Optical micromanipulation using a Bessel light beam,” Opt. Commun. 197(4-6), 239–245 (2001). [CrossRef]

–20

K. C. Neuman and S. M. Block, “Optical trapping,” Rev. Sci. Instrum. 75(9), 2787–2809 (2004). [CrossRef] [PubMed]

]. Among those we found the Bessel beams.

Single Bessel beams furnish an easy way for multiple trapping, although their extended focus makes a full three-dimensional trap unachievable. Several methods have been developed for theoretically predicting the behavior of dielectric spherical particles, under the influence of a Bessel beam, in the ray optics regime, where the wavelength λ is much smaller than the diameter d of the particle (λ << d), and in the case where λ >> d, for which a dipole model can be assumed. Also, in the Rayleigh-Gans regime, simple formulas can be used to predict the optical trapping forces with relative success [21

A. N. Rubinov, A. A. Afanas’ev, I. E. Ermolaev, Y. A. Kurochkin, and S. Y. Mikhnevich, “Localization of spherical particles under the action of gradient forces in the field of a zero-order Bessel beam. Rayleigh-Gans approximation,” J. Appl. Spectrosc. 70(4), 565–572 (2003). [CrossRef]

]. Even recent works that use the GLMT for studying transverse dynamics of silicon particles under the influence of Bessel beams still formulate their theory based on slightly modified versions of quadratures [22

G. Milne, K. Dholakia, D. McGloin, K. Volke-Sepulveda, and P. Zemánek, “Transverse particle dynamics in a Bessel beam,” Opt. Express 15(21), 13972–13987 (2007). [CrossRef] [PubMed]

], explicitly emphasizing the time-demanding character of such approach as a function of the particle radius. The inclusion of the ILA for multi-ringed beams in the GLMT may fulfill a bottleneck in this theory and, therefore, really seems to deserve some special attention.

In this way, this paper is devoted to the analysis of the beam shape coefficients in order to describe an ordinary (zero-order) Bessel beam using the GLMT with the integral localized approximation and is organized as follows. In Section 2 we theoretically deduce closed-form solutions of the BSC’s for linearly polarized Bessel beams, which can be automatically extended to circularly polarized beams by applying simple symmetry considerations. In Section 3, numerical results are presented and compared with quadrature methods, which are exact numerical integrations directly derived from the GLMT based on orthogonality conditions. This reveals the efficient and fast character of the ILA algorithm, enabling us to describe, with enough accuracy, the electromagnetic fields of an ordinary Bessel beam. Section 4 shows that the ILA, applied to Bessel beams, makes optical trapping force calculations over spherical micro particles in the Mie regime very simple and accurate when compared to previous works in ray optics, Rayleigh-Gans and Mie regimes. In section 5, our conclusions are presented.

2. The GLMT and the ILA applied to ordinary Bessel beams

Consider an ideal monochromatic zero-order Bessel beam (BB) propagating with speed c and angular frequency ω ₀ parallel to +z, with its optical axis displaced ρ ₀ = (x ₀,y ₀) from the z-axis and making an angle ϕ ₀ relative to the x-axis, according to Fig. 1 . Although we could have shifted the beam in a three-dimensional fashion, we must take into account that, as long as the ideal definition is considered and that the corresponding shift in z be constricted to the focal length (thus keeping the radial Bessel profile still valid), this implies in a more complete but unnecessary mathematical formulation. It is well known that the longitudinal and transverse wave numbers for a BB are given by k_z = kcosθ_a and k_ρ = ksinθ_a , respectively, θ_a being the associated axicon angle and k = ω ₀/c the wavenumber [23

J. Durnin, J. J. Miceli Jr., and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58(15), 1499–1501 (1987). [CrossRef] [PubMed]

]. If this beam is linearly polarized in x or in y, we can write its electric field as

Fig. 1 Geometrical description of an ordinary Bessel beam propagating parallel to + z (out of the page). The optical axis makes an angle ϕ ₀ relative to the x-axis and is displaced ρ ₀ from the origin O of the coordinate system.

E (ρ, ϕ, z) = {\begin{cases} \hat{x} \\ \hat{y} \end{cases}} E_{0} J_{0} (k_{ρ} \sqrt{ρ^{2} + ρ_{0}^{2} - 2 ρ ρ_{0} \cos (ϕ - ϕ_{0})}) e^{- i k_{z} z}

(1)

where J ₀(.) is the ordinary Bessel function and E ₀ is the electric field strength. A factor exp(iω0t) has been omitted for convenience. The field profile in Eq. (1) may be considered to be valid as long as β = sinθa/(1 + cosθa) << 1 [22

G. Milne, K. Dholakia, D. McGloin, K. Volke-Sepulveda, and P. Zemánek, “Transverse particle dynamics in a Bessel beam,” Opt. Express 15(21), 13972–13987 (2007). [CrossRef] [PubMed]

]; i.e., the results presented here will represent good descriptions of real experiments if the axicon angle is small enough that a paraxial approximation can be assumed.

In the framework of the GLMT, there are several ways in which the beam-shape coefficients associated with Eq. (1) can be computed, all of them demanding solely a previous knowledge of the radial components E_r and H_r of the beam as functions of spherical coordinates (r,θ,ϕ) relative to the origin O in Fig. 1 (in the case of optical trapping, for instance, O would coincide with the center of a spherical scatterer, as will be the case in Section 4). The radial component E_r of the electric field for this BB is readily found by putting a multiplicative factor sinθcosϕ (sinθsinϕ) to the right of Eq. (1) in the case of an x-polarized (y-polarized) beam, and replacing ρ by rsinθ and z by rcosθ. Then, for the radial component H_r of the magnetic field, we recall Faraday’s law in the frequency domain and in its differential form, so that we can write E_r and H_r , after some simple algebra, as

E_{r}_{{\begin{cases} x \\ y \end{cases}}} = E_{0} J_{0} [\sin θ_{a} \sqrt{{(k r)}^{2} \sin^{2} θ + ρ_{0}^{2} k^{2} - 2 (k r) ρ_{0} \sin θ \cos (ϕ - ϕ_{0})}] e^{- i (k r) \cos θ_{a} \cos θ} \sin θ {\begin{cases} \cos ϕ \\ \sin ϕ \end{cases}}

(2)

H_{r}_{{\begin{cases} x \\ y \end{cases}}} \approx H_{0} \cos θ_{a} J_{0} [\sin θ_{a} \sqrt{{(k r)}^{2} \sin^{2} θ + ρ_{0}^{2} k^{2} - 2 (k r) ρ_{0} \sin θ \cos (ϕ - ϕ_{0})}] e^{- i (k r) \cos θ_{a} \cos θ} \sin θ {\begin{matrix} \sin ϕ \\ - \cos ϕ \end{matrix}}

(3)

where

H_{0} = E_{0} / η = E_{0} \sqrt{ε / μ}

and η is the intrinsic impedance of the propagating medium, ε and μ being, respectively, the permittivity and permeability of this medium. The slashes represent either an x- or a y-polarized beam, with its corresponding sinϕ or cosϕ factor.

From Eqs. (2) and (3), the beam shape coefficients

g_{n, T E}^{m}

and

g_{n, T M}^{m}

in the integral localized approximation are readily found from the following relations [7

K. F. Ren, G. Gouesbet, and G. Gréhan, “Integral localized approximation in generalized lorenz-mie theory,” Appl. Opt. 37(19), 4218–4225 (1998). [CrossRef] [PubMed]

g_{n, T E}^{m} = \frac{Z_{n}^{m}}{2 π H_{0}} \int_{0}^{2 π} \hat{G} [H_{r} (r, θ, ϕ)] \exp (- i m ϕ) d ϕ

(4)

g_{n, T M}^{m} = \frac{Z_{n}^{m}}{2 π E_{0}} \int_{0}^{2 π} \hat{G} [E_{r} (r, θ, ϕ)] \exp (- i m ϕ) d ϕ

(5)

where n and m are integer numbers ranging from 1 ≤ n < ∞ and –n < m < n, respectively, corresponding to a particular spherical harmonic function

Y_{n}^{m} (θ, ϕ) = P_{n}^{m} (\cos θ) \exp (i m ϕ)

P_{n}^{m} (\cos θ)

being associated Legendre functions. In the integral localized approximation description of

g_{n, T E}^{m}

and

g_{n, T M}^{m}

, the localization operator

\hat{G}

changes the factor kr to (n + 1/2) and θ to π/2, while

Z_{n}^{m}

are normalization factors that read as [7

K. F. Ren, G. Gouesbet, and G. Gréhan, “Integral localized approximation in generalized lorenz-mie theory,” Appl. Opt. 37(19), 4218–4225 (1998). [CrossRef] [PubMed]

]

Z_{n}^{0} = \frac{2 n (n + 1) i}{2 n + 1}, m = 0

(6)

Z_{n}^{m} = {(\frac{- 2 i}{2 n + 1})}^{| m | - 1}, m \neq 0.

(7)

By substituting Eq. (2) into Eq. (5) for m = 0, we find

g_{n, T M}^{0}_{{\begin{cases} x \\ y \end{cases}}} = i \frac{2 n (n + 1)}{2 π (2 n + 1)} \int_{0}^{2 π} J_{0} [\sin θ_{a} \sqrt{{(n + 1 / 2)}^{2} + ρ_{0}^{2} k^{2} - 2 (n + 1 / 2) ρ_{0} \cos (ϕ - ϕ_{0})}] {\begin{matrix} \cos ϕ \\ \sin ϕ \end{matrix}} d ϕ

(8)

where we have used the normalization factor

Z_{n}^{0}

from Eq. (6). To analytically solve the above integral, we perform the following expansion for a zero-order Bessel function [24

G. N. Watson, A Treatise on the Theory of Bessel Functions (Cambridge, University Press, 1944), p. 358.

]

J_{0} [\sqrt{ϖ^{2} + ξ^{2} - 2 ϖ ξ \cos (ϕ - ϕ_{0})}] = \sum_{p = 0}^{\infty} ε_{p} J_{p} (ϖ) J_{p} (ξ) \cos [p (ϕ - ϕ_{0})]

(9)

ε_p being 1 if p = 0 and 2 otherwise. Thus, letting

ϖ = \sin θ_{a} (n + 1 / 2)

and

ξ = ρ_{0} k \sin θ_{a}

, we have from Eq. (8)

\begin{matrix} g_{n, T M}^{0}_{{\begin{cases} x \\ y \end{cases}}} = i \frac{2 n (n + 1)}{2 π (2 n + 1)} \int_{0}^{2 π} \sum_{p = 0}^{\infty} ε_{p} J_{p} (\sin θ_{a} (n + 1 / 2)) J_{p} (ρ_{0} k \sin θ_{a}) \cos [p (ϕ - ϕ_{0})] {\begin{matrix} \cos ϕ \\ \sin ϕ \end{matrix}} d ϕ \\ = i \frac{2 n (n + 1)}{2 π (2 n + 1)} \sum_{p = 0}^{\infty} ε_{p} J_{p} (\sin θ_{a} (n + 1 / 2)) J_{p} (ρ_{0} k \sin θ_{a}) \int_{0}^{2 π} \cos [p (ϕ - ϕ_{0})] {\begin{matrix} \cos ϕ \\ \sin ϕ \end{matrix}} d ϕ \\ = i \frac{2 n (n + 1)}{2 π (2 n + 1)} \sum_{p = 0}^{\infty} ε_{p} J_{p} (\sin θ_{a} (n + 1 / 2)) J_{p} (ρ_{0} k \sin θ_{a}) {\begin{cases} \cos p ϕ_{0} \int_{0}^{2 π} {\begin{matrix} \frac{\cos (p - 1) ϕ + \cos (p + 1) ϕ}{2} \\ 0 \end{matrix}} d ϕ + \\ \sin p ϕ_{0} \int_{0}^{2 π} {\begin{matrix} 0 \\ \frac{\cos (p - 1) ϕ - \cos (p + 1) ϕ}{2} \end{matrix}} d ϕ \end{cases}} \\ = i \frac{2 n (n + 1)}{(2 n + 1)} J_{1} (\sin θ_{a} (n + 1 / 2)) J_{1} (ρ_{0} k \sin θ_{a}) {\begin{matrix} \cos ϕ_{0} \\ \sin ϕ_{0} \end{matrix}}, \end{matrix}

(10)

which will always be zero whenever ρ ₀ = 0 (on-axis case), because

J_{1} (ρ_{0} k \sin θ_{a}) = 0

For m ≠ 0, the analysis is similar to the previous one and we have, from Eqs. (5), (7) and (9),

g_{n, T M}^{m \neq 0}_{{\begin{cases} x \\ y \end{cases}}} = \frac{1}{2 π} {(\frac{- 2 i}{2 n + 1})}^{| m | - 1} \sum_{p = 0}^{\infty} ε_{p} J_{p} (\sin θ_{a} (n + 1 / 2)) J_{p} (ρ_{0} k \sin θ_{a}) [\begin{matrix} \cos p ϕ_{0} \int_{0}^{2 π} {\begin{cases} \cos p ϕ \cos m ϕ \cos ϕ \\ - i \cos p ϕ \sin m ϕ \sin ϕ \end{cases}} d ϕ + \\ \sin p ϕ_{0} \int_{0}^{2 π} {\begin{cases} - i \sin p ϕ \sin m ϕ \cos ϕ \\ \sin p ϕ \cos m ϕ \sin ϕ \end{cases}} d ϕ \end{matrix}]

\begin{array}{l} = \frac{1}{2 π} {(\frac{- 2 i}{2 n + 1})}^{| m | - 1} \sum_{p = 0}^{\infty} ε_{p} J_{p} (\sin θ_{a} (n + 1 / 2)) J_{p} (ρ_{0} k \sin θ_{a}) [\begin{matrix} \cos p ϕ_{0} \int_{0}^{2 π} {\begin{cases} \cos p ϕ \cos | m | ϕ \cos ϕ \\ \mp i \cos p ϕ \sin | m | ϕ \sin ϕ \end{cases}} d ϕ + \\ \sin p ϕ_{0} \int_{0}^{2 π} {\begin{cases} \mp i \sin p ϕ \sin | m | ϕ \cos ϕ \\ \sin p ϕ \cos | m | ϕ \sin ϕ \end{cases}} d ϕ \end{matrix}] \\ = \frac{1}{2} {(\frac{- 2 i}{2 n + 1})}^{| m | - 1} [\begin{array}{l} {\begin{cases} 1 \\ \mp i \end{cases}} J_{| m | - 1} (\sin θ_{a} (n + 1 / 2)) J_{| m | - 1} (ρ_{0} k \sin θ_{a}) [\cos (| m | - 1) ϕ_{0} \mp i \sin (| m | - 1) ϕ_{0}] + \\ {\begin{cases} 1 \\ \pm i \end{cases}} J_{| m | + 1} (\sin θ_{a} (n + 1 / 2)) J_{| m | + 1} (ρ_{0} k \sin θ_{a}) [\cos (| m | + 1) ϕ_{0} \mp i \sin (| m | + 1) ϕ_{0}] \end{array}], \end{array}

(11)

Notice that, in deriving Eqs. (10) and (11), we have made use of some basic orthogonality conditions for trigonometric functions. Also, because J _| _m _{| ± 1}(ρ ₀ ksinθ_a ) = 0 whenever ρ ₀ = 0 and m ≠ 1, as expected for an on-axis (z-axis) beam, the only non-zero beam-shape coefficients for ρ ₀ = 0 are

g_{n, T E}^{m}

and

g_{n, T M}^{m}

, both of which with a magnitude of

(1 / 2) J_{0} (\sin θ_{a} (n + 1 / 2))

that decreases with increasing n.

Similarly, the evaluation of the BSCs

g_{n, T E}^{m}

is analogous to the previous derivation, Eq. (4) being used instead of Eq. (5). Alternatively, one may also realize the similitude between Eqs. (2) and (3) with only an interchange between sines and cosines inside the slashes, thus leading to an immediate determination of the

g_{n, T E}^{m}

g_{n, T E}^{0}_{{\begin{cases} x \\ y \end{cases}}} = i \frac{2 n (n + 1)}{(2 n + 1)} J_{1} (\sin θ_{a} (n + 1 / 2)) J_{1} (ρ_{0} k \sin θ_{a}) {\begin{matrix} \sin ϕ_{0} \\ - \cos ϕ_{0} \end{matrix}},

(12)

g_{n, T E}^{m \neq 0}_{{\begin{cases} x \\ y \end{cases}}} = \frac{1}{2} {(\frac{- 2 i}{2 n + 1})}^{| m | - 1} [\begin{array}{l} {\begin{cases} \mp i \\ - 1 \end{cases}} J_{| m | - 1} (\sin θ_{a} (n + 1 / 2)) J_{| m | - 1} (ρ_{0} k \sin θ_{a}) [\cos (| m | - 1) ϕ_{0} \mp i \sin (| m | - 1) ϕ_{0}] + \\ {\begin{cases} \pm i \\ - 1 \end{cases}} J_{| m | + 1} (\sin θ_{a} (n + 1 / 2)) J_{| m | + 1} (ρ_{0} k \sin θ_{a}) [\cos (| m | + 1) ϕ_{0} \mp i \sin (| m | + 1) ϕ_{0}] \end{array}],

(13)

These are the beam-shape coefficients necessary to fully describe an x- or y- polarized (on- or off-axis) ordinary Bessel beam in the GLMT using the integral localized approximation. They may be greatly simplified if, for example, ϕ ₀ = 0 for x-polarization, or ϕ ₀ = π/2 together with an y-polarized BB, or other analogous particular cases. Notice that we can readily extend the set of equations (10)–(13) to circularly polarized beams by means of simple symmetry relations and writing each of the new TM and TE BSC’s as functions of both TM and TE BSC’s in Eqs. (10)–(13) [15

H. Polaert, G. Gréhan, and G. Gouesbet, “Forces and torques exerted on a multilayered spherical particle by a focused Gaussian beam,” Opt. Commun. 155(1-3), 169–179 (1998). [CrossRef]

3. Numerical results for the BSC’s of an ordinary Bessel beam

In the framework of the generalized Lorenz-Mie theory, the beam shape coefficients can be primarily determined by means of finite series or using quadratures. In the latter case, double and triple integrations over the spherical coordinates (r,θ,ϕ) can be a really time-consuming task, as already pointed out by previous authors [5

–9

G. Gouesbet, “Generalized Lorenz-Mie theories, the third decade: A perspective,” J. Quant. Spectrosc. Radiat. Transf. 110(14-16), 1223–1238 (2009). [CrossRef]

,12

,13

]. For double integration, if the radial components of the electromagnetic fields of a laser beam are given, we can numerically evaluate the BSC’s

g_{n, T E}^{m}

and

g_{n, T M}^{m}

by solving the relations

g_{n, T E}^{m} = \frac{- 1}{4 π} (i^{n - 1}) \frac{R}{j_{n} (R)} \frac{(n - | m |)!}{(n + | m |)!} \int_{0}^{π} \sin θ d θ \int_{0}^{2 π} d ϕ P {}_{n}^{| m |}{(\cos θ)} \exp (- i m ϕ) \frac{H_{r} (R, θ, ϕ)}{H_{0}}

(14)

g_{n, T M}^{m} = \frac{- 1}{4 π} (i^{n - 1}) \frac{R}{j_{n} (R)} \frac{(n - | m |)!}{(n + | m |)!} \int_{0}^{π} \sin θ d θ \int_{0}^{2 π} d ϕ P {}_{n}^{| m |}{(\cos θ)} \exp (- i m ϕ) \frac{E_{r} (R, θ, ϕ)}{E_{0}}

(15)

P_{n}^{m} (\cos θ)

being the associated Legendre polynomials, R = kr an auxiliary coordinate and j_n (R) are spherical Bessel functions of integer order n. If E_r and H_r represent an exact solution to Maxwell’s equations, then R is arbitrary because, in such a case, the integrals over the spherical angles θ and ϕ are proportional to j_n (R)/R [12

,13

,25

K. R. Fen, “Diffusion des Faisceaux Feuille Laser par une Particule Sphérique et Applications aux Ecoulements Diphasiques,” Ph.D thesis (Faculté des Sciences de L’Université de Rouen, 1995).

,26

C. F. Bohren and D. R. Huffmann, Absorption and Scattering of Light by Small Particles (Wiley-Interscience, 1983).

]. But some care must be exercised in choosing a specific R when numerically evaluating Eqs. (14) and (15). It is recommended, because of the van de Hulst principle, and we have adopted it here, to impose R = n + 1/2 [11

H. C. van de Hulst, Light Scattering by Small Particles (Wiley, New York, 1957).

]. This holds because the most significant contribution of the integrands in Eqs. (14) and (15) in the case of a zero-order Bessel beam also happens around R = n + 1/2, as one can verify by plotting several curves of them and performing the same analysis as those for a focused Gaussian beam [25

]. By recalling orthogonality relations of the spherical Bessel functions, the above equations may be written as triple integrals, thus getting rid of the R dependence, with a corresponding undesirable computational cost [12

,13

To compare the numerical calculation performance of the BSCs obtained in section 2 by means of the integral localized approximation with double (Eqs. (14), (15)) and triple integration, a Fortran code was developed and is available under request. The code was then run on a PC (AMD 64 Athlon processor 3500+, 805MHz, with 768Mb of RAM) with windows XP. Table 1 compares the coefficients

g_{n, T M}^{1}

calculated via integral localized approximation (called ILA) and quadratures with double (F1) and triple (F2) integration for an x-polarized Bessel beam with wavelength λ = 1064 nm and an axicon angle θ_a = 0.0141 rad (thus, with a transverse spot of Δρ = 2.405v/(ω ₀sinθ_a ) = 28.89 μm, v being the velocity of the wave in the propagating medium), with its optical axis coincident with the z-axis, i.e., ρ ₀ = 0, in air. Only the non-zero terms calculated by the ILA method are shown (|m| = 1). Analogous results can be achieved for off-axis beams. For m < 0, the coefficients are immediately calculated by considering symmetry relations [27

K. F. Ren, G. Gréhan, and G. Gouesbet, “Symmetry relations in generalized Lorenz-Mie theory,” J. Opt. Soc. Am. A 11(6), 1812–1817 (1994). [CrossRef]

], although in our code we have preferred to calculate them directly by applying the set of Eqs. (10)-(13). This obviously requires more computational time for applications where these (m < 0) coefficients are needed, such as in optical force calculations for optical trapping.

Table 1 Beam-shape coefficients

g_{n, T M}^{1}

evaluated by quadratures (methods F1 and F2) and using the integral localized approximation ILA for an on-axis (ρ ₀ = ϕ ₀ = 0) zero-order Bessel beam with λ = 1064 nm and θ_a = 0.0141 rad ^a

n	Numerical method
n	ILA	F1	F2
1	0.499975	0.499976 + i0.962191 x 10⁻¹⁵	0.499450 + i0.117501 x 10⁻¹⁴
2	0.499900	0.499901 + i0.719626 x 10⁻¹⁵	0.497927 + i0.860414 x 10⁻¹⁵
3	0.499776	0.499773 - i0.360131 x 10⁻¹⁵	0.494984 + i0.326233 x 10⁻¹⁵
4	0.499602	0.499593 - i0.307063 x 10⁻¹⁵	0.495807 + i0.406808 x 10⁻¹⁵
5	0.499379	0.499362 + i0.160106 x 10⁻¹⁵	0.495963 + i0.351898 x 10⁻¹⁵
10	0.497518	0.497454 + i0.453892 x 10⁻¹⁵	0.491495 + i0.838867 x 10⁻¹⁴
15	0.494424	0.494301 - i0.576601 x 10⁻¹⁵	0.487111 + i0.781222 x 10⁻¹⁵
20	0.490109	0.489921 - i0.253856 x 10⁻¹⁵	0.484364 + i0.153648 x 10⁻¹⁴
50	0.439779	0.439164 - i0.162183 x 10⁻¹⁵	0.432009 - i0.869916 x 10⁻¹⁵
100	0.280726	0.279540 - i0.327642 x 10⁻¹⁵	0.266963 + i0.933362 x 10⁻¹⁶
200	−0.965629 x 10⁻¹	0.111491 - i0.995058 x 10⁻¹³	—
400	0.201049 x 10⁻¹	0.928932 x 10⁻³ + i0.336547 x 10⁻¹³	—

^a We have chosen to present numbers with six digits, according to the convention adopted by previous authors [7

K. F. Ren, G. Gouesbet, and G. Gréhan, “Integral localized approximation in generalized lorenz-mie theory,” Appl. Opt. 37(19), 4218–4225 (1998). [CrossRef] [PubMed]

The number of steps for numerically evaluating the BSCs of Table 1, using methods F1 and F2, was set to 200 for integrating both in ϕ and θ. Furthermore, for method F2, the integral over R was performed from 0 to 200 for n = 1 and up to 10, from 0 to 500 for n = 15 and 20, and from 0 to 1500 for n = 50 and 100, all of them with 500 integration points. This ensures a reasonable convergence of both quadrature methods [25

As expected,

g_{n, T M}^{m}

= 0 whenever m ≠ 1 for ρ ₀ = 0.This happens because, as already mentioned in the previous section, J _| _m _{| ± 1}(ρ ₀ ksinθ_a ) = 0. It is easily verified that, for off-axis BBs, significant values for the BSC’s appear for m ≠ 0, as also expected [12

,25

,27

K. F. Ren, G. Gréhan, and G. Gouesbet, “Symmetry relations in generalized Lorenz-Mie theory,” J. Opt. Soc. Am. A 11(6), 1812–1817 (1994). [CrossRef]

From Table 1, one can see that the magnitudes of

g_{n, T M}^{m}

(or, equivalently,

g_{n, T E}^{m}

, as

g_{n, T E}^{m} = - i g_{n, T M}^{m}

for ρ ₀ = 0) from the ILA method are very close to those obtained by means of quadratures, when we take into account only the most relevant BSCs in magnitude (low n). In fact, for all the BSCs presented, whenever the relative magnitude

g_{n, T M}^{1} / g_{1, T M}^{1}

is small there is a dissimilitude between ILA, F1 and F2 methods (for example,

{g_{100, T M}^{1} |}_{I L A}

= 0.280726,

{g_{100, T M}^{1} |}_{F 1}

≈0.279240 and

{g_{100, T M}^{1} |}_{F 2}

≈0.266963, representing a relative difference of 0.53% and 5.16%, respectively) which, overall, does not contribute significantly to the description of the field intensity profile of the Bessel beam with the parameters chosen. But we must point out that this does not represent a failure of the ILA method for high n, being solely a matter of raising the number of numerical points in θ, ϕ (for F1 and F2) and R (for F2) to ensure an adequate numerical convergence. For example, using 500 points in θ and ϕ leads to

{g_{100, T M}^{1} |}_{F 1}

≈0.280659 and

{g_{100, T M}^{1} |}_{F 2}

≈0.268462, and increasing n necessarily imposes additional refinements in our numerical integration for F1 and F2.

As we have shown that the ILA can furnish the values of the most significant BSCs with great accuracy, we now turn our attention to the most outstanding advantage of using the integral localized approximation, viz., the computation time of these coefficients. Due to the absence of numerical integrations, calculating

g_{n, T E}^{m}

and

g_{n, T M}^{m}

using the ILA is much faster than using the quadrature methods F1 and F2. Table 2 shows the average time taken for numerically evaluating the BSCs presented in Table 1. These values were obtained by calculating each ten times and then taking the time average. Notice that, as n increases, so does the average time for F1 and F2. This happens in our simulations because our subroutine uses recursive relations for the spherical Bessel functions, so that, for example, for n = 200, 200 iterations are performed to compute j ₂₀₀(R) in these methods.

Table 2 Elapsed time (in seconds) for computing the BSC’s

g_{n, T M}^{1}

of an on-axis (ρ ₀ = ϕ ₀ = 0) zero-order Bessel beam with λ = 1064 nm and θ_a = 0.0141 rad ^a

n	Numerical method
n	ILA	F1b	F2b
1	0.549968 x 10⁻³	0.512970 (0.933)	0.210406 x 10³ (0.383)
2	0.449984 x 10⁻³	0.531625 (1.181)	0.215234 x 10³ (0.478)
3	0.675186 x 10⁻³	0.546375 (0.809)	0.220765 x 10³ (0.327)
4	0.574975 x 10⁻³	0.593000 (1.031)	0.222328 x 10³ (0.387)
5	0.599977 x 10⁻³	0.593521 (0.989)	0.224562 x 10³ (0.374)
10	0.349988 x 10⁻³	0.624250 (1.784)	0.236734 x 10³ (0.676)
15	0.599975 x 10⁻³	0.702250 (1.170)	0.252047 x 10³ (0.420)
20	0.524983 x 10⁻³	0.734125 (1.398)	0.280725 x 10³ (0.535)
50	0.274985 x 10⁻³	0.968750 (3.523)	0.422434 x 10³ (1.536)
100	0.599973 x 10⁻³	0.139013 x 10¹ (2.317)	0.570141 x 10³ (0.950)
200	0.349975 x 10⁻³	0.223350 x 10¹ (6.382)	—
400	0.268730 x 10⁻³	0.398363 x 10¹ (14.824)	—

^a We have chosen to present numbers with six digits.

^b The numbers in parentheses represent the ratio (x10³ for F1 and x10⁶ for F2) between the elapsed time for the method and the corresponding one from the integral localized approximation.

Thus, the integral localized approximation proves to be a computationally efficient method for computing the beam-shape coefficients necessary to describe an ordinary BB in the framework of the generalized Lorenz-Mie theory. Although this approximation is known since de 90’s, it is noteworthy and remarkable that so few laser beams have ever been described with it [28

F. Corbin, G. Gréhan, and G. Gouesbet, “Top-hat beam technique: improvements and application to bubble measurements,” Part. Part. Sys. Characterization. 8(1-4), 222–228 (1991). [CrossRef]

–30

G. Gouesbet, “Validity of the localized approximations for arbitrary shaped beams in the generalized Lorenz-Mie theory for spheres,” J. Opt. Soc. Am. A 16(7), 1641–1650 (1999). [CrossRef]

]. It can be shown that good results are also achieved if we replace the x-polarized zero-order Bessel beam by a circularly polarized (xy plane) with the same parameters.

Now, suppose that we are given a set of beam-shape coefficients

g_{n, T E}^{m}

and

g_{n, T M}^{m}

, regardless of the numerical method used to evaluate them. The electric field components in the GLMT can be written as [13

E_{r} (R, θ, ϕ) = - i E_{0} \sum_{n = 1}^{\infty} {(- i)}^{n} (2 n + 1) \frac{j_{n} (R)}{R} \sum_{m = - n}^{n} g_{n, T M}^{m} π_{n}^{| m |} (θ) \sin θ \exp (i m ϕ)

\begin{array}{l} E_{θ} (R, θ, ϕ) = - E_{0} \sum_{n = 1}^{\infty} {(- i)}^{n} \frac{(2 n + 1)}{n (n + 1)} {\begin{cases} j_{n} (R) \sum_{m = - n}^{n} g_{n, T E}^{m} i m π_{n}^{| m |} (θ) \exp (i m ϕ) + \\ i [j_{n - 1} (R) - \frac{n}{R} j_{n} (R)] \sum_{m = - n}^{n} g_{n, T M}^{m} τ_{n}^{| m |} (θ) \exp (i m ϕ) \end{cases}} \\ E_{ϕ} (R, θ, ϕ) = - E_{0} \sum_{n = 1}^{\infty} {(- i)}^{n} \frac{(2 n + 1)}{n (n + 1)} {\begin{cases} - j_{n} (R) \sum_{m = - n}^{n} g_{n, T E}^{m} τ_{n}^{| m |} (θ) \exp (i m ϕ) + \\ i [j_{n - 1} (R) - \frac{n}{R} j_{n} (R)] \sum_{m = - n}^{n} g_{n, T M}^{m} i m π_{n}^{| m |} (θ) \exp (i m ϕ) \end{cases}} \end{array}

(16)

with similar expressions for the magnetic field components. In Eq. (16),

π_{n}^{| m |} (θ) = \sin θ P_{n}^{m} (\cos θ)

and

τ_{n}^{| m |} (θ) = (1 / \sin θ) d P_{n}^{m} (\cos θ) / d θ

. Furthermore, consider an x-polarized zero-order Bessel beam propagating along + z. The dominant electric field component E_x is written, based on the above equation, as

E_{x} (x, y, z) = E_{r} \sin θ \cos ϕ + E_{θ} \cos θ \cos ϕ - E_{ϕ} \sin ϕ

but, but, if we further assume an observation point ρ ₀ along the x-axis, we can express the electric field E_x in the same way as Ref. [13

E_{x} (x, y, z) = {\begin{cases} E_{r} (r = | x | = | ρ_{0} |, θ = \frac{π}{2}, ϕ = 0), i f x > 0 \\ - E_{r} (r = | x | = | ρ_{0} |, θ = \frac{π}{2}, ϕ = π), i f x < 0 \end{cases}

(17)

Figures 2(a)-(d) are plots of the intensity profiles of Eq. (17) for ϕ = 0 and ρ ₀ = 0, 30, 60 and 90 μm for different ranges of n and with m_max = 15. There is a clear compromise between x and the maximum n entering superposition Eq. (16) for E_r from which a good prediction of the Bessel profile is achieved. For the parameters chosen, a maximum n = 700 is sufficient, as long as the physical analysis is limited to 0 < x < 120 μm. If x is to be kept below approximately 60 μm (i.e., in cases where the physical analysis does not demand the knowledge of the field for distances above x = 60 μm), then 1 ≤ n < 500 could be used to make field calculations faster. Note that this is in contradiction with some previous analysis for focused Gaussian beams, where off-axis descriptions using the ILA may possess serious limitations due to the approximation of the beam itself, i.e., the truncation of the Davis’ series description of this type of beam [13

]. In fact, the integral localized approximation, as originally developed, is based on a Gaussian description of order L– or L [6

G. Gouesbet, G. Gréhan, and B. Maheu, “Localized interpretation to compute all the coefficients gnm in the generalized Lorenz-Mie theory,” J. Opt. Soc. Am. A 7(6), 998–1007 (1990). [CrossRef]

]. For higher ρ ₀’s, m_max may also have to assume higher values.

Fig. 2 E_x -field intensity profile for a Bessel beam with λ = 1064 nm, θ_a = 0.0141 rad and displaced (a) ρ ₀ = 0, (b) ρ ₀ = 30 μm, (c) ρ ₀ = 60 μm and (d) ρ ₀ = 90 μm along x. The accuracy of the ILA for increasing x depends on increasing the number of BSCs entering superposition (18) and (19). In all cases, m_max = 15.

Thus, for ordinary Bessel beams, the ILA is capable of reproducing the original electric field with great accuracy and reliability. We have also analyzed the intensity profile and relative difference for all the other electromagnetic field components using the ILA formulation. All of them can be reproduced by an adequate choice of n_max and m_max . Similar results were then obtained.

Finally, Fig. 3 shows the three-dimensional electric field intensity profiles for an x-polarized BB with λ = 500 nm and spot Δρ = 10.0 μm for the on-axis case, whereas Fig. 4 represents an x-polarized BB with λ = 532 nm and Δρ = 2.336 μm for an off-axis (ρ ₀ = 20 μm and ϕ ₀ = 0) situation. Excellent agreement is achieved with the ideal situation (paraxial approximation) for the transverse fields and also for the longitudinal fields when compared with vector Bessel beams presented by previous authors [22

G. Milne, K. Dholakia, D. McGloin, K. Volke-Sepulveda, and P. Zemánek, “Transverse particle dynamics in a Bessel beam,” Opt. Express 15(21), 13972–13987 (2007). [CrossRef] [PubMed]

Fig. 3 (a) 3D and (b) 2D Ey-field intensity profile for an on-axis x-polarized Bessel beam with λ = 500 nm, Δρ = 10.0 μm (θa ≈ 1.1°). The Fortran code was run by imposing nmax = 1800 (non-zero BSCs occur only for |m| = 1). (c) and (d) are the corresponding Ez-field.

Fig. 4 .a) 3D and (b) 2D Ey-field intensity profile for an off-axis (ρ0 = 5 μm and ϕ0 = π/2) x-polarized Bessel beam with λ = 532 nm, Δρ = 2.336 μm (θa = 5°, representing a limiting angle for the paraxial approximation). The Fortran code was run by imposing nmax = and mmax = 15. (c) and (d) show the corresponding Ez-field.

4. Application to optical force calculations

The success of the integral localized approximation for describing a BB can be used advantageously in biomedical applications, such as in the determination of the optical forces exerted over biological particles under the influence of such beams.

As a first example, let a dielectric spherical particle with radius a and refractive index n_p be immersed in a liquid medium with refractive index n_m , so that the relative refractive index is given by n_rel = n_p / n_m . Consider that the product k_ρa = 3.5, i.e., for a Bessel beam with λ = 1064 nm and θ_a = 0.0141 rad (thus with a transverse wave number k_ρ = ksinθ_a = (2π/λ)sinθ_a = 8.33x10⁻⁴ m⁻¹) propagating along the liquid medium, our particle has a radius a ≈42.04 μm. The beam is again assumed to be x-polarized. These parameters were chosen in such a manner as to fit those of Ref. [21

]. In the Mie regime, expressions for the x-component of the radiation pressure cross-section C_pr _, _x (which is intrinsically related and directly proportional to the x-component of the optical force) can be found elsewhere and will not be reproduced here [25

Figure 5 reveals the normalized axial force exerted on the particle when the BB is shifted along x for n_rel = 1.1. The force intensity profile for k_ρa = 0.1, i.e., a ≈1.20 μm, is also shown for comparison. In this normalized plot, F_x _,max| _kρa _{= 3.5}/ F_x _,max| _kρa _{= 0.1} = 2858.9. To reproduce these slopes, we have set n_max = 700 and m_max = 15, thus ensuring the adequate description of the Bessel beam profile, as already pointed out in previous sections.

Fig. 5 Normalized force profile exerted on a dielectric simple particle of relative refractive index nrel = 1.1. The field intensity of the zero-order Bessel beam is also shown as a solid line. All parameters of the incident beam are the same as those of Fig. 2. Positive Fx/Fmax means an attractive force towards the optical axis of the beam.

By comparing Fig. 5 with Figs. 2(a) and (b) of Ref. [21

], we see that the ILA can be successfully applied in optical force calculations, predicting the same stable equilibrium positions and magnitudes of other previous methods.

Notice, however, that according to Fig. 2, as long as the particle is assumed to be fixed at x = 0, n_max can be significantly reduced because at this point few n-iterations are necessary in calculating C_pr _, _x , to achieve an adequate description of the fields and, consequently, of the optical forces. In this way, for plotting Fig. 5, for instance, we could have safely chosen n_max = 200 and still get the same force profiles as shown.

As a second example, Fig. 6 shows the x component of the radiation pressure cross-section profile C_pr _, _x for silicon spheres (n_p = 1.4496) of four different radii immersed in water (n_m = 1.33). The incident beam is a Bessel beam with Δρ ≈2.35 μm (θ_a ≈7.51°) and a wavelength λ = 802.7 nm, both in water. These parameters were chosen to best fit the theoretical data of Ref. [22

G. Milne, K. Dholakia, D. McGloin, K. Volke-Sepulveda, and P. Zemánek, “Transverse particle dynamics in a Bessel beam,” Opt. Express 15(21), 13972–13987 (2007). [CrossRef] [PubMed]

], so that negative C_pr _, _x means an attractive force towards the optical axis of the beam.

Fig. 6 Radiation pressure cross-section C_pr,x (solid) for an x-polarized Bessel beam displaced along x using the ILA. The beam has λ = 802.7 nm and Δρ ≈2.35 μm in water (n_m = 1.33). The beam intensity is shown as a dotted line. The silicon spheres have a refractive index n_p = 1.4496 and radii a = 1.15 μm (a), 2.15 μm (b), 2.50 μm (c) and 3.42 μm (d). Points of stable equilibrium are close to those predicted in Ref. [21

], where a quadrature scheme was adopted for numerically implementing the GLMT.

Comparing our results with the above mentioned reference, we see approximately the same points of stable equilibrium. As for the simulation time, each C_pr _, _x curve took, in average, 512 s for n_max = 1300, m_max = 20 and 100 radial positions (we point out that, if we had implemented symmetry relations for calculating the BSCs for m < 0 in our code, or adopted a lower n_max , as pointed out above in deriving Fig. 5, this time would have been reduced further). These upper limits for n and m are again necessary to ensure that the set of BSCs reproduce the electromagnetic fields for this particular BB within the range 0 ≤ x ≤ 20 μm with negligible difference, thus providing the correct picture for the radiation pressure cross-section. This simulation time is radically in contrast to the elapsed time observed in Ref. [22

G. Milne, K. Dholakia, D. McGloin, K. Volke-Sepulveda, and P. Zemánek, “Transverse particle dynamics in a Bessel beam,” Opt. Express 15(21), 13972–13987 (2007). [CrossRef] [PubMed]

], where it explicitly depended on the radius of the silicon spheres and computational calculations could last more than 3.5 hours for a silicon sphere of radius a = 10 μm (for the same PC configuration, this would represent a multiplicative factor of 25.2. In practice, however, this factor is even higher, as numerical simulations of Ref. [22

G. Milne, K. Dholakia, D. McGloin, K. Volke-Sepulveda, and P. Zemánek, “Transverse particle dynamics in a Bessel beam,” Opt. Express 15(21), 13972–13987 (2007). [CrossRef] [PubMed]

]. were performed on a higher performance PC). Obviously, here it is the number of BSCs which ultimately establishes the average time, as increasing n_max and m_max turn the evaluation of Eq. (16) and C_pr _, _x more time consuming and, furthermore, more and more BSCs have to be calculated.

In fact, the use of the ILA in the generalized Lorenz-Mie theory provides an efficient tool not only for force calculations, but also for other physical entities such as the electromagnetic field components themselves, as already shown in this paper, or angular momentum, scattering and absorbing cross-sections.

5. Conclusions

For the first time in the literature, the beam-shape coefficients for ordinary Bessel beams were numerically evaluated using the integral localized approximation in the context of the generalized Lorenz-Mie theory. Closed-form solutions were presented, extending the range of applicability of the ILA to beams other than Gaussian and laser sheets.

The method presented here furnishes almost the same results of those based on quadratures for the magnitude of the BSCs necessary to adequately describe the electromagnetic fields of zero-order Bessel beams. But the fundamental advantage lies on the incredibly reduced elapsed time to compute each of these coefficients. This ultimately justifies the adoption of this method for subsequent works.

We have also applied our formalism to compute optical forces in optical trapping systems with incident ordinary Bessel beams, comparing their profiles with those already obtained by previous authors. The similitude of results, together with the lower computational time, reinforce the reliability of our approach and strongly suggest the extension of the integral localized approximation to a wider variety of laser beams, including BB of higher orders. For this particular class of multi-ringed beams, however, it is expectable that some of the symmetry relations in the original development of the ILA method will not hold true. Of course, this can certainly pose difficulties in our attempt to eliminate the angular integration presented in the formulation of the ILA for the BSCs, thus making closed-form solutions more challengeable.

Extending the range of applicability of the ILA certainly would represent a significant gain in the study of optical properties other than forces. We could use this method, for example, in optical torque calculations or simply to find absorbing, scattering and extinction cross-sections for a particle under the influence of a particular laser beam such as Laguerre-Gaussian, for example.

Acknowledgments

The authors wish to thank FAPESP—Fundação de Amparo à Pesquisa do Estado de São Paulo—under contracts 2009/54494-9 (L. A. Ambrosio’s post doctorate grant) and 2005/51689-2 (CePOF, Optics and Photonics Research Center), for supporting this work.

References and links

1.	G. Mie, “Beiträge zur Optik Trüber Medien, Speziell Kolloidaler Metallösungen,” Ann. Phys. 330(3), 377–452 (1908). [CrossRef]
2.	G. Gouesbet and G. Gréhan, “Sur la généralisation de la théorie de Lorenz-Mie,” J. Opt. (Paris) 13, 97–103 (1982).
3.	B. Maheu, G. Gouesbet, and G. Gréhan, “A concise presentation of the generalized Lorenz-Mie theory for arbitrary incident profile,” J. Opt. (Paris) 19, 59–67 (1988).
4.	G. Gouesbet, G. Gréhan, and B. Maheu, “Expressions to compute the coefficients gnm in the generalized Lorenz-Mie theory using finite series,” J. Opt. (Paris) 19, 35–48 (1988).
5.	G. Gouesbet, G. Gréhan, and B. Maheu, “On the generalized Lorenz-Mie theory: first attempt to design a localized approximation to the computation of the coefficients gnm,” J. Opt. (Paris) 20, 31–43 (1989).
6.	G. Gouesbet, G. Gréhan, and B. Maheu, “Localized interpretation to compute all the coefficients gnm in the generalized Lorenz-Mie theory,” J. Opt. Soc. Am. A 7(6), 998–1007 (1990). [CrossRef]
7.	K. F. Ren, G. Gouesbet, and G. Gréhan, “Integral localized approximation in generalized lorenz-mie theory,” Appl. Opt. 37(19), 4218–4225 (1998). [CrossRef] [PubMed]
8.	J. A. Lock and G. Gouesbet, “Generalized Lorenz-Mie theory and applications,” J. Quant. Spectrosc. Radiat. Transf. 110(11), 800–807 (2009). [CrossRef]
9.	G. Gouesbet, “Generalized Lorenz-Mie theories, the third decade: A perspective,” J. Quant. Spectrosc. Radiat. Transf. 110(14-16), 1223–1238 (2009). [CrossRef]
10.	G. Gouesbet, J. A. Lock, and G. Gréhan, “Generalized Lorenz-Mie theories and description of electromagnetic arbitrary shaped beams: localized approximations and localized beam models,” J. Quant. Spectrosc. Radiat. Transf. 112(1), 1–27 (2011). [CrossRef]
11.	H. C. van de Hulst, Light Scattering by Small Particles (Wiley, New York, 1957).
12.	G. Gouesbet and J. A. Lock, “Rigorous justification of the localized approximation to the beam-shape coefficients in generalized Lorenz-Mie theory. I. On-axis beams,” J. Opt. Soc. Am. A 11(9), 2503–2515 (1994). [CrossRef]
13.	G. Gouesbet and J. A. Lock, “Rigorous justification of the localized approximation to the beam-shape coefficients in generalized Lorenz-Mie theory. II. Off-axis beams,” J. Opt. Soc. Am. A 11(9), 2516–2525 (1994). [CrossRef]
14.	K. F. Ren, G. Gréhan, and G. Gouesbet, “Radiation pressure forces exerted on a particle arbitrarily located in a Gaussian beam by using the generalized Lorenz-Mie theory, and associated resonance effects,” Opt. Commun. 108(4-6), 343–354 (1994). [CrossRef]
15.	H. Polaert, G. Gréhan, and G. Gouesbet, “Forces and torques exerted on a multilayered spherical particle by a focused Gaussian beam,” Opt. Commun. 155(1-3), 169–179 (1998). [CrossRef]
16.	L. A. Ambrosio and H. E. Hernández-Figueroa, “Fundamentals of negative refractive index optical trapping: forces and radiation pressures exerted by focused Gaussian beams using the generalized Lorenz-Mie theory,” Biomed. Opt. Express 1(5), 1284–1301 (2010). [CrossRef] [PubMed]
17.	J. Arlt, V. Garces-Chavez, W. Sibbett, and K. Dholakia, “Optical micromanipulation using a Bessel light beam,” Opt. Commun. 197(4-6), 239–245 (2001). [CrossRef]
18.	V. Garcés-Chávez, D. McGloin, H. Melville, W. Sibbett, and K. Dholakia, “Simultaneous micromanipulation in multiple planes using a self-reconstructing light beam,” Nature 419(6903), 145–147 (2002). [CrossRef] [PubMed]
19.	V. Garcés-Chávez, D. Roskey, M. D. Summers, H. Melville, D. McGloin, E. M. Wright, and K. Dholakia, “Optical levitation in a Bessel light beam,” Appl. Phys. Lett. 85(18), 4001–4003 (2004). [CrossRef]
20.	K. C. Neuman and S. M. Block, “Optical trapping,” Rev. Sci. Instrum. 75(9), 2787–2809 (2004). [CrossRef] [PubMed]
21.	A. N. Rubinov, A. A. Afanas’ev, I. E. Ermolaev, Y. A. Kurochkin, and S. Y. Mikhnevich, “Localization of spherical particles under the action of gradient forces in the field of a zero-order Bessel beam. Rayleigh-Gans approximation,” J. Appl. Spectrosc. 70(4), 565–572 (2003). [CrossRef]
22.	G. Milne, K. Dholakia, D. McGloin, K. Volke-Sepulveda, and P. Zemánek, “Transverse particle dynamics in a Bessel beam,” Opt. Express 15(21), 13972–13987 (2007). [CrossRef] [PubMed]
23.	J. Durnin, J. J. Miceli Jr., and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58(15), 1499–1501 (1987). [CrossRef] [PubMed]
24.	G. N. Watson, A Treatise on the Theory of Bessel Functions (Cambridge, University Press, 1944), p. 358.
25.	K. R. Fen, “Diffusion des Faisceaux Feuille Laser par une Particule Sphérique et Applications aux Ecoulements Diphasiques,” Ph.D thesis (Faculté des Sciences de L’Université de Rouen, 1995).
26.	C. F. Bohren and D. R. Huffmann, Absorption and Scattering of Light by Small Particles (Wiley-Interscience, 1983).
27.	K. F. Ren, G. Gréhan, and G. Gouesbet, “Symmetry relations in generalized Lorenz-Mie theory,” J. Opt. Soc. Am. A 11(6), 1812–1817 (1994). [CrossRef]
28.	F. Corbin, G. Gréhan, and G. Gouesbet, “Top-hat beam technique: improvements and application to bubble measurements,” Part. Part. Sys. Characterization. 8(1-4), 222–228 (1991). [CrossRef]
29.	G. Gouesbet, J. A. Lock, and G. Gréhan, “Partial-wave representations of laser beams for use in light-scattering calculations,” Appl. Opt. 34(12), 2133–2143 (1995). [CrossRef] [PubMed]
30.	G. Gouesbet, “Validity of the localized approximations for arbitrary shaped beams in the generalized Lorenz-Mie theory for spheres,” J. Opt. Soc. Am. A 16(7), 1641–1650 (1999). [CrossRef]

OCIS Codes
(080.0080) Geometric optics : Geometric optics
(170.4520) Medical optics and biotechnology : Optical confinement and manipulation
(290.4020) Scattering : Mie theory
(350.4855) Other areas of optics : Optical tweezers or optical manipulation
(290.5825) Scattering : Scattering theory

ToC Category:
Optical Traps, Manipulation, and Tracking

History
Original Manuscript: April 4, 2011
Revised Manuscript: June 7, 2011
Manuscript Accepted: June 7, 2011
Published: June 9, 2011

Citation
Leonardo A. Ambrosio and Hugo E. Hernández-Figueroa, "Integral localized approximation description of ordinary Bessel beams and application to optical trapping forces," Biomed. Opt. Express 2, 1893-1906 (2011)
http://www.opticsinfobase.org/boe/abstract.cfm?URI=boe-2-7-1893

Sort: Author | Year | Journal | Reset

References

G. Mie, “Beiträge zur Optik Trüber Medien, Speziell Kolloidaler Metallösungen,” Ann. Phys. 330(3), 377–452 (1908). [CrossRef]
G. Gouesbet and G. Gréhan, “Sur la généralisation de la théorie de Lorenz-Mie,” J. Opt. (Paris) 13, 97–103 (1982).
B. Maheu, G. Gouesbet, and G. Gréhan, “A concise presentation of the generalized Lorenz-Mie theory for arbitrary incident profile,” J. Opt. (Paris) 19, 59–67 (1988).
G. Gouesbet, G. Gréhan, and B. Maheu, “Expressions to compute the coefficients gnm in the generalized Lorenz-Mie theory using finite series,” J. Opt. (Paris) 19, 35–48 (1988).
G. Gouesbet, G. Gréhan, and B. Maheu, “On the generalized Lorenz-Mie theory: first attempt to design a localized approximation to the computation of the coefficients gnm,” J. Opt. (Paris) 20, 31–43 (1989).
G. Gouesbet, G. Gréhan, and B. Maheu, “Localized interpretation to compute all the coefficients gnm in the generalized Lorenz-Mie theory,” J. Opt. Soc. Am. A 7(6), 998–1007 (1990). [CrossRef]
K. F. Ren, G. Gouesbet, and G. Gréhan, “Integral localized approximation in generalized lorenz-mie theory,” Appl. Opt. 37(19), 4218–4225 (1998). [CrossRef] [PubMed]
J. A. Lock and G. Gouesbet, “Generalized Lorenz-Mie theory and applications,” J. Quant. Spectrosc. Radiat. Transf. 110(11), 800–807 (2009). [CrossRef]
G. Gouesbet, “Generalized Lorenz-Mie theories, the third decade: A perspective,” J. Quant. Spectrosc. Radiat. Transf. 110(14-16), 1223–1238 (2009). [CrossRef]
G. Gouesbet, J. A. Lock, and G. Gréhan, “Generalized Lorenz-Mie theories and description of electromagnetic arbitrary shaped beams: localized approximations and localized beam models,” J. Quant. Spectrosc. Radiat. Transf. 112(1), 1–27 (2011). [CrossRef]
H. C. van de Hulst, Light Scattering by Small Particles (Wiley, New York, 1957).
G. Gouesbet and J. A. Lock, “Rigorous justification of the localized approximation to the beam-shape coefficients in generalized Lorenz-Mie theory. I. On-axis beams,” J. Opt. Soc. Am. A 11(9), 2503–2515 (1994). [CrossRef]
G. Gouesbet and J. A. Lock, “Rigorous justification of the localized approximation to the beam-shape coefficients in generalized Lorenz-Mie theory. II. Off-axis beams,” J. Opt. Soc. Am. A 11(9), 2516–2525 (1994). [CrossRef]
K. F. Ren, G. Gréhan, and G. Gouesbet, “Radiation pressure forces exerted on a particle arbitrarily located in a Gaussian beam by using the generalized Lorenz-Mie theory, and associated resonance effects,” Opt. Commun. 108(4-6), 343–354 (1994). [CrossRef]
H. Polaert, G. Gréhan, and G. Gouesbet, “Forces and torques exerted on a multilayered spherical particle by a focused Gaussian beam,” Opt. Commun. 155(1-3), 169–179 (1998). [CrossRef]
L. A. Ambrosio and H. E. Hernández-Figueroa, “Fundamentals of negative refractive index optical trapping: forces and radiation pressures exerted by focused Gaussian beams using the generalized Lorenz-Mie theory,” Biomed. Opt. Express 1(5), 1284–1301 (2010). [CrossRef] [PubMed]
J. Arlt, V. Garces-Chavez, W. Sibbett, and K. Dholakia, “Optical micromanipulation using a Bessel light beam,” Opt. Commun. 197(4-6), 239–245 (2001). [CrossRef]
V. Garcés-Chávez, D. McGloin, H. Melville, W. Sibbett, and K. Dholakia, “Simultaneous micromanipulation in multiple planes using a self-reconstructing light beam,” Nature 419(6903), 145–147 (2002). [CrossRef] [PubMed]
V. Garcés-Chávez, D. Roskey, M. D. Summers, H. Melville, D. McGloin, E. M. Wright, and K. Dholakia, “Optical levitation in a Bessel light beam,” Appl. Phys. Lett. 85(18), 4001–4003 (2004). [CrossRef]
K. C. Neuman and S. M. Block, “Optical trapping,” Rev. Sci. Instrum. 75(9), 2787–2809 (2004). [CrossRef] [PubMed]
A. N. Rubinov, A. A. Afanas’ev, I. E. Ermolaev, Y. A. Kurochkin, and S. Y. Mikhnevich, “Localization of spherical particles under the action of gradient forces in the field of a zero-order Bessel beam. Rayleigh-Gans approximation,” J. Appl. Spectrosc. 70(4), 565–572 (2003). [CrossRef]
G. Milne, K. Dholakia, D. McGloin, K. Volke-Sepulveda, and P. Zemánek, “Transverse particle dynamics in a Bessel beam,” Opt. Express 15(21), 13972–13987 (2007). [CrossRef] [PubMed]
J. Durnin, J. J. Miceli, and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58(15), 1499–1501 (1987). [CrossRef] [PubMed]
G. N. Watson, A Treatise on the Theory of Bessel Functions (Cambridge, University Press, 1944), p. 358.
K. R. Fen, “Diffusion des Faisceaux Feuille Laser par une Particule Sphérique et Applications aux Ecoulements Diphasiques,” Ph.D thesis (Faculté des Sciences de L’Université de Rouen, 1995).
C. F. Bohren and D. R. Huffmann, Absorption and Scattering of Light by Small Particles (Wiley-Interscience, 1983).
K. F. Ren, G. Gréhan, and G. Gouesbet, “Symmetry relations in generalized Lorenz-Mie theory,” J. Opt. Soc. Am. A 11(6), 1812–1817 (1994). [CrossRef]
F. Corbin, G. Gréhan, and G. Gouesbet, “Top-hat beam technique: improvements and application to bubble measurements,” Part. Part. Sys. Characterization. 8(1-4), 222–228 (1991). [CrossRef]
G. Gouesbet, J. A. Lock, and G. Gréhan, “Partial-wave representations of laser beams for use in light-scattering calculations,” Appl. Opt. 34(12), 2133–2143 (1995). [CrossRef] [PubMed]
G. Gouesbet, “Validity of the localized approximations for arbitrary shaped beams in the generalized Lorenz-Mie theory for spheres,” J. Opt. Soc. Am. A 16(7), 1641–1650 (1999). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper