Imagine you're trying to forge a signature by layering simple geometric stamps — circles, rectangles, straight lines. For most signatures, you can get reasonably close with a modest stack. But some signatures have a property — a particular curlicue, an irreducible flourish — that forces you to use exponentially more stamps as you scale up. No clever arrangement helps. This paper proves that non-Gaussianity is exactly that irreducible flourish for quantum states: if a state isn't Gaussian, then approximating tensor-power copies of it from Gaussian building blocks requires an exponentially growing number of terms, period. The committed claim: for every pure non-Gaussian state (with definite parity for fermions), the approximate border Gaussian rank of n copies grows exponentially in n at any fixed approximation error below one. This is not a result about specific exotic states — it is universal across both fermionic and bosonic systems. The authors prove this for every such state on finitely many modes, which is a substantially stronger statement than prior work that handled only specific examples or weaker rank notions. The proof architecture is elegant in its modularity. For fermions, the authors reduce the problem to four modes via a clever projection argument, then deploy Majorana spectral bounds — eigenvalue constraints on the state's representation in the Majorana operator basis — to establish the exponential floor. For bosons, the strategy is different: Gaussian postselection combined with an entropy-based rank bound. The bosonic result is notable because it requires no energy assumption on the target state, sidestepping a common limitation in continuous-variable quantum information. Concrete examples anchor the abstract universality. The four-mode fermionic GHZ state and the bosonic single-photon state both receive explicit exponential lower bounds. These are canonical non-Gaussian states, so the examples serve as proof-of-concept that the general theorem bites on states people actually care about in quantum computing and quantum optics. The integrity profile is that of a mathematical proof paper, not an experimental or computational one. The results are theorems, not simulations. The validation is internal mathematical rigor — reduction steps, spectral arguments, entropy inequalities — not benchmarking against external datasets. This is the appropriate validation mode for this kind of result, but it means the strength rests entirely on the correctness of the proofs, which awaits independent verification by the community. Why does this matter beyond the math? Decomposition-based classical simulation is one of the primary strategies for simulating quantum circuits on classical hardware. Every time you want to classically simulate a non-Gaussian quantum computation, you need to decompose states into sums of Gaussian ("free") states. The number of terms in that sum — the Gaussian rank — directly controls how long your classical simulation takes. This paper says that cost is inescapably exponential for non-Gaussian resources, no matter how clever your decomposition strategy. It draws a hard line: non-Gaussianity universally implies classical simulation hardness via this route. The result is a direct analog of the qubit stabilizer rank problem, which has been notoriously difficult to resolve with tight lower bounds. By establishing robust exponential lower bounds in the Gaussian setting — for both fermions and bosons, approximately and not just exactly, and for border rank not just standard rank — this paper sets a new high-water mark for what can be proven about decomposition complexity. The open question now is whether these techniques can be adapted back to the qubit stabilizer setting or extended to mixed states.
