{"id":"6e4742bd-747c-4ddd-bf46-d96b4c389135","arxiv_id":"2607.05343","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":8.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":3,"one_line_summary":"A unified Magic Rényi Entropy measure for spins, bosons, and fermions is shown to have a universal critical contribution determined by the Affleck-Ludwig boundary entropy.","lead":"This paper introduces a unified measure called the Magic Rényi Entropy (MRE) that quantifies quantum computational resources across spins, bosons, and fermions. It uses conformal field theory to show that the universal part of this measure at criticality is governed by the Affleck-Ludwig boundary entropy, and confirms predictions with numerical calculations.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"Conformal invariance of the boundary state in the rotated-bulk theory is the right concern, but it is well-addressed for the TLL; the more practical risk is the very small system sizes used for numerical verification of non-perturbative predictions.","rationale":"The reader identified the correct load-bearing assumption — conformal invariance of the boundary state in the rotated-bulk theory. However, for the TLL (the paper's main concrete application), this concern is addressed through a chain of arguments: (1) the unitary rotation preserves the bulk spectrum, (2) the Gaussian part H_0 explicitly satisfies the conformal boundary condition, (3) the non-Gaussian perturbation is exactly marginal, and (4) the bulk-boundary OPE analysis systematically identifies all relevant boundary channels. The assumption about the conformal manifold (footnote [81]) is verified for the TLL but would need separate justification for other models.\n\nThe analytical framework is internally consistent: the Liouville-space path integral formulation is exact, the BCFT extraction of the g-factor follows standard techniques, the perturbative calculation involves a clean torus correlator computation with a factorization that is verified through charge neutrality, and the duality K↔K⁻¹ provides a non-trivial consistency check. The cancellation of all γ²L/β terms (noted after Eq. 151) is a useful internal check that the result is a genuine boundary quantity.\n\nThe numerical verification is the weakest link, but it does not undermine the central analytical claim. The perturbative prediction near K=1 is confirmed, and the boundary transition prediction comes from an independent analytical result (the scaling dimension h_min). The small system sizes (L ≤ 14) limit confidence in the non-perturbative regime, but the paper is transparent about finite-size effects and the agreement with perturbation theory where it should hold is excellent.\n\nThe paper makes a substantive contribution: it provides the first unified field-theoretical framework for many-body magic across spins, bosons, and fermions, with a concrete analytical calculation and numerical verification for the fermionic case. The ACCEPT verdict with HIGH confidence is appropriate — the central claim is well-supported for the TLL, and the limitations (small numerics, scope of the conformal manifold assumption) are honestly discussed.","tokens_in":54588,"tokens_out":5427,"duration_ms":228099,"concrete_test":"Perform DMRG or tensor-network calculations for the interacting fermion chain (Eq. 109) at L = 32, 64, 128 to extract the universal constant s(K) and locate the boundary transition. If the crossing point K_c shifts by more than 0.1 from the L≤12 value of 3.05, or if s(K) away from K≈1 deviates by more than 20% from the L≤12 extrapolation, the non-perturbative predictions would need revision. Near K=1, the perturbative formula s = 3γ²/2 should remain stable, providing an internal consistency check.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader correctly identifies the conformal invariance of |J⟩⟩ in the rotated-bulk picture as the load-bearing assumption. However, the concern is less severe than it first appears for the paper's main concrete application (the TLL). The rotation is unitary, so the bulk spectrum and universality class are unchanged (Sec. V.B). The paper decomposes the rotated Hamiltonian as H_0 + H_Δ^rot, explicitly verifies (T−T̄)|J⟩⟩=0 for the Gaussian part H_0 (Eq. 136), shows H_Δ^rot is exactly marginal via the SU(2) current algebra (Appendix E), and performs a systematic bulk-boundary OPE analysis to identify all potentially relevant boundary operators (Appendix D, Eq. 163). The assumption that H_0, H_R, and H_R^rot lie on the same conformal manifold (footnote [81]) is verified for the TLL but remains an unverified assumption for the broader class of density-density models in Sec. V — though the paper is appropriately scoped to the TLL for its detailed results.\n\nThe more practical concern is the numerical verification. The exact diagonalization uses L ≤ 14 sites, with the universal constant s extracted from three-point fits at L = 8, 10, 12 (Eq. 165). The boundary transition at K_c = 3 is located via a crossing analysis using only L₁, L₂, L₃ = 12, 10, 8 (Eq. 168-169), giving K_c = 3.05. At these sizes, higher-order corrections to the scaling form (1/L² terms, etc.) could shift both the extracted s values and the crossing point. The perturbative prediction s = 3γ²/2 near K=1 is confirmed numerically, which validates the analytical framework. But the non-perturbative predictions — the boundary transition location and the behavior of s away from K=1 — rest on very small systems where asymptotic scaling is not guaranteed.","agreement_with_reader":"partial"},"referee_report":{"model":"glm-5.2","summary":"This paper introduces the magic Renyi entropy (MRE), a unified measure of quantum computational resources (nonstabilizerness for spins, non-Gaussianity for bosons/fermions) based on convolution of replicated states. The authors formulate the MRE within a Euclidean path-integral framework and show that, at one-dimensional critical points, its universal size-independent contribution is determined by the Affleck-Ludward boundary entropy (g-factor) of a conformal boundary condition in the replicated theory. As a concrete application, they analyze the Tomonaga-Luttinger liquid (TLL) of interacting spinless fermions: they derive perturbatively that the universal constant s = (3/2)gamma^2 + O(gamma^4) near the free-fermion point (K=1), predict boundary phase transitions at K=1/3 and K=3 via a bulk-boundary OPE analysis, and confirm these predictions with exact diagonalization (finding K_c ~ 3.05 vs. the predicted 3). The paper also proves key resource-theoretic properties (faithfulness, additivity, monotonicity under Gaussian protocols for bosons) and extends the qubit/qudit MRE equivalence with the SRE to general coprime (n,d).","tokens_in":55144,"tokens_out":1509,"duration_ms":199226,"significance":"The paper makes a substantial contribution by providing the first unified field-theoretical framework for quantum magic across spins, bosons, and fermions. The derivation of the universal MRE constant from the Affleck-Ludward g-factor is a parameter-free result: the perturbative formula s = (3/2)gamma^2 is derived analytically from BCFT data and then compared against independent numerics, not fitted. The systematic bulk-boundary OPE analysis (Appendix D) that predicts the boundary transitions at K=1/3 and K=3 is a falsifiable, non-trivial prediction. The proof of monotonicity under adaptive bosonic Gaussian protocols (Theorem 11 in Appendix A) and the qudit MRE-SRE equivalence (Appendix B) are additional strengths. The cancellation of gamma^2 L/beta terms in the perturbative calculation (Sec. VI.C, Eq. 151) provides a strong internal consistency check. The main limitation is the small system sizes (L <= 14) used for numerical verification, which leaves the non-perturbative regime and the boundary transition location only qualitatively confirmed.","major_comments":[{"comment":"Sec. VI.E, Eqs. (165)-(169): The numerical verification uses very small system sizes (L = 8, 10, 12 for the three-point fits; L1, L2, L3 = 12, 10, 8 for the crossing analysis). At these sizes, subleading corrections to the scaling form M = mL - s + d/L (i.e., 1/L^2 and higher-order terms) could be significant and may shift both the extracted s values and the crossing point K_c. The agreement K_c ~ 3.05 vs. 3 is encouraging but could be fortuitous given the finite-size effects. The authors should discuss the potential impact of these corrections more explicitly and, if possible, provide an estimate of the systematic error. This is load-bearing because the numerical confirmation of the boundary transition is a central claim of the paper.","section":null},{"comment":"Sec. V.B, footnote [81]: The assumption that H_0, H_R, and H_R^rot lie on the same conformal manifold is verified for the TLL but stated as an unverified assumption for the broader class of density-density models in Sec. V. While the paper is appropriately scoped to the TLL for its detailed results, the general framework in Sec. V is presented as applicable to this broader class. The authors should clarify the scope of their general claims: are the results of Sec. V (the three scenarios for boundary RG flows) presented as general principles that should hold when the conformal-manifold condition is met, or are they claimed to hold for all density-density models? A brief clarifying statement would strengthen the presentation.","section":null}],"minor_comments":[{"comment":"Fig. 5 caption: The notation s(K, L_0) is introduced in the text but the figure caption could benefit from explicitly stating that the plotted quantity is the fitted universal constant at central size L_0.","section":null},{"comment":"Eq. (110): The expressions for v and K are given in terms of arccos(V/2t). It would help the reader to note that these follow from the standard Bethe ansatz solution of the model (109), or to cite the relevant reference.","section":null},{"comment":"Sec. VI.E: The boundary conditions for ED are stated as antiperiodic (periodic) when N = L/2 is even (odd). It would be useful to briefly explain why this choice is made (presumably to avoid degeneracies and access the ground state in the NS sector).","section":null},{"comment":"Appendix F, Eq. (F3): The gluing matrix G is an 8x8 matrix. For readability, it may help to indicate the block structure (e.g., that it decomposes into 4x4 blocks) more explicitly in the text.","section":null},{"comment":"Sec. III.C, Eq. (56): The bosonic MRE is written as a phase-space integral. It may be worth noting that this integral may require regularization for certain states, or stating the class of states for which it is well-defined.","section":null},{"comment":"Typo in Sec. VI.C, Eq. (147): The expression uses both E_2(q) and vartheta_3(q); the notation should be consistent with standard conventions (e.g., specifying whether E_2 is the normalized or unnormalized Eisenstein series).","section":null},{"comment":"Sec. VII (Outlook): The discussion of the SYK model as a future direction is interesting but somewhat speculative. A brief mention of what specific non-Gaussianity signature the MRE might reveal in the SYK model would make this outlook more concrete.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The reader's report and stress-test note correctly identify the conformal invariance of the boundary state in the rotated-bulk picture as the load-bearing assumption. On reading the paper, I find this concern is well-addressed for the TLL: the rotation is unitary (preserving the spectrum), the Gaussian part is explicitly verified to be conformally invariant (Eq. 136), the perturbation is shown to be exactly marginal (Appendix E), and a systematic OPE analysis identifies all relevant boundary operators (Appendix D). The more practical concern about small ED system sizes is valid and should be addressed in revision, but it does not undermine the central analytical results. The paper is a strong contribution to the emerging field of many-body magic and is appropriate for a high-impact condensed matter or quantum information journal."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful and constructive report. The referee correctly identifies the central results of the paper and raises two substantive points: (1) the potential impact of subleading finite-size corrections on the numerical extraction of the universal constant s and the boundary transition point K_c, and (2) the scope of the general framework in Sec. V regarding the conformal-manifold assumption. Both points are well-taken, and we will address them in a revised manuscript.","responses":[{"response":"The referee raises a valid concern. We agree that at the system sizes accessible to exact diagonalization (L <= 14), subleading 1/L^2 corrections to the scaling form M = mL - s + d/L can be non-negligible and may shift both the extracted s values and the crossing point K_c. We will revise the manuscript to discuss this limitation more explicitly. Specifically, we will add a paragraph in Sec. VI.E acknowledging that: (i) the 1/L^2 and higher-order corrections are not included in the three-point fit, and their omission introduces a systematic error that is difficult to quantify at these sizes; (ii) the agreement K_c ~ 3.05 vs. the predicted K_c = 3 is encouraging and consistent with the field-theoretical prediction, but we cannot rule out that finite-size effects shift the crossing point; (iii) the perturbative formula s = (3/2)gamma^2 is confirmed most robustly near K = 1, where the universal constant is largest relative to the finite-size corrections, and the agreement there provides the strongest numerical support for the theory. We will also note that the duality collapse in Fig. 6 (using L = 10, 12, 14) provides an additional consistency check that is less sensitive to the absolute value of s. We emphasize that the central analytical results — the perturbative formula, the boundary stability criterion, and the predicted transition points — are derived from BCFT data and do not depend on the numerics. The numerics serve as a qualitative confirmation, and we will make this framing clearer in the revised text.","revision_made":"partial","referee_comment":"Sec. VI.E, Eqs. (165)-(169): The numerical verification uses very small system sizes (L = 8, 10, 12 for the three-point fits; L1, L2, L3 = 12, 10, 8 for the crossing analysis). At these sizes, subleading corrections to the scaling form M = mL - s + d/L (i.e., 1/L^2 and higher-order terms) could be significant and may shift both the extracted s values and the crossing point K_c. The agreement K_c ~ 3.05 vs. 3 is encouraging but could be fortuitous given the finite-size effects. The authors should discuss the potential impact of these corrections more explicitly and, if possible, provide an estimate of the systematic error."},{"response":"We agree that the scope of the general claims in Sec. V should be stated more precisely. The intention is that the three scenarios for boundary RG flows (Sec. V.B) are presented as general principles that apply when the conformal-manifold condition is met — that is, when H_0, H_R, and H_R^rot lie on the same conformal manifold. This condition is verified for the TLL (Sec. VI) but is not established for the broader class of density-density models in general. We will add a clarifying statement at the beginning of Sec. V.B making explicit that the three scenarios are conditional on the conformal-manifold assumption, and that the detailed results of the paper (perturbative calculations, boundary transition predictions) are derived specifically for the TLL where this assumption is verified. This scoping is already implicit in the structure of the paper but should be stated more directly.","revision_made":"yes","referee_comment":"Sec. V.B, footnote [81]: The assumption that H_0, H_R, and H_R^rot lie on the same conformal manifold is verified for the TLL but stated as an unverified assumption for the broader class of density-density models in Sec. V. While the paper is appropriately scoped to the TLL for its detailed results, the general framework in Sec. V is presented as applicable to this broader class. The authors should clarify the scope of their general claims: are the results of Sec. V (the three scenarios for boundary RG flows) presented as general principles that should hold when the conformal-manifold condition is met, or are they claimed to hold for all density-density models? A brief clarifying statement would strengthen the presentation."}],"tokens_in":54650,"tokens_out":971,"duration_ms":176367,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"The main thing to know: this paper constructs a single convolution-based magic measure (the MRE) that reduces to the SRE for qubits/qudits and quantifies non-Gaussianity for bosons and fermions, then shows that at 1D criticality the universal size-independent term is the Affleck-Ludwig boundary g-factor. That is a genuine unification and a clean field-theoretic result. The TLL application is concrete: the perturbative formula s = (3/2)γ² near K=1 is derived analytically with no free parameters, the γ² L/β terms cancel as a consistency check, and the boundary stability analysis predicting transitions at K=1/3 and K=3 follows from a systematic bulk-boundary OPE. The resource-theoretic properties (faithfulness, additivity, monotonicity under Gaussian protocols for the Helmert convolution) are proven carefully in Appendix A, including the adaptive measurement case. The qudit generalization in Appendix B, showing MRE = SRE when n and d are coprime, is a nice bonus. The reader's concern about conformal invariance of the boundary state in the rotated-bulk picture is the right question to ask, but the paper handles it well for the TLL: the rotation is unitary so the bulk spectrum is unchanged, (T−T̄)|J⟩⟩=0 is verified for the Gaussian part, the perturbation is shown exactly marginal via the SU(2) current algebra, and the bulk-boundary OPE systematically identifies all relevant boundary operators. The assumption that H₀, H_R, and H_R^rot lie on the same conformal manifold is verified for the TLL but stated as an assumption for the broader class in Section V — this is appropriately scoped. The real soft spot is the numerics. ED uses L ≤ 14, with the universal constant extracted from three-point fits at L = 8, 10, 12. The boundary transition at K_c = 3 is located via a crossing analysis on those same three sizes, giving K_c ≈ 3.05. At these system sizes, 1/L² corrections could shift both the extracted s values and the crossing point. The perturbative prediction near K=1 is confirmed convincingly, which validates the analytical framework. But the non-perturbative predictions — the transition location and the behavior of s away from K=1 — rest on very small systems where asymptotic scaling is not guaranteed. This is a limitation but not a fatal one, because the analytical results stand on their own; the numerics are confirmation, not load-bearing. This paper is for researchers working on quantum resource theories, many-body magic, and CFT methods. It deserves a serious referee who can check the BCFT construction and the perturbative calculation in detail.","headline":"Unified MRE measure for spins/bosons/fermions with BCFT derivation of universal term; numerics are small but the analytics carry the paper.","tokens_in":55527,"tokens_out":673,"would_cite":true,"duration_ms":116453,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.67.Mn","11.25.Hf","71.10.Pm"],"model":"glm-5.2","headline":"Quantum magic unifies spins, bosons, and fermions via boundary entropy","keywords":[],"falsifier":"If exact-diagonalization or other numerical methods at larger system sizes showed that the universal constant s does not follow the predicted (3/2)gamma^2 scaling near K=1, or that no boundary transition occurs near K=3, the central field-theoretical predictions would be falsified. Alternatively, if the conformal invariance of the boundary condition were found to be broken by the inter-replica couplings in the rotated bulk, the identification of s with the Affleck-Ludwig g-factor would fail.","tokens_in":54868,"feed_emoji":"🎲","tokens_out":1145,"duration_ms":85154,"temperature":0.7,"pith_summary":"This paper introduces the magic Renyi entropy (MRE), a single measure that quantifies the computational resource content of quantum many-body states regardless of whether the underlying degrees of freedom are spins, bosons, or fermions. For spins, the MRE reduces to the known stabilizer Renyi entropy, which measures nonstabilizerness. For bosons and fermions, it measures non-Gaussianity. The central claim is that when a one-dimensional quantum system is at a critical point described by a conformal field theory, the universal, size-independent contribution to the MRE is exactly the Affleck-Ludwig boundary entropy, a known quantity characterizing conformal boundary conditions. The authors construct the MRE via a convolution operation: multiple copies of a quantum state are mixed by a unitary, most copies are discarded, and the purity of the remaining copy is measured. Free states (stabilizer states for spins, Gaussian states for bosons and fermions) remain pure under this operation; resourceful states become mixed. Rewriting this construction as a Euclidean path integral on replicated copies of the system, the convolution becomes a boundary condition in a folded geometry. The resulting boundary entropy directly gives the universal term in the MRE. The paper then applies this framework to interacting spinless fermions described by a Tomonaga-Luttinger liquid. It shows that non-Gaussianity (deviation from the free-fermion point K=1) continuously renormalizes the boundary entropy through the identity channel of the bulk-boundary operator product expansion, yielding a perturbative formula s = (3/2)gamma^2 + O(gamma^4) where gamma = (K-1)/(K+1). Beyond a critical interaction strength, a relevant nonidentity boundary operator is generated, driving a boundary phase transition at K=1/3 and K=3. Exact-diagonalization numerics confirm the perturbative formula near K=1, the duality under K to 1/K, and locate the transition at K approximately 3.05, close to the predicted K=3.","feed_headline":"Quantum magic unifies spins, bosons, and fermions via boundary entropy","feed_subtitle":"A single measure of computational resource is governed by conformal boundary data across all three platforms, with fermionic interactionsdr ","key_machinery":"The magic Renyi entropy (MRE) is defined by mixing n copies of a quantum state via a convolution unitary, discarding n-1 copies, and measuring the purity loss of the remaining copy. When recast as a replicated Euclidean path integral and folded, the convolution becomes a boundary condition whose Affleck-Ludwig g-factor determines the universal, size-independent term of the MRE. For fermions in the Tomonaga-Luttinger liquid, the bulk-boundary OPE of the exactly marginal perturbation induced by convolution determines whether the boundary is stable (only identity and irrelevant channels) or undergoes a transition (relevant nonidentity channel generated).","core_discovery":"The universal contribution to quantum magic in one-dimensional critical many-body states is governed by the Affleck-Ludwig boundary entropy of a conformal boundary condition induced by a convolution operation on replicated states. This holds uniformly across spins, bosons, and fermions, and non-Gaussianity in fermionic systems can either continuously renormalize this entropy or drive a boundary phase transition at specific interaction strengths.","pith_inferences":[],"forward_implications":["The MRE provides a common language for comparing computational resource content across different physical platforms, enabling direct quantitative comparison of magic in spin-chain experiments, bosonic continuous-variable systems, and fermionic condensed-matter systems.","Boundary phase transitions in the MRE at K=1/3 and K=3 represent a new class of critical phenomena intrinsic to quantum computational resources, distinct from conventional entanglement-driven or symmetry-breaking transitions.","The connection between convolution-based free-state testing and boundary CFT data suggests that other resource-theoretic quantities defined through similar replica-mixing operations may also admit universal field-theoretical descriptions at criticality.","The framework can be applied to other critical fermionic systems, including the Sachdev-Ye-Kitaev model, where non-Gaussianity measured by the MRE may reveal aspects of quantum magic complementary to those captured by spin-based nonstabilizerness."],"fun_headline_variants":["Quantum computational resources tied to boundary entropy across spin, boson, fermion syste","Affleck-Ludwig boundary entropy governs many-body magic in spins, bosons, and fermions","Magic Rényi entropy unifies nonstabilizerness and non-Gaussianity via CFT boundary data","Conformal boundary entropy controls quantum magic in critical fermionic and bosonic chains","Unified magic measure reveals boundary phase transitions in Tomonaga-Luttinger liquids"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The analysis assumes that the rotated bulk Hamiltonian, which contains inter-replica couplings generated by the convolution unitary, lies on the same conformal manifold as the original decoupled theory, so that the boundary CFT machinery applies. If the rotated Hamiltonian flows to a different bulk universality class, the boundary entropy analysis breaks down.","fun_headline_variants_meta":{"raw":{"variants":["Quantum computational resources tied to boundary entropy across spin, boson, fermion systems","Affleck-Ludwig boundary entropy governs many-body magic in spins, bosons, and fermions","Magic Rényi entropy unifies nonstabilizerness and non-Gaussianity via CFT boundary data","Conformal boundary entropy controls quantum magic in critical fermionic and bosonic chains","Unified magic measure reveals boundary phase transitions in Tomonaga-Luttinger liquids"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":768,"prompt_tokens":652,"completion_tokens":116,"prompt_tokens_details":null},"tokens_in":652,"tokens_out":116,"duration_ms":18356,"temperature":1.0,"reasoning_tokens":null,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-07T16:09:41.859790+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If exact-diagonalization or other numerical methods at larger system sizes showed that the universal constant s does not follow the predicted (3/2)gamma^2 scaling near K=1, or that no boundary transition occurs near K=3, the central field-theoretical predictions would be falsified. Alternatively, if the conformal invariance of the boundary condition were found to be broken by the inter-replica couplings in the rotated bulk, the identification of s with the Affleck-Ludwig g-factor would fail.","supporting_citations":[],"review_version":1}