Pith. sign in

REVIEW 3 major objections 5 minor 49 references

Giant Bulk-Rashba Splitting in Polar Topological Insulator BiSbTeSe$_2$

T0 review · 3 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash

Pith's one-line read 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.

desk verdict 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. read the letter →

arxiv 2607.26311 v2 pith:ZLDGD4RO submitted 2026-07-28 cond-mat.mes-hall cond-mat.mtrl-sci

classification cond-mat.mes-hallcond-mat.mtrl-sci
keywords BiSbTeSe2topologicalinsulatorbulkRashbasplittingspin-momentumlockingC3vsymmetryk·pmodelWannierfunctionsZ2invariant
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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

What would settle it

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.

Watch

Extended reading notes

Core claim

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

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

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.

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 (3)
  1. [Section III.A and Ref. [10]] 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
  2. [Section IV.B] 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.
  3. [Section III.C, Eq. (18)] 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.
minor comments (5)
  1. [Section III.C.1] The equal-weight radial-annulus procedure is not fully reproducible; please specify the radial bin width and how “nonempty” annuli are defined.
  2. [Fig. 7 and Fig. S4] The surface spectral functions use an arbitrary color intensity; adding a colorbar and stating the normalization would improve reproducibility.
  3. [Data Availability] “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.
  4. [Section I] 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.
  5. [Table I] 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”.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the Rashba parameters are explicit fits to DFT band splittings, and the topological indices come from independent Wannier/Wilson-loop calculations.

full rationale

The central claims are self-contained against first-principles data. The Rashba coefficients α_CB and α_VB are obtained by least-squares fitting the DFT half-splitting d_D,DFT(k) to the symmetry-constrained two-band k·p Hamiltonian (Eqs. 10–14), with explicit kmax-plateau analysis (Fig. 5c) and comparison to linear-only fits; this is parameter extraction, not a prediction that assumes its own target. The topological classification is established by Wannier interpolation, Wilson-loop Z2 indices, and semi-infinite surface spectral functions derived from the same DFT Hamiltonian, with no target value inserted. The C3v model form follows from symmetry and is used to interpret, not to generate, the DFT spin textures. The only self-citation is Ref. [9], used to support the empirical context that BSTS is already used for transport studies; it is not load-bearing for any derived result. The choice of the ordered Se–Bi–Se–Sb–Te R3m structure is supported by the paper's own total-energy comparison and by independent phonon/Raman work (Ref. [10]). Whether real samples realize long-range polar order is an external-validity question, not a circularity in the derivation.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The central claim rests on standard DFT approximations and on the assumed realization of a specific ordered stacking; the Rashba coefficients are fitted characterizations rather than parameter-free derivations. No new physical entities are introduced.

free parameters (4)
  • α_CB (conduction-band linear Rashba coefficient) = 2.666 ± 0.005 eV·Å
    Fitted to the DFT conduction-band half-splitting using the C3v two-band model over k_max windows 0.015–0.050 Å^-1; window-to-window variation quoted as systematic uncertainty.
  • β_CB (isotropic cubic radial Rashba correction) = -140 ± 10 eV·Å^3
    Fitted simultaneously with α_CB; required to reproduce the reduction of effective splitting at finite momentum.
  • |λ_CB| (cubic C3v warping coefficient) = ~90 eV·Å^3 (finite-window estimate)
    Fitted from λ²k^6 contribution; statistically unresolved in small windows and reported only as a finite-window estimate.
  • α_VB (valence-band linear Rashba coefficient) = 0.345 ± 0.005 eV·Å
    Fitted to the DFT valence-band half-splitting over intermediate windows 0.018–0.030 Å^-1; cubic coefficients for VB do not converge.
assumptions (5)
  • domain assumption The Se–Bi–Se–Sb–Te ordered stacking is the structure realized in bulk BiSbTeSe2.
    The paper selects this sequence because it is lowest in energy among three considered and cites ref [10]; no experimental determination of long-range cation order in the bulk is cited. This is the load-bearing structural premise.
  • domain assumption PBE-GGA with DFT-D3 and fully relativistic pseudopotentials accurately describes band inversion, gaps, and Rashba splitting.
    Standard approximation; no hybrid-functional or experimental band-gap benchmark is provided, and the gap is noted to be underestimated.
  • domain assumption The two-band C3v Hamiltonian of Eq. (4), with only α, β and λ, fully describes the near-Γ conduction- and valence-band doublets.
    Fit quality (RMS 0.21 meV for CB) supports sufficiency in the fitted window, but any omitted coupling to other bands could bias the extracted α.
  • domain assumption Time-reversal symmetry holds; the system is nonmagnetic.
    Used in Eq. (1) to relate E(k) and spin at ±k, and to classify Kramers doublets.
  • domain assumption The Wannier projection onto the p-orbital manifold captures the low-energy physics, so Wilson-loop and surface Green-function results are valid.
    Validated by Fig. 6 agreement between Wannier-interpolated and DFT bands near the gap.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Giant Bulk-Rashba Splitting in Polar Topological Insulator BiSbTeSe$_2$." pith.science (2026). https://pith.science/paper/ZLDGD4RO

@misc{pith2026260726311,
  author       = {Pith},
  title        = {Pith review of: Giant Bulk-Rashba Splitting in Polar Topological Insulator BiSbTeSe$_2$},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZLDGD4RO}},
  note         = {Machine review of arXiv:2607.26311}
}
abstract

Bulk-Rashba spin splitting is forbidden in tetradymite topological insulators like Bi$_2$Se$_3$ or Bi$_2$Te$_3$, since their quintuple-layer stacking preserves inversion symmetry. We show that BiSbTeSe$_2$ escapes this restriction: in the Se-Bi-Se-Sb-Te sequence, the structure loses its inversion center, reducing the point group symmetry at $\Gamma$ from $D_{3d}$ to $C_{3v}$. First-principles density functional calculations with spin-orbit coupling show that this ordered structure retains bulk band inversion and a linearly dispersive surface state of a strong topological insulator. Additionally, its bulk bands acquire a pronounced linear-in-$k$ spin splitting away from $\Gamma$. Fitting the conduction- and valence-band doublets to symmetry-constrained two-band $k\cdot p$ Hamiltonians, we extract intrinsic linear Rashba coefficients of $\alpha_{\mathrm{CB}}\approx2.66~\mathrm{eV\,\text{\AA}}$ and $\alpha_{\mathrm{VB}}\approx0.35~\mathrm{eV\,\text{\AA}}$. The conduction-band value places ordered BiSbTeSe$_2$ among the strongest bulk-Rashba topological-insulator systems reported to date and approaches the coupling found in the benchmark polar Rashba semiconductor BiTeI. Sublattice ordering thus provides a route to giant bulk spin--momentum locking that coexists with protected topological surface states, offering a platform in which bulk-Rashba and topological surface contributions to spin and charge transport can be investigated within the same material.

Figures

Figures reproduced from arXiv: 2607.26311 by the authors.

Figure 1
Figure 1. FIG. 1: Crystal structure and bulk Brillouin zone of ordered BiSbTeSe [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Spin-resolved electronic band structure of ordered BiSbTeSe [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Calculated physical-spin texture of the two Rashba-split conduction-band branches in the momentum plane [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: FIG. 4: Calculated physical-spin texture of the two Rashba-split valence-band branches in the momentum plane [PITH_FULL_IMAGE:figures/full_fig_p006_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5: Conduction-band fit to the symmetry-constrained [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]
Figure 7
Figure 7. Figure 7: FIG. 7: Surface spectral function of semi-infinite [PITH_FULL_IMAGE:figures/full_fig_p009_7.png]
Figure 6
Figure 6. Figure 6: FIG. 6: Comparison between the first-principles DFT [PITH_FULL_IMAGE:figures/full_fig_p009_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

49 extracted references

  1. [10]

    Y. L. Chen, J. G. Analytis, J.-H. Chu, Z. K. Liu, S.-K. Mo, X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, S. C. Zhang, I. R. Fisher, Z. Hussain, and Z.-X. Shen, Ex- perimental realization of a three-dimensional topological insulator, Bi2Te3,Science325, 178 (2009)

  2. [1]

    Since any linear combination of the two degenerate states is also an allowed eigenstate, an individual branch helicity is not uniquely defined exactly at Γ

    Conduction-band doublet At Γ, the two conduction-band states form a degener- ate Kramers pair. Since any linear combination of the two degenerate states is also an allowed eigenstate, an individual branch helicity is not uniquely defined exactly at Γ. At finite in-plane momentum, however, the degen- eracy is lifted and the lower and upper energy-ordered b...

  3. [2]

    Valence-band doublet The two valence-band branches display the same qual- itative signatures of Rashba spin–momentum locking as the conduction-band (see Fig. 4). Close to Γ, their in- plane spin components wind approximately tangentially around the momentum-space origin, with opposite he- licities for the two energy-ordered branches. The spin component pa...

  4. [3]

    (10) using nonlinear least squares

    Fitting procedure At each momentum point, the two DFT bands belong- ing to a given doublet were ordered by energy, EDFT D,− (k)< EDFT D,+ (k).(12) They were then decomposed into a band-center energy, ED,DFT(k) = EDFT D,+ (k) +E DFT D,− (k) 2 ,(13) and a half-splitting, dD,DFT(k) = EDFT D,+ (k)−E DFT D,− (k) 2 .(14) The Rashba parametersα D,β D, and|λ D|we...

  5. [4]

    Conduction-band doublet For the conduction-band doublet, the fitting-window dependence was examined over 0.015≤k max ≤0.050 ˚A −1 .(17) The fitted linear coefficient remains between 2.661 and 2.670 eV ˚A throughout this interval. Taking the window- to-window variation as the dominant systematic uncer- tainty gives αCB = 2.666±0.005 eV ˚A.(18) The weak dep...

  6. [5]

    Since the valence splitting is much smaller than the conduction-band splitting, a denser momentum grid was used in a smaller region around Γ

    Valence-band doublet The same symmetry-constrained Hamiltonian was ap- plied independently to the valence-band-edge doublet. Since the valence splitting is much smaller than the conduction-band splitting, a denser momentum grid was used in a smaller region around Γ. The fitting-window dependence was examined over 0.008≤k max ≤0.030 ˚A −1 .(24) For the sma...

  7. [6]

    M. Z. Hasan and C. L. Kane, Colloquium: Topological insulators,Rev. Mod. Phys.82, 3045 (2010)

  8. [7]

    Qi and S.-C

    X.-L. Qi and S.-C. Zhang, Topological insulators and su- perconductors,Rev. Mod. Phys.83, 1057 (2011). 11

Show all 49 references
  1. [8]

    Zhang, C.-X

    H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. Zhang, Topological insulators in Bi 2Se3, Bi 2Te3, and Sb2Te3 with a single Dirac cone on the surface,Nat. Phys.5, 438 (2009)

  2. [9]

    Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, R. J. Cava, and M. Z. Hasan, Observation of a large-gap topological-insulator class with a single Dirac cone on the surface,Nat. Phys. 5, 398 (2009)

  3. [11]

    Z. Ren, A. A. Taskin, S. Sasaki, K. Segawa, and Y. Ando, Large bulk resistivity and surface quantum oscillations in the topological insulator Bi 2Te2Se,Phys. Rev. B82, 241306 (2010)

  4. [12]

    Ando, Topological insulator materials,J

    Y. Ando, Topological insulator materials,J. Phys. Soc. Jpn.82, 102001 (2013)

  5. [13]

    Lohani, P

    H. Lohani, P. Mishra, A. Banerjee, K. Majhi, R. Gane- san, U. Manju, D. Topwal, P. S. Anil Kumar, and B. R. Sekhar, Band structure of topological insulator BiSbTe1.25Se1.75,Sci. Rep.7, 4567 (2017)

  6. [14]

    Shyaga, R

    N. Shyaga, R. Sarkar, L. Nanda, R. Chakraborty, S. K. Negi, A. Baral, J. Rajendran, S. Narasimhan, and D. Suri, Multi-channel conduction in electron transport dynamics of three-dimensional topological insulators,J. Phys.: Condens. Matter37, 275701 (2025)

  7. [15]

    German, E

    R. German, E. V. Komleva, P. Stein, V. G. Mazurenko, Z. Wang, S. V. Streltsov, Y. Ando, and P. H. M. van Loosdrecht, Phonon mode calculations and Raman spectroscopy of the bulk-insulating topological insulator BiSbTeSe2,Phys. Rev. Mater.3, 054204 (2019)

  8. [16]

    Ishizaka, M

    K. Ishizaka, M. S. Bahramy, H. Murakawa, M. Sakano, T. Shimojima, T. Sonobe, K. Koizumi, S. Shin, H. Miya- hara, A. Kimura,et al., Giant Rashba-type spin splitting in bulk BiTeI,Nat. Mater.10, 521 (2011)

  9. [17]

    M. S. Bahramy, B.-J. Yang, R. Arita, and N. Nagaosa, Emergence of non-centrosymmetric topological insulat- ing phase in BiTeI under pressure,Nat. Commun.3, 679 (2012)

  10. [18]

    Y. L. Chen, M. Kanou, Z. K. Liu, H. J. Zhang, J. A. Sobota, D. Leuenberger, S. K. Mo, B. Zhou, S.-L. Yang, P. S. Kirchmann, D. H. Lu, R. G. Moore, Z. Hussain, Z. X. Shen, X. L. Qi, and T. Sasagawa, Discovery of a single topological Dirac fermion in the strong inversion asymmet...

  11. [19]

    Landolt, S

    G. Landolt, S. V. Eremeev, Y. M. Koroteev, B. Slom- ski, S. Muff, T. Neupert, M. Kobayashi, V. N. Stro- cov, T. Schmitt, Z. S. Aliev, M. B. Babanly, I. R. Ami- raslanov, E. V. Chulkov, J. Osterwalder, and J. H. Dil, Disentanglement of Surface and Bulk Rashba Spin Split- tings ...

  12. [20]

    J.-J. Zhou, W. Feng, Y. Zhang, S. A. Yang, and Y. Yao, Engineering Topological Surface States and Gi- ant Rashba Spin Splitting in BiTeI/Bi 2Te3 Heterostruc- tures,Sci. Rep.4, 3841 (2014)

  13. [21]

    Di Sante, P

    D. Di Sante, P. Barone, R. Bertacco, and S. Picozzi, Elec- tric Control of the Giant Rashba Effect in Bulk GeTe, Adv. Mater.25, 509 (2013)

  14. [22]

    Liebmann, C

    M. Liebmann, C. Rinaldi, D. Di Sante, J. Kellner, C. Pauly, R. N. Wang, J. E. Boschker, A. Gius- sani, S. Bertoli, M. Cantoni, L. Baldrati, M. Asa, I. Vobornik, G. Panaccione, D. Marchenko, J. S´ anchez- Barriga, O. Rader, R. Calarco, S. Picozzi, R. Bertacco, and M. Morgenster...

  15. [23]

    Plekhanov, P

    E. Plekhanov, P. Barone, D. Di Sante, and S. Picozzi, En- gineering relativistic effects in ferroelectric SnTe,Phys. Rev. B90, 161108 (2014)

  16. [24]

    T. H. Hsieh, H. Lin, J. Liu, W. Duan, A. Bansil, and L. Fu, Topological crystalline insulators in the SnTe ma- terial class,Nat. Commun.3, 982 (2012)

  17. [25]

    Giannozziet al., QUANTUM ESPRESSO: A modu- lar and open-source software project for quantum simula- tions of materials,J

    P. Giannozziet al., QUANTUM ESPRESSO: A modu- lar and open-source software project for quantum simula- tions of materials,J. Phys.: Condens. Matter21, 395502 (2009)

  18. [26]

    Giannozziet al., Advanced capabilities for materials modelling with QUANTUM ESPRESSO,J

    P. Giannozziet al., Advanced capabilities for materials modelling with QUANTUM ESPRESSO,J. Phys.: Con- dens. Matter29, 465901 (2017)

  19. [27]

    Vanderbilt, Soft self-consistent pseudopotentials in a generalized eigenvalue formalism,Phys

    D. Vanderbilt, Soft self-consistent pseudopotentials in a generalized eigenvalue formalism,Phys. Rev. B41, 7892 (1990)

  20. [28]

    Dal Corso, Pseudopotentials periodic table: From H to Pu,Comput

    A. Dal Corso, Pseudopotentials periodic table: From H to Pu,Comput. Mater. Sci.95, 337 (2014)

  21. [29]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple,Phys. Rev. Lett.77, 3865 (1996)

  22. [30]

    Grimme, J

    S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, A con- sistent and accurateab initioparametrization of density functional dispersion correction (DFT-D) for the 94 ele- ments H–Pu,J. Chem. Phys.132, 154104 (2010)

  23. [31]

    C. G. Broyden, J. E. Dennis, Jr., and J. J. Mor´ e, On the local and superlinear convergence of quasi-Newton methods,IMA J. Appl. Math.12, 223 (1973)

  24. [32]

    Fletcher, A new approach to variable metric algo- rithms,Comput

    R. Fletcher, A new approach to variable metric algo- rithms,Comput. J.13, 317 (1970)

  25. [33]

    Goldfarb, A family of variable-metric methods derived by variational means,Math

    D. Goldfarb, A family of variable-metric methods derived by variational means,Math. Comp.24, 23 (1970)

  26. [34]

    D. F. Shanno, Conditioning of quasi-Newton methods for function minimization,Math. Comp.24, 647 (1970)

  27. [35]

    H. J. Monkhorst and J. D. Pack, Special points for Brillouin-zone integrations,Phys. Rev. B13, 5188 (1976)

  28. [36]

    S. P. Ong, W. D. Richards, A. Jain, G. Hautier, M. Kocher, S. Cholia, D. Gunter, V. L. Chevrier, K. A. Persson, and G. Ceder, Python Materials Genomics (py- matgen): A robust, open-source Python library for ma- terials analysis,Comput. Mater. Sci.68, 314 (2013)

  29. [37]

    Marzari and D

    N. Marzari and D. Vanderbilt, Maximally localized gen- eralized Wannier functions for composite energy bands, Phys. Rev. B56, 12847 (1997)

  30. [38]

    Souza, N

    I. Souza, N. Marzari, and D. Vanderbilt, Maximally lo- calized Wannier functions for entangled energy bands, Phys. Rev. B65, 035109 (2001)

  31. [39]

    A. A. Mostofi, J. R. Yates, Y.-S. Lee, I. Souza, D. Van- derbilt, and N. Marzari, wannier90: A tool for obtaining maximally-localised Wannier functions,Comput. Phys. Commun.178, 685 (2008)

  32. [40]

    A. A. Mostofi, J. R. Yates, G. Pizzi, Y.-S. Lee, I. Souza, D. Vanderbilt, and N. Marzari, An updated version of wannier90: A tool for obtaining maximally-localised Wannier functions,Comput. Phys. Commun.185, 2309 (2014)

  33. [41]

    Q. Wu, S. Zhang, H.-F. Song, M. Troyer, and A. A. Soluyanov, WannierTools: An open-source software package for novel topological materials,Comput. Phys. 12 Commun.224, 405 (2018)

  34. [42]

    Fu and C

    L. Fu and C. L. Kane, Time reversal polarization and a Z2 adiabatic spin pump,Phys. Rev. B74, 195312 (2006)

  35. [43]

    Yu, X.-L

    R. Yu, X.-L. Qi, A. Bernevig, Z. Fang, and X. Dai, Equiv- alent expression ofZ 2 topological invariant for band in- sulators using the non-Abelian Berry connection,Phys. Rev. B84, 075119 (2011)

  36. [44]

    A. A. Soluyanov and D. Vanderbilt, Computing topolog- ical invariants without inversion symmetry,Phys. Rev. B 83, 235401 (2011)

  37. [45]

    M. P. L´ opez Sancho, J. M. L´ opez Sancho, and J. Rubio, Highly convergent schemes for the calculation of bulk and surface Green functions,J. Phys. F: Met. Phys.15, 851 (1985)

  38. [46]

    Fu, Hexagonal warping effects in the surface states of the topological insulator Bi 2Te3,Phys

    L. Fu, Hexagonal warping effects in the surface states of the topological insulator Bi 2Te3,Phys. Rev. Lett.103, 266801 (2009)

  39. [47]

    Vajna, E

    S. Vajna, E. Simon, A. Szilva, K. Palot´ as, B. Uj- falussy, and L. Szunyogh, Higher-order contributions to the Rashba–Bychkov effect with application to the Bi/Ag(111) surface alloy,Phys. Rev. B85, 075404 (2012)

  40. [48]

    Zhang, C

    F. Zhang, C. L. Kane, and E. J. Mele, Surface states of topological insulators,Physical Review B86, 081303 (2012)

  41. [49]

    Ghosh, N

    S. Ghosh, N. Stoji´ c, and N. Binggeli, Rashba-induced spin texture and spin-layer-locking effects in the antifer- romagnetic CrI 3 bilayer,Journal of Physics and Chem- istry of Solids173, 111100 (2023). SUPPLEMENT AL MA TERIAL Crystal Structure and Band Inversion We considere...

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.