Skip to content

Homodyne Detection

Homodyne detection measures one quadrature of an optical or microwave field by interfering the signal with a strong phase reference called the local oscillator. In balanced homodyne detection, the two outputs of a beam splitter are sent to photodetectors and the photocurrents are subtracted. The large local-oscillator amplitude converts a field quadrature into an intensity difference.

The lossless two-port unitary and its phase conventions are developed at Beam Splitters. This page takes that mixer as part of the detector and owns the resulting continuous measurement record and conditional dynamics. Homodyne and Heterodyne Detection owns the practical optical receiver, mode selection, vacuum calibration, quadrature distributions, and laboratory nonidealities.

In continuous-measurement language, homodyne detection produces a diffusive record:

dYt=conditional quadrature signal dt+dWt.dY_t = \text{conditional quadrature signal}\,dt + dW_t.

The Wiener increment dWtdW_t is not a technical afterthought. It is the shot-noise or vacuum-noise part of the record, and it is what makes homodyne-conditioned states follow stochastic trajectories rather than deterministic expectation-value curves.

For a single mode with annihilation operator aa, define a quadrature convention

Xϕ=e−iϕa+eiϕa†2.X_\phi = \frac{ e^{-i\phi}a+e^{i\phi}a^\dagger }{\sqrt2}.

The phase ϕ\phi is the local-oscillator phase. Changing ϕ\phi rotates the measured quadrature:

Xϕ+π/2=−ie−iϕa+ieiϕa†2,X_{\phi+\pi/2} = \frac{ -i e^{-i\phi}a+i e^{i\phi}a^\dagger }{\sqrt2},

which is conjugate to XϕX_\phi in the usual single-mode normalization.

Continuous-output records often absorb factors of 2\sqrt2 into the definition of dYtdY_t and the coupling operator. Always check the convention before comparing formulas. The operational statement is convention independent: a homodyne detector measures one phase-selected quadrature of the output field.

Treat the local oscillator as a large coherent amplitude

αLO=∣αLO∣eiϕ.\alpha_{\mathrm{LO}} = |\alpha_{\mathrm{LO}}|e^{i\phi}.

Let asa_s be the signal mode entering the other input port of a balanced beam splitter. With one common phase convention, the two output modes are

a+=as+αLO2,a−=as−αLO2.a_+ = \frac{a_s+\alpha_{\mathrm{LO}}}{\sqrt2}, \qquad a_- = \frac{a_s-\alpha_{\mathrm{LO}}}{\sqrt2}.

The photon-number difference is

n+−n−=a+†a+−a−†a−=αLO∗as+αLOas†.\begin{aligned} n_+ - n_- &= a_+^\dagger a_+ - a_-^\dagger a_- \\ &= \alpha_{\mathrm{LO}}^*a_s + \alpha_{\mathrm{LO}}a_s^\dagger . \end{aligned}

Dividing by the local-oscillator amplitude gives

n+−n−∣αLO∣=e−iϕas+eiϕas†.\frac{n_+ - n_-}{|\alpha_{\mathrm{LO}}|} = e^{-i\phi}a_s + e^{i\phi}a_s^\dagger .

Thus the difference current measures the signal quadrature selected by ϕ\phi. The sum current is dominated by the local-oscillator intensity and is usually used for normalization and diagnostics rather than as the quantum signal.

This derivation also explains why balance matters. Subtracting the two detector outputs cancels the common large local-oscillator intensity to leading order, while preserving the interference term linear in the signal field.

For a Markovian output port, Input–Output Theory writes

bout(t)=bin(t)+L(t),b_{\mathrm{out}}(t) = b_{\mathrm{in}}(t)+L(t),

where LL is the system coupling operator for that port. With vacuum input and known coherent amplitudes subtracted, a common homodyne record convention is

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

Here η\eta is the total detection efficiency, including collection loss, propagation loss, mode mismatch, detector efficiency, and any added-noise penalty that is represented as loss in the model. The conditional mean is the output quadrature signal; the innovation is

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

The record is a property of the traveling output field. It is usually not correct to replace boutb_{\mathrm{out}} by an intracavity operator without the input–output relation and the relevant phase convention.

The same monitored output channel gives a diffusive stochastic master equation. In a common normalized convention,

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

where

H[C]ρ=Cρ+ρC†−Tr⁡ ⁣[(C+C†)ρ]ρ.\mathcal H[C]\rho = C\rho+\rho C^\dagger - \operatorname{Tr} \!\left[ (C+C^\dagger)\rho \right]\rho.

The unconditional generator L\mathcal L contains the Hamiltonian and all monitored and unmonitored Lindblad channels. Averaging over homodyne records recovers the unconditional master equation.

This page explains the detection model. The general stochastic-equation framework is the canonical home for Stochastic Master Equations, and the broader family of noisy-current trajectories is treated in Diffusive Trajectories.

Photon counting measures output intensity by registering discrete detection events. Homodyne detection measures a field quadrature through interference with a phase reference. The difference is not cosmetic:

  • photon counting gives a point process dNtdN_t;
  • homodyne detection gives a diffusive current dYtdY_t;
  • photon counting is naturally sensitive to L†LL^\dagger L;
  • homodyne detection is sensitive to e−iϕL+eiϕL†e^{-i\phi}L+e^{i\phi}L^\dagger;
  • the same unmonitored dissipator D[L]ρ\mathcal D[L]\rho can underlie either record.

Thus a Lindblad operator LL is not automatically a detector click. How the output field is measured determines the unraveling.

Homodyne detection measures one quadrature selected by a phase reference. Heterodyne detection records both quadratures by introducing an additional mode or frequency offset, which brings extra vacuum noise. Homodyne is therefore preferable when one quadrature contains the signal of interest and the orthogonal quadrature is not needed.

In phase-space language, homodyne samples a marginal quadrature distribution, while Heterodyne Detection gives a noisy complex-amplitude record. The two methods answer different experimental questions even when they monitor the same output port.

In resonance fluorescence, homodyne detection of emitted light can monitor a dipole quadrature of a driven atom. The resulting trajectories are diffusive rather than jump-like: the state is steered by a noisy photocurrent instead of by isolated photon detections.

In cavity QED, homodyne detection of a transmitted or reflected field measures a quadrature of the cavity output. If the cavity field depends on the state of an atom or qubit, the homodyne record becomes an indirect continuous measurement of that internal degree of freedom.

In circuit QED, microwave homodyne or near-homodyne detection after amplification is the basis of dispersive qubit readout. After eliminating the resonator in an appropriate regime, the record is often written schematically as

dYt=2ηΓm ⟨σz⟩c dt+dWt.dY_t = 2\sqrt{\eta\Gamma_{\mathrm m}}\, \langle\sigma_z\rangle_c\,dt + dW_t.

This formula hides considerable hardware and calibration work: resonator response, amplifier noise, demodulation phase, integration windows, and leakage outside the simplified qubit model.

A realistic homodyne record depends on more than the ideal stochastic equation. Important experimental ingredients include:

  • local-oscillator phase and phase noise;
  • mode matching or visibility between signal and local oscillator;
  • detector quantum efficiency and propagation loss;
  • electronic noise and amplifier added noise;
  • finite detector bandwidth;
  • saturation or nonlinearity in the photodetectors or amplifiers;
  • imperfect subtraction of the two photocurrents;
  • offsets and slow drifts in the recorded current.

Many of these effects can be folded into an effective efficiency η\eta only if they behave like loss or independent Gaussian noise over the bandwidth of interest. Phase drift, detector saturation, and filtering generally require more explicit modeling.

  • Treating the homodyne current as the quadrature expectation value rather than signal plus noise.
  • Forgetting that the local-oscillator phase selects which quadrature is measured.
  • Confusing the traveling output field with an intracavity mode.
  • Comparing two formulas without checking quadrature normalization.
  • Describing homodyne detection as photon counting with smaller bins.
  • Ignoring the known coherent input amplitude before interpreting the output record.
  • Treating detector inefficiency as though it only weakens the plotted signal while leaving the conditional state equally pure.
  • Assuming a balanced detector removes quantum shot noise; it mainly removes common classical intensity noise from the local oscillator.
  • D. F. Walls and G. J. Milburn, Quantum Optics, 2nd ed., Springer, 2008.
  • C. W. Gardiner and P. Zoller, Quantum Noise, 3rd ed., Springer, 2004.
  • H. M. Wiseman and G. J. Milburn, Quantum Measurement and Control, Cambridge University Press, 2010.
  • H. J. Carmichael, An Open Systems Approach to Quantum Optics, Springer, 1993.
  • U. Leonhardt, Measuring the Quantum State of Light, Cambridge University Press, 1997.
  • H.-A. Bachor and T. C. Ralph, A Guide to Experiments in Quantum Optics, 2nd ed., Wiley-VCH, 2004.
  1. Starting from
a±=as±αLO2,a_\pm = \frac{a_s\pm\alpha_{\mathrm{LO}}}{\sqrt2},

derive the difference signal n+−n−n_+ - n_- used in balanced homodyne detection.

Solution

Compute

n+=12(as†+αLO∗)(as+αLO),n_+ = \frac12 \left( a_s^\dagger+\alpha_{\mathrm{LO}}^* \right) \left( a_s+\alpha_{\mathrm{LO}} \right),

and

n−=12(as†−αLO∗)(as−αLO).n_- = \frac12 \left( a_s^\dagger-\alpha_{\mathrm{LO}}^* \right) \left( a_s-\alpha_{\mathrm{LO}} \right).

Subtracting cancels as†asa_s^\dagger a_s and ∣αLO∣2|\alpha_{\mathrm{LO}}|^2, leaving

n+−n−=αLO∗as+αLOas†.n_+ - n_- = \alpha_{\mathrm{LO}}^*a_s + \alpha_{\mathrm{LO}}a_s^\dagger.

With αLO=∣αLO∣eiϕ\alpha_{\mathrm{LO}}=|\alpha_{\mathrm{LO}}|e^{i\phi}, the normalized difference is

n+−n−∣αLO∣=e−iϕas+eiϕas†.\frac{n_+ - n_-}{|\alpha_{\mathrm{LO}}|} = e^{-i\phi}a_s+e^{i\phi}a_s^\dagger.
  1. Show that changing the local-oscillator phase by π\pi reverses the sign of the measured quadrature.
Solution

Using

Xϕ=e−iϕa+eiϕa†2,X_\phi = \frac{ e^{-i\phi}a+e^{i\phi}a^\dagger }{\sqrt2},

one has

Xϕ+π=e−i(ϕ+π)a+ei(ϕ+π)a†2=−e−iϕa−eiϕa†2=−Xϕ.\begin{aligned} X_{\phi+\pi} &= \frac{ e^{-i(\phi+\pi)}a + e^{i(\phi+\pi)}a^\dagger }{\sqrt2} \\ &= \frac{ -e^{-i\phi}a - e^{i\phi}a^\dagger }{\sqrt2} = -X_\phi. \end{aligned}
  1. For the homodyne record
dYt=η⟨C+C†⟩cdt+dWt,dY_t = \sqrt{\eta} \langle C+C^\dagger\rangle_c dt + dW_t,

what is the conditional mean and variance of dYtdY_t?

Solution

Since E[dWt∣ρc(t)]=0\mathbb E[dW_t|\rho_c(t)]=0 and dWt2=dtdW_t^2=dt,

E[dYt∣ρc(t)]=η⟨C+C†⟩cdt,\mathbb E[dY_t|\rho_c(t)] = \sqrt{\eta} \langle C+C^\dagger\rangle_c dt,

and

Var⁡(dYt∣ρc(t))=dt.\operatorname{Var}(dY_t|\rho_c(t)) = dt.

The mean is order dtdt, while the noise scale is order dt\sqrt{dt}.

  1. Why does inefficient homodyne detection generally produce less pure conditional states than ideal homodyne detection, even when the same unconditional dissipator is present?
Solution

Efficiency less than one means that only part of the output information reaches the observer. A useful model splits the channel into a monitored part and an unmonitored loss part. The monitored part contributes to the stochastic innovation, while the unmonitored part still contributes to decoherence in the unconditional generator. The observer therefore lacks some information about the backaction, so the conditional state is generally more mixed than in an ideal measurement.