Tailoring Optical Gradient Force and Optical Scattering and Absorption Force

The introduction of the concept of gradient force and scattering and absorption force is an important milestone in optical trapping. However the profiles of these forces are usually unknown, even for standard setups. Here, we successfully calculated them analytically via multipole expansion and numerically via Mie theory and fast Fourier transform. The former provides physical insight, while the latter is highly accurate and efficient. A recipe to create truly conservative energy landscapes is presented. These may open up qualitatively new features in optical manipulation.

non-conservative force amounts to about 5% of the maximum transverse conservative force. While this serves as a perturbation in many cases, we remark that the energy associated with the non-conservative force can accumulate while that of conservative force cannot. An example where the non-conservative force plays a significant role is a particle trapped by a circularly polarized Gaussian beam, where the non-conservative force would rotate the particle about the beam axis. Incident light flows in one direction, this induces non-conservative scattering and absorption force that cannot be described by a potential energy approach. The subtle distinction between conservative and non-conservative forces is that conservative/non-conservative forces can/cannot be derived from a potential energy and the associated work done is path independent/dependent. Consequently, they are of different attributes and have different applications. Generally speaking, F g and F k are, respectively, responsible for optical trapping 12 and particle transportation 40 . They can also be combined to achieve other functionalities. For a conservative force, the textbook conservative classical mechanics and equilibrium statistical mechanics may be applied, which is in general significantly simpler than their non-conservative analog 43,44 . Furthermore, a colloid in a conservative force field can also be used to simulate other systems where experiments are more difficult 45,46 .
Here, we present an analytical and a numerical approach to calculate these forces. With these tools, we created a recipe to produce a fairly general class of conservative optical force field characterized by F k = 0. In general, particles immersed in an optical force field do not obey equilibrium statistical mechanics, making the analysis complicated 43,44 . With conservative forces, these issues are resolved. Such conservative optical force field is required in many applications of optical micromanipulation.

Results
Analytical Method. Based on a previously derived multipolar expression for optical forces acting on a spherical particle 19,36,47 , we derived in the supplemental material the analytical expression of F g and F k for the first few leading multipoles: where α′, β′, γ′, γ′ m , Ω′ and α″ are the multipole moments obtainable from Mie theory and are tabulated in the supplemental material, k is the wavenumber, a is the particle radius, and E in and B in are the arbitrary incident electromagnetic fields. In principle, one may keep adding higher order terms (such as octopole moment and beyond) into the multipole expansion to derive more accurate expression, but in practice the mathematics can be prohibitive, especially when Toroidal moments are involved. We note that Eq. (1) goes beyond the previous dipolar theory and reaches into the Mie regime.
Numerical Method. Nevertheless, optically trapped particles often have sizes beyond the validity of Eq. (1).
The real advantage of Eq. (1) lies in its transparent physics and insight, rather than its ability to compute the numerical values of F g and F k . Here, we adopted an independent efficient numerical method based on fast Fourier transform (FFT) to treat particles with arbitrary sizes:  2)) approaches is demonstrated. Gradient force (red) and scattering and absorption force (blue) calculated by the analytical expression Eq. (1) (symbols) and numerical approach Eq. (2) (solid lines), respectively. The 300 nm diameter particle of refractive index 1.59 is immersed in water (refractive index = 1.33). The wavelength is 1064 nm. 3 x is the Fourier transform of the total optical force F(x). We remark that the use of FFT in Eq. (2) significantly increases the computational speed and also the wavelet analysis can be an alternative to the FFT 48 . Eq. (2) can be applied to treat sufficiently fast decaying fields or periodic fields. The total optical force needed in the FFT is computed by particle surface is the time averaged Maxwell stress tensor, with the required E and H calculated by the generalized Mie theory [22][23][24]26 , which at least within classical electrodynamics, is exact. It can be readily verified that in both Eqs (1) and (2), F = F g + F k while ∇ × F g = 0 and ∇⋅F k = 0, so they are indeed the gradient force and the scattering and absorption force. Figure 1 plotted F g and F k for a dielectric spherical particle with a diameter of 300 nm illuminated by an x-polarized fundamental Gaussian beam with a wavelength of 1064 nm, calculated analytically using Eq. (1) (symbols) and numerically using Eq. (2) (solid lines). The strongly focused Gaussian beam, which is known as optical tweezers, is modeled by using the highly accurate generalized vector Debye integral 31,49,50 , which is known to generate results that can be directly compared with experiments for particles larger than the wavelength. We remark that for particles smaller than wavelength, one must take into account the astigmatism to obtain an accurate result 51 . In short, the incident unfocused laser beam is illuminated on an objective lens with high numerical aperture. Since the lens is macroscopic in size, the focusing of light can be treated using geometrical optics with negligible errors. Then the vector Debye integral maps the geometrical optics solution to the field in the focal region. As far as the 300nm-diameter dielectric particle is concerned, excellent agreement is achieved. This validates both analytical and numerical approaches.
In Fig. 2, we plotted F g and F k using the numerical approach (Eq. (2)), for the widely employed linearly polarized fundamental Gaussian beam. To obtain converged calculation for the strongly focused Gaussian beam, the unit cell for FFT is chosen such that the forces near the edges are at least two to three orders of magnitude smaller than that of the center. We also repeated our calculations with different unit cell sizes. For sufficiently large unit cell sizes, the results converge very well.

Generation of a Conservative Force Field.
We now return to the generation of a conservative force field characterized by , which is equivalent to having an incident standing wave. In other words, a standing wave can generate a conservative force field in three dimensions for particles with diameter less than roughly half a wavelength. Figure 3(a-c), plotted, respectively, the potential energy U (where F g = −∇U), |F g |, and |F k | for a 1-micron diameter particle illuminated by a standing wave generated by interfering four plane waves (each with an intensity of 10 4 W/cm 2 ). Clearly, a conservative periodic potential is generated, since |F k | in Fig. 3(c) is very small, on the same order as the numerical noise. We note that a 1-micron diameter particle is beyond the validity of Eq. (1), and therefore the observed conservative force is not a prediction of Eq. (1). We repeated the calculation with 5 microns diameter particle (data not shown), still F k is comparable with the numerical noise. These findings are consequences of a more general theorem we analytically proved in the supplemental material: for a spherical particle illuminated by TE or TM standing wave in two dimensions (i.e. all incident wave vectors lie on the z = 0 plane), the optical force is conservative, just as shown in Fig. 3. We remark that the vanishing of the scattering and absorption force is a consequence of symmetry (see supplemental material). We can generate a conservative force even for a lossy, absorptive particle (data not shown), according to our analytical proof in the supplemental material.
For comparison, the forces acting on a 1-micron diameter particle when illuminated by three plane waves are plotted in Fig. 4. These three plane waves do not form a standing wave, but their incident momentum do cancel each other completely, i.e. they have the same amplitude and ∑ = = k 0 . When the waves are coherent, both F k and F g are non-zero due to interference. The force field is clearly non-conservative, as the maximum value of F k is actually greater than that of F g . This highlights the importance of having a standing wave.

Discussion
In summary, we devised an analytical and a numerical approach to calculate the gradient force and the scattering and absorption force for the experimentally accessible micro-particles. The profile of these forces associated with the widely employed optical tweezers is presented. This will enable more detailed analysis and precise control on optical micromanipulation. As an example, we provided a recipe and a sufficient condition to generate or tailor a conservative force field for particles smaller than half a wavelength using arbitrary standing wave. For TE or TM incident standing wave in two dimensions, the induced force field is conservative for spherical particle of any size. This will allow us to create a truly conservative force field, paving the way to mimic a wide variety of phenomena in equilibrium statistical mechanics using optical micromanipulation system. Finally, we remark that our approach can also be used in other non-conservative force field, such as the acoustic force.