pith. sign in

arxiv: 2606.22014 · v1 · pith:FQYT3YZSnew · submitted 2026-06-20 · ❄️ cond-mat.mtrl-sci

Stacking-Directed Polarization and Excitonic Engineering in MoS₂/MoSe₂ van der Waals Heterostructures

Mohammed El Amine Miloudi , Oliver K\"uhn This is my paper

Pith reviewed 2026-06-26 12:04 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords van der Waals heterostructuresMoS2MoSe2stacking orderexcitonspolarizationmany-body perturbation theorytrilayers
0
0 comments X

The pith

Stacking order in MoS2/MoSe2 trilayers controls photogenerated electron location via internal electric fields and 60-70 meV band shifts.

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

The paper establishes that the sequence of atomic layers in MoS2 and MoSe2 structures dictates the direction of internal electric fields and the resulting shifts in electronic band edges. In trilayer structures with two MoS2 layers and one MoSe2, different stackings move photogenerated electrons predictably between the central and bottom MoS2 layers. Calculations show these shifts reach 60 to 70 millielectronvolts and produce a 36 millielectronvolt change in interlayer exciton energy that matches experiments. A sympathetic reader would care because this provides a way to engineer light-matter interactions at the atomic scale just by sliding the layers relative to each other.

Core claim

Using GW+BSE many-body perturbation theory on MoS2/MoSe2 van der Waals heterostructures, the stacking sequence in 2L-MoS2/MoSe2 trilayers enables deterministic control of photogenerated electrons between the central and bottom MoS2 layers. This control is governed by internal electric fields and quasiparticle band-edge shifts of 60-70 meV. The calculations predict a 36 meV interlayer excitonic shift that agrees with recent experiments. Homobilayer MoS2 shows a switchable interlayer dipole from registry-induced symmetry breaking, while the MoS2/MoSe2 hetero-interface remains pinned by the chemical potential mismatch between sulfur and selenium.

What carries the argument

Stacking sequence that induces registry-dependent symmetry breaking and internal electric fields controlling band edges and excitons in the trilayer structures.

If this is right

  • Homobilayer MoS2 exhibits a switchable interlayer dipole driven by registry-induced symmetry breaking.
  • The MoS2/MoSe2 hetero-interface remains pinned by the intrinsic chemical potential mismatch between sulfur and selenium.
  • Trilayer stacking enables deterministic control of photogenerated electrons between MoS2 layers via 60-70 meV band-edge shifts.
  • A 36 meV interlayer excitonic shift is predicted in agreement with experiments.
  • TMD trilayers form a platform for sliding ferroelectricity and programmable optoelectronic functionalities.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same stacking mechanism could be tested in other transition metal dichalcogenide combinations to check for comparable electron control.
  • Mechanical sliding of layers in a device might enable gate-free switching of optical responses without external electric fields.
  • The microscopic link between registry and many-body effects could inform design of thicker multilayer stacks or mixed-material systems for similar polarization engineering.

Load-bearing premise

The GW+BSE many-body perturbation theory calculations on the chosen structural models accurately capture the stacking-dependent polarization and excitonic energies without substantial errors from approximations or incomplete convergence.

What would settle it

Spatially resolved spectroscopy or scanning tunneling microscopy that measures whether photogenerated electrons switch between the central and bottom MoS2 layers when the stacking order is changed in a trilayer sample.

Figures

Figures reproduced from arXiv: 2606.22014 by Mohammed El Amine Miloudi, Oliver K\"uhn.

Figure 1
Figure 1. Figure 1: Top and side views of the atomic configurations for (a) AB-stacked and (b) [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Planar-averaged charge density difference (∆ [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Planar-averaged electrostatic potentials, quasiparticle band structures, and [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Many-body excitonic absorption spectra and stacking-dependent photolumi [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Atomic configurations of 2L-MoS2/MoSe2 heterostructures for different stacking sequences. Panels (a–d) show the ABA, BAA, ABB, and BAB configurations, respectively, in both top and side views. The interlayer distances dint1 and dint2 denote the separations between L1–L2 and L2–L3, highlighting the structural variation induced by stacking order. AAB, and BBA exhibit larger c values (approximately 36.3–36.5 … view at source ↗
Figure 6
Figure 6. Figure 6: Charge redistribution in 2L-MoS2/MoSe2 heterostructures for different stacking sequences. Panels (a–d) show the charge density difference (∆ρ) isosurfaces for ABA, BAA, ABB, and BAB configurations, respectively, where yellow and cyan regions represent electron accumulation and depletion (isovalue: 0.001 e/˚A3 ). Panels (e–h) present the corresponding planar-averaged ∆ρ(z) profiles across the trilayer, reve… view at source ↗
Figure 7
Figure 7. Figure 7: Quasiparticle electronic structures and band alignments of trilayer 2L [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Excitonic absorption of trilayer 2L-MoS2/MoSe2 heterostructures.(e)–(h) Many-body excitonic absorption spectra (black solid lines, Gaussian broadening with width 0.1 eV) and oscillator strengths (blue bars) calculated at the EV G0W0+BSE level. The variations in the low-energy excitonic peaks reflect the modulation of interlayer coupling and spatial charge separation induced by the different stacking regist… view at source ↗
Figure 9
Figure 9. Figure 9: Schematic representation of the energy level shifts for ILXs in AB- and [PITH_FULL_IMAGE:figures/full_fig_p018_9.png] view at source ↗
read the original abstract

The stacking-dependent polarization and excitonic response of MoS$_2$/MoSe$_2$ heterostructures were investigated using GW+BSE many-body perturbation theory. While homobilayer MoS$_2$ exhibited a switchable interlayer dipole driven by registry-induced symmetry breaking, the MoS$_2$/MoSe$_2$ hetero-interface remained pinned by the intrinsic chemical potential mismatch between sulfur and selenium. In 2L-MoS$_2$/MoSe$_2$ trilayers, the stacking sequence enabled a deterministic control of photogenerated electrons between the central and bottom MoS$_2$ layers, governed by internal electric fields and quasiparticle band-edge shifts of 60--70~meV. Our calculations predicted a 36~meV interlayer excitonic shift, in remarkable agreement with recent experiments. These results elucidate the microscopic link between atomic registry and many-body interactions, establishing transition metal dichalcogenide trilayers as a potential platform for sliding ferroelectricity and programmable optoelectronic functionalities.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript applies GW+BSE many-body perturbation theory to MoS2/MoSe2 van der Waals heterostructures, reporting that homobilayer MoS2 exhibits a registry-driven switchable interlayer dipole while the hetero-interface is pinned by S/Se chemical-potential mismatch. In 2L-MoS2/MoSe2 trilayers, stacking sequence is claimed to enable deterministic control of photogenerated electrons between central and bottom MoS2 layers via internal fields and 60-70 meV quasiparticle band-edge shifts; the calculations predict a 36 meV interlayer excitonic shift in agreement with experiment.

Significance. If the numerical results are robust, the work establishes a direct microscopic connection between atomic registry, polarization, and many-body excitonic energies in TMD trilayers, supporting potential applications in sliding ferroelectricity and programmable optoelectronics. The explicit use of GW+BSE for both polarization and excitons, together with the reported experimental agreement, constitutes a concrete strength.

major comments (2)
  1. [Methods] Methods/Computational Details: The reported 60-70 meV band-edge shifts and 36 meV interlayer excitonic shift are presented without convergence data on k-point sampling, vacuum padding, or sensitivity to GW approximations (plasmon-pole model, static screening). These parameters directly affect 2D polarization and excitonic energies; their omission is load-bearing for the quantitative claims in the abstract and trilayer results.
  2. [Results] Results, trilayer section: The assertion of 'deterministic control' of electron localization relies on the internal-field picture and the 60-70 meV shifts, yet no explicit test is shown for how these values change under altered structural relaxations or different registries; this weakens the link between stacking sequence and the predicted electron routing.
minor comments (2)
  1. [Abstract] Abstract: The phrase 'remarkable agreement' with experiment for the 36 meV shift should be accompanied by a brief statement of the experimental reference value and uncertainty.
  2. [Results] Notation: Consistent use of 'quasiparticle band-edge shifts' versus 'internal electric fields' would clarify whether the 60-70 meV values are purely electrostatic or include chemical-potential contributions.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments and positive evaluation of the work's significance. We address each major comment point by point below. We agree that additional documentation of convergence and robustness tests will strengthen the manuscript and have revised accordingly.

read point-by-point responses

  1. Referee: [Methods] Methods/Computational Details: The reported 60-70 meV band-edge shifts and 36 meV interlayer excitonic shift are presented without convergence data on k-point sampling, vacuum padding, or sensitivity to GW approximations (plasmon-pole model, static screening). These parameters directly affect 2D polarization and excitonic energies; their omission is load-bearing for the quantitative claims in the abstract and trilayer results.

    Authors: We agree that explicit convergence data are important to support the quantitative claims. Although our calculations employed parameters consistent with prior GW+BSE studies on TMD heterostructures (dense k-grids, 20 Å vacuum, and the plasmon-pole approximation), we did not report dedicated convergence tests in the original submission. We have now performed additional calculations varying k-point sampling (up to 24×24×1), vacuum thickness (up to 30 Å), and GW screening models. The 60–70 meV shifts and 36 meV excitonic shift remain stable within ~5 meV. A new subsection and supplementary figure documenting these tests will be added to the revised Methods section. revision: yes

  2. Referee: [Results] Results, trilayer section: The assertion of 'deterministic control' of electron localization relies on the internal-field picture and the 60-70 meV shifts, yet no explicit test is shown for how these values change under altered structural relaxations or different registries; this weakens the link between stacking sequence and the predicted electron routing.

    Authors: The referee is correct that direct sensitivity tests to structural variations would reinforce the connection between stacking and electron routing. Our original calculations used fully relaxed geometries for each registry, with the reported shifts arising from the stacking-dependent polarization. To address the concern, we have tested small deviations in interlayer spacing (±0.1 Å) around the equilibrium values and confirmed that the band-edge shifts vary by less than 10 meV. We will add a short paragraph in the revised Results section (and a supplementary table) summarizing this robustness check, thereby strengthening the evidence for deterministic control via stacking sequence. revision: yes

Circularity Check

0 steps flagged

No circularity: GW+BSE results presented as independent predictions

full rationale

The paper reports quasiparticle band-edge shifts and a 36 meV interlayer excitonic shift obtained from GW+BSE calculations on chosen MoS2/MoSe2 stackings. These are framed as predictions that happen to agree with external experiments, with no indication that the 36 meV value or the 60-70 meV shifts were fitted to those experiments or derived by re-expressing input parameters. No self-citations, ansatzes smuggled via prior work, or self-definitional steps appear in the provided abstract or derivation description. The central claim therefore rests on the external validity of the many-body method rather than on any internal reduction to its own fitted inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on the accuracy of the GW+BSE method and the representativeness of the modeled atomic registries; no free parameters or invented entities are mentioned in the abstract.

axioms (1)
  • domain assumption GW+BSE many-body perturbation theory accurately describes stacking-dependent polarization and excitonic response in these heterostructures
    Invoked as the computational framework without stated limitations or validation steps in the abstract.

pith-pipeline@v0.9.1-grok · 5722 in / 1181 out tokens · 30454 ms · 2026-06-26T12:04:21.787943+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references

  1. [1]

    Zhao K, He D, Fu S, Bai Z, Miao Q, Huang M, Wang Y and Zhang X 2022 Nanomaterials123418

    2022

  2. [2]

    Xue F, He J H and Zhang X 2021Applied Physics Reviews8021316

  3. [3]

    Vizner Stern M, Waschitz Y, Cao W, Nevo I, Watanabe K, Taniguchi T, Sela E, Urbakh M, Hod O and Ben Shalom M 2021Science3721462–1466

  4. [4]

    Chen X, Xuan X, Guo W and Zhang Z 2025npj 2D Materials and Applications9 10

  5. [5]

    Niu Y, Li L, Qi Z, Aung H H, Han X, Tenne R, Yao Y, Zak A and Guo Y 2023 Nature Communications145578

    2023

  6. [7]

    Tran K, Choi J and Singh A 20212D Materials8022002

  7. [8]

    Amin B, Singh N and Schwingenschl¨ ogl U 2015Physical Review B92075439

  8. [9]

    Berkelbach T C, Hybertsen M S and Reichman D R 2013Physical Review B88 045318

  9. [10]

    Mak K F, Lee C, Hone J, Shan J and Heinz T F 2010Physical Review Letters105 136805

  10. [11]

    Wang F, Wang J, Guo S, Zhang J, Hu Z and Chu J 2017Scientific Reports744712

  11. [12]

    Xiao D, Liu G B, Feng W, Xu X and Yao W 2012Physical Review Letters108 196802

  12. [13]

    Wang G, Chernikov A, Glazov M M, Heinz T F, Marie X, Amand T and Urbaszek B 2018Reviews of Modern Physics90021001

  13. [14]

    Hong X, Kim J, Shi S F, Zhang Y, Jin C, Sun Y, Tongay S, Wu J, Zhang Y and Wang F 2014Nature Nanotechnology9682–686

  14. [15]

    Ceballos F, Bellus M Z, Chiu H Y and Zhao H 2014ACS Nano812717–12724

  15. [16]

    Tran K, Moody G, Wu F, Lu X, Choi J and Li X 2019Nature56771–75 REFERENCES21

  16. [17]

    Zhang N, Surrente A, Baranowski M, Maude D K, Gant P C, Castellanos-Gomez A and Plochocka P 2018Nano Letters187651–7657

  17. [18]

    Choi J, Florian M, Steinhoff A, Erben D, Tran K, Kim D S and Li X 2021Physical Review Letters126047401

  18. [19]

    Jiang Y, Chen S, Zheng W and Pan A 2021Light: Science & Applications1072

  19. [20]

    Li S, Wei K, Liu Q, Tang Y and Jiang T 2024Frontiers of Physics1942501

  20. [21]

    Jung E, Park J C, Seo Y S, Kim J H, Hwang J and Lee Y H 2022Scientific Reports 124543

  21. [22]

    Lin Q, Fang H, Kalaboukhov A, Liu Y, Zhang Y, Fischer M, Li J, Hagel J, Brem S, Malic E, Stenger N, Sun Z, Wubs M and Xiao S 2024Nature Communications 158762

  22. [23]

    Zhang Q, Fan A, Wang Y, Wu F, Li L, Meng H and Geng D 2025npj 2D Materials and Applications976

  23. [24]

    Cao W, Hod O and Urbakh M 2022ACS Applied Materials & Interfaces1457492– 57499

  24. [25]

    Pan J, Wang H, Zou L, Wang X, Zhang L, Dong X, Xie H B, Ding Y, Zhang Y, Watanabe K, Taniguchi T, Wang S and Wang Z 2025Physical Review Letters135 126603

  25. [26]

    Yang D, Liang J, Wu J and Ye Z 2024Nature Communications151389

  26. [27]

    Jiang Ying import C S, Zheng W, Zheng B and Pan A 2021Light: Science & Applications1072

  27. [28]

    Hou S, Xie X, Li S, Chen J, Zhang X and Liu Y 2025ACS Nano1934954–34965

  28. [29]

    Schwandt-Krause J, Miloudi M E A, Blundo E, Deb S, Heidkamp J N, Watanabe K, Taniguchi T and Korn T 2026Nano Letters26214–221

  29. [30]

    Bl¨ ochl P E 1994Physical Review B5017953–17979

  30. [31]

    Kresse G and Joubert D 1999Physical Review B591758–1775

  31. [32]

    Kresse G and Furthm¨ uller J 1996Physical Review B5411169–11186

  32. [33]

    Kresse G and Furthm¨ uller J 1996Computational Materials Science615–50

  33. [34]

    Perdew J P, Burke K and Ernzerhof M 1996Physical Review Letters773865–3868

  34. [35]

    Hybertsen M S and Louie S G 1986Physical Review B345390–5413

  35. [36]

    Wilhelm J, Seewald P and Golze D 2021Journal of Chemical Theory and Computation171662–1677

  36. [37]

    Hanke W and Sham L J 1980Physical Review B214656–4673

  37. [38]

    Strinati G 1988Rivista del Nuovo Cimento111–86

  38. [39]

    Thygesen K S 20172D Materials4022004

  39. [40]

    S1) • Bilayer electrostatic potentials (Fig

    Rivera P, Schaibley J R, Jones A M, Ross J S, Wu S, Aivazian G, Klement P, Seyler K, Clark G, Ghimire N J, Yan J, Mandrus D G and Xu X 2015Nature Communications66242 Supporting Material: Stacking-Directed Polarization and Excitonic Engineering in MoS 2/MoSe2 van der Waals Heterostructures Mohammed El Amine Miloudi and Oliver K¨ uhn Institute of Physics, U...