Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2605.18986 · Europhysics Letters (EPL)

Non-Gaussianity of random quantum states

Filiberto Ares, Sara Murciano, Pasquale Calabrese

cond-mat.stat-mechquant-phSubmitted May 18, 2026• Announced May 20, 2026• 0 citations
arXivHTMLPDF

Abstract

We study the fermionic non-Gaussianity in typical quantum states, focusing on Haar random states of qubits with or without a global $U(1)$ symmetry. Using the Weingarten calculus, we derive analytical predictions for the non-Gaussianity, defined as the relative entropy between the reduced density matrix and its Gaussianized counterpart. We identify two regimes controlled by the ratio between the subsystem and the system size, $\ell/L$. For $\ell/L < 1/2$, the non-Gaussianity vanishes in the absence of symmetries, because typical reduced density matrices are exponentially close to the maximally mixed state. In the presence of a global $U(1)$ symmetry, instead, it remains small but finite. By contrast, in the regime $\ell/L > 1/2$, the non-Gaussianity becomes extensive. These results establish the typical scaling of fermionic non-Gaussianity in random states and analyze how this is modified by the presence of global symmetries.

Journal reference
EPL 156 (2026) 38001
DOI
10.1209/0295-5075/ae9219

Related papers

Universal Bound and Phase Transition in Many-Body Fermionic Non-Gaussianity

Masahiro Hoshino, Ryota Matsuda, Yuto Ashida

Fermionic non-Gaussianity is a resource for universal quantum computation that can be generated by interactions in quantum many-body systems. Using the magic Rényi entropy (MRE) as a measure of non-Gaussianity, we derive its universal upper bound and Haar mean, proving that typical Haar-random states attain the maximal MRE density of $\ln(4/3)$ per Majorana in the large-system limit. We also show that the MRE can differ extensively between states with identical covariance matrices. To investigate how non-Gaussianity grows toward this maximal density, we evolve a Gaussian state in imaginary and real time under the Sachdev-Ye-Kitaev Hamiltonian. By varying the imaginary- and real-time durations, we identify a first-order transition marked by a kink in the MRE density and spontaneous breaking of the permutation symmetry among the four copies used to evaluate the MRE. This transition represents a qualitative change in the non-Gaussianity of the state that is not reflected in the thermal free energy.

quant-pharXiv:2610.02075•Oct 2, 2026•0 citations

Pauli spectrum and nonstabilizerness of random fermionic Gaussian states

Xhek Turkeshi, Piotr Sierant, Poetri Sonya Tarabunga

We characterize the magic of random fermionic Gaussian states through their Pauli spectrum and their stabilizer entropies. For the Haar ensemble of Majorana Gaussian states we derive closed finite-$N$ expressions for the average stabilizer purities and reconstruct the Pauli spectrum exactly, as a mixture of products of beta-distributed variables resolved by Majorana weight. The stabilizer entropies, and their filtered versions, freeze at Rényi index $q_c=2$: for $q>2$ the magic density is $1/(q-1)$, controlled by rare weight-two Majorana strings. Resolving these contributions requires exponentially many characteristic samples, making direct sampling estimates of the stabilizer entropies inefficient in this frozen regime. For Haar-random Gaussian states at fixed filling, we express the averaged stabilizer purities at positive integer $q$ as coefficient integrals over $\mathrm{SU}(2q)$, with dimension independent of $N$, and evaluate the $q=2$ case exactly. We establish that the frozen density describes typical states in the Majorana and half-filled number-conserving ensembles and relate these ensembles to the corresponding $\mathrm{SYK}_2$ ground states. Finally, we show that typical fermionic Gaussian states retain certified magic after discarding any fixed fraction $κ<κ_\star^{\rm G}\simeq0.7613$ of their modes in the thermodynamic limit, remaining magical even beyond the $2/3$ threshold for Haar-random states in the full Hilbert space.

quant-pharXiv:2609.15739•Sep 15, 2026•0 citations

Random matrix prediction of average entanglement entropy in non-Abelian symmetry sectors

Anwesha Chakraborty, Lucas Hackl, Mario Kieburg

We study the average bipartite entanglement entropy of Haar-random pure states in quantum many-body systems with global $\mathrm{SU}(2)$ symmetry, constrained to fixed total spin $J$ and magnetization $J_z = 0$. Focusing on spin-$\tfrac12$ lattices and subsystem fractions $f < \frac{1}{2}$, we derive a asymptotic expression for the average entanglement entropy up to constant order in the system volume $V$. In addition to the expected leading volume law term, we prove the existence of a $\frac{1}{2}\log V$ finite-size correction resulting from the scaling of the Clebsch-Gordon coefficients and compute explicitly the $O(1)$ contribution reflecting angular-momentum coupling within magnetization blocks. Our analysis uses features of random matrix ensembles and provides a fully analytical treatment for arbitrary spin densities, thereby extending Page type results to non-Abelian sectors and clarifying how $\mathrm{SU}(2)$ symmetry shapes average entanglement.

quant-pharXiv:2512.22942•Dec 30, 2025•0 citations

Non-Local Magic Resources for Fermionic Gaussian States

Daniele Iannotti, Beatrice Magni, Riccardo Cioli +2 more

Entanglement and magic are fundamental resources that capture the complexity of quantum many-body systems. Non-local magic isolates the irreducible nonstabilizerness intrinsically tied to entanglement. However, evaluating this quantity generally requires a prohibitive minimization over the full Hilbert space, making it computationally inaccessible beyond a few qubits. Here, we overcome this bottleneck by establishing a closed-form expression for the non-local stabilizer entropies of fermionic Gaussian states over local Gaussian unitaries, which we prove at Rényi index $α=2$ for arbitrary subsystem size, and which can be evaluated in polynomial time directly from the eigenvalues of the reduced Majorana covariance matrix. We apply this framework to characterize fermionic non-local magic across diverse physical regimes: we derive an exact Page-like curve for typical random states, reveal logarithmic scaling at the quantum critical point of the XY model, and establish a quasiparticle picture for magic generation during out-of-equilibrium quantum quenches. Crucially, because our result relies solely on two-point correlation functions, it provides a scalable route for the experimental estimation of fermionic non-local magic in large-scale quantum processors via fermionic shadow tomography.

quant-pharXiv:2604.27049•May 1, 2026•0 citations

One-Body Purity, Non-Gaussianity, and Entanglement in Interacting Integrable Models

Rafał Świętek, Maksymilian Kliczkowski, Lev Vidmar +1 more

When describing entanglement in typical midspectrum eigenstates of many-body lattice Hamiltonians, two paradigms have emerged that capture the behavior observed in integrable and nonintegrable systems, Haar-random fermionic Gaussian states and Haar-random pure states, respectively. Remarkably, the former capture the behavior of interacting integrable systems, whose eigenstates are non-Gaussian. We argue that the paradigm that captures both the entanglement properties and the lack of Gaussianity in integrable systems is that of random superpositions of polynomially many Gaussian states. In contrast, eigenstates of nonintegrable systems are consistent with being described by random superpositions of exponentially many Gaussian states. We gain this understanding by comparing analytical and numerical results for the one-body purity, the non-Gaussianity, and the entanglement entropy of the random superpositions and the Hamiltonian eigenstates.

quant-pharXiv:2607.01326•Jul 3, 2026•0 citations

An asymmetry lower bound on fermionic non-Gaussianity

Filiberto Ares, Michele Mazzoni, Sara Murciano +3 more

Fermionic Gaussian states are a fundamental tool in many-body physics, faithfully representing non-interacting quantum systems and allowing for efficient numerical simulations. Given a many-body wave function, it is therefore interesting to ask how much it differs from that of a Gaussian state, as quantified by the notion of non-Gaussianity. In this work, we relate measures of non-Gaussianity with the Shannon entropy of the particle-number distribution, coinciding with the particle-number asymmetry for pure states. We derive a lower bound on the relative entropy of non-Gaussianity in terms of the exponential of the Shannon entropy, and study numerically its tightness for large system sizes. Our bound is non-trivial for large values of the asymmetry and relies on the concentration of the particle-number distribution of (mixed) fermionic Gaussian states. Since the Shannon entropy of the particle-number distribution is often efficient to compute or experimentally measure, our results can be viewed as a practical way to lower bound non-Gaussianity, highlighting a non-trivial interplay with particle-number asymmetry.

quant-pharXiv:2603.16762•Mar 18, 2026•0 citations
We study the fermionic non-Gaussianity in typical quantum states, focusing on Haar random states of qubits with or without a global U(1)U(1)U(1) symmetry. Using the Weingarten calculus, we derive analytical predictions for the non-Gaussianity, defined as the relative entropy between the reduced density matrix and its Gaussianized counterpart. We identify two regimes controlled by the ratio between the subsystem and the system size, ℓ/L\ell/Lℓ/L. For ℓ/L<1/2\ell/L < 1/2ℓ/L<1/2, the non-Gaussianity vanishes in the absence of symmetries, because typical reduced density matrices are exponentially close to the maximally mixed state. In the presence of a global U(1)U(1)U(1) symmetry, instead, it remains small but finite. By contrast, in the regime ℓ/L>1/2\ell/L > 1/2ℓ/L>1/2, the non-Gaussianity becomes extensive. These results establish the typical scaling of fermionic non-Gaussianity in random states and analyze how this is modified by the presence of global symmetries.
Fermionic non-Gaussianity is a resource for universal quantum computation that can be generated by interactions in quantum many-body systems. Using the magic Rényi entropy (MRE) as a measure of non-Gaussianity, we derive its universal upper bound and Haar mean, proving that typical Haar-random states attain the maximal MRE density of ln⁡(4/3)\ln(4/3)ln(4/3) per Majorana in the large-system limit. We also show that the MRE can differ extensively between states with identical covariance matrices. To investigate how non-Gaussianity grows toward this maximal density, we evolve a Gaussian state in imaginary and real time under the Sachdev-Ye-Kitaev Hamiltonian. By varying the imaginary- and real-time durations, we identify a first-order transition marked by a kink in the MRE density and spontaneous breaking of the permutation symmetry among the four copies used to evaluate the MRE. This transition represents a qualitative change in the non-Gaussianity of the state that is not reflected in the thermal free energy.
We characterize the magic of random fermionic Gaussian states through their Pauli spectrum and their stabilizer entropies. For the Haar ensemble of Majorana Gaussian states we derive closed finite-NNN expressions for the average stabilizer purities and reconstruct the Pauli spectrum exactly, as a mixture of products of beta-distributed variables resolved by Majorana weight. The stabilizer entropies, and their filtered versions, freeze at Rényi index qc=2q_c=2qc​=2: for q>2q>2q>2 the magic density is 1/(q−1)1/(q-1)1/(q−1), controlled by rare weight-two Majorana strings. Resolving these contributions requires exponentially many characteristic samples, making direct sampling estimates of the stabilizer entropies inefficient in this frozen regime. For Haar-random Gaussian states at fixed filling, we express the averaged stabilizer purities at positive integer qqq as coefficient integrals over SU(2q)\mathrm{SU}(2q)SU(2q), with dimension independent of NNN, and evaluate the q=2q=2q=2 case exactly. We establish that the frozen density describes typical states in the Majorana and half-filled number-conserving ensembles and relate these ensembles to the corresponding SYK2\mathrm{SYK}_2SYK2​ ground states. Finally, we show that typical fermionic Gaussian states retain certified magic after discarding any fixed fraction κ<κ⋆G≃0.7613κ<κ_\star^{\rm G}\simeq0.7613κ<κ⋆G​≃0.7613 of their modes in the thermodynamic limit, remaining magical even beyond the 2/32/32/3 threshold for Haar-random states in the full Hilbert space.
We study the average bipartite entanglement entropy of Haar-random pure states in quantum many-body systems with global SU(2)\mathrm{SU}(2)SU(2) symmetry, constrained to fixed total spin JJJ and magnetization Jz=0J_z = 0Jz​=0. Focusing on spin-12\tfrac1221​ lattices and subsystem fractions f<12f < \frac{1}{2}f<21​, we derive a asymptotic expression for the average entanglement entropy up to constant order in the system volume VVV. In addition to the expected leading volume law term, we prove the existence of a 12log⁡V\frac{1}{2}\log V21​logV finite-size correction resulting from the scaling of the Clebsch-Gordon coefficients and compute explicitly the O(1)O(1)O(1) contribution reflecting angular-momentum coupling within magnetization blocks. Our analysis uses features of random matrix ensembles and provides a fully analytical treatment for arbitrary spin densities, thereby extending Page type results to non-Abelian sectors and clarifying how SU(2)\mathrm{SU}(2)SU(2) symmetry shapes average entanglement.
Entanglement and magic are fundamental resources that capture the complexity of quantum many-body systems. Non-local magic isolates the irreducible nonstabilizerness intrinsically tied to entanglement. However, evaluating this quantity generally requires a prohibitive minimization over the full Hilbert space, making it computationally inaccessible beyond a few qubits. Here, we overcome this bottleneck by establishing a closed-form expression for the non-local stabilizer entropies of fermionic Gaussian states over local Gaussian unitaries, which we prove at Rényi index α=2α=2α=2 for arbitrary subsystem size, and which can be evaluated in polynomial time directly from the eigenvalues of the reduced Majorana covariance matrix. We apply this framework to characterize fermionic non-local magic across diverse physical regimes: we derive an exact Page-like curve for typical random states, reveal logarithmic scaling at the quantum critical point of the XY model, and establish a quasiparticle picture for magic generation during out-of-equilibrium quantum quenches. Crucially, because our result relies solely on two-point correlation functions, it provides a scalable route for the experimental estimation of fermionic non-local magic in large-scale quantum processors via fermionic shadow tomography.
When describing entanglement in typical midspectrum eigenstates of many-body lattice Hamiltonians, two paradigms have emerged that capture the behavior observed in integrable and nonintegrable systems, Haar-random fermionic Gaussian states and Haar-random pure states, respectively. Remarkably, the former capture the behavior of interacting integrable systems, whose eigenstates are non-Gaussian. We argue that the paradigm that captures both the entanglement properties and the lack of Gaussianity in integrable systems is that of random superpositions of polynomially many Gaussian states. In contrast, eigenstates of nonintegrable systems are consistent with being described by random superpositions of exponentially many Gaussian states. We gain this understanding by comparing analytical and numerical results for the one-body purity, the non-Gaussianity, and the entanglement entropy of the random superpositions and the Hamiltonian eigenstates.
Fermionic Gaussian states are a fundamental tool in many-body physics, faithfully representing non-interacting quantum systems and allowing for efficient numerical simulations. Given a many-body wave function, it is therefore interesting to ask how much it differs from that of a Gaussian state, as quantified by the notion of non-Gaussianity. In this work, we relate measures of non-Gaussianity with the Shannon entropy of the particle-number distribution, coinciding with the particle-number asymmetry for pure states. We derive a lower bound on the relative entropy of non-Gaussianity in terms of the exponential of the Shannon entropy, and study numerically its tightness for large system sizes. Our bound is non-trivial for large values of the asymmetry and relies on the concentration of the particle-number distribution of (mixed) fermionic Gaussian states. Since the Shannon entropy of the particle-number distribution is often efficient to compute or experimentally measure, our results can be viewed as a practical way to lower bound non-Gaussianity, highlighting a non-trivial interplay with particle-number asymmetry.