---
title: Finite-Temperature Stabilizer Rényi Entropy
url: https://www.emergentmind.com/topics/finite-temperature-stabilizer-renyi-entropy
type: topic
---

# Finite-Temperature Stabilizer Rényi Entropy

Searching arXiv for recent papers on finite-temperature stabilizer Rényi entropy and related formulations.
Finite-temperature stabilizer Rényi entropy extends stabilizer Rényi entropy from pure states to Gibbs states \(\rho(\beta)=e^{-\beta H}/\mathrm{Tr}\,e^{-\beta H}\) by combining a participation-type Rényi entropy of Pauli-string or Majorana-string expectation values with the ordinary Rényi entropy of the thermal density matrix. Across recent work, the object appears in several closely related notational conventions: for spin systems, \(S_n^{(\mathrm{stab})}(\beta)\equiv M_n(\rho)-S_n(\rho)\) for Rényi index \(n\) [2405.19577]; for Gibbs states at \(\alpha=2\), \(\widetilde M_2(\rho)=M_2(\rho)-S_2(\rho)\) [2501.12146]; for coupled SYK, the second-Rényi SRE is defined from the Majorana-string expansion of \(\rho_\beta\) [2509.17417]; and for the open critical transverse-field Ising chain, the mixed-state stabilizer Rényi-\(\tfrac12\) entropy is written \(S_{\rm stab}^{(1/2)}(T)=M_{1/2,L}(\beta)\) [2606.08606]. In all cases, the finite-temperature quantity is designed to isolate nontrivial string-structure beyond ordinary thermodynamics, although recent papers draw sharply different conclusions about its status as a mixed-state magic measure.

## 1. Definitions and operator-string formulations

For a general \(N\)-qubit density matrix \(\rho\), one definition used in spin systems is
\[
M_{\alpha}(\rho)=\frac{1}{1-\alpha}\ln\!\biggl[\frac{1}{2^N}\sum_{P\in\mathcal P_N}\bigl|\mathrm{Tr}(\rho P)\bigr|^{2\alpha}\biggr],
\]
with \(\mathcal P_N=\{I,X,Y,Z\}^{\otimes N}\!/\{\pm1,\pm i\}\), and for Gibbs states one extends the construction by subtracting the usual Rényi entropy. In particular, for \(\alpha=2\),
\[
\widetilde M_2(\rho)=M_2(\rho)-S_2(\rho),\qquad S_2(\rho)=-\ln \mathrm{Tr}\,\rho^2.
\]
Introducing \(Q=\sum_{P\in\mathcal P_N}\bigl|\mathrm{Tr}[e^{-\beta H}P]\bigr|^4\) and \(Z_2=\mathrm{Tr}[e^{-2\beta H}]\), one obtains
\[
\widetilde M_2(\rho)=-\ln\!\Bigl[\tfrac{1}{2^N}\,\tfrac{Q}{Z^4}\Bigr]+\ln\!\Bigl[\tfrac{Z_2}{Z^2}\Bigr].
\]
This representation separates ordinary partition functions from a generalized four-replica object \(Q\) [2501.12146].

In the sign-problem-free spin-Hamiltonian framework, the \(n\)-th finite-temperature stabilizer Rényi entropy is written
\[
S_n^{(\mathrm{stab})}(\beta)\equiv M_n(\rho)-S_n(\rho),
\]
where
\[
M_n(\rho)=\frac{1}{1-n}\ln\!\Bigl[\sum_{\sigma\in\mathcal P_N}\frac{1}{2^N}\mathrm{Tr}(\rho\sigma)^{2n}\Bigr],\qquad
S_n(\rho)=\frac{1}{1-n}\ln \mathrm{Tr}\,\rho^n.
\]
Within that convention, \(S_n^{(\mathrm{stab})}(\beta)\ge 0\), vanishes if and only if \(\rho\) is a stabilizer (Clifford) state, and is additive under tensor-product states [2405.19577].

For the coupled SYK model, Zhang et al. specialize to the second-Rényi case and define the participation Rényi entropy of the Majorana-string expansion of the thermal density matrix by
\[
M_2(\rho_\beta)\equiv-\ln\!\Bigl[2^{-N}\sum_{\mathcal v}c_{\mathcal v}(\rho_\beta)^4\Bigr],\qquad
\rho_\beta=2^{-N}\sum_{\mathcal v}c_{\mathcal v}(\rho_\beta)\Psi_{\mathcal v},
\]
with \(\Psi_{\mathcal v}\) an orthonormal operator basis of Majorana strings. The stabilizer Rényi entropy is then
\[
S_{\mathrm{SRE}}(\beta)\equiv M_2(\rho_\beta)-S_2(\rho_\beta).
\]
Equivalently, writing \(S_{\mathrm{SRE}}=-\ln Z_{\mathrm{SRE}}\), one has the exact identity
\[
M_2(\rho_\beta)=S_{\mathrm{SRE}}(\beta)-4S_\beta+N\ln 2,
\]
with \(S_\beta=-\ln Z_\beta\) [2509.17417].

A distinct but related formulation appears at Rényi index \(\tfrac12\) for the open critical transverse-field Ising chain:
\[
S_{\rm stab}^{(1/2)}(T)\equiv M_{1/2,L}(\beta)
=2\ln\!\Bigl[2^{-L}\sum_{P\in\mathcal P_L}\bigl|\mathrm{Tr}(\rho_\beta P)\bigr|\Bigr]-2\ln \mathrm{Tr}\,\rho_\beta^2.
\]
Here the numerator is the first moment of thermal Pauli-string expectation values, and Wick’s theorem reduces it to a sum over absolute values of all square minors of the finite-temperature correlation matrix \(G(\beta)\) [2606.08606].

## 2. Replica and path-integral constructions

A common structural feature is replica enlargement. In the generic spin-system construction, one introduces \(2n\) replicas and uses the identity
\[
\sum_{\sigma\in\mathcal P_N}\sigma^{\otimes 2n}=\bigotimes_{s=1}^N \mathcal T,
\]
where \(\mathcal T\) is a site-local connection tensor of rank \(4n\). For any subset \(B\subseteq\{1,\dots,N\}\), the interpolating operators
\[
T(B)\equiv (1/2^{|B|})\bigotimes_{s=1}^N \bigl[B_s\mathcal T+(1-B_s)I^{\otimes 2n}\bigr]
\]
define partition functions
\[
\mathcal Z_B\equiv \mathrm{Tr}\bigl[(e^{-\beta H/2})^{\otimes 2n}T(B)(e^{-\beta H/2})^{\otimes 2n}\bigr].
\]
Then \(\mathcal Z_{[N]}\) gives the string moment and \(\mathcal Z_\emptyset=(\mathrm{Tr}\,e^{-\beta H})^{2n}\), so that
\[
M_n(\beta)=\frac{1}{1-n}\ln\!\bigl[\mathcal Z_{[N]}/\mathcal Z_\emptyset\bigr].
\]
An interpolation \(\mathcal Z(\lambda)=\sum_{B\subseteq[N]}\lambda^{|B|}(1-\lambda)^{N-|B|}\mathcal Z_B\) converts the problem into a free-energy difference, and Jarzynski’s equality yields \(e^{-(1-n)M_n(\beta)}=\langle e^{-W}\rangle_{\rm paths}\) [2405.19577].

For \(\alpha=2\) in spin models, Ding, Wang, and Yan formulate the four-replica quantity \(Q\) in stochastic-series-expansion form and identify the source of sign cancellations in direct Pauli-string sampling. They replace \(Y\) by \(\tilde Y=\begin{pmatrix}0&-1\\1&0\end{pmatrix}\) and restrict the sum to a reduced configuration space in which, at each site and at the imaginary-time boundary, an even number of replicas carry a spin-up. This yields a reduced Pauli-string set \(\tilde{\mathcal P}_N\) and the exact rewriting
\[
Q=\sum_{\tilde P\in\tilde{\mathcal P}_N}\sum_{\mathrm{valid}\;S_\Lambda^k,\alpha_0^k}\prod_{k=1}^4
W\bigl(\tilde P;S_\Lambda^k,\ket{\alpha_0^k}\bigr),
\]
with all weights strictly nonnegative [2501.12146].

In the coupled SYK model, the finite-temperature SRE is represented by four replicas \(\alpha=1,\dots,4\) and a quartic insertion at imaginary time \(\tau=0\),
\[
\prod_m\bigl[1+4\psi_m^{(1)}\psi_m^{(2)}\psi_m^{(3)}\psi_m^{(4)}\bigr].
\]
Equivalently, one decouples the insertion by auxiliary Ising spins \(\sigma_m=\pm1\) and site-local SWAP operators
\[
S_m^{\alpha\beta}[\sigma_m]=\exp\!\bigl[(\pi/2)\sigma_m\,\psi_m^{(\alpha)}\psi_m^{(\beta)}\bigr],
\]
leading to
\[
Z_{\mathrm{SRE}}
=\sum_{\{\sigma_m\}}\mathrm{Tr}\!\Bigl[(e^{-\beta H})^{\otimes 4}\prod_m S_m^{12}[\sigma_m]\,S_m^{34}[\sigma_m]\Bigr].
\]
After disorder averaging and bilocal Hubbard-Stratonovich fields, the large-\(N\) limit becomes a self-consistent saddle-point problem [2509.17417].

The open-chain Ising analysis replaces stochastic sampling or large-\(N\) saddle methods by an exact algebraic reduction. The exponentially large sum of absolute values of all square minors of the correlation matrix is exactly reducible to a single Pfaffian, which in the open-chain basis can be rewritten as a finite-size block Toeplitz–Hankel Pfaffian. This exact reformulation places the stabilizer quantity within a class of Pfaffian determinants amenable to asymptotic methods [2606.08606].

## 3. Large-\(N\) finite-temperature SRE in the coupled SYK model

The coupled SYK analysis establishes a general framework for analyzing the SRE in solvable SYK models in the large-\(N\) limit, enabling the application of the saddle-point approximation [2509.17417]. After disorder averaging in the Maldacena–Qi model and introducing bilocal fields \(G_{ss'}^{(\alpha\beta)}(\tau,\tau')\) and \(\Sigma_{ss'}^{(\alpha\beta)}(\tau,\tau')\), with \(s,s'=L/R\) labeling the two SYK clusters, the saddle equations are
\[
\Sigma_{ss'}^{(\alpha\beta)}(\tau,\tau')=J^2\,G_{ss'}^{(\alpha\beta)}(\tau,\tau')^3,
\]
together with
\[
G_{ss'}^{(\alpha\beta)}
=
\frac{\sum_{\sigma_L,\sigma_R} g_{ss',\sigma_L\sigma_R}^{(\alpha\beta)}
\,[\det g_{\sigma_L\sigma_R}^{(\alpha\beta)}]^{-1/2}}
{\sum_{\sigma_L,\sigma_R}[\det g_{\sigma_L\sigma_R}^{(\alpha\beta)}]^{-1/2}},
\]
where the fixed-boundary Green’s functions \(g\) are defined by inverting
\[
(\partial_\tau-\mu(\sigma^y)-\Sigma)g=1
\]
subject to twisted anti-periodic boundary conditions \(\sigma_s=\pm1\) on each cluster. The saddle action is
\[
S_{\mathrm{SRE}}[G,\Sigma]
=
-N\ln\!\Bigl[\sum_{\sigma_L,\sigma_R}\det[g_{\sigma_L\sigma_R}^{(\alpha\beta)}]^{-1/2}\Bigr]
+\frac{3N}{8}\sum_{\alpha\beta,ss'}\int_0^\beta d\tau\,d\tau'\,
G_{ss'}^{(\alpha\beta)}(\tau,\tau')\Sigma_{ss'}^{(\alpha\beta)}(\tau,\tau').
\]
Together with the usual SYK free energy \(S_\beta[N,G_0,\Sigma_0]\), this yields \(M_2\) and hence \(S_{\mathrm{SRE}}(\beta)=M_2(\beta)+4S_\beta-N\ln 2\) [2509.17417].

The analysis identifies controlled limits. At high temperature, \(\beta J\ll1\) and \(\mu\ll J\), one finds perturbatively
\[
M_2(\beta)\simeq N\bigl[\ln 2-c(\mu\beta)^2+\cdots\bigr],
\]
so the SRE starts at zero at \(\beta\to0\) and grows quadratically. At low temperature, in the decoupled limit \(\mu=0\), one exactly shows \(M_2(\beta)=N\ln 2\) for all \(\beta\). For \(\mu>0\) small and \(\beta\to\infty\),
\[
M_2(\infty)=N\bigl[2\ln 2-4s_0\bigr]<N\ln 2,
\]
where \(s_0\approx0.2324\) is the SYK zero-\(T\) entropy density, so \(S_{\mathrm{SRE}}(\beta)=M_2(\beta)-S_2(\beta)\) saturates to a constant less than \(N(\ln 2-2s_0)\) [2509.17417].

For \(\mu/J=0.1\), numerical solution of the saddle-point equations yields three first-order transitions in \(M_2(\beta)\), and hence in SRE, as temperature is tuned. The first is \(\beta^*_{HP/2}\approx12\), described as a “half” Hawking–Page transition and inherited from a discontinuity in \(S_{2\beta}\) at \(\beta=2\beta_{HP}\). The second is \(\beta^*_{\mathrm{SRE}}\approx19\), an “intrinsic” SRE transition with no counterpart in the thermal free energy, arising from a change in the dominant \(\sigma\)-boundary-twist sector. The third is \(\beta^*_{HP}\approx27\), the usual wormhole\(\leftrightarrow\)black-hole transition in \(S_\beta\), which also feeds into \(M_2\) through the \(-4S_\beta\) term. In each case, the order parameter is the dominant connectivity of the four replicas, and the critical point \(\beta^*_{\mathrm{SRE}}\) is characterized by two coexisting large-\(N\) saddles of equal action \(S_{\mathrm{SRE}}\) but different replica-connectivity pattern [2509.17417].

A central point is that \(\beta^*_{\mathrm{SRE}}\) is invisible to thermodynamics: neither the thermal free energy \(S_\beta\) nor any standard correlator shows a singularity there, yet \(M_2(\beta)\) and thus \(S_{\mathrm{SRE}}(\beta)\) jump. The paper therefore states that SRE serves as a genuine order parameter for a new class of finite-\(T\) transitions which cannot be detected by conventional thermodynamic observables, and describes the change of SWAP-twist boundary conditions as a switch from a “disconnected” to a “connected” phase, analogous to replica-wormhole saddles in gravitational path integrals [2509.17417].

## 4. Exact finite-size scaling and hidden boundary data in the open critical Ising chain

For the open critical transverse-field Ising chain, the mixed-state stabilizer Rényi-\(\tfrac12\) entropy admits an exact finite-size treatment [2606.08606]. Writing
\[
\sum_{P\in\mathcal P_L}\bigl|\mathrm{Tr}(\rho_\beta P)\bigr|
=
S_L(\beta)
=
\sum_{\substack{A,B\subset\{1,\dots,L\}\\ |A|=|B|}}
\bigl|\det G_{A,B}(\beta)\bigr|,
\]
the numerator becomes the sum of absolute values of all square minors of the \(L\times L\) Majorana correlation matrix \(G(\beta)\). Introducing the \(2L\times2L\) antisymmetric lift \(\mathcal A(G)\) and the selector matrix \(\mathcal J_{2L}\), with a staggered Jordan–Wigner gauge conjugation by
\[
D_s=\mathrm{diag}\bigl((-1)^1,(-1)^1,\dots,(-1)^L,(-1)^L\bigr),
\]
one obtains the exact Pfaffian identity
\[
S_L(\beta)=(-1)^L\operatorname{Pf}
\begin{pmatrix}
\mathcal A(G(\beta)) & I_{2L}\\
-I_{2L} & -D_s\mathcal J_{2L}D_s
\end{pmatrix}.
\]
After block-wise site reordering, this becomes
\[
S_L(\beta)=(-1)^L\Pf\bigl[T_L(\Phi)+H_L(\Psi)\bigr],
\]
a finite-size block Toeplitz–Hankel Pfaffian [2606.08606].

The exact reformulation yields asymptotics in four thermal regimes. For fixed \(T\) and \(L\to\infty\),
\[
\ln S_L(\beta)=L f_0(\beta)+c_0(\beta)+o(1),
\]
with
\[
f_0(\beta)=\frac1{4\pi}\int_0^{2\pi}\ln\!\bigl[1+2\,t_\beta(\theta)/\sin(\theta/2)+t_\beta(\theta)^2\bigr]\,d\theta,\qquad
t_\beta(\theta)=\tanh\!\bigl(2\beta J\sin(\theta/2)\bigr).
\]
The purity part has an analogous expansion \(\Pi_L(\beta)=Lp_0(\beta)+p_1(\beta)+o(1)\), so
\[
M_{1/2,L}(\beta)=L(2f_0-2p_0)+(2c_0-2p_1)+o(1).
\]
In the high-\(T\) regime, \(\beta\to0\) at fixed \(L\),
\[
\ln S_L(\beta)=L f_{\rm HT}(\beta)+c_{\rm HT}(\beta)+\dots,
\]
with
\[
f_{\rm HT}(\beta)=2\beta-3\beta^2+\tfrac{16}{3}\beta^3+\cdots,\qquad
c_{\rm HT}(\beta)=-\beta+\tfrac52\beta^2-\tfrac{20}{3}\beta^3+\cdots.
\]
In the saturated low-\(T\) regime, \(L\to\infty\) with \(\beta/L\to\infty\),
\[
\ln S_L(\infty)=L f_\infty-\tfrac18\ln L+c_\infty+o(1),
\]
where
\[
f_\infty=\frac1{2\pi}\int_0^\pi \ln\!\bigl(2+2/\sin x\bigr)\,dx,
\]
and the \(-\tfrac18\ln L\) term is the boundary Fisher–Hartwig logarithm [2606.08606].

In the finite-size thermal crossover \(\beta=\tau L\), the stabilizer quantity factorizes:
\[
\ln S_L(\tau L)=\ln S_L(\infty)+\ln\mathcal F(\tau)+o(1),
\]
or equivalently
\[
M_{1/2,L}(\tau L)=2L(f_\infty-\ln2)-\tfrac14\ln L+2c_\infty+\ln\mathcal G(\tau)+o(1),
\]
where \(\mathcal G(\tau)=\mathcal F(\tau)^2/\mathcal P(\tau)^2\). The universal factor is
\[
\mathcal F(\tau)=
\frac{\eta(i\tau/2)\,\eta(2i\tau)^{7/4}}
{\eta(i\tau)^2\,\eta(4i\tau)^{1/2}},
\]
a level-eight eta quotient, rather than the ordinary free-boundary Majorana thermal factor
\[
\mathcal F_{ff}(\tau)=\prod_{p=1}^\infty (1+q^{2p-1})^{-1},\qquad q=e^{-\pi\tau}.
\]
As \(\tau\to0\),
\[
\ln\mathcal F(\tau)=-\frac{\pi}{16\tau}-\frac18\ln\tau+\frac18\ln2+o(1),
\]
which the paper interprets in terms of an effective “Pauli-weight depletion central charge”
\[
c_{\rm dep}^P=\tfrac34.
\]
The reported conclusion is that finite-temperature stabilizer entropy reveals hidden defect-like conformal data invisible to ordinary thermodynamic probes [2606.08606].

## 5. Quantum Monte Carlo evaluation in spin systems

Two complementary Monte Carlo strategies have been developed for many-body SRE calculations in sign-problem-free spin systems. The first is a non-equilibrium QMC algorithm based on the path integral of the work between two partition-function ensembles, implemented within stochastic-series-expansion QMC and applicable to all spatial dimensions and temperatures [2405.19577]. The second is a reduced-Pauli-string sampling method that treats \(\alpha\)-SRE as partition-function ratios and eliminates the sign problem in the imaginary-time path integral by sampling reduced Pauli strings within a reduced configuration space [2501.12146].

In the non-equilibrium construction, one thermalizes \(2n\) independent replicas at \(\lambda=0\), evolves through a schedule of \(\lambda\)-values while updating connection topology and SSE operator configurations, accumulates the work
\[
W=-\int d\lambda\,\partial_\lambda \ln g(\lambda,|B|),\qquad g(\lambda,k)=\lambda^k(1-\lambda)^{N-k},
\]
and finally estimates
\[
M_n(\beta)=-\frac{1}{1-n}\ln\langle e^{-W}\rangle.
\]
The reported time per path is \(\mathcal O(\beta N)\), with number of \(\lambda\)-steps chosen \(\sim10^4\ldots10^3\) independent of \(N\). Empirically, the signal-to-noise ratio for \(e^{-W}\) scales as \(N^{-\alpha}\) with \(\alpha\simeq1.20\), so \(R=\mathcal O(N^\alpha)\) samples suffice and the total CPU cost is \(\sim \mathcal O(R\beta N)=\mathcal O(\beta N^{1+\alpha})\approx \mathcal O(\beta N^{2.2})\) [2405.19577].

Benchmarks in the transverse-field Ising model show quantitative agreement with tensor-network-based algorithms. In a 1D ring with \(N=20\) and periodic boundary conditions, ground-state SRE \( \widetilde m_2\equiv S_2/N \) peaks near \(h/J=1\), and QMC agrees to within \(10^{-3}\) of exact MPS-contraction data. In a 2D cylinder \(4\times4\), ground-state and finite-\(T\) SRE \( \widetilde m_2 \) versus \(h/J\) or \(\beta\) agree with MPS with bond dimension \(\chi\) up to \(640\) to within \(1\%\). In a 2D square \(L\times L\) with \(L\) up to \(16\), the finite-\(T\) SRE density at \(h/J=2.75\) versus \(\beta\) displays two extremal points at \(\beta\approx0.36\) and \(0.78\), and a feature near the critical \(\beta_c\approx1.0874\) [2405.19577].

The reduced-configuration-space approach of Ding, Wang, and Yan addresses a different bottleneck. By replacing direct sampling of \(\{X,Y,Z\}\) strings with a reduced set \(\tilde{\mathcal P}_N\) and valid SSE configurations obeying the even-up constraint across replicas, it renders the four-replica generalized partition function \(Q\) sign-problem free and enables efficient classical computations of \(\alpha\)-SRE and its derivatives in previously inaccessible \(2\)D and higher-dimensional systems [2501.12146]. This paper emphasizes that the method gives scalable access not only to the full \(\widetilde M_2\) but also to its separate free-energy and characteristic-function contributions.

## 6. Interpretation, limitations, and points of tension

A central tension in the literature concerns whether finite-temperature stabilizer Rényi entropy should be regarded as a mixed-state magic measure. Ding, Wang, and Yan state that \(2\)-SRE fails to characterize magic in mixed states, yielding nonphysical results [2501.12146]. In the 2D transverse-field Ising test at \((J=1,h=2.5)\), the true finite-\(T\) critical point is \(\beta_c\simeq0.7851\), but the measured \(d\widetilde M_2/d\beta\) shows a clear singularity at a lower “ghost” temperature \(\beta^*\approx0.687\). Their decomposition of
\[
-\beta^2\frac{d^2\widetilde m_2}{d\beta^2}
=
C_{\beta,Q}-2C_{\beta,Z}-C_{\beta,Z_2}
\]
shows that the \(Z\)-part diverges exactly at \(\beta_c\), the \(Z_2\)-part diverges at \(\beta_c/2\), and the \(Q\)-part diverges at \(\beta^*<\beta_c\). Because \(\widetilde M_2\) mixes these three contributions with fixed coefficients \(\{1,-2,-1\}\), the total exhibits a nonphysical singularity at \(\beta^*\). The same paper concludes that no meaningful volume-law expansion survives at nonzero \(T\), and that the finite-\(T\) \(\widetilde M_2\) neither scales simply nor locates phase transitions correctly [2501.12146].

By contrast, the coupled SYK work and the open-chain Ising work assign finite-temperature stabilizer observables a different role. In the coupled SYK model, SRE is presented as an order parameter for a new class of finite-\(T\) transitions which cannot be detected by conventional thermodynamic observables, including an intrinsic SRE-only jump at \(\beta^*_{\mathrm{SRE}}\approx19\) [2509.17417]. In the open critical Ising chain, stabilizer Rényi-\(\tfrac12\) entropy reveals hidden defect-like conformal data invisible to ordinary thermodynamic probes, encoded in a level-eight eta quotient rather than the ordinary free-boundary Majorana factor [2606.08606].

Taken together, these results distinguish two uses of finite-temperature stabilizer entropy. One use is as a mixed-state magic monotone, where the criticism of \(\widetilde M_2\) for Gibbs states is explicit [2501.12146]. The other use is as an informational order parameter or probe of replica structure, boundary data, and saddle reorganization beyond standard thermodynamics [2509.17417; 2606.08606]. This suggests that the phrase “finite-temperature stabilizer Rényi entropy” now covers a mathematically coherent family of replica observables whose physical interpretation depends strongly on context, Rényi index, and the distinction between resource-theoretic monotonicity and diagnostic sensitivity.

## 7. Relation to entanglement, thermodynamics, and many-body structure

The motivation for finite-temperature SRE is rooted in the separation between quantum entanglement and quantum magic as distinct resources [2509.17417]. In the coupled SYK analysis, entanglement-type thermodynamic quantities can remain smooth while the SRE jumps, so the stabilizer observable isolates structure in the replicated boundary-twist sector that ordinary free energy does not register. In the open critical Ising chain, the relevant extra information is boundary or defect-like conformal data hidden from standard thermal factors. In spin-model Monte Carlo studies, the decomposition into \(Q\), \(Z\), and \(Z_2\) likewise separates a characteristic-function contribution from ordinary free-energy terms [2501.12146].

At zero temperature, Ding, Wang, and Yan fit the SRE on an \(L^d\) cluster as
\[
\widetilde M_2(L)=a_dL^d+b_d,
\]
and report that the subleading constant \(b_d\), the “volume-law correction,” jumps discontinuously as \(J\) crosses the critical point \(J_c\). They interpret \(b_d\neq0\) as genuine nonlocal many-body magic in correlations and propose that \(b_d\) is a sharper diagnostic of criticality than the full SRE density \(a_d=L^{-d}\widetilde M_2\) [2501.12146]. For finite temperature, however, that paper rejects the same observable as a faithful magic measure, whereas the SYK and conformal-boundary analyses emphasize precisely its ability to register structures invisible to ordinary thermodynamics.

The resulting research program is therefore two-sided. On one side are exact and large-\(N\) results showing that finite-temperature stabilizer observables can be computed analytically, reduced to Pfaffians, or solved by saddle points, and that they detect replica-connectivity changes or hidden boundary data [2509.17417; 2606.08606]. On the other side are quantum Monte Carlo results showing both that such quantities can be evaluated efficiently in large sign-problem-free systems and that, for Gibbs states, the resulting number need not behave as a physically faithful mixed-state magic measure [2405.19577; 2501.12146]. Within current literature, finite-temperature stabilizer Rényi entropy is therefore best understood as a technically rich family of replica-based observables whose interpretive status remains model- and purpose-dependent.

Source: https://www.emergentmind.com/topics/finite-temperature-stabilizer-renyi-entropy