Skip to content

Quantum illumination is entanglement-assisted discrimination between two lossy, noisy channels: one describing an absent target and one describing a weakly reflecting target. A source prepares correlated signal–idler mode pairs. The signal interrogates a target region, the idler is retained, and a receiver jointly processes the collected return and retained idler to decide which channel acted.

Its best-known result is striking but narrow. For many independent low-brightness mode pairs, a weak specular reflector of known phase, and a bright thermal background, an ideal two-mode-squeezed transmitter with an optimum joint receiver has an error-probability exponent four times that of an optimum-reception coherent-state transmitter at the same transmitted signal energy. This factor of four is conventionally called a 6 dB error-exponent advantage. It can survive even when loss and noise make the received return–idler state separable.

That theorem does not establish a six-decibel improvement in received power, range, cross section, or arbitrary radar performance. It assumes a declared hypothesis pair, known mode matching, a particular background model, an ideal idler, and an optimum receiver. Unknown phase, speckle, range search, Doppler, clutter, receiver loss, idler storage, and wall-clock duty cycle can change or remove the advantage. Quantum illumination is therefore best understood as a precise result in quantum channel discrimination, not as a synonym for every sensor using correlated radiation.

Required background. Quantum Measurement as Estimation introduces priors and decision losses; Density Operators for Quantum Information supplies mixed states and trace distance; Two-Mode Entanglement develops the signal–idler state; POVMs defines general measurement outcomes.

Target absence and presence as quantum channels

Section titled “Target absence and presence as quantum channels”

For one matched temporal-frequency-spatial mode, let a^S\hat a_S denote the transmitted signal, a^I\hat a_I the retained idler, and a^R\hat a_R the collected return. A standard specular-target model is

H0:a^R=b^0,H1:a^R=κ eiϕa^S+1−κ b^1.\begin{aligned} H_0 &: \quad \hat a_R = \hat b_0, \\[3pt] H_1 &: \quad \hat a_R = \sqrt{\kappa}\,e^{i\phi}\hat a_S + \sqrt{1-\kappa}\,\hat b_1. \end{aligned}

Here H0H_0 means target absent, H1H_1 means target present, κ≪1\kappa\ll1 is the round-trip transmissivity into the collected mode, and ϕ\phi is the return phase. The bath modes are thermal and independent of the source. In the usual no-passive-signature convention, they are chosen so that the received background brightness is NBN_B under either hypothesis:

⟨b^0†b^0⟩=NB,⟨b^1†b^1⟩=NB1−κ.\langle \hat b_0^\dagger\hat b_0\rangle=N_B, \qquad \langle \hat b_1^\dagger\hat b_1\rangle = \frac{N_B}{1-\kappa}.

Consequently,

⟨a^R†a^R⟩H0=NB,⟨a^R†a^R⟩H1=NB+κNS,\begin{aligned} \langle \hat a_R^\dagger\hat a_R\rangle_{H_0} &=N_B, \\ \langle \hat a_R^\dagger\hat a_R\rangle_{H_1} &=N_B+\kappa N_S, \end{aligned}

where NS=⟨a^S†a^S⟩N_S=\langle\hat a_S^\dagger\hat a_S\rangle is signal photons per mode. This convention prevents an artificial detector from deciding between H0H_0 and H1H_1 merely because inserting the target changes the background level. A different physical background model is legitimate, but it defines a different discrimination problem and must be used for both quantum and classical systems.

Let Π0\Pi_0 and Π1\Pi_1 be the receiver’s two decision operators, with Π0+Π1=I\Pi_0+\Pi_1=\mathbb I. The false-alarm and miss probabilities are

PFA=Tr⁡(Π1ρ0),PM=Tr⁡(Π0ρ1).P_{\mathrm{FA}} = \operatorname{Tr}(\Pi_1\rho_0), \qquad P_{\mathrm{M}} = \operatorname{Tr}(\Pi_0\rho_1).

For prior probabilities π0\pi_0 and π1\pi_1 and equal costs for the two error types, the Bayes error is

Pe=π0PFA+π1PM.P_e = \pi_0P_{\mathrm{FA}} + \pi_1P_{\mathrm{M}}.

Equal priors are mathematically convenient, but many sensing systems operate at low target prevalence and impose a stringent false-alarm constraint. The appropriate output is then a receiver-operating characteristic or a Neyman–Pearson miss probability at fixed PFAP_{\mathrm{FA}}, not only an equal-prior error number. Threshold tuning after seeing test data invalidates either comparison.

A broadband pulse or observation interval contains many approximately orthogonal mode pairs. In the ideal independent-and-identically-distributed model,

ρj(M)=ρj⊗M,j∈{0,1}.\rho_j^{(M)} = \rho_j^{\otimes M}, \qquad j\in\{0,1\}.

For duration TT and usable bandwidth WW, the order of magnitude is M∼TWM\sim TW, up to real-mode versus complex-mode conventions and polarization factors. The transmitted photon budget is

Ntx=MNS,Etx=ℏωMNSN_{\mathrm{tx}}=MN_S, \qquad E_{\mathrm{tx}}=\hbar\omega MN_S

for a narrowband carrier. Reporting NSN_S without MM, or MM without source brightness, does not specify the exposure.

The canonical Gaussian transmitter uses MM two-mode squeezed-vacuum pairs. For each pair,

∣ψ⟩SI=∑n=0∞NSn(NS+1)n+1∣n⟩S∣n⟩I.\lvert\psi\rangle_{SI} = \sum_{n=0}^{\infty} \sqrt{\frac{N_S^n}{(N_S+1)^{n+1}}} \lvert n\rangle_S\lvert n\rangle_I.

Each marginal is thermal,

⟨a^S†a^S⟩=⟨a^I†a^I⟩=NS,\langle \hat a_S^\dagger\hat a_S\rangle = \langle \hat a_I^\dagger\hat a_I\rangle =N_S,

but the pair has phase-sensitive cross correlation

⟨a^Sa^I⟩=NS(NS+1)\langle \hat a_S\hat a_I\rangle = \sqrt{N_S(N_S+1)}

in a convenient phase convention. A classical signal–idler state with a positive Glauber–Sudarshan representation obeys the Cauchy–Schwarz bound

∣⟨a^Sa^I⟩∣2≤NSNI.\bigl\lvert \langle \hat a_S\hat a_I\rangle \bigr\rvert^2 \leq N_SN_I.

For equal signal and idler brightness NI=NSN_I=N_S, a two-mode squeezed vacuum has correlation squared NS(NS+1)N_S(N_S+1) rather than at most NS2N_S^2. The relative gap is especially large when NS≪1N_S\ll1. Quantum illumination encodes target presence in the weak remnant of this initially nonclassical correlation.

Use X=(a+a†)/2X=(a+a^\dagger)/\sqrt2 and P=(a−a†)/(i2)P=(a-a^\dagger)/(i\sqrt2), so [X,P]=i[X,P]=i and vacuum variance is 1/21/2. Order the covariance as (XR,PR,XI,PI)(X_R,P_R,X_I,P_I) and define

Z=(100−1).\mathbf Z = \begin{pmatrix} 1&0\\ 0&-1 \end{pmatrix}.

For known phase ϕ=0\phi=0, the return–idler covariance matrices are

VRI(0)=((NB+12)I200(NS+12)I2)V_{RI}^{(0)} = \begin{pmatrix} (N_B+\tfrac12)\mathbb I_2 & 0\\ 0 & (N_S+\tfrac12)\mathbb I_2 \end{pmatrix}

and

VRI(1)=((NB+κNS+12)I2κNS(NS+1) ZκNS(NS+1) Z(NS+12)I2).V_{RI}^{(1)} = \begin{pmatrix} (N_B+\kappa N_S+\tfrac12)\mathbb I_2 & \sqrt{\kappa N_S(N_S+1)}\,\mathbf Z \\ \sqrt{\kappa N_S(N_S+1)}\,\mathbf Z & (N_S+\tfrac12)\mathbb I_2 \end{pmatrix}.

When NB≫κNSN_B\gg\kappa N_S, the return marginal changes only slightly. The main signature available to a joint receiver is the off-diagonal block. Under H0H_0 it vanishes; under H1H_1 it is proportional to κNS(NS+1)\sqrt{\kappa N_S(N_S+1)}.

In this convention, the thermal-loss channel is entanglement breaking when NB≥κN_B\geq\kappa. Its added quadrature noise is Y=[NB+(1−κ)/2]I2Y=[N_B+(1-\kappa)/2]I_2; the phase-insensitive attenuator is entanglement breaking when Y≥(1+κ)I2/2Y\geq(1+\kappa)I_2/2, which gives this threshold. The returned mode and stored idler are then separable even under H1H_1. This does not make the two conditional states identical. The initial entangled source created a phase-sensitive correlation stronger than any equal-energy classical source could create; attenuation leaves a weak but measurable remnant of that correlation.

The operational resource statement is therefore:

Initial signal–idler entanglement enables a stronger transmitter correlation, while target detection uses distinguishability between the two final return–idler states.

It is misleading to say that entanglement itself travelled out and back intact. It is also unnecessary to identify one surviving correlation measure, such as discord, as the unique cause of the advantage. The performance claim is settled by a matched channel-discrimination bound.

For arbitrary priors, the minimum Bayes error over all two-outcome POVMs is

Pe⋆=12[1−∥π1ρ1(M)−π0ρ0(M)∥1].P_e^\star = \frac12 \left[ 1- \left\lVert \pi_1\rho_1^{(M)} - \pi_0\rho_0^{(M)} \right\rVert_1 \right].

For equal priors this becomes

Pe⋆=12[1−12∥ρ1(M)−ρ0(M)∥1].P_e^\star = \frac12 \left[ 1- \frac12 \left\lVert \rho_1^{(M)}-\rho_0^{(M)} \right\rVert_1 \right].

The corresponding Helstrom measurement projects onto the positive and negative eigenspaces of π1ρ1(M)−π0ρ0(M)\pi_1\rho_1^{(M)}-\pi_0\rho_0^{(M)}. This defines the optimum, but directly diagonalizing a many-mode mixed-state difference does not by itself give a practical optical receiver.

For equal priors and independent copies, define

Q:=inf⁡0≤s≤1Tr⁡(ρ0sρ11−s),ξQ:=−ln⁡Q.Q := \inf_{0\leq s\leq1} \operatorname{Tr} \left( \rho_0^s\rho_1^{1-s} \right), \qquad \xi_Q:=-\ln Q.

The quantum Chernoff bound is

Pe⋆(M)≤12QM=12e−MξQ.P_e^\star(M) \leq \frac12Q^M = \frac12e^{-M\xi_Q}.

More importantly, the exponent is asymptotically exact:

lim⁡M→∞−1Mln⁡Pe⋆(M)=ξQ.\lim_{M\to\infty} -\frac1M\ln P_e^\star(M) = \xi_Q.

An exponent is not an exact finite-sample error probability. Prefactors, non-identical modes, calibration uncertainty, and threshold estimation matter before the asymptotic regime is reached.

The familiar comparison assumes

NS≪1,κ≪1,NB≫1,M≫1,N_S\ll1, \qquad \kappa\ll1, \qquad N_B\gg1, \qquad M\gg1,

with known return phase and ideal mode matching. Let

γ:=κNSNB.\gamma := \frac{\kappa N_S}{N_B}.

For a coherent-state transmitter with optimum reception, the leading error exponent is

ξcoh≃γ4,Pe,coh=exp⁡[−Mξcoh+o(M)].\xi_{\mathrm{coh}} \simeq \frac{\gamma}{4}, \qquad P_{e,\mathrm{coh}} =\exp[-M\xi_{\mathrm{coh}}+o(M)].

This is not merely a comparison with direct detection. A coherent state is the reference classical transmitter for the specified Gaussian channel, and it is paired here with optimum quantum reception. The original positive-PP analysis also bounds arbitrary classical signal–idler transmitters, but that general bound is looser than the exact coherent-state result. An experiment claiming the best classical benchmark must therefore state its allowed transmitter, idler, receiver, amplifier, and reference resources rather than infer that claim from a comparison with one convenient classical apparatus.

Structured pairwise quantum-illumination receivers such as ideal OPA and phase-conjugate receivers attain

ξpair≃γ2.\xi_{\mathrm{pair}} \simeq \frac{\gamma}{2}.

The optimum joint quantum receiver attains

ξQI≃γ.\xi_{\mathrm{QI}} \simeq \gamma.

Thus

ξpairξcoh≃2,ξQIξcoh≃4.\frac{\xi_{\mathrm{pair}}}{\xi_{\mathrm{coh}}} \simeq2, \qquad \frac{\xi_{\mathrm{QI}}}{\xi_{\mathrm{coh}}} \simeq4.

Expressed in decibels,

10log⁡102≃3.01 dB,10log⁡104≃6.02 dB.10\log_{10}2\simeq3.01\ \mathrm{dB}, \qquad 10\log_{10}4\simeq6.02\ \mathrm{dB}.

The logarithm applies to a ratio of error exponents. It does not mean that PeP_e itself is smaller by 6 dB. Because PeP_e is exponential in MM, a fourfold exponent can yield a much larger error-probability ratio after many modes. Conversely, outside the asymptotic assumptions there need not be a fourfold ratio at all.

For the ideal known-phase Gaussian model, the coherent-state Chernoff exponent before the bright-background approximation is

ξcoh=κNS(NB+1−NB)2.\xi_{\mathrm{coh}} = \kappa N_S \left( \sqrt{N_B+1}-\sqrt{N_B} \right)^2.

Since

(NB+1−NB)2≃14NB\left( \sqrt{N_B+1}-\sqrt{N_B} \right)^2 \simeq \frac{1}{4N_B}

for NB≫1N_B\gg1, this reduces to the expression above. An accessible comparison is phase-matched homodyne detection. Each return quadrature is Gaussian with variance NB+1/2N_B+1/2 and means separated by 2κNS\sqrt{2\kappa N_S}. A midpoint threshold on the sum of MM measurements gives the exact error

Pe,hom(M)=12erfc⁡ ⁣(MκNS4NB+2).P_{e,\rm hom}(M) =\frac12\operatorname{erfc}\!\left( \sqrt{\frac{M\kappa N_S}{4N_B+2}}\right).

Its exponent κNS/(4NB+2)\kappa N_S/(4N_B+2) approaches the optimal coherent exponent in bright background, but homodyne is not the optimal receiver at every NBN_B. The complementary error function also makes the nonconstant finite-MM prefactor explicit.

Quoting the 6 dB result at optical frequencies with negligible thermal occupancy, without recomputing the relevant bound, silently changes regime.

The Standard Quantum Limit is a matched independent-probe benchmark for an estimation problem, commonly expressed as a variance or root-mean-square error. Quantum illumination is usually a finite-separation hypothesis test, and its central figure of merit is an error probability, ROC, or asymptotic error exponent. Calling ξcoh\xi_{\mathrm{coh}} an SQL obscures the loss function and invites an invalid comparison between 1/N1/\sqrt N precision scaling and exponential decision error. A reflectivity-estimation version can be formulated with Fisher information, but it is a different task from target absence versus presence.

Compare equal transmitted energy, mode bandwidth and duration, target and background models, priors or false-alarm constraints, and allowed receiver resources. A classical reference can use optimum quantum reception even though its transmitter has a positive PP representation. Beating a split-thermal source or coherent direct detection does not automatically beat that stronger reference. Idler storage, references, loss, and preparation failures must be included when comparing complete instruments.

The source creates an opportunity; the receiver decides how much of it is realized.

A low-gain optical parametric amplifier mixes each return–idler pair. In one phase convention an output mode is

c^m=G a^Im+G−1 a^Rm†,G−1≪1.\hat c_m = \sqrt{G}\,\hat a_{I_m} + \sqrt{G-1}\,\hat a_{R_m}^{\dagger}, \qquad G-1\ll1.

Writing CRI=⟨aRaI⟩C_{RI}=\langle a_Ra_I\rangle, the target-present mean is

⟨cm†cm⟩1=GNS+(G−1)(NB+κNS+1)+2G(G−1)Re⁡CRI.\langle c_m^\dagger c_m\rangle_1 =G N_S+(G-1)(N_B+\kappa N_S+1) +2\sqrt{G(G-1)}\operatorname{Re}C_{RI}.

Under H0H_0, both CRIC_{RI} and the κNS\kappa N_S term vanish. The last term is the useful phase-sensitive signal, but a weaker return-intensity term also remains. Photon counts accumulated over MM modes therefore shift between the hypotheses. With optimized gain in the canonical regime, this receiver reaches the 3 dB exponent advantage. It requires stable phase and good return–idler mode matching, and it remains suboptimal even with ideal components.

A phase-conjugate receiver converts the return so that its phase-sensitive correlation with the idler becomes a phase-insensitive interference signal. The converted return and idler interfere on a balanced beam splitter, and the difference of output counts is the test statistic. Its ideal asymptotic exponent also reaches the 3 dB level.

Physical phase conjugation adds vacuum noise. A digital implementation that heterodynes the return and idler and conjugates one classical record is not the same measurement as lossless optical joint detection: separate heterodyne measurements have already added measurement noise and destroyed access to a later collective POVM.

Sum-frequency generation is the inverse nonlinear process to down-conversion. For many weak return–idler pairs with common phase-sensitive correlation, the interaction

H^SFG=ℏg∑m=1M(b^†a^Rma^Im+b^a^Rm†a^Im†)\hat H_{\mathrm{SFG}} = \hbar g \sum_{m=1}^{M} \left( \hat b^\dagger \hat a_{R_m}\hat a_{I_m} + \hat b \hat a_{R_m}^\dagger\hat a_{I_m}^\dagger \right)

can convert pair correlations into a sum-frequency mode b^\hat b, for real gg in a fixed pump-phase convention. Initially, d⟨b⟩/dt=−ig∑m⟨aRmaIm⟩d\langle b\rangle/dt=-ig\sum_m\langle a_{R_m}a_{I_m}\rangle. For MM identical aligned pairs this is −igMCRI-igM C_{RI}, not a M\sqrt M amplitude at fixed gg and interaction time. A normalized sum M−1/2∑maRmaImM^{-1/2}\sum_m a_{R_m}a_{I_m} instead has mean MCRI\sqrt M C_{RI}. Conversion time, depletion, and noise must be included to obtain receiver performance; initial coherent growth is not an M2M^2 information law. A staged SFG receiver can asymptotically attain the quantum Chernoff exponent in the ideal low-brightness model (Zhuang, Zhang, and Shapiro). Feedforward nulling between stages can approach the Helstrom limit in the low-brightness regime.

The ideal construction requires controlled weak nonlinear interactions, mode-resolved processing, idler storage, and adaptive photon counting; losses in those operations change the predicted error.

Correlation-to-displacement receivers conditionally map broadband return–idler correlation onto a coherent displacement of a selected mode. The final problem can then be attacked with coherent-state receivers such as Kennedy-like nulling or photon counting. This architecture is theoretically capable of optimal performance for several known-phase entanglement-assisted tasks and may reduce the nonlinear and memory burden of staged SFG.

Its practical performance still depends on conversion efficiency, amplifier noise, idler loss, phase knowledge, mode selection, and the final detector. Ideal noiseless amplification and actual microwave amplification are not interchangeable assumptions.

Time-tagged photon-pair experiments often estimate a coincidence excess or a photon-number covariance,

CRI=⟨n^Rn^I⟩−⟨n^R⟩⟨n^I⟩.C_{RI} = \langle \hat n_R\hat n_I\rangle - \langle \hat n_R\rangle \langle \hat n_I\rangle.

Such receivers can reject uncorrelated background and give accessible proof-of-principle demonstrations. Their performance depends on coincidence window, detector jitter, dark counts, multipair emission, accidental coincidences, and whether no-click records enter the likelihood. Discarding noncoincident trials without counting their acquisition cost can make a conditional statistic look better while reducing information per transmitted photon or per second.

Loss, Imperfections, and Nuisance Parameters

Section titled “Loss, Imperfections, and Nuisance Parameters”

Signal attenuation is already represented by κ\kappa and is the regime in which quantum illumination was designed to operate. Idler loss is different: it directly erases the local reference needed to read the return correlation. With return efficiency ηR\eta_R, idler-memory efficiency ηI\eta_I, mode overlap amplitude μ\mu, and Gaussian phase jitter of variance σϕ2\sigma_\phi^2, a useful schematic correlation budget is

CRIeff=μηRηIκ e−σϕ2/2NS(NS+1).C_{RI}^{\mathrm{eff}} = \mu \sqrt{\eta_R\eta_I\kappa} \,e^{-\sigma_\phi^2/2} \sqrt{N_S(N_S+1)}.

This is an averaged correlation, not a complete noisy-state specification. The expression assumes independent, zero-mean phase jitter and uncorrelated loss environments. Observed phase tags can allow compensation; a shared unknown phase across many pairs produces correlated records. Noise variances and the full likelihood are still needed before converting ∣CRIeff∣2\lvert C_{RI}^{\mathrm{eff}}\rvert^2 into a signal-to-noise or error ratio. A long idler delay can therefore erase the ideal factor-of-four advantage even though the outward signal path is supposed to be lossy. Any proposal for long range must specify how the idler is stored for the round-trip time without excessive attenuation, added noise, or mode drift.

At NS≪1N_S\ll1, two-mode squeezing offers its largest correlation advantage per transmitted photon, but useful decisions require many modes. Increasing pump power raises NSN_S and total flux but moves the source away from the regime in which the canonical exponent ratios are derived. In photon-pair language it also raises multipair events. Source purity, unwanted thermal population, spectral Schmidt number, pump phase noise, and collection loss must be included in the state delivered to the transmit and idler reference planes.

The 6 dB Gaussian result assumes a specular return with known phase. A rough target may produce speckle with random amplitude and nearly uniform phase. For one pair, the target-present state is then a mixture,

ρˉ1=∫dκ dϕ p(κ,ϕ)ρ1(κ,ϕ).\bar\rho_1 = \int d\kappa\,d\phi\, p(\kappa,\phi)\rho_1(\kappa,\phi).

If the same unknown κ,ϕ\kappa,\phi persist throughout a pulse, however, the MM-pair state is

ρˉ1(M)=∫dκ dϕ p(κ,ϕ)ρ1(κ,ϕ)⊗M,\bar\rho_1^{(M)} =\int d\kappa\,d\phi\,p(\kappa,\phi) \rho_1(\kappa,\phi)^{\otimes M},

which generally differs from ρˉ1⊗M\bar\rho_1^{\otimes M}. The latter describes fresh independent fading in every mode. Averaging and taking tensor powers do not commute: common phase can leave relative information across the block even when each marginal loses its mean correlation. The cited Rayleigh-fading analysis uses a common fading coefficient over the return bandwidth and observation interval.

Uniform phase averaging gives

∫02πdϕ2πeiϕ=0,\int_0^{2\pi}\frac{d\phi}{2\pi} e^{i\phi} =0,

so the mean phase-sensitive correlation used by an OPA receiver vanishes, although its small return-intensity contribution need not vanish. The OPA’s known-phase 3 dB advantage can disappear. More collective receivers can retain an advantage for some fading models, but its scaling may become subexponential rather than a simple 6 dB exponent improvement. Phase tracking with a reference beam consumes power and changes the comparison.

The basic protocol asks whether a target occupies one known mode. Radar and lidar usually search over delay, angle, and Doppler cells. If the return could occupy one of KK modes, the receiver faces a multiple-hypothesis problem,

H0,H1,1,H1,2,…,H1,K,H_0, H_{1,1}, H_{1,2}, \ldots, H_{1,K},

with a search penalty and a mode bank. Waveform ambiguity, sidelobes, multiple comparisons, moving-target phase, and the cost of a quantum memory must be included. A one-cell target-detection exponent is not a complete ranging or tracking analysis.

For a bosonic mode at temperature TT,

Nth(ω,T)=1exp⁡(ℏω/kBT)−1.N_{\mathrm{th}}(\omega,T) = \frac{1}{ \exp(\hbar\omega/k_{\mathrm B}T)-1 }.

Room-temperature microwave modes can have Nth≫1N_{\mathrm{th}}\gg1, while optical thermal occupancy is ordinarily negligible. Optical experiments often inject background light to access the bright-noise regime. Microwave systems must contend with cryogenic-to-room-temperature interfaces and amplifier noise before and after the target path.

Thermal Gaussian noise is not a universal model for terrain clutter, multipath, atmospheric fading, interference, or an adaptive jammer. Robustness to an entanglement-breaking thermal channel does not imply immunity to spoofing, interception, or adversarial waveform knowledge.

Consider an idealized known-phase test with

NS=10−2,κ=10−2,NB=20,M=2×106.N_S=10^{-2}, \qquad \kappa=10^{-2}, \qquad N_B=20, \qquad M=2\times10^6.

The dimensionless accumulated scale is

Mγ=MκNSNB=10.M\gamma = \frac{M\kappa N_S}{N_B} =10.

Using only the leading asymptotic exponents gives

systemMξ12e−Mξcoherent, optimum2.54.10×10−2OPA or phase-conjugate53.37×10−3QI, optimum joint102.27×10−5\begin{array}{c|c|c} \text{system} & M\xi & \tfrac12e^{-M\xi} \\ \hline \text{coherent, optimum} &2.5 &4.10\times10^{-2} \\ \text{OPA or phase-conjugate} &5 &3.37\times10^{-3} \\ \text{QI, optimum joint} &10 &2.27\times10^{-5} \end{array}

These are exponent-level estimates, not certified finite-MM errors. They show why saying “6 dB lower error probability” is wrong: the two approximate errors differ by orders of magnitude here, not by a fixed factor of four.

For a leading-envelope target of 10−310^{-3}, define

L:=ln⁡ ⁣(12Pe)=ln⁡500≃6.215.L := \ln\!\left(\frac{1}{2P_e}\right) = \ln 500 \simeq6.215.

The corresponding leading mode requirements are

Mcoh≃4NBLκNS≃4.97×106,Mpair≃2NBLκNS≃2.49×106,MQI≃NBLκNS≃1.24×106.\begin{aligned} M_{\mathrm{coh}} &\simeq \frac{4N_BL}{\kappa N_S} \simeq4.97\times10^6, \\ M_{\mathrm{pair}} &\simeq \frac{2N_BL}{\kappa N_S} \simeq2.49\times10^6, \\ M_{\mathrm{QI}} &\simeq \frac{N_BL}{\kappa N_S} \simeq1.24\times10^6. \end{aligned}

The ideal joint receiver uses one quarter as many modes as the coherent benchmark at the same exponent target. Whether it finishes sooner depends on source bandwidth, idler storage, receiver cycle time, calibration, and accepted-shot fraction.

Three asymptotic error envelopes reach the same error at different accumulated signal-to-background exposures

Leading envelopes 12e−rx\tfrac12e^{-rx} with x=MκNS/NBx=M\kappa N_S/N_B and r=1/4,1/2,1r=1/4,1/2,1 for coherent, pairwise, and optimum-QI reception. At envelope 10−310^{-3}, the required xx values are approximately 24.924.9, 12.412.4, and 6.26.2. These illustrate the asymptotic mode-count ratio; they are not exact finite-mode error probabilities. The assumptions are low brightness, weak reflection, bright background, known phase, and ideal reception.

Suppose each fixed-duration cycle TcT_c delivers McM_c useful pairs with probability pp, the same under both hypotheses. The failure flag carries no target information, and successful cycles are independent. If a useful pair has Chernoff coefficient Q=e−ξQ=e^{-\xi}, the full cycle has

Qcycle=(1−p)+pQMc,Rξ=−ln⁡[1−p+pe−Mcξ]Tc.Q_{\rm cycle}=(1-p)+pQ^{M_c}, \qquad \mathcal R_\xi= \frac{-\ln[1-p+p e^{-M_c\xi}]}{T_c}.

The direct-sum success and failure records give the first expression; the same Chernoff minimizer applies to the identical accepted modes. It applies to quantum states with optimum reception, or to classical outcome distributions with a specified receiver. For Mcξ≪1M_c\xi\ll1, Rξ≃pMcξ/Tc\mathcal R_\xi\simeq pM_c\xi/T_c. Simply multiplying by the mean accepted fraction is not exact for a fixed number of attempted cycles: unusually many failures contribute to the average error. As Mcξ→∞M_c\xi\to\infty, the exponent per attempted cycle approaches −ln⁡(1−p)-\ln(1-p).

Stopping after a fixed number of successful cycles is a different design. Its expected duration is increased by 1/p1/p, and its exponent per expected duration is pMcξ/TcpM_c\xi/T_c; the completion time is then random. If acceptance depends on the target, include its informative flag in the joint model.

When comparing instruments, include photons transmitted on failed cycles, bandwidth, idler storage, receiver efficiency, and acquisition time. Pump, cryogenic, and local-oscillator power may matter for deployment even when a theoretical comparison fixes only photons incident on the sample.

Spontaneous parametric down-conversion naturally supplies broadband signal–idler pairs, low thermal source noise, mature photon counting, and optical nonlinear processing. Optical backgrounds are not normally bright thermal baths per mode, so experiments often inject broadband noise. Free-space optical targets introduce turbulence, rough-surface speckle, detector dark counts, daylight background, and mode-matching challenges. At long range, an optical idler can be delayed in fiber, but fiber attenuation and dispersion grow with round-trip time.

Microwave frequencies naturally provide bright room-temperature thermal backgrounds and are relevant to conventional radar bands. Josephson parametric devices can generate two-mode-squeezed microwave fields, but they operate cryogenically. Sending one mode into a room-temperature path while retaining and later measuring the idler requires low-loss coupling, amplification, isolation, and often a quantum memory or transduction stage.

Conventional microwave receivers are exceptionally mature. Their coherent transmit power, antenna aperture, low-noise amplification, waveform agility, and digital processing form a demanding baseline. A laboratory quantum advantage in a low-brightness channel is scientifically meaningful without yet establishing superior range, scan rate, robustness, or cost over a field radar.

Optical experiments with OPA reception (Zhang and colleagues) and microwave experiments with idler storage and joint measurement (Assouly and colleagues) provide examples of the detection mechanism under controlled channels. They test particular transmitters, receivers, and noise models. Their performance cannot be converted into a range or search-rate advantage without analyzing those additional tasks. The broader experimental questions are discussed in Quantum Metrology Frontiers.

For a classical signal–idler state represented as a probability mixture of coherent states with complex amplitudes αS\alpha_S and αI\alpha_I, prove

∣⟨a^Sa^I⟩∣2≤NSNI.\bigl\lvert\langle\hat a_S\hat a_I\rangle\bigr\rvert^2 \leq N_SN_I.

Compare the result with a two-mode squeezed vacuum for NI=NSN_I=N_S.

Solution

For a positive PP representation, normally ordered moments are classical averages. Therefore

⟨a^Sa^I⟩=E(αSαI),NS=E∣αS∣2,NI=E∣αI∣2.\langle\hat a_S\hat a_I\rangle = \mathbb E(\alpha_S\alpha_I), \qquad N_S=\mathbb E\lvert\alpha_S\rvert^2, \qquad N_I=\mathbb E\lvert\alpha_I\rvert^2.

Cauchy–Schwarz gives

∣E(αSαI)∣2≤E∣αS∣2E∣αI∣2=NSNI.\bigl\lvert\mathbb E(\alpha_S\alpha_I)\bigr\rvert^2 \leq \mathbb E\lvert\alpha_S\rvert^2 \mathbb E\lvert\alpha_I\rvert^2 =N_SN_I.

For a two-mode squeezed vacuum with equal marginal brightness,

∣⟨a^Sa^I⟩∣2=NS(NS+1)>NS2\bigl\lvert\langle\hat a_S\hat a_I\rangle\bigr\rvert^2 =N_S(N_S+1)>N_S^2

for every NS>0N_S>0. The fractional gap is largest at low brightness.

Show that choosing ⟨b^1†b^1⟩=NB/(1−κ)\langle\hat b_1^\dagger\hat b_1\rangle=N_B/(1-\kappa) in the target-present channel gives return brightness NB+κNSN_B+\kappa N_S.

Solution

Independence and zero bath mean remove cross terms. Hence

⟨a^R†a^R⟩H1=κ⟨a^S†a^S⟩+(1−κ)⟨b^1†b^1⟩=κNS+(1−κ)NB1−κ=NB+κNS.\begin{aligned} \langle\hat a_R^\dagger\hat a_R\rangle_{H_1} &= \kappa \langle\hat a_S^\dagger\hat a_S\rangle + (1-\kappa) \langle\hat b_1^\dagger\hat b_1\rangle \\ &= \kappa N_S + (1-\kappa)\frac{N_B}{1-\kappa} \\ &=N_B+\kappa N_S. \end{aligned}

Thus the target does not create a change in the bath contribution itself; it adds only the attenuated signal energy.

Suppose two systems have exponents ξA\xi_A and ξB=4ξA\xi_B=4\xi_A. Compute the advantage in decibels and the asymptotic ratio of mode counts needed to reach the same target error.

Solution

The exponent advantage is

10log⁡10ξBξA=10log⁡104≃6.02 dB.10\log_{10}\frac{\xi_B}{\xi_A} = 10\log_{10}4 \simeq6.02\ \mathrm{dB}.

At leading order Pe∼12e−MξP_e\sim\tfrac12e^{-M\xi}, so fixed PeP_e requires MξM\xi to be fixed. Therefore

MBMA=ξAξB=14.\frac{M_B}{M_A} = \frac{\xi_A}{\xi_B} = \frac14.

This mode-count result, not a fixed fourfold probability ratio, is the direct meaning of the exponent comparison.

Under H1H_1, suppose the return phase is uniform on [0,2π)[0,2\pi). Show that the ensemble-averaged phase-sensitive return–idler correlation vanishes. Why does this invalidate the known-phase OPA argument but not prove that every quantum receiver is useless?

Solution

Conditioned on ϕ\phi,

⟨a^Ra^I⟩ϕ=κ eiϕNS(NS+1).\langle\hat a_R\hat a_I\rangle_{\phi} = \sqrt{\kappa}\,e^{i\phi} \sqrt{N_S(N_S+1)}.

Averaging gives

Eϕ⟨a^Ra^I⟩ϕ=κNS(NS+1)∫02πdϕ2πeiϕ=0.\mathbb E_\phi \langle\hat a_R\hat a_I\rangle_\phi = \sqrt{\kappa N_S(N_S+1)} \int_0^{2\pi}\frac{d\phi}{2\pi}e^{i\phi} =0.

The OPA’s leading mean-count shift is linear in a chosen quadrature of this correlation, so averaging removes that contribution. The weaker intensity shift (G−1)κNS(G-1)\kappa N_S remains; zero average correlation does not make the OPA output distributions identical. The mixed state ρˉ1\bar\rho_1 can still differ from ρ0\rho_0 through higher-order or phase-insensitive structure. A receiver designed for the mixture may retain some advantage, but the known-phase exponent cannot simply be reused.

Take ηR=0.8\eta_R=0.8, ηI=0.5\eta_I=0.5, κ=10−3\kappa=10^{-3}, NS=0.02N_S=0.02, ∣μ∣=0.9\lvert\mu\rvert=0.9, and σϕ=0.3\sigma_\phi=0.3. Estimate ∣CRIeff∣\lvert C_{RI}^{\mathrm{eff}}\rvert.

Solution

Using the schematic budget,

∣CRIeff∣=0.9(0.8)(0.5)(10−3)e−0.32/2(0.02)(1.02)≃2.46×10−3.\begin{aligned} \lvert C_{RI}^{\mathrm{eff}}\rvert &= 0.9 \sqrt{(0.8)(0.5)(10^{-3})} e^{-0.3^2/2} \sqrt{(0.02)(1.02)} \\ &\simeq 2.46\times10^{-3}. \end{aligned}

The source correlation before these losses is 0.02×1.02≃0.143\sqrt{0.02\times1.02}\simeq0.143. The small returned correlation illustrates why a receiver must coherently accumulate many matched modes.

Estimate the room-temperature thermal occupancy at 10 GHz10\ \mathrm{GHz} and at an optical frequency of 200 THz200\ \mathrm{THz}. Use T=300 KT=300\ \mathrm K and kBT/h≃6.25 THzk_{\mathrm B}T/h\simeq6.25\ \mathrm{THz}.

Solution

Since ℏω=hν\hbar\omega=h\nu,

Nth=1ehν/kBT−1.N_{\mathrm{th}} = \frac{1}{e^{h\nu/k_{\mathrm B}T}-1}.

At 10 GHz10\ \mathrm{GHz},

hνkBT≃0.0106.25=1.6×10−3,\frac{h\nu}{k_{\mathrm B}T} \simeq \frac{0.010}{6.25} =1.6\times10^{-3},

so Nth≃1/(1.6×10−3)≃625N_{\mathrm{th}}\simeq1/(1.6\times10^{-3})\simeq625. At 200 THz200\ \mathrm{THz} the exponent is about 3232, giving Nth≃e−32∼10−14N_{\mathrm{th}}\simeq e^{-32}\sim10^{-14}. Bright natural thermal noise is therefore intrinsic in the microwave regime but not usually in the optical regime.

A target is present with prior probability π1=10−3\pi_1=10^{-3}. A receiver has PFA=10−2P_{\mathrm{FA}}=10^{-2} and PM=0.1P_{\mathrm M}=0.1. Compute its equal-cost Bayes error. Which term dominates, and why can equal-prior optimization be misleading?

Solution

With π0=0.999\pi_0=0.999,

Pe=(0.999)(10−2)+(10−3)(0.1)=0.00999+0.00010=0.01009.\begin{aligned} P_e &= (0.999)(10^{-2}) + (10^{-3})(0.1) \\ &=0.00999+0.00010 \\ &=0.01009. \end{aligned}

False alarms dominate because the target-absent hypothesis is overwhelmingly common. A receiver optimized for equal priors may choose a threshold with an unacceptable false-alarm rate even if its symmetric error is small. Here the receiver even loses to always declaring the target absent, whose Bayes error is π1=0.001\pi_1=0.001. Thus the low-prevalence example demands a different decision threshold, not merely a smaller equal-prior error. The ROC or a fixed-false-alarm comparison reveals the tradeoff.

Exercise 8: Failed trials change the error exponent

Section titled “Exercise 8: Failed trials change the error exponent”

A cycle succeeds with probability p=0.1p=0.1, independent of the target. Successful records have Chernoff coefficient e−5e^{-5}, and a failed record contains no target information. Find the exponent per attempted cycle and compare it with multiplying the successful exponent by pp.

Solution

The flagged outcome coefficient is Qcycle=0.9+0.1e−5Q_{\rm cycle}=0.9+0.1e^{-5}, so

ξcycle=−ln⁡(0.9+0.1e−5)≃0.1046.\xi_{\rm cycle}=-\ln(0.9+0.1e^{-5})\simeq0.1046.

The naive product p×5=0.5p\times5=0.5 is nearly five times too large. For nn attempts, the all-failure event alone has probability 0.9n0.9^n and leaves equal-prior error 1/21/2, preventing an exponent above −ln⁡0.9≃0.1054-\ln0.9\simeq0.1054. If instead one waits for a fixed number of successes, the error and random duration must be described using that different stopping rule.

  • J. Angeletti, H. Shi, T. Lakshmanan, D. Vitali, and Q. Zhuang, “Microwave Quantum Illumination with Correlation-to-Displacement Conversion,” Physical Review Applied 20, 024030 (2023), doi:10.1103/PhysRevApplied.20.024030.

  • R. Assouly, R. Dassonneville, T. Peronnin, A. Bienfait, and B. Huard, “Quantum Advantage in Microwave Quantum Radar,” Nature Physics 19, 1418–1422 (2023), doi:10.1038/s41567-023-02113-4.

  • S. Barzanjeh, S. Pirandola, D. Vitali, and J. M. Fink, “Microwave Quantum Illumination Using a Digital Receiver,” Science Advances 6, eabb0451 (2020), doi:10.1126/sciadv.abb0451.

  • S. Guha and B. I. Erkmen, “Gaussian-State Quantum-Illumination Receivers for Target Detection,” Physical Review A 80, 052310 (2009), doi:10.1103/PhysRevA.80.052310.

  • C. W. Helstrom, Quantum Detection and Estimation Theory, Academic Press (1976), doi:10.1016/C2013-0-10310-8.

  • A. Karsa, A. Fletcher, G. Spedalieri, and S. Pirandola, “Quantum Illumination and Quantum Radar: A Brief Overview,” Reports on Progress in Physics 87, 094001 (2024), doi:10.1088/1361-6633/ad6279.

  • R. Nair and M. Gu, “Fundamental Limits of Quantum Illumination,” Optica 7, 771–774 (2020), doi:10.1364/OPTICA.391335.

  • H. Shi, B. Zhang, J. H. Shapiro, Z. Zhang, and Q. Zhuang, “Optimal Entanglement-Assisted Electromagnetic Sensing and Communication in the Presence of Noise,” Physical Review Applied 21, 034004 (2024), doi:10.1103/PhysRevApplied.21.034004.

  • S.-H. Tan, B. I. Erkmen, V. Giovannetti, S. Guha, S. Lloyd, L. Maccone, S. Pirandola, and J. H. Shapiro, “Quantum Illumination with Gaussian States,” Physical Review Letters 101, 253601 (2008), doi:10.1103/PhysRevLett.101.253601.

  • M. M. Wilde, M. Tomamichel, S. Lloyd, and M. Berta, “Gaussian Hypothesis Testing and Quantum Illumination,” Physical Review Letters 119, 120501 (2017), doi:10.1103/PhysRevLett.119.120501.

  • Z. Zhang, S. Mouradian, F. N. C. Wong, and J. H. Shapiro, “Entanglement-Enhanced Sensing in a Lossy and Noisy Environment,” Physical Review Letters 114, 110506 (2015), doi:10.1103/PhysRevLett.114.110506.

  • Q. Zhuang, Z. Zhang, and J. H. Shapiro, “Optimum Mixed-State Discrimination for Noisy Entanglement-Enhanced Sensing,” Physical Review Letters 118, 040801 (2017), doi:10.1103/PhysRevLett.118.040801.

  • Q. Zhuang, Z. Zhang, and J. H. Shapiro, “Quantum Illumination for Enhanced Detection of Rayleigh-Fading Targets,” Physical Review A 96, 020302(R) (2017), doi:10.1103/PhysRevA.96.020302.