Imagine you're trying to learn the shape of a bell curve by drawing samples from it. In the classical world, if the distribution is a 2n-dimensional Gaussian, you need roughly n² samples to pin down its covariance matrix — that's a well-known result from statistics. Now imagine the same bell curve, but drawn on a quantum canvas: a bosonic Gaussian state of n modes, which has a covariance matrix of the same size but lives in a Hilbert space governed by the uncertainty principle. The question this paper answers is: does the quantum canvas make the learning problem harder, and if so, when? The committed claim: there exists a sharp, temperature-controlled crossover in the sample complexity of learning bosonic Gaussian states. When the state is 'cold' — its covariance matrix is close to the vacuum noise floor of ½I — you provably need Ω(n³) copies under any single-copy (non-entangled) measurement strategy, which is strictly more than the Θ(n²) samples a classical learner would need for the same-dimensional Gaussian. When the state is 'warm' — thermal fluctuations exceed the vacuum noise by any constant amount ν — the sample complexity drops to Θ(n²), matching the classical case exactly, achievable with simple heterodyne measurements. The architectural move here is information-theoretic, not algorithmic. This is not a paper proposing a new tomography protocol and benchmarking it against competitors. It's a paper proving matching upper and lower bounds on sample complexity as a function of the thermal parameter ν. The lower bounds use quantum information-geometric techniques (quantum Fisher information, Le Cam-style reductions). The upper bound is achieved by heterodyne detection — arguably the simplest possible measurement on a bosonic mode — making the achievability result practically significant. The key structural insight is that the parameter ν controls a smooth interpolation: the sample complexity is Θ(n² · min(n, 1 + ν⁻¹)), which is n³ when ν → 0 and n² when ν = Ω(1). On the ladder, this paper advances a line of work on quantum state tomography sample complexity. Prior results established that general quantum state tomography scales exponentially, and that structure (like Gaussianity) helps enormously. For bosonic Gaussian states specifically, the Θ(n²) classical baseline was known from classical statistics (learning a 2n × 2n covariance matrix). The n³ lower bound for cold states and the tight characterization of the crossover are new. The closest prior work includes results on finite-dimensional Gaussian state tomography (Holevo, Serafini, and more recently Ohliger et al.), but the tight ν-parameterized crossover and the matching bounds across the full temperature range appear to be a genuine first. Integrity is strong for this category of work. These are mathematical proofs — upper and lower bounds — not empirical benchmarks. The validation is as non-circular as it gets in theoretical CS/physics: you either prove matching bounds or you don't. There's no cherry-picking of benchmarks, no hyperparameter tuning, no test-set leakage. The few-copy entangled measurement result (showing hardness persists even with limited entanglement) strengthens the lower bound side. The main caveat is that the results assume idealized measurements and noiseless state preparation — but that's standard for information-theoretic bounds and is clearly scoped. The milestone question is interesting because this paper is more of a foundational characterization than a stepping stone on a hardware roadmap. The practical implication is for quantum sensing and continuous-variable quantum information processing: experimentalists doing optical or microwave quantum tomography now have a precise thermometer telling them when their problem is fundamentally quantum-hard versus classically easy. The next concrete milestone would be an experimental demonstration of the n³ → n² crossover — take an n-mode optical system, cool it to near vacuum, show the extra sample cost, then warm it and watch the cost drop. That would validate the theory in a real lab setting. The obvious experiment not run: an adaptive measurement strategy for the cold regime. The paper shows that non-adaptive, non-entangled measurements require n³ copies in the cold regime, and that even few-copy entanglement doesn't help. But what about fully adaptive strategies where each measurement is chosen based on prior outcomes? The authors likely didn't close this gap because it requires different proof techniques (adaptive information-theoretic lower bounds are notoriously harder), and the current result is already a clean, complete story for the non-adaptive case. This feels like (c) — a deliberate scope choice leaving the adaptive question for the next paper.