Skip to content

Repeatability and QND Measurement

Repeatability and quantum nondemolition measurement are related but not identical ideas.

A measurement is repeatable when, after obtaining an outcome, an immediate repetition of the same ideal measurement gives the same outcome with probability one. This is a statement about the post-measurement state and the next measurement.

A measurement is quantum nondemolition (QND) when an observable can be monitored over time without the measurement process randomizing the later values of that same observable. This is a dynamical statement: it involves the system Hamiltonian, the measurement interaction, environmental loss, and the time scale on which repeated predictions are claimed.

Helpful background. Lüders rule supplies immediate repeatability, while measurement backaction explains how other observables can still be disturbed.

An ideal projective measurement can be repeatable and still disturb phases, conjugate variables, or unresolved internal labels. A QND measurement is allowed to disturb other degrees of freedom; it is designed so that the disturbance avoids the monitored observable.

For the finite-dimensional repeatability argument, let a sharp measurement be represented by projectors {Πa}\{\Pi_a\}, with

ΠaΠb=δabΠa,∑aΠa=I.\Pi_a\Pi_b=\delta_{ab}\Pi_a, \qquad \sum_a\Pi_a=I.

For the ideal Lüders operation,

Ia(ρ)=ΠaρΠa,\mathcal I_a(\rho)=\Pi_a\rho\Pi_a,

the conditional state after outcome aa is

ρa=ΠaρΠaTr⁡(Πaρ).\rho_a = \frac{\Pi_a\rho\Pi_a} {\operatorname{Tr}(\Pi_a\rho)}.

A second immediate measurement of the same projective decomposition gives aa with probability one:

Tr⁡(Πaρa)=Tr⁡(ΠaΠaρΠa)Tr⁡(Πaρ)=1.\operatorname{Tr}(\Pi_a\rho_a) = \frac{\operatorname{Tr}(\Pi_a\Pi_a\rho\Pi_a)} {\operatorname{Tr}(\Pi_a\rho)} =1.

This is the simplest repeatability theorem for ideal projective measurements. It assumes no intervening dynamics, no noise, and a post-measurement system in the same Hilbert space.

For a more general outcome operation Ia\mathcal I_a, repeatability relative to the sharp question {Πa}\{\Pi_a\} requires the output of the aa operation to lie entirely in the aa subspace:

ΠaIa(ρ)Πa=Ia(ρ)\Pi_a\mathcal I_a(\rho)\Pi_a = \mathcal I_a(\rho)

for all inputs ρ\rho for which outcome aa can occur. Equivalently, when Tr⁡Ia(ρ)≠0\operatorname{Tr}\mathcal I_a(\rho)\ne0,

Tr⁡[ΠaIa(ρ)]Tr⁡Ia(ρ)=1.\frac{ \operatorname{Tr}[\Pi_a\mathcal I_a(\rho)] }{ \operatorname{Tr}\mathcal I_a(\rho) } =1.

The operation may still disturb other observables. Repeatability fixes only the next answer to the same question.

Repeatability is often misread as “the system was not disturbed.” That is too strong.

An ideal measurement of {Πa}\{\Pi_a\} can be repeatable while changing the state from ρ\rho to ρa\rho_a. If the input had coherence between different outcome subspaces, the selective update removes all components outside the observed subspace. If the outcome is ignored, the nonselective channel

ρ⟼∑aΠaρΠa\rho \longmapsto \sum_a\Pi_a\rho\Pi_a

removes coherence between different eigenspaces.

Degeneracy makes the point sharper. A Lüders measurement of a degenerate observable preserves coherence inside the measured eigenspace, while a refined apparatus can secretly measure additional labels and then report only the coarse outcome. Both devices can be repeatable for the coarse eigenvalue, but they leave different states for later compatible measurements. See Lüders Rule.

Destructive detection is different again. A photodetector click may be a sharp event in the detector record while absorbing the photon. It is not a QND measurement of the photon number remaining in the same mode; the carrier has been removed.

Suppose the first measurement prepares the system in the Πa\Pi_a subspace and the system then evolves unitarily under HSH_S for time tt:

U(t)=exp⁡(−iHSt/ℏ).U(t)=\exp(-iH_St/\hbar).

The outcome aa remains certain at the later time if the subspace is invariant under the evolution:

U(t)ΠaU†(t)=Πa.U(t)\Pi_aU^\dagger(t)=\Pi_a.

A sufficient time-independent condition is

[HS,Πa]=0[H_S,\Pi_a]=0

for every measured projector. Then

Tr⁡ ⁣[ΠaU(t)ρaU†(t)]=1.\operatorname{Tr} \!\left[ \Pi_aU(t)\rho_aU^\dagger(t) \right] =1.

If the measured observable is

A=∑aaΠa,A=\sum_a a\Pi_a,

then [HS,A]=0[H_S,A]=0 implies the spectral projectors are conserved, in finite dimension. For unbounded operators, conservation requires the corresponding spectral commutation or invariance statement; a formal commutator on an unspecified domain is insufficient. This is why energy measurements are the textbook example of repeatability over time in a closed system: ideal energy eigenprojectors commute with the Hamiltonian.

For comparison, take A=σzA=\sigma_z and HS=ℏΩσx/2H_S=\hbar\Omega\sigma_x/2. After obtaining +z+z, the state evolves to cos⁡(Ωt/2)∣0⟩−isin⁡(Ωt/2)∣1⟩\cos(\Omega t/2)|0\rangle-i\sin(\Omega t/2)|1\rangle. The next +z+z probability is cos⁡2(Ωt/2)\cos^2(\Omega t/2): immediate repeatability is exact, but finite-time certainty is lost.

The condition can be relaxed. In stroboscopic protocols, one may only require the measured Heisenberg-picture observables at selected times to be mutually compatible:

[A(ti),A(tj)]=0[A(t_i),A(t_j)]=0

for the times actually probed. This allows QND-style monitoring of selected oscillator quadratures at special phases, even when AA is not conserved at every instant.

Position readouts separated by half a cycle

Section titled “Position readouts separated by half a cycle”

Consider an isolated harmonic oscillator with HS=P2/(2m)+mω2X2/2H_S=P^2/(2m)+m\omega^2X^2/2, m,ω>0m,\omega>0, and instantaneous position-measurement pulses. Between pulses,

X(t)=X(0)cos⁡ωt+P(0)mωsin⁡ωt,[X(t),X(t′)]=iℏmωsin⁡[ω(t′−t)].X(t)=X(0)\cos\omega t+\frac{P(0)}{m\omega}\sin\omega t, \qquad [X(t),X(t')]=\frac{i\hbar}{m\omega}\sin[\omega(t'-t)].

At times tn=nπ/ωt_n=n\pi/\omega, the free evolution gives X(tn)=(−1)nX(0)X(t_n)=(-1)^nX(0). The sign reverses each half cycle, but these sampled observables commute. An instantaneous position measurement changes momentum; that momentum disturbance contributes nothing to position at the next exact half cycle. This is the stroboscopic example discussed by Braginsky, Vorontsov, and Thorne (1980), Eq. (4a). It requires the harmonic dynamics and the specified measurement coupling, not only a choice of sampling times.

To quantify the backaction without preparing an unnormalizable position eigenstate, use the Gaussian position instrument with resolution σ>0\sigma>0. Its unread channel is equivalent to an independent zero-mean momentum kick ξ\xi with variance ℏ2/(4σ2)\hbar^2/(4\sigma^2). A pulse at t=0t=0 therefore adds, relative to evolution without that pulse,

δVar⁡X(t)=ℏ24m2ω2σ2sin⁡2ωt.\delta\operatorname{Var}X(t) =\frac{\hbar^2}{4m^2\omega^2\sigma^2}\sin^2\omega t.

This added variance vanishes at every tnt_n, although it is nonzero between them. The same channel adds mean energy ℏ2/(8mσ2)\hbar^2/(8m\sigma^2): nondemolition of the sampled variable does not mean no heating. These are ensemble statements for inputs with finite second moments. Finite-resolution readouts still fluctuate, and the state conditioned on a record still updates; QND does not promise identical raw numbers from a noisy detector. Finite pulse duration, damping, anharmonicity, or a timing error can feed backaction into the next sampled position.

A QND observable is one whose value can be inferred repeatedly without the measurement backaction demolishing the future value being inferred. In a common idealized system-meter model, the interaction has the form

HI(t)=ℏg(t)A⊗PM,H_I(t)=\hbar g(t) A\otimes P_M,

where PMP_M is a meter variable conjugate to the pointer readout. If

[A,HS]=0,[A⊗IM,HI(t)]=0,[A,H_S]=0, \qquad [A\otimes I_M,H_I(t)]=0,

then the interaction can correlate the meter with AA without directly inducing transitions between different eigenspaces of AA. The measurement may dephase superpositions of different AA values, but it does not change the populations in the AA basis in the ideal model.

These commutation relations are best read as sufficient idealizations, not as a complete laboratory definition. Real QND claims must also specify:

  1. the Hilbert space and observable being monitored;
  2. the Hamiltonian terms retained in the approximation;
  3. leakage, loss, relaxation, and detector inefficiency;
  4. the time scale over which nondemolition behavior is asserted;
  5. which later predictions are protected and which variables are allowed to be disturbed.

A measurement can be QND for AA while being strongly demolishing for a conjugate phase or quadrature.

Continuous QND monitoring often appears as measurement-induced dephasing in the monitored basis. For a Hermitian dimensionless observable AA, an unconditional master equation may contain

ρ˙=−iℏ[H,ρ]+ΓmD[A]ρ,\dot\rho = -\frac{i}{\hbar}[H,\rho] + \Gamma_m\mathcal D[A]\rho,

with

D[A]ρ=AρA−12{A2,ρ}.\mathcal D[A]\rho = A\rho A - \frac{1}{2} \{A^2,\rho\}.

If A∣a⟩=a∣a⟩A|a\rangle=a|a\rangle and HH is diagonal in the same basis, then

ρ˙ab=−iℏ(Ea−Eb)ρab−Γm2(a−b)2ρab.\dot\rho_{ab} = -\frac{i}{\hbar}(E_a-E_b)\rho_{ab} - \frac{\Gamma_m}{2}(a-b)^2\rho_{ab}.

For a=ba=b, the dephasing term vanishes. The populations ρaa\rho_{aa} are unchanged by the measurement dephasing, while coherences between different AA values decay. This is the unconditional evolution produced by an ideal nondemolition monitor. The same equation can also describe unobserved environmental or classical noise, so dephasing alone does not demonstrate that a usable measurement record was obtained. In a monitored realization, conditional populations change as the record reveals information, even though the ensemble populations remain fixed.

The conditioned version is described by a stochastic master equation. The measurement record gradually reveals information about AA, while the conditioned state is updated by the observed noise. For the trajectory formalism, see Stochastic Master Equations and Diffusive Trajectories.

Dispersive superconducting-qubit readout is a central approximate QND example. In the simplest two-level dispersive model,

Hdispℏ=(ωr+χσz)a†a+12(ωq+χ)σz.\frac{H_{\mathrm{disp}}}{\hbar} = \left(\omega_r+\chi\sigma_z\right)a^\dagger a + \frac{1}{2}(\omega_q+\chi)\sigma_z.

This Hamiltonian commutes with σz\sigma_z in the idealized qubit subspace:

[Hdisp,σz]=0.[H_{\mathrm{disp}},\sigma_z]=0.

A microwave pulse sent through the resonator acquires a qubit-state-dependent phase. The outgoing field therefore carries information about σz\sigma_z, while the dominant ideal backaction dephases superpositions of ∣g⟩\lvert g\rangle and ∣e⟩\lvert e\rangle.

The word “approximate” matters. Dispersive readout ceases to be perfectly QND when relaxation, Purcell decay, leakage outside the qubit subspace, nonadiabatic pulse effects, dressed-state transitions, residual thermal photons, or detector-induced excitation become relevant. The platform discussion is Circuit QED.

In cavity QED, photon-number QND measurements use a probe whose phase depends on the number of photons in a cavity mode without absorbing those photons in the ideal limit. A dispersive atom-field interaction can have the schematic form

Hintℏ=χa†a σz.\frac{H_{\mathrm{int}}}{\hbar} = \chi a^\dagger a\,\sigma_z.

Since

[a†a,Hint]=0,[a^\dagger a,H_{\mathrm{int}}]=0,

the interaction can imprint photon-number information on the probe atom while preserving the photon number of the cavity field in the idealized interaction. The state may still lose phase coherence between different number states, and real cavities have loss. Thus the QND statement is always limited by photon lifetime, imperfect probes, and unwanted absorption or emission channels.

See Cavity QED for the broader open-system setting.

Gravitational-Wave and Optomechanical Contexts

Section titled “Gravitational-Wave and Optomechanical Contexts”

QND ideas entered gravitational-wave detection because a continuous position measurement of a free mass is not automatically nondemolition. For a free mass,

X(t)=X(0)+P(0)tm,X(t)=X(0)+\frac{P(0)t}{m},

so

[X(t),X(t′)]=iℏm(t′−t).[X(t),X(t')] = \frac{i\hbar}{m}(t'-t).

Backaction on momentum at one time can feed into position at a later time. The standard quantum limit for simple interferometric position readout reflects this tradeoff between imprecision and radiation-pressure backaction.

QND and backaction-evading strategies change the measurement design: they may monitor a variable whose later value is not contaminated by the backaction, correlate imprecision and backaction noise, use squeezed input states, or measure selected mechanical quadratures. These strategies do not remove quantum backaction from nature; they route it away from the signal variable for a specified task. See Optomechanics for the platform-level context.

In the oscillator example, the next pulse arrives at t=π/ω+ϵt=\pi/\omega+\epsilon with ∣ωϵ∣≪1|\omega\epsilon|\ll1. Find the leading added position variance from the first pulse. Compare it with the detector variance σ2\sigma^2 and obtain a condition on ∣ϵ∣|\epsilon| for this contribution to be negligible. Why does making the first readout sharper make the timing requirement stricter?

Solution

The kick contributes δX=ξsin⁡(π+ωϵ)/(mω)\delta X=\xi\sin(\pi+\omega\epsilon)/(m\omega). Averaging its square gives

δVar⁡X=ℏ2ϵ24m2σ2[1+O((ωϵ)2)].\delta\operatorname{Var}X =\frac{\hbar^2\epsilon^2}{4m^2\sigma^2} \left[1+O((\omega\epsilon)^2)\right].

The extra variance is much smaller than σ2\sigma^2 when ∣ϵ∣≪2mσ2/ℏ|\epsilon|\ll2m\sigma^2/\hbar, in addition to ∣ωϵ∣≪1|\omega\epsilon|\ll1. This comparison isolates the first pulse’s backaction; it is not the total prediction error, which also includes the initial state uncertainty and the next detector’s imprecision. Reducing σ\sigma increases the first pulse’s momentum disturbance and reduces the error budget against which it is compared, so the timing scale falls as σ2\sigma^2.

At exact half-cycle spacing this contribution is zero even though the mean energy added by the first pulse remains positive. From P↦P+ξP\mapsto P+\xi, zero mean and independence remove the cross term, giving δ⟨HS⟩=Var⁡(ξ)/(2m)=ℏ2/(8mσ2)\delta\langle H_S\rangle=\operatorname{Var}(\xi)/(2m)=\hbar^2/(8m\sigma^2).

Let U(t)=exp⁡(−iHSt/ℏ)U(t)=\exp(-iH_St/\hbar) and suppose [Πa,HS]=0[\Pi_a,H_S]=0. Show that if ρa\rho_a satisfies ΠaρaΠa=ρa\Pi_a\rho_a\Pi_a=\rho_a, then the later state U(t)ρaU†(t)U(t)\rho_aU^\dagger(t) is still supported in the same subspace.

Solution

The commutator condition implies U(t)Πa=ΠaU(t)U(t)\Pi_a=\Pi_aU(t). Then

ΠaU(t)ρaU†(t)Πa=U(t)ΠaρaΠaU†(t)=U(t)ρaU†(t).\Pi_aU(t)\rho_aU^\dagger(t)\Pi_a = U(t)\Pi_a\rho_a\Pi_aU^\dagger(t) = U(t)\rho_aU^\dagger(t).

Thus a later ideal measurement of the same projector gives outcome aa with probability one.

For a qubit with A=σzA=\sigma_z, consider the unconditional master equation

ρ˙=Γm(σzρσz−ρ).\dot\rho = \Gamma_m \left( \sigma_z\rho\sigma_z-\rho \right).

Write the equations for the density-matrix elements in the σz\sigma_z basis.

Solution

Let

ρ=(ρ00ρ01ρ10ρ11).\rho = \begin{pmatrix} \rho_{00} & \rho_{01}\\ \rho_{10} & \rho_{11} \end{pmatrix}.

Since

σzρσz=(ρ00−ρ01−ρ10ρ11),\sigma_z\rho\sigma_z = \begin{pmatrix} \rho_{00} & -\rho_{01}\\ -\rho_{10} & \rho_{11} \end{pmatrix},

the equation gives

ρ˙00=0,ρ˙11=0,ρ˙01=−2Γmρ01,ρ˙10=−2Γmρ10.\dot\rho_{00}=0, \qquad \dot\rho_{11}=0, \qquad \dot\rho_{01}=-2\Gamma_m\rho_{01}, \qquad \dot\rho_{10}=-2\Gamma_m\rho_{10}.

The populations in the measured basis are unchanged, while coherence between the two measured values decays. This is the unconditional dephasing associated with ideal σz\sigma_z monitoring; an observed record is additional information not specified by this master equation alone.

Why is ordinary absorption photodetection not a QND measurement of photon number in the same mode?

Solution

An absorption click removes a photon from the mode. In quantum-jump language the click operation is proportional to

aρa†.a\rho a^\dagger.

This changes the photon number of the remaining field. The detector record can be sharp and useful, but it demolishes the photon whose presence was detected. A photon-number QND measurement instead tries to infer a†aa^\dagger a through a dispersive probe while preserving the photon number, up to cavity loss and other imperfections.

  • A. Blais, A. L. Grimsmo, S. M. Girvin, and A. Wallraff, “Circuit quantum electrodynamics,” Reviews of Modern Physics 93, 025005 (2021).
  • V. B. Braginsky, Yu. I. Vorontsov, and K. S. Thorne, “Quantum nondemolition measurements,” Science 209, 547–557 (1980), Eq. (4a) and the discussion of stroboscopic measurements. Author-hosted article.
  • V. B. Braginsky and F. Y. Khalili, Quantum Measurement, Cambridge University Press (1992).
  • C. M. Caves, K. S. Thorne, R. W. P. Drever, V. D. Sandberg, and M. Zimmermann, “On the measurement of a weak classical force coupled to a quantum-mechanical oscillator. I. Issues of principle,” Reviews of Modern Physics 52, 341–392 (1980).
  • S. Gleyzes et al., “Quantum jumps of light recording the birth and death of a photon in a cavity,” Nature 446, 297–300 (2007).
  • S. Haroche and J.-M. Raimond, Exploring the Quantum: Atoms, Cavities, and Photons, Oxford University Press (2006).
  • K. Jacobs, Quantum Measurement Theory and its Applications, Cambridge University Press (2014).
  • H. M. Wiseman and G. J. Milburn, Quantum Measurement and Control, Cambridge University Press (2010).