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.
Abstract The symmetry-resolved Rényi entanglement entropy is the Rényi entanglement entropy of each symmetry sector of a density matrix
ρ ρ . This experimentally relevant quantity is known to have rich theoretical connections to conformal field theory (CFT). For a family of critical free-fermion chains, we present a rigorous lattice-based derivation of its scaling properties using the theory of Toeplitz determinants. We consider a class of critical quantum chains with a microscopic U(1) symmetry; each chain has a low energy description given byNmassless Dirac fermions. For the density matrix,
ρA ρA , of subsystems ofLneighbouring sites we calculate the leading terms in the largeLasymptotic expansion of the symmetry-resolved Rényi entanglement entropies. This follows from a largeLexpansion of the charged moments of
ρA ρA ; we derive
tr(eiαQAρAn) = aeiα⟨QA⟩(σL)−x(1+O(L−μ)) tr(eiαQAρAn)=aeiα⟨QA⟩(σL)-x(1+O(L-μ)) , wherea, xand
μ μ are universal and
σ σ depends only on theNFermi momenta. We show that the exponentxcorresponds to the expectation from CFT analysis. The error term
O(L−μ) O(L-μ) is consistent with but weaker than the field theory prediction
O(L−2μ) O(L-2μ) . However, using further results and conjectures for the relevant Toeplitz determinant, we find excellent agreement with the expansion over CFT operators.
The entanglement spectra for a subsystem in a spin chain fine-tuned to a quantum-critical point contains signatures of the underlying quantum field theory that governs its low-energy properties. For an open chain with given boundary conditions described by a 2D conformal field theory~(CFT), the entanglement spectrum of the left/right half of the system coincides with a boundary CFT spectrum, where one of the boundary conditions arise due to the `entanglement cut'. The latter has been argued to be conformal and has been numerically found to be the `free' boundary condition for Ising, Potts and free boson theories. For these models, the `free' boundary condition for the lattice degree of freedom has a counterpart in the continuum theory. However, this is not true in general. Here, this question is analyzed for the unitary minimal models of 2D CFTs using the density matrix renormalization group technique. The entanglement spectra are computed for blocks of spins in open chains of A-type restricted solid-on-solid models with identical boundary conditions at the ends. The imposed boundary conditions are realized exactly for these lattice models due to their integrable nature. The obtained entanglement spectra are in good agreement with certain boundary CFT spectra. The boundary condition for the entanglement cut is found to be conformal and to coincide with the one with the highest boundary entropy. This identification enables determination of the exponents governing the unusual corrections to the entanglement entropy from the CFT partition functions. These are compared with numerical results.
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.
Symmetry and entanglement are two fundamental concepts in quantum many-body physics. Their interplay is captured by symmetry-resolved entanglement, which decomposes the total entanglement into contributions from different symmetry sectors. Computing symmetry-resolved entanglement in strongly interacting higher-dimensional quantum systems remains challenging. Here, we formulate and implement an estimator-based quantum Monte Carlo (QMC) framework for computing symmetry-resolved Rényi entropies (SRRE) in sign-problem-free interacting lattice systems by measuring disorder (symmetry-twisted) operators in ordinary and replica ensembles and reconstructing SRRE from the corresponding charged moments. We validate the framework in two controlled one-dimensional settings: the transverse-field Ising model (TFIM), for which exact conformal-field-theory predictions are available, and the interacting Heisenberg chain, which tests the U(1) U ( 1 ) symmetry-sector reconstruction and its finite-size behavior. We then apply the method to the two-dimensional TFIM. Within the accessible system sizes and a phenomenological finite-size extrapolation, our data provide numerical evidence consistent with entanglement equipartition at the (2+1) ( 2 + 1 ) D Ising critical point. Our work establishes a practical numerical route to symmetry-resolved entanglement in interacting lattice models and provides a framework for future studies beyond one dimension.
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.