{"id":"4f53fbe0-75ce-44c2-87b3-11abc3689417","arxiv_id":"2502.07884","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The type of magnetic order induced in local moments by a Kitaev-type spin liquid is controlled by whether the spin operators commute with the spin liquid's flux operators, yielding short-range Kitaev, long-range RKKY-like, or dimerized phases.","lead":"This paper asks what happens when a layer of ordinary magnetic atoms is placed next to a quantum spin liquid, a state of matter with no ordinary magnetism. Using three exactly solvable models, the authors show that the spin liquid's ghostly Majorana particles can make the magnetic atoms order in different patterns, from short-range spin-liquid order to long-range antiferromagnetism.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Two-impurity-derived Jeff is assumed pairwise additive and unmodified by the order it induces; without a multi-spin/self-consistency check, the YL AFM and partial-gapping claims are not fully established.","rationale":"The reader's weakest assumption identifies pairwise additivity and 16-site finite-size effects; I agree partially. The pairwise-additivity concern is the most load-bearing because it underpins both the Yao-Lee and square-lattice central claims, whereas the 16-site issue only affects the square-lattice phase boundaries and is less likely to change the deep-DQP conclusion at Jw = 0. The self-consistency/back-action point is an under-appreciated part of the same assumption: the induced order modifies the mediating QSL, so the effective exchange kernel should in principle be recomputed in the ordered state. The paper has real strengths: exact solvability of the QSL models, a benchmark of Kitaev perturbation against exact diagonalization, and a clean fit to RKKY forms. However, the missing many-spin check means the conditional verdict is appropriate. I recommend keeping the reader's CONDITIONAL verdict; the proposed concrete test would either support or falsify the central inference. No ad hominem is intended; this is purely a question about the argument's quantitative closure.","tokens_in":21743,"tokens_out":18983,"duration_ms":184951,"concrete_test":"Perform a self-consistent classical-spin ground-state search on a finite Yao-Lee cluster (e.g., 6x6 with a few local moments): treat each local moment as a classical unit vector, diagonalize the full Majorana Hamiltonian including JK Si·sigma_i for the current spin configuration, compute the total energy, and iterate by minimizing over spin orientations. Compare the converged ground-state spin pattern and the Majorana spectrum with the Heisenberg AFM predicted from Heff = Sum Jeff Si·Sj using the paper's Jeff values. If the converged pattern differs, or if the energy curvature for a pair of spins in the converged background deviates from the two-impurity Jeff by more than 20%, the pairwise-additive assumption is falsified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central mechanism rests on converting two-impurity energy differences into a many-spin effective Hamiltonian Heff = Sum_ij Jeff Si·Sj (Eq. 14), but this step is uncontrolled. In the Yao-Lee model, Jeff is extracted from exact diagonalization of the Majorana Hamiltonian with two fixed impurity spins; nothing rules out multi-spin interactions of comparable size once all local moments are present, and no small-parameter estimate for their strength is given. Moreover, the AFM order predicted from Heff produces a staggered field that gaps two of the three Majorana flavors; this gap suppresses the gapless spin susceptibility responsible for the 1/R^3 exchange, so the bare two-impurity Jeff is not obviously the correct kernel in the ordered state. The back-action paragraph in Sec. III.B is a one-way mean-field statement, not a self-consistent solution. The square-lattice phase diagram (Fig. 7) adds a finite-size layer: the DQP-AFM boundaries at Jxy/Jz ≈ 1.5/2.4 come from a single 16-site cluster with truncated long-range Jz, so the quoted Jw interval [-2.81, 0.79] is not established in the thermodynamic limit. The qualitative classification by flux commutation is plausible and well motivated; the weak link is the quantitative inference to the many-spin ground state.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies how local moments coupled by a Kondo term to three exactly solvable Kitaev-type spin liquids acquire effective exchange interactions. For the honeycomb Kitaev model, second-order perturbation theory gives a short-range Kitaev Hamiltonian for the local moments, and higher-order terms add plaquette interactions while preserving a gapless QSL. For the Yao-Lee model, exact diagonalization of two impurity spins in the Majorana representation gives a long-range 1/R^3 exchange with graphene-like oscillations; the authors conclude that the many-spin ground state is a Heisenberg antiferromagnet that partially gaps two of the three Majorana flavors. For a square-lattice gamma-matrix model, a short-range XY bond coexists with a long-range Ising 1/R^3 interaction; exact diagonalization of the effective model on 16 sites yields a dimerized quantum paramagnet in a window of the flux-tuning parameter Jw and an Ising antiferromagnet outside it. The paper's central organizing principle is that the range and sign structure of the induced exchange are controlled by whether the spin operators commute with the flux operators of the spin liquid.","tokens_in":1867,"tokens_out":2021,"duration_ms":114284,"significance":"If the central claims hold, the paper identifies a useful organizing principle: flux commutation distinguishes short-range Kitaev-like exchange from long-range Majorana-mediated RKKY exchange. The perturbative treatment of the Kitaev model is careful and benchmarked against 16-site exact diagonalization, and the idea that Yao-Lee Majorana fermions can mediate a graphene-like 1/R^3 interaction is falsifiable and potentially relevant to heterostructures. However, the decisive steps from two-impurity energy differences to many-spin order, and from the bare exchange kernel to the ordered state, are not fully controlled. As it stands, the paper is a promising proposal with strong numerical evidence for the two-impurity couplings, but the quantitative phase diagram and the partial-gapping claim require additional self-consistency and finite-size analysis.","major_comments":[{"comment":"The central step from two-impurity energy differences to the many-spin Hamiltonian H_eff = sum_ij J_eff(R_i-R_j) S_i dot S_j assumes pairwise additivity. The paper computes J_eff from exact diagonalization of the Majorana Hamiltonian with exactly two impurity spins (Fig. 4) and never checks whether three- and four-spin terms generated by the Kondo coupling are negligible when all local moments are present. Since J_eff itself is approximately two orders of magnitude weaker than the Kitaev short-range exchange, there is no explicit small parameter that guarantees multi-spin terms are subdominant. A three-impurity calculation, or a controlled estimate of the next-order multi-spin corrections, is needed before the Heisenberg-antiferromagnet ground-state claim can be accepted.","section":"Sec. III.B, Eq. (14)"},{"comment":"The claim that the antiferromagnetic state partially gaps two of the three Majorana flavors is not self-consistent with the RKKY kernel used to predict that state. The 1/R^3 graphene form of Eq. (13) is the static susceptibility of gapless Dirac fermions; once the staggered field J_K(-1)^n sigma^z opens a gap in two of the Majorana flavors, the same susceptibility is suppressed at long distances, so the bare two-impurity J_eff is not obviously the correct interaction kernel in the ordered phase. The paragraph asserts a one-way back action but does not solve the coupled order-parameter/gap problem. A self-consistent mean-field treatment, or an estimate showing that the induced gap is negligible on the relevant length scales, is required to support the AFM and partial-gapping claims.","section":"Sec. III.B, final paragraph"},{"comment":"The DQP-to-Ising-AFM boundaries at Jxy/Jz of approximately 1.5 (pi-flux) and 2.4 (0-flux) are obtained from a single 16-site cluster with the long-range Jz truncated at fourth-nearest-neighbor distances, and the quoted Jw interval [-2.81, 0.79] in Table I inherits these values without any estimate of finite-size or truncation error. Critical ratios in frustrated quantum paramagnets can shift significantly with cluster size, and the thermodynamic-limit phase diagram is not established. The authors should either provide finite-size scaling on larger clusters or state explicitly that the interval is a cluster-size-dependent estimate rather than a thermodynamic-limit result.","section":"Sec. III.C.3, Fig. 7 and Table I"},{"comment":"The identification of the Yao-Lee J_eff with the charge-neutrality graphene RKKY form, Eq. (13), is not fully derived. Appendix C maps the z-component Kondo problem to the two graphene-like Hamiltonians H_p and H_h of Eqs. (C2)-(C3), which contain chemical-potential terms of +/-J_K, whereas Eq. (13) is the mu equals 0 result of Ref. [38]. The paper neither shows that the J_K chemical potential is negligible for the fitted distances nor derives the 1/R^3 form directly from the mapped Hamiltonian; it fits the overall prefactor C. The prefactor is therefore a fit parameter, not a prediction, and the statement that the YL result is identical to graphene RKKY up to prefactors is stronger than the presented derivation supports.","section":"Sec. III.B and Appendix C"}],"minor_comments":[{"comment":"The text after Eq. (12) defines k1 = k dot b1 and then k1 = k dot b2; the second should read k2 = k dot b2.","section":"Eq. (12) and Appendix B"},{"comment":"The sentence 'The SU(2) symmetry in Heff energes from...' contains a typo; 'energes' should be 'emerges'.","section":"Sec. III.B"},{"comment":"The phrase 'even tough the z-component has long-ranged unfrustrated interactions' contains a typo; 'tough' should be 'though'.","section":"Sec. III.C.1"},{"comment":"The Jw interval in Table I should specify its units (presumably |K|) and use consistent notation for 'not in' rather than the TeX fragment shown in the table.","section":"Table I"},{"comment":"The main text states that at the fixed value Jxy/Jz = 44 the dimer overlap is 0.9997, while the Fig. 7 caption quotes Jxy/Jz = 17 for |<psi_D|psi>|^2 >= 0.99; the relation between these two numbers should be clarified.","section":"Fig. 7 caption and Sec. III.C.3"},{"comment":"The relabeling of the b^4 and b^5 Majorana operators into the c^z and c^x flavors is described only in words; a small table or explicit substitution would improve readability.","section":"Sec. II, Eqs. (6)-(8)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope and is likely to interest the spin-liquid and heterostructure communities. My main concern is not the two-impurity numerics, which appear sound, but the uncontrolled leap from those two-impurity couplings to the many-spin ground state, compounded by the one-way back-action argument and the single-cluster phase diagram. These issues are fixable in revision, so I do not recommend rejection. The use of a fitted prefactor C should also be clearly labeled as a fit rather than a parameter-free prediction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a genuinely useful addition to the Kitaev-spin-liquid literature. The new results are the Yao-Lee and square-lattice calculations: the Yao-Lee model gives a graphene-like 1/R^3 RKKY interaction that favors Heisenberg antiferromagnetism, and the square-lattice model produces a dimerized quantum paramagnet or an Ising antiferromagnet depending on a flux-gap tuning term. The Kitaev second-order exchange was already in Ref. [16], but the higher-order terms and the overall commutation-based classification (spin operators commuting with flux operators => long range; anticommuting => short range) are new and well motivated. The Kitaev perturbative calculation is clean, benchmarked against 16-site exact diagonalization, and the mapping of the Yao-Lee model to a graphene-like complex-fermion problem is a nice step. The soft spots are real, and the stress-test note hits the main one: the effective Hamiltonian Heff = sum_{ij} Jeff S_i.S_j is extracted from two-impurity energy differences, and pairwise additivity is assumed without a control parameter. In the Yao-Lee case the induced order gaps two Majorana flavors, which will feed back on the exchange kernel, and the paper's back-action discussion is a one-way mean-field statement rather than a self-consistent solution. For the square-lattice model the phase boundaries come from a single 16-site cluster, so the quoted J_w interval [-2.81, 0.79] is hostage to finite-size effects. These issues do not kill the qualitative message--the commutation principle is compelling and likely correct--but they do undercut the quantitative claims about partial gapping and the precise DQP/AFM boundary. The fitted prefactor C in the Yao-Lee RKKY form is not derived, but that is a minor caveat since it only sets the overall scale and does not affect the sign or phase. Who is this for? Anyone working on magnetic impurities or heterostructures involving Kitaev-type spin liquids will get value from it, and the classification table is a nice summary. It deserves a serious referee; I would engage with it, but the referee should push for either a multi-spin check (e.g., classical or mean-field treatment of the full coupled system) or a clear statement that the quantitative phase boundaries are provisional. My recommendation: send it to peer review, with revision.","headline":"Solid theory paper with a compelling organizing principle, but the step from two-impurity exchange to many-spin order needs work before the quantitative claims can be trusted.","tokens_in":671,"tokens_out":1352,"would_cite":true,"duration_ms":32632,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The commutation of spins with fluxes decides which magnetic order a spin liquid induces.","keywords":["quantum spin liquids","Kitaev model","Yao-Lee model","Majorana fermions","RKKY interaction","Kondo lattice","magnetic order","dimerized quantum paramagnet"],"falsifier":"The central claim would be settled by solving the full Kondo-lattice Hamiltonian with all local moments present, not just two impurities, in the perturbative regime $J_K/\\Delta\\ll 1$: if the ground state of the full system disagrees with the ground state of the pairwise-derived Hamiltonian (14) or (16), the pairwise-additivity assumption fails. For the square-lattice model specifically, exact diagonalization of Hamiltonian (16) on clusters larger than 16 sites would show whether the transition thresholds $J_{xy}/J_z\\approx 1.5$ ($\\pi$-flux) and $\\approx 2.4$ (0-flux) persist; a shift larger than the finite-size scatter would refute the phase boundary.","tokens_in":21506,"feed_emoji":"🧲","tokens_out":8581,"duration_ms":67050,"temperature":0.7,"pith_summary":"The paper asks whether local magnetic moments can order when they interact only through a nearby Kitaev-type spin liquid, with no direct exchange among themselves. It claims that the decisive property is how each spin component of the spin liquid commutes with the conserved flux operators. A single Majorana flavor (Kitaev) produces only short-range induced Kitaev couplings and keeps the spin liquid intact; three flavors on the honeycomb (Yao-Lee) produce a long-range $1/R^3$ Heisenberg antiferromagnet that gaps two of the three Majorana flavors; the two-flavor square-lattice model produces mostly a dimerized quantum paramagnet, tunable into an Ising antiferromagnet by a flux-gap term. A sympathetic reader would care because the result sketches a design rule: which fractionalized excitations mediate which magnetic order, and how to choose a spin liquid to engineer a desired ordered phase.","feed_headline":"Magnetic order follows from spin-flux commutation in Kitaev liquids","feed_subtitle":"Yao-Lee liquid makes a 1/R^3 antiferromagnet; the square-lattice liquid mostly dimerizes instead.","key_machinery":"The central object is the commutation relation between the spin operator $\\sigma_i^\\alpha$ on the Kondo-coupled site and the conserved plaquette flux operators $W_p$ of the Kitaev-type spin liquid. In the Kitaev model $\\{\\sigma_i^\\alpha, u_{ij}^\\alpha\\}=0$, so a spin flip creates two visons and the induced exchange is short-ranged; in the Yao-Lee model $[\\sigma_i^\\alpha, W_p]=0$, so no flux is created and the itinerant Majorana fermions mediate long-range RKKY exchange; in the square-lattice model $\\sigma^z$ commutes with $W_p$ while $\\sigma^x,\\sigma^y$ anticommute, producing the mixed short-range XY plus long-range Ising Hamiltonian. The second object that carries the argument is the two-impurity energy difference $J_{\\mathrm{eff}}=(E_{\\mathrm{FM}}-E_{\\mathrm{AFM}})/2$ computed from aligned versus anti-aligned local moments in the Majorana representation, fitted to the graphene-neutrality form $J_{\\mathrm{eff}}(R)\\propto [1+\\cos(2k_D\\cdot R)]/|R|^3$.","core_discovery":"On its own terms, the central discovery is that the range and sign structure of the induced exchange are controlled by the commutation of the spin operators with the flux (Wilson-loop) operators of the host spin liquid. When spin operators anticommute with flux operators, each Kondo term creates vison excitations, so the lowest-order induced interaction is short-ranged: for the Kitaev model it remains a Kitaev Hamiltonian with renormalized couplings, and the spin-liquid ground state survives higher-order terms. When spin operators commute with flux, the Majorana fermions act as the carriers and the induced interaction is RKKY-like: for the Yao-Lee model it is $J_{\\mathrm{eff}}(R)\\propto (J_K^2/K)\\,[1+\\cos(2k_D\\cdot R)]/|R|^3$, ferromagnetic within a sublattice and antiferromagnetic between sublattices, so the local moments form a Heisenberg antiferromagnet and back-react on the spin liquid by gapping two of the three Majorana flavors. The square-lattice model mixes both behaviors: strong short-range XY couplings on one bond type together with long-range Ising interactions, giving a dimerized quantum paramagnet for $J_{xy}/J_z^{\\max}\\gtrsim 1.5$ ($\\pi$-flux) or $\\gtrsim 2.4$ (0-flux), with an Ising antiferromagnet outside a window set by the added flux-gap term $J_w$.","pith_inferences":["A testable extension of the commutation-rule criterion is a four-flavor Kitaev-type model on the honeycomb lattice: the one-, two-, and three-flavor comparison in the paper suggests the induced exchange becomes richer, but the criterion alone does not predict which ground state wins.","An experimental consequence, if the Yao-Lee result is realized in a heterostructure with a spin-orbital liquid, is an antiferromagnetic response decaying as $1/R^3$ with Dirac-point oscillations, distinguishable from metallic RKKY by the absence of Fermi-surface nesting features and by the partial gapping of the Majorana spectrum.","A caution: the square-lattice phase diagram rests on pairwise $J_{\\mathrm{eff}}$ values and a 16-site diagonalization, so larger clusters could shift the $J_{xy}/J_z$ thresholds and the $J_w$ window for the dimerized paramagnet."],"forward_implications":["If the central claim is right, coupling a Mott-insulating layer to a Yao-Lee spin liquid produces a Heisenberg antiferromagnet on the local moments, and the ordered moments partially gap the Majorana spectrum of the spin liquid.","For the Kitaev model, the local-moment layer remains a Kitaev spin liquid with renormalized couplings; higher-order terms preserve exact solvability and only renormalize the Fermi velocity without opening a gap.","For the square-lattice model, the default induced Hamiltonian sits deep in the dimerized quantum paramagnet phase ($J_{xy}/J_z\\approx 44$ at $J_w=0$), with nearly perfect overlap with the dimer product state.","Adding a flux-gap term $J_w$ tunes the ratio $J_{xy}/J_z$ through a $\\pi$-flux to 0-flux transition, so the same effective model can be switched between a dimerized paramagnet and an Ising antiferromagnet by changing flux-sector energetics.","Because both $J_{xy}$ and $J_z$ scale with $J_K^2/K$, their ratio is independent of the Kondo coupling at small $J_K$, making the phase assignment stable against the precise coupling strength."],"supporting_citations":[{"why":"Supplies the exactly solvable honeycomb Kitaev model and its Majorana representation, the base case of the study.","marker":"[21]"},{"why":"Defines the Yao-Lee model, the three-flavor spin-orbital liquid whose commuting spins give long-range RKKY exchange.","marker":"[22]"},{"why":"Defines the square-lattice Gamma-matrix generalization with two Majorana flavors and a pi-flux ground state.","marker":"[23]"},{"why":"Provides the graphene-at-charge-neutrality RKKY formula used to fit the induced exchange decay and oscillations.","marker":"[38]"},{"why":"Fixes the ground-state flux sectors (0-flux for honeycomb, pi-flux for square lattice) on which the flux-commutation argument depends.","marker":"[37]"},{"why":"Gives the exact nearest-neighbor spin correlation used for the second-order Kitaev and square-lattice couplings.","marker":"[39]"},{"why":"Previous result for the effective exchange of magnetic impurities in the Kitaev spin liquid that the paper benchmarks its perturbative Kitaev coupling against.","marker":"[16]"}],"fun_headline_variants":["Spin-flux commutation decides magnetic order in Kitaev liquids","Kitaev liquids order magnetically when spins commute with flux","Flux commutation governs mediated magnetic order in Kitaev models","From Kitaev short-range to Yao-Lee long-range order via flux commutation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes that the exchange between two local moments, extracted from the aligned versus anti-aligned energy difference, is pairwise additive and determines the many-moment ground state, without solving the full system of all coupled moments together; for the square-lattice model it also assumes the 16-site phase boundaries survive in the thermodynamic limit.","fun_headline_variants_meta":{"raw":{"variants":["Spin-flux commutation decides magnetic order in Kitaev liquids","Kitaev liquids order magnetically when spins commute with flux","Flux commutation governs mediated magnetic order in Kitaev models","From Kitaev short-range to Yao-Lee long-range order via flux commutation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000786,"raw_usage":{"total_tokens":3543,"prompt_tokens":1093,"completion_tokens":2450,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":709,"completion_tokens_details":{"reasoning_tokens":2376}},"tokens_in":709,"tokens_out":2450,"duration_ms":17940,"temperature":1.0,"reasoning_tokens":2376,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T11:31:45.860693+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"The central claim would be settled by solving the full Kondo-lattice Hamiltonian with all local moments present, not just two impurities, in the perturbative regime $J_K/\\Delta\\ll 1$: if the ground state of the full system disagrees with the ground state of the pairwise-derived Hamiltonian (14) or (16), the pairwise-additivity assumption fails. For the square-lattice model specifically, exact diagonalization of Hamiltonian (16) on clusters larger than 16 sites would show whether the transition thresholds $J_{xy}/J_z\\approx 1.5$ ($\\pi$-flux) and $\\approx 2.4$ (0-flux) persist; a shift larger than the finite-size scatter would refute the phase boundary.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the exactly solvable honeycomb Kitaev model and its Majorana representation, the base case of the study."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Yao-Lee model, the three-flavor spin-orbital liquid whose commuting spins give long-range RKKY exchange."},{"cited_title":"Zheng and A","cited_arxiv_id":null,"evidence_quote":"Defines the square-lattice Gamma-matrix generalization with two Majorana flavors and a pi-flux ground state."},{"cited_title":"Poliakov, W.-H","cited_arxiv_id":null,"evidence_quote":"Provides the graphene-at-charge-neutrality RKKY formula used to fit the induced exchange decay and oscillations."},{"cited_title":"Microscopic Roadmap to a Yao-Lee Spin-Orbital Liquid","cited_arxiv_id":"2410.21389","evidence_quote":"Gives the exact nearest-neighbor spin correlation used for the second-order Kitaev and square-lattice couplings."}],"review_version":1}