{"id":"009ab42f-03d0-4860-9e17-e3e14c630886","arxiv_id":"1908.04554","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":6,"one_line_summary":"RKKY spin coupling between opposite surfaces of a Weyl semimetal slab is strong, non-monotonic in thickness, and peaks at a critical thickness set by the Weyl node separation.","lead":"This paper studies how magnetic impurity atoms on the two surfaces of a thin Weyl semimetal film influence each other through conduction electrons. The coupling is unexpectedly strong and strongest at a particular film thickness, which may guide experiments on magnetic ordering in thin topological films.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Optimum-thickness scaling tested at a single lateral separation; peak shifts with R, so Lc ~ 1/k0 may not be intrinsic.","rationale":"The reader's weakest assumption, representativeness of the minimal four-node model for real materials, is a legitimate external caveat but not an internal flaw. The stress-test pass instead identifies a more concrete internal gap: the paper's strongest quantitative claim about the thin-film limit, the Lc ~ 1/k0 scaling, is verified at a single lateral separation R = 40a, while the authors themselves show that the position of the peak shifts with R (Fig. 5, bottom). Because R·k0 varies substantially across the three values of m/λ used in Appendix B, the observed collapse could be coincidental. A complementary gap is that the non-monotonic thickness dependence is only shown at µ = 0; at finite chemical potential, bulk Fermi-surface contributions could mask or shift the maximum, and the paper does not provide a thickness scan at µ ≠ 0. These concerns do not undermine the qualitative existence of a maximum within the model, nor the internal consistency of the calculations, but they do affect how much weight can be placed on the specific 'thin-film limit' and the TaAs thickness estimate. A conditional accept is therefore appropriate: the central physics is credible, but the quantitative scaling claim should be either backed by multi-R and finite-µ scans or stated with the appropriate caveats.","tokens_in":13787,"tokens_out":47505,"duration_ms":478042,"concrete_test":"Recompute the opposite-surface Jij thickness scan for m/λ = 0.2, 0.5, and 0.95 at lateral separations R = 20a, 40a, and 60a, and also at a finite chemical potential such as µ = 0.1λ, using the same recursive Green's function method (or the analytical n = 0 wavefunctions). For each case, extract the critical thickness Lc at which the coupling is maximal. If Lc·k0 is not constant across R and µ, the claimed 1/k0 thin-film scaling is a fixed-R, µ = 0 artifact and the paper should present Lc as a function of both R and k0 rather than collapsing it onto k0 alone.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Appendix B verifies Lc·k0 ≈ constant using only R = 40a for all three values of m/λ (Fig. 8). Yet the paper's own Fig. 5 (bottom) shows that the peak of the opposite-surface RKKY coupling shifts with lateral separation R, so Lc depends on at least two length scales: R and 1/k0. Across the three k0 values tested, R·k0 changes by roughly a factor of four (about 55 for m/λ = 0.2 versus about 13 for m/λ = 0.95), so the apparent collapse of Lc onto 1/k0 could be an artifact of the single fixed R rather than evidence that k0 is the only relevant scale. If Lc actually depends on R (or on the dimensionless product R·k0), then the 'thin-film limit' is not a well-defined material property, and the estimated thickness for real systems such as TaAs (several tens to a hundred lattice spacings) is not supported. The qualitative non-monotonic thickness dependence is not in question; the load-bearing quantitative part is the claimed intrinsic 1/k0 scaling and the associated material prediction. In addition, the thickness dependence is only demonstrated at µ = 0; the finite-µ case (Fig. 7) is shown at a single thickness, so the persistence of the optimum at realistic doping is also untested.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies RKKY interactions between impurity spins placed on the surfaces of a thin-film Weyl semimetal. The authors use a minimal four-node time-reversal-symmetric tight-binding model and a low-energy Dirac model with infinite-mass boundary conditions. They show that spins on the same surface couple through Fermi arc states with pronounced anisotropy, while spins on opposite surfaces couple through both Fermi arcs and bulk states. The central claim is that the opposite-surface coupling is unexpectedly strong and non-monotonic in film thickness, with a maximum at an optimum thickness. Results are obtained both from analytical wavefunctions of the slab and from a recursive Green's function computation on the tight-binding model, and the authors further claim that the critical thickness scales approximately as 1/k0, where k0 is the separation of the Weyl nodes.","tokens_in":14131,"tokens_out":25282,"duration_ms":259619,"significance":"If the central claims hold, the paper identifies a concrete thin-film regime in which intersurface RKKY coupling is strong and controllable by thickness, with possible consequences for magnetic ordering on WSM surfaces and for the interpretation of future experiments. The paper has real strengths: it provides two independent calculational methods that agree qualitatively, the recursive Green's function scheme is efficient and well suited to the two-surface problem, and the symmetry analysis in Appendix A is a useful reference for the structure of the spin susceptibility matrix. The analytical derivation in Sec. II is self-contained, and the comparison between the low-energy wavefunction approach and the tight-binding numerics gives confidence in the qualitative behavior. The main quantitative claim, however, namely the Lc ~ 1/k0 scaling and the associated material estimate for TaAs, rests on a single lateral separation and therefore needs additional support before it can be regarded as established.","major_comments":[{"comment":"The claim that the critical thickness scales as Lc ~ 1/k0 is tested for only one lateral separation, R = 40a, in Fig. 8, while Fig. 5 (bottom) shows that the peak thickness shifts with R. Since R·k0 changes by roughly a factor of four across the three values of m/λ studied (about 55 for m/λ = 0.2 versus about 13 for m/λ = 0.95), the apparent collapse of Lc·k0 onto a constant could be a coincidental consequence of the fixed R rather than evidence for an intrinsic 1/k0 scale. Please provide a two-parameter study varying both R and k0, or at least show that the extracted Lc is independent of R at fixed k0, before using this scaling to define the thin-film limit and to estimate the critical thickness for TaAs in Section V.","section":"Appendix B, Fig. 8; Fig. 5 (bottom)"},{"comment":"The non-monotonic thickness dependence and the existence of an optimum thickness are demonstrated only at µ = 0. The finite-µ results in Fig. 7 are shown for a single thickness (Nz = 33), so the statement in Section V that the envelopes 'behave rather similarly' to the µ = 0 case is not supported by a thickness scan. Since the abstract states the thickness maximum without a µ qualifier, and since real materials may not be at the nodal energy, the persistence of the optimum thickness at finite doping is an open point. Either a finite-µ thickness sweep or an explicit qualification of the claim is needed.","section":"Section V and Fig. 7"},{"comment":"The quantitative material prediction for TaAs in Section V assumes that the minimal four-node model with straight, spin-polarized Fermi arcs captures the essential intersurface physics of a real WSM. This is a reasonable starting point, but it is a nontrivial assumption: real Fermi arcs can have different connectivity, curvature, and spin texture, and the penetration behavior near the Weyl node projections may differ. Because the TaAs thickness estimate is the most actionable claim, the authors should either soften it to an explicit model-based estimate or test the sensitivity of Lc to the arc structure, for example by using a model with curved arcs or a TaAs-derived surface spectrum.","section":"Eq. (1) and Section V"}],"minor_comments":[{"comment":"The bottom-right entry of the 2x2 matrix G' should be G_{N_z N_z}, not G_{1 N_z}; as written the matrix is inconsistent with the stated goal of computing both surface-to-surface propagators.","section":"Eq. (24)"},{"comment":"The text contains an unresolved placeholder '[REF]' in the discussion of boundary conditions for the Dirac equation; a proper citation is needed.","section":"Section II.B"},{"comment":"The phrase 'as well using a two-surface recursive Green's function analysis' should read 'as well as using'.","section":"Abstract"},{"comment":"Reference 15 lists the year as '20015'; it should be 2015.","section":"References"},{"comment":"The caption says 'The energy values decreases exponentially'; this should be 'decrease'.","section":"Fig. 2 caption"},{"comment":"The argument of the cosine contains 'π y' with y appearing to be a dimensioned coordinate; the notation should specify that lattice spacing a is set to unity or otherwise define the dimensionless combination.","section":"Eq. (35)"}],"recommendation":"major_revision","confidential_remarks":"This is a solid paper and the qualitative central claim is convincing. The main obstacle is the insufficiently supported quantitative scaling Lc ~ 1/k0, which is fixable with additional computations and does not require a change of scope. I would not push the model-realism concern beyond what the authors already acknowledge; a clearly stated caveat around the TaAs estimate should be sufficient on that point."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a solid theory paper worth a serious referee. The genuinely new thing is the opposite-surface RKKY coupling in a finite-thickness WSM slab, and the finding that it is non-monotonic in thickness with a maximum. The two-surface recursion plus the analytical slab wavefunctions agree, and the symmetry structure in Appendix A is a nice check.\n\nWhat it does well: the authors extend the semi-infinite single-surface treatments (Refs. 19-21) to a slab, and they identify a qualitative effect that is absent there: the intersurface coupling mediated by Fermi arcs and bulk states. The same-surface Jxx that decays with thickness is also a plausible surface-related effect, and they show it disappears when the surfaces are perpendicular to the node separation. That control calculation is good practice.\n\nWhere I'd push back: the claimed Lc ~ 1/k0 scaling is weaker than the paper suggests. Appendix B tests it at a single lateral separation R = 40a for all three k0 values, but Fig. 5 (bottom) shows the peak position shifts with R. So R*k0 changes by roughly a factor of four across the three m/λ values, and the collapse onto 1/k0 could be accidental. To call 1/k0 the material-relevant scale, they need at least one additional R for each k0, or an argument for why the R dependence drops out. This doesn't kill the qualitative non-monotonic result, but it does mean the 'thin-film limit' defined by Lc is not yet established as an intrinsic property. The finite-µ results are similarly thin: Fig. 7 shows 2kF oscillations at one thickness, so the optimum thickness at realistic doping is untested.\n\nThe minimal four-node model is a standard caveat, but the TaAs estimate (tens to a hundred lattice spacings) depends on the straight-arc model; real Fermi-arc connectivity could shift it. I would not call that a flaw, just a limit.\n\nBottom line: the paper deserves peer review. I would recommend 'accept with revisions' conditional on the authors either testing Lc at a second R or softening the scaling claim to 'consistent with 1/k0.' The central physics of inter-surface coupling is solid.","headline":"Solid theory paper on RKKY in WSM slabs with a genuine new intersurface effect; the Lc ~ 1/k0 scaling is only single-separation data and needs a second test or softer wording.","tokens_in":14615,"tokens_out":1959,"would_cite":true,"duration_ms":20756,"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":"The paper claims that impurity spins on opposite surfaces of a Weyl-semimetal thin film couple strongly through Fermi arcs and bulk states, with the coupling peaking at an optimum thickness set by the Weyl-node separation.","keywords":["RKKY interaction","Weyl semimetal thin films","Fermi arcs","spin-spin coupling","recursive Green's function","thin-film limit","Weyl nodes","magnetic impurities"],"falsifier":"Measure the thickness dependence of the magnetic correlation between impurities on opposite surfaces of a Weyl-semimetal slab, for example in TaAs thin films at low temperature; the claim fails if the coupling decays monotonically with thickness without any maximum, or if the opposite-surface coupling is orders of magnitude weaker than same-surface coupling at all thicknesses. A less demanding check is to compute the same RKKY tensor in a material-specific model with curved or differently connected Fermi arcs and see whether the non-monotonic peak survives.","tokens_in":13623,"feed_emoji":"🧲","tokens_out":6427,"duration_ms":67001,"temperature":0.7,"pith_summary":"This paper asks how magnetic impurity spins on the surfaces of a thin Weyl-semimetal film couple to one another through the conduction electrons. It claims that when the two spins sit on opposite surfaces, the coupling is strong and non-monotonic in film thickness, reaching a maximum at a critical thickness set by the separation of the Weyl nodes in momentum space. It also claims that even same-surface couplings acquire a thickness-dependent contribution from the opposite surface, so a proper thin-film treatment must keep the Fermi-arc states of both surfaces. The authors establish this with two independent methods that agree: analytical slab wavefunctions for the low-energy states and a recursive Green's function computation on the tight-binding model.","feed_headline":"Opposite-surface spin coupling peaks at a critical film thickness","feed_subtitle":"Fermi arcs and bulk states make the coupling strong and non-monotonic, set by Weyl-node separation.","key_machinery":"The argument is carried by two complementary tools. Analytically, the authors solve the low-energy Dirac Hamiltonian of each time-reversed block in a slab with infinite-mass boundary conditions at the two surfaces; the resulting transcendental equation labels a family of $n=0$ bands that contain both the Fermi-arc surface states and the low-energy bulk states, and the penetration length of the arc states diverges as the surface momentum approaches the Weyl-node projection. Numerically, a recursive Green's function scheme built from the layer structure of the tight-binding Hamiltonian computes the surface-to-surface Green's function elements $G_{1,N_z}$ and $G_{N_z,1}$, as well as the intrasurface ones, essentially exactly, allowing the full RKKY tensor to be evaluated for a slab of arbitrary thickness. The matching of these two methods is what lets the paper attribute the non-monotonic thickness dependence to Fermi-arc penetration and bulk-state competition.","core_discovery":"On the paper's own terms, the central discovery is that the RKKY coupling between impurity spins on opposite surfaces of a Weyl-semimetal slab is neither weak nor monotonic: it grows as the film thickens, peaks at an optimum thickness $L_c$, and then decays rapidly as bulk-state dominance sets in. The maximum exists because the Fermi-arc states localize more firmly on their surfaces as the thickness increases, raising the surface density of states, while the number of conduction-electron states sensitive to both surfaces decreases; the competition defines a thin-film limit with $L_c \\sim 1/k_0$, where $k_0$ is the momentum-space separation of the Weyl nodes. For fixed thickness the opposite-surface coupling falls off roughly as $(R^2+L_z^2)^{-\\alpha/2}$ with $\\alpha$ between 3 and 4, and for large thickness as $L_z^{-\\gamma}$ with $\\gamma\\approx 5$, the same falloff as bulk Weyl RKKY. A related result is that the same-surface component $J_{xx}$, which vanishes for straight Fermi arcs in a single-surface model, remains nonzero in a finite slab and disappears with increasing thickness, showing that the second surface is responsible.","pith_inferences":["The divergent penetration of Fermi-arc states near Weyl-node projections suggests that the optimum thickness and peak coupling should depend on where the impurity spins sit laterally on the surface relative to the node projections; mapping $L_c$ as a function of position could image the arcs' real-space extent.","If a real material has curved Fermi arcs or a different spin texture, the non-monotonic peak may shift to a different thickness, split among spin components, or weaken; a material-specific tight-binding calculation would show which of these happens.","Because the coupling strength is controlled by thickness, a WSM slab could act as a tunable magnetic coupler: changing the film thickness by a few lattice spacings near $L_c$ would switch the effective intersurface spin interaction from small to maximal."],"forward_implications":["For films thinner than the critical thickness $L_c\\sim 1/k_0$, RKKY calculations that model only a single surface are incomplete: the opposite surface contributes even to same-surface spin couplings like $J_{xx}$.","Opposite-surface RKKY coupling can be strong and long-ranged enough to favor ferromagnetic ordering between the two surfaces, with parallel surface magnetizations in the ground state.","In materials such as TaAs, with Weyl-node separation $k_0\\approx 0.1\\pi/a$, the predicted critical thickness is tens to a hundred lattice spacings, placing the effect in experimentally accessible thin films.","At finite chemical potential the RKKY coupling acquires the usual $2k_F$ oscillations while the non-monotonic thickness envelope survives, so the optimum-thickness phenomenon is not an artifact of the nodal limit."],"supporting_citations":[{"why":"Defines the RKKY interaction formula, Eq. (18), that the entire calculation evaluates.","marker":"2"},{"why":"Provides the prior single-surface Fermi-arc RKKY results and asymptotic Green's function that the thick-slab limit must recover.","marker":"19"},{"why":"Supplies the time-reversal-symmetric four-node tight-binding Hamiltonian underlying the slab and numerical Green's function computations.","marker":"22"},{"why":"Supplies the recursive surface Green's function scheme used for the numerical results.","marker":"23"},{"why":"Gives the real-material Weyl-node separation estimate used to translate the critical thickness into experimentally accessible film sizes.","marker":"24"}],"fun_headline_variants":["Opposite-surface RKKY peaks at critical film thickness","Spin coupling across Weyl slab peaks at optimal thickness","Fermi arcs drive strong RKKY across Weyl thin films","Critical thickness maximizes spin coupling in Weyl thin films"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the minimal tight-binding model with two straight, spin-polarized Fermi arcs per surface, and the low-energy states derived from it, faithfully represent a real Weyl-semimetal thin film; if real Fermi-arc connectivity, spin texture, or penetration behavior differs substantially, the predicted strong opposite-surface coupling and its optimum thickness could change qualitatively.","fun_headline_variants_meta":{"raw":{"variants":["Opposite-surface RKKY peaks at critical film thickness","Spin coupling across Weyl slab peaks at optimal thickness","Fermi arcs drive strong RKKY across Weyl thin films","Critical thickness maximizes spin coupling in Weyl thin films"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001259,"raw_usage":{"total_tokens":5138,"prompt_tokens":909,"completion_tokens":4229,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":525,"completion_tokens_details":{"reasoning_tokens":4161}},"tokens_in":525,"tokens_out":4229,"duration_ms":29914,"temperature":1.0,"reasoning_tokens":4161,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:39:42.545387+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the thickness dependence of the magnetic correlation between impurities on opposite surfaces of a Weyl-semimetal slab, for example in TaAs thin films at low temperature; the claim fails if the coupling decays monotonically with thickness without any maximum, or if the opposite-surface coupling is orders of magnitude weaker than same-surface coupling at all thicknesses. A less demanding check is to compute the same RKKY tensor in a material-specific model with curved or differently connected Fermi arcs and see whether the non-monotonic peak survives.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the RKKY interaction formula, Eq. (18), that the entire calculation evaluates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the prior single-surface Fermi-arc RKKY results and asymptotic Green's function that the thick-slab limit must recover."},{"cited_title":"Quantum transport simulation of non-local response in Weyl semimetals","cited_arxiv_id":"1812.05504","evidence_quote":"Supplies the time-reversal-symmetric four-node tight-binding Hamiltonian underlying the slab and numerical Green's function computations."}],"review_version":1}