{"id":"52d49bd6-b65d-4ac5-9684-16a6f3fce862","arxiv_id":"2506.12380","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"An exact beyond-relaxation-time Boltzmann calculation for a Rarita-Schwinger-Weyl node predicts opposite curvature in the two conduction bands and a sign flip controlled by interband scattering.","lead":"This paper calculates how electrical conductivity along a magnetic field behaves for a fourfold degenerate band crossing called a Rarita-Schwinger-Weyl node, solving the transport equations more carefully than earlier work. The predicted band-resolved magnetoconductivity curves bend in opposite directions and change sign when interband scattering is added, offering a fingerprint of this node.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Collision-term normalization appears to double-count scattering strengths: Eq. (23) already carries the β factors whose prefactor reappears in Eq. (26), so A N = Υ is not uniquely defined.","rationale":"The paper's central claim is conditional on solving a finite linear system. I find the polynomial ansatz concern less threatening than the reader did: since the angular kernel in Eq. (23) is a polynomial of degree three in cosθ and cosθ′, the collision integral maps any distribution to a cubic polynomial in cosθ, so Eq. (27) is an exact solution space, not a truncation. The truly load-bearing problem is the normalization of that kernel: the printed equations do not specify whether the β factors appear once or twice in the collision integral. The appendix matrix is linear in β, but the displayed T in Eq. (23) plus the prefactor in Eq. (26) is consistent only with β². Until this is resolved, the numerical curves and threshold claims cannot be reproduced. This keeps the manuscript in CONDITIONAL status: the framework is plausible and the derivation is mostly transparent, but the central quantitative signature rests on an ambiguity in the printed equations. The absence of code/data and unshown parameter scans reinforce the need for a normalization check, though they are secondary.","tokens_in":16326,"tokens_out":17606,"duration_ms":223735,"concrete_test":"Independently re-derive Eq. (26) by substituting Eq. (18) and Eq. (19) into Eq. (24); check whether the product of the prefactor and the β factors printed in Eq. (23) yields terms linear or quadratic in β in the A matrix. Then recompute σ_{zz}^{s}(B) for the Fig. 3 parameter sets using the normalization that matches the appendix (linear in β) and the naive double-counted version (quadratic in β). If the downward-curving/negative σ_{zz}^{3/2} above threshold does not persist under the corrected linear-in-β normalization, the central distinguishing claim is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (19) defines |V_{s,\\tilde{s}}|^2 = (16×2π/ρ_imp) β_{...}. Eq. (23) then displays T_{s,\\tilde{s}}(θ,θ′) with factors β_{1/2,1/2}^{intra}, β_{3/2,3/2}^{intra}, and β_{inter} already included. Eq. (26) multiplies this same T by ρ_imp|V_{s,\\tilde{s}}|^2/(16×4π) = β_{...}/2. Unless the β symbols in Eq. (23) are meant to denote something different from the β parameters of Eq. (19), every scattering amplitude enters the Boltzmann kernel twice. The appendix matrix A is linear in β (entries such as 3 c_2^1 β_{1/2,1/2}^{intra}), which is consistent only if T in Eq. (26) is a pure spinor-overlap angular function, not the β-weighted object printed in Eq. (23). This ambiguity controls the numerical solution of A N = Υ, and hence controls the claimed sign/curvature change of σ_{zz}^{3/2} as a function of β_{inter}/β_{intra}. The phrase \"exact computation\" cannot be validated from the printed equations until this normalization is settled.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the longitudinal magnetoconductivity of an isolated isotropic Rarita-Schwinger-Weyl node, using the semiclassical Boltzmann equation with Berry curvature and orbital magnetic moment included, beyond the constant relaxation-time approximation. Impurity scattering is treated with intraband and interband channels characterized by parameters β_intra and β_inter. The authors reduce the linearized Boltzmann equation to an 8×8 linear system via a cubic-polynomial ansatz for the mean-free path, impose charge conservation, and compute band-resolved δσ_zz(B). Their central claim is that the s=3/2 and s=1/2 bands respond in opposite ways, with interband scattering flipping the s=3/2 curve into the negative domain, which would distinguish RSWNs from ordinary Weyl nodes.","tokens_in":16581,"tokens_out":10118,"duration_ms":128006,"significance":"If the calculation is correct, the predicted sign and curvature reversal of the band-resolved longitudinal magnetoconductivity is a concrete and falsifiable distinction between RSWNs and conventional Weyl nodes, going beyond the relaxation-time approximation. The paper's strengths are the explicit spinor overlap computation, the inclusion of Berry curvature and orbital magnetic moment, and the attempt to solve the full linearized Boltzmann equation. However, the printed equations contain a normalization inconsistency that directly controls the numerical solution of A N = Υ, and the appendix is not self-contained. These issues must be resolved before the central claim can be accepted.","major_comments":[{"comment":"The normalization of the collision kernel is internally inconsistent. Equation (19) defines |V_{s,tilde s}|^2 = (16×2π/ρ_imp) β_{...}, and Eq. (23) then prints T_{s,tilde s}(θ,θ') as a β-weighted combination of angular functions. In Eq. (26), the prefactor ρ_imp |V_{s,tilde s}|^2/(16×4π) supplies an additional factor β_{...}/2. If Eq. (23) is used literally, every scattering amplitude enters the Boltzmann equation twice, and the Appendix matrix A, which is linear in β (e.g., entries like 3 c_2^1 β_{1/2,1/2}^{intra}), is inconsistent with Eq. (26) combined with Eq. (23). Since the numerical solution of A N = Υ determines the sign and curvature claims in Sec. III C, this ambiguity is load-bearing: either Eq. (23) should be the pure spinor-overlap polynomial with the β factors removed, or Eq. (26) and the Appendix must carry β² terms.","section":"Eqs. (19), (23), (26), and Appendix"},{"comment":"The definitions of c^n_{αs} and hc^n_{αs} in Eq. (32) are printed identically, yet the Υ vector of the Appendix contains terms such as 3 hc^2_2 β_inter + 3 hc^0_2 β_inter, which must involve an extra factor of the non-polynomial function h_s(θ) in the integrand. Without a correct definition of hc^n_{αs}, the reader cannot verify the right-hand side of A N = Υ, and the claimed 'exact' solution is not reproducible. The authors should provide the explicit integrands for both c^n and hc^n and show at least one row of A and Υ in full.","section":"Appendix, Eq. (32)"},{"comment":"The paper asserts that the eight-dimensional cubic ansatz spans the exact solution space because the collision kernel is polynomial in cosθ, but it does not demonstrate that the right-hand side of Eq. (26) contains no powers of cosθ higher than 3 after the non-polynomial factors τ_s, D_s, and k_F(θ) are integrated. The phrase 'observing the powers of cosθ in (23)' is insufficient for the 'exact computation' claim made in the Abstract and Sec. IV. Either the closure proof should be supplied, or the claim should be softened to state that the solution is obtained within a truncated polynomial basis.","section":"Sec. III B, Eq. (27)"}],"minor_comments":[{"comment":"There is a typo: 'bbtained' should be 'obtained'.","section":"Sec. III C"},{"comment":"The notation β_{1,1}^{intra}, β_{3,3}^{intra} is inconsistent with β_{1/2,1/2}^{intra}, β_{3/2,3/2}^{intra} used elsewhere; please unify the notation.","section":"Eq. (23)"},{"comment":"The horizontal axis is labeled 'B (in eV 2)' in the figures and text; the units of magnetic field in natural units should be specified clearly, and the axis label should be consistent with the definition σ_zz(B)/σ_zz(0)−1.","section":"Sec. III C and figure captions"},{"comment":"The sentence stating that rank deficiency 'prevents the system from becoming overdetermined' is confusing: a rank-7 8×8 system is underdetermined unless an additional independent constraint is added. Please clarify that charge conservation supplies the missing equation and state how linear independence is verified.","section":"Sec. III B"},{"comment":"Several typographical errors appear, including 'degenrate', 'form' for 'from', 'focussing', and 'quasipaticle'; a careful proofread is recommended.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is well aligned with the journal's scope in mesoscopic and topological transport. The main concern is not the physics idea but the reproducibility of the numerical solution: the normalization ambiguity in Eqs. (19), (23), and (26) and the defective definitions in Eq. (32) prevent the reader from verifying the central result. If the author can confirm that the Appendix was generated with pure spinor overlaps (linear β) rather than with the β-weighted T of Eq. (23), and can fix the appendix definitions, the paper could become publishable after a moderate revision. The 'exact' claim in the Abstract should be toned down or justified if the cubic ansatz is not proven to close."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, it does deliver something genuinely new: the first beyond-RTA, band-resolved longitudinal magnetoconductivity for an isolated Rarita-Schwinger-Weyl node, with Berry curvature and orbital magnetic moment both included. Second, the central qualitative claim—that interband scattering above a small threshold flips the s=3/2 response negative while s=1/2 stays positive—is plausible but rests on printed equations that contain a real notation ambiguity about the beta factors, so I would not sign off on the word \"exact\" until that is cleaned up.\n\nWhat is actually new: the formalism is taken from the author's earlier Kramers-Weyl and pseudospin-1 work plus Knoll-Timm-Meng, but the RSWN case with its s=1/2 and s=3/2 Fermi surfaces is a legitimate extension, and the sign-change prediction is a concrete, falsifiable fingerprint that distinguishes multifold from ordinary Weyl nodes. The paper is readable and gives enough equations to reconstruct the calculation, which is more than most transport papers do.\n\nThe main soft spot is the beta normalization. Equation (19) defines |V|^2 proportional to beta, but Eq. (23) prints T with beta factors already inside, and Eq. (26) multiplies that T by rho_imp|V|^2/(16 pi), which is again beta/2. If both are taken literally, every scattering amplitude enters squared. The appendix matrix A is linear in beta, so the actual implementation must be using T as a pure overlap function, not the beta-weighted object in Eq. (23). That is a fixable notation issue, but it is exactly the kind of thing that makes the \"exact computation\" claim impossible to verify from the manuscript alone.\n\nThe other concern in the reader report—that the cubic-polynomial ansatz is unproven because D_s and k_F are non-polynomial—does not hold up. The collision integral maps a function of theta' through T(theta,theta'), which is polynomial of degree three in cos(theta). So the RHS is exactly a cubic polynomial in cos(theta) regardless of the non-polynomial measure factors. The ansatz closes exactly; the integrals in the appendix absorb all the non-polynomial physics. So that weakest assumption is actually not weak.\n\nWhat is missing: the parameter scan behind the threshold claim is not shown, only representative plots. That is a reproducibility gap, not a fatal one.\n\nRecommendation: send it to peer review. A solid referee can sort out the beta notation and check the numerics. The qualitative prediction is worth putting on record.","headline":"A genuinely new RSWN transport calculation with a plausible qualitative prediction, but a beta-factor notation ambiguity in the printed equations must be fixed before the 'exact' claim is verifiable.","tokens_in":17163,"tokens_out":5617,"would_cite":true,"duration_ms":65960,"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":"Interband scattering flips the s=3/2 band's longitudinal magnetoconductivity negative in an isolated Rarita-Schwinger-Weyl node, giving a transport fingerprint that distinguishes it from ordinary Weyl nodes.","keywords":["Rarita-Schwinger-Weyl node","longitudinal magnetoconductivity","semiclassical Boltzmann equation","Berry curvature","orbital magnetic moment","interband scattering","multifold fermions","chiral semimetal"],"falsifier":"Recompute $\\sigma_{zz}(B)$ from the same linearized Boltzmann equation without imposing the cubic ansatz -- by expanding $\\Lambda_z^s(\\mu,\\theta)$ in Legendre polynomials up to high order for the paper's parameter values -- and check whether coefficients beyond $\\cos^3\\theta$ vanish and whether the $s=3/2$ band still turns negative; any significant higher harmonic or a different sign would falsify the central claim.","tokens_in":16054,"feed_emoji":"🧲","tokens_out":16004,"duration_ms":265852,"temperature":0.7,"pith_summary":"The paper aims to establish that the longitudinal magnetoconductivity of an isolated Rarita-Schwinger-Weyl node -- the current response when electric and magnetic fields are applied parallel to each other -- carries a band-resolved signature that separates the fourfold node from ordinary Weyl nodes. It solves the linearized semiclassical Boltzmann equation without the relaxation-time approximation, keeping the full angle-dependent spinor overlaps, Berry curvature, and orbital magnetic moment. The central result is that once interband scattering between the $s=1/2$ and $s=3/2$ bands exceeds a small threshold, the $s=3/2$ band's magnetoconductivity turns negative while the $s=1/2$ band remains positive. If correct, this sign flip gives a transport fingerprint of the Rarita-Schwinger-Weyl node and shows that momentum-independent relaxation-time estimates can miss the qualitative physics.","feed_headline":"Sign flip in magnetoconductivity fingerprints a Rarita-Schwinger-Weyl node","feed_subtitle":"For an isolated fourfold node, interband scattering turns the s=3/2 band's response negative while s=1/2 stays positive.","key_machinery":"The central object is the cubic-polynomial ansatz for the mean-free path, Eq. (27): $\\Lambda_z^s(\\mu,\\theta) = \\tau_s(\\mu,\\theta)[\\lambda_s - h_s + a_s\\cos\\theta + b_s\\cos^2\\theta + c_s\\cos^3\\theta]$. The argument relies on the claim that the collision operator closes on this four-dimensional space for each band, because the spinor-overlap functions $T_{s,\\tilde{s}}(\\theta,\\theta')$ are polynomials of degree three in $\\cos\\theta$. This turns the linearized Boltzmann equation into an $8\\times 8$ linear system for the coefficients $\\{\\lambda_s, a_s, b_s, c_s\\}$; the matrix has rank 7, and the missing independent equation is provided by electron-number conservation. The $\\theta$-dependent Fermi-surface radii, the phase-space factor $D_s$, and the field-induced Fermi-surface displacements enter through the integrals $c_{\\alpha s}^n$ that build the matrix and source vector, which is how the magnetic field enters the solution.","core_discovery":"On the paper's own terms, the discovery is that the two occupied bands of an isolated Rarita-Schwinger-Weyl node respond oppositely to a collinear magnetic field once interband scattering is allowed. Starting from $H = v_F\\,\\mathbf{k}\\cdot\\mathbf{J}$ with bands $s = \\pm 1/2, \\pm 3/2$, the Berry curvature $\\Omega_s$ enters through the phase-space factor $D_s = [1 + e\\mathbf{B}\\cdot\\Omega_s]^{-1}$, and the orbital magnetic moment distorts the Fermi surfaces by $\\varepsilon_s^{(m)} = -\\mathbf{B}\\cdot\\mathbf{m}_s$. Solving the linearized Boltzmann equation with the cubic-in-$\\cos\\theta$ ansatz for the mean free path and imposing charge conservation yields the band-resolved longitudinal magnetoconductivity. The numerical solutions show that for $\\beta_{\\mathrm{inter}}/\\beta_{\\mathrm{intra}}$ above a small threshold the $s=3/2$ contribution to $\\delta\\sigma_{zz}$ flips negative while the $s=1/2$ contribution curves upward and stays positive; with the orbital magnetic moment switched off the two curves reverse their behavior, showing that the OMM is the ingredient that pushes $s=1/2$ positive and pulls $s=3/2$ down. The paper presents this exact computation as correcting the earlier relaxation-time-approximation results and as the distinguishing feature of a RSWN.","pith_inferences":["If the cubic closure is exact, the sign-flip criterion should be a generic property of isolated Rarita-Schwinger-Weyl nodes, and a measurement that isolates the two Fermi pockets at positive chemical potential would be a sharper test than measuring the total conductivity.","The 'exact' designation depends on the ansatz closing; a numerical solution of the full linearized Boltzmann equation expanded in higher angular harmonics at representative parameters would settle whether coefficients beyond $\\cos^3\\theta$ are genuinely absent or merely small.","The same machinery with angle-dependent spinor overlaps could be applied to other multifold fermions, such as pseudospin-1 and sixfold nodes, where interband scattering may also produce sign changes that a relaxation-time treatment would miss.","In real materials such as RhSi the isolated-node assumption is tied to the chemical potential sitting near the node; doping away from that region would bring other Fermi pockets and could wash out the predicted band-resolved sign flip."],"forward_implications":["The exact solution shows that a momentum-independent relaxation time is not adequate for Rarita-Schwinger-Weyl nodes: keeping the angular dependence of the spinor overlaps and the mean free path changes the band-resolved magnetoconductivity, including its sign.","For an isolated Rarita-Schwinger-Weyl node with interband scattering above a small threshold, the $s=3/2$ band contribution to $\\delta\\sigma_{zz}$ turns negative while the $s=1/2$ band contribution remains positive and curved upward, providing a transport signature absent in two-fold Weyl nodes.","The orbital magnetic moment contributes with opposite signs in the two bands: positive for $s=1/2$, enough to flip the Berry-curvature-only response positive, and negative for $s=3/2$, pulling the response downward.","$\\sigma_{zz}(B)$ contains only even powers of $B$, consistent with Onsager-Casimir reciprocity; odd-in-$B$ terms are absent because the RSWN Hamiltonian has no tilt term.","The decoupled limit $\\beta_{\\mathrm{inter}}=0$ obeys charge conservation band by band, whereas $\\beta_{\\mathrm{inter}}\\neq 0$ conserves only the total charge; this difference is why turning on interband scattering qualitatively reorganizes the response rather than merely renormalizing it."],"supporting_citations":[{"why":"Supplies the linearized-Boltzmann-equation method that the paper adapts and generalizes to the RSWN.","marker":"[12]"},{"why":"Gives the earlier relaxation-time-approximation calculation for RSWNs that the exact computation is meant to correct.","marker":"[33]"},{"why":"Gives the earlier internode-scattering relaxation-time treatment whose insufficiency motivates the present work.","marker":"[34]"},{"why":"Applies the same exact formalism to Kramers-Weyl nodes, providing the comparison case of how interband scattering reshapes band-resolved responses.","marker":"[35]"},{"why":"Provides the material parameters (including v_F = 1.23 eV Å) used in the plotted magnetoconductivity curves.","marker":"[58]"},{"why":"Identifies the isolated RSWN in RhSi and the ~0.4 eV energy separation that justifies modeling a single node.","marker":"[43]"}],"fun_headline_variants":["Sign flip in magnetoconductivity singles out Rarita-Schwinger-Weyl node","Opposite band responses expose fourfold RSWN in magnetoconductivity","Exact Boltzmann solution reveals band-resolved sign change at RSWN","Interband scattering flips s=3/2 response, s=1/2 stays positive at RSWN"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the exact mean-free path has no angular dependence beyond $\\cos^3\\theta$, so the cubic ansatz closes under the collision operator; if higher angular harmonics are generated by the scattering, the predicted sign and curvature of the magnetoconductivity could change.","fun_headline_variants_meta":{"raw":{"variants":["Sign flip in magnetoconductivity singles out Rarita-Schwinger-Weyl node","Opposite band responses expose fourfold RSWN in magnetoconductivity","Exact Boltzmann solution reveals band-resolved sign change at RSWN","Interband scattering flips s=3/2 response, s=1/2 stays positive at RSWN"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001021,"raw_usage":{"total_tokens":4351,"prompt_tokens":1034,"completion_tokens":3317,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":650,"completion_tokens_details":{"reasoning_tokens":3224}},"tokens_in":650,"tokens_out":3317,"duration_ms":29790,"temperature":1.0,"reasoning_tokens":3224,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:52:24.731025+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute $\\sigma_{zz}(B)$ from the same linearized Boltzmann equation without imposing the cubic ansatz -- by expanding $\\Lambda_z^s(\\mu,\\theta)$ in Legendre polynomials up to high order for the paper's parameter values -- and check whether coefficients beyond $\\cos^3\\theta$ vanish and whether the $s=3/2$ band still turns negative; any significant higher harmonic or a different sign would falsify the central claim.","supporting_citations":[{"cited_title":"Flicker, F","cited_arxiv_id":null,"evidence_quote":"Gives the earlier relaxation-time-approximation calculation for RSWNs that the exact computation is meant to correct."},{"cited_title":"Mandal, Chiral anomaly and internode scatterings in multifold semimetals, Phys","cited_arxiv_id":null,"evidence_quote":"Applies the same exact formalism to Kramers-Weyl nodes, providing the comparison case of how interband scattering reshapes band-resolved responses."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies the isolated RSWN in RhSi and the ~0.4 eV energy separation that justifies modeling a single node."}],"review_version":1}