Quantum information · Proposal

Distributed dual-comb spectroscopy

Haowei Shi; Kartikeya Arora; Zheshen Zhang; Quntao Zhuang

arXiv:2610.30462v1 ↗PDF ↗Submitted 21 October 2026

Question
Quantum-enhanced dual-comb spectroscopy uses squeezed or entangled frequency combs to beat shot noise over a broad band, and distributed quantum sensing uses entanglement across sites to estimate global quantities. Can one combine both, estimating weighted averages of absorption across many samples and many frequencies, with a single quantum light source?
Result
Yes. Each sample sits in a Mach–Zehnder interferometer with a bright classical comb in one port and a share of a single two-mode-squeezed comb, split \(K\) ways by a passive beam splitter, in the dark port. With splitter weights matched to the target weights, the variance of every weighted-average estimator falls by the squeezing gain \(G\). Splitters that do not depend on frequency force a trade-off between gases with different targets: for orthogonal targets the guaranteed advantage is \(2G/(G+1)\), at most 3 dB.

The paper in detail

In dual-comb spectroscopy a signal comb probes the sample and beats against a local-oscillator comb with a slightly different repetition rate. The beat notes at \(m\Delta\omega_r\) compress an optical spectrum into radio frequencies. Heterodyne with a strong uniform oscillator gives \(\langle\hat d_m\rangle=A_m+A_{-m}^*\) and noise \(\sum_n(\hat{\mathcal A}_{n,m}+\hat{\mathcal A}^\dagger_{n,-m})\), so squeezing the appropriate pairs of sidebands reduces the noise. Earlier work displaced an entangled comb with a classical one on a nearly transparent splitter, which wasted most of the classical power.

The single-sensor fix interferes a zero-mean squeezed comb \(\hat A^Q\) with a bright classical comb \(\hat A^C\) on a 50:50 splitter and sends one arm through the sample. At the second splitter the dark port carries \(\hat A^Q+(\mathcal L[\hat A_+]-\hat A_+)/\sqrt2\): the noise is that of the squeezed comb, and the signal \(\langle\hat d_m\rangle=\tfrac12A_m(\tau_m-1)\) uses a quarter of the classical comb energy rather than a vanishing fraction. A division receiver, which takes ratios of photocurrents, works the same way in the weak-absorption limit, up to which sidebands are squeezed and a constant factor.

For \(K\) samples the target is \(\theta_m=\sum_\ell w^{(\ell)}(\tau^{(\ell)}_m-1)\) at every comb line. The estimator \(\hat O_m=\sum_\ell v^{(\ell)}\hat d^{(\ell)}_m\) is unbiased when \(v^{(\ell)}A^{(\ell)}_m=2w^{(\ell)}\). Its noise is a weighted sum of the sensors’ quantum combs, which is fully squeezed if a single squeezed comb is split with weights \(q^{(\ell)}\propto w^{(\ell)}\). Under a total power budget, optimal comb amplitudes \(|A^{(\ell)}|^2\propto|w^{(\ell)}|\) give \(\mathrm{Var}=\tfrac{2\hbar\Omega_c}{P_mT}\bigl(\sum_n G_n^{-1}\bigr)\bigl(\sum_\ell|w^{(\ell)}|\bigr)^2\). For a gas-density gradient this sends all power to the two end sensors; with equal power per sensor the optimum is the least-squares slope, at most a factor 3 worse.

When several gas species need different functionals, such as a background average and a leak gradient, one frequency-independent splitter \(\mathbf q\) must serve them all. The advantage for species \(g\) is \(\mathcal Q_g=[1-(1-1/G)(\tilde{\mathbf p}_g\cdot\mathbf q)^2]^{-1}\), so it depends on one direction cosine. The reachable alignments lie in an ellipse fixed by the overlap of the targets. The symmetric point gives \(\mathcal Q_\star=G\) for collinear targets but only \(2G/(G+1)\) for orthogonal ones, however strong the squeezing. Methods extend the analysis to division receivers, whose optimum is the same configuration with twice the heterodyne variance because they need two squeezed sources. They also show that the outputs can be recombined before a single detector.

Transferable methods

  1. Match the entanglement splitter to the estimated linear functional

    The noise of a weighted sum of sensor outputs is a projection of the network’s quantum noise onto the weight vector. Distributing one squeezed mode along exactly that direction squeezes the estimator, and the other ports contribute only vacuum noise orthogonal to it.

    With unused ports in vacuum, \(\mathrm{cov}(\hat S^{(\ell)},\hat S^{(\ell')})=M\delta_{\ell\ell'}+q^{(\ell)}q^{(\ell')}(V-M)\), \(V=\sum_nG_n^{-1}\), so \(\mathrm{Var}\,\hat O=4[M\|\mathbf p\|^2-(M-V)(\mathbf p\cdot\mathbf q)^2]\) with \(\mathbf p=(w^{(\ell)}/A^{(\ell)})_\ell\). By Cauchy–Schwarz the minimum is at \(\mathbf q=\mathbf p/\|\mathbf p\|\), giving \(4V\|\mathbf p\|^2\), an advantage \(M/V=G\) for uniform gain.

    Gradient example (editorial check of Eqs. 27 and 30): endpoint weights give variance \(\propto64/(K-1)^2\) and least-squares weights under equal power give \(\propto192/(K^2-1)\). Their ratio \(3(K-1)/(K+1)\) is 1.8 at \(K=4\) and approaches 3 as \(K\to\infty\).

    How to reuse it

    For any distributed estimate of a linear combination of local parameters, write the estimator noise as a quadratic form, identify the direction the target projects onto, and split one nonclassical resource along it. The advantage is set by the squared overlap between splitter and target direction.

  2. Bound multi-target advantage by a shared-resource geometry

    When one fixed splitter must serve several targets, each advantage depends only on the squared overlap of the splitter with that target’s direction. Unit vectors cannot overlap strongly with two nearly orthogonal directions at once, which caps the advantage that can be guaranteed to all targets.

    For two targets, \(c_g=(\tilde{\mathbf p}_g\cdot\mathbf q)^2\) obey \(c_1+c_2-2\sqrt{c_1c_2}\,|\tilde{\mathbf p}_1\cdot\tilde{\mathbf p}_2|\le1-(\tilde{\mathbf p}_1\cdot\tilde{\mathbf p}_2)^2\). The symmetric optimum \(c_1=c_2=(1+|\tilde{\mathbf p}_1\cdot\tilde{\mathbf p}_2|)/2\) gives \(\mathcal Q_\star=[1-(1-1/G)(1+|\tilde{\mathbf p}_1\cdot\tilde{\mathbf p}_2|)/2]^{-1}\).

    Editorial check at \(G=15\) dB: \(\mathcal Q_\star\) is 15.0 dB for collinear targets, 5.6 dB for overlap 0.5 and 2.9 dB for orthogonal targets. A random search over unit \(\mathbf q\in\mathbb R^4\) for the orthogonal pair (average, least-squares gradient) finds a best guaranteed squared alignment of 0.499, consistent with the bound 1/2.

    How to reuse it

    When a single quantum resource must serve several estimation tasks, express each task’s gain through its overlap with the resource direction. Then the achievable trade-off is a question about unit vectors and a Gram matrix, and it tells you in advance when frequency- or task-dependent routing is worth building.

Source: Sections II–IV (Eqs. 1–21); Section V (Eqs. 22–34); Methods VIII.1 (Eqs. 35–48). Based on the linked paper version. Derivations marked as editorial, and suggested uses in other problems, are editorial guidance rather than results of the paper.

Context & limitations

Why it matters. It shows that spatial and spectral quantum advantage can be obtained together from one engineered comb, where a separable strategy would need \(K\) independent quantum sources. The interferometric design also fixes the low-brightness problem of earlier entangled-comb proposals, in which most of the classical comb power was wasted.

Limitations. The analysis is in the weak-absorption limit, with ideal lossless distribution and uniform squeezing. Loss, phase noise and finite detector efficiency are deferred to supplementary material or future work. Full advantage for several gas species at once needs frequency-dependent beam splitters that do not yet exist. No experiment is reported.

Description based on the full text of the linked arXiv version, read 2026-10-24; equations follow the paper’s notation.