Evaluating many-body stabilizer Rényi entropy by sampling reduced Pauli strings: singularities, volume law, and nonlocal magic

Abstract

We present a novel quantum Monte Carlo method for evaluating the $\hat{I}\pm$-stabilizer Rényi entropy (SRE) for any integer $\hat{I}\pm \ge 2$. By interpreting $\hat{I}\pm$-SRE as partition function ratios, we eliminate the sign problem in the imaginary-time path integral by sampling \emph{reduced Pauli strings} within a \emph{reduced configuration space}, which enables efficient classical computations of $\hat{I}\pm$-SRE and its derivatives to explore magic in previously inaccessible 2D/higher-dimensional systems. We first isolate the free energy part in $2$-SRE, which is a trivial term. Notably, at quantum critical points in 1D/2D transverse field Ising (TFI) models, we reveal nontrivial singularities associated with the \emph{characteristic function} contribution, directly tied to magic. Their interplay leads to complicated behaviors of $2$-SRE, avoiding extrema at critical points generally. In contrast, analyzing the volume-law correction to SRE reveals a discontinuity tied to criticalities, suggesting that it is more informative than the full-state magic. For conformal critical points, we claim it could reflect nonlocal magic residing in correlations. Finally, we verify that $2$-SRE fails to characterize magic in mixed states (e.g. Gibbs states), yielding nonphysical results. This work provides a powerful tool for exploring the roles of magic in large-scale many-body systems, and reveals intrinsic relation between magic and many-body physics.

15 pages, 12 figures