Elliptical orbits of microspheres in an evanescent field

Lulu Liu; Simon Kheifets; Vincent Ginis; Andrea Di Donato; Federico Capasso

doi:10.1073/pnas.1714953114

New Research In

Contributed by Federico Capasso, August 30, 2017 (sent for review January 6, 2017; reviewed by Vincent Daria and Kishan Dholakia)

Significance

Using a highly sensitive particle-tracking scheme, we report an observation of elliptical motion of microparticles driven by a single evanescent field. We show that this behavior is highly tunable and predictable using a theoretical model that accounts for Mie scattering and hydrodynamic drag. The significance is twofold. First, our work represents an important step in understanding the detailed dynamics of microparticles in a fluidic environment—especially near a surface—a prerequisite for effective application of optically driven microparticles as functional elements in optofluidic devices. Second, our method could complement current structured-light approaches for inducing orbital motion in microparticles, with readily tunable parameters of motion including orbital frequency, radius, and ellipticity.

Abstract

We examine the motion of periodically driven and optically tweezed microspheres in fluid and find a rich variety of dynamic regimes. We demonstrate, in experiment and in theory, that mean particle motion in 2D is rarely parallel to the direction of the applied force and can even exhibit elliptical orbits with nonzero orbital angular momentum. The behavior is unique in that it depends neither on the nature of the microparticles nor that of the excitation; rather, angular momentum is introduced by the particle’s interaction with the anisotropic fluid and optical trap environment. Overall, we find this motion to be highly tunable and predictable.

Recently, much work has gone into the investigation of optical forces on micro- and nanoparticles near surfaces, primarily in the context of an electric field localized by a microstructured surface (1 ⇓⇓⇓–5).

Surface-based geometries have raised theoretical excitement due to, for instance, their ability to significantly enhance optical forces (6 ⇓–8) as well as the emergence of lateral forces due to the extraordinary momentum and spin in evanescent waves (9 ⇓⇓–12). From an applied point of view, such geometries can enable miniaturization and parallelization of efficient optical traps enabling integration into optofluidic devices (13, 14). Moreover, light-controlled microspheres near surfaces have found applications as force transducers (15 ⇓⇓–18) and pumps and switches (19, 20), and in general, their collective manipulation is an advancing field (21 ⇓–23).

However, in the case where a microparticle is within several diameters of a surface, optical and hydrodynamic surface effects cannot be neglected. Effects arising from optical coupling or reflections can complicate trapping and detection schemes (24), and hydrodynamic interactions can cause the motion of a microparticle to become highly nontrivial (25). Despite this, little quantitative study has been done on the dynamics of optically driven particles near a surface. In studies introducing new schemes to optically manipulate matter, the particle’s response function is often ignored (11, 12, 26 ⇓⇓–29).

In this work, we investigate the dynamics of a system with both near field optical forces and surface-induced hydrodynamic effects (see Fig. 1). By driving the particle with an oscillating force from a modulated evanescent field and tracking the particle’s motion closely in two dimensions, we map out a range of dynamics that can arise from the interplay of optical and hydrodynamic surface effects. We find that the magnitude and direction of the optical force depends on particle size. We also observe that the trajectory of the particle in general does not follow the direction of the force. Instead, the shape and the orientation of the trajectory vary with modulation frequency and distance of the particle from the surface, a result of the anisotropy in both the hydrodynamic drag and the optical trap spring constant.

Fig. 1.

The optically trapped microsphere in its anisotropic environment. The optical trap (660 nm, 12 mW) has a Rayleigh length, b, larger than the beam waist, w, resulting in a restoring force approximately 5 times stronger in the lateral direction. Additionally, near-wall hydrodynamic coupling results in anisotropic drag that is consistently larger in the vertical direction. As a result, the particle has a different response function in the two directions. The optical force of the pump beam has components in z and x. At a driving frequency between the two cutoff frequencies, the particle’s response can have large phase differences between the two directions, generating elliptical orbits. Fk represents the restoring force of the trap. The dissipative drag and the stochastic force are not visualized. For a detailed description of the setup, see Fig. S1.

Fig. S1.

Setup. An optical tweezer is formed by the 660-nm laser beam focused by a high-N.A. water immersion objective (OBJ). The trap confines polystyrene microspheres in water near an AR-coated glass surface; the back-scattered trap beam is used to determine its lateral position via a balanced detector (BAL). The bead-surface separation is controlled by vertical translation of the objective. A low-intensity evanescent wave used for vertical position detection is produced by total internal reflection of a 637-nm (probe) beam. Probe light scattered by the particle is sent to a low-noise photodiode (PD1). A second, high-intensity, 785-nm (pump) evanescent wave exerts an optical force on the particle in the x–z plane. A chopper modulates the intensity of the pump beam, which is monitored by PD2. This signal as well as the signals from BAL and PD1 are digitized and recorded by a computer for lock-in signal processing.

In particular, we find that under certain conditions, controlled, elliptical motion of a microsphere can be achieved in this configuration, without the use of a beam that carries either orbital or spin angular momentum. This motion is of theoretical and practical interest as it arises due to the particle’s interaction with the inherently anisotropic mechanical environment nearby a solid interface as it moves in response to a periodic driving force. It has several distinguishing features that may make it an applicable technique for generating orbital angular momentum, alongside existing methods predominantly using structured beams carrying orbital or spin angular momentum (30 ⇓⇓⇓⇓⇓–36). For a detailed discussion, please see The Evanescent Wave Technique as a Method of Angular Momentum Generation, following Experiment.

Forces on an Optically Driven and Trapped Microsphere

An optically trapped microsphere in fluid is well-modeled as a stochastically driven, damped harmonic oscillator in three dimensions (37, 38):m⋆xi¨(t)=−κijxj(t)−γijxj˙(t)+Bi(t)+Fi(t).[1]The terms on the right of Eq. 1 correspond, in order, to the restoring force of the optical trap, the dissipative drag in fluid, the stochastic force due to thermal fluctuations, and an externally applied force. Here, we apply a force and consider motion only in the two measured (x and z) dimensions.

Due to the symmetries in our system, the trap stiffness, κij, is a diagonal tensor. For a single beam optical trap, its diagonal elements are not equal, κx≠κz; that is, the harmonic potential is anisotropic. Usually, the magnitude of the spring constant in the lateral direction is several times larger than the spring constant in the axial direction (39). In our model describing the particle motion, we approximate that κi is a constant in z, an assumption that is generally valid until the separation becomes smaller than the Debye length in water (see Height-Dependent Trap Spring Constant).

In the absence of boundary effects, viscous drag, γ, is isotropic, given by Stokes’ Law, γ0=6πηR, where η is the viscosity of the fluid and R is the radius of a sphere. The presence of the boundary results in an anisotropic, height-dependent drag coefficient (40, 41). The predicted drag in the direction perpendicular to the surface (z) is always larger than that in the lateral direction (x) (see Height-Dependent Damping Coefficient). Due to the spherical symmetry of the particle, in low Reynolds number hydrodynamics (Re≪1), we may assume that the drag forces parallel and perpendicular to a plane boundary are uncoupled, such that γij is diagonal with unequal elements (42).

B→(t) is the fluctuating thermal force due to collisions between the microparticle and molecules of the fluid. In each dimension, the fluctuations are expected to be uncorrelated, with ⟨Bi(t)⟩=0, where B∼(ω) is spectrally flat. Its mean-squared magnitude is predicted by the fluctuation-dissipation theorem to be ⟨Bi(t)2⟩=2γikBT, where kB is the Boltzmann constant. It is important to note that thermal noise increases not only with temperature but also with increased viscous drag and is the dominant source of noise in our system (18).

Lastly, the external force, which drives the periodic motion of the particle, is an optical force produced by the interaction of the microsphere with the evanescent field of a totally internally reflected TM-polarized plane wave. In the dipolar approximation, the interaction can be broadly understood to have two components: one proportional to the gradient of the field, which attracts the particle toward the surface, and the other proportional to the linear momentum of the wave, which pushes the particle along the surface (see Optical Force on a Dipolar Particle in an Evanescent Field).

For small, lossless particles (Im[α]≪Re[α]), the optical force is predominantly oriented along the z direction. But as the particle increases in size, its interaction with the evanescent wave becomes more complex (43). In this regime, called the Mie regime, light is scattered into every direction in the xz plane. As a result, the direction of the net applied force is observed to rotate toward the horizontal direction. A detailed analysis of the forces acting on the particle in the Mie regime is shown in Optical Force on a Mie Particle in an Evanescent Field.

Height-Dependent Trap Spring Constant

In predicting the trajectories in the main text, the spring constant was assumed to be a constant independent of height. This assumption is roughly valid in the lateral direction but begins to break down in the axial direction near the surface.

At very close separations less than 350 nm from the surface, the spring constant in the axial direction of the beam increases sharply (see Fig. S5). This is due to the particle’s electrostatic repulsion from the surface arising as an additional confining force (discussed in ref. 18). The trap potential narrows and the particle becomes more tightly confined.

Fig. S5.

Fitted trap stiffnesses, κx and κz, in the horizontal and vertical directions as a function of height for a 2-μm-radius PS particle. The line is to guide the eye.

This effect is neglected in the trajectory figures in the main text. For this reason, the trajectory predictions very close to the surface are expected to disagree somewhat from the measured results.

Height-Dependent Damping Coefficient

For slow motion of a microsphere in fluid in close proximity to a hard surface with no-slip boundary conditions, Brenner and Goldman derived expressions for the height-dependent viscous drag in the direction perpendicular to the wall, γz(z) (40), and in the direction parallel to the wall, γx(z) (41), as correction to the Stokes law result γ0=6πηR.

The drag is generally larger than the case with no boundary present due to hydrodynamic coupling between the sphere and the wall and generally larger in the perpendicular direction than the parallel direction (see Fig. S4). On contact with the wall, the correction coefficient diverges for the perpendicular direction—that is, γz→∞—while the correction coefficient for the parallel direction reaches a constant value γx→C. Far from the wall, the boundary’s influence diminishes and γx(z) and γz(z) approach each other and the isotropic result of Stokes law. The diffusion coefficient depends inversely on the drag as D(z)=kBT/γ(z).

Fig. S4.

The measured perpendicular damping coefficient as a function of height above surface γ(z) for a PS microsphere is fit to theoretical hydrodynamic predictions (40). The fitted radius is then used to predict the damping coefficient in the lateral direction, which is generally smaller. Fitted results yield a radius of 1.1 μm for the sphere at a temperature of 22 °C.

Although the height-dependent damping complicates force measurement, since it results in height dependence of the mechanical susceptibility and the thermal noise, we take advantage of the relationship between damping in the two directions to perform the calibration between signal intensity and position in the x direction and to better estimate the z=0 point or the point-of-contact in the z direction (18).

Optical Force on a Dipolar Particle in an Evanescent Field

In the dipole approximation, the force on a small, polarizable particle can be calculated analytically. The time-averaged Lorentz force acting upon a dipolar particle with radius a and permittivity ϵp, surrounded by a background with permittivity ϵb, consists of three terms:F→=Re[α]4∇|E0→|2+Im[α]2k→|E0→|2−Im[α]2Im[E0→∇E0→],[S4]where E0 is the amplitude of the electric field with wave vector k and α is the polarizability of the particle in the fluid. For small particles (ka≪1), this polarizability can be approximated by α=α0(1+2/3ik3|α0|2) with α0=a3(ϵr−1)/(ϵr+2) and ϵr=ϵp/ϵb. The first term in Eq. 4 corresponds to the gradient force, whereas the second and the third term are scattering forces. In our setup, we consider an evanescent field generated by an internally reflected TM-polarized plane wave. In this case, the first term pulls the particle toward the surface, the second one pushes the particle along the surface in the direction of propagation, and the final term vanishes.

The ratio of these forces perpendicular to and longitudinal along the surface is thus given by:|FzFx|=1Re[α]∇|E0→|22Im[α]kx|E0→|2,[S5]where E0→=E01→Eexp[ikxx−qz]. |Fz/Fx| is thus entirely determined by the ratio of the real and the imaginary parts of the particle’s polarizability α and the evanescent field’s wave vector:|FzFx|=qkxRe[α]Im[α].[S6]For a given angle of incidence, q and kz remain fixed. The polarizability α of a small particle is related to its static polarizability α0 using Draine’s radiation correction:α=α0(1−23ik3α0)−1≈α0(1+23ik3α0∗).[S7]For lossless particles, the ratio of the real and the imaginary parts of the polarizability thus equals:Re[α]Im[α]=(23k3α0)−1.[S8]The ratio of the force in the z and the x direction grows to infinity for decreasing particle radius a, as the static polarizability α0 is proportional to a3. This effect is shown in the small particle regime (a<0.1μm) of Fig. S7B.

Fig. S7.

(A) In the dipole regime (a<0.1μm), the amplitude of both the lateral and the perpendicular force monotonically increase with increasing particle size. In the Mie regime, the forces continue to grow, with signature Mie oscillations superimposed. (B) The relative amplitude of the attractive force versus the longitudinal force grows to infinity in the dipole regime. In the Mie regime, the forces in both directions are comparable in amplitude.

Optical Force on a Mie Particle in an Evanescent Field

To calculate the optical force on dielectric particles with sizes on the order of the wavelength of light, one needs to solve the Mie scattering problem. As shown in ref. 43, the integral of the Maxwell stress tensor can be expressed as a function of the Mie scattering coefficients. The optical force that acts on a dielectric particle in an evanescent field can thus be evaluated using an algebraic combination of the scattering coefficients. Using this algorithm, the time-averaged attractive force and longitudinal forces was calculated as:⟨Fz⟩a2ϵ0E02=Re[iα24∑l=1∞∑m=−ll(l+m+2)(l+m+1)(2l+1)(2l+3)l(l+2)×(2nw2al,mal+1,m+1∗+nw2al,mAl+1,m+1∗+nw2Al,mal+1,m+1∗+2bl,mbl+1,m+1∗+bl,mBl+1,m+1∗+Bl,mbl+1,m+1∗)+(l−m+1)(l−m+2)(2l+1)(2l+3)l(l+2)×(2nw2al+1,m−1al,m∗+nw2al+1,m−1Al,m∗+nw2Al+1,m−1al,m∗+2bl+1,m−1bl,m∗+bl+1,m−1Bl,m∗+Bl+1,m−1bl,m∗)−(l+m+1)(l−m)nw×(−2al,mbl,m+1∗+2bl,mal,m+1−al,mBl,+1∗+bl,mAl,m+1∗+Bl,mal,m+1∗−Al,mbl,m+1∗)],[S9]⟨Fx⟩a2ϵ0E02=Im[−α22∑l=1∞∑m=−ll(l−m+1)(l+m+1)(2l+1)(2l+3)l(l+2)×(2nw2al+1,mal,m∗+nw2al+1,mAl,m∗+nw2Al+1,mal,m∗+2bl+1,mbl,m∗+bl+1,mBl,m∗+Bl+1,mbl,m∗+nwm(2al,mbl,m∗+al,mBl,m∗+Al,mbl,m∗))].[S10]The Mie scattering coefficients Al,m, Bl,m, al,m, and bl,m in the previous equation are related to the incident (superscript (i)) and the scattered (superscript (s)) fields and are defined by:Er(i)=a2E0r2∑l=1∞∑m=−lll(l+1)Al,mψl(nwk0r)Yl,m(θ,ϕ),[S11]Hr(i)=a2H0r2∑l=1∞∑m=−lll(l+1)Bl,mψl(nwk0r)Yl,m(θ,ϕ),[S12]Er(s)=a2E0r2∑l=1∞∑m=−lll(l+1)al,mξl(1)(nwk0r)Yl,m(θ,ϕ),[S13]Hr(s)=a2H0r2∑l=1∞∑m=−lll(l+1)bl,mξl(nwk0r)Yl,m(θ,ϕ),[S14]where a is the radius of the spherical particle, Yl,m(θ,ϕ) are the spherical harmonics, ψl(r)=rjl(r) with jl(r) the spherical Bessel function, and ξl(r)=rhl(1)(r) with hl(1)(r) the spherical Hankel function of the first kind.

The value of these scattering coefficients was determined using a traditional Mie scattering algorithm, and the infinite series was truncated at a value l for which the relative magnitude of the optical force contribution dropped below 10−7. The result is shown in Fig. S7A. It is interesting to note that in the Mie regime, the longitudinal force along the surface approaches the same amplitude as the attractive force toward the surface. This is clear in the large particle regime (a>0.5μm) in Fig. S7 B and Inset.

Microparticle Mechanical Response

The mechanical response of our microsphere to an applied force is frequency-dependent and can be found by solving Eq. 1. We can neglect the particle’s inertia while working at time scales much larger than the momentum relaxation time, τp=m⋆/γ, which, in our system, is below 1 μs (44, 45). Inserting a time-harmonic solution, we then find thatx∼i(ω)=χij(ω)F∼j(ω),[2]in which we identified the mechanical susceptibility χij(ω) of the system as (χ−1(ω))ij=κij+iωγij.

The motion is separable, with two different cutoff frequencies, fc=κ/2πγ, in the x and z directions. Eq. 2 implies that though Fx(t) and Fz(t) may be in phase, as is the case with our optically driven system, x(t) and z(t) need not. The maximum phase lag between motion in the two directions is π/2, making it possible to induce elliptical trajectories of the particle. Furthermore, the linearity of the system ensures that Eq. 2 applies whether the force is random or periodic. Thus, in our case, where the external force on our particle is composed of both a driven and a stochastic component, the thermal motion of the microsphere may be considered separately from the time evolution of the particle’s mean position.

Experiment

A schematic of the experiment is shown in Fig. 1 (also see Fig. S1 for details). All experiments are performed at room temperature in a water-filled, closed, 25-μm-deep microfluidic chamber. An antireflection (AR) coating is applied to the bottom glass–water interface of the chamber to eliminate standing-wave modulation of the trap laser beam (46). Each polystyrene microsphere is optically trapped by a focused CW laser beam (660-nm-wavelength) via a high numerical aperture (N.A. = 1.2) water-immersion objective. The vertical position of the focus, and consequently the height of the trapped bead, is adjusted using a piezo—by lifting and lowering the objective relative to the chamber. In addition to the trap beam, which enters the chamber from above, two lasers, both p-polarized, are incident from below at greater than critical angle (θc = 61.4^∘): a low-power (<1 mW), 637-nm beam acts as the probe for the vertical (z) position of the particle (15, 47), and a second laser (785 nm, around 100 mW average power before modulation) generates the periodic optical force that pulls on the microparticle. The second beam, which we refer to as the pump beam, is modulated by a chopper set to a desired frequency. A detailed explanation of the analysis may be found in Analysis Flow.

Fig. 2 plots the predicted and observed 2D trajectories of a 2-μm PS microsphere under the influence of a periodic optical force. Predictions were made based on the optical force, mean trap stiffnesses, and bead radius fitted to experimental data. The measured trajectories are mean trajectories, extracted from 100-s measurements where each point is the averaged 2D position of the particle at a certain phase in the chopper actuation cycle. As the frequency of the modulation is increased, several distinct behaviors are observed that can be understood in relation to the cutoff frequencies in the two spatial directions.

Fig. 2.

(A) Predicted and measured 2D motion of a 1-μm-radius bead when driven by a periodic optical force from an evanescent field. The blue arrows indicate the direction of the optical force. Each tile in A corresponds to a particular height (z) and modulation frequency (ω0) and is numbered corresponding to its location on the grid (B). The details of the measurement parameters are given in Table S1. (B) Grid of predicted particle trajectories where the horizontal axis represents the driving frequency and the vertical axis is the height above the surface. The displacement amplitude in nanometers is defined as the length of the semimajor axis of each orbit. Calculated cutoff frequencies in the x and z directions are drawn as blue lines. At all heights, the cutoff frequency in the z direction is lower than that in the x direction, due to a larger trap stiffness in the lateral direction. (C) Measured amplitude and phase of frequency-dependent mechanical susceptibility, χ(f), at a height of 400 nm, compared with predictions.

View this table:

Table S1.

Details of experimental parameters in Fig. 2

Below both cutoff frequencies, the position of the particle is in phase with the applied force. In this case, the particle’s motion in x and z is in phase with one another, and it undergoes linear oscillation with a direction determined primarily by the ratio of the trap stiffnesses κx/κz. In our case, the ratio is about 5. Above both cutoff frequencies, the motion is again linear as the oscillations are now out of phase with the applied force in both directions. However, the direction of motion is now determined by the ratio of the damping coefficients γx/γz. As the particle nears the surface, γz diverges, while γx approaches a constant value, and the motion becomes more and more parallel to the wall. But between the two cutoff frequencies, the motions in the two spatial directions go out of phase with one another, and stable, elliptical orbits are established.

In Fig. 3 we compare our measurements against Mie theory predictions for a fixed modulation frequency. As the microparticle radius increases, the net optical force from the evanescent field increases for a given field strength and configuration, and the direction of the force rotates slowly toward the horizontal. In addition, since hindered diffusion theory predicts drag near a surface to increase as a function of the ratio z/R, where z is the separation and R is the radius, increasing the radius has similar effects as decreasing separation. Therefore, large beads tend to move parallel to the surface, regardless of their direction of excitation.

Fig. 3.

(A) Predicted and measured motion of 0.6- to 2.5-μm-radius beads when driven by a periodic optical force from an evanescent field that is switched on and off at 20 Hz. The blue arrows indicate the direction of the optical force. Each tile in A corresponds to a particular height and bead radius and is numbered corresponding to its location on the grid (B). The details of the measurement parameters are given in Table S2. (B) Grid of predicted particle trajectories based on Mie theory calculations, where the horizontal axis represents bead radius and the vertical axis is the height above the surface. The displacement amplitude in nanometers is defined as the length of the semimajor axis of each orbit. The trajectories labeled 7 and 8 are for a 2.85-μm-radius bead, outside the scope of our Mie theory calculations. (C) Ratio of measured optical forces in the z and x directions compared with Mie theory. Shaded area shows the uncertainty in the microsphere index of refraction (n = 1.575 ± 0.005). Generally, as microparticle radius increases, the direction of the optical force rotates toward the horizontal.

View this table:

Table S2.

Details of experimental parameters in Fig. 3

Analysis Flow

The unknown parameters in our experiment (calibrated positions, lock-in amplitudes, damping coefficients, and trap stiffnesses in the two spatial directions) are determined sequentially in the procedure diagrammed above. Calibrated positions in z follow a logarithmic relation to scattered light intensity and are determined according to the procedure detailed in ref. 46. In Fig. S2, 0 the parameters are defined, with the known values in black and the unknown values in red. As the analysis progresses, parameters determined in each step turn black in color; in this way, the logical progression of our process is revealed.

Fig. S2.

Analysis flow for determination of unknown parameters: 0 describes what is known at the start of this analysis; 1 and 5 refer to the lock-in algorithm detailed in Lock-In Algorithm; 2 and 4 refer to the best fit to mean-squared-displacement (MSD) of particle position described in Fitting the MSD; and 3 references the fit to hindered diffusion model described in Height-Dependent Damping Coefficient.

The Evanescent Wave Technique as a Method of Angular Momentum Generation

The observed phenomenon may find application as an efficient generator of small-radius orbital motion. While we do not suggest that this method supplants the study and use of structured light beams in colloidal systems, here we provide a detailed quantitative comparison and discuss the distinguishing features of our system.

One manner of quantifying angular motion is to consider the torque applied by the particle to the fluid, τapp(t)=r→(t)×γv→(t), where γ is the tensor described in Eq. 1. For a particle undergoing steady-state, elliptical motion in fluid, this quantity is on average the same as the torque applied to the particle. Fig. 4 shows τapp accumulating over 100 s for a bead driven at two different frequencies, which agree well with predictions. Another method of quantifying this angular motion is to calculate the constant, nonzero intrinsic orbital angular momentum in the steady state, given by ⟨Lin⟩=ω|x∼||z∼|sin(ϕz−ϕx), where ϕx and ϕz are the phases of the complex mechanical susceptibilities χx and χz in the x and z directions, respectively, and ⟨Lin⟩ is normalized by mass (see Angular Momentum Analysis). This quantity is useful in allowing direct comparison with existing reports in literature, where we do not have access to values of the damping parameter, γ, but can estimate velocity and orbital radius. In our unoptimized set-up, the largest Lin was measured to be 4.9×10−15m2/s, where ϕz−ϕx was closest to π/2. The configuration for this measurement was a 2-μm-diameter bead driven at 4.6 Hz at a separation of around 200 nm from the surface (corresponding to point 5 in Fig. 2). In comparing this value with previous reports for dielectric particles driven by structured light, such as by optical vortices (30 ⇓⇓–33), we find our value comparable when normalized by laser power. With optimization of driving frequency and bead radius or changes in the power or wavelength of the pump beam, further gains may be expected. For details regarding this calculation, see Angular Momentum Analysis and Comparing the Generation of Angular Momentum with Common Existing Techniques.

Fig. 4.

Measured and predicted torque applied to fluid integrated over time for a 2-μm PS bead driven at 11 Hz and 5 Hz frequency at a height of about 200 nm, corresponding to points 5 and 6 in Fig. 2. The force on the bead is the same in both curves.

It is relevant, at this point, to bring up another body of previous work on the orbital motion of colloids, distinct from structured light studies. One specific scheme involves a particle positioned in between a pair of misaligned beams (48). It is a part of a general class of studies in which particles are observed to undergo circulation while immersed in a static force field with nonzero curl (25, 49). In addition, Angelsky et al. observed orbital angular momentum transfer in circularly polarized Gaussian beams as a manifestation of the macroscopic “spin energy flow,” thus confirming the theory of inhomogeneously polarized paraxial beams (34, 50).

We are careful to distinguish the nature of the orbital motion described in our current report from this class of studies. The motion observed in our work is not fully determined by a particle’s passive interaction with a static optical force field but rather by a periodically driven microsphere’s mechanical response. Since our system exhibits driven rather than passive motion, an array of such trapped particles can be made to move in phase, enabling unique collective behaviors that are thus far unexplored.

In addition, while a particle or a collection of particles may be easily made to orbit in a plane parallel to a nearby surface using existing techniques, methods for obtaining orbits in a plane perpendicular to such a surface are less readily available. Particularly in on-chip geometries where space may be an issue, our method of actuation may provide a useful alternative.

Finally, since the motion observed in our work is actively controlled, main parameters of orbital motion—such as rotational rate and radius—are decoupled and can be independently tuned. For instance, to increase the orbital radius in our experiment without affecting the rotation rate, one can simply raise the pump beam’s intensity without changing its modulation frequency. It should be noted, however, that in our scheme the maximum radius of rotation cannot exceed the limits of the optical trap. If complete flexibility in orbital radii is desired, structured light approaches are more suitable. Additionally, if the goal is near-circular motion, a nearby surface or some other symmetry-breaking element is required, making the work described here potentially unsuitable for manipulation in bulk fluid or in vacuum. Also, larger radii of rotation can be achieved in general by the use of structured light, as the requirement that a particle remains in a single-beam Gaussian harmonic trap may prove too limiting for certain applications.

Thus, our technique can be considered as complementary to existing methods for generating orbital angular momentum, rather than substitutive. The ideas underlying both can be combined to extend the parameter range of orbital angular momentum generation in microparticles using optical forces.

Angular Momentum Analysis

In this section, we discuss how we extract the angular momentum from the experimental data and how we compare it to our analytical model. The intrinsic angular momentum—the angular momentum per unit kilogram—was calculated directly by evaluating L=r→×v→, both for the experimental data and the analytical model.

However, one needs to be cautious when extracting the mean angular momentum from the experimental data. The differentiation of the displacement to calculate the velocity causes a blue-shift in the noise spectrum, amplifying high-frequency noise. The data, as a result, need to be low pass-filtered to be interpreted correctly. This is achieved using a triangular windowing of the velocity data in the time-domain. The resulting calculated angular momentum as a function of the windowing width is shown in Fig. S8.

Fig. S8.

The angular momentum as a function of the windowing width of the velocity data. At small windowing widths, the angular momentum rapidly drops as the high-frequency noise is filtered out. The angular momentum of the experimental data (blue circles) is extracted by fitting the curve in the region of large windowing widths to the analytical curve of a windowed noiseless ellipse (red curve) and by evaluating that fit at zero windowing width. The result is in excellent agreement with the value predicted from our analytical model (orange curve). (A) Frequency, 4.6 Hz; height, 170 nm. (B) Frequency, 4.6 Hz; height, 230 nm. (C) Frequency, 11 Hz; height, 160 nm. (D) Frequency, 11 Hz; height, 360 nm. The aforementioned results were obtained with the particle having a radius equal to 1.1μm.

In these graphs, we can clearly distinguish different regimes. For low windowing widths (high frequency), noise dominates the signal. As the windowing width increases, this noise is filtered out, but undersampling of the ellipse begins to affect the angular momentum resolved, causing a systematic underestimate of the true angular momentum.

To resolve this issue, we fit the curve in the low-frequency regime—where the high-frequency noise is filtered out—to an equivalently undersampled analytical orbit that does not contain any noise. We then extract the angular momentum of the experimental data as the angular momentum of the fitted ellipse at zero windowing width. As shown in Fig. S8, the predicted values of the angular momenta (red curves) are in excellent agreement with the experimentally fitted ones (yellow curves).

In the main text we provide an analytical estimate for the intrinsic angular momentum expected from such an analysis, which depends only on the angular frequency, ω; amplitudes of oscillation, |x∼| and |z∼|; and the response phase difference in the two directions, ϕz−ϕx. A brief derivation is provided below:⟨L/m⟩=⟨r→×v→⟩[S15]where the position and velocity vectors are given by:r→=|x∼|cos(ωt−ϕx)x^+|z∼|cos(ωt−ϕz)z^,andv→=dr→dt=−|x∼|ωsin(ωt−ϕx)x^−|z∼|ωsin(ωt−ϕz)z^.[S16]Taking their cross-product, simplifying, and taking the time average over each cycle then gives:⟨L/m⟩=ω|x∼||z∼|sin(ϕz−ϕx)[S17]We note that the calculations performed in this paper assume that the particle motion, together with the gradient of spring and damping constants, at each fixed height are sufficiently small that the parameters γ, κ, as well as the amplitude and direction of applied optical force may be treated as constants and the equations of motion considered linear.

Comparing the Generation of Angular Momentum with Common Existing Techniques

In our setup, we measure an intrinsic angular momentum, Lin=r→×v→, of 4.89×10−15m2/s. We cite two results in the literature for comparison (31, 33). The first citation reports using a 700-mW beam to generate rotation with angular frequency 0.07 rad/s; the beads were seen to rotate within an ellipse around 5 μm in radius. Scaling down to the powers used in our experiment (80 mW), the intrinsic angular momentum equals 2.28×10−14m2/s, which is about 4.7 times our result.

The second citation reports particles rotating at 0.06 Hz around an inner ring of a Bessel beam of 2.9-μm radius at a power of 600 mW at the sample surface. Rescaling to our power levels, this intrinsic angular momentum equals 5.64×10−14m2/s, which is about 11 times our result.

Thus, our unoptimized system produced angular momentum in optically driven microparticles within an order of magnitude of previous structured-light results. If desired, the angular momentum generation in our system can be greatly enhanced by performing some simple optimizations. For instance, two orders of magnitude improvement can be expected, if, simply, the pump beam is focused so that the spot size is shrunk by a factor of 3. Currently, the low-intensity beam used for excitation has a semimajor axis of around 100 μm. Focusing to a spot size of around 30 μm would increase the field intensity by a factor of 10, and the angular momentum generated by a factor of 100.

In our geometry it is possible to independently tune the radius and angular velocity through the manipulation of the chopping frequency and the power of the pump beam.

The angular frequency at which the particle moves is directly given by the chopping frequency of the pump (evanescent) beam. However, to obtain near circular motion, one needs to operate this frequency between the two corner frequencies that determine the particle’s dynamics (blue lines in Fig. 2B). The corner frequencies are given by fcornerx,z=κx,z/γx,z, where κx,z and γx,z are the spring constants and the damping coefficients in the x and the z direction, respectively. In the vertical (z) direction, the corner frequency can be made smaller than 0.5 Hz. This defines a certain lower limit of the chopping, or orbital, frequency. In the longitudinal (x) direction, the maximal corner frequency is tunable by adjusting the trap stiffness (most directly by increasing the trap laser power). Although this parameter cannot be increased indefinitely, we can easily tune it such that the upper corner frequency lies above 50 Hz (as shown in Fig. S9). The range of angular rotational rates can thus be adjusted approximately between 0.5 Hz and 50 Hz.

Fig. S9.

PSD in the x direction, where the corner frequency is up to ∼50 Hz. This higher frequency was obtained by increasing the trap stiffness. Raw data (light red), running average (dark red), and predicted PSW (black line) are based on a fit of the MSD.

The radius of the orbital motion, on the other hand, can be controlled by adjusting the power (intensity) of the pumping beam. Indeed, the radius can be estimated as r=F/(iω∗γ+κ), where ω is the angular chopping frequency, γ is the damping constant, κ is the spring constant, and F has an amplitude proportional to the intensity of the pump beam. This shows that the orbital radius will roughly scale with the pump power.

Conclusion

Using a highly sensitive detection scheme, we investigated the dynamics of trapped microspheres under the influence of a temporally modulated force. We find good agreement between our model and the measured motion validating our method of generating elliptical orbits with sustained nonzero orbital angular momentum.

Experimental Set-Up

For details about other portions of the experimental set-up, including the microfluidic chambers, calibration, AR-coating methods, and imaging optics, please see the supplemental materials of Liu et al. (18).

Lock-In Algorithm

The lock-in algorithm is the same as the one used in ref. 18 and is described in detail in the supplemental materials of that paper and is broadly summarized below.

In general, a lock-in measurement is used to estimate the amplitude and phase of a sinusoidal signal obscured by noise using a reference signal. In this case, the output of the monitor photodiode that records the intensity of the chopper-modulated pump beam serves as the reference signal.

The algorithm generates a sin wave with unity amplitude that is in-phase with the square wave with a fundamental frequency equal to the frequency of modulation.

The noisy (discrete) signal zi is multiplied by the complex reference signal, exp[iϕi], resulting in the complex lock-in signal li=2ziexp[iϕi].

To correct for the phase delay caused by differences in the signal pathways between the pump monitor photodiode and the TIRM detection photodiode, the delay was directly measured and compensated in the data analysis (18).

The output of the lock-in algorithm corresponds to the amplitude of the signal at the fundamental frequency. The coefficient of the first term of the Fourier series expansion of a square wave is larger than the amplitude of the square wave, by a factor of π/4. The forces we report in our results corresponds to the amplitude of this fundamental component. It is larger than the amplitude of the square wave but smaller than the peak-to-peak height of the square wave, which corresponds to the DC force of an unmodulated beam, by a factor of 2/π.

Fitting the MSD

The position power spectral density (PSD), Sx(w), pictured in Fig. S3, is defined as the Fourier transform of the autocorrelation function ⟨x(t)x(0)⟩:Sx(w)=∫−∞∞dteiωt⟨x(t)x(0)⟩[S1]The position autocorrelation function is then related to the MSD by:MSD=⟨(x(t)−x(0))2⟩=2⟨x2⟩−2⟨x(t)x(0)⟩[S2]where ⟨x2⟩ is the variance of the function x(t).

Fig. S3.

PSD in the z and x directions. Pictured are raw data, running average, and the predicted PSD based on fit to MSD.

For a Brownian particle in an optical trap driven by an oscillating force at frequency ωL, ignoring inertial effects, the MSD is given by:MSD(t)=2kBTκ[1−e−t/τc]+(2A)2[1−cos(ωLt)],[S3]where τc=γ/κ.

The damping coefficient γ and trap stiffness κ are extracted by a least-squares fit to the MSD of the Brownian particle, with lock-in amplitude A and frequency ωL set to the values determined using the lock-in algorithm described in the previous section. One fit is performed per measurement and per spatial dimension at a given height.

Net Force on Microparticles for Various Trajectories

Fig. S6 shows the computed net force acting on a 2.2-μm-radius PS microsphere under periodic forcing at various stages of its orbit. The figure shows the physical origin of these elliptical trajectories.

Fig. S6.

Net force normalized by effective mass of particle in nanometers per second squared for a few select trajectories of a 2.2-μm-radius particle actuated at frequencies between 2 and 100 Hz. Force direction and relative magnitude are plotted on top of particle position at each point in the particle’s orbit.

Acknowledgments

We thank the groups of Evelyn Hu and David Weitz at Harvard University for shared equipment and laboratory access. All fabrication was done in the Harvard Center for Nanoscale Sciences (CNS) clean room facility. We acknowledge the support of NSF GFRP Grant DGE1144152 and Research Foundation Flanders Grant 12O9115N.

Footnotes

↵¹To whom correspondence may be addressed. Email: capasso{at}seas.harvard.edu or ginis{at}seas.harvard.edu.

Author contributions: L.L., S.K., V.G., and F.C. designed research; L.L., S.K., and A.D.D. performed research; F.C. oversaw research; L.L. and V.G. contributed new reagents/analytic tools; L.L., S.K., V.G., and A.D.D. analyzed data; and L.L., S.K., V.G., A.D.D., and F.C. wrote the paper.
Reviewers: V.D., Australian National University; and K.D., University of St. Andrews.
The authors declare no conflict of interest.
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1714953114/-/DCSupplemental.

Freely available online through the PNAS open access option.

References

↵
1. Wang K,
2. Schonbrun E,
3. Steinvurzel P,
4. Crozier KB
(2011) Trapping and rotating nanoparticles using a plasmonic nano-tweezer with an integrated heat sink. Nat Commun 2:469.
.
OpenUrl CrossRef PubMed
↵
1. Righini M,
2. Volpe G,
3. Girard C,
4. Petrov D,
5. Quidant R
(2008) Surface plasmon optical tweezers: Tunable optical manipulation in the femtonewton range. Phys Rev Lett 100:8–11.
.
OpenUrl
↵
1. Min C, et al.
(2013) Focused plasmonic trapping of metallic particles. Nat Commun 4:2891.
.
OpenUrl
↵
1. Wang K,
2. Schonbrun E,
3. Crozier KB
(2009) Propulsion of gold nanoparticles with surface plasmon polaritons: Evidence of enhanced optical force from near-field coupling between gold particle and gold film. Nano Lett 9:2623–2629.
.
OpenUrl PubMed
↵
1. Grigorenko A,
2. Roberts N,
3. Dickinson M,
4. Zhang Y
(2008) Nanometric optical tweezers based on nanostructured substrates. Nat Photon 2:365–370.
.
OpenUrl
↵
1. Povinelli ML, et al.
(2005) Evanescent-wave bonding between optical waveguides. Opt Lett 30:3042–3044.
.
OpenUrl PubMed
↵
1. Volpe G,
2. Quidant R,
3. Badenes G,
4. Petrov D
(2006) Surface plasmon radiation forces. Phys Rev Lett 96:238101.
.
OpenUrl CrossRef PubMed
↵
1. Woolf D,
2. Kats MA,
3. Capasso F
(2014) Spoof surface plasmon waveguide forces. Opt Lett 39:517–520.
.
OpenUrl
↵
1. Bliokh KY,
2. Bekshaev AY,
3. Nori F
(2013) Extraordinary momentum and spin in evanescent waves. Nat Commun 5:14.
.
OpenUrl
↵
1. Wang SB,
2. Chan CT
(2014) Lateral optical force on chiral particles near a surface. Nat Commun 5:1–8.
.
OpenUrl CrossRef PubMed
↵
1. Rodríguez-Fortuño FJ,
2. Engheta N,
3. Martínez A,
4. Zayats AV
(2015) Lateral forces on circularly polarizable particles near a surface. Nat Commun 6:8799, and erratum (2015) 6:10263.
.
↵
1. Hayat A,
2. Mueller JPB,
3. Capasso F
(2015) Lateral chirality-sorting optical forces. Proc Natl Acad Sci USA 112:13190–13194.
.
OpenUrl Abstract/FREE Full Text
↵
1. Grier DG
(2003) A revolution in optical manipulation. Nature 424:810–816.
.
OpenUrl CrossRef PubMed
↵
1. Neuman KC,
2. Block SM
(2004) Optical trapping. Rev Sci Instrum 75:2787–2809.
.
OpenUrl CrossRef PubMed
↵
1. Prieve DC,
2. Bike SG,
3. Frej NA
(1990) Brownian motion of a single microscopic sphere in a colloidal force field. Faraday Discuss Chem Soc 90:209–222.
.
OpenUrl CrossRef
↵
1. Florin EL,
2. Pralle A,
3. Hörber JH,
4. Stelzer EH
(1997) Photonic force microscope based on optical tweezers and two-photon excitation for biological applications. J Struct Biol 119:202–211.
.
OpenUrl CrossRef PubMed
↵
1. Geraci AA,
2. Papp SB,
3. Kitching J
(2010) Short-range force detection using optically cooled levitated microspheres. Phys Rev Lett 105:101101.
.
OpenUrl CrossRef PubMed
↵
1. Liu L,
2. Kheifets S,
3. Ginis V,
4. Capasso F
(2016) Subfemtonewton force spectroscopy at the thermal limit in liquids. Phys Rev Lett 116:228001.
.
OpenUrl
↵
1. Leach J,
2. Mushfique H,
3. di Leonardo R,
4. Padgett M,
5. Cooper J
(2006) An optically driven pump for microfluidics. Lab A Chip 6:735–739.
.
OpenUrl
↵
1. Terray A,
2. Oakey J,
3. Marr DWM
(2002) Microfluidic control using colloidal devices. Science 296:1841–1844.
.
OpenUrl Abstract/FREE Full Text
↵
1. MacDonald M,
2. Spalding G,
3. Dholakia K
(2003) Microfluidic sorting in an optical lattice. Nature 426:421–424.
.
OpenUrl CrossRef PubMed
↵
1. Lutz C,
2. Kollmann M,
3. Bechinger C
(2004) Single-file diffusion of colloids in one-dimensional channels. Phys Rev Lett 93:026001.
.
OpenUrl CrossRef PubMed
↵
1. Garces-Chavez V, et al.
(2006) Extended organization of colloidal microparticles by surface plasmon polariton excitation. Phys Rev B 73:085417.
.
OpenUrl
↵
1. Jonas A,
2. Zemanek P,
3. Florin EL
(2001) Single-beam trapping in front of reflective surfaces. Opt Lett 26:1466–1468.
.
OpenUrl CrossRef PubMed
↵
1. Roichman Y,
2. Sun B,
3. Stolarski A,
4. Grier DG
(2008) Influence of nonconservative optical forces on the dynamics of optically trapped colloidal spheres: The fountain of probability. Phys Rev Lett 101:1–5.
.
OpenUrl CrossRef
↵
1. Wang S,
2. Chan CT
(2016) Strong optical force acting on a dipolar particle over a multilayer substrate. Opt Express 24:2235–2241.
.
OpenUrl
↵
1. Salary MM,
2. Mosallaei H
(2016) Strong optical force acting on a dipolar particle over a multilayer substrate. Phys Rev B 24:2235–2241.
.
OpenUrl
↵
1. Bekshaev AY,
2. Bliokh KY,
3. Nori F
(2013) Mie scattering and optical forces from evanescent fields: A complex-angle approach. Opt Express 21:7082–7095.
.
OpenUrl CrossRef PubMed
↵
1. Yang Y,
2. Zang WP,
3. Zhao ZY,
4. Tian JG
(2012) Optical forces on Mie particles in an airy evanescent field. Opt Express 20:25681–25692.
.
OpenUrl PubMed
↵
1. O’neil A,
2. MacVicar I,
3. Allen L,
4. Padgett M
(2002) Intrinsic and extrinsic nature of the orbital angular momentum of a light beam. Phys Rev Lett 88:053601.
.
OpenUrl CrossRef PubMed
↵
1. Volke-Sepulveda K,
2. Garces-Chavez V,
3. Chavez-Cerda S,
4. Arlt J,
5. Dholakia K
(2002) Orbital angular momentum of a high-order Bessel light beam. J Opt B Quantum Semiclassical Opt 4:S82.
.
OpenUrl
↵
1. Garces-Chavez V, et al.
(2003) Observation of the transfer of the local angular momentum density of a multiringed light beam to an optically trapped particle. Phys Rev Lett 91:093602.
.
OpenUrl PubMed
↵
1. Loxpez-Mariscal C,
2. Gutierrez-Vega JC,
3. Milne G,
4. Dholakia K
(2006) Orbital angular momentum transfer in helical Mathieu beams. Opt Express 14:4182–4187.
.
OpenUrl PubMed
↵
1. Angelsky OV, et al.
(2012) Circular motion of particles suspended in a Gaussian beam with circular polarization validates the spin part of the internal energy flow. Opt Express 20:11351–11356.
.
OpenUrl PubMed
↵
1. Angelsky OV, et al.
(2012) Orbital rotation without orbital angular momentum: Mechanical action of the spin part of the internal energy flow in light beams. Opt Express 20:3563–3571.
.
OpenUrl PubMed
↵
1. Xiao G,
2. Yang K,
3. Luo H,
4. Chen X,
5. Xiong W
(2016) Orbital rotation of trapped particle in a transversely misaligned dual-fiber optical trap. IEEE Photon J 8:1–8.
.
OpenUrl
↵
1. Uhlenbeck GE,
2. Ornstein LS
(1930) On the theory of the Brownian motion. Phys Rev 36:823–841.
.
OpenUrl CrossRef
↵
1. Berg-Sorensen K,
2. Flyvbjerg H
(2004) Power spectrum analysis for optical tweezers. Rev Sci Instrum 75:594–612.
.
OpenUrl CrossRef
↵
1. Nieminen TA, et al.
(2014) Optical tweezers: Theory and modelling. J Quant Spectrosc Radiat Transf 146:59–80.
.
OpenUrl
↵
1. Brenner H
(1961) The slow motion of a sphere through a viscous fluid towards a plane surface. Chem Eng Sci 16:242–251.
.
OpenUrl CrossRef
↵
1. Goldman A,
2. Cox RG,
3. Brenner H
(1967) Slow viscous motion of a sphere parallel to a plane wall i: Motion through a quiescent fluid. Chem Eng Sci 22:637–651.
.
OpenUrl CrossRef
↵
1. Nägele G
(1996) On the dynamics and structure of charge-stabilized suspensions. Phys Rep 272:215–372.
.
OpenUrl
↵
1. Almaas E,
2. Brevik I
(1995) Radiation forces on a micrometer-sized sphere in an evanescent field. J Opt Soc Am B 12:2429–2438.
.
OpenUrl
↵
1. Huang R, et al.
(2011) Direct observation of the full transition from ballistic to diffusive Brownian motion in a liquid. Nat Phys 7:576–580.
.
OpenUrl CrossRef
↵
1. Kheifets S,
2. Simha A,
3. Melin K,
4. Li T,
5. Raizen MG
(2014) Observation of Brownian motion in liquids at short times: Instantaneous velocity and memory loss. Science 343:1493–1496.
.
OpenUrl Abstract/FREE Full Text
↵
1. Liu L,
2. Woolf A,
3. Rodriguez AW,
4. Capasso F
(2014) Absolute position total internal reflection microscopy with an optical tweezer. Proc Natl Acad Sci USA 111:E5609–E5615.
.
OpenUrl Abstract/FREE Full Text
↵
1. Prieve DC,
2. Walz JY
(1993) Scattering of an evanescent surface wave by a microscopic dielectric sphere. Appl Opt 32:1629–1641.
.
OpenUrl CrossRef PubMed
↵
1. Xiao G,
2. Yang K,
3. Luo H,
4. Chen X,
5. Xiong W
(2016) Orbital rotation of trapped particle in a transversely misaligned dual-fiber optical trap. IEEE Photon J 8:1–8.
.
OpenUrl
↵
1. Pesce G,
2. Volpe G,
3. De Luca AC,
4. Rusciano G,
5. Volpe G
(2009) Quantitative assessment of non-conservative radiation forces in an optical trap. Europhys Lett 86:38002.
.
OpenUrl
↵
1. Bekshaev A,
2. Bliokh KY,
3. Soskin M
(2011) Internal flows and energy circulation in light beams. J Opt 13:053001.
.
OpenUrl

Article Alerts

Email Article

Citation Tools

Request Permissions

Proceedings of the National Academy of Sciences: 114 (42)

Table of Contents

Submit

Article Classifications

Physical Sciences
Engineering

[1] ↵
Wang K,
Schonbrun E,
Steinvurzel P,
Crozier KB
(2011) Trapping and rotating nanoparticles using a plasmonic nano-tweezer with an integrated heat sink. Nat Commun 2:469.
.
OpenUrl CrossRef PubMed

[2] Wang K,

[3] Schonbrun E,

[4] Steinvurzel P,

[5] Crozier KB

[6] ↵
Righini M,
Volpe G,
Girard C,
Petrov D,
Quidant R
(2008) Surface plasmon optical tweezers: Tunable optical manipulation in the femtonewton range. Phys Rev Lett 100:8–11.
.
OpenUrl

[7] Righini M,

[8] Volpe G,

[9] Girard C,

[10] Petrov D,

[11] Quidant R

[12] ↵
Min C, et al.
(2013) Focused plasmonic trapping of metallic particles. Nat Commun 4:2891.
.
OpenUrl

[13] Min C, et al.

[14] ↵
Wang K,
Schonbrun E,
Crozier KB
(2009) Propulsion of gold nanoparticles with surface plasmon polaritons: Evidence of enhanced optical force from near-field coupling between gold particle and gold film. Nano Lett 9:2623–2629.
.
OpenUrl PubMed

[15] Wang K,

[16] Schonbrun E,

[17] Crozier KB

[18] ↵
Grigorenko A,
Roberts N,
Dickinson M,
Zhang Y
(2008) Nanometric optical tweezers based on nanostructured substrates. Nat Photon 2:365–370.
.
OpenUrl

[19] Grigorenko A,

[20] Roberts N,

[21] Dickinson M,

[22] Zhang Y

[23] ↵
Povinelli ML, et al.
(2005) Evanescent-wave bonding between optical waveguides. Opt Lett 30:3042–3044.
.
OpenUrl PubMed

[24] Povinelli ML, et al.

[25] ↵
Volpe G,
Quidant R,
Badenes G,
Petrov D
(2006) Surface plasmon radiation forces. Phys Rev Lett 96:238101.
.
OpenUrl CrossRef PubMed

[26] Volpe G,

[27] Quidant R,

[28] Badenes G,

[29] Petrov D

[30] ↵
Woolf D,
Kats MA,
Capasso F
(2014) Spoof surface plasmon waveguide forces. Opt Lett 39:517–520.
.
OpenUrl

[31] Woolf D,

[32] Kats MA,

[33] Capasso F

[34] ↵
Bliokh KY,
Bekshaev AY,
Nori F
(2013) Extraordinary momentum and spin in evanescent waves. Nat Commun 5:14.
.
OpenUrl

[35] Bliokh KY,

[36] Bekshaev AY,

[37] Nori F

[38] ↵
Wang SB,
Chan CT
(2014) Lateral optical force on chiral particles near a surface. Nat Commun 5:1–8.
.
OpenUrl CrossRef PubMed

[39] Wang SB,

[40] Chan CT

[41] ↵
Rodríguez-Fortuño FJ,
Engheta N,
Martínez A,
Zayats AV
(2015) Lateral forces on circularly polarizable particles near a surface. Nat Commun 6:8799, and erratum (2015) 6:10263.
.

[42] Rodríguez-Fortuño FJ,

[43] Engheta N,

[44] Martínez A,

[45] Zayats AV

[46] ↵
Hayat A,
Mueller JPB,
Capasso F
(2015) Lateral chirality-sorting optical forces. Proc Natl Acad Sci USA 112:13190–13194.
.
OpenUrl Abstract/FREE Full Text

[47] Hayat A,

[48] Mueller JPB,

[49] Capasso F

[50] ↵
Grier DG
(2003) A revolution in optical manipulation. Nature 424:810–816.
.
OpenUrl CrossRef PubMed

[51] Grier DG

[52] ↵
Neuman KC,
Block SM
(2004) Optical trapping. Rev Sci Instrum 75:2787–2809.
.
OpenUrl CrossRef PubMed

[53] Neuman KC,

[54] Block SM

[55] ↵
Prieve DC,
Bike SG,
Frej NA
(1990) Brownian motion of a single microscopic sphere in a colloidal force field. Faraday Discuss Chem Soc 90:209–222.
.
OpenUrl CrossRef

[56] Prieve DC,

[57] Bike SG,

[58] Frej NA

[59] ↵
Florin EL,
Pralle A,
Hörber JH,
Stelzer EH
(1997) Photonic force microscope based on optical tweezers and two-photon excitation for biological applications. J Struct Biol 119:202–211.
.
OpenUrl CrossRef PubMed

[60] Florin EL,

[61] Pralle A,

[62] Hörber JH,

[63] Stelzer EH

[64] ↵
Geraci AA,
Papp SB,
Kitching J
(2010) Short-range force detection using optically cooled levitated microspheres. Phys Rev Lett 105:101101.
.
OpenUrl CrossRef PubMed

[65] Geraci AA,

[66] Papp SB,

[67] Kitching J

[68] ↵
Liu L,
Kheifets S,
Ginis V,
Capasso F
(2016) Subfemtonewton force spectroscopy at the thermal limit in liquids. Phys Rev Lett 116:228001.
.
OpenUrl

[69] Liu L,

[70] Kheifets S,

[71] Ginis V,

[72] Capasso F

[73] ↵
Leach J,
Mushfique H,
di Leonardo R,
Padgett M,
Cooper J
(2006) An optically driven pump for microfluidics. Lab A Chip 6:735–739.
.
OpenUrl

[74] Leach J,

[75] Mushfique H,

[76] di Leonardo R,

[77] Padgett M,

[78] Cooper J

[79] ↵
Terray A,
Oakey J,
Marr DWM
(2002) Microfluidic control using colloidal devices. Science 296:1841–1844.
.
OpenUrl Abstract/FREE Full Text

[80] Terray A,

[81] Oakey J,

[82] Marr DWM

[83] ↵
MacDonald M,
Spalding G,
Dholakia K
(2003) Microfluidic sorting in an optical lattice. Nature 426:421–424.
.
OpenUrl CrossRef PubMed

[84] MacDonald M,

[85] Spalding G,

[86] Dholakia K

[87] ↵
Lutz C,
Kollmann M,
Bechinger C
(2004) Single-file diffusion of colloids in one-dimensional channels. Phys Rev Lett 93:026001.
.
OpenUrl CrossRef PubMed

[88] Lutz C,

[89] Kollmann M,

[90] Bechinger C

[91] ↵
Garces-Chavez V, et al.
(2006) Extended organization of colloidal microparticles by surface plasmon polariton excitation. Phys Rev B 73:085417.
.
OpenUrl

[92] Garces-Chavez V, et al.

[93] ↵
Jonas A,
Zemanek P,
Florin EL
(2001) Single-beam trapping in front of reflective surfaces. Opt Lett 26:1466–1468.
.
OpenUrl CrossRef PubMed

[94] Jonas A,

[95] Zemanek P,

[96] Florin EL

[97] ↵
Roichman Y,
Sun B,
Stolarski A,
Grier DG
(2008) Influence of nonconservative optical forces on the dynamics of optically trapped colloidal spheres: The fountain of probability. Phys Rev Lett 101:1–5.
.
OpenUrl CrossRef

[98] Roichman Y,

[99] Sun B,

[100] Stolarski A,

[101] Grier DG

[102] ↵
Wang S,
Chan CT
(2016) Strong optical force acting on a dipolar particle over a multilayer substrate. Opt Express 24:2235–2241.
.
OpenUrl

[103] Wang S,

[104] Chan CT

[105] ↵
Salary MM,
Mosallaei H
(2016) Strong optical force acting on a dipolar particle over a multilayer substrate. Phys Rev B 24:2235–2241.
.
OpenUrl

[106] Salary MM,

[107] Mosallaei H

[108] ↵
Bekshaev AY,
Bliokh KY,
Nori F
(2013) Mie scattering and optical forces from evanescent fields: A complex-angle approach. Opt Express 21:7082–7095.
.
OpenUrl CrossRef PubMed

[109] Bekshaev AY,

[110] Bliokh KY,

[111] Nori F

[112] ↵
Yang Y,
Zang WP,
Zhao ZY,
Tian JG
(2012) Optical forces on Mie particles in an airy evanescent field. Opt Express 20:25681–25692.
.
OpenUrl PubMed

[113] Yang Y,

[114] Zang WP,

[115] Zhao ZY,

[116] Tian JG

[117] ↵
O’neil A,
MacVicar I,
Allen L,
Padgett M
(2002) Intrinsic and extrinsic nature of the orbital angular momentum of a light beam. Phys Rev Lett 88:053601.
.
OpenUrl CrossRef PubMed

[118] O’neil A,

[119] MacVicar I,

[120] Allen L,

[121] Padgett M

[122] ↵
Volke-Sepulveda K,
Garces-Chavez V,
Chavez-Cerda S,
Arlt J,
Dholakia K
(2002) Orbital angular momentum of a high-order Bessel light beam. J Opt B Quantum Semiclassical Opt 4:S82.
.
OpenUrl

[123] Volke-Sepulveda K,

[124] Garces-Chavez V,

[125] Chavez-Cerda S,

[126] Arlt J,

[127] Dholakia K

[128] ↵
Garces-Chavez V, et al.
(2003) Observation of the transfer of the local angular momentum density of a multiringed light beam to an optically trapped particle. Phys Rev Lett 91:093602.
.
OpenUrl PubMed

[129] Garces-Chavez V, et al.

[130] ↵
Loxpez-Mariscal C,
Gutierrez-Vega JC,
Milne G,
Dholakia K
(2006) Orbital angular momentum transfer in helical Mathieu beams. Opt Express 14:4182–4187.
.
OpenUrl PubMed

[131] Loxpez-Mariscal C,

[132] Gutierrez-Vega JC,

[133] Milne G,

[134] Dholakia K

[135] ↵
Angelsky OV, et al.
(2012) Circular motion of particles suspended in a Gaussian beam with circular polarization validates the spin part of the internal energy flow. Opt Express 20:11351–11356.
.
OpenUrl PubMed

[136] Angelsky OV, et al.

[137] ↵
Angelsky OV, et al.
(2012) Orbital rotation without orbital angular momentum: Mechanical action of the spin part of the internal energy flow in light beams. Opt Express 20:3563–3571.
.
OpenUrl PubMed

[138] Angelsky OV, et al.

[139] ↵
Xiao G,
Yang K,
Luo H,
Chen X,
Xiong W
(2016) Orbital rotation of trapped particle in a transversely misaligned dual-fiber optical trap. IEEE Photon J 8:1–8.
.
OpenUrl

[140] Xiao G,

[141] Yang K,

[142] Luo H,

[143] Chen X,

[144] Xiong W

[145] ↵
Uhlenbeck GE,
Ornstein LS
(1930) On the theory of the Brownian motion. Phys Rev 36:823–841.
.
OpenUrl CrossRef

[146] Uhlenbeck GE,

[147] Ornstein LS

[148] ↵
Berg-Sorensen K,
Flyvbjerg H
(2004) Power spectrum analysis for optical tweezers. Rev Sci Instrum 75:594–612.
.
OpenUrl CrossRef

[149] Berg-Sorensen K,

[150] Flyvbjerg H

[151] ↵
Nieminen TA, et al.
(2014) Optical tweezers: Theory and modelling. J Quant Spectrosc Radiat Transf 146:59–80.
.
OpenUrl

[152] Nieminen TA, et al.

[153] ↵
Brenner H
(1961) The slow motion of a sphere through a viscous fluid towards a plane surface. Chem Eng Sci 16:242–251.
.
OpenUrl CrossRef

[154] Brenner H

[155] ↵
Goldman A,
Cox RG,
Brenner H
(1967) Slow viscous motion of a sphere parallel to a plane wall i: Motion through a quiescent fluid. Chem Eng Sci 22:637–651.
.
OpenUrl CrossRef

[156] Goldman A,

[157] Cox RG,

[158] Brenner H

[159] ↵
Nägele G
(1996) On the dynamics and structure of charge-stabilized suspensions. Phys Rep 272:215–372.
.
OpenUrl

[160] Nägele G

[161] ↵
Almaas E,
Brevik I
(1995) Radiation forces on a micrometer-sized sphere in an evanescent field. J Opt Soc Am B 12:2429–2438.
.
OpenUrl

[162] Almaas E,

[163] Brevik I

[164] ↵
Huang R, et al.
(2011) Direct observation of the full transition from ballistic to diffusive Brownian motion in a liquid. Nat Phys 7:576–580.
.
OpenUrl CrossRef

[165] Huang R, et al.

[166] ↵
Kheifets S,
Simha A,
Melin K,
Li T,
Raizen MG
(2014) Observation of Brownian motion in liquids at short times: Instantaneous velocity and memory loss. Science 343:1493–1496.
.
OpenUrl Abstract/FREE Full Text

[167] Kheifets S,

[168] Simha A,

[169] Melin K,

[170] Li T,

[171] Raizen MG

[172] ↵
Liu L,
Woolf A,
Rodriguez AW,
Capasso F
(2014) Absolute position total internal reflection microscopy with an optical tweezer. Proc Natl Acad Sci USA 111:E5609–E5615.
.
OpenUrl Abstract/FREE Full Text

[173] Liu L,

[174] Woolf A,

[175] Rodriguez AW,

[176] Capasso F

[177] ↵
Prieve DC,
Walz JY
(1993) Scattering of an evanescent surface wave by a microscopic dielectric sphere. Appl Opt 32:1629–1641.
.
OpenUrl CrossRef PubMed

[178] Prieve DC,

[179] Walz JY

[180] ↵
Xiao G,
Yang K,
Luo H,
Chen X,
Xiong W
(2016) Orbital rotation of trapped particle in a transversely misaligned dual-fiber optical trap. IEEE Photon J 8:1–8.
.
OpenUrl

[181] Xiao G,

[182] Yang K,

[183] Luo H,

[184] Chen X,

[185] Xiong W

[186] ↵
Pesce G,
Volpe G,
De Luca AC,
Rusciano G,
Volpe G
(2009) Quantitative assessment of non-conservative radiation forces in an optical trap. Europhys Lett 86:38002.
.
OpenUrl

[187] Pesce G,

[188] Volpe G,

[189] De Luca AC,

[190] Rusciano G,

[191] Volpe G

[192] ↵
Bekshaev A,
Bliokh KY,
Soskin M
(2011) Internal flows and energy circulation in light beams. J Opt 13:053001.
.
OpenUrl

[193] Bekshaev A,

[194] Bliokh KY,

[195] Soskin M

Main menu

User menu

Search

New Research In

Elliptical orbits of microspheres in an evanescent field

Significance

Abstract

Forces on an Optically Driven and Trapped Microsphere

Height-Dependent Trap Spring Constant

Height-Dependent Damping Coefficient

Optical Force on a Dipolar Particle in an Evanescent Field

Optical Force on a Mie Particle in an Evanescent Field

Microparticle Mechanical Response

Experiment

Analysis Flow

The Evanescent Wave Technique as a Method of Angular Momentum Generation

Angular Momentum Analysis

Comparing the Generation of Angular Momentum with Common Existing Techniques

Conclusion

Experimental Set-Up

Lock-In Algorithm

Fitting the MSD

Net Force on Microparticles for Various Trajectories

Acknowledgments

Footnotes

References

Citation Manager Formats

Article Classifications

You May Also be Interested in

Similar Articles

Articles

PNAS Portals

Information

Main menu

User menu

Search

New Research In

Physical Sciences

Featured Portals

Articles by Topic

Social Sciences

Featured Portals

Articles by Topic

Biological Sciences

Featured Portals

Articles by Topic

Elliptical orbits of microspheres in an evanescent field

Significance

Abstract

Forces on an Optically Driven and Trapped Microsphere

Height-Dependent Trap Spring Constant

Height-Dependent Damping Coefficient

Optical Force on a Dipolar Particle in an Evanescent Field

Optical Force on a Mie Particle in an Evanescent Field

Microparticle Mechanical Response

Experiment

Analysis Flow

The Evanescent Wave Technique as a Method of Angular Momentum Generation

Angular Momentum Analysis

Comparing the Generation of Angular Momentum with Common Existing Techniques

Conclusion

Experimental Set-Up

Lock-In Algorithm

Fitting the MSD

Net Force on Microparticles for Various Trajectories

Acknowledgments

Footnotes

References

Citation Manager Formats

Sign up for Article Alerts

Article Classifications

Jump to section

You May Also be Interested in

Similar Articles

Articles

PNAS Portals

Information