{"id":"fea2d282-8a17-4277-b5b1-5eac789c686e","arxiv_id":"1908.08575","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Current-induced spin polarization combined with spin-orbit defect scattering generates a bilinear magnetoresistance that scales as j*b*sinθ/|εF|^3 and dominates over the warping mechanism at low Fermi energies.","lead":"This paper derives a new mechanism for 'bilinear magnetoresistance' in topological insulator surfaces, where resistance depends on both current and magnetic field direction. It predicts this effect grows strongly at low Fermi energies and could explain experiments in materials without the usual hexagonal warping.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (14) has an internal ℏ-factor error: the printed BMR amplitude is inconsistent with Eq. (13); the quoted numerical values match the corrected (ℏ^2 numerator) form, so the central mechanism is not invalidated but the headline formula needs correction.","rationale":"The reader's verdict of CONDITIONAL is appropriate. The most concrete load-bearing problem is the internal inconsistency in the printed BMR formula: Eq. (14) and Eq. (17) place ℏ^{-2} where the preceding derivation of Eq. (13) requires ℏ^2. This is not a conceptual flaw in the mechanism, since the quoted numbers and the comparison in Fig. 2e use the correct prefactor, but it is a serious typo in the central quantitative result that must be corrected before the paper is usable. The reader's weakest_assumption, the field-independence of the current-induced spin polarization, is less damaging than stated: any correction to S_y linear in B is also linear in j and would contribute to the BMR only at order B^2 j, not at the leading order B j. Still, the paper's S1 justification for this assumption is thin, relying on a citation to a related calculation rather than a derivation in the same model, and a public derivation would strengthen the paper. The sign of the BMR term also appears convention-dependent but was not the main issue. Overall, the central mechanism is plausible and supported by a closed Green's-function calculation; the correct remedy is a minor revision fixing Eq. (14) and clarifying the CISP assumption, not rejection.","tokens_in":11348,"tokens_out":29649,"duration_ms":325312,"concrete_test":"Recompute the BMR amplitude from Eq. (13) by substituting S_y = ℏ^2 j_x/(2 e v), J = -8π v_F/k_F, and B = g μ_B b, then evaluate the numerical estimates for Bi2Se3 at ε_F = 0.256 eV and 0.02 eV with v_F = 5×10^5 m/s, g = 2, b = 5 T, j = 10 A/m. If the amplitude is taken with ℏ^2 in the numerator, the results match the quoted 4.40×10^{-5} and 9.2×10^{-2}; with the printed ℏ^{-2} denominator, the results differ by a factor of roughly ℏ^4 and are unphysical. This single check settles whether Eq. (14) is merely a typo.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is the BMR amplitude in Eq. (14)/(17). Tracing the coefficient from Eq. (13) with S_y = ℏ^2 j_x/(2 e v), J = -8π v_F/k_F, B = g μ_B b, and k_F = |ε_F|/(ℏ v_F) gives a prefactor proportional to ℏ^2 v_F/|ε_F|^3, not ℏ^{-2} v_F/|ε_F|^3. Concretely, the linear term in Eq. (13) is (15 h/|e|) j_x B sinθ/(k_F ε_F^2); substituting h=2πℏ and k_F=|ε_F|/(ℏ v_F) yields 30π ℏ^2 g μ_B v_F j b sinθ/(|e| |ε_F|^3). The printed Eq. (14) has the same factor divided by ℏ^2 instead of multiplied, i.e., off by ℏ^4. The numerical values quoted in the text (4.40×10^{-5} for ε_F=0.256 eV and 9.2×10^{-2} for ε_F=0.02 eV) are reproduced by the ℏ^2-numerator expression, so the error appears to be typographical rather than a failure of the derivation. This matters because the printed formula is the paper's headline result and would mislead any reader using it directly. Separately, the S1 justification for taking the current-induced spin polarization at its zero-field value is only a citation to a related calculation; however, a B-linear correction to S_y would enter the BMR only at order B^2 j, so this assumption does not threaten the leading bilinear coefficient.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a new mechanism for bilinear magnetoresistance (BMR) in topological insulator surface states. Using the minimal Dirac-cone model without hexagonal warping, the authors include the current-induced spin polarization (CISP) as an effective exchange field J S⋅σ, perform a gauge transformation that removes the combined in-plane field B + J S from the clean Hamiltonian, and study scattering by random spin-orbit defects in the Born approximation with ladder vertex corrections. The central result is the longitudinal conductivity (Eq. 13), from which they extract a BMR amplitude (Eq. 14) that grows linearly in current and magnetic field, varies as sinθ, scales as v_F/|ε_F|^3, and is independent of disorder strength. They compare this with the hexagonal-warping mechanism and argue that the CISP mechanism dominates at low Fermi energies, with numerical estimates for Bi2Se3. The derivation is presented in the main text and in supplementary sections S1–S3.","tokens_in":11719,"tokens_out":24089,"duration_ms":240137,"significance":"The proposed mechanism is physically interesting and timely: it offers an explanation for BMR in systems without hexagonal warping (e.g., α-Sn(001)) and makes falsifiable predictions, namely linear scaling in current and magnetic field, sinθ angular dependence, 1/|ε_F|^3 Fermi-energy scaling, and a disorder-strength-independent amplitude. The calculation is a standard Green's-function treatment with Born self-energies and ladder vertex corrections, and no BMR data are used to set constants; the coupling J and the spin polarization S_y are derived within the model, making the amplitude parameter-free apart from band-structure parameters. However, the headline formula Eq. (14) contains an internal ℏ-factor error that must be fixed, and the sign of the effect needs to be clarified. With those corrections, the central claim appears sound and the paper would be a valuable contribution to the field.","major_comments":[{"comment":"Eq. (14) (and the identical prefactor in Eq. (17)) is inconsistent with Eq. (13). Taking the coefficient 15h/(|e| k_F ε_F^2) j B sinθ from Eq. (13) and substituting k_F = |ε_F|/(ℏ v_F) and B = g μ_B b gives 30π g μ_B ℏ^2 v_F/(|e| |ε_F|^3) j b sinθ, i.e., a prefactor proportional to ℏ^2, not ℏ^{-2}. The numerical values quoted in the Discussion (4.40×10^{-5} for ε_F = 0.256 eV and 9.2×10^{-2} for ε_F = 0.02 eV) reproduce the ℏ^2-numerator expression, so the error is likely typographical rather than a failure of the derivation. Nevertheless, since this is the paper's headline quantitative prediction, the formula must be corrected before publication.","section":"Discussion, Eq. (14) and Eq. (17)"},{"comment":"The sign of the BMR in Eq. (14) appears inconsistent with Eq. (13) under the paper's own conventions. With J = -8π v_F/k_F from S1 and S_y = ℏ^2 j_x/(2 e v) from S1, the linear term in Eq. (13) is positive for j_x > 0 and B_y > 0, so ρ(B)/ρ(0) - 1 ≈ -c j_x B_y with c > 0. Using the paper's definition BMR = [MR(j_x = j) - MR(j_x = -j)]/2 then gives BMR = -c j b sinθ, opposite to the sign in Eq. (14). The authors should check this sign, state their chirality and charge-sign conventions explicitly, and correct Eq. (14) accordingly.","section":"Discussion, Eq. (14) versus Eq. (13)"}],"minor_comments":[{"comment":"The assumption that S_y retains its zero-field value in the presence of B and J S is justified only by a citation to prior work; since any B-linear correction to S_y would enter the BMR at order B^2 j, not at the leading bilinear order, this assumption is safe for the leading coefficient, but it would be helpful to state this argument explicitly in S1.","section":"S1, Eq. (22)"},{"comment":"In S2 the relaxation rate Γ0 is written with λ^2, whereas the scattering strength is denoted α elsewhere; this is presumably a typographical inconsistency and should be corrected.","section":"S2, definition of Γ0"},{"comment":"There are several typographical and formatting errors: 'm odel' in the title, 'antysymmetric' for 'antisymmetric', inconsistent use of 'bilinear' and 'bi-linear', and garbled accents in the author affiliation.","section":"Title and text"},{"comment":"Reference [30] is cited as a private communication for the experimental observation in α-Sn(001); a published reference would be more appropriate, and reference [48] is cited as 'to be published' and should be either published or described in enough detail for the reader to verify the claim.","section":"References"},{"comment":"The numerical comparison at ε_F = 0.02 eV is made deep in the low-energy regime where the minimal Dirac model may not be quantitatively reliable; the authors should state explicitly how the validity condition |ε_F| ≫ Γ0 is met for the disorder parameters used in the estimates.","section":"Discussion, numerical estimates"}],"recommendation":"major_revision","confidential_remarks":"The ℏ-factor error in Eq. (14) is almost certainly typographical, because the numerical estimates in the Discussion use the corrected ℏ^2 numerator. The sign discrepancy between Eq. (13) and Eq. (14) is potentially more substantive and should be resolved carefully; if the authors intend e in Eq. (13) to be the signed electron charge, they must state this explicitly, since Eq. (14) uses |e|. The reliance on an unpublished reference and a private communication is a fit-to-journal concern but not a reason to reject. In my assessment, the mechanism and the central derivation are defensible, and the paper is publishable after the requested corrections."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe paper identifies a new mechanism for bilinear magnetoresistance in TI surface states, and the central derivation holds up. The mechanism relies on current-induced spin polarization plus scattering off spin-orbit defects, and it works without hexagonal warping. That's new, and it explains why BMR appears in α-Sn where warping is absent.\n\nWhat the paper does well: the Green's function calculation in Born approximation with ladder vertex corrections is standard but competently executed. The final conductivity formula, Eq. (13), has the expected structure, and the comparison with the Zhang–Vignale warping mechanism gives a clean, testable prediction: the CISP mechanism dominates at low Fermi energies, with a 1/|ε_F|^3 enhancement. The derivation is honest, not fitted.\n\nSoft spots, in order of severity. First, Eq. (14) contains a concrete factor error: the printed second expression has ℏ^{-2} in the denominator, while tracing the derivation from Eq. (13) gives ℏ^2 in the numerator. The numerical values quoted later match the corrected form, so this is almost certainly a typographical error, but it will mislead anyone who copies the formula directly. The authors should fix it.\n\nSecond, the assumption that the current-induced spin polarization S_y remains at its zero-field value is justified only by a citation to the authors' own pending work. Annoying, but not fatal: a B-linear correction to S_y would enter the BMR only at order B^2 j, so the leading bilinear coefficient is unaffected.\n\nThird, the model includes only spin-orbit defects, not scalar disorder. Real TI surfaces have both, and if scalar disorder dominates, the predicted amplitude could be smaller. The paper should at least comment on this.\n\nThe reliance on a private communication for the α-Sn observation is weak, but it's a theory letter; the mechanism stands on its own.\n\nWho this is for: researchers working on nonlinear magnetotransport in topological insulators, Rashba systems, or spintronics. The paper deserves a serious referee; it should not be desk-rejected. A referee should require the Eq. (14) correction, a public derivation or reference for the S_y(B=0) assumption, and a short discussion of scalar disorder. With those changes, this is a solid contribution.\n\nMy recommendation: engage with it. The central mechanism is plausible, and the main quantitative claim holds up once the typo is fixed.\n\nBest,\n[Your name]","headline":"A mostly sound new mechanism for bilinear magnetoresistance in TI surface states, with a factor-ℏ^4 typo in the headline formula that must be corrected.","tokens_in":12287,"tokens_out":5698,"would_cite":true,"duration_ms":51561,"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":"In the minimal model of topological-insulator surface states, bi-linear magnetoresistance arises from current-induced spin polarization acting through magnetic-field-dependent spin-orbit-defect scattering, and this mechanism dominates at…","keywords":["bilinear magnetoresistance","topological insulator surface states","current-induced spin polarization","spin-orbit defect scattering","unidirectional magnetoresistance","Dirac cone","hexagonal warping","relaxation rate anisotropy"],"falsifier":"Gate-tune the Fermi energy in a topological insulator with a circular, unwarped surface Fermi surface and measure the antisymmetric magnetoresistance as a function of $b$, $j$, and Fermi energy. If the mechanism holds, the bilinear amplitude should follow $\\left(\\frac{30\\pi g\\mu_B}{\\hbar^2 |e|}\\right)\\frac{v_F}{|\\varepsilon_F|^3} j b \\sin\\theta$, growing steeply as $|\\varepsilon_F|$ is reduced within the regime $|\\varepsilon_F| \\gg B, \\Gamma_0$; if instead the signal scales with the hexagonal-warping parameter and grows with $|\\varepsilon_F|$, the proposed mechanism is not the operative one.","tokens_in":11102,"feed_emoji":"🧲","tokens_out":18190,"duration_ms":155424,"temperature":0.7,"pith_summary":"The paper proposes and calculates a new microscopic mechanism for bi-linear magnetoresistance (BMR) in topological-insulator surface states: a resistance change linear in both the applied current and the magnetic field, so that reversing the field is equivalent to reversing the current. The mechanism needs no hexagonal Fermi-surface warping. Instead, a current-induced spin polarization $S_y = \\hbar^2 j_x/(2ev)$ combines with the magnetic field inside the scattering rate for random spin-orbit defects, making relaxation depend on the sign of $\\mathbf{B}\\cdot\\mathbf{j}$. Within the minimal Dirac-cone model the longitudinal conductivity becomes $\\sigma_{xx} = \\sigma_{0xx}\\left[1 - \\frac{15}{2}\\frac{B_y J S_y}{\\varepsilon_F^2} - 3\\frac{B^2}{\\varepsilon_F^2}\\right]$, giving a BMR amplitude $\\left(\\frac{30\\pi g\\mu_B}{\\hbar^2 |e|}\\right)\\frac{v_F}{|\\varepsilon_F|^3} j b \\sin\\theta$. Because this amplitude grows as $|\\varepsilon_F|^{-3}$, the mechanism dominates the warping-based one at low Fermi energies, consistent with observations of BMR in materials without hexagonal symmetry.","feed_headline":"Current-driven spin polarization creates one-way magnetoresistance","feed_subtitle":"The resistance term grows steeply near the Dirac point and needs no Fermi-surface warping.","key_machinery":"The central object is the effective in-plane field $\\mathbf{B}_{\\mathrm{eff}} = \\mathbf{B} + J\\mathbf{S}$, where $\\mathbf{S}$ is the current-induced spin polarization $S_y = \\hbar^2 j_x/(2ev)$. A gauge transformation $\\mathbf{k}\\to\\mathbf{q}-\\frac{e}{\\hbar}\\boldsymbol\\Lambda$ removes $\\mathbf{B}_{\\mathrm{eff}}$ from the Dirac Hamiltonian, but the spin-orbit defect scattering potential retains a residual term proportional to $B_x$ and to $B_y+J S_y$. The Born self-energy of that residual term makes the quasiparticle relaxation rate $\\Gamma(\\varphi)$ depend on the direction of the momentum on the Fermi circle, and the ladder vertex correction converts that angular dependence into a term in $\\sigma_{xx}$ linear in $\\mathbf{B}$ and $\\mathbf{j}$. The mechanism is therefore field-dependent relaxation in the presence of a non-equilibrium spin polarization, not Fermi-surface warping.","core_discovery":"The paper establishes that in the minimal Dirac-cone model of a topological-insulator surface, described by $\\hat H_0 = v(\\mathbf{k}\\times\\hat z)\\cdot\\boldsymbol\\sigma + \\mathbf{B}\\cdot\\boldsymbol\\sigma + J\\mathbf{S}\\cdot\\boldsymbol\\sigma$, scattering from local spin-orbit fluctuations becomes magnetic-field dependent after the gauge transformation that removes the effective field $\\mathbf{B}+J\\mathbf{S}$ from the band Hamiltonian. The Born self-energy then yields an angle-dependent relaxation rate $\\Gamma(\\varphi) = \\Gamma_0\\left[1 + \\frac{2B_y J S_y + B^2}{\\varepsilon^2} + \\frac{2}{|\\varepsilon|}\\left((B_y+J S_y)\\cos\\varphi - B_x\\sin\\varphi\\right)\\right]$ at the Fermi surface. This angular dependence enters the ladder vertex correction and produces the longitudinal conductivity $\\sigma_{xx} = \\sigma_{0xx}\\left[1 - \\frac{15}{2}\\frac{B_y J S_y}{\\varepsilon_F^2} - 3\\frac{B^2}{\\varepsilon_F^2}\\right]$. The corresponding magnetoresistance contains a bilinear term $\\mathrm{BMR} = \\left(\\frac{30\\pi g\\mu_B}{\\hbar^2 |e|}\\right)\\frac{v_F}{|\\varepsilon_F|^3} j b \\sin\\theta$ (with the sign set by the current direction), linear in current and field and antisymmetric under reversal of either, plus a symmetric quadratic term $\\mathrm{sMR} = 3g^2\\mu_B^2 b^2/\\varepsilon_F^2$. Since the bilinear amplitude decays as $|\\varepsilon_F|^{-3}$ while the warping-based amplitude grows linearly with $|\\varepsilon_F|$, the new mechanism dominates at low Fermi energies.","pith_inferences":["The paper keeps $S_y$ at its zero-field value even when $\\mathbf{B}$ and the exchange field are present; a fully self-consistent treatment that lets $\\mathbf{B}$ renormalize $S_y$ is the natural next step, and could correct the prefactor without changing the symmetry of the effect.","Because the mechanism needs only current-induced spin polarization and spin-orbit disorder, it likely transfers to other spin-momentum-locked conductors, such as spin-split interface states, where the same gauge-shift logic would apply with a suitably generalized scattering potential.","The steep $1/|\\varepsilon_F|^3$ growth suggests BMR could serve as a quantitative local probe of current-induced spin polarization in gated devices, an application the paper points toward but does not develop.","A direct numerical calculation of the vertex correction beyond the ladder approximation would test whether the angular relaxation rate is the only carrier of the effect, or whether higher-order impurity scattering contributes comparably."],"forward_implications":["BMR should appear on any topological-insulator surface with spin-orbit disorder even when the Fermi surface is perfectly circular.","The bilinear signal should scale linearly with current density $j$ and magnetic field $b$ and vary as $\\sin\\theta$, vanishing when the field is parallel to the current.","Gating the Fermi energy toward the Dirac point should strongly enhance the BMR amplitude as $|\\varepsilon_F|^{-3}$, while the warping contribution would move in the opposite direction.","Because the symmetric quadratic magnetoresistance is independent of field orientation, current-reversal measurements separate the two terms and can in principle extract band parameters such as $k_F$, $v_F$, and the $g$-factor.","The mechanism explains why BMR persists in topological-insulator systems that lack hexagonal symmetry and thus cannot be described by warping alone."],"supporting_citations":[{"why":"Defines the current-induced spin polarization that the mechanism couples to.","marker":"[12]"},{"why":"Reports bilinear magnetoresistance in a three-dimensional topological insulator, the experimental fact the mechanism is built to explain.","marker":"[26]"},{"why":"Supplies the warping-based BMR theory whose amplitude is the comparison baseline.","marker":"[28]"},{"why":"Gives the gauge transformation used to move the effective in-plane field out of the band Hamiltonian and into the scattering potential.","marker":"[38]"},{"why":"Provides the random spin-orbit disorder potential used to model the defects.","marker":"[41]"},{"why":"Provides the Green-function and Born self-energy formalism on which the relaxation-rate and conductivity calculations rest.","marker":"[44]"}],"fun_headline_variants":["Bilinear magnetoresistance from current-induced spin polarization and spin-orbit defects","Spin-orbit defects give one-way magnetoresistance near Dirac point","Low-energy BMR from spin-orbit scattering, no warping needed","Magnetoresistance linear in current and field thanks to spin-orbit defects","One-way magnetoresistance emerges from spin-polarized scattering"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the current-induced spin polarization keeps its zero-field value even after the magnetic field and the polarization's own feedback field are switched on; if those fields substantially change the polarization, the predicted amplitude changes, and if the polarization is not proportional to current, the effect vanishes.","fun_headline_variants_meta":{"raw":{"variants":["Bilinear magnetoresistance from current-induced spin polarization and spin-orbit defects","Spin-orbit defects give one-way magnetoresistance near Dirac point","Low-energy BMR from spin-orbit scattering, no warping needed","Magnetoresistance linear in current and field thanks to spin-orbit defects","One-way magnetoresistance emerges from spin-polarized scattering"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00051,"raw_usage":{"total_tokens":2528,"prompt_tokens":1034,"completion_tokens":1494,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":650,"completion_tokens_details":{"reasoning_tokens":1400}},"tokens_in":650,"tokens_out":1494,"duration_ms":12107,"temperature":1.0,"reasoning_tokens":1400,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:38:32.082892+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Gate-tune the Fermi energy in a topological insulator with a circular, unwarped surface Fermi surface and measure the antisymmetric magnetoresistance as a function of $b$, $j$, and Fermi energy. If the mechanism holds, the bilinear amplitude should follow $\\left(\\frac{30\\pi g\\mu_B}{\\hbar^2 |e|}\\right)\\frac{v_F}{|\\varepsilon_F|^3} j b \\sin\\theta$, growing steeply as $|\\varepsilon_F|$ is reduced within the regime $|\\varepsilon_F| \\gg B, \\Gamma_0$; if instead the signal scales with the hexagonal-warping parameter and grows with $|\\varepsilon_F|$, the proposed mechanism is not the operative one.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the current-induced spin polarization that the mechanism couples to."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports bilinear magnetoresistance in a three-dimensional topological insulator, the experimental fact the mechanism is built to explain."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the warping-based BMR theory whose amplitude is the comparison baseline."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the gauge transformation used to move the effective in-plane field out of the band Hamiltonian and into the scattering potential."},{"cited_title":"Winkler, Spin-Orbit Interaction Effects in Two- Di- mensional Electron and Hole Systems (Springer- Verlag, Berlin, 2003)","cited_arxiv_id":null,"evidence_quote":"Provides the random spin-orbit disorder potential used to model the defects."},{"cited_title":"Strom, H","cited_arxiv_id":null,"evidence_quote":"Provides the Green-function and Born self-energy formalism on which the relaxation-rate and conductivity calculations rest."}],"review_version":1}