On optical forces in spherical whispering gallery mode resonators
In this paper we discuss the force exerted by the field of an optical cavity on a polarizable dipole. We show that the modification of the cavity modes due to interaction with the dipole significantly alters the properties of the force. In particular, all components of the force are found to be non-conservative, and cannot, therefore, be derived from a potential energy. We also suggest a simple generalization of the standard formulas for the optical force on the dipole, which reproduces the results of calculations based on the Maxwell stress tensor.
Department of Physics, Queens College of the City University of New York, Flushing, NY 11367, USA & The Graduate Center of CUNY, 365 5th Avenue, New York, NY 10016, USA \email@example.com
(350.4855) Optical tweezers or optical manipulation, (140.3945) Microcavities (140.7010) Laser trapping
Appendix A Introduction
Starting with the pioneering works of Ashkin [1, 2], where trapping and manipulation of objects by optical forces was first demonstrated, there has been an explosion of interest in optical forcing of particles. Optical tweezers have since been developed into a standard tool used in many applications, while on the fundamental side optical cooling of atoms has opened unique opportunities for exploring various quantum mechanical many-body phenomena. Recently a great deal of interest has been devoted to the possibility of optical cooling of macroscopic objects such as mechanical nanoresonators [3, 4, 5, 6], membranes [7, 8, 9] or particles oscillating in optical traps [10, 11, 12, 13, 14, 15]. While traditional optical trapping experiments involve freely propagating laser beams, optical microcavities, which confine light in a small volume, have emerged as a candidate for the source of optical forces [16, 17]. The spatial confinement of light within the small mode volume of the cavity results in an increase of the strength of the force. Additionally, a material object interacts with the cavity mode, resulting in a shifted modal frequency and altered field distribution. As a result, the dynamics of mechanical degrees of freedom become coupled with those of the field giving rise to so-called dynamical backaction  responsible for effects such as cavity cooling or heating of mechanical degrees of freedom [5, 19, 12].
Understanding cavity optomechanical phenomena depends on correct representation of the optical force exerted by the cavity modes. In the case of free-propagating optical fields (i.e. laser beams), the force on a subwavelength object (dipole) is naturally separated into gradient and scattering components . The gradient component is analogous to the force on a static dipole, which tends to draw a particle into regions of greater field intensity. It can be presented as the gradient of the electromagnetic energy of the polarized particle and is therefore conservative. The scattering component results from radiation pressure and is expressed in terms of the momentum flux impinging on the particle per unit time. This force is non-conservative because it results from the process of irreversible exchange of momentum and energy between the particle and optical field. Due to its conceptual simplicity and apparent universality, this paradigm has become firmly engrained in the current literature and has been accepted as a framework for calculating optical forces also due to cavity modes. The effects of spatial confinement are taken into account by using cavity modes to represent electromagnetic field while allowing the resonant frequencies of the cavity to depend on mechanical degrees of freedom.
However, in our recent paper  we pointed out that the validity of the conservative/gradient, non-conservative/scattering paradigm explicitly depends on the assumption that the particle itself does not change sources of the incident field. While this condition is fulfilled for free-propagating electromagnetic fields, it is violated for fields confined within a cavity. The position-dependence of the cavity resonance frequencies is just one manifestation of this phenomenon. In order to elucidate all consequences of the particle-induced modification of the cavity field, we considered the interaction between whispering-gallery-modes of a spherical microresonator and a small dielectric particle. In this system forces can be calculated by a rigorous analytical approach based on the Maxwell stress tensor. These calculations show that the force cannot be described within the standard paradigm . In particular, no vector component of the optical force is conservative, i.e. derivable as a gradient of potential energy. Furthermore, the force tangential to the surface of the cavity, which is responsible for the ”carousel” effect observed in Ref. , was found to have a contribution proportional to the real part of the particle’s polarizability. This contradicts the assertion of Ref.  that the tangential optical force, which is in the direction of the momentum flux of the unmodified by the particle cavity mode, is of strictly ”scattering” origin. In Ref.  these results were obtained using formal ab initio approach, which while providing accurate expressions for the force, does not allow for the simple physical interpretation of the results and for finding connections between them and heuristic approaches used in works of other authors. The objective of this paper is to show how the traditional gradient-scattering force approach can be generalized to derive the results of Ref.  without relying on complicated Maxwell stress tensor calculations. It can be conjectured that the ”pseudo-gradient” formalism of this paper can be applied even in situations in which rigorous treatment is not possible. In addition, the results presented here allow elucidating limitations of previous heuristic approaches to optical forces.
To achieve this objective in the most efficient way, we organize this paper in the following manner. We begin by reviewing the derivations of the gradient and scattering forces paying particular attention to the assumptions involved. Based on this discussion we propose a psuedo-gradient procedure as a way to extend standard gradient/scattering approach to optical cavities. We use this idea to determine the force on a particle due to a spherical whispering gallery mode resonator and compare it with the results of stress tensor based calculations.
Appendix B Optical force on a small polarizable particle
b.1 Gradient force: thermodynamic derivation
We begin by recalling the thermodynamic approach to deriving the electrostatic force on a small dielectric particle. This approach emphasizes one of the key requirements for the validity of the gradient paradigm. Following the classical textbook by Landau & Lifshitz  the electrostatic component of the free energy of the polarizable particle is
where and the integral runs over the volume of the particle. here refers to the total field in and around the particle, and is its polarization. This total energy contains the energy, , of the external field, , which would have existed in the absence of the particle. Subtracting , the energy required to polarize the particle is
Using this expression one can derive the change in due to a small variation of as
The force is determined by substituting in the change in the field, , due to an infinitesimal translation of the particle, , which is given by: . Equation 2 then takes the following form:
Taking into account that this change in energy can be related to the work of the force acting on the particle: , one obtains expression for the force as
If is approximately constant over the dimensions of the particle (dipole approximation), it can be taken out of the integral, giving
where the dipole moment dV. Assuming linear polarizability of the particle, , one derives the final gradient expression for the force:
It should be noted, however, that this derivation depends critically on the assumption that a displacement of the particle does not affect the distribution of the sources of the external field . Without this assumption one would not be able to equate the change of due to small variations of the sources (Eq. (2)) with its change due to particle’s displacement making the subsequent equations invalid.
b.2 Gradient force: direct derivation
An alternative derivation, which is also commonly encountered in the textbooks, is based on a model of an electric dipole as a system of equal and opposite charges , separated by some small, ultimately infinitesimal distance . The same approach can be used to describe forces on an induced dipole characterized by polarizability . The dipole is assumed to be placed in some external field , so that . We make no demands on other than that be defined self consistently with it. In particular, may be dependent on the dipole itself. For example, if the external field is due to charges on a conductor, the presence of the dipole will alter the charge distribution and thus . We make this explicit by writing , where and are respectively a field point and the position vector of the particle. The total force is derived by considering the Coulomb forces at each charge comprising the dipole: (see Fig. 1). By taking the limit , keeping constant, the force derived in this case is:
where refers to a gradient with respect to field coordinates . While this result looks similar to Eq. (5), there is an important difference between them. The force given by Eq. (5) implies that before taking the spatial derivative of the field with respect to field coordinates , one sets coordinate of the particle to coincide with . On the other hand, Eq. (7) requires that this procedure is reversed. Physically, it reflects the fact that the electric force on a dipole results from the spatial variation of the electric field across it. The two equations, Eq. (5) and Eq. (7), produce identical results only if the field exerting the force does not depend upon particle’s position.
Equation (7) can also be transformed into a ”gradient” form
which, however, differs from Eq. (6). Unlike the latter, Eq. (8) involves taking the gradient of the function of two variables, and, therefore, the expression cannot be interpreted as a potential energy unless is independent of . Correspondingly, the force calculated according to Eq. (6) does not have to be conservative.
To simplify terminology and notations in the subsequent consideration we will call the operation presented in Eq. (8) a ”pseudo-gradient” and will use notation to represent it.
b.3 Total force on a dipole
In this subsection we generalize the previous results to the case of a dipole interacting with a time-dependent harmonic electromagnetic field. The oscillating dipole moment creates a current density , giving rise to a magnetic force. The total Lorentz force on the particle for incident fields independent of the dipole’s position have the standard form . The time averaged expression for this force can be rewritten as :
where now is a complex valued (with radiative corrections included) polarizability of the dipole. Based upon analysis of the previous sub-section we conjecture that the expression for the force in the case of the field dependent on the particle’s position can be obtained from Eq. (9) by substituting with respectively. For a standard dipole particle with radius , refractive index , and polarizability
Equation (9) for the force can be re-written as
where is the average polarization energy of the dipole , is its scattering cross section , and is the average momentum density of the field , and is the speed of light in vacuum.
In the subsequent sections of the paper we will apply Eq. (10) to the case of a small dielectric particle interacting with a whispering gallery mode (WGM) of a spherical resonator. Comparing the obtained expressions for the force with results of Ref. , where the optical force in this system was calculated using Maxwell stress tensor approach, we will be able to substantiate validity of Eq. (10) and shed additional light on physical properties of the optical forces due to cavity-confined electromagnetic field.
Appendix C Optical force of a WGM resonator
c.1 WGMs of a single spherical resonator
An optical whispering gallery mode is a long living excitation of a spherical resonator, which can be thought of as a ray of light propagating along the equator of the sphere and trapped in it due to total internal reflection. In a spherical coordinate system with polar axis perpendicular to the plane of the propagation of the mode ( system in Fig. 2), its field is described by a single vector spherical harmonic (VSH) characterized by polar, azimuthal, and radial indexes respectively and polarization, TE or TM. For concreteness we focus only on TE polarized modes with radial index defined as (time harmonic factor is assumed and suppressed). Here is a normalization factor, is the scattering amplitude describing response of the resonator to an incident radiation with frequency , is the spherical Hankel function of the first kind, and is the magnitude of the wavevector . In the close vicinity of a chosen resonant frequency the scattering amplitude can be approximated as
with and being respectively the resonant frequency and decay rate of the mode. The angular portion of the VSH is defined as , where is the dimensionless angular momentum operator , and is the scalar spherical harmonic. We focus here on so called fundamental modes with . WGMs with long life-time are characterized by and in expressions that follow we neglect terms of order compared to those or order unity.
For calculation of optical forces due to this WGM it is also convenient to find an expression for its field in a coordinate system with polar axis connecting the center of the sphere and the point of observation ( system in Fig. 2.) Such an expression can be obtained using rotational properties of VSH  expressed as
where is the Wigner D function and are the Euler angles specifying the rotation from the unprimed to primed coordinate system. The D-functions are defined by
where the function for can be written as
In order to find an expression for the WGM in the primed system one has to apply inverse transformation to the VSH defined in the unprimed coordinate system. Rotation by the angle is equivalent to shifting the coordinate, so we only consider transformations whith . Applying this transformation to a WGM with orbital number we find its representation in the rotated system in the following form:
Here we introduced dimensionless variable representing the relative detuning of the external frequency from the resonance of the WGM with respect to its width. At any point on the polar axis (), all vanish except for . Assuming , and dispensing with the prime on , the field can be written explicitly as:
where and and are the spherical coordinate unit vectors referred to the global, unprimed system. The magnetic field is given by:
where the prime denotes differentiation with respect to argument . We can express the field at any point in space in the form of Eq. (16) and Eq. (17) by changing the Euler angles appearing in . It is important to note that angles and correspond to the respective angular coordinates and of the axis as viewed from the unprimed coordinate system.
c.2 Calculation of the force neglecting particle-resonator coupling
In order to elucidate the effects of the particle-induced modification of WGMs on the properties of the optical force, we first compute this force with this modification neglected. In this case, assuming that a particle with coordinates as defined in the coordinate system, lies on the axis of the system, we can substitute the results of the previous sub-section, Eq. (16) and (17), into Eq. (10) for the force, while replacing the operator of pseudo-gradient with the regular gradient. In the limit function can be approximated as
where we took into account that . Consequently, the field coefficients for can be presented in terms of the coefficient as . It is clear, therefore, that products of the form in this case are purely real, so that any terms in the force proportional to vanish. Since the field strength decays quickly when the particle moves out of the equatorial plane , we only consider , and expand in terms of . However, since for , we shall only expand the terms, which do not contain . Keeping linear terms in the expansion of , and taking into account that for the imaginary part of the Hankel function is much greater than its real part in the region , we obtain the gradient portion of the optical force :
and is the spherical Neumann function.
The leading term in is of the second order in and can be neglected compared to the scattering force, , which takes the form:
where . Thus, in this approximation, the gradient force draws the particle toward the resonator radially, and maintains it in the equatorial plane via its component. It is important to note that its dependence on the particle’s position follows the behavior of the Neumann function , which is almost exponential in the considered range of parameters. At the same time, the azimuthal component of the force is of purely scattering nature and proportional to the moment density, which can be presented in the simple form , where is the electromagnetic energy density (the use of symbol for the field emphasizes that it is not affected by the presence of the particle).
c.3 Effect of the particle-induced modification of the WGM on the optical force
Modification of the WGM by the particle
The problem of determination of the electromagnetic field of the coupled resonator-dipole system is analytically tractable and was solved in Ref. [26, 27]. The dipole is modeled as a small sphere with radius , where , and refractive index , and the field is found in the form of a general expansion in terms VSHs with all and polarizations. The particle is found to modify the WGM in two significant ways. First, it creates an additional resonance at frequency with width in addition to the resonance of a single sphere, with and depending only on . Second, the steady state of the resonator field associated with the resonance is significantly modified compared to the field distribution of the initial WGM. Initially isotropic field turns into a highly directional distribution oriented predominantly toward the particle. Thus, a displacement of the particle in the direction causes the resonator’s field to move with the particle (see for details Ref. [26, 27]). In addition, the interaction with the particle excites in the resonator WGMs with different and polarization. However, these contributions are small, and can be neglected. In this approximation, the scattered field of the resonator can again be presented in the form of Eq. (14), but with expansion coefficients, which are no longer given by Eq. (15). They have the following form
where is the VSH translation coefficient , which arises when the field scattered by one sphere is expressed in terms of VSH centered about the other, and is given by
Calculation of the force with particle-modified field
In this part of the paper we assess the effects of the particle-induced shift of the resonance frequency and of the changes in the spatial configuration of the field of the resonator on the optical forces exerted by it. To this end we shall analyze the expressions for the force obtained by evaluating Eq. (10) with the field at the location of the particle given by Eq. (14). The role of the pseudo-gradient operator in this equation is to distinguish between field coordinates and particle coordinates even though we calculate the force at the point . Taking into account that dependence on is only contained in the expansion coefficients , this procedure becomes rather trivial: one needs to find the gradient of all respective expressions treating these coefficients as constants, and after that equate . Calculating the required gradients we obtain for the pseudo-gradient component of the force:
The scattering component of the force, , takes in this case the form of
A contribution from the remaining term in Eq. (10), which is proportional to remains negligible and will not be considered any further.
One can see that Eq. (24) and Eq. (27) differ significantly from the respective Eq. (19) and Eq. (21) obtained under the assumption of the unmodified WGM. Further analysis of the obtained expression will be performed in two steps. Since it is often assumed that the main effect of the particle on the cavity mode consists in changing the resonance frequency, we first separate this effect. To achieve this we allow the particle to shift and broaden the resonance according to Eq. (23), but will assume that the field coefficients are given by the unmodified Eq. (15) with replacement of by . In this case the particle can modify the amplitude of the resonator’s field, but does not change its spatial distribution. Then, one immediately sees that in Eq. (24), , which contains terms proportional to with vanishes. For the same reason the -component of the scattering force also vanishes.
In order to present the radial component in the form easily related to expressions used in papers of other authors such as those of Ref. [10, 13, 28], we introduce the power radiated by those modes of the resonator that interact with the particle as,
and the number of photons in them as
The expression for the radial component of the force can be written down now as
which should be contrasted with the expression obtained from the usual derivative of interaction energy
The last term in this expression is spurious as was first noticed in Ref. , where it was shown by direct numerical simulations that the force on a dipole in Fabry-Perot resonator must be given by Eq. (28) rather than by Eq. (29).
This result is a clear demonstration of the fact that the polarization energy of the particle cannot be considered a true potential energy even if one neglects the spatial modification of the cavity mode due its interaction with the particle. However, taking this modification into account results even in more drastic changes in the optical force yielding a non-zero azimuthal component of the pseudo-gradient force and a non-zero polar component of the scattering force. The expression for the total force in this case is found by using the correct set of the field coefficients as defined by Eq. (22). The radial component of the force does not change from Eq. (28), while its polar and azimuthal components now can be presented as
Terms proportional to come from the psuedo-gradient force, while terms come from the scattering force. These expressions demonstrate significant deviation of the force from both completely unmodified and frequency-only modified WGM approximations. First, let us note that the radial dependence of the force is determined by the factor in addition to the Neumann function in and . The role of this factor can be seen as follows. The condition is satisfied for some at which driving frequency coincides with the particle-induced resonance. If one linearizes about this point as , where
If is sufficiently small, then in the region where can be considered constant, the factor in is also constant. Therefore the spatial profile of the force has Lorentzian shape peaked at with width . Let us recall that in the unmodified WGM approximation the magnitude of the force monotonically (essentially exponentially) decreases with .
The azimuthal force is no longer solely due to the scattering contribution. Two different limits are of interest based upon choice of the external driving frequency . In the limit , . The magnitude of the pseudo-gradient term exceeds that of the scattering term unless . This can only happen for very small values of , when the magnitude of the force is also very small. When , the pseudo-gradient contribution exceeds the scattering force, and the tangential component can be written as , where is the scattering force in the unmodified WGM approximation as given by Eq. (21). If satisfies , the tangential force exceeds . When the scattering contribution to becomes negligible as well. In this case the force can be written , which also exceeds the magnitude that would have been obtained in the unmodified WGM approximation. These results show that the force propelling the particle in the experiments like the one of Ref.  is not necessarily of scattering origin and might have a pseudo-gradient contribution. The two can be distinguished by their dependence on : while the pseudo-gradient force is linear in this parameter, the scattering force is quadratic.
Even in the range of parameters where the pseudo-gradient contribution to the azimuthal component of the force dominates, it remains non-conservative since it imparts net kinetic energy to the particle moving along a closed orbit around the resonator. This occurs because there is a field gradient which pushes the particle in the direction. When the particle moves to a new point, the field re-adjusts so that there is again a field gradient in the direction. Implicit in this analysis is the assumption that the particle moves slowly enough to consider the field always remaining in the quasi-steady state. Velocity dependent effects can become significant when the time scale of particle motion is comparable to the relaxation time of the resonator, .
The relative magnitudes of the force components can be analyzed by comparing to . From the asymptotic expansions for the spherical Nuemann functions and their derivatives in the region , , we have  where prime denotes differentiation with respect to argument and is defined by . Since resonances are in the region , where is the refractive index of the resonator, we have and thus . Thus, assuming we are near a particle induced resonance so that , the relative magnitude of the forces will be determined by the factor in and . If the system is driven at a frequency close to the ideal Mie frequency, so that is of order unity, then the radial and azimuthal forces will be of comparable magnitude, while the polar force will be smaller by a factor . If on the other hand, then the azimuthal force will be smaller than the radial force by a factor , while the polar force is smaller by a factor .
It can also be seen that the scattering contribution to the azimuthal and polar forces is in general smaller than the pseudo-gradient contribution due to the fact that (where ). In the limit where both the driving frequency and particle induced resonance frequencies are very close to the ideal Mie resonance, so that , the scattering contribution to becomes appreciable, while it vanishes in . This is to be expected given that is the limit where the particle induced modification of the cavity mode becomes vanishingly small, and accordingly the behavior of the forces approaches that of their unmodified forms of Eq. (19) and Eq. (21). This is likely the regime encountered in the experiments of Ref. .
The results of the calculation of the force within the pseudo-gradient approach can be compared with calculations carried out by integrating the Maxwell stress tensor over a surface of the particle, which is assumed to have a small, but finite size. For a field represented by a VSH expansion, the stress tensor integral over a spherical region can be performed analytically and the force given in terms of the VSH expansion coefficients . For the present case these coefficients are given in Ref.  while full details of the calculations of the force can be found in Ref. . In the large limit, the forces obtained agree exactly with those calculated from the pseudo-gradient approach validating the latter.
It is interesting to note that the large limit of the stress tensor calculations is necessary to maintain consistency with the assumed point-like nature of the particle in the pseudo-gradient approach. To see this, note that for a given resonator of radius and refractive index , the lowest order approximation to the resonant frequency is the geometric optical condition . At the same time, a point dipole is defined by the the limit with electromagnetic size parameter kept constant. Combining these two conditions we have , which implies that taking requires that .
Appendix D Conclusion
We have presented here a generalization of the theory of optical forces on a dipole for the case when it interacts with the electromagnetic field of an optical cavity. The traditional gradient/scattering paradigm is shown to be invalid when the dipole can modify the source of the field. In particular, all vector components of the force are found to be non-conservative and, consequently, no component can be derived from a gradient of electromagnetic polarization energy. We have further shown that, when the particle-induced modification of the resonator field is taken into account, the force in the direction of the energy and momentum flux of the wave cannot be interpreted as a ’scattering’ force. In place of the gradient/scattering paradigm we have proposed a pseudo-gradient framework which is conceptually simpler and computationally more efficient than the exact Maxwell stress tensor approach. Using an example of a small dielectric particle interacting with whispering-gallery-modes of a spherical resonator we demonstrated that the suggested pseudo-gradient formalism reproduces all results of the calculations based on the Maxwell stress tensor.
The results of this work have important implications for the quantum theory of optomechanical interaction, which is commonly based on the assumed potential nature of the gradient force. These results are also of importance for proposed optofluidic sensors which rely on a tangential force to drive the particle in orbit around the resonator.
- A. Ashkin, “Acceleration and trapping of particles by radiation pressure,” Phys. Rev. Lett. 24, 156–159 (1970).
- A. Ashkin, J. M. Dziedzic, J. E. Bjorkholm, and S. Chu, “Observation of a single-beam gradient force optical trap for dielectric particles,” Opt. Lett. 11, 288–290 (1986).
- G. Anetsberger, O. Arcizet, Q. P. Unterreithmeier, R. Riviere, A. Schliesser, E. M. Weig, J. P. Kotthaus, and T. J. Kippenberg, “Near-field cavity optomechanics with nanomechanical oscillators,” Nature Phys. 5, 909–914 (2009).
- A. Schliesser, O. Arcizet, R. Riviere, G. Anetsberger, and T. J. Kippenberg, “Resolved-sideband cooling and position measurement of a micromechanical oscillator close to the heisenberg uncertainty limit,” Nature Phys. 5, 509–514 (2009).
- A. Schliesser, R. Riviere, G. Anetsberger, O. Arcizet, and T. J. Kippenberg, “Resolved-sideband cooling of a micromechanical oscillator,” Nature Phys. 4, 415–419 (2008).
- S. Groeblacher, K. Hammerer, M. R. Vanner, and M. Aspelmeyer, “Observation of strong coupling between a micromechanical resonator and an optical cavity field,” Nature 460, 724–727 (2009).
- O. Arcizet, C. Molinelli, T. Briant, P.-F. Cohadon, A. Heidmann, J.-M. Mackowski, C. Michel, L. Pinard, O. Francais, and L. Rousseau, “Experimental optomechanics with silicon micromirrors,” N. J. Phys. 10 (2008).
- A. M. Jayich, J. C. Sankey, B. M. Zwickl, C. Yang, J. D. Thompson, S. M. Girvin, A. A. Clerk, F. Marquardt, and J. G. E. Harris, “Dispersive optomechanics: a membrane inside a cavity,” N. J. Phys. 10 (2008).
- D. J. Wilson, C. A. Regal, S. B. Papp, and H. J. Kimble, “Cavity optomechanics with stoichiometric sin films,” Phys. Rev. Let. 103 (2009).
- D. E. Chang, C. A. Regal, S. B. Papp, D. J. Wilson, J. Ye, O. Painter, H. J. Kimble, and P. Zoller, “Cavity opto-mechanics using an optically levitated nanosphere,” Proc. Natl. Acad. Scie. 107, 1005–1010 (2010).
- O. Romero-Isart, and M. L. Juan, and R. Quidant, and J. I. Cirac, “Toward quantum superposition of living organisms,” N. J. Phys. 12, 033015 (2010).
- P. F. Barker and M. N. Shneider, “Cavity cooling of an optically trapped nanoparticle,” Phys. Rev. A 81, 023826 (2010).
- O. Romero-Isart, A. C. Pflanzer, M. L. Juan, R. Quidant, N. Kiesel, M. Aspelmeyer, and J. I. Cirac, “Optically levitating dielectrics in the quantum regime: Theory and protocols,” Phys. Rev. A 83, 013803 (2011).
- Z.-q. Yin, T. Li, and M. Feng, “Three-dimensional cooling and detection of a nanosphere with a single cavity,” Phys. Rev. A 83, 013816 (2011).
- T. Li, S. Kheifets, and M. G. Raizen, “Millikelvin cooling of an optically trapped microsphere in vacuum,” Nature Phys. 7, 527 – 530 (2011).
- T. J. Kippenberg and K. J. Vahala, “Cavity optomechanics: Back-action at the mesoscale,” Science 321, 1172–1176 (2008).
- A. Schliesser and T. J. Kippenberg, “Cavity optomechanics with whispering-gallery mode optical micro-resonators,” in Advances In Atomic Molecular and Optical Physics 58, (2010), 207–323.
- V. Braginsky and A. Manukin, Measurment of Weak Forces in Physics Experiments (University of Chicago Press, Chicago, 1977).
- R. J. Schulze, C. Genes, and H. Ritsch, “Optomechanical approach to cooling of small polarizable particles in a strongly pumped ring cavity,” Phys. Rev. A 81 (2010).
- M. Nieto-Vesperinas, P. Chaumet, and A. Rahmani, “Near-field photonic forces,” Phil. Trans. R. Soc. Lond. A 362, 719–737 (2004).
- J. Rubin and L. Deych, “Optical forces due to spherical microresonators and their manifestation in optically induced orbital motion of nanoparticles,” Phys. Rev. A 84, 023844 (2011).
- S. Arnold, D. Keng, S. I. Shopova, S. Holler, W. Zurawsky, and F. Vollmer, “Whispering gallery mode carousel – a photonic mechanism for enhanced nanoparticle detection in biosensing,” Opt. Express 17, 6230–6238 (2009).
- L. Landau, E. Lifshitz, and L. Pitaevskiĭ, Electrodynamics of continuous media, Course of theoretical physics (Butterworth-Heinemann, 1984).
- V. Wong and M. A. Ratner, “Explicit computation of gradient and nongradient contributions to optical forces in the discrete-dipole approximation,” J. Opt. Soc. Am. B 23, 1801–1814 (2006).
- M. I. Mishchenko, L. D. Travis, and A. A. Lacis, Scattering, Absorption and Emission of Light by Small Particles (Cambridge University Press, Cambridge, 2002).
- L. Deych and J. Rubin, “Rayleigh scattering of whispering gallery modes of microspheres due to a single dipole scatterer,” Phys. Rev. A 80, 061805 (2009).
- J. T. Rubin and L. Deych, “Ab initio theory of defect scattering in spherical whispering-gallery-mode resonators,” Phys. Rev. A 81, 053827 (2010).
- J. Hu, S. Lin, L. C. Kimerling, and K. Crozier, “Optical trapping of dielectric nanoparticles in resonant cavities,” Phys. Rev. A 82, 053819 (2010).
- M. Abramowitz and I. Stegun, Handbook of Mathematical Functions with Formulas, Graphs and Mathematical Tables (National Bureau of Standards, Washington, DC, 1972).
- J. Chen, J. Ng, S. Liu, and Z. Lin, “Analytical calculation of axial optical force on a rayleigh particle illuminated by gaussian beams beyond the paraxial approximation,” Phys. Rev. E 80, 026607 (2010).