A Yardstick for Quantum Advantage: Exact Simulation Frontiers Mapped Across Four Photonic Architectures

Open-source GPU library benchmarks exactly where classical simulation of photonic quantum computers breaks down — 38 single photons on one H100 — giving the field a shared ruler instead of competing claims.

Imagine you're timing how long it takes to shuffle a deck of cards into every possible order. With 10 cards it's fast; with 52 it outlasts the universe. The question isn't whether it's hard — it's exactly where the crossover happens between 'my laptop handles this' and 'no classical machine on Earth can keep up.' That's the question this paper answers for photonic quantum computing, and the answer is more precise than anything the field has had before. The committed claim: there is a single, well-defined 'simulability frontier' for photonic quantum devices — the device size at which the fastest exact classical sampler on a fixed GPU produces one sample per second. The authors implement this for four distinct light sources (single photons, entangled photons, vacuum-single-photon superpositions, and squeezed light) in an open-source GPU library called photonic-engines, validated against independent simulators sharing no code. On a single NVIDIA H100, the frontier sits at 38.3 single photons, 31.0 entangled photons from one emitter, 41.5 squeezed photonic inputs (mean photon number 0.5), and — critically — 50.0 for coherent vacuum-photon superpositions versus 64.2 for the same sources mixed incoherently. Coherence makes simulation harder, not easier. The ladder here matters. Previous claims about when boson sampling beats classical computation range wildly — 50 to 98 photons, depending on who you ask and what they assumed. The divergence is not about disagreement over physics; it's about disagreement over the metric. Some groups measure approximate sampling (whose difficulty depends on required accuracy), others use different hardware baselines, still others compare to algorithms that were never implemented. This paper cuts through the noise by defining a single, implementable measure: exact sampling time per sample, on real hardware, with real code. It doesn't claim quantum advantage — it tells you exactly how far away it is for each architecture. Architecturally, the samplers exploit the permanent/hafnian computation family that underlies boson sampling theory, implemented as GPU-accelerated exact algorithms. The key structural insight is that simulation cost scales exponentially in detected photon number, so doubling classical compute only pushes the frontier by about one photon. This is the central bottleneck: better hardware buys almost nothing, while better algorithms change the exponential growth rate itself. The library covers Fock-state (single photon), Gaussian (squeezed), and hybrid regimes, making it the first apples-to-apples comparison across photonic architectures. Integrity is a genuine strength. Each sampler is cross-validated against a simulator written by a different codebase — not the same team running the same code twice. The code is open-source. The benchmarks use Haar-random interferometers (the standard hard-instance distribution), twice as many modes as photonic inputs, and photon-number-resolving detectors. The authors are transparent that their frontier numbers will move as algorithms improve, and they explicitly parameterize the frontier as a function of hardware and algorithmic advances. The one gap: no independent group has yet replicated the measurements, and the squeezed-light frontier relies on a timing model fitted to measurements rather than direct computation at the boundary. The milestone framing is unusually clean. Today's frontier: 38.3 single photons on one H100. The field needs roughly 50–100 detected photons to demonstrate unambiguous quantum advantage for any practically interesting sampling task. Each doubling of classical compute buys only ~1 photon, so the frontier moves slowly from the classical side. The real race is whether photonic hardware can scale detected photon counts faster than algorithmic improvements push the frontier upward. The paper gives both sides of the race a shared scoreboard. The obvious experiment not run: scaling to multi-GPU or cluster-scale classical simulation. The paper benchmarks a single H100; a 256-GPU cluster would push the frontier up by about 8 photons (log₂ of the speedup). The authors almost certainly have access to multi-GPU setups — this omission likely reflects a deliberate choice to define the frontier on a reproducible single-card baseline, saving the multi-GPU extension for a follow-up that would let them argue for even higher classical thresholds.