Spherical Expansion of an Arbitrary Electromagnetic Field

Optical trapping schemes that exploit radiation forces, such as optical tweezers, are well understood and widely used to manipulate microparticles; however, these are typically effective only on short (sub-millimeter) length scales. In the past decade, manipulating micron sized objects over large distances (∼0.5 meters) using photophoretic forces has been experimentally established. Photophoresis, discovered by Ehrenhaft in the early twentieth century, is the force a small particle feels when exposed to radiation while immersed in a gas. The anisotropic heating caused by the radiation results in a net momentum transfer on one side with the surrounding gas. To date, there is no theoretical evaluation of the photophoretic force in the case of an arbitrary illumination profile (i.e. not a plane wave) incident on a dielectric sphere, starting from Maxwell’s equations. Such a treatment is needed for the case of recently published photophoretic particle manipulation schemes that utilize complicated wavefronts such as diverging Laguerre-Gaussian-Bessel beams. Here the equations needed to determine the expansion coefficients for electromagnetic fields when represented as a superposition of spherical harmonics are derived. The algorithm of Driscoll and Healy for the efficient numerical integration of functions on the 2-sphere is applied and validated with the plane wave, whose analytic expansion is known. The expansion coefficients of the incident field are related to the field inside the sphere, from which the distribution of heat deposition can be evaluated. The incident beam is also related to the scattered field, from which the scattering forces may be evaluated through the Maxwell stress tensor. In future work, these results will be combined with heat diffusion/convection simulations within the sphere and a surrounding gas to predict the total forces on the sphere, which will be compared against experimental observations that so far remain unexplained.