Stabilizer Rényi entropy has emerged as a probe of nonstabilizerness in quantum many-body systems, but its scaling structure beyond critical points remains poorly understood compared with entanglement entropy. Recent field-theory approaches indicate that stabilizer entropy contains universal critical data and boundary-sensitive terms, raising the question of how these structures extend into massive and crossover regimes. We address this problem for a broad class of finite-range spin chains at Rényi index one-half. We derive exact finite-size formulas for both full periodic chains and finite intervals of the infinite chain, making the universal crossover from critical to noncritical behavior analytically accessible. In periodic geometry, the entropy obeys a volume law away from criticality and exhibits a universal finite-size crossover controlled by the competition between system size and correlation length. We also show that the large-scale SRE density develops a cusp across the field-tuned critical line, while the XX endpoint is governed by a distinct scaling regime associated with the saturation point. In the subsystem geometry, the interval entropy separates bulk critical behavior from boundary contributions generated by the way the finite region cuts the infinite chain. The crossover from critical to massive behavior is then encoded in boundary constants and universal functions controlled by the correlation length. Through exact stabilizer-entropy correspondences, the scaling theory extends to internal XY reductions, Finite-range spin chains, and Cluster--Ising representatives. Our results provide an exact lattice benchmark for the emerging QFT description of stabilizer entropy beyond isolated conformal points.
Shannon-Rényi and stabilizer entropies are key diagnostics of structure, "nonstabilizerness," phase transitions, and universality in quantum many-body states. We establish an exact correspondence for quadratic fermions: for any Gaussian eigenstate, the stabilizer Rényi entropy equals the Shannon-Rényi entropy of a number-conserving free-fermion eigenstate on a doubled system, evaluated in the computational basis. Specializing to the transverse-field Ising (TFI) chain, the TFI ground-state stabilizer entropies map to the Shannon-Rényi entropies of the XX-chain ground state of length 2L. Building on this correspondence, together with other exact identities we prove, we derive closed expressions for the stabilizer entropy at indices α=1/2,2, 4 for a broad class of critical closed free-fermion systems. Each of these can be written with respect to the universal functions of the TFI chain. We further derive conformal-field-theory scaling laws for the stabilizer entropy at arbitrary Rényi index under both periodic and open boundary conditions. At α=4, these scaling forms display a discontinuity for both open and periodic boundary conditions.
We study the finite-temperature stabilizer Rényi entropy of the open critical quantum spin chains. At Rényi index one half, this observable probes the distribution of thermal Pauli-string expectation values and can be written as a sum over absolute values of all square minors of a finite-temperature correlation matrix for the transverse-field Ising chain. We show that this exponentially large sum is exactly reducible to a single Pfaffian. The Pfaffian representation reveals a block Toeplitz--Hankel structure and allows us to extract the large-size scaling in several thermal regimes. In the crossover window where the inverse temperature is proportional to the system size, the stabilizer entropy factorizes into a saturated extensive contribution and a universal finite-size scaling function. We find that this scaling function is a level-eight eta quotient, rather than the ordinary free-boundary Majorana thermal factor. The deviation is exponentially hidden at low temperature but controls the high-temperature crossover, where it gives a Cardy-like asymptotic for the Pauli-string expectation-weight spectrum. These results show that finite-temperature stabilizer entropy reveals hidden defect-like conformal data invisible to ordinary thermodynamic probes.
While nonstabilizerness (''magic'') is a key resource for universal quantum computation, its behavior in many-body quantum systems, especially near criticality, remains poorly understood. We develop a spectral transfer-matrix framework for the stabilizer Rényi entropy (SRE) in infinite matrix product states, showing that its spectrum contains universal subleading information. In particular, we identify an SRE correlation length -- distinct from the standard correlation length -- which diverges at continuous phase transitions and governs the spatial response of the SRE to local perturbations. We derive exact SRE expressions for the bond dimension
χ=2 MPS ''skeleton'' of the cluster-Ising model, and we numerically probe its universal scaling along the
Z2 critical lines in the phase diagram. These results demonstrate that nonstabilizerness captures signatures of criticality and local perturbations, providing a new lens on the interplay between computational resources and emergent phenomena in quantum many-body systems.
Understanding universal aspects of many-body systems is one of the central themes in modern physics. Recently, the stabilizer Rényi entropy (SRE) has emerged as a computationally tractable measure of nonstabilizerness, a crucial resource for fault-tolerant universal quantum computation. While numerical results suggested that the SRE in critical states can exhibit universal behavior, its comprehensive theoretical understanding has remained elusive. In this work, we develop a field-theoretical framework for the SRE in a ( 1 + 1 )-dimensional many-body system and elucidate its universal aspects by using boundary conformal field theory. We demonstrate that the SRE is equivalent to a participation entropy in the Bell basis of a doubled Hilbert space, which can be calculated from the partition function of a replicated field theory with the interlayer line defect created by the Bell-state measurements. This identification allows us to characterize the universal contributions to the SRE on the basis of the data of conformal boundary conditions imposed on the replicated theory. We find that the SRE of the entire system contains a universal size-independent term determined by the noninteger ground-state degeneracy known as the g factor. In contrast, we show that the mutual SRE exhibits the logarithmic scaling with a universal coefficient given by the scaling dimension of a boundary-condition-changing operator, which elucidates the origin of universality previously observed in numerical results. As a concrete demonstration, we present a detailed analysis of the Ising criticality, where we analytically derive the universal quantities at arbitrary Rényi indices and numerically validate them with high accuracy by employing tensor-network methods. These results establish a field-theoretical approach to understanding the universal features of nonstabilizerness in quantum many-body systems.
Quantum magic---the resource that separates universal quantum computation from efficiently simulable Clifford circuits---has emerged as a diagnostic of many-body quantum states, yet the stabilizer Rényi entropy (SRE) that quantifies it remains largely inaccessible in interacting fermion systems. Estimating the global SRE requires a sum over exponentially many Majorana strings, which in determinant quantum Monte Carlo can be sampled without a sign problem only for the interacting density matrix obtained after averaging the auxiliary fields: sampling strings and fields simultaneously incurs a sign problem even for sign-problem-free models, so that the known sign-free alternative is a nested Monte Carlo that does not scale. We propose instead the two-point SRE---a practical stabilizer-Rényi correlator built from rank-2 SREs of one- and two-site reduced density matrices, which are fixed exactly by single-particle Green functions and density correlations already measured in standard simulations---and show that it obeys universal finite-size scaling at fermionic quantum critical points. In the one-dimensional half-filled spinless
t-
V chain, the correlator distinguishes algebraic and exponential decay regimes and tracks the inverse-logarithmic finite-size drift characteristic of the Berezinskii--Kosterlitz--Thouless transition. On the honeycomb lattice, sign-problem-free quantum Monte Carlo yields Gross--Neveu--Ising scaling with finite anomalous dimension at zero temperature and two-dimensional Ising collapse at the thermal transition. Our results suggest stabilizer-Rényi microscopy as a spatially resolved, quantitatively universal probe of critical correlations in interacting fermionic matter, on par with conventional order-parameter correlators and accessible to quantum-simulator measurements.
Stabilizer Rényi entropy quantifies the nonstabilizerness of a quantum state through Rényi moments of its Pauli expectation-value distribution. Its standard form sums over Pauli-string degrees and retains only the total Rényi weight. We introduce a fugacity-resolved partition function for the critical transverse-field Ising chain that resolves this degree in the balanced Majorana representation and generates its full counting statistics. This Ising problem extends beyond a single model: exact decimation identities and the correspondence with the half-filled
XX chain establish it as a common finite-size building block for stabilizer and computational-basis Shannon--Rényi entropies. We map the fugacity-resolved all-minors sum, for every positive Rényi index, to a checkerboard-weighted discrete Selberg ensemble on a half-filled doubled root lattice. For positive integer indices, it reduces to a finite aliased Dyson constant term. At
α=21,1,2, determinant and Pfaffian compressions yield product formulas. At
α=4, the generic-fugacity problem admits an exact inverse Jack--Kostka representation, while at unit fugacity a complementary middle-minor identity relates it to the square of the
α=2 result. These generating functions determine the balanced-degree statistics. Complement symmetry makes the distribution symmetric about
k=L/2, with an index-dependent width. It is exactly binomial at
α=1, has variance proportional to
L at
α=21, and develops an
LlogL enhancement at
α=2. After variance rescaling, the centered distributions converge to Gaussian limits at all three indices. Beyond these solvable cases, finite-size numerics reveal an evolution from a central peak to a symmetric bimodal profile and eventually to endpoint dominance as the Rényi index increases.