Skip to content

Optomechanics

Cavity optomechanics studies the coupling between electromagnetic fields and mechanical motion. The mechanical system may be a suspended mirror, membrane, nanobeam, microwave drumhead, levitated particle, phononic mode, or collective elastic coordinate. The common ingredient is that the optical or microwave resonance frequency depends on a mechanical displacement.

This page is an application map. It does not replace the canonical harmonic-oscillator pages, a full device-physics treatment, or the general derivation of input–output theory. Its purpose is to show how radiation-pressure coupling, continuous position measurement, measurement backaction, sideband cooling, and the standard quantum limit appear in open-system language.

A minimal single-mode optomechanical Hamiltonian is

Hℏ=ωca†a+Ωmb†b−g0a†a(b+b†),\frac{H}{\hbar} = \omega_c a^\dagger a + \Omega_m b^\dagger b - g_0 a^\dagger a \left( b+b^\dagger \right),

where aa annihilates a cavity photon, bb annihilates a mechanical phonon, ωc\omega_c is the cavity frequency, Ωm\Omega_m is the mechanical frequency, and g0g_0 is the single-photon optomechanical coupling rate.

The mechanical displacement operator is

x=xzpf(b+b†),xzpf=ℏ2meffΩm.x = x_{\mathrm{zpf}} \left( b+b^\dagger \right), \qquad x_{\mathrm{zpf}} = \sqrt{\frac{\hbar}{2m_{\mathrm{eff}}\Omega_m}} .

If the cavity frequency depends on displacement, then

g0=−xzpf∂ωc∂x.g_0 = - x_{\mathrm{zpf}} \frac{\partial\omega_c}{\partial x}.

The sign convention is not important by itself; the measurable quantities are frequency shifts, damping rates, spectra, and phases. The essential physics is that photon number exerts a force on the mechanical oscillator, while mechanical displacement shifts the phase and frequency of the outgoing field.

Most experiments use a strong coherent drive. In a frame rotating at the drive frequency ωd\omega_d, write the cavity operator as

a=α+d,a = \alpha+d,

where α\alpha is the coherent intracavity amplitude and dd represents fluctuations. Keeping terms up to first order in fluctuations gives the linearized Hamiltonian

Hlinℏ=−Δd†d+Ωmb†b−(G∗d+Gd†)(b+b†),\frac{H_{\mathrm{lin}}}{\hbar} = - \Delta d^\dagger d + \Omega_m b^\dagger b - \left( G^*d+Gd^\dagger \right) \left( b+b^\dagger \right),

with detuning

Δ=ωd−ωc,G=g0α.\Delta = \omega_d-\omega_c, \qquad G = g_0\alpha .

The drive changes the problem qualitatively. The single-photon coupling g0g_0 can be small, while the enhanced coupling GG is large enough to cool, amplify, measure, or hybridize mechanical motion. The price is that drive noise, heating, and classical technical fluctuations may enter the same ports as the quantum field.

A standard Markovian model includes cavity loss at rate κ\kappa and mechanical damping at rate Γm\Gamma_m:

ρ˙=−iℏ[H,ρ]+κD[a]ρ+Γm(nˉm+1)D[b]ρ+ΓmnˉmD[b†]ρ.\dot\rho = - \frac{i}{\hbar} [H,\rho] + \kappa\mathcal D[a]\rho + \Gamma_m(\bar n_m+1)\mathcal D[b]\rho + \Gamma_m\bar n_m\mathcal D[b^\dagger]\rho .

For optical cavities at ordinary temperatures the thermal photon occupation is usually negligible. For microwave cavities it may not be negligible unless the mode is well thermalized and cold. The mechanical bath occupation is

nˉm=1exp⁡(ℏΩm/kBT)−1.\bar n_m = \frac{1}{ \exp(\hbar\Omega_m/k_{\mathrm B}T)-1 }.

This master equation is not the whole story in every device. Strong colored noise, non-Markovian supports, internal two-level fluctuators, excess laser noise, photothermal forces, gas damping, and nonlinearities may require a more detailed bath model. Still, the Lindblad model is the useful baseline against which those corrections are judged.

In the linear regime the outgoing field carries information about mechanical displacement. Schematically,

dout(t)=din(t)+κext d(t),d_{\mathrm{out}}(t) = d_{\mathrm{in}}(t) + \sqrt{\kappa_{\mathrm{ext}}}\,d(t),

and an appropriate output quadrature contains a contribution proportional to x(t)x(t). Homodyne detection therefore implements a continuous weak measurement of mechanical position or of a filtered mechanical quadrature.

The measurement record may be written abstractly as

I(t)=ηΓmeas ⟨Xm⟩t+ξ(t),I(t) = \sqrt{\eta\Gamma_{\mathrm{meas}}}\, \langle X_m\rangle_t + \xi(t),

where η\eta is the detection efficiency, Γmeas\Gamma_{\mathrm{meas}} is a convention-dependent measurement rate, Xm=b+b†X_m=b+b^\dagger, and ξ(t)\xi(t) is measurement noise. The conditional mechanical state follows a stochastic master equation when the record is retained. If the record is discarded, the same interaction produces unconditional diffusion and decoherence.

This is the platform version of Measurement Backaction and Diffusive Trajectories.

An optomechanical position measurement has two quantum noise terms:

  • imprecision noise, because the optical phase quadrature has shot noise;
  • backaction force noise, because amplitude fluctuations shake the mechanical oscillator through radiation pressure.

A compact convention-independent statement is

Sxximp(ω)SFFba(ω)≥ℏ24,S_{xx}^{\mathrm{imp}}(\omega) S_{FF}^{\mathrm{ba}}(\omega) \ge \frac{\hbar^2}{4},

for an ideal linear detector with no useful cross-correlation between imprecision and force noise. Different one-sided and two-sided spectral-density conventions move factors of 22, but not the tradeoff.

Increasing probe power usually lowers SxximpS_{xx}^{\mathrm{imp}} while raising SFFbaS_{FF}^{\mathrm{ba}}. At low power the measurement is too noisy to resolve the motion. At high power the measurement itself drives the motion. The standard quantum limit is the optimum that results when these two contributions are balanced for a specified force, displacement, or quadrature measurement.

The standard quantum limit is not a universal ban on better measurements. It assumes a particular observable, detector model, and absence of helpful correlations. Backaction-evading measurements, squeezed input light, variational readout, two-tone schemes, and quantum nondemolition quadrature measurements can beat the simplest SQL for selected tasks, while still respecting the full quantum noise constraints.

The mechanical oscillator modifies the intracavity field with a delay set by κ−1\kappa^{-1}. That delayed radiation-pressure force changes the mechanical damping and spring constant. In linear response one writes an effective mechanical susceptibility

χeff−1(ω)=χm−1(ω)−Σopt(ω),\chi_{\mathrm{eff}}^{-1}(\omega) = \chi_m^{-1}(\omega) - \Sigma_{\mathrm{opt}}(\omega),

where Σopt\Sigma_{\mathrm{opt}} is an optical self-energy. Its real part shifts the mechanical frequency; its imaginary part changes the damping.

This is called dynamical backaction. It can cool a mechanical mode, amplify it, drive self-oscillation, or entangle it with the optical field. It is backaction in the literal dynamical sense, not merely a philosophical synonym for measurement disturbance.

For a red-detuned drive near

Δ≃−Ωm,\Delta \simeq - \Omega_m,

the dominant resonant process converts a mechanical phonon plus a drive photon into a higher-frequency cavity photon. This removes mechanical energy. In the resolved-sideband regime,

Ωm>κ,\Omega_m \gt \kappa,

the anti-Stokes scattering process can be stronger than the Stokes heating process.

With the detuning convention Δ=ωd−ωc\Delta=\omega_d-\omega_c, useful approximate rates are

A−=∣G∣2κ(κ/2)2+(Δ+Ωm)2,A+=∣G∣2κ(κ/2)2+(Δ−Ωm)2.A_- = \frac{|G|^2\kappa}{ (\kappa/2)^2+(\Delta+\Omega_m)^2 }, \qquad A_+ = \frac{|G|^2\kappa}{ (\kappa/2)^2+(\Delta-\Omega_m)^2 }.

Here A−A_- removes phonons and A+A_+ adds phonons. The optical damping rate is

Γopt=A−−A+.\Gamma_{\mathrm{opt}} = A_- - A_+ .

The approximate final occupation is

nˉfinal≃Γmnˉm+A+Γm+Γopt.\bar n_{\mathrm{final}} \simeq \frac{ \Gamma_m\bar n_m + A_+ }{ \Gamma_m+\Gamma_{\mathrm{opt}} }.

Ground-state cooling requires more than a large optical damping rate. One also needs sufficiently low classical heating, a sideband-resolved or otherwise optimized spectrum, stable drive power, and a bath that does not add excess noise faster than the cooling removes phonons.

Beam-Splitter and Two-Mode-Squeezing Limits

Section titled “Beam-Splitter and Two-Mode-Squeezing Limits”

Under a rotating-wave approximation, the red-detuned interaction becomes approximately

Hbsℏ=−(Gd†b+G∗db†),\frac{H_{\mathrm{bs}}}{\hbar} = - \left( Gd^\dagger b + G^*db^\dagger \right),

which swaps excitations between light and mechanics. This is the interaction behind sideband cooling, state transfer, and optomechanically induced transparency.

For a blue-detuned drive near Δ≃+Ωm\Delta\simeq+\Omega_m, the resonant interaction is instead

Htmsℏ=−(Gd†b†+G∗db),\frac{H_{\mathrm{tms}}}{\hbar} = - \left( Gd^\dagger b^\dagger + G^*db \right),

which creates correlated photon-phonon pairs. It can amplify motion, generate squeezing, or become unstable when gain exceeds damping.

These approximations are powerful only when the rotating terms are well separated from the counter-rotating terms. Outside the resolved-sideband regime, both processes can contribute appreciably.

The mechanical mode is rarely coupled to only one clean thermal bath. Practical noise sources include:

  • substrate phonons and support losses;
  • gas collisions or clamping loss;
  • absorption-induced heating;
  • photothermal forces;
  • frequency noise of the optical or microwave drive;
  • charge, dielectric, or two-level fluctuators;
  • recoil heating in levitated systems.

The open-system task is to identify which noise source contributes at which frequency and through which operator. A force-noise spectrum near Ωm\Omega_m drives energy exchange. Low-frequency frequency noise causes phase diffusion of the mechanical oscillator. Nonlinear damping and mode coupling can make the effective bath depend on amplitude.

The relation between damping and thermal noise is governed by fluctuation–dissipation assumptions when the bath is in equilibrium. See Fluctuation–Dissipation Relation and Caldeira–Leggett Model for the canonical background.

Different optomechanical platforms emphasize different limits:

  • optical microcavities often have strong gradients, low optical thermal occupation, and challenging absorption heating;
  • microwave drumhead circuits naturally interface with superconducting electronics and microwave input–output measurement;
  • membrane-in-the-middle and Fabry–Perot systems provide clean textbook radiation-pressure geometry;
  • levitated particles can have extremely weak mechanical contact with a substrate but require careful control of recoil, gas, and trap noise;
  • phononic-crystal devices engineer mechanical and optical localization together.

The same symbols g0g_0, κ\kappa, Γm\Gamma_m, Ωm\Omega_m, and nˉm\bar n_m do not by themselves determine the experiment. Geometry, bath temperature, drive noise, optical absorption, coupling efficiency, and detection efficiency decide which open-system model is trustworthy.

  • Treating a homodyne record as a nondisturbing readout of a preexisting trajectory.
  • Quoting g0g_0 without also quoting κ\kappa, Γm\Gamma_m, Ωm\Omega_m, drive strength, and thermal occupation.
  • Calling any red-detuned cooling “ground-state cooling” without checking the final occupation and added heating.
  • Confusing optical damping with ordinary mechanical damping from the support.
  • Forgetting that increasing probe power can increase radiation-pressure backaction and absorption heating.
  • Using the sideband-cooling rate formulas outside their weak-coupling, linearized, and Markovian limits.
  • Treating the standard quantum limit as a universal impossibility theorem rather than a detector- and task-dependent bound.

Suppose a cavity resonance changes with displacement as ωc(x)=ωc(0)+Gxx\omega_c(x)=\omega_c(0)+G_xx. Express g0g_0 in terms of GxG_x and xzpfx_{\mathrm{zpf}} using the sign convention in this page.

Solution

The page uses

g0=−xzpf∂ωc∂x.g_0 = - x_{\mathrm{zpf}} \frac{\partial\omega_c}{\partial x}.

Since ∂ωc/∂x=Gx\partial\omega_c/\partial x=G_x,

g0=−xzpfGx.g_0 = - x_{\mathrm{zpf}}G_x .

Only the relative sign matters after the phases of aa and bb are fixed; rates depend on ∣g0∣|g_0| or ∣G∣|G|.

Using the sideband-cooling rates in this page, show which rate is resonantly enhanced when Δ=−Ωm\Delta=-\Omega_m.

Solution

The cooling rate is

A−=∣G∣2κ(κ/2)2+(Δ+Ωm)2.A_- = \frac{|G|^2\kappa}{ (\kappa/2)^2+(\Delta+\Omega_m)^2 }.

At Δ=−Ωm\Delta=-\Omega_m, the second term in the denominator vanishes, so A−A_- is maximized:

A−=4∣G∣2κ.A_- = \frac{4|G|^2}{\kappa}.

The heating rate is

A+=∣G∣2κ(κ/2)2+(Δ−Ωm)2=∣G∣2κ(κ/2)2+4Ωm2.A_+ = \frac{|G|^2\kappa}{ (\kappa/2)^2+(\Delta-\Omega_m)^2 } = \frac{|G|^2\kappa}{ (\kappa/2)^2+4\Omega_m^2 }.

Thus red detuning favors anti-Stokes scattering when the sidebands are spectrally resolved.

Assume A+≪A−A_+\ll A_- and Γopt=A−−A+≫Γm\Gamma_{\mathrm{opt}}=A_--A_+\gg\Gamma_m. What does the approximate final occupation reduce to?

Solution

Start from

nˉfinal≃Γmnˉm+A+Γm+Γopt.\bar n_{\mathrm{final}} \simeq \frac{ \Gamma_m\bar n_m + A_+ }{ \Gamma_m+\Gamma_{\mathrm{opt}} }.

If Γopt≫Γm\Gamma_{\mathrm{opt}}\gg\Gamma_m, the denominator is approximately Γopt\Gamma_{\mathrm{opt}}. If A+≪A−A_+\ll A_-, then Γopt≃A−\Gamma_{\mathrm{opt}}\simeq A_-. Therefore

nˉfinal≃ΓmnˉmA−+A+A−.\bar n_{\mathrm{final}} \simeq \frac{\Gamma_m\bar n_m}{A_-} + \frac{A_+}{A_-}.

The first term is residual heating from the original mechanical bath. The second term is the quantum backaction limit from Stokes scattering.

In a simplified detector model, suppose Sxximp=C/PS_{xx}^{\mathrm{imp}}=C/P and SFFba=DPS_{FF}^{\mathrm{ba}}=DP, where PP is probe power and C,D>0C,D\gt0. Find the power-independent product and explain why there is an optimum total noise.

Solution

The product is

SxximpSFFba=CPDP=CD.S_{xx}^{\mathrm{imp}}S_{FF}^{\mathrm{ba}} = \frac{C}{P}DP = CD .

Increasing PP lowers imprecision but raises backaction. For a mechanical displacement or force measurement, the total noise contains an imprecision contribution plus a backaction contribution filtered by the mechanical susceptibility. Since one decreases and the other increases with PP, their sum has an optimum.

  • V. B. Braginsky and F. Y. Khalili, Quantum Measurement, Cambridge University Press, 1992.
  • T. J. Kippenberg and K. J. Vahala, “Cavity optomechanics: back-action at the mesoscale,” Science 321, 1172-1176 (2008).
  • A. A. Clerk, M. H. Devoret, S. M. Girvin, F. Marquardt, and R. J. Schoelkopf, “Introduction to quantum noise, measurement, and amplification,” Reviews of Modern Physics 82, 1155-1208 (2010).
  • J. D. Teufel et al., “Sideband cooling of micromechanical motion to the quantum ground state,” Nature 475, 359-363 (2011).
  • J. Chan et al., “Laser cooling of a nanomechanical oscillator into its quantum ground state,” Nature 478, 89-92 (2011).
  • M. Aspelmeyer, T. J. Kippenberg, and F. Marquardt, “Cavity optomechanics,” Reviews of Modern Physics 86, 1391-1452 (2014).