Research ArticlePHYSICS

Rotational Doppler cooling and heating

See allHide authors and affiliations

Science Advances  06 Jan 2021:
Vol. 7, no. 2, eabd6705
DOI: 10.1126/sciadv.abd6705


Doppler cooling is a widely used technique to laser cool atoms, molecules, and nanoparticles by exploiting the Doppler shift associated with translational motion. The rotational Doppler effect arising from rotational coordinate transformation should similarly enable optical manipulation of the rotational motion of nanosystems. Here, we show that rotational Doppler cooling and heating (RDC and RDH) effects embody rich and unexplored physics, including an unexpected strong dependence on particle morphology. For geometrically constrained particles, cooling and heating are observed at red- or blue-detuned laser frequencies relative to particle resonances. In contrast, for nanosystems that can be modeled as solid particles, RDH appears close to resonant illumination, while detuned frequencies produce cooling of rotation. We further predict that RDH can lead to optomechanical spontaneous chiral symmetry breaking, where an achiral particle under linearly polarized illumination starts spontaneously rotating. Our results open up new exciting possibilities to control the rotational motion of nanosystems.


The advent of lasers gave birth to various techniques for trapping and manipulating atoms, small particles, and macroscopic objects using optical forces (15). An important class of optical manipulation methods, such as Doppler cooling, relies on control over the dynamics of a target by exploiting the dependence of optical pressure on its velocity. As a result of the Doppler shift in the laser frequency associated with translational motion, changes in the target velocity lead to increases or reductions in the optical pressure, depending on the frequency detuning of the laser relative to the intrinsic resonances of the target. In consequence, in a particle subject to illumination by two coherent and monochromatic counter-propagating light beams, red- and blue-detuned laser frequencies can decelerate and accelerate the particle motion, respectively, thus cooling and heating the translational temperature of an ensemble of particles. The theoretical description of Doppler cooling and heating is simply a textbook result derived from classical theory. In experiments, Doppler cooling was first achieved for atoms (6, 7) and recently for molecules (812), and it was also generalized to cool the mechanical motion of a macroscopic cavity (13, 14). Compared with cooling, laser heating is less explored, although it can give rise to exciting effects such as mechanical instabilities of illuminated optical cavities (14).

Optical forces can also affect the rotational motion of nanosystems with rotational degrees of freedom, such as molecules (1519) and nanoparticles (2025). More precisely, circularly polarized light with intrinsic angular momentum can exert a torque on the particles, which is the analog of optical pressure resulting from linear momentum transfer. Compared with the optical manipulation of translational motion reviewed above, it is striking to see that the very fundamental question regarding how rotational motion affects the torque exerted by light remains unanswered. Similar to the translational counterpart, such an important piece of knowledge can be used to generalize the physics of Doppler cooling and heating to the rotational degrees of freedom. We note that, in contrast to the relatively simple instances of translationally moving inertial frames, and although rotational cooling has been achieved using specific molecules (2628) and molecular ions (2931), a quantum-mechanical theoretical explanation of rotational cooling involves a complex geometry-dependent coupling between optical excitations, radiation, internal rotations, and vibrational states in noninertial rotating frames (3234); understandably, the general physics underlying rotational cooling is not yet well understood. It is thus crucial to develop a compact classical model to unveil the general physics underlying the rotational cooling effect without resorting to the internal rotational and vibrational quantum states, similar to the theoretical description of its translational counterpart. For a molecular ensemble (1519) that is not tightly confined close to the rotational quantum ground state, the statistical average of the rotation dynamics can be approximately characterized by a continuous rotation velocity (18), using the quantum-classical correspondence. Such a classical model is rigorous and illuminating for studying novel optical manipulations of nano- or microparticles in the rapidly evolving field of nanophotonics (2125, 3538).

Here, we generalize the theory of Doppler cooling and heating, usually discussed for translational motion as schematically illustrated in Fig. 1 (A and B), to deal with rotational degrees of freedom (see Fig. 1C), namely, rotational Doppler cooling and heating (RDC and RDH). We present a theory that is valid for both small particles and molecules rotating with classical velocities, which allows us to predict a marked dependence of rotational laser cooling and heating on the geometry of those particles or molecules. Specifically, for nanosystems that can be modeled as round solid particles, RDC and RDH appear in different frequency regimes compared with its translational Doppler counterpart. More notably, because of RDH, an achiral particle under linearly polarized illumination starts rotating spontaneously, manifesting a spontaneous chiral symmetry breaking (SCSB) effect. We also explore the dynamics of particles under linearly polarized illumination and find a series of stable and metastable configurations in which the particles maintain steady rotation. Our results are of fundamental importance for understanding the optical properties of rotating nanosystems, paving the way for the development of new techniques to optically manipulate rotations of nanosystems.

Fig. 1 Comparison between translational and rotational Doppler cooling and heating.

(A and B) Illustration of translational Doppler cooling. In the laboratory frame [(A), top], a particle is moving with velocity v in the presence of two counter-propagating light waves (wavy curves) of frequency ω. In the frame moving with the particle [(A), bottom], the extinction cross section (σext) of these two waves [see (B)] is affected by Doppler shifts to ω(1 ± v/c), leading to an optical force that pushes the particle to the left, thus causing particle deceleration in the laboratory frame. (C) We consider a particle rotating with angular velocity Ω and illuminated by linearly polarized light of frequency ω (top). The incident light can be decomposed into right circularly polarized (RCP) and left circularly polarized (LCP) components, which are Doppler-shifted to ω ∓ Ω in the frame rotating with the particle (bottom). Black circular arrows denote the directions of the electric field rotation on the plane perpendicular to the wave vector.


Comparing translational and rotational scenarios

In Fig. 1A, we show a particle moving with velocity v and illuminated by two counter-propagating light waves of equal intensity and frequency ω (Fig. 1A, top). In the frame moving with the particle (Fig. 1A, bottom), the two light waves propagating parallel or antiparallel with v are red- or blue-Doppler–shifted, respectively, to ω(1 ± v/c). Assuming a red-detuned laser frequency ω relative to a dominant intrinsic particle resonance ω0 (Fig. 1B), compared with the red-shifted light wave, the blue-shifted antiparallel wave is closer to ω0 and thus undergoes stronger scattering, resulting in deceleration and cooling of the particle. Similarly, heating can be achieved through blue-detuned illumination.

To generalize this scheme to rotational motion, we consider a particle trapped by a linearly polarized laser beam (Fig. 1C, top). Right and left circularly polarized (RCP and LCP) components of the incident light have then equal intensity, thus resembling the two counter-propagating waves in Fig. 1A. Then, they experience different rotational Doppler shifts to ω = ω ∓ Ω (hereafter, upper and lower signs denote quantities for RCP and LCP, respectively), when transforming the system to the frame rotating with the particle (Fig. 1C, bottom), namely, a rotational Doppler effect (15, 39). Similar to translational dynamics, the direction of the torque acting on the particle is then determined by the relative magnitude of the torque exerted by RCP and LCP components. One might naïvely conclude that RDC and RDH should also appear at red- and blue-detuned ω. However, as we show below, this conclusion is only valid for particles having certain geometries. Rotational transformations are noninertial, so that the optical response of the rotating particle is more complicated and thus displays richer physics. The translational trapping of particles at the antinodes of standing waves due to gradient forces cannot be extrapolated to the rotational scenario, which, as we show below, can lead to instabilities of the particle at rest (Ω = 0) for laser frequencies in the heating regime.

Circular polarizabilities of rotating nanosystems

A detailed analysis of the torque exerted by linearly polarized light on rotating small particles or molecules with different geometries requires the knowledge of their optical polarizabilities. Quantum theoretical treatments of this problem commonly assume a rigid body to discuss rotation dynamics (33). However, the coupling between optical excitations and vibrations in molecules cannot be neglected, and importantly, this coupling has a strong dependence on geometry (34), which further complicates a general quantum description. Besides, a naïve application of first-order quantum perturbation theory fails to comply with the optical theorem (40) and therefore neglects radiative corrections that are relevant to the optical torque. To avoid this complication, here we adopt a classical model that effectively captures all relevant physical processes and retains a tutorial character to explain the dependence of RDC and RDH on particle geometry. Our model describes the dipolar optical mode of the particle as a harmonically oscillating effective charge (mass m, charge Q) driven by the electric field of light. Such classical description is capable of adequately addressing transition dipole moments, such as those of molecules, as well as resonant modes in particles. Reassuringly, this model satisfies the optical theorem (see Supplementary text) and, therefore, it correctly accounts for scattering processes. In terms of its radial position vector r, the classical equation of motion of the effective charge becomesr¨=ω02rr̂γ(r·Ωrφ̂)+τr+QmE+Freactm(1)where ω0 is an intrinsic oscillator resonance frequency, we incorporate an internal dissipation force Fdis=γm(r·Ωrφ̂) proportional to a phenomenological damping rate γ times the center-of-mass velocity, the Abraham-Lorentz force mτr with τ = 2Q2/3mc3 introduces corrections due to radiation reaction, and Freact is the force imposed by the boundary defined by the particle geometry. The radial component of Fdis describes the coupling to vibrational phonons, while Freact and the azimuthal components of Fdis account for the coupling between electron motion and particle rotations.

In a thin nanorod (Fig. 2A), the oscillating dipole p produced by the bounded charges is oriented along the rod axis, which rotates with the particle. Likewise, in a nanocross (Fig. 2B), two dipoles (one per branch) are induced, also rotating with the particle. In contrast, for isotropic nanoparticles containing freely moving electrons, such as the nanodisk shown in Fig. 2C, we need to consider the coupling between the two orthogonal charge oscillators through the Coriolis force. At this point, it is useful to consider the electric field E±=(x̂±iŷ)E±eiωt/2 of RCP (+) and LCP (−) light and relate the induced elastic dipole moment p±ω=α±E± (oscillating at frequency ω in the laboratory frame) to the circular polarizability of the particle α±. In this work, the particle rotation axis is assumed to be along the z direction and only the electric field component in the x-y plane is relevant to the physics, so we define the chirality of the light field (i.e., LCP and RCP) according to the rotation direction of Epm (see Fig. 1C), regardless of the direction of the photon wave vector. We note that p±ω is formed by two orthogonal degenerate components. Applying Eq. 1 to these three representative types of particles (see Methods, Supplementary text, and fig. S1), we findα±(ω)=Q2/mA(2)withArod=2[ω02Ω2ω2iω(γ+τω2+3τΩ2)]Across=ω02Ω2ω2i(γω+τω3)Adisk=ω02Ω2ω2i(γω+τω3)where Q denotes the total effective charge in the nanorod, while the cross and disk have a charge of Q oscillating along each of the two orthogonal directions. Such theory can be rigorous in classical scenarios, such as nanoparticles and molecules consisting of larger numbers of atoms, with rotation energy not tightly confined close to the rotational ground state.

Fig. 2 Optical torque acting on rotating particles.

(A to C) Optical response of rotating particles of different geometries, namely, nanorod (A), nanocross (B), and nanodisk (C), to RCP (red) and LCP (blue) incident light. We plot the imaginary part of the polarizability, which is related to the extinction cross section according to σext = (4πω/c) Im {α}. A solid particle such as the nanodisk (C) exhibits different circular dichroism relative to particles with confined electron motion, such as the nanorod (A) and nanocross (B). (D to F) Time-averaged torque Mdr acting on the rotating particles considered in (A) to (C) under linearly polarized illumination. Insets: Mdr versus particle rotating frequency Ω for two typical light frequencies in the cooling and heating regimes. In all cases, damping rates are assumed to be γ = 0.2ω0 and τ−1 = 0.02ω0. Apart from the insets, particles are rotating with angular velocity Ω = 0.1ω0. All frequencies are normalized to the particle resonance frequency ω0, the polarizability is normalized to α0=Q2/mω02, and the torque is normalized to M0 = α0E±2/2. We assume a large resonance linewidth for illustration.

This classical model treats the rotation of particles as a continuum of states characterized by a rotation frequency Ω, which we remark can always be considered as a rigorous approach for nanoparticles and molecules consisting of a large number of atoms. For an ensemble of small molecules containing only a few atoms (e.g., diatomic molecules) and rotating with high angular velocity, according to the quantum-classical correspondence guaranteed by the Ehrenfest theorem, the statistical average of the rotational dynamics can then also be characterized by the rotation frequency Ω, so the molecular dynamics is adequately described through our classical model. For example, the rotational Raman scattering observed in molecules (15, 39) with a high rotation frequency is captured by our model without resorting to quantized rotational modes (see results below). When small molecules are gradually cooled down to extremely low temperatures and the rotation energy is lowered near the quantum rotational ground state, the quantization of the rotational states becomes increasingly more important. However, since the general physical description should be consistent in both classical and quantum scenarios, as already demonstrated in the translational Doppler cooling, our model and the results obtained from it below still offer qualitative understanding to estimate the behavior of electron-phonon/rotations interaction in the quantum limit. As reviewed in Introduction, the complicated morphology-dependent coupling between electrons and phonon/rotations in molecules conceals the essential physics underlying rotational cooling. As we show below, our simple model is capable of intuitively capturing the geometry dependence of the coupling between material electrons and rotational modes.

Optical torque exerted by linearly polarized light

The polarizabilities given by Eq. 2 can fully describe the optical response of small rotating particles and the optical torque produced under external illumination. According to the optical theorem, the extinction cross sections of these particles are given by σext±=4πkIm{α±}, where k = ω/c. However, elastic scattering processes, as described by the cross section σω±=8πk4α±2/3, maintain the angular momentum of light and thus do not lead to a torque on the particle. In contrast, each absorbed circularly polarized photon directly transfers angular momentum ħ to the particle, so we are interested in the absorption cross section σabs±=σext±σω±σω2Ω±. Here, σω2Ω± is the cross section of inelastic scattering at output frequency ω ∓ 2Ω [the so-called rotational Doppler scattering or rotational Raman scattering (15, 39); see Supplementary text], which is only present in rotating anisotropic particles and comes accompanied by an exchange of angular momentum 2ħ and energy 2ħΩ for each scattered photon. For the nanorod, we find σω2Ω±=8π(ω2Ω)4α±2/3c4, while this cross section vanishes for the cross and the disk. Last, the total torque exerted on the particle by RCP and LCP components reduces toM±=±(2ησω2Ω±+σabs±)I/ω(3)where I = cE±2/8π is the incident light intensity, and η = ω/(ω ∓ 2Ω) accounts for the energy shifts of inelastically scattered photons.

Figure 2 (A to C) shows the imaginary part of the circular polarizabilities Im{α±} for the three types of rotating particles described by Eq. 2, whose internal losses dominate over radiative losses (γ ≫ τ−1). We observe a strong circular dichroism (CD) in the rotating thin nanorod and nanocross, characterized by a splitting of 2Ω in the resonance peaks. Similar dichroism can also be observed in particles without internal dissipation (i.e., for γ = 0, see fig. S2). In contrast, the polarizability of the rotating disk with freely moving electrons inside it does not exhibit any resonance splitting. However, the factor Adisk in Eq. 2 predicts a difference in the decay rate for different circular polarizations, which leads to a weak CD, manifested by the discrepancy in the linewidth and magnitude of the resonance peaks in Fig. 2C.

Provided with the CD response of these rotating particles, we can readily conclude that linearly polarized illumination should exert optical torques on the particles because this type of light can be decomposed into RCP and LCP components with equal amplitudes, which contribute with opposite and imbalanced torques. Figure 2 (D to F) shows rigorous results for the total optical torque experienced by the particles described in Fig. 2 (A to C) under linearly polarized light, calculated from Eq. 3 according to Mdr = M+ + M. For particles with γ ≫ τ−1, such as those considered in Fig. 2, Eq. 3 can be approximated as M± ≈ 2 Im {α±}∣E±2, so that the torques in Fig. 2 (D to F) are proportional to Im{α±}, as shown in Fig. 2 (A to C). For example, for rotating nanorods and nanocrosses under red-detuned laser illumination (ω < ω0; Fig. 2, A and B), absorption of the LCP component should be strong compared with that of the RCP component (Im{α} > Im {α+}), so that the total torque exerted by linearly polarized light should decelerate the particle rotation, leading to an RDC effect (blue shaded area in Fig. 2, D and E; see Supplementary text for discussions on the RDC limit). Similarly, rotational acceleration (red shaded area in Fig. 2, D and E) is observed for blue-detuned illumination (ω > ω0). The RDC and RDH effects are also confirmed by the dependence of Mdr on Ω (insets in Fig. 2, D and E), which is similar to the velocity dependence of optical pressure in optical molasses. The conditions for RDC and RDH in rotating nanorods and nanodisks are also similar to those for their translational counterparts. This similarity originates in the fact that the polarizability observed in the rotating frame is equal to that of the motionless particle, α++) = α) ≈ α(ω, Ω = 0), so that the discussions for translational cooling based on Fig. 1B are equivalently applied to rotating particles with these types of geometries. However, for a rotating solid particle, such as the disk shown in Fig. 2F, because of its different CD shown in Fig. 2C, an RDH effect is observed for a nearly resonant laser frequency ω, while RDC is found at off-resonance laser frequencies in both blue- and red-detuned regimes. We also note that the dissipationless limit (γ = 0), in which the torque acting on isotropic particles disappears according to Eq. 3, is generally unphysical because dissipation channels are generally present in molecules and nanoparticles that can quench optical excitations. Even in the limit γ = 0 (rigid molecules), RDC and RDH effects can be observed in anisotropic particles because of the torque arising in rotational Raman scattering (see Eq. 3 and fig. S2).

Rotational stability and dynamics

Besides the optical torque discussed in Fig. 2, light absorption can also increase the temperature of a particle or the vibration energy of a molecule, and the subsequent thermal emission can produce a frictional torque on the particle if it is rotating (4143). This thermal frictional torque arises because the two circularly polarized dipoles p± have different thermal populations and, therefore, there is an imbalance in the angular momentum released through thermal emission, leading to a frictional torque (see Methods)Mfr=cħΣν=±ν0dωρ0(ω)σabsνNν(ω)(4)where N±(ω) = n1) − n0(ω) is the thermal imbalance of particle modes with vacuum, and ρ0(ω) = ω2/3π2c3 is the projected local density of optical states in free space. For rotating particles, when the light frequency falls into the cooling regime, both thermal friction and optical cooling lead to slowing down of rotation; consequently, provided the light frequency is in the heating regime, the driving torque Mdr exerted by the external illumination needs to exceed the thermal friction Mfr to produce acceleration.

An interesting phenomenon can be intuitively foreseen for a particle at rest under linearly polarized illumination when the light frequency is in the heating regime: If the condition ∂ΩMdr > ∂ΩMfr is satisfied at Ω = 0, a small particle rotation induced by any fluctuation can be amplified by RDH. This effect implies an instability of the particle at rest, and considering the chiral symmetry of the system Hamiltonian, such instability manifests as an SCSB process. We note that the SCSB term relates to the dynamics of a single particle, while for an ensemble, the statistical average cancels the macroscopic chirality. To further analyze the particle stability at Ω = 0, we need to find the steady state of the particle temperature under light irradiation, which is reached when the power absorbed from the laser, Pabs=Σ±I±σabs±, is exactly compensated by the thermal-emission power. At arbitrary Ω, the latter is given by (42) (see Methods)Pems=cħΣν=±0ωdωρ0(ω)σabsνNν(ω)(5)

The stability of the particle can be investigated by taking the following aspects into account: (i) given the environment temperature T0, under stationary conditions, the laser intensity is uniquely determined by the particle temperature I(T1); (ii) the driving torque ∂ΩMdr is then uniquely related to T1 through I(T1); (iii) ∂ΩMfr depends on T0 and T1, and its magnitude compared with ∂ΩMdr determines the stability. Following these considerations, at a given incident frequency ω, the rotational stability of the particle can be mapped into a T0 versus T1 plot.

In Fig. 3A, we choose the nanocross of Fig. 2B at rest (Ω = 0) as an example to illustrate the stability in the (T0, T1) plane, and in fact, this diagram of stability also represents the nanorod in Fig. 2A (see discussion on the nanodisk in Fig. 2C and fig. S3). This plot constitutes a universal phase diagram, considering that phase transitions are featured by spontaneous symmetry breaking—in the normal phase, the motionless state is stable; when increasing the intensity of illumination I to heat the particle above the critical temperature dictated by the phase boundary in Fig. 3A (black curves), the particle starts rotating spontaneously and the system enters into an SCSB phase.

Fig. 3 Rotational dynamics of particles under linearly polarized illumination.

(A) Stability of the nanocross considered in Fig. 2 (B and E) at rest under linearly polarized illumination of frequency ω = 1.1ω0 (solid curve) and 1.3ω0 (dashed curve), as a function of vacuum and particle temperatures T0 and T1 normalized to Θ0 = ħω0/kB. For each value of T0, a steady particle temperature T1 is reached at a laser intensity I(T1). The black solid and dashed curves denote phase boundaries for frequencies ω = 1.1ω0 and 1.3ω0, respectively. (B and C) Evolution of the nanocross at arbitrary initial T1 and Ω for light frequency ω = 1.1ω0, with vacuum temperature T0 = 0.4Θ0 and laser intensities I(0.41Θ0) (B) or I(0.5Θ0) (C), which correspond to the black dots I and II in (A). (D) Illustration of the particle state in different phases. Top: The equilibrium state at Ω = 0 is stable [black dot, also in (B)]. Middle: The equilibrium state at Ω = 0 is unstable for high light intensity [gray dot, also in (C)]. Bottom: A metastable configuration (see below). (E) Driving torque acting on the particle rotating at different velocities Ω for light frequencies ω = 1.1ω0 and 1.3ω0. (F) Same as (C) for light frequency ω = 1.3ω0. Two metastable configurations are observed, yielding the equilibrium at Ω = 0 stable [point II in (A) is in the normal phase for ω = 1.3ω0], which is intuitively illustrated in the bottom panel of (D).

To clearly reveal the features of the SCSB phase and find the final stable configuration, we simulate the evolution of Ω and T1 as governed by the dynamical equations of motion Ω·=Mtot/J and T·1=(PabsPemsMtotΩ)/C, where Mtot = MdrMfr, and J and C are the moment of inertia and thermal capacity of the particle, respectively. The evolution of the system for a vacuum temperature T0 = 0.4Θ0, where Θ0 = ħω0/kB, is shown in Fig. 3 (B and C), taking the laser intensities to be either I(0.41Θ0) or I(0.5Θ0), corresponding to points I and II in Fig. 3A, respectively. Here, we discuss the general behavior of rotational heating based on normalized quantities, so we can easily obtain absolute values for specific particles based on Fig. 3 (an example is shown in fig. S4).

For an incident frequency ω = 1.1ω0, the two points in Fig. 3A fall into normal and SCSB phases, respectively (see boundary shown as a black solid curve in Fig. 3A). When the system is in the normal phase, as shown in Fig. 3B, regardless of the initial conditions for Ω and T1, the system evolves toward a trivial stable equilibrium point (black dot) at Ω = 0 and T1 = 0.41Θ0, since the laser intensity used is I(0.41Θ0). By increasing the laser intensity to I(0.5Θ0), the equilibrium point located at Ω = 0 is shifted to higher temperature at T1 = 0.5Θ0 (gray dot). As expected, this equilibrium point becomes unstable, and the particle then starts rotating toward a random direction, which eventually reaches one of the two new stable equilibrium states (black dots, Fig. 3C). The characteristics of the normal and SCSB phases revealed by these dynamical evolutions are intuitively illustrated in Fig. 3D (top and middle), where the black and gray dots correspond to the stable and unstable equilibrium states in Fig. 3 (B and C).

In Fig. 3A, the phase boundary for ω = 1.3ω0 is also indicated through a dashed curve, which lies above the phase boundary for ω = 1.1ω0. As shown in Fig. 3E, given a fixed laser intensity, the particle acquires a large optical torque from a laser at ω = 1.1ω0 at a small rotation frequency Ω ≈ 0 compared to ω = 1.3ω0, so a laser with frequency ω = 1.1ω0 can more easily break the stability at Ω = 0. Point II for ω = 1.3ω0 lies within the normal phase, and the corresponding dynamical evolution is shown in Fig. 3F. Although the equilibrium state at Ω = 0 in Fig. 3F is stable, we also observe two metastable configurations (red dots) at relatively high rotation frequency, which arise because a larger driving torque can be exerted at a far-detuned frequency ω = 1.3ω0 for higher Ω (Fig. 3E). In the metastable state, the particle can maintain its rotation, and only a large perturbation comparable to the energy barrier surrounding the metastable region can break the stability and induce the particle to return to equilibrium at Ω = 0, as illustrated in Fig. 3D (bottom). The stable and metastable states for finite rotation velocity revealed in Fig. 3 (C and F) should also be observed in a nanocross with small anisotropy, since the time-averaged torque acting on the rotating particle should then vanish. However, for Ω ≈ 0, such torque tends to align the long axis of the particle with the direction of linear polarization, which can increase the threshold temperature T1 shown by the curve in Fig. 3A and, consequently, also the corresponding laser intensity.


In Figs. 2 and 3, we choose a large linewidth to illustrate the mechanisms involved in RDC and RDH. For small particles with γ > τ−1, according to ∂ΩMdr ∝ ∂Ω Im {α(ω)} (near Ω = 0), the predicted phenomena are more easily attainable in particles with sharp resonances. We also note that a high rotation frequency Ω comparable with ω0 is adopted here also for intuitive illustrations, whereas RDC and RDH are readily observed for lower rotation frequency, as can be seen in Fig. 2 (insets in D to F). Detailed calculations also show that the predicted RDC and RDH in this work can be achieved for plasmonic nanoparticles, where all necessary conditions are feasible using currently available experimental techniques (see fig. S5).

Regarding various new materials emerging in nanophotonics, we argue that graphene particles with high electron mobility are good candidates to test the predicted phenomena, since they sustain ultranarrow plasmon resonances at long wavelength, which also help reduce thermal friction. Optomechanical SCSB should be achievable in graphene, since the electron temperature of this material (i.e., T1) can easily reach ∼103 K and exceed Θ0 under external illumination. Sparse particles or molecules trapped in a high vacuum could be used to avoid gas friction and coupling between rotational and translational degrees of freedom originating in these scattering events. In experiments, concepts and techniques that are well developed in current optical cooling setups, such as Zeeman splitting and chirping of the light frequency, could be combined with the mechanism here revealed to explore actual applications. We also note that the optical response of a rotating nanoring is similar to a nanocross (34). Considering a ring particle with high electron mobility, with its lattice fixed in the laboratory frame and exposed to linearly polarized illumination, we expect a spontaneous electron current to arise, mimicking the mechanical rotation in our model. However, we argue that such spontaneous electron current cannot be achieved in a nanodisk according to our model dealing only with the linear response (44), because CD as in Fig. 2C cannot be observed in such case. Our work unveils fundamental mechanisms enabling new approaches toward optical trapping and also offers unexplored insights into the optical response of out-of-equilibrium rotating systems.


Modeling circular polarizabilities

Our model extends a previously formulated approach (see details in Supplementary text) (45). To find the circular polarizability of a rotating nanosystem, we need to solve the motion of the effective charge inside the particle (position vector r±), which is governed by Eq. 1 for a circularly polarized optical electric field E±=(x̂±iŷ)E±eiωt/2. This directly leads to an induced electric dipole moment p± = Qr±. In the laboratory frame, the electric dipole moment in a rotating particle of arbitrary symmetry should include in general two Fourier components, p±=p±ω+p±ω2Ω, where p±ω and p±ω2Ω are associated with frequencies ω and ω ∓ 2Ω, respectively. Here, we define the circular polarizability through the elastic electric dipole relation p±ω=α±E±.

In a thin nanorod (Fig. 2A), the oscillating dipole p produced by the bounded charge is oriented along the rod axis, which rotates with the particle. We assume that the rod axis is fixed to x in the Cartesian coordinate system rotating with the nanorod (xy, see fig. S1A), so that Eq. 1 reduces to one-dimensional motion in the rotating frame. For RCP (+) and LCP (−) light, the electric field E±(ω) in the rotating frame is Doppler-shifted to frequency ω ∓ Ω (see Fig. 1C), which forces the charge to oscillate also at the same frequency ω ∓ Ω. However, when observed in the laboratory frame, the charge oscillation in the nanorod includes Fourier components at frequencies ω and ω ∓ 2Ω, which correspond to elastic and inelastic scattering processes, respectively. For an optically isotropic nanoparticle, the total dipole moment p± is the sum of two orthogonal degenerate rotating dipole moments, p±,1 and p±,2. An example of an optically isotropic particle is formed by connecting two orthogonal nanorods at their centers with two branches fixed along x and y (see Fig. 2B and fig. S1B), where we assume a charge Q oscillating along each of the two orthogonal directions. Another example of an optically isotropic particle is a nanodisk (see Fig. 2C and fig. S1C), which contains freely moving electrons. In these optically isotropic particles, the superposition of the elastic dipole moments from the two orthogonal oscillations p±,1ω+p±,2ω add up constructively, while the inelastic dipole moments p±,1ω2Ω+p±,2ω2Ω cancel out. In contrast, inelastic scattering is present in optically anisotropic nanorods, which manifests as a discrepancy in the losses captured by the imaginary parts of Arod and Across, disk in Eq. 2. The real parts of Arod and Across are the same because the electron motions in these two types of particles are both confined to the radial directions, thus leading to similar resonance frequency shifts (see Fig. 2, A and B). However, a disk shows different resonance behavior, as shown by the real part of Adisk and Fig. 2C, because its internal electrons can move freely in the absence of a reaction force Freact in Eq. 1.

Optical cross sections

The time-averaged powers associated with the different physical processes, described by the force terms in Eq. 1, can be easily calculated as P=F·r·. By inserting here the force terms of Eq. 1, including the light field force QE, the Abraham-Lorentz force mτr, the dissipation force Fdis, and the reaction force Freact (solved according to the force balance in Eq. 1), we can find the extinction (Pext), scattering (Psca), dissipation (Pdis), and mechanical power (Preact) contributions to P, respectively. The extinction and scattering powers are consistent with conventional definitions for motionless particles, which describe the powers of the total input energy and the re-emission of the dipole moments. The dissipation and mechanical power are associated with the energy converted to heat and rotational energy, respectively, and these two components constitute the absorption energy that is stored in the particle, Pabs = Psca + Pdis. Likewise, these four force terms also allow us to find the corresponding time-averaged optical torques according to M = ⟨r × F⟩ (see details in Supplementary text). In particular, the sum of the torques exerted by Fdis and Freact stand for the net torque acting on the particle.

As mentioned in previous sections, the time-averaged powers and optical torques permit us to define the optical cross section as σ = P/I. From the partial powers found above, we further obtain the extinction (σext), scattering (σsca), dissipation (σdis), and mechanical (σreact) cross sections. As the absorption power Pabs is the sum of Psca and Pdis, the absorption cross section corresponds to the sum σabs = σdis + σreact. Energy conservation is thus stated in terms of these cross sections as σext = σsca + σabs. In addition, by separating the scattering power into elastic and inelastic components, corresponding to the emission powers associated with pω and pω ∓ p2Ω, we can define the elastic and inelastic scattering cross sections σω and σω ∓ 2Ω. Now, we can recast the net optical torque acting on the particle in a compact form using these cross sections, as shown in Eq. 3.

Torque due to exchanges with the thermal field

To determine the optical torque acting on rotating particles due to thermal radiation, we consider optical absorption from the vacuum thermal field, as well as radiation emission (i.e., the reciprocal process). Vacuum photon modes can be represented by a complete set of plane waves, each of them characterized by a wave vector k and a polarization vector êi, where êi can be chosen to run over two orthogonal linear polarizations. For the thin nanoparticles considered in Fig. 2, only the electric field components in the x-y plane perpendicular to the rotation axis (see fig. S1) are relevant to light absorption and emission. By further projecting the polarization vector onto the in-plane circular polarization unit vectors û±=(x̂±iŷ)/2, we find the effective number of photons absorbed from each incident plane wave mode to be cσabs±êi·û±2/V, where V is the quantization volume; each of these photons transfers an angular momentum ±ħ to the particle. Now, the number of photons in each plane wave mode of frequency ω is given by n0(ω), the Bose-Einstein distribution evaluated at the vacuum temperature T0. By summing contributions from modes with all possible wave vectors k and taking into account the photon population, we find the total absorption rate of thermal photons to beΓ±abs=cVΣi,kσabs± n0(ω)[n1(ω)+1] êi·û±2=13π2c20ω2dω σabs± n0(ω)[n1(ω)+1]where we have carried out the sum over modes using the relations V1Σi,kêi·û±2=(2π)3dkΣiêi·û±2=ω2dω3π2c3, and we have introduced a factor n1) + 1 to represent the fact that a (bosonic) particle polarization mode of Doppler-shifted frequency ω has increased its population from n1) (the Bose-Einstein distribution at the particle temperature T1) to n1) + 1. Similarly, the thermal photon emission rate only differs from the absorption rate in the changes incurred on the thermal distributions, with a particle mode and a photon mode decreasing and increasing their populations by 1, respectively. We findΓ±emi=13π2c20ω2dω σabs± [n0(ω)+1]n1(ω).

The absorption and emission of each photon are accompanied by a transfer of angular momentum ±ħ, so the net torque on the rotating particle isMfri=Σ±(±ħ)(Γ±absΓ±emi)which readily leads to Eq. 4. Likewise, the absorption and emission of each photon are accompanied by an energy exchange ħω, so we find that the net emission power from a rotating particle reduces toPemi=Σ±ħω(Γ±emiΓ±abs)which directly leads to Eq. 5. Incidentally, there is no contribution to the net power that does not experience a rotational Doppler shift because the particle is assumed to be thin along the rotation direction and, therefore, the emission associated with polarization along that axis can be neglected.


Supplementary material for this article is available at

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: Funding: This work has been supported, in part, by the European Research Council (ERC Advanced Grant 789104-eNANO), the Spanish MINECO (MAT2017-88492-R and no. SEV2015-0522), the Catalan CERCA Program, Fundació Privada Cellex, and Fundació Mir-Puig. Author contributions: D.P. conceived the research and carried out the calculations. D.P. and F.J.G.d.A. developed the theory and wrote the paper. H.X. discussed the results. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article