{"id":"5a089109-25d0-426d-98db-672c4cabec22","arxiv_id":"1908.08058","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In the transverse field XY spin chain, single-qubit magic turns on just after the critical point, peaks at the factorizable ground state, and two-qubit magic persists over long distances and traces the quantum critical region.","lead":"This paper calculates how much magic, a resource required for fault-tolerant quantum computation, appears in the spins of a magnetic chain near its quantum phase transition. It finds a special point where the ground state is a product of optimal magic states, and shows that magic between distant spins can detect the quantum critical region.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the missing |±i> optimality proof in Eq. (8) is supplied by a simple shift argument.","rationale":"The paper's central single-qubit claims all rest on Eq. (8): magic vanishes in the disordered phase, appears at an MPP near λ = 1.00015, and reaches √2 - 1 at the factorizable point γ = 1/3, λ = √(9/8). The reader correctly identified that Appendix A assumes, without proof, that decompositions over only the four xz-plane stabilizer states are optimal. I checked this and found the missing argument: for y-magnetization zero, the weights on |+i> and |-i> must be equal, and adding half that common weight to each of the four xz coefficients preserves normalization and all Pauli expectations while not increasing L1. Therefore the restriction is harmless and Eq. (8) is exact for the positive-order-parameter branch. The factorizable-point result is independently supported by the mean-field orientation x = √(2γ/(1+γ)), z = √((1-γ)/(1+γ)), which equals (1/√2,1/√2) exactly at γ = 1/3, giving pure H states. What remains are presentation and numerical-rigor issues: Eq. (6) is notationally inconsistent with Eq. (7) by a factor of two, fitted exponents μ, ν, κ and the δλ_c scaling are reported without error bars, and the claim that two-qubit QR peaks exactly at MPP is supported by DMRG and the derivative structure of QR rather than by an analytic derivation. None of these threatens the central construction, so I do not change the reader's CONDITIONAL verdict.","tokens_in":19390,"tokens_out":27065,"duration_ms":258789,"concrete_test":"Run the LP of Eq. (7) over all six single-qubit stabilizer states for the H-state Bloch vector (x,z)=(1/√2,1/√2) and for a grid of points with x,z ≥ 0 and x²+z² ≤ 1, verifying that the optimum equals max[|x|+|z|-1,0]; the shift argument predicts exact agreement, so any counterexample would refute Eq. (8).","verdict_should_be":"UNCHANGED","load_bearing_attack":"No significant objection identified. The reader's weakest assumption is that Appendix A restricts the RoM pseudo-mixture to the four xz-plane stabilizer states without proving that |+i> and |-i> cannot lower the L1 norm. This gap is real but the concern does not land. For any state with ⟨σy⟩=0, the weights on |+i> and |-i> in any pseudo-mixture must be equal; call this common weight t. Replace each weight a_i on |0>, |1>, |+>, |-> by a_i + t/2. The normalization is preserved because the new sum is (Σa_i) + 4(t/2) = (1 - 2t) + 2t = 1. The x and z expectations are preserved because the differences a1 - a2 and a3 - a4 are unchanged, and the y expectation remains zero. By the triangle inequality, Σ|a_i + t/2| ≤ Σ|a_i| + 2|t|, which is exactly the L1 contribution of the removed |±i> pair. Hence every six-state decomposition can be converted to a four-state decomposition with no larger RoM. For the symmetry-broken branch with x,z ≥ 0, the four-state LP gives R = max[|x| + |z| - 1, 0] = max[x + z - 1, 0]. Thus Eq. (8) is correct; only the presentation of Appendix A is incomplete. Remaining caveats are numerical rather than load-bearing: fitted exponents lack error bars, and the exact MPP-peak statement for two-qubit QR rests on DMRG and a derivative-kink argument rather than an analytic proof.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the robustness of magic (RoM) in single- and two-qubit reduced states of the transverse-field anisotropic XY spin chain. For the symmetry-broken ground state, it claims a simple closed-form expression for single-qubit RoM, R(λ) = max[⟨σx⟩ + ⟨σz⟩ − 1, 0], valid when ⟨σy⟩ = 0, and uses this to show that magic vanishes in the disordered phase, emerges at a 'magic pseudocritical point' (MPP) just inside the ordered phase, and reaches the maximal equatorial value √2−1 at the factorizable point (γ0 = 1/3, λ0 ≈ 1.06), where the ground state is a product of H states. The paper further investigates two-qubit RoM, demonstrating long-range persistence, a sharp maximum of a 'global magic' quantity at the MPP, scaling laws near criticality, finite-size scaling behavior, thermal sudden death of magic, and detection of the quantum critical region via the Gruneisen parameter. The results are obtained through a combination of analytic formulas, linear programming, and DMRG numerics.","tokens_in":19776,"tokens_out":13511,"duration_ms":113923,"significance":"If the results hold, the paper gives an operational resource-theoretic perspective on a classic integrable model, showing that factorizable ground states can serve as sources of pure H-type magic states for fault-tolerant quantum computation and that magic can act as a long-range probe of quantum criticality, complementing entanglement and discord. The analytic formula Eq. (8) is a useful simplification, and the identification of the factorizable point as a magic-state factory is a concrete, falsifiable prediction. The paper also provides extensive numerical evidence, including machine-checkable linear programs and DMRG data, for the scaling and thermal properties. However, the proof of the central formula in Appendix A is incomplete as written, and several numerical exponents are extracted without quantitative error estimates, which tempers the strength of the quantitative claims.","major_comments":[{"comment":"The proof of Eq. (8) assumes that the optimal RoM pseudo-mixture for states with ⟨σy⟩ = 0 uses only the four stabilizer states |0⟩, |1⟩, |+⟩, |−⟩, without justifying why the |+i⟩ and |−i⟩ states cannot yield a lower L1 norm. This is a genuine gap in the derivation of the central formula. The result is nevertheless correct: for any decomposition that includes |±i⟩ with equal weight t (required by ⟨σy⟩ = 0), replacing them by t/2 of each of the four xz-plane states preserves the state and, by the triangle inequality, does not increase the L1 norm. Please add this argument or a reference to a complete proof.","section":"Appendix A, Eq. (A.1)"},{"comment":"The claim that the global magic QR attains its maximum exactly at the MPP for all inter-site distances r is stated as an exact result, but it is supported only by DMRG numerics and a qualitative derivative-kink observation. Since this is a central result of the paper, please provide either an analytic argument for the coincidence or a more extensive numerical verification with error estimates and a scan over anisotropy γ.","section":"Section 3.2, Figs. 4(b) and 4(c)"},{"comment":"The finite-size scaling exponents μ and ν, as well as the infinite-chain exponents μ and the power law δλc ∼ γ^5.55, are determined by visual fits without quantitative uncertainties. Because these exponents are used to support the analytic prediction μ ≈ 1 − βx (Eq. (11)), please provide a systematic collapse analysis with error bars or at least quantitative fit-quality measures.","section":"Section 3.1.3 and Figs. 1(d), 2(b), 2(c), 3"},{"comment":"The 'global magic' QR = log(1 + R(ρ12)) − log(1 + R(ρ1⊗ρ2)) is proposed as a correlation measure, but the paper does not establish that QR is nonnegative for correlated states or discuss conditions under which it could be negative. If QR can be negative, its interpretation as a measure of correlation needs qualification. Please clarify this point, either by proving nonnegativity for the states considered or by discussing the sign behavior.","section":"Section 3.2, Eq. (13)"}],"minor_comments":[{"comment":"The phrase 'the the factorizable ground state' contains a typo and should read 'the factorizable ground state'.","section":"Abstract"},{"comment":"The sentence 'We assume thst, the optimal decomposition...' contains a typo; it should read 'We assume that'.","section":"Appendix A"},{"comment":"In the argument for reducing a one-negative-coefficient decomposition, the stated RoM of the new decomposition is '2|1−µ| + 2ϵ', which appears to omit a factor of |a4|. The conclusion is unaffected, but the expression should be corrected for accuracy.","section":"Appendix A, case (ii)"},{"comment":"The function g(λ) is defined as 1 in the ordered phase and 0 in the disordered phase; explicitly writing g(λ) = Θ(λ − 1) would improve clarity.","section":"Section 2.1, Eq. (3)"},{"comment":"The sudden death temperature scaling Tc ∝ r^κ is presented with fitted κ values but no error bars; adding uncertainties would strengthen the quantitative claim.","section":"Section 4.2, Eq. (16)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a serious contribution with a clear main result and extensive numerical support. The primary issue is the incomplete proof of Eq. (8) in Appendix A; the gap is easily fixable by the shift argument described in the report. The numerical fitting concerns are secondary but should be addressed for the quantitative claims. Overall, the paper fits the journal's scope and is publishable after major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: this is a solid early map of magic (robustness of magic) in the transverse-field anisotropic XY chain. The central formula, Eq. (8), is right, and the factorizable ground state as a large-scale source of H magic states is a genuinely useful observation. The proof gap in Appendix A is real but easily fixed, so the main conclusions stand.\n\nWhat's new: the paper gives the first systematic study of RoM in an exactly solvable many-body model. It identifies a 'magic pseudocritical point' slightly inside the ordered phase where single-qubit magic turns on, with a linear ramp near that point and a derivative that scales as (λ−λ_c)^{β_x−1}, β_x=1/8. That exponent is checked against the known longitudinal magnetization exponent rather than extracted from the magic data, so the scaling law is not forced by fitting. The factorizable point (γ=1/3, λ≈1.06) yields a product state of H states, each with maximal single-qubit RoM √2−1, and the state is robust to a few percent mistuning. The two-qubit RoM persists over long distances, unlike bipartite entanglement, and the 'global magic' quantity Q_R peaks at the MPP for all distances studied. The thermal part shows sudden death and a crossover that can detect the quantum critical region via two distant qubits.\n\nSoft spots, in order. Appendix A assumes the optimal pseudo-mixture uses only the four stabilizer states in the xz-plane, without justifying the exclusion of |±i⟩. That is a real gap in the written proof, but it is patchable: if ⟨σ_y⟩=0, the weights on |+i⟩ and |−i⟩ are equal, and replacing that pair by half-weight on each of the other four states preserves the state and does not increase the L1 norm. So Eq. (8) survives; the appendix just needs the argument. Second, several fitted exponents (β_z, μ, ν, κ, and the tanh-fit constants) come with no error bars. The fits are plausible but not self-validating, and the quoted numbers should carry uncertainties. Third, the claim that Q_R peaks exactly at the MPP for arbitrary distance is supported by DMRG plus a derivative-kink argument, not an analytic proof, and 'exactly' is stronger than the numerics warrant. None of these are load-bearing flaws.\n\nBottom line: this is a paper for people working on magic state resources and on quantum-information probes of many-body criticality. The literature coverage is adequate, and the central claims look correct. It deserves a serious referee. If I were the editor, I'd send it out, asking the authors to fix the Appendix A argument and report uncertainties on fitted exponents.","headline":"Solid first map of magic in the XY chain: Eq. (8) is correct, the FGS H-state source is the main prize, but Appendix A needs a patch and fitted exponents need error bars.","tokens_in":20312,"tokens_out":5029,"would_cite":true,"duration_ms":48180,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In the XY spin chain, magic appears just past the critical point and reaches its ceiling at a factorizable point, where the ground state is a product of H states.","keywords":["robustness of magic","magic states","XY spin chain","quantum phase transition","factorizable ground state","H states","quantum critical region","thermal sudden death"],"falsifier":"Compute the RoM of a single-site reduced density matrix of the XY ground state near $\\lambda^*_c$ by solving the full linear program over all six single-qubit stabilizer states, with no restriction to the $xz$-plane. If for any $\\gamma$ and $\\lambda$ the full optimum is strictly below $\\max[\\langle\\sigma_x\\rangle+\\langle\\sigma_z\\rangle-1,0]$, then Eq. (8) is false and the location of the magic pseudocritical point and the extracted exponents change accordingly.","tokens_in":19151,"feed_emoji":"✨","tokens_out":17936,"duration_ms":170323,"temperature":0.7,"pith_summary":"This paper studies robustness of magic (RoM), the resource behind magic-state quantum computation, in the ground and thermal states of the transverse-field anisotropic XY spin chain. It establishes that in the symmetry-broken ground state, the single-qubit RoM is $\\max\\left[\\langle\\sigma_x\\rangle+\\langle\\sigma_z\\rangle-1,0\\right]$, so magic is absent throughout the disordered phase and turns on sharply at a magic pseudocritical point just inside the ordered phase, at $\\lambda\\approx1.00015$ for the Ising chain. At the factorizable point $\\gamma=1/3$, $\\lambda\\approx1.06$, the ground state becomes a product of $H$ states, each carrying the maximal single-qubit RoM $\\sqrt{2}-1$; this means a critical spin system can supply many unencoded magic states for injection schemes. The paper also shows that two-qubit magic survives over long distances, unlike bipartite entanglement, and that its temperature behaviour locates the quantum critical region. The result turns an operational resource for fault-tolerant computation into a long-range probe of critical many-body physics.","feed_headline":"A spin-chain point yields maximal magic states","feed_subtitle":"At γ=1/3, λ≈1.06 the ground state is a product of H states, the raw resource for fault-tolerant quantum computation.","key_machinery":"The load-bearing object is the robustness of magic (RoM), defined as the minimal $\\ell^1$ excess over a stabilizer-state pseudo-mixture: $R(\\rho)=\\inf\\{\\sum_k|X_k|-1 : AX=B\\}$ after fixing the Bloch-vector data. The paper's central device is a reduction valid for qubit states with $\\langle\\sigma_y\\rangle=0$: the optimal pseudo-mixture is assumed to involve only the four stabilizer states lying in the $xz$-plane, turning the RoM into the closed form $R_\\gamma(\\lambda)=\\max[\\langle\\sigma_x\\rangle+\\langle\\sigma_z\\rangle-1,0]$. That identity, combined with the exact Toeplitz-determinant solutions for the one- and two-point functions, converts the resource-theoretic problem into a one-dimensional scaling analysis. For two qubits the same reduced density matrix is fed into the RoM linear program, and the difference $Q_R=\\log(1+R(\\rho_{12}))-\\log(1+R(\\rho_1\\otimes\\rho_2))$ isolates the correlation-borne part of the magic.","core_discovery":"On its own terms, the paper's central claim is that the single-qubit robustness of magic in the symmetry-broken ground state of the XY chain is governed by $R_\\gamma(\\lambda)=\\max[\\langle\\sigma_x\\rangle+\\langle\\sigma_z\\rangle-1,0]$. Because $\\langle\\sigma_x\\rangle$ is the order parameter and vanishes in the disordered phase, all magic sits on the ordered side of the transition and starts at a pseudocritical point $\\lambda^*_c(\\gamma)$ that moves toward $\\lambda_c=1$ as $\\gamma\\to0$. Using the exact correlation functions, the paper locates the maximal single-qubit magic at $\\gamma_0=1/3$, $\\lambda_0\\approx1.06$, exactly the factorizable point, where the ground state is a product of $H$ states with RoM $\\sqrt{2}-1$. For two qubits it defines a global magic $Q_R=\\log(1+R(\\rho_{12}))-\\log(1+R(\\rho_1\\otimes\\rho_2))$, whose peak is at $\\lambda^*_c$ for every inter-site distance, and it shows that the two-qubit RoM persists over long distances. In the symmetry-unbroken thermal state, the first derivative of two-qubit RoM diverges logarithmically at criticality, the magic undergoes finite-temperature sudden death with $T_c\\propto r^\\kappa$, and the crossover at $T^*=a\\,T_{\\mathrm{cross}}$ delineates the quantum critical region.","pith_inferences":["The restricted-decomposition assumption behind the single-qubit formula is checkable with a small semidefinite program; if it holds, the same closed-form device may apply to any qubit-state model whose Bloch vector lies in a single plane, not only the XY chain.","Because the factorizable point supplies pure $H$ states without any entanglement, the mechanism gives a concrete preparation recipe: tune an interacting spin chain to its factorizing point to mass-produce magic-state ancillae, a purpose for which such points were previously avoided.","The distance-independent location of the global-magic maximum hints that $Q_R$ behaves like a connected correlation function of the order parameter; testing its scaling in finite-size systems could connect magic to standard critical exponents.","The finite-temperature sudden death of magic implies a hard temperature ceiling for using condensed-matter sources as magic-state factories, and the crossover scaling $T^*=a\\,T_{\\mathrm{cross}}$ gives a quantitative operating bound for such sources."],"forward_implications":["At the factorizable point $\\gamma_0=1/3$, $\\lambda_0\\approx1.06$, a spin chain in its symmetry-broken ground state yields an extensive number of unencoded $H$-type magic states, each with the maximal single-qubit RoM $\\sqrt{2}-1$, and small mistuning of the Hamiltonian costs less than about 0.1 percent in fidelity.","Magic does not flag the quantum critical point itself but a nearby magic pseudocritical point $\\lambda^*_c$ slightly inside the ordered phase; the offset falls with anisotropy roughly as $\\gamma^{5.55}$, so the effect is most pronounced near the Ising limit.","The correlation-borne part of two-qubit magic, $Q_R$, peaks exactly at $\\lambda^*_c$ for every inter-qubit distance and decays only slowly with distance, providing a long-range quantum-correlation probe where bipartite entanglement dies after the second neighbour.","In the symmetry-unbroken thermal state, the first derivative of two-qubit RoM diverges logarithmically at criticality, two-qubit magic suffers sudden death at finite temperature with $T_c\\propto r^\\kappa$, and derivatives of the RoM collapse as functions of $T/T_{\\mathrm{cross}}$, tracing out the quantum critical region."],"supporting_citations":[{"why":"Provides the exact thermodynamic-limit magnetization and correlation functions used to build the reduced density matrices.","marker":"[60]"},{"why":"Supplies the Toeplitz-determinant form of the two-point correlation functions used in Eq. (2).","marker":"[61]"},{"why":"Gives the algebraic scaling of the order parameter near criticality that underlies the RoM scaling laws.","marker":"[14]"},{"why":"Establishes the linearized minimization formulation of RoM that the paper optimizes and reduces.","marker":"[36]"},{"why":"Sets the single-qubit equatorial-plane maximum for H states and the two-qubit maximum used as benchmarks.","marker":"[37]"},{"why":"Specifies the single-qubit stabilizer octahedron and the equatorial-plane bound used to identify H states as maximal.","marker":"[65]"},{"why":"Identifies the factorizable ground-state point at which the paper finds maximal single-qubit magic.","marker":"[24]"},{"why":"Provides the quantum-critical scaling and crossover-temperature ansatz used to map the finite-temperature critical region.","marker":"[80]"}],"fun_headline_variants":["Magic peaks at factorizable point in XY chain","Critical spin chain: maximal magic at distant qubits","Magic in correlations spikes near phase transition","XY model ground state is a sea of magical H states","Two-qubit magic survives long distances near criticality"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that for a qubit with $\\langle\\sigma_y\\rangle=0$, the optimal RoM pseudo-mixture never needs the $\\left|+i\\right\\rangle$ and $\\left|-i\\right\\rangle$ stabilizer states; if including them lowered the $\\ell^1$ norm, the closed-form RoM and every single-qubit scaling law built on it would have to be revised.","fun_headline_variants_meta":{"raw":{"variants":["Magic peaks at factorizable point in XY chain","Critical spin chain: maximal magic at distant qubits","Magic in correlations spikes near phase transition","XY model ground state is a sea of magical H states","Two-qubit magic survives long distances near criticality"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000596,"raw_usage":{"total_tokens":2830,"prompt_tokens":1024,"completion_tokens":1806,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":640,"completion_tokens_details":{"reasoning_tokens":1733}},"tokens_in":640,"tokens_out":1806,"duration_ms":481964,"temperature":1.0,"reasoning_tokens":1733,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:53:14.623911+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the RoM of a single-site reduced density matrix of the XY ground state near $\\lambda^*_c$ by solving the full linear program over all six single-qubit stabilizer states, with no restriction to the $xz$-plane. If for any $\\gamma$ and $\\lambda$ the full optimum is strictly below $\\max[\\langle\\sigma_x\\rangle+\\langle\\sigma_z\\rangle-1,0]$, then Eq. (8) is false and the location of the magic pseudocritical point and the extracted exponents change accordingly.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the exact thermodynamic-limit magnetization and correlation functions used to build the reduced density matrices."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Toeplitz-determinant form of the two-point correlation functions used in Eq. (2)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the algebraic scaling of the order parameter near criticality that underlies the RoM scaling laws."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the linearized minimization formulation of RoM that the paper optimizes and reduces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Sets the single-qubit equatorial-plane maximum for H states and the two-qubit maximum used as benchmarks."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Specifies the single-qubit stabilizer octahedron and the equatorial-plane bound used to identify H states as maximal."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies the factorizable ground-state point at which the paper finds maximal single-qubit magic."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the quantum-critical scaling and crossover-temperature ansatz used to map the finite-temperature critical region."}],"review_version":1}