{"id":"586fb984-f455-40cb-b03c-26efd6b98f61","arxiv_id":"2411.17231","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"The continuum entropic c-function of the 2+1 dimensional Z2 gauge theory is reported to show a power-law short-distance regime and an exponential large-distance decay, with a crossover near l m_g = 1.","lead":"The authors use the Kramers-Wannier duality to turn a hard entanglement calculation in a 2+1 dimensional Z2 gauge theory into a spin-model simulation, and they extract the entropic c-function in the continuum limit. The result is compared with a holographic prediction about how entanglement signals confinement.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fit (15) allows a free decay rate alpha=0.360(19), while the predicted form (12) has alpha=1; without testing alpha=1 the claim of numerical confirmation of (12) is not supported.","rationale":"The reader's weakest assumption is Eq. (11), the duality preservation of the entropic c-function. This is indeed load-bearing and is not derived in the present paper; it is taken from the companion paper [21] and from [30]. However, a more immediate and internal problem is the fitted value of alpha in Eq. (15). Eq. (12), which the paper claims to confirm, has no free alpha; setting m equal to the mass gap fixes the decay rate exactly. The best fit gives alpha = 0.360(19), which is incompatible with 1 at the quoted precision. This means either the data are not in the asymptotic large-l regime where (12) applies, or the relevant mass scale is not m_g, or the functional form is right but the prefactor or kinematics differ. In all three cases, the sentence 'our analysis is the first numerical confirmation of the prediction done in [16]' is too strong without an additional argument. The correct statement would be that the data are consistent with an exponential decay at some scale, with the extracted scale 0.36/m_g. A test with alpha fixed to 1 would settle this. The reader's verdict of CONDITIONAL is therefore appropriate: the numerical work is valuable and the duality approach is well motivated, but the interpretation must be revised before the confirmation claim is accepted. This is a precise, checkable mismatch between the fitted model and the prediction, not a dispute about consensus or a personal critique.","tokens_in":8291,"tokens_out":3151,"duration_ms":29493,"concrete_test":"Re-fit the continuum-extrapolated C2 data shown in Fig. 2 with alpha fixed to 1 and A free, reporting chi2/dof and residual plots; independently extract the large-l decay rate from the raw data and compare it with 2 m_g. If the fixed-alpha fit is rejected or the extracted mass differs from m_g by more than about two sigma, the claim of numerical confirmation of (12) must be withdrawn or qualified as confirmation only of an exponential form with a different mass scale.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim is that the continuum-limit entropic c-function of the (2+1)-dimensional Z2 gauge theory constitutes 'the first numerical confirmation' of the Klebanov-Kutasov-Murugan prediction (12). That prediction is the specific function C2 = l * integral dk exp(-2 sqrt(m^2+k^2) l), with m equal to the mass gap of the theory. In Eq. (15) the authors instead fit with an extra parameter alpha multiplying the mass in the exponent, finding alpha = 0.360(19). Since (12) corresponds to alpha = 1, the quoted result is about 34 sigma away from the predicted decay rate. A fit with alpha free tests only the exponential envelope, not the predicted mass scale; many functional forms with a free scale parameter can fit a smooth decaying curve. The paper does not report a fit with alpha fixed to 1, nor does it explain why the effective mass scale should be 0.36 m_g rather than m_g. The claimed crossover at l m_g = 1 is also inferred from where the power-law and exponential fits break down; with the fitted alpha the exponential's natural length scale is 1/(2 alpha m_g) ~ 1.4/m_g, so the crossover location is not an independent check of the prediction. In addition, the entire route relies on Eq. (11), C_n^Ising = C_n^gauge, imported from the companion paper [21] and not derived here; if that equality holds only for a particular center choice or fails in the replicated slab geometry, the gauge-theory interpretation collapses. But the alpha mismatch is the immediate, internally checkable inconsistency: the reported data are consistent with a one-parameter family of exponential decays, not specifically with (12).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript is a Lattice 2024 proceedings contribution in which the authors use Kramers-Wannier duality to compute the entropic c-function C2 of the (2+1)-dimensional Z2 gauge theory from Monte Carlo simulations of the dual Ising model. After reviewing the need for an unambiguous treatment of entanglement in gauge theories, they recall the dual replica construction of [21] and the statement that the Kramers-Wannier duality preserves C_n. They then simulate a two-replica Ising system with a Jarzynski-based estimator, perform thermodynamic and continuum extrapolations, and fit the extrapolated C2(l) to a power law at short distances and to an exponential-integral form at large distances. The fitted large-distance form has amplitude A=0.33(3) and exponent parameter alpha=0.360(19). The authors conclude that these data provide the first numerical confirmation of the Klebanov-Kutasov-Murugan prediction (12) and that the power-law to exponential crossover occurs at l m_g approximately 1.","tokens_in":8627,"tokens_out":9112,"duration_ms":79146,"significance":"The computational approach is sound and the continuum extrapolation is a valuable step: using a spin-model dual to define a replica geometry, combined with a non-equilibrium Monte Carlo estimator, is a promising way to access entanglement quantities that are otherwise ambiguous in lattice gauge theories. The paper also makes good use of external inputs for the CFT normalization [44] and the mass gap [45], and the reported chi-squared values indicate that the quoted fits are internally stable. However, the central claim of confirmation of Eq. (12) is not supported by the fit actually reported, because the fitted exponent alpha deviates from the predicted value 1 by roughly 34 standard deviations. With the free-exponent fit, the result is better described as a measurement of an effective decay rate; the headline conclusion needs to be corrected or supported by a fit with alpha fixed to 1.","major_comments":[{"comment":"The large-distance fit does not test the prediction (12), because Eq. (15) introduces a free parameter alpha multiplying the mass in the exponent. Equation (12) is the special case alpha=1; the fitted value alpha=0.360(19) is about 34 standard deviations from 1. With alpha free, the fit only checks that the data fall on an exponential envelope and that a smooth decaying curve can be described with an effective scale; it does not confirm that the mass-gap scale appears in the exponent as predicted. The manuscript reports no fit with alpha fixed to 1, no chi-squared for that constrained fit, and no discussion of why the effective decay rate should be 0.36 m_g. Since the 'first numerical confirmation' claim rests on exactly this comparison, the analysis must be redone or the claim must be weakened.","section":"Section 4, Eq. (15)"},{"comment":"The crossover at l m_g approximately 1 is presented as a physical finding, but it is inferred from the breakdown of two phenomenological fits. The fitted alpha affects the natural length scale of the exponential, 1/(2 alpha m_g) approximately 1.4/m_g, so the location of the crossover is not an independent confirmation of the mass-gap scale. To make the crossover claim meaningful, the authors should connect it to a fit in which the predicted mass scale is used, or present it as a property of the specific fitted forms.","section":"Section 4, Eqs. (14)-(15)"},{"comment":"The gauge-theory interpretation of the Ising c-function is imported from the companion paper [21] without derivation. Because the whole simulation route equates C_2^Ising with C_2^gauge for the replicated slab geometry in the continuum limit, the manuscript should either reproduce the key steps of the proof or state explicitly the conditions under which Eq. (11) holds, including the choice of center and boundary operators. As written, a reader cannot assess whether the measured curve is the gauge-theory c-function or a dual, center-dependent quantity.","section":"Section 3, Eq. (11)"}],"minor_comments":[{"comment":"There is a typo in the opening sentence: 'One one hand' should be 'On the one hand'.","section":"Section 1"},{"comment":"After Eq. (10), 'Karamers-Wannier' should be 'Kramers-Wannier'.","section":"Section 3"},{"comment":"The phrase 'the inverse mass-gap of the theory m_g' is confusing: m_g is later used as a mass scale in the dimensionless combination l m_g. It should be written as 'the mass gap m_g' or 'the inverse mass gap 1/m_g'.","section":"Section 4"},{"comment":"The sentence 'First, this study [21] is the first numerical confirmation of the prediction (12)' is ambiguous, since [21] is the companion paper and the introduction states that [21] already performed a high-precision test of the conjectures of [16]. The authors should clarify whether the current work is the first continuum-extrapolated confirmation and how it relates to the earlier numerical results in [21].","section":"Section 4"},{"comment":"The figure would be easier to assess if the continuum-extrapolated data points were shown with their errors and if the fit covariance, or at least the chi-squared for the constrained alpha=1 case, were reported.","section":"Figure 2"}],"recommendation":"major_revision","confidential_remarks":"As a proceedings contribution, the paper can be concise, but the discrepancy between the reported fit and the claim needs to be addressed before publication. The authors should add a fit with alpha fixed to 1 and report the corresponding chi-squared; if that fit is poor, the conclusion should be revised. The reliance on [21] for Eq. (11) is acceptable if the companion paper is cited and the conditions are stated, but the current text is too terse for a standalone claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a well-executed proceedings summary of [21], but the central claim is not supported by the reported fit. The stress-test note gets it right. The authors fit Eq. (15) with a free alpha and get alpha=0.360(19), while the predicted form (12) has alpha=1. That's about 34 sigma away. They never report a fit with alpha fixed to 1, so 'first numerical confirmation' is an overstatement. Many smooth decaying curves can be fit by an exponential with a free scale; without testing the predicted mass scale, the data do not specifically confirm (12).\n\nThat said, there is real substance. The Kramers-Wannier route around the replica-trick ambiguity is sensible, and the numerical work looks careful: thermodynamic and continuum extrapolations, reduced chi-squared near 1, and external inputs for C_CFT and m_g from [44,45] avoid circularity. The paper is honest enough to report alpha itself, so the problem is in the interpretation, not the data.\n\nThe other soft spot: this is a proceedings report, and the text attributes the first confirmation to [21]. So the standalone novelty is limited. The equality (11) is also taken from [21]; it is parameter-free and based on standard duality, so I don't see it as circular, but a reader would need [21] to judge it.\n\nIf the authors fixed alpha to 1 and showed the fit fails, or explained why the effective scale is 0.36 m_g, the claim would be credible. As written, the conclusion should be softened to 'consistent with an exponential decay' rather than 'confirmation of (12).' I'd engage with the work, but I would not cite the proceedings version; the companion paper is the citable source. For a journal, I'd send to peer review only if the alpha issue is fixable; otherwise the overclaim justifies a desk reject. For a proceedings, it is acceptable as a summary once the claim is revised.","headline":"Solid numerics, but the 'first numerical confirmation' claim fails against the paper's own fit, which finds alpha=0.36 rather than the predicted alpha=1.","tokens_in":9178,"tokens_out":2273,"would_cite":false,"duration_ms":22470,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Kramers-Wannier duality lets the entropic c-function of (2+1)-dimensional Z2 gauge theory be computed from dual Ising simulations, yielding a power-law to exponential crossover at the mass-gap scale.","keywords":["entanglement entropy","lattice gauge theory","Kramers-Wannier duality","entropic c-function","Rényi entropy","confinement","Monte Carlo simulation","Ising model"],"falsifier":"A direct lattice calculation of the R\\'enyi entropy ratio in the Z2 gauge theory, using an alternative replica construction (e.g., electric- or magnetic-center subalgebra), could be compared with the Ising-side result; a mismatch would signal that duality does not preserve the c-function. Additionally, computing $C_3$ from a three-replica simulation of the Ising model and comparing it with a direct gauge-theory computation would test the equality (11) beyond the second R\\'enyi entropy.","tokens_in":8078,"feed_emoji":"⚛️","tokens_out":4602,"duration_ms":37886,"temperature":0.7,"pith_summary":"This paper aims to show that the entanglement structure of a lattice gauge theory can be read off from its dual spin model. The authors exploit Kramers-Wannier duality to compute the entropic c-function of the (2+1)-dimensional Z2 gauge theory in the continuum limit, by simulating the dual Ising model. They report that the resulting c-function follows a power law at short distances and crosses over to an exponential decay when the slab thickness is comparable to the inverse mass gap, matching the behavior predicted by holographic models. If correct, this provides the first numerical confirmation of that prediction in a non-holographic confining theory, and establishes a general strategy for entanglement calculations in Abelian gauge theories.","feed_headline":"Mass gap sets entanglement crossover in Z2 gauge theory","feed_subtitle":"Simulating the dual Ising model yields the first numerical check of the holographic prediction in a non-holographic confining theory.","key_machinery":"The argument rests on Kramers-Wannier duality, which maps the partition function of the 3D Ising model to that of the Z2 lattice gauge theory, and on the companion result [21] that this duality preserves the entropic c-function, $C_n^{\\text{Ising}} = C_n^{\\text{gauge}}$. The entropic c-function is defined as $C_n = \\frac{l^{D-1}}{|\\partial A|} \\frac{\\partial S_n}{\\partial l}$, where $S_n$ is the R\\'enyi entropy of a slab subsystem of thickness $l$, and the ratio of partition functions entering it is computed with a Jarzynski-based non-equilibrium Monte Carlo algorithm for two replicas. Thermodynamic and continuum limits are taken using scale setting from [42], mass gap from [45], and the conformal value $C_2^{\\text{CFT}}$ from [44] for normalization.","core_discovery":"The central claim is that the entropic c-function of the (2+1)-dimensional Z2 gauge theory, extracted through duality from Monte Carlo simulations of the 3D Ising model, is described by a power law $C_2 \\sim B/(l m_g)^c$ at short distances and by the exponential decay form (12) at large distances, with a crossover at $l m_g = 1$. After thermodynamic and continuum extrapolations, the fits give $B = 0.360(9)$, $c = 0.48(2)$, and $A = 0.33(3)$, $\\alpha = 0.360(19)$, with reduced chi-squares near unity. The authors state that this is the first numerical confirmation of the Klebanov-Kutasov-Murugan prediction in a (2+1)-dimensional, non-holographic theory, and that the crossover scale is set by the mass gap of the theory.","pith_inferences":["The same duality route could be applied to the $U(1)$ gauge theory, whose dual spin model is well known, to test whether the c-function crossover at $l m_g = 1$ is a universal feature of confining Abelian theories.","If the equality $C_n^{\\text{Ising}} = C_n^{\\text{gauge}}$ holds for higher R\\'enyi orders, the method would allow extraction of higher-order entanglement measures in gauge theories, which are otherwise extremely difficult to compute directly.","The fitted power-law exponent $c = 0.48(2)$ near $l m_g = 1$ may be related to the scaling dimension of the mass operator, a connection the paper does not explore but which could be tested against analytic predictions.","Combining this duality approach with flow-based sampling could extend the same strategy to continuous gauge groups, where non-equilibrium algorithms currently face performance limitations."],"forward_implications":["A universal strategy emerges for entanglement in Abelian gauge theories that admit a spin-model dual: simulate the dual model and map the result back to the gauge theory.","The results confirm that the holographic prediction of exponential decay of the entropic c-function holds beyond holography, in a genuinely non-holographic theory.","The crossover scale $l m_g = 1$ shows that the mass gap, rather than the number of colors, controls the transition from a short-distance power law to a long-distance exponential behavior.","The continuum-extrapolated c-function provides a benchmark for direct entanglement calculations in gauge theories and for tests of other duality-based approaches."],"supporting_citations":[{"why":"Supplies the holographic prediction of the exponential decay form (12) for the entropic c-function in confining theories, which the paper tests numerically.","marker":"[16]"},{"why":"Establishes the duality mapping of the replica partition-function ratio and the preservation of the entropic c-function under Kramers-Wannier duality, the core identity used here.","marker":"[21]"},{"why":"Provides the non-equilibrium Monte Carlo algorithm used to compute the ratio of partition functions with high precision.","marker":"[26]"},{"why":"Introduced the out-of-equilibrium protocol for R\\'enyi entropies via the Jarzynski equality that the algorithm of [26] generalizes.","marker":"[41]"},{"why":"Supplies the scale-setting procedure for the Z2 gauge theory used to convert lattice units into physical units.","marker":"[42]"},{"why":"Gives the conformal value $C_2^{\\text{CFT}}$ used to normalize the entropic c-function.","marker":"[44]"},{"why":"Provides the inverse mass gap $m_g$ used as the reference scale for the plots and fits.","marker":"[45]"},{"why":"The original Kramers-Wannier duality transformation that underlies the whole method.","marker":"[24,25]"}],"fun_headline_variants":["Z2 gauge entanglement crossover pinned by mass gap","Duality unlocks entanglement in Z2 lattice gauge theory","First numerical check of holographic prediction in Z2 gauge","Mass gap governs Z2 gauge entanglement across scales","Entropic c-function crossover from dual Ising model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire numerical route depends on the equality $C_n^{\\text{Ising}} = C_n^{\\text{gauge}}$ (Eq. 11) holding for the replicated slab geometry in the continuum limit; if Kramers-Wannier duality does not preserve the entropic c-function, the reported curve is not the gauge theory c-function.","fun_headline_variants_meta":{"raw":{"variants":["Z2 gauge entanglement crossover pinned by mass gap","Duality unlocks entanglement in Z2 lattice gauge theory","First numerical check of holographic prediction in Z2 gauge","Mass gap governs Z2 gauge entanglement across scales","Entropic c-function crossover from dual Ising model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000502,"raw_usage":{"total_tokens":2412,"prompt_tokens":864,"completion_tokens":1548,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":480,"completion_tokens_details":{"reasoning_tokens":1473}},"tokens_in":480,"tokens_out":1548,"duration_ms":11359,"temperature":1.0,"reasoning_tokens":1473,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:22:56.056509+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct lattice calculation of the R\\'enyi entropy ratio in the Z2 gauge theory, using an alternative replica construction (e.g., electric- or magnetic-center subalgebra), could be compared with the Ising-side result; a mismatch would signal that duality does not preserve the c-function. Additionally, computing $C_3$ from a three-replica simulation of the Ising model and comparing it with a direct gauge-theory computation would test the equality (11) beyond the second R\\'enyi entropy.","supporting_citations":[{"cited_title":"Universal divergence of the Renyi entropy of a thinly sliced torus at the Ising fixed point","cited_arxiv_id":"1904.08955","evidence_quote":"Gives the conformal value $C_2^{\\text{CFT}}$ used to normalize the entropic c-function."},{"cited_title":"The Spectrum of the 2+1 Dimensional Gauge Ising Model","cited_arxiv_id":"hep-lat/9607029","evidence_quote":"Provides the inverse mass gap $m_g$ used as the reference scale for the plots and fits."}],"review_version":1}