simulation

Computable measures of fermionic non-Gaussianity from the covariance matrix

Curator's Take

This article delivers the first practical resource‑theoretic tools for quantifying fermionic “magic,” turning a previously abstract notion into computable monotones that can be read directly from a state’s covariance matrix. By linking occupation‑number and natural‑orbital entropies to gate counts and classical simulation costs, it bridges the gap between many‑body physics and algorithmic complexity—a connection that complements recent work on Gaussian circuit simulability and magic‑state distillation for qubits. The results give researchers a concrete way to assess when fermionic systems become classically intractable, which could guide both hardware design and the development of hybrid quantum‑classical algorithms. As with any resource theory, the measures are most informative for pure states; extending them to mixed or noisy settings will be an important next step.

— Mark Eatherly

Summary

Fermionic non-Gaussianity, or fermionic magic, is a key resource underlying the computational complexity of fermionic quantum systems, yet tractable and operationally meaningful ways to quantify it remain limited. We address this challenge by developing a convex resource theory of fermionic non-Gaussianity and introducing two families of computable measures for pure fermionic states, both derived from the Williamson normal form of the covariance matrix. The first family, occupation number entropies, is defined as the Tsallis-$α$ entropy of the occupation numbers. We prove that one member of this family is monotonic under Gaussian protocols, establishing it as a computable convex resource monotone. It consequently lower bounds the number of non-Gaussian gates needed for state preparation. The second family, natural-orbital participation entropies, is given by the Rényi-$α$ entropy of the squared amplitudes of the state in the natural-orbital basis, defined by the eigenvectors of the covariance matrix. These measures quantify state compressibility in this basis and thus upper bound the classical simulation cost in an orthonormal Gaussian basis. We analyze both families for stabilizer and translation-invariant states, where they simplify and reveal additional structure. We further study representative examples, including random SWAP-doped matchgate circuits and the bond-modulated XXZ model, highlighting the role of non-Gaussianity in many-body phenomena. Our work establishes a resource-theoretic framework for computable fermionic non-Gaussianity that unifies notions arising across quantum information, condensed-matter physics, and quantum chemistry, opening new directions for studying the complexity of quantum many-body systems and providing practical tools to assess the classical simulability of fermionic states relevant for quantum advantage.