We study the non-stabilizerness (quantum magic) content of the Hubbard dimer, an analytically solvable, yet completely non-trivial, model of strongly correlated fermions. We can access zero- and finite-temperature properties as well as the time evolution in a quantum quench protocol. We evaluate local and nonlocal non-stabilizerness using both the robustness of magic and the stabilizer Renyi entropy, demonstrating how the latter often fails in detecting the mixed stabilizer states that are typically found in this kind of systems. Finally, we compare the non-stabilizerness with other genuine resources of quantum-state complexity, i.e., the fermionic non-Gaussianity and the superselected two-site entanglement. Our findings corroborate the notion of non-stabilizerness as a fundamentally different quantum resource, able to give profound insights that are missed by more traditional information-theoretic quantities.
We study the non-stabilizerness (quantum magic) content of the Hubbard dimer, an analytically solvable, yet completely non-trivial, model of strongly correlated fermions. We can access zero- and finite-temperature properties as well as the time evolution in a quantum quench protocol. We evaluate local and nonlocal non-stabilizerness using both the robustness of magic and the stabilizer Renyi entropy, demonstrating how the latter often fails in detecting the mixed stabilizer states that are typically found in this kind of systems. Finally, we compare the non-stabilizerness with other genuine resources of quantum-state complexity, i.e., the fermionic non-Gaussianity and the superselected two-site entanglement. Our findings corroborate the notion of non-stabilizerness as a fundamentally different quantum resource, able to give profound insights that are missed by more traditional information-theoretic quantities.
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.
We investigate the structure of quantum magic in interacting disordered fermionic systems, quantifying non-stabilizerness via the fermionic stabilizer Rényi entropy (SRE). To resolve the distribution of magic across different scales, we introduce a multipartite non-local magic functional, constructed from an inclusion-exclusion combination of subsystem contributions. This measure serves as a fine-grained diagnostic, isolating genuinely global contributions and revealing nontrivial interactions between local and collective supports of magic. We illustrate the measure on paradigmatic multipartite states and apply these diagnostics to the Sachdev-Ye-Kitaev model and its variants. Crucially, for thermal/typical ensembles, we observe a marked disparity between Thermal Pure Quantum (TPQ) states and the thermal density matrix. This reveals a concealed complexity: the immense computational hardness characterizing the unitary evolution is encoded in the specific microstructure of the black hole microstates, while being washed out in the coarse-grained thermodynamic description. Furthermore, in N=2 supersymmetric SYK, we show that while fortuitous BPS states exhibit intermediate stabilizer complexity, the multipartite measure unveils a rich, sector-dependent pattern of global correlations, distinguishing them from generic chaotic states.
We introduce an efficient method to quantify nonstabilizerness in fermionic Gaussian states, overcoming the long-standing challenge posed by their extensive entanglement. Using a perfect sampling scheme based on an underlying determinantal point process, we compute the Stabilizer Renyi Entropies (SREs) for systems with hundreds of qubits. Benchmarking on random Gaussian states with and without particle conservation, we reveal an extensive leading behavior equal to that of Haar random states, with logarithmic subleading corrections. We support these findings with analytical calculations for a set of related quantities, the participation entropies in the computational (or Fock) basis, for which we derive an exact formula. We also investigate the time evolution of non-stabilizerness in a random unitary circuit with Gaussian gates, observing that it converges in a time that scales logarithmically with the system size. Applying the sampling algorithm to a two-dimensional free-fermionic topological model, we uncover a sharp transition in non-stabilizerness at the phase boundaries, highlighting the power of our approach in exploring different phases of quantum many-body systems, even in higher dimensions.
Characterizing a quantum state through the lens of quantum resources provides an information-theoretic perspective on many-body systems. While quantum entanglement serves as the paradigmatic example of a quantum resource, recent studies have shown that quantum magic, a resource for universal quantum computation, captures aspects of many-body states complementary to those described by entanglement. For instance, in spin systems, conformal field theory (CFT) analysis of the stabilizer Rényi entropy has revealed universal features of nonstabilizerness qualitatively distinct from entanglement. In bosonic and fermionic systems, however, a comparable formulation for their computational resource, non-Gaussianity, has yet to be established. In this work, we introduce a unified measure, the magic Rényi entropy (MRE), to quantify computational resources in spins, bosons, and fermions on an equal footing. We show that the MRE is a resource monotone under stabilizer and Gaussian protocols involving measurements and feedforward operations. The MRE reveals common universal aspects of nonstabilizerness and non-Gaussianity in critical many-body states. In particular, our CFT analysis shows that the universal contribution to the MRE appears as the size-independent term determined by the Affleck-Ludwig boundary entropy. We find that non-Gaussianity can continuously renormalize this universal contribution or drive a boundary transition through bulk-induced boundary renormalization-group flows. As a concrete example, we present a CFT analysis of non-Gaussianity in interacting spinless fermions described by the Tomonaga-Luttinger liquid, showing boundary transitions at the Luttinger parameters K=1/3 and K=3. Our field-theoretical predictions are confirmed by numerical calculations. These results provide a unified field-theoretical understanding of many-body magic across spins, bosons, and fermions.
Understanding how collective phenomena such as pairing and deformation emerge in interacting fermionic systems, and how the correlations underlying them shape the resources required for their simulation, is a central challenge in many-body physics. We investigate various aspects of quantum complexity in a fermionic system exhibiting two main forms of collectivity, specifically deformation and pairing, through entanglement, non-stabilizerness, and fermionic non-Gaussianity. In the deformed phase, high entanglement coexists with low non-stabilizerness, while the superfluid pairing phase, by contrast, exhibits substantial multipartite entanglement and non-stabilizerness, together with maximal fermionic anti-flatness. In a second step we examine how to describe the phase diagram of the system in a classically scalable and symmetry-preserving way, by combining stabilizer-state expansion and variational basis optimization. A physics-guided decomposition of the Hamiltonian identifies the relevant stabilizer states, which can be constructed from graph states. A single stabilizer state captures the deformed phase to a large extent, and, although the number of stabilizer states required to describe the pairing phase grows exponentially with system size, a symmetry-preserving basis rotation yields a high-fidelity approximation using far fewer states. This rotated state can then serve as a practical starting point for quantum refinement.
Understanding the computational complexity of quantum states is a central challenge in quantum many-body physics. In qubit systems, fermionic Gaussian states can be efficiently simulated on classical computers and hence can be employed as a natural baseline for evaluating quantum complexity. In this work, we develop a framework for quantifying fermionic magic resources, also referred to as fermionic non-Gaussianity, which constitutes an essential resource for universal quantum computation. We leverage the algebraic structure of the fermionic commutant to define the fermionic antiflatness (FAF)—an efficiently computable and experimentally accessible measure of non-Gaussianity, with a clear physical interpretation in terms of Majorana fermion correlation functions. Studying systems in equilibrium, we show that FAF detects phase transitions, reveals universal features of critical points, and uncovers special solvable points in many-body systems. Extending the analysis to out-of-equilibrium settings, we demonstrate that fermionic magic resources become more abundant in highly excited eigenstates of many-body systems. We further investigate the growth and saturation of FAF under ergodic many-body dynamics, highlighting the roles of conservation laws and locality in constraining the increase of non-Gaussianity during unitary evolution. This work provides a framework for probing quantum many-body complexity from the perspective of fermionic Gaussian states and opens up directions for investigating fermionic magic resources in many-body systems. Our results establish fermionic non-Gaussianity, alongside entanglement and nonstabilizerness, as a resource relevant not only to foundational studies but also to experimental platforms aiming to achieve quantum advantage.