Squeezed Light and Gaussian Boson Sampling

Boson sampling's Achilles heel is its fuel. It needs n identical single photons in the same clock cycle, and a heralded SPDC source delivers one only a few percent of the time — so the n-fold coincidence rate collapses as p_{\text{source}}^{\,n}. Twenty photons: essentially never. The escape, found by Hamilton and colleagues in 2017, is a change of fuel rather than of engine. The same nonlinear crystal that makes heralded pairs makes, when you simply don't herald, squeezed light — a state that pours out deterministically, every single pump pulse. Feed squeezed light straight into the mesh and count photons at the output: this is Gaussian boson sampling (GBS), the design behind the first photonic quantum-advantage claims — China's Jiuzhang and Xanadu's Borealis. To understand it we need this course's one excursion into continuous-variable quantum optics: field quadratures, the vacuum's noise floor, and what it means to squeeze below it.

Quadratures: the harmonic oscillator returns

A single optical mode is, mathematically, a quantum harmonic oscillator — the photon number is the oscillator's level. Its "position" and "momentum" are the field's quadratures,

\hat X = \tfrac{1}{\sqrt2}\bigl(\hat a + \hat a^\dagger\bigr), \qquad \hat P = \tfrac{1}{i\sqrt2}\bigl(\hat a - \hat a^\dagger\bigr), \qquad [\hat X, \hat P] = i ,

the cosine and sine components of the oscillating field — what a radio engineer calls I and Q. Being conjugate, they obey an uncertainty relation \Delta X\,\Delta P \ge \tfrac12. The vacuum saturates it symmetrically: \Delta X^2 = \Delta P^2 = \tfrac12. This vacuum fuzz is not an instrument artefact — it is the shot noise that Module 5 met as the fundamental noise floor of optical detection: measure any quadrature of "nothing" and you get a Gaussian jitter of variance ½. Laser light (a coherent state) is just displaced vacuum, carrying the same circle of noise around a nonzero mean.

Squeezing: reshaping the vacuum

The uncertainty product is a floor on the product — not on either factor. A squeezed state trades the two off. Pumping a χ² crystal applies the squeezing operator S(r) and turns the noise circle into an ellipse:

\Delta X^2 = \tfrac12\,e^{-2r}, \qquad \Delta P^2 = \tfrac12\,e^{+2r},

quieter than the vacuum along one axis, correspondingly louder along the other. The squeezing strength is quoted in decibels below shot noise, -10\log_{10}\!\bigl(e^{-2r}\bigr) \approx 8.69\,r dB; today's best crystals reach about 15 dB (r \approx 1.7). Two facts make squeezed vacuum the perfect GBS fuel. First, it is deterministic: every pump pulse squeezes, no herald, no lottery. Second, expand it in photon number and only even terms appear,

|r\rangle \;=\; \frac{1}{\sqrt{\cosh r}} \sum_{n=0}^{\infty} \frac{\sqrt{(2n)!}}{2^n n!}\,(-\tanh r)^{\,n}\,|2n\rangle ,

because the crystal creates photons strictly in pairs — squeezed vacuum is the un-heralded face of SPDC, with mean photon number \sinh^2 r. That pair structure is about to leave a fingerprint on the mathematics.

GBS: sampling from hafnians

The Gaussian boson sampler squeezes several input modes of an m-mode mesh, lets the unitary interfere them, and counts photons at every output. The output probabilities are again given by a matrix function of the mesh unitary — but the permanent, which pairs inputs with outputs one-to-one, is replaced by the function that pairs output photons with each other, respecting their pairwise birth:

P(S) \;\propto\; \bigl|\operatorname{Haf}(A_S)\bigr|^2 ,

where A is a matrix built from U and the squeezing strengths, and the hafnian Haf sums over all perfect matchings — all ways of splitting the detected photons into pairs. (The permanent is recoverable as a special case, so the hafnian is the harder object; it is likewise #P-hard, with the best classical algorithms scaling roughly as 2^{n/2}.) The complexity story then runs parallel to Aaronson–Arkhipov: efficient classical GBS simulation would collapse the polynomial hierarchy, with analogous conjectures covering the noisy case. What changes is the economics: with deterministic fuel in every input port, experiments jumped from boson sampling's handful of photons to hundreds, essentially overnight.

Hence the record books. Jiuzhang (USTC, December 2020): 50 squeezed states into a 100-mode free-space interferometer, up to 76 photons detected, and a sampling task estimated at the time to need billions of years classically — the first quantum-advantage claim on light. Borealis (Xanadu, 2022): the same mathematics rebuilt as engineering — a single squeezed source firing into a loop of fibre delay lines, 216 time-bin modes interfered by three dynamically programmed beamsplitters, and quantum advantage claimed on a programmable, fibre-networked machine you could queue jobs to over the cloud. Between them sits the pattern this course keeps finding: the physics fits on one line; the leap is in the packaging.

Loss: the advantage-eating channel

Now the villain. Every photonic system loses photons, but for squeezed light loss is crueller than a dimmer switch. Transmission \eta acts as a beamsplitter that swaps a fraction 1-\eta of your state for vacuum — in the language of quantum channels, a Gaussian channel mixing in fresh noise. A quadrature variance V becomes

V \;\longrightarrow\; \eta\,V + (1 - \eta)\cdot\tfrac12\cdot 2 \;\;\text{(in shot-noise units: } V \to \eta V + 1 - \eta\text{)} ,

so the squeezed variance e^{-2r} is dragged back toward the vacuum's 1. The floor is brutal: even infinite squeezing behind transmission \eta can never show a variance below 1-\eta. At 50% transmission, no source in the universe looks better than 3 dB. Play with the slider — watch how loss crushes the squeezed curve toward vacuum while barely denting the noisy one:

For GBS machines this is existential. Squeezing is the quantum resource; as loss mounts, the output state slides toward a classically simulable thermal state, and classical "spoofing" algorithms exploit exactly this — several later analyses narrowed or contested the early advantage claims by simulating the lossy experiment rather than the ideal one. The lesson generalises far beyond sampling: in photonic quantum computing, loss is not one error among many but the error, the quantity every architecture is shaped around. The final lesson takes that as its founding constraint.

Worked example: a squeezing budget

Your crystal produces 10 dB of squeezing: e^{-2r} = 0.1. The route to the detector — coupling, propagation, detection efficiency — has total transmission \eta = 0.5. The observed variance is

V \;=\; 0.5 \times 0.1 + 0.5 \;=\; 0.55 \quad\Longrightarrow\quad -10 \log_{10} 0.55 \approx 2.6\ \text{dB} .

Ten decibels of hard-won laboratory squeezing arrive as 2.6. Run the budget the other way and the engineering targets write themselves: to observe even 6 dB (V = 0.25) from a perfect source you need \eta > 0.75; state-of-the-art squeezing experiments fight for every percent of detection efficiency, and GBS machines live or die on end-to-end transmission in a way qubit machines — which can retry a lost photon — do not.

Jiuzhang is named for the Jiuzhang Suanshu — "Nine Chapters on the Mathematical Art" — the Han-dynasty text that anchored two millennia of Chinese mathematics, a pointed answer to Google naming its rival chip after a computing pioneer. The machine itself was gloriously un-chip-like: a room-sized table of free-space optics, 25 crystals pumped in lockstep, and a 100-mode interferometer built from mirrors and beamsplitters aligned to sub-wavelength tolerances — phase stability maintained not by a foundry but by heroic optomechanics. Its 2020 run detected up to 76 photons across 100 modes. The advantage claim aged the way such claims do: classical spoofers exploiting loss and photon-number structure steadily shrank the gap, USTC answered with bigger runs (Jiuzhang 2.0 and 3.0), and the true legacy settled elsewhere — GBS at the hundred-photon scale was proven physically real, and the pressure it created produced the tensor-network and matchgate spoofing algorithms that now define the classical frontier. Advantage, it turns out, is not a finish line but an arms race.

For a laser beam, 50% loss means half the power and an unchanged signal-to-noise character; for a squeezed state, 50% loss means most of the quantumness is gone. The mechanism is the part to internalise: the lost fraction is replaced by vacuum, and vacuum is noisy — variance 1 in shot-noise units. The channel therefore averages your carefully quietened e^{-2r} with 1, and the average is dominated by whichever is larger: beyond mild squeezing, more r buys almost nothing (\eta e^{-2r} + 1 - \eta \to 1 - \eta). Two traps follow. First, "we'll just squeeze harder" cannot beat a lossy channel — only better transmission can. Second, the same arithmetic read backwards is a diagnostic: measure 3 dB where your source makes 10, and you have measured your losses. And because the anti-squeezed quadrature keeps growing as \eta e^{2r}, a lossy strongly-pumped source is worse than useless for sensing — all the excess noise, little of the quiet.

Where this goes next

GBS closed the argument that photons can do classically impossible things at scale — and underlined, twice, what stands between here and a useful machine: loss, and the absence of error correction in sampling architectures. The final lesson of this module assembles everything — dual-rail qubits, HOM interference, KLM-style heralded entanglement, and loss-as-erasure — into the architecture that industry has actually bet on: fusion-based quantum computing, where small photonic resource states are stitched into a fault-tolerant whole by Bell measurements, on silicon photonics stamped out in a semiconductor fab.