Skip to content

Photon Counting

Photon counting is direct detection of quanta in an optical, microwave, or other bosonic output channel. Instead of a continuous quadrature current, the record is a sequence of events: detector clicks, time tags, or counts in finite bins. Calibrated optical number-resolving and threshold POVMs, efficiency, dark counts, recovery, and detector tomography are developed on the quantum-optics Photon Counting page. This page takes the detector model as given and concentrates on the conditional record produced by monitoring an open-system output.

In the ideal Markovian limit, photon counting is modeled by a counting process NtN_t with increment

dNt∈{0,1},dNt2=dNt.dN_t\in\{0,1\}, \qquad dN_t^2=dN_t.

The conditional mean of the increment is set by the instantaneous detected flux. For an output channel with coupling operator LL,

E[dNt∣ρc(t)]=η ⟨L†L⟩c dt\mathbb E[dN_t|\rho_c(t)] = \eta\, \langle L^\dagger L\rangle_c\,dt

in a common ideal convention, where η\eta is the total detection efficiency. This page explains the record and detector model. The conditioned state dynamics associated with clicks and no-click intervals is the canonical subject of Quantum Jump Trajectories.

Photon counting is sensitive to intensity, not field quadrature. In input–output language, a monitored port satisfies

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

For vacuum input with known coherent amplitudes subtracted or accounted for, the emitted photon flux associated with the system channel is

Φ(t)=⟨L†L⟩.\Phi(t) = \langle L^\dagger L\rangle.

The operator LL is often a cavity-loss operator such as κa\sqrt{\kappa}a, or a spontaneous-emission operator such as Γ σ−\sqrt{\Gamma}\,\sigma_-. The same LL contributes the unconditional dissipator D[L]ρ\mathcal D[L]\rho when the output is not retained.

Photon counting asks a different question from homodyne detection. It samples emitted quanta. Homodyne detection interferes the output with a strong local oscillator and samples a quadrature. Both may monitor the same physical port, but they produce different records and different conditioned trajectories.

The cumulative count

NtN_t

is a nondecreasing integer-valued process. In a small interval [t,t+dt][t,t+dt], dNt=1dN_t=1 means one detected event and dNt=0dN_t=0 means no detected event. For a time bin of finite width Δt\Delta t,

ΔNk=Ntk+1−Ntk\Delta N_k = N_{t_{k+1}}-N_{t_k}

may be 00, 11, or larger if the bin is wide enough or the flux is high enough.

Equivalently, an ideal time-tagging detector reports event times

0<t1<t2<⋯<tm≤T.0\lt t_1\lt t_2\lt\cdots\lt t_m\le T.

The two descriptions are related by binning. Fine time tags preserve more information than coarse bins, but real detector timing jitter and dead time set a practical limit on time resolution.

The conditional intensity or hazard rate λt\lambda_t is the rate predicted from the current conditional state:

E[dNt∣ρc(t)]=λt dt.\mathbb E[dN_t|\rho_c(t)] = \lambda_t\,dt.

For one ideal detected channel,

λt=η Tr⁡(L†Lρc(t)).\lambda_t = \eta\,\operatorname{Tr}(L^\dagger L\rho_c(t)).

With dark counts, a simple model is

λt=η Tr⁡(L†Lρc(t))+λd,\lambda_t = \eta\,\operatorname{Tr}(L^\dagger L\rho_c(t)) + \lambda_{\mathrm d},

where λd\lambda_{\mathrm d} is a background count rate. The dark-count branch carries little or no information about the system, so a click should not always be interpreted as a system quantum being emitted.

If several detectors monitor distinct channels LjL_j, the record has several counting processes Nj,tN_{j,t} with intensities

λj,t=ηj Tr⁡(Lj†Ljρc(t))+λj,d.\lambda_{j,t} = \eta_j\, \operatorname{Tr}(L_j^\dagger L_j\rho_c(t)) + \lambda_{j,\mathrm d}.

The detector label is part of the record. Losing it is a coarse graining that changes the conditioned state update.

A click is an outcome, but so is no click. In an interval dtdt, the two ideal alternatives are:

click: dN_t = 1
no click: dN_t = 0

If a click occurs in channel LL, the normalized state update has the jump form

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

If no click occurs, the state is also updated because the absence of an expected event is information. The no-click update involves the effective non-Hermitian Hamiltonian and normalization by the probability of no event. The detailed formulas, waiting-time distributions, and ensemble recovery are treated in Quantum Jump Trajectories.

The important record-level lesson is that a data file containing no clicks is not the same as no measurement. It may be strong evidence that the emitting system was dark, in the ground state, outside the collection mode, or hidden by detector inefficiency.

Real photon counters require more than the ideal point-process model. Common effects include:

  • finite efficiency η\eta;
  • dark counts or background light;
  • dead time after a detection event;
  • timing jitter;
  • finite time-bin width;
  • afterpulsing or detector recovery artifacts;
  • number resolution or lack of number resolution;
  • saturation at high flux;
  • loss before the detector;
  • mode mismatch and imperfect filtering.

Some imperfections can be represented by adding unmonitored loss channels or background counting branches. Others, such as dead time and afterpulsing, create detector memory and can invalidate a simple Markovian counting-process model unless the detector state is included explicitly.

Ideal photodetection is naturally connected with normally ordered field correlations. For a positive-frequency field operator E(+)(t)E^{(+)}(t), the single-time counting rate is proportional to

G(1)(t,t)=⟨E(−)(t)E(+)(t)⟩.G^{(1)}(t,t) = \left\langle E^{(-)}(t)E^{(+)}(t) \right\rangle.

Two-count correlations involve

G(2)(t1,t2)=⟨E(−)(t1)E(−)(t2)E(+)(t2)E(+)(t1)⟩,G^{(2)}(t_1,t_2) = \left\langle E^{(-)}(t_1)E^{(-)}(t_2) E^{(+)}(t_2)E^{(+)}(t_1) \right\rangle,

up to detector response and normalization conventions. This is the origin of photon-correlation measurements such as photon antibunching and Hanbury Brown–Twiss interferometry.

In Markovian input–output notation, E(+)E^{(+)} is represented by the relevant output field. For a vacuum-input system channel, the normally ordered flux reduces to the system expectation ⟨L†L⟩\langle L^\dagger L\rangle after the input-output relation is applied.

This page only gives the preview. A quantum-optics treatment owns the full theory of optical coherence functions.

Prepared optical Fock-state distributions, loss-thinned number statistics, phase properties, and source-quality diagnostics are developed on Photon Number States.

For a two-level atom with spontaneous emission operator

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

the ideal click rate is

λt=Γ⟨σ+σ−⟩c=Γpe(t).\lambda_t = \Gamma \langle\sigma_+\sigma_-\rangle_c = \Gamma p_e(t).

A click indicates an emitted photon and updates the atom toward the ground state. A long no-click interval makes excitation less likely, provided the detector efficiency is high and background counts are low.

For a leaking cavity with

L=κa,L=\sqrt{\kappa}a,

the ideal count rate is

λt=κ⟨a†a⟩c.\lambda_t = \kappa\langle a^\dagger a\rangle_c.

The count record probes the output photon flux, not a nondestructive record of the intracavity photon number at every instant. Cavity loss removes photons from the mode.

  • Treating no-click intervals as absence of measurement rather than no-click outcomes.
  • Interpreting every click as a system photon when dark counts or background light are present.
  • Forgetting that inefficiency leaves unobserved emissions that still decohere the system.
  • Using a jump trajectory when the experiment actually measures a homodyne or heterodyne current.
  • Confusing detected output flux with intracavity photon number.
  • Ignoring detector dead time at high count rates.
  • Applying infinitesimal dNt∈{0,1}dN_t\in\{0,1\} rules to large time bins without allowing multiple counts.
  • Reading physical meaning into a collapse operator without specifying the monitored port and detector.
  • R. J. Glauber, “The quantum theory of optical coherence,” Physical Review 130, 2529–2539, 1963.
  • L. Mandel and E. Wolf, Optical Coherence and Quantum Optics, Cambridge University Press, 1995.
  • H. J. Carmichael, An Open Systems Approach to Quantum Optics, Springer, 1993.
  • 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.
  • M. B. Plenio and P. L. Knight, “The quantum-jump approach to dissipative dynamics in quantum optics,” Reviews of Modern Physics 70, 101–144, 1998.
  1. For a constant counting intensity λ\lambda, compute the probability of no clicks in time TT.
Solution

For a Poisson process with constant rate λ\lambda, the no-count probability is

P0(T)=e−λT.P_0(T)=e^{-\lambda T}.

Equivalently, multiplying no-click probabilities over short bins gives

lim⁡N→∞(1−λTN)N=e−λT.\lim_{N\to\infty} \left( 1-\lambda\frac{T}{N} \right)^N = e^{-\lambda T}.
  1. A two-level atom has L=Γσ−L=\sqrt{\Gamma}\sigma_- and conditional excited-state population pep_e. With detection efficiency η\eta and dark-count rate λd\lambda_{\mathrm d}, write the observed click intensity.
Solution

Since

L†L=Γσ+σ−,L^\dagger L = \Gamma\sigma_+\sigma_-,

one has

Tr⁡(L†Lρc)=Γpe.\operatorname{Tr}(L^\dagger L\rho_c) = \Gamma p_e.

The observed intensity in the simple dark-count model is

λobs=ηΓpe+λd.\lambda_{\mathrm{obs}} = \eta\Gamma p_e + \lambda_{\mathrm d}.
  1. Suppose a detector reports only the total count from two unresolved channels L1L_1 and L2L_2. What information is lost compared with a detector that labels which channel clicked?
Solution

The labeled record distinguishes the alternatives L1ρL1†L_1\rho L_1^\dagger and L2ρL2†L_2\rho L_2^\dagger. If the detector reports only the total count, the click operation is the coarse-grained sum

Iclick(ρ)=L1ρL1†+L2ρL2†,\mathcal I_{\mathrm{click}}(\rho) = L_1\rho L_1^\dagger + L_2\rho L_2^\dagger,

with the appropriate rates absorbed into L1L_1 and L2L_2. Coarse graining loses information about which channel carried the emission and can leave the conditional state more mixed.

  1. Why is a large-bin count ΔNk\Delta N_k not always well represented by an infinitesimal increment dNtdN_t?
Solution

The infinitesimal rule keeps only terms through first order in dtdt, so the probability of two or more counts is O(dt2)O(dt^2) and is neglected. In a finite bin Δt\Delta t, especially at high flux, multiple events can have appreciable probability. Then ΔNk\Delta N_k should be treated as a finite-bin count with Poisson or detector-specific statistics, not as a binary increment.