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.