Skip to content

Quantum Optics

Quantum optics is the natural laboratory for open quantum systems. Photons leave cavities, atoms emit into radiation modes, detectors produce clicks and currents, and measured output fields condition the state of the source. The same experiment can be described as a master equation, a Langevin equation, an input–output relation, or a quantum trajectory depending on which degrees of freedom and records are retained.

This page is an application map. It does not try to replace the full quantum-optics treatment or the AMO treatment of light–matter interaction. Its goal is to show how the measurement and open-system language used throughout this volume appears in optical systems.

A typical quantum-optical setup has:

  • a localized system: atom, ion, molecule, cavity mode, mechanical mode, or nonlinear optical mode;
  • one or more electromagnetic reservoirs;
  • coherent drives supplied by lasers;
  • detectors that monitor selected output modes;
  • unmonitored loss channels that produce decoherence and damping.

The unconditional state often obeys a quantum optical master equation,

ρ˙=−iℏ[H,ρ]+∑jD[Lj]ρ,\dot\rho = - \frac{i}{\hbar}[H,\rho] + \sum_j \mathcal D[L_j]\rho,

where

D[L]ρ=LρL†−12{L†L,ρ}.\mathcal D[L]\rho = L\rho L^\dagger - \frac12\{L^\dagger L,\rho\}.

The collapse operators LjL_j are not just mathematical terms. In quantum optics they often correspond directly to physical output channels: spontaneous emission, cavity leakage, absorption, scattering, or detector-monitored fluorescence.

For a single cavity mode with annihilation operator aa and loss rate κ\kappa, zero-temperature damping is described by

ρ˙=−iℏ[H,ρ]+κD[a]ρ.\dot\rho = - \frac{i}{\hbar}[H,\rho] + \kappa\mathcal D[a]\rho.

If the reservoir has thermal occupation nˉ\bar n, the damping terms become

κ(nˉ+1)D[a]ρ+κnˉ D[a†]ρ.\kappa(\bar n+1)\mathcal D[a]\rho + \kappa\bar n\,\mathcal D[a^\dagger]\rho.

At optical frequencies and ordinary laboratory temperatures, nˉ\bar n is usually negligible. At microwave frequencies it may not be. This distinction matters when comparing optical cavity QED with circuit QED.

The same cavity coupling appears in input–output theory:

bout(t)=bin(t)+κ a(t),b_{\mathrm{out}}(t) = b_{\mathrm{in}}(t) + \sqrt{\kappa}\,a(t),

up to phase convention. The output field is what detectors actually see. See Input–Output Theory.

For a two-level atom with lowering operator σ−=∣g⟩⟨e∣\sigma_-=\lvert g\rangle\langle e\rvert, spontaneous emission into unobserved modes is modeled by

ρ˙=−iℏ[H,ρ]+ΓD[σ−]ρ.\dot\rho = - \frac{i}{\hbar}[H,\rho] + \Gamma\mathcal D[\sigma_-]\rho.

The excited-state population decays as

ρ˙ee=−Γρee,\dot\rho_{ee} = -\Gamma\rho_{ee},

while the optical coherence decays at rate Γ/2\Gamma/2 before adding any extra pure dephasing. This is the optical prototype behind the Amplitude Damping Master Equation.

Adding a coherent near-resonant drive gives the standard Optical Bloch Equations for saturation, resonance fluorescence, and damped Rabi dynamics.

If emitted photons are monitored, the same channel can be unraveled into quantum jumps. If the detector sees a photon, the conditioned state updates by

ρc⟼LρcL†Tr⁡(L†Lρc).\rho_c \longmapsto \frac{ L\rho_c L^\dagger }{ \operatorname{Tr}(L^\dagger L\rho_c) }.

If no photon is detected, the state evolves under a non-Hermitian no-click evolution plus normalization. Averaging over all records returns the unconditional master equation.

For a monitored channel with operator LL, the probability of a click in a short interval dtdt is

pclick=Tr⁡(L†Lρc) dt.p_{\mathrm{click}} = \operatorname{Tr}(L^\dagger L\rho_c)\,dt.

Direct detection is therefore naturally modeled by a Poisson increment dNtdN_t satisfying

dNt2=dNt,E[dNt∣ρc]=Tr⁡(L†Lρc) dt.dN_t^2=dN_t, \qquad \mathbb E[dN_t\mid\rho_c] = \operatorname{Tr}(L^\dagger L\rho_c)\,dt.

Photon counting is the operational basis for antibunching, resonance fluorescence records, photon-correlation measurements, and many single-emitter experiments. See Photon Counting and Quantum Jump Trajectories.

Homodyne detection mixes the output field with a strong local oscillator and measures one field quadrature. For a monitored channel cc and detection efficiency η\eta, a common normalized stochastic master equation is

dρc=Lρc dt+η H[e−iϕc]ρc dWt,d\rho_c = \mathcal L\rho_c\,dt + \sqrt{\eta}\, \mathcal H[e^{-i\phi}c]\rho_c\,dW_t,

where ϕ\phi is the local-oscillator phase and

H[c]ρ=cρ+ρc†−Tr⁡[(c+c†)ρ]ρ.\mathcal H[c]\rho = c\rho+\rho c^\dagger - \operatorname{Tr}[(c+c^\dagger)\rho]\rho.

The measured current has the schematic form

dYt=η ⟨e−iϕc+eiϕc†⟩c dt+dWt.dY_t = \sqrt{\eta}\, \langle e^{-i\phi}c+e^{i\phi}c^\dagger\rangle_c\,dt + dW_t.

Heterodyne detection measures both quadratures at the cost of added vacuum noise. Both descriptions are diffusive rather than jump-like; the measurement record is a noisy continuous current. See Homodyne Detection, Heterodyne Detection, and Diffusive Trajectories.

A coherent optical drive is often treated as a classical amplitude plus quantum vacuum fluctuations:

bin(t)=βin(t)+ξin(t).b_{\mathrm{in}}(t) = \beta_{\mathrm{in}}(t) + \xi_{\mathrm{in}}(t).

For a driven cavity, the mean field obeys a damped oscillator equation. In a frame rotating at the drive frequency,

α˙=−(iΔc+κ2)α+E,\dot\alpha = - \left( i\Delta_c+\frac{\kappa}{2} \right)\alpha + \mathcal E,

where Δc\Delta_c is the drive-cavity detuning and E\mathcal E is a drive amplitude in the chosen convention.

Coherent states are pointer-like for linear loss: damping maps a coherent state to another coherent state with decaying amplitude. In covariance language this is a Gaussian attenuation channel; see Gaussian Channels for the open-system convention. This is one reason classical optical fields emerge naturally from lossy quantum optical modes. Nonlinearities, single-photon sources, photon blockade, and conditional measurements are what make the field strongly nonclassical.

For a single mode, define the quadrature

Xθ=12(ae−iθ+a†eiθ).X_\theta = \frac12 \left( ae^{-i\theta} + a^\dagger e^{i\theta} \right).

In the vacuum state,

(ΔXθ)2=14(\Delta X_\theta)^2 = \frac14

for every θ\theta. A squeezed state has reduced noise in one quadrature,

(ΔXθ0)2<14,(\Delta X_{\theta_0})^2 \lt \frac14,

with increased noise in the conjugate quadrature. Squeezed light is therefore both a nonclassical state and a practical measurement resource: it can reduce imprecision in one quadrature while increasing backaction or noise elsewhere.

Open-system descriptions of squeezed reservoirs modify the usual damping master equation by adding phase-sensitive noise correlations. The physical message is that “the bath” need not be thermal or vacuum; engineered optical reservoirs can carry squeezing, correlations, and directionality.

Quantum optics often characterizes output fields through normally ordered correlation functions. The second-order coherence is

g(2)(τ)=⟨a†(0)a†(τ)a(τ)a(0)⟩⟨a†a⟩2g^{(2)}(\tau) = \frac{ \langle a^\dagger(0)a^\dagger(\tau)a(\tau)a(0) \rangle }{ \langle a^\dagger a\rangle^2 }

for a stationary single mode. Standard reference values are:

Fieldg(2)(0)g^{(2)}(0)Interpretation
coherent state11Poissonian photon statistics
thermal state22bunching
ideal single-photon stream00antibunching

The quantum regression theorem connects such multitime correlations to the same Liouvillian used in the master equation, under Markovian assumptions. When those assumptions fail, correlation functions can reveal memory, filtering, detector bandwidth, and mode mismatch.

Quantum optics forces a useful distinction:

  • the intracavity or atomic state is a model variable;
  • the output field is the propagating quantum system that reaches the detector;
  • the detector record is a classical stochastic process derived from the measured field.

Confusing these levels leads to mistakes. A click record is not the same object as the intracavity photon number. A homodyne current is not a projective measurement of a field quadrature at a single instant. Loss into unmonitored modes is still a measurement by the environment, even if no experimentalist records the outcome.

  • Treating a laser drive as if it removed vacuum fluctuations from the input field.
  • Forgetting detector efficiency and dark counts when interpreting photon records.
  • Confusing normally ordered optical correlations with symmetrized noise spectra.
  • Assuming every master-equation collapse operator is directly monitored.
  • Treating quantum jumps as unique physical histories rather than one unraveling selected by the measurement scheme.
  • Ignoring finite detector bandwidth and mode filtering.
  • Applying the quantum regression theorem outside its Markovian domain.
  • Calling any sub-Poissonian statistic “squeezing”; squeezing refers to quadrature variance.

For a monitored spontaneous-emission channel L=Γσ−L=\sqrt{\Gamma}\sigma_-, write the conditioned state immediately after a detected photon.

Solution

The jump update is

ρc⟼LρcL†Tr⁡(L†Lρc).\rho_c \longmapsto \frac{ L\rho_c L^\dagger }{ \operatorname{Tr}(L^\dagger L\rho_c) }.

Substituting L=Γσ−L=\sqrt{\Gamma}\sigma_- gives

ρc⟼σ−ρcσ+Tr⁡(σ+σ−ρc).\rho_c \longmapsto \frac{ \sigma_-\rho_c\sigma_+ }{ \operatorname{Tr}(\sigma_+\sigma_-\rho_c) }.

For an ideal two-level atom, this leaves the atom in the ground state whenever the denominator is nonzero.

A cavity mode obeys ρ˙=κD[a]ρ\dot\rho=\kappa\mathcal D[a]\rho with no drive. What happens to an initial coherent amplitude α(0)\alpha(0)?

Solution

The mean amplitude obeys

ddt⟨a⟩=−κ2⟨a⟩.\frac{d}{dt}\langle a\rangle = - \frac{\kappa}{2} \langle a\rangle.

Therefore

α(t)=α(0)e−κt/2.\alpha(t) = \alpha(0)e^{-\kappa t/2}.

For linear damping into vacuum, an initial coherent state remains coherent while its amplitude decays. The photon number decays as e−κte^{-\kappa t}.

Why does an ideal single-photon source have g(2)(0)=0g^{(2)}(0)=0?

Solution

The quantity g(2)(0)g^{(2)}(0) measures the normalized probability of detecting two photons at the same time. An ideal single-photon source emits at most one photon in the relevant mode at a time, so the two-photon coincidence probability at zero delay vanishes:

g(2)(0)=0.g^{(2)}(0)=0.

This is antibunching and has no classical analogue for an ordinary positive-intensity random field.

In a homodyne measurement with channel cc, what quadrature appears in the measurement signal when the local-oscillator phase is ϕ\phi?

Solution

The signal term is proportional to

⟨e−iϕc+eiϕc†⟩c.\langle e^{-i\phi}c + e^{i\phi}c^\dagger \rangle_c.

Changing ϕ\phi rotates the measured quadrature. The stochastic noise term remains a Wiener increment in the ideal Markovian model.

  • R. J. Glauber, “The quantum theory of optical coherence,” Physical Review 130, 2529-2539 (1963).
  • M. O. Scully and M. S. Zubairy, Quantum Optics, Cambridge University Press (1997).
  • D. F. Walls and G. J. Milburn, Quantum Optics, Springer (2008).
  • C. W. Gardiner and P. Zoller, Quantum Noise, Springer (2004).
  • H. J. Carmichael, An Open Systems Approach to Quantum Optics, Springer (1993).
  • J. Dalibard, Y. Castin, and K. Molmer, “Wave-function approach to dissipative processes in quantum optics,” Physical Review Letters 68, 580-583 (1992).
  • H. M. Wiseman and G. J. Milburn, Quantum Measurement and Control, Cambridge University Press (2010).
  • H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford University Press (2002).