Fugacity-Resolved Stabilizer Entropy in Critical Quantum Chains: Discrete Selberg Sums and Exactly Solvable Rényi Indices

AI summary

Overview

Research area: Quantum information and quantum many-body physics, specifically the resource theory of nonstabilizerness ("magic") in critical free-fermion chains, combined with algebraic combinatorics (Vandermonde/Selberg sums, Dyson constant terms, Pfaffians, Jack–Kostka theory).

Technical level: Advanced. The paper relies on Majorana fermionization, Cauchy determinant identities, de Bruijn–type Pfaffian formulas, root-of-unity Fourier aliasing, and symmetric-function theory. The physics motivation is accessible; the derivations are not.

Scope (one sentence): The paper introduces a fugacity-resolved generating function for the stabilizer Rényi entropy of the critical transverse-field Ising chain and solves it exactly at several Rényi indices, obtaining explicit product formulas and the full counting statistics of the balanced Majorana degree.

What This Paper Is About

The stabilizer Rényi entropy measures how far a quantum state is from the efficiently simulable "stabilizer" sector by taking Rényi moments of the state's Pauli expectation values. For Gaussian (free-fermion) states, each Pauli expectation value is a determinant or Pfaffian of a submatrix, but the entropy requires summing over all "balanced" minors, whose number grows as a central binomial coefficient — so even an exactly solvable chain does not automatically give a tractable entropy. This paper asks which Rényi indices allow that all-minors sum to be compressed exactly at finite size, and it answers by keeping track of how many Pauli strings contribute at each balanced Majorana degree through an extra fugacity variable.

Key Contributions

1. A fugacity-resolved partition function for the critical TFI chain. The authors define a polynomial in a fugacity variable that weights each balanced Majorana sector of degree k, so that the standard stabilizer Rényi moment is recovered at unit fugacity. This upgrades the single number of the usual entropy into a full probability-generating function.

2. An exact term-by-term mapping to a discrete Selberg (Vandermonde) gas. Balanced minors of the half-shift correlation matrix are shown to be in bijection with half-filled configurations on a doubled root-of-unity lattice of 2L sites, with the fugacity acting as a checkerboard chemical potential distinguishing the two interlaced sublattices. This representation holds for every real α > 0.

3. A finite aliased Dyson–Morris constant-term representation for positive integer α. The Vandermonde weight has finite Fourier bandwidth, so the root-of-unity measure truncates to a finite alias kernel inside constant-term extraction, giving an exact representation built from the Dyson Laurent polynomial plus a one-body checkerboard kernel.

4. Exact product formulas at α = 1/2, 1, 2 and an exact (non-product) representation at α = 4, together with the resulting degree-resolved full counting statistics and its Gaussian limits.

Main Findings

• Fugacity generates the full statistics. Because the fugacity-resolved partition function is the generating function of the balanced Majorana degree, the exact formulas determine the entire degree distribution, normalized as P_{α,L}(k).

• Complement symmetry. The doubled-root representation gives Z_{α,L}(u) = u^L Z_{α,L}(u^{-1}) and W_{α,L}(k) = W_{α,L}(L-k), so the degree distribution is symmetric about k = L/2. The endpoint sectors, containing only one complete sublattice, satisfy W_{α,L}(0) = W_{α,L}(L) = 1.

• α = 1/2, closed product. The fugacity-resolved result is a product over j = 1, …, L/2 of [1 + u² + 2u sec(π(2j−1)/(2L))]. The unit-fugacity value reduces to 2^{L/2} times a product of [1 + sec(π(2j−1)/(2L))], attributed to earlier work. Derivation routes: an ordered-root de Bruijn Pfaffian with checkerboard 2×2 block factorization, plus an independent direct 4L×4L balanced-all-minors Pfaffian representing the complete fugacity polynomial.

• α = 1, exactly binomial. Z_{1,L}(u) = (1+u)^L, and Z_{1,L}(1) = 2^L, obtained via a Cauchy–Binet determinant and root orthogonality for the discrete β = 2 ensemble. The degree distribution is therefore exactly binomial.

• α = 2, closed product with L log L broadening. Z_{2,L}(u) is a product over r = 1, …, L/2 of [(1+u)² − 4u((2r−1)/L)²], with unit-fugacity value (2L)!/(L! L^L). The derivation uses a multiplicity-two confluent de Bruijn Pfaffian for the discrete β = 4 ensemble. The variance develops an L log L enhancement.

• α = 4, structurally different. At generic fugacity the problem admits an exact neutral seven-charge shifted-Dyson and rectangular inverse Jack–Kostka representation, which is exact but neither a closed product nor an ordinary-Pfaffian compression. At unit fugacity a complementary middle-minor identity collapses it to a product: 2^L/L^{2L} [(2L−1)!!]², equal to 2^{−L} Z_{2,L}(1)².

• α = 0 as a formal counting endpoint. The table lists the formal counting case with Σ_k C(L,k)² u^k and unit-fugacity value C(2L,L).

• Fluctuation scales differ but limits are Gaussian. At α = 1/2 the variance is proportional to L; at α = 2 it shows the L log L enhancement; at α = 1 the distribution is exactly binomial. Despite these different scales, the centered distributions rescaled by their standard deviations converge to Gaussian limits at all three product-solvable indices.

• Beyond the solvable indices, profile reorganization. Finite-size numerics show that increasing the Rényi tilt evolves the degree-resolved Pauli weight from a central peak to a symmetric bimodal profile and eventually to endpoint dominance. The paper does not report the system sizes used for these numerics in the excerpted content.

• General structural hierarchy. Perfect-square indices form a confluent hyperpfaffian hierarchy, while positive integer indices admit complementary-minor and shifted-Dyson organizations. The authors explicitly distinguish exact structural representations from polynomial-size determinant and Pfaffian evaluations.

• The TFI result is a common building block. Exact reductions are collected for: the TFI–XX relation M_α^TFI(L) = H_α^XX(2L) − L log 2; stabilizer decimation M_α^{(z^n+1)}(L) = n M_α^TFI(L/n) for n | L and M_α^{(z^m+z^{-m})}(L) = 2m M_α^TFI(L/(2m)) for 2m | L; Shannon factorization H_α^{(z^m+z^{-m})}(L) = m H_α^XX(L/m); the combined relation H_α^{(z^m+z^{-m})}(2mL) = m M_α^TFI(L) + mL log 2; and the Haldane–Shastry reduction H_α^HS(2L) = L log 2 + [(1−2α)M_{2α}^TFI(L) + α M_2^TFI(L)]/(1−α).

• New probability-level factorization proved. For L = mℓ, the full computational-basis distribution of the range-m XX chain factorizes into m squeezed XX chains of length ℓ, giving an exact Shannon–Rényi decimation for the whole family.

• Haldane–Shastry enters as an order-two escort. Since p_HS(U) ∝ |Δ(U)|⁴ and p_XX(U) ∝ |Δ(U)|², the Haldane–Shastry distribution is the order-two escort of the half-filled XX distribution. The authors note this applies to computational-basis wave-function probabilities and does not imply a corresponding stabilizer Rényi identity.

• Ordinary Dyson is not enough. The ordinary identity CT D_α = (αL)!/(α!)^L applies to the unrefined problem only when the visible-alias kernel is constant — which happens for α = 1 and α = 2 at unit fugacity, whereas for α = 3 and α = 4 the unit-fugacity kernel is 1 + x^{2L} + x^{−2L}.

Methodology in Plain English

The authors start from the standard definition of stabilizer Rényi entropy as a sum over Pauli strings of powers of expectation values. Using a Majorana representation, each Pauli string is labeled by two equal-size index sets S and T, and its expectation value is a minor of an L×L matrix whose entries are cosecants with a half-integer shift — the "half-shift matrix."

They then attach a variable u to the size k of those index sets. This turns one number into a polynomial whose coefficients record how much weight sits at each balanced degree, and whose derivatives at u = 1 give moments of the degree distribution.

The key technical move is to recognize that the half-integer shift makes the matrix a Cauchy matrix. The Cauchy determinant formula converts each minor into a Vandermonde product over a set of roots of unity living on a lattice of twice the size, half of whose sites are occupied. The fugacity becomes a staggered weight favoring one of the two interlaced sublattices. The all-minors sum thus becomes a discrete version of a Selberg integral — a gas of particles with logarithmic repulsion at inverse temperature set by α.

For integer α, the Vandermonde weight is a Laurent polynomial with bounded exponents, so modes of the root-of-unity Fourier comb beyond a finite range cannot be cancelled and drop out exactly. What remains is a constant-term problem: extract the zero-mode coefficient of the Dyson product times a small one-body kernel. At the special indices the constant term can be evaluated in closed form, either through determinant identities and root orthogonality (α = 1), through de Bruijn-type Pfaffians for the β = 2 and β = 4 ensembles (α = 1/2, 2), or through symmetric-function expansions in inverse Jack–Kostka form (α = 4). Standard probability tools then convert the generating polynomial into variance, limiting-shape, and Gaussianity statements, supplemented by finite-size numerical evaluation away from the solvable points.

Why This Matters

The paper supplies a template for turning a hard "sum over exponentially many minors" problem into a structured statistical-mechanical sum that can sometimes be done exactly. It clarifies that integrability of the chain is not the bottleneck; the bottleneck is the algebraic mechanism available at a given Rényi index. It also shows that a single free-fermion building block underlies a wide family of stabilizer and Shannon–Rényi observables, which organizes a fragmented literature into one reduction table.

• Quantum error correction and classical simulation: the paper frames stabilizer states and Clifford operations as the efficiently tractable sector underlying quantum error correction and classical simulation, so quantitative measures of departure from that sector are directly relevant to assessing simulability.
• Benchmarking magic measures: exact finite-size results provide ground truth against which approximate, tensor-network, Pauli-sampling, and determinantal-sampling methods for nonstabilizerness can be validated.
• Critical-point diagnostics: earlier work cited in the paper argues magic is enhanced near quantum critical points and can carry universal data, so exact lattice results help separate universal from finite-size contributions.
• Sampling-based estimation protocols: the degree-resolved distribution identifies which Pauli-string sectors dominate, which is the kind of structural information that sampling schemes need.

Industry relevance: The paper reports no experimental, hardware, or industrial demonstration, and no concrete algorithmic speedup claims for practical tasks. Its relevance is upstream and methodological — it strengthens the theory base for quantifying nonstabilizerness in fermionic systems, which is the same quantity used to reason about classical simulability of quantum circuits. Any industrial connection runs through that broader program, not through results stated in this paper.

Future Directions

• Extending the constant-term machinery beyond the solvable indices. The paper establishes exact representations for all positive integer α but closed products only at α = 1/2, 1, 2 (and at unit fugacity for α = 4). Whether other indices admit comparably compact evaluations, or whether the hyperpfaffian and shifted-Dyson hierarchies can be pushed to explicit evaluations, is left open.

• Turning structural representations into efficient algorithms. The authors explicitly distinguish exact structural representations from polynomial-size determinant and Pfaffian evaluations. A natural next step is to determine which of the generic-fugacity constructions (the rectangular inverse Jack–Kostka form, the seven-charge shifted-Dyson form) can be evaluated at a cost polynomial in L.

• Sharpening the non-solvable fluctuation profiles. The central-peak to bimodal to endpoint-dominated evolution is reported from finite-size numerics. Whether these regimes have sharp transitions in α, and what sets their scaling with L, is not settled.

• Transfers to other models and to continuum descriptions. Since the TFI result is a common building block for the range-m XX chain, the Haldane–Shastry chain, and paired free-fermion chains, the fugacity-resolved formulas can be carried over through the reduction table; connecting the lattice full counting statistics to conformal-boundary and defect descriptions is a further open direction.

Target Audience

Quantum information theorists and condensed-matter physicists working on nonstabilizerness, magic, and entanglement in free-fermion and critical systems; mathematical physicists and combinatorialists interested in Selberg-type sums, Dyson constant terms, de Bruijn Pfaffians, and Jack–Kostka theory; and researchers developing classical-simulation or Pauli-sampling algorithms who need exact benchmarks for Gaussian states. The paper assumes comfort with Majorana fermionization and symmetric-function techniques, so it is best suited to readers at the graduate-student level and above in theoretical quantum many-body physics.

Authors’ abstract

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 \(α=\tfrac12,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 \(α=\tfrac12\), and develops an \(L\log L\) 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.