{"id":"c820b6fd-4fdb-461f-a8ca-7c44a9fbdbc3","arxiv_id":"2608.01782","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":1,"one_line_summary":"In the warped surface states of Bi2Te3 the linear spin Hall response is purely out-of-plane, and a second-order intrinsic in-plane spin current polarized along the field survives, but the claimed extrinsic collinear spin current is exactly zero by time-reversal symmetry.","lead":"This paper calculates how electron spin can flow as a pure spin current in the surface states of the topological insulator Bi2Te3 when an electric field is applied, extending the calculation beyond the usual linear response. The authors claim two second-order spin-current channels, but one of them, a current with spin pointing along the flow direction, is exactly zero by symmetry, so the headline application built on it disappears.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The extrinsic second-order spin current in Eq. (35) vanishes identically under time-reversal; the collinear extrinsic response (Fig. 4(b), Sec. IV) is therefore unsupported.","rationale":"I independently checked the symmetry of the central extrinsic integral. The reader's weakest assumption is exactly the load-bearing point: in a time-reversal-invariant band structure, the spin-weighted Berry curvature is T-even while the Fermi-surface derivative is T-odd, so the dipole integral cannot be finite. This is confirmed by the paper's own equations: d(-k)=-d(k), ε(-k)=ε(k), and Eq. (28) show the integrand is odd. I also note the secondary fragility (relaxation-time approximation) is real but independent; the T-symmetry cancellation holds even if the Boltzmann equation is accepted. The intrinsic second-order response, involving BCP and BC terms integrated against f, is structurally distinct and appears consistent. Thus a rejection of the paper's central claim is warranted, but the concern is specific to the extrinsic contribution. I agree with the reader's verdict and do not see a different more-load-bearing objection.","tokens_in":11980,"tokens_out":6210,"duration_ms":70695,"concrete_test":"Evaluate Eq. (35) numerically at μ = 0.5ε0 using the explicit Ξ^x, Ξ^y from Eqs. (32)–(33) and the dispersion (27). Discretize kx,ky on a uniform grid (k0 = 0.1 Å⁻¹, step ~0.002 k0, cutoff ~5 k0) and compute ∫ d²k Ξ^l ∂f/∂k_x with a low-temperature Fermi function (T=0.01 ε0). If the result is zero to numerical precision for both l=x and l=y, while a replica of the paper's method produces the nonzero values in Fig. 4(b), the extrinsic response is an artifact of the angular reduction in Appendix B. Repeat for μ = 0.2 and 1.0 ε0 to confirm.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The extrinsic SOSC in Eq. (35) is the load-bearing premise for the CPSC claim. In this model, Ξ^l_ab = Ω^ab S^l is even under k→−k: Eq. (26) gives d(−k)=−d(k), so S^l(−k)=−S^l(k) (spin expectation is ±d/|d|) and Ω^ab(−k)=−Ω^ab(k) from Eq. (28); hence Ξ^l(−k)=Ξ^l(k). The equilibrium occupation f(ε(k)) is also even because ε(−k)=ε(k) in Eq. (27), so ∂f/∂k_c is odd. The integrand of Eq. (35) is an odd function of k, and its integral over the full 2D Brillouin zone is identically zero for every l and c, for both bands (the two bands give identical Ξ^l). The same conclusion follows from the mirror symmetries displayed in Fig. 1(a): the in-plane Ξ components (Eqs. 32,33) are odd under one of the mirror reflections while ∂f/∂k_x is even under that reflection, so the contribution vanishes already before the T sum. Therefore Γ^{l,ext}=0, and the nonzero curve in Fig. 4(b), the CPSC discussion in Sec. IV, and the abstract's extrinsic in-plane spin-current claim cannot be reproduced from the paper's own equations. This is an internal symmetry inconsistency, not a disagreement with an external consensus. The intrinsic BCP/BC second-order response is not affected by this objection.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies linear and second-order pure spin currents in the hexagonal-warping surface states of Bi2Te3 within the conserved-spin-current formalism. Starting from the Fu model, the authors derive the k-resolved Berry curvature, BCP, and spin-weighted Berry curvature, and combine first-order perturbation theory with a Boltzmann relaxation-time description. They report (i) a purely out-of-plane linear spin Hall conductivity, (ii) a second-order intrinsic in-plane spin current polarized along the applied electric field, arising from BCP and Berry-curvature terms, and (iii) a second-order extrinsic in-plane collinear spin current attributed to a spin-weighted Berry curvature dipole. The paper also discusses experimental probes and potential applications of the collinear spin current.","tokens_in":12198,"tokens_out":7747,"duration_ms":92174,"significance":"If correct, the paper would identify a new extrinsic nonlinear spin response in topological insulator surfaces. The linear out-of-plane result and the intrinsic second-order field-polarized response are indeed consistent with the symmetries of the model, and the derivations in Eqs. (32)-(34) and Eqs. (23)-(24) have the expected parity properties. The work usefully shows that the conserved-spin-current formalism can be applied to the warped TI surface without introducing a mass term. However, the extrinsic second-order response is the paper's main new predictive claim, and it vanishes identically under time reversal. Since the abstract, Sec. IV, and Fig. 4(b) all rely on this extrinsic result, the central claim of the paper is unsupported by the paper's own equations.","major_comments":[{"comment":"The extrinsic SOSC vanishes identically. From Eq. (26), d(-k) = -d(k), so the dispersion is even, ε(-k)=ε(k); the spin expectation is odd, S^l(-k)=-S^l(k); and Eq. (28) gives Ω^xy(-k)=-Ω^xy(k). Therefore Ξ^l_yx(-k)=Ω^xy(-k)S^l(-k)=Ξ^l_yx(k), i.e. the spin-weighted Berry curvature is even under k→-k. Since f^(0)(k) is even, ∂f/∂k_x is odd, so the integrand in Eq. (35) is an odd function of k. The integral over the inversion-symmetric Brillouin zone is therefore zero for every spin component l and for both bands. The same conclusion follows from the Fermi-surface form in Eq. (B2): k_F and the denominator are even under φ→φ+π, while ∂ε/∂k_x is odd, so the angular integral vanishes. Thus Γ^{l,ext}=0; the nonzero Γ^{y,ext} in Fig. 4(b) cannot be reproduced from the paper's own formulas.","section":"III.B, Eq. (35) and Appendix B"},{"comment":"The abstract's claim that the second-order response exhibits an extrinsic in-plane spin current 'induced by ... Berry curvature dipole (BCD)', and Sec. IV's conclusion that 'The extrinsic SOSC ... leads to a collinearly polarized spin current', rest entirely on Γ^{l,ext} in Eq. (35). Since that quantity is zero by time-reversal symmetry, these statements are unsupported. This is not a matter of numerical magnitude or of the relaxation-time approximation; it is a symmetry property of the model. A revision that wishes to retain an extrinsic response would need to break time-reversal symmetry (e.g. by magnetic doping) or change the definition of the distribution function, neither of which is present in the manuscript.","section":"Abstract and Sec. IV"}],"minor_comments":[{"comment":"The terminology 'spin-dependent BCD' is misleading: the quantity in Eq. (35) is the product of Berry curvature and spin expectation, and its dipole is not a standard Berry curvature dipole. The vanishing of this dipole is a direct consequence of time-reversal symmetry, and the text should reflect that.","section":"Sec. IV"},{"comment":"The figure shows a nonzero extrinsic y-polarized conductivity, but the analytic integral in Eq. (B2) is odd under φ→φ+π and hence integrates to zero. The figure should be re-examined; if it was obtained numerically, the integration grid or the inclusion of the Jacobian may have broken the symmetry.","section":"Fig. 4(b)"},{"comment":"The relaxation-time approximation in Eq. (13) is applied to the surface states without modeling surface-bulk coupling. This is a secondary issue for the intrinsic response, but it would need discussion if the extrinsic response were retained.","section":"Sec. II.A"}],"recommendation":"reject","confidential_remarks":"To the editor: The manuscript contains a clear symmetry-based error in the extrinsic second-order spin current. The intrinsic results may be publishable as a shorter paper after removing the extrinsic claim, but as submitted the headline CPSC result is not valid. I recommend rejection. I see no indication of misconduct; the error appears to be an oversight in the numerical integration."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Fair warning up front: the paper's headline extrinsic result is wrong, but the intrinsic half is real. The second-order response in Eq. (35) is the integral of Ξ^l = Ω S^l against ∂f/∂k. Under time reversal in this model, d(-k) = -d(k), so ε is even, Ω and S^l are odd, making Ξ^l even. The Fermi function derivative is odd, so the integrand is odd over the full BZ: the extrinsic SOSC is identically zero. The nonzero collinear curve in Fig. 4(b) cannot come from the stated equations. The paper's own mirror-symmetry plots make the same point: the in-plane Ξ components are odd under one mirror while ∂f/∂k_x is even, so they vanish before the T-sum. So the abstract's 'both intrinsic and extrinsic in-plane spin currents' and the CPSC application are unsupported.\n\nWhat is worth keeping: the linear response is correctly shown to be purely out-of-plane, and the intrinsic second-order response is field-polarized along x, with the symmetry selection rules verified by the parity of Eqs. (23), (24), (29). The numerical estimate of ~0.4 e µm/V for the intrinsic conductivity is within the range of other predictions. The application of the conserved-spin-current formalism to the Fu model is new, and the presentation is clear.\n\nSoft spots beyond the fatal one: Appendix B has dimensionally inconsistent intermediate steps (the derivative expressions mix k and dimensionless k). The relaxation-time approximation for surface states without surface-bulk coupling is a simplification, but that is secondary. The self-citation to Ref. 45 is fine—it is the companion formalism paper.\n\nIn short, this is a paper with one good result and one false claim. The intrinsic second-order spin Hall effect in Bi2Te3 surface states is a legitimate contribution that could stand alone. The extrinsic collinear channel, however, is not there. Given the solid intrinsic result, the paper deserves a serious referee rather than a desk rejection. But the current version is not publishable; the authors need to drop the extrinsic claim and resubmit as a study of the intrinsic second-order spin Hall effect. If they cannot do that, reject.","headline":"The intrinsic second-order spin current in Bi2Te3 surface states is a real result, but the extrinsic collinear spin current vanishes by time-reversal symmetry, so the paper's main new claim is unsupported.","tokens_in":12869,"tokens_out":4094,"would_cite":false,"duration_ms":43654,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"On warped topological-insulator surfaces, in-plane pure spin currents appear at second order in the applied electric field, with directions fixed by crystal symmetry.","keywords":["conserved spin current","nonlinear spin Hall effect","spin-weighted Berry curvature","Berry connection polarizability","Berry curvature dipole","hexagonal warping","topological insulator surface states","time-reversal symmetry"],"falsifier":"Compute the Fermi-surface integral in Eq. (35) at a single Fermi energy with the paper's own Ξ^l(k). The extrinsic second-order conductivity survives only if this integral is nonzero; a second-harmonic spin-current measurement on Bi2Te3 can then check whether the predicted collinear component actually appears.","tokens_in":11712,"feed_emoji":"⚡","tokens_out":10266,"duration_ms":119770,"temperature":0.7,"pith_summary":"The paper claims that on the warped surface of a three-dimensional topological insulator such as Bi2Te3, pure spin currents that are absent in the linear response appear at second order in an applied electric field. Working within the conserved-spin-current formalism, it derives linear and nonlinear spin Hall conductivities for the Dirac surface Hamiltonian with hexagonal warping, whose k^3 warping term gives a finite Berry curvature without breaking time-reversal symmetry. The second-order response decomposes into intrinsic and extrinsic channels: the intrinsic channel is governed by the Berry connection polarizability dipole and a spin-dependent Berry-curvature term, while the extrinsic channel comes from the Fermi-surface dipole of the spin-weighted Berry curvature. The C3v anisotropy of the dispersion acts as a selection rule, leaving a field-polarized x-spin current in the intrinsic channel and a collinearly polarized y-spin current in the extrinsic channel. A sympathetic reader would care because these nonlinear in-plane spin currents could enable pure spin injection and field-free magnetization switching without an external magnetic field.","feed_headline":"In-plane pure spin currents turn on at second order in Bi2Te3","feed_subtitle":"Theory predicts both field-polarized and collinear pure spin currents, detectable as edge spin accumulation.","key_machinery":"The central object is the spin-weighted Berry curvature, Ξ^l_ab = Ω^ab S^l, the product of the band Berry curvature and the spin polarization. At linear order it produces the out-of-plane spin Hall conductivity; at second order its Fermi-surface dipole generates the extrinsic in-plane spin current. The intrinsic second-order channel is carried by the Berry connection polarizability (BCP) dipole and by a Berry-curvature/spin-matrix-element term that arise from first-order corrections to the Bloch states. The enabling ingredient is the k^3 hexagonal warping term in the surface Hamiltonian, which makes the Berry curvature finite and anisotropic while leaving time-reversal symmetry intact.","core_discovery":"On the paper's own terms, the central discovery is that the second-order response of the conserved spin current is not empty for time-reversal-invariant hexagonal-warped surface states. The linear spin Hall conductivity is purely out-of-plane; the in-plane spin components vanish not because of time-reversal symmetry alone, but because the C3 anisotropy makes their momentum integrals cancel. At order E^2, the same anisotropy selects nonzero components: the intrinsic conductivity Γ^x is polarized along the applied field, and the extrinsic Fermi-surface conductivity Γ^y is collinear with the current. Both are expressed through band-geometric quantities—spin-weighted Berry curvature, Berry conne","pith_inferences":["Extension: strain or a tilted surface that lowers the C3v symmetry to C_s should unlock additional spin-conductivity tensor components, since the selection rules derived here are purely symmetry-driven.","Extension: the same spin-weighted-Berry-curvature machinery could be applied to other hexagonal-warped Dirac materials, such as Bi2Se3 or Sb2Te3, to predict the sign and magnitude of their nonlinear spin currents without new formalism.","Extension: a first-principles band-structure calculation of Bi2Te3's surface states could test whether the low-energy C3v model captures the dominant contributions or whether higher-order warping terms change the surviving components."],"forward_implications":["Linear response on these surfaces carries only out-of-plane spin polarization; in-plane spin currents are a second-order phenomenon, so they appear at twice the driving frequency in ac experiments.","The intrinsic second-order spin conductivity is predicted to be on the order of 0.4 e µm/V, large enough to produce measurable edge spin accumulation or spin-orbit torque.","The extrinsic contribution is a Fermi-surface effect and is much smaller than the intrinsic one, so it should be distinguishable by its relaxation-time scaling.","The C3v symmetry acts as a selection switch: among the three spin components, only the components even under the mirror symmetries survive the momentum integration.","A collinear spin current, with spins parallel to the travel direction, emerges from the extrinsic channel and could apply field-free torques to an adjacent magnet."],"supporting_citations":[{"why":"Establishes the conserved spin-current formalism whose k-resolved spin current the paper extends to second order.","marker":"[14]"},{"why":"Shows the Berry curvature is the source of the proper conserved spin current, the object carried into the nonlinear regime.","marker":"[17]"},{"why":"Introduces the hexagonal-warping k^3 term used in the surface Hamiltonian; this term is what makes the Berry curvature finite without breaking time-reversal symmetry.","marker":"[25]"},{"why":"Supplies the low-energy model Hamiltonian for Bi2Te3 surface states that the paper starts from.","marker":"[21]"},{"why":"Shows the charge Hall response in these surfaces is suppressed up to third order, which motivates seeking a second-order pure spin current.","marker":"[24]"},{"why":"Defines the intrinsic nonlinear spin Hall effect and the collinear spin-current concept used to classify the extrinsic channel.","marker":"[37]"}],"fun_headline_variants":["Second-order spin currents turn in-plane in Bi2Te3","Hexagonal warping drives nonlinear pure spin currents","Bi2Te3: in-plane spin currents emerge at second order","Beyond linear: conserved spin transport goes transverse","Nonlinear effect yields in-plane spin currents in TIs"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The central calculation depends on the spin-weighted Berry curvature having a nonzero Fermi-surface dipole; if symmetry forces that dipole to vanish, the extrinsic second-order spin current disappears and only the intrinsic terms remain.","fun_headline_variants_meta":{"raw":{"variants":["Second-order spin currents turn in-plane in Bi2Te3","Hexagonal warping drives nonlinear pure spin currents","Bi2Te3: in-plane spin currents emerge at second order","Beyond linear: conserved spin transport goes transverse","Nonlinear effect yields in-plane spin currents in TIs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000141,"raw_usage":{"total_tokens":986,"prompt_tokens":714,"completion_tokens":272,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":458,"completion_tokens_details":{"reasoning_tokens":193}},"tokens_in":458,"tokens_out":272,"duration_ms":3769,"temperature":1.0,"reasoning_tokens":193,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T21:06:23.110666+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the Fermi-surface integral in Eq. (35) at a single Fermi energy with the paper's own Ξ^l(k). The extrinsic second-order conductivity survives only if this integral is nonzero; a second-harmonic spin-current measurement on Bi2Te3 can then check whether the predicted collinear component actually appears.","supporting_citations":[],"review_version":1}