{"id":"b5fd5f48-ac9b-4454-9e60-e54c8cee826a","arxiv_id":"2607.26311","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Ordered BiSbTeSe2 in the Se–Bi–Se–Sb–Te stacking is predicted to be a strong topological insulator with a giant bulk Rashba spin splitting (α_CB ≈ 2.67 eV·Å) approaching that of BiTeI.","lead":"Using computer simulations, this paper predicts that a chemically ordered form of the topological insulator BiSbTeSe2 develops a large spin-split bulk band structure while still hosting protected surface states. The work suggests a single material could combine bulk spin–momentum locking with topological surface transport, relevant for spintronics and spin–charge conversion.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central prediction is credible, but it hinges on real BiSbTeSe2 adopting the ordered Se–Bi–Se–Sb–Te R3m phase; the paper offers only total-energy comparison and Raman evidence, not long-range structural confirmation.","rationale":"The reader's weakest assumption identifies the same structural realizability concern, and I agree that this is the most load-bearing point. The paper's internal calculations — Wannier interpolation, Wilson-loop Z2, surface spectral functions, and the stable fitting-window analysis for α_CB — are credible and mutually consistent. No internal contradiction in the k·p extraction or the topological classification undermines the central claim. The single point at which the argument can fail without a numerical error is whether the ordered R3m phase is what is actually meant by 'BiSbTeSe2' in the experimental record. The paper cites one Raman/phonon study for energetic favorability, but that is indirect and does not demonstrate long-range order. Therefore the reader's CONDITIONAL verdict is appropriate; my read does not change it.","tokens_in":15607,"tokens_out":6821,"duration_ms":78503,"concrete_test":"Perform single-crystal X-ray diffraction and atomic-resolution STEM on the same BSTS crystals used for transport/ARPES and search for the superstructure reflections / atomic-column ordering corresponding to the Se–Bi–Se–Sb–Te R3m stacking. If no long-range polar order is observed (or if only short-range domains are present), the bulk Rashba splitting and topological-surface-state coexistence are not properties of real BiSbTeSe2, and the paper should be read purely as a prediction for a hypothetical ordered phase.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the identification of 'BiSbTeSe2' with the ordered Se–Bi–Se–Sb–Te R3m structure. The electronic-structure part of the paper is internally strong: the α_CB plateau over kmax (Fig. 5c), the necessity of the cubic term (Fig. 5d), the Wannier interpolation (Fig. 6), Wilson-loop ν0=1, and Dirac surface states on both terminations. But all of this applies to a specific ordered polymorph. Section III.A chooses this polymorph by comparing three single-QL stacking sequences and citing Ref. [10] for its energetic favorability. Ref. [10] is a phonon/Raman study; Raman probes local coordination and can be consistent with short-range order, but it does not establish a long-range polar R3m lattice across the crystals used in transport/ARPES. No diffraction, STEM, or other direct structural evidence is cited. The experiments quoted in the Introduction concern conventionally grown BSTS, which is usually a Bi/Sb and Te/Se solid solution. If actual samples are disordered or contain only nanometer-scale ordered domains, inversion symmetry is restored on average, the bulk Rashba splitting cancels, and the 'same material' platform does not exist. This is an external-validity risk, not an internal numerical inconsistency.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports first-principles DFT calculations on a proposed ordered polymorph of BiSbTeSe2 with layer sequence Se–Bi–Se–Sb–Te in space group R3m. The authors find that this structure lacks inversion symmetry, reducing the point group at Γ from D3d to C3v, and that the bulk conduction and valence band edges acquire Rashba-type spin splitting. They fit the spin-split Kramers doublets to a symmetry-constrained two-band k·p Hamiltonian and extract α_CB = 2.666 ± 0.005 eV Å and α_VB = 0.345 ± 0.005 eV Å, with a negative cubic radial correction for the conduction band. Using Wannier interpolation, Wilson loops, and semi-infinite surface spectral functions, they argue that the same ordered phase remains a strong topological insulator, with Te- and Se-terminated surfaces hosting Dirac-like surface states at different energies. The central claim is that sublattice ordering offers a bulk-Rashba topological-insulator platform in an already bulk-insulating tetradymite material.","tokens_in":15911,"tokens_out":10338,"duration_ms":98641,"significance":"The computational core of the paper is credible and internally consistent. The α_CB extraction is supported by a plateau over fitting radii (Fig. 5c), by a large reduction in RMS error when cubic terms are included (Fig. 5d), and by fitting-window checks in the SM. The spin textures are calculated directly from fully relativistic wavefunctions, and the topological diagnostics are independent of the k·p fitting, so there is no circularity. If the ordered R3m phase can be realized, the predicted coexistence of a giant bulk Rashba splitting with a topological surface state would be of significant interest for spin transport. However, the physical relevance of the entire prediction is tied to the experimental realization of long-range polar order, which the manuscript does not establish. The paper would be strengthened by either direct structural evidence or a clear re-scoping as a prediction for a specific ordered polymorph.","major_comments":[{"comment":"The identification of BiSbTeSe2 with the ordered Se–Bi–Se–Sb–Te R3m phase is load-bearing for the paper’s central claim. The evidence offered is a total-energy comparison among three layer sequences plus a Raman/phonon study. Raman spectroscopy is sensitive to local coordination and does not prove long-range polar order; conventional BiSbTeSe2 is typically described as a Bi/Sb and Te/Se solid solution in an averaged structure. If real samples are disordered or contain only nanoscale ordered domains, inversion symmetry is restored on average and the predicted bulk Rashba splitting and the “same material” platform no longer exist. The authors should either provide direct experimental structural evidence (diffraction, STEM, or equivalent) or explicitly re-scope the abstract, title, and conclusions to present this as a prediction for a specific ordered polymorph, with a candid discussion of","section":"Section III.A and Ref. [10]"},{"comment":"The strong topological invariant is asserted in a single sentence: “The resulting indices have strong component ν0 = 1.” Because the R3m structure lacks inversion symmetry, parity-based classification is unavailable, so the Wilson-loop data are the primary evidence for the strong topological phase. The manuscript should show the Wilson-loop (Wannier-charge-center) evolution or provide equivalent numerical data, including the number of occupied bands used and the convergence with respect to the k-mesh. Without this, the topological claim is not independently checkable from the manuscript.","section":"Section IV.B"},{"comment":"The quoted uncertainty α_CB = 2.666 ± 0.005 eV Å is based on window-to-window variation, but the fitted value also depends on the order of the k·p expansion. The same data set yields α_linear = 2.428 eV Å in the linear-only model. While the RMS drop from 3.04 meV to 0.21 meV justifies the cubic model, the authors should report the fit-covariance uncertainty and a sensitivity test to higher-order terms (e.g., k^4/k^5 corrections) so the reader can judge how much of the “intrinsic” coefficient is model-dependent.","section":"Section III.C, Eq. (18)"}],"minor_comments":[{"comment":"The equal-weight radial-annulus procedure is not fully reproducible; please specify the radial bin width and how “nonempty” annuli are defined.","section":"Section III.C.1"},{"comment":"The surface spectral functions use an arbitrary color intensity; adding a colorbar and stating the normalization would improve reproducibility.","section":"Fig. 7 and Fig. S4"},{"comment":"“Available from the author upon reasonable request” is a weak reproducibility statement; depositing relaxed structures, Wannier Hamiltonians, and fitting scripts in a public repository is recommended.","section":"Data Availability"},{"comment":"The sentence introducing Ref. [10] (“Raman spectroscopy and first-principles calculations identify the Se–Bi–Se–Sb–Te sequence as energetically favorable”) could be misread as direct structural confirmation; please phrase it more carefully.","section":"Section I"},{"comment":"For |λ_CB|, consider reporting the largest-window fitted value with its uncertainty (90.1 ± 5.2 eV Å^3 from Fig. S3) rather than the approximate symbol “∼90”.","section":"Table I"}],"recommendation":"major_revision","confidential_remarks":"The paper is a single-author computational study with a sound DFT/k·p core. The main risk is external validity: the ordered R3m phase is chosen by energetics and a Raman study, but no direct experimental evidence is cited for long-range polar order in real BiSbTeSe2. This is not an internal inconsistency, and it can be addressed by explicit re-scoping as a prediction for a specific polymorph plus a discussion of synthesis feasibility. The Wilson-loop section also needs to show its data rather than a one-sentence assertion. I would not reject the manuscript; these issues are fixable within revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know upfront: this is a competent DFT paper whose central number, α_CB ≈ 2.67 eV·Å, is probably right for the ordered Se–Bi–Se–Sb–Te R3m phase. The real question is whether that phase exists in the lab.\n\nWhat's genuinely new: the specific prediction that this particular ordered phase is simultaneously a strong TI and a bulk Rashba system with a coefficient approaching BiTeI. I checked the fitting procedure—it's careful. The plateau in α_CB across fitting windows (2.661–2.670), the clear need for the cubic β term (RMS error drops from 3.04 to 0.21 meV), and the separate handling of λ via k^6 scaling all show the authors know how to extract an intrinsic k→0 coefficient rather than a window-dependent slope. The Wilson-loop ν0=1 and the surface spectral functions on both terminations are consistent and independently computed from the Wannier model. That's the standard toolkit, applied properly.\n\nThe load-bearing assumption is the identification of BiSbTeSe2 with this specific ordered stacking. The paper's evidence is a total-energy comparison among three model stackings and a citation to a phonon/Raman study [10]. Raman and energetics can indicate short-range order, but they don't establish a long-range polar R3m lattice in the crystals used for transport or ARPES. Conventionally grown BSTS is a disordered solid solution; if inversion symmetry is restored on average, the bulk Rashba splitting cancels and the 'same material' platform is gone. The authors acknowledge disorder only as future work, not as a threat to the current claim. That's the honest soft spot, and it's external validity rather than internal inconsistency.\n\nAlso minor: data availability is 'upon reasonable request'—for a computational paper with Wannier Hamiltonians and scripts, that's weaker than shipping them, though not a fatal flaw.\n\nWho's this for? People working on tetradymite topological insulators, Rashba physics, and spin-charge conversion. The paper deserves a serious referee: the electronic-structure analysis is solid and the prediction is crisp. A referee should push on the structural-order question, maybe ask for a discussion of disorder and a realistic assessment of whether long-range order can be stabilized (e.g., growth conditions, substrate, annealing). That's a revision, not a rejection.","headline":"A clean first-principles prediction of coexisting topological surface state and giant bulk Rashba splitting in one ordered polymorph of BiSbTeSe2 — but the paper never establishes that real BSTS actually adopts that ordering.","tokens_in":16453,"tokens_out":1858,"would_cite":true,"duration_ms":18844,"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":"Ordered BiSbTeSe2 is claimed to be a polar strong topological insulator whose chemical stacking produces a bulk Rashba spin splitting close to that of BiTeI while retaining a protected surface state.","keywords":["BiSbTeSe2","topological insulator","bulk Rashba splitting","spin-momentum locking","C3v symmetry","k·p model","Wannier functions","Z2 invariant"],"falsifier":"Measure the near-Γ bulk dispersion of a single-domain BiSbTeSe2 crystal with spin-resolved photoemission: if the conduction-band branches do not split with opposite helicities at a rate compatible with α ≈ 2.67 eV·Å, or if no Dirac-like surface state appears on either Te- or Se-terminated surface, the central claim fails. Similarly, diffraction or atomic-resolution imaging showing that real BiSbTeSe2 lacks long-range Se–Bi–Se–Sb–Te order would remove the premise.","tokens_in":15430,"feed_emoji":"⚛️","tokens_out":6356,"duration_ms":52510,"temperature":0.7,"pith_summary":"BiSbTeSe2 is a tetradymite usually treated with inversion symmetry; this paper claims that when its atomic layers are ordered in the Se–Bi–Se–Sb–Te sequence, the quintuple layer becomes polar. The loss of inversion symmetry (point group at Γ reduces from D3d to C3v) allows a bulk Rashba spin splitting, while the band inversion and surface Dirac cone of a strong topological insulator survive. From first-principles calculations, the paper extracts a conduction-band linear Rashba coefficient of about 2.67 eV·Å—close to the benchmark polar semiconductor BiTeI—and a valence-band coefficient of about 0.35 eV·Å, with sizeable cubic corrections. The appeal is that one material would then host both bulk spin–momentum locking and a protected topological surface channel, making it a platform for studying and tuning bulk and surface spin transport together.","feed_headline":"Giant bulk-Rashba splitting in polar topological insulator BiSbTeSe2","feed_subtitle":"Layer ordering gives a tetradymite a giant bulk spin splitting without losing its protected surface state.","key_machinery":"The load-bearing object is the ordered Se–Bi–Se–Sb–Te quintuple layer, whose chemical inequivalence on its two sides removes the inversion center while preserving the threefold rotation and vertical mirrors, reducing the Γ point group from D3d to C3v. The quantitative workhorse is a symmetry-constrained two-band k·p Hamiltonian applied separately to the conduction- and valence-band Kramers doublets, H_D = ε_D(k) σ_0 + (α_D + β_D k²)(k_x σ_y − k_y σ_x) + λ_D(3k_x² k_y − k_y³) σ_z. This form separates the intrinsic linear Rashba coefficient α_D from the isotropic cubic correction β_D and the anisotropic C3v warping λ_D. Fitting the DFT half-splitting d_D(k) directly to this Hamiltonian—rather","core_discovery":"The central claim is that chemical sublattice ordering in BiSbTeSe2—the Se–Bi–Se–Sb–Te stacking—breaks inversion symmetry in the bulk tetradymite, and that this noncentrosymmetric polar phase is simultaneously a strong topological insulator and a bulk Rashba system. The evidence presented includes bulk band inversion with a direct gap near 0.54 eV; a Z2 index ν0 = 1 from Wilson-loop evolution; a Dirac-like surface state on both Te- and Se-terminated surfaces; and spin-resolved band structures whose conduction- and valence-band doublets split linearly away from Γ with opposite helicities. Fitting the DFT dispersion near Γ to a symmetry-constrained two-band k·p Hamiltonian yields α_CB = 2.666","pith_inferences":["If ordered BiSbTeSe2 can be grown with sufficient domain size, spin- and angle-resolved photoemission should resolve the predicted conduction-band doublet with opposite helicities; the sin(3φ) out-of-plane spin modulation offers a sharp fingerprint of the cubic C3v term.","The near-BiTeI coupling suggests the bulk Edelstein effect could be sizable in this material; measurements of current-induced spin polarization might show a bulk contribution that is distinguishable from the topological surface channel by its dependence on Fermi-level position.","The strength of the claim rests on the physical realization of long-range Se–Bi–Se–Sb–Te order; if disorder or competing stackings dominate, the effective bulk splitting in real crystals could be much smaller than the ordered-phase value.","A natural extension is to test whether strain or gating can move the Fermi level across the Dirac crossing and the Rashba-split bands, allowing controlled switching between surface-dominated and bulk-Rashba-dominated transport in one sample."],"forward_implications":["If the ordered polar phase is realized, bulk spin–momentum locking and topological surface states coexist at ambient conditions in a single bulk-insulating tetradymite.","The conduction-band coefficient α_CB ≈ 2.67 eV·Å places ordered BiSbTeSe2 among the strongest bulk-Rashba topological insulators reported, approaching BiTeI.","The negative cubic correction β_CB ≈ −140 eV·Å³ means the effective in-plane Rashba coupling decreases by roughly 14% at k ≈ 0.05 Å⁻¹, so finite-momentum measurements will see a smaller slope than the intrinsic k→0 value.","Because the Te- and Se-terminated surfaces place the Dirac crossing at different energies (about 0.30 eV and 0.10 eV below the Fermi level), termination choice offers a tuning knob for surface carrier density.","The same ordered structure provides a platform to separate bulk Rashba and topological surface contributions to spin–charge conversion and current-induced spin polarization."],"fun_headline_variants":["Sublattice ordering unlocks giant bulk Rashba in topological insulator","BiSbTeSe2: topological insulator with a giant bulk Rashba split","How layer order gives a tetradymite giant spin splitting","Polar topological insulator now boasts bulk Rashba splitting","Giant bulk Rashba states coexist with topological surface states"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"Real BiSbTeSe2 must actually form the ordered Se–Bi–Se–Sb–Te stacking used in all calculations; if the physical compound is the disordered or centrosymmetric alloy typically reported, the predicted bulk Rashba splitting and the same-material coexistence would not occur.","fun_headline_variants_meta":{"raw":{"variants":["Sublattice ordering unlocks giant bulk Rashba in topological insulator","BiSbTeSe2: topological insulator with a giant bulk Rashba split","How layer order gives a tetradymite giant spin splitting","Polar topological insulator now boasts bulk Rashba splitting","Giant bulk Rashba states coexist with topological surface states"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00022,"raw_usage":{"total_tokens":1335,"prompt_tokens":846,"completion_tokens":489,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":590,"completion_tokens_details":{"reasoning_tokens":403}},"tokens_in":590,"tokens_out":489,"duration_ms":5040,"temperature":1.0,"reasoning_tokens":403,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T00:10:56.480714+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the near-Γ bulk dispersion of a single-domain BiSbTeSe2 crystal with spin-resolved photoemission: if the conduction-band branches do not split with opposite helicities at a rate compatible with α ≈ 2.67 eV·Å, or if no Dirac-like surface state appears on either Te- or Se-terminated surface, the central claim fails. Similarly, diffraction or atomic-resolution imaging showing that real BiSbTeSe2 lacks long-range Se–Bi–Se–Sb–Te order would remove the premise.","supporting_citations":[],"review_version":1}