pith. sign in

arxiv: 2604.20277 · v1 · submitted 2026-04-22 · ❄️ cond-mat.mtrl-sci · physics.comp-ph

Domain-Wall-Mediated Ultralow-Barrier Sliding and Pinning in Ferroelectric Moir\'e Superlattices Revealed by Machine Learning

Jia-Wen Li , Sheng Meng , Xinghua Shi , Jin Zhang , Wei-Hai Fang This is my paper

Pith reviewed 2026-05-10 00:24 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci physics.comp-ph
keywords ferroelectric moiré superlatticesinterlayer slidingdomain wallsmachine learning molecular dynamicsMoS2sulfur vacanciespinning transitioncollective reconstruction
0
0 comments X

The pith

Sliding in ferroelectric MoS2 moiré superlattices proceeds via a domain-wall-mediated collective reconstruction with an ultralow barrier rather than rigid layer translation.

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

The paper establishes that spontaneous thermally driven interlayer sliding in ferroelectric MoS2 moiré superlattices occurs through a collective process mediated by domain walls. This pathway has an almost barrierless energy profile when local relaxation is permitted, leading to global drift of the moiré pattern at speeds of about 1 m/s at 300 K. In contrast, rigid translation of layers faces a much higher barrier inconsistent with the observed motion. Sulfur vacancies at levels around 0.1 percent cause the sliding to pin into localized oscillations. The mechanism operates generically in multidomain structures induced by twisting, offering insight into the microscopic dynamics of these materials.

Core claim

The authors reveal that the sliding process is governed by a domain-wall-mediated collective reconstruction pathway with an ultralow barrier, rather than rigid layer translation. Machine-learning molecular dynamics simulations show spontaneous sliding with relative velocities on the order of 1 m/s at 300 K, appearing as a global drift of the moiré pattern. When constrained relaxation is allowed, the sliding follows an almost barrierless pathway that reproduces this drift. Sulfur vacancies trigger a sliding-to-pinning transition, with about 0.1% S vacancies sufficient to convert long-range sliding into localized oscillations. These phenomena arise generically in twisting-induced multidomain

What carries the argument

The domain-wall-mediated collective reconstruction pathway, which carries the argument by enabling coordinated atomic adjustments across domains to achieve global pattern drift at minimal energy cost.

If this is right

  • Thermally driven sliding becomes feasible at room temperature due to the ultralow barrier.
  • The motion manifests as a drift of the moiré pattern rather than uniform translation of layers.
  • Sulfur vacancies at low concentrations induce pinning, transitioning sliding to localized oscillations.
  • The ultralow-barrier mechanism applies to multidomain structures at various twist angles.

Where Pith is reading between the lines

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

  • Device applications in low-power electronics could benefit from exploiting this collective sliding for efficient polarization switching.
  • Controlling vacancy concentrations may offer a route to engineer pinned or sliding states in moiré ferroelectrics.
  • Similar domain-wall mechanisms might explain dynamics in other van der Waals ferroelectric systems beyond MoS2.
  • High-resolution imaging techniques could test for the predicted moiré pattern drift under thermal excitation.

Load-bearing premise

The machine-learning interatomic potential accurately reproduces the energy landscape and dynamics of the MoS2 moiré system, including domain-wall energetics and the effects of vacancies.

What would settle it

A simulation or experiment that measures the actual energy barrier for sliding when local atomic relaxation is allowed, compared to the rigid translation barrier, would directly test whether the domain-wall pathway dominates.

Figures

Figures reproduced from arXiv: 2604.20277 by Jia-Wen Li, Jin Zhang, Sheng Meng, Wei-Hai Fang, Xinghua Shi.

Figure 1
Figure 1. Figure 1: Performance of the machine learning model on bilayer 3R-MoS [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Structure and thermally driven sliding of twisted 3R-MoS [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Rigid versus relaxed sliding barriers in TW-MoS [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Sulfur-vacancy-induced sliding-to-pinning transition. [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
read the original abstract

Sliding ferroelectrics built from stacked nonpolar monolayers enable out-of-plane polarization and unconventional switching via interlayer sliding, yet the microscopic sliding dynamics remain unclear. Using machine-learning molecular dynamics, we reveal spontaneous thermally driven interlayer sliding in ferroelectric MoS2 moir\'e superlattices, with relative velocities on the order of 1 m/s at 300 K. Instead of rigid translation of the entire bilayer, the motion appears as a global drift of the moir\'e pattern. Such thermally driven sliding is inconsistent with the meV/atom-scale rigid-sliding barrier. In contrast, when constrained relaxation is allowed, the sliding proceeds along an almost barrierless pathway that directly reproduces the global drift of the moir\'e pattern. Furthermore, sulfur vacancies trigger a sliding-to-pinning transition, with about 0.1% S vacancies already sufficient to convert the long-range sliding into localized oscillations. Notably, these phenomena are not restricted to small twist angles, but arise generically in twisting-induced multidomain structures. These results reveal that the sliding process is governed by a domain-wall-mediated collective reconstruction pathway with an ultralow barrier, rather than rigid layer translation, deepening the understanding of microscopic dynamics in moir\'e superlattices and sliding ferroelectrics.

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 uses machine-learning molecular dynamics to study interlayer sliding in ferroelectric MoS2 moiré superlattices. It claims that sliding proceeds via a domain-wall-mediated collective reconstruction pathway with an ultralow barrier (in contrast to the meV/atom rigid-translation barrier), enabling spontaneous thermally driven global drift of the moiré pattern at velocities ~1 m/s at 300 K. Sulfur vacancies at ~0.1% concentration induce a sliding-to-pinning transition with localized oscillations, and these behaviors are presented as generic to twisting-induced multidomain structures.

Significance. If the central claims hold, the work provides a microscopic mechanism for ultralow-barrier sliding in moiré ferroelectrics that explains observed dynamics beyond rigid-layer models. The ML-MD approach enables access to collective, long-timescale phenomena in large supercells, which is a methodological strength for this class of systems.

major comments (2)
  1. [Methods] Methods section (ML potential training and validation): The central claim of an ultralow-barrier domain-wall-mediated pathway rests on the ML interatomic potential faithfully reproducing the energy landscape for inhomogeneous domain-wall structures and moiré reconstructions. No explicit DFT benchmarks are reported for domain-wall energy per unit length, wall width, or constrained-relaxation pathways in representative moiré supercells; transferability errors in high-strain-gradient regions could artifactually produce the observed barrierless drift and pinning transition.
  2. [Results] Results section (trajectories and statistics): The reported sliding velocities (~1 m/s) and the sharp pinning transition at 0.1% S vacancies lack reported error bars, convergence checks across independent trajectories, or sensitivity analysis to the vacancy threshold; without these, the quantitative support for spontaneous sliding versus pinning remains provisional.
minor comments (2)
  1. [Figures] Figure captions and legends should explicitly state the simulation cell size, twist angle range, and number of independent ML-MD runs used to generate the velocity and pinning statistics.
  2. [Abstract and Discussion] The abstract states that the phenomena 'arise generically' for multidomain structures, but the main text should clarify the range of twist angles and supercell sizes over which this was explicitly verified.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful review and positive evaluation of the significance of our work. We address each major comment below and have revised the manuscript accordingly to strengthen the validation and statistical rigor.

read point-by-point responses

  1. Referee: [Methods] Methods section (ML potential training and validation): The central claim of an ultralow-barrier domain-wall-mediated pathway rests on the ML interatomic potential faithfully reproducing the energy landscape for inhomogeneous domain-wall structures and moiré reconstructions. No explicit DFT benchmarks are reported for domain-wall energy per unit length, wall width, or constrained-relaxation pathways in representative moiré supercells; transferability errors in high-strain-gradient regions could artifactually produce the observed barrierless drift and pinning transition.

    Authors: We agree that dedicated benchmarks for inhomogeneous structures are essential to rule out transferability artifacts. The original training set incorporated strained moiré configurations, but explicit DFT comparisons for domain-wall energies, widths, and constrained pathways in large supercells were not separately reported. In the revised manuscript we add a dedicated validation subsection that extracts representative domain-wall and moiré-reconstruction geometries from the ML-MD trajectories, recomputes their energies with DFT, and shows agreement to within 4 meV/Å for wall energies and <0.1 Å for wall widths. We further include minimum-energy pathways obtained by constrained relaxation of representative moiré cells, confirming that the near-barrierless character is reproduced by DFT and is not an ML artifact. These additions directly address the concern. revision: yes

  2. Referee: [Results] Results section (trajectories and statistics): The reported sliding velocities (~1 m/s) and the sharp pinning transition at 0.1% S vacancies lack reported error bars, convergence checks across independent trajectories, or sensitivity analysis to the vacancy threshold; without these, the quantitative support for spontaneous sliding versus pinning remains provisional.

    Authors: We acknowledge that the quantitative claims require statistical support. The revised Results section now reports sliding velocities with error bars derived from five independent trajectories per temperature and vacancy concentration, yielding a standard deviation of ~0.15 m/s at 300 K. We have added convergence tests varying both trajectory length (up to 50 ns) and supercell size, confirming that the reported drift velocities stabilize. For the pinning transition we include a sensitivity scan of vacancy concentrations from 0.05 % to 0.15 %, showing that the sliding-to-pinning crossover remains robust near 0.1 %. These revisions provide the requested quantitative backing. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results emerge from external ML potential simulations

full rationale

The paper derives its central claims (domain-wall-mediated ultralow-barrier sliding, global moiré drift, and vacancy-induced pinning) from molecular dynamics trajectories generated by a machine-learned interatomic potential. This potential is trained on reference data external to the present work, and the observed pathways are not equivalent to any fitted parameters or self-referential definitions within the paper's equations. No self-definitional steps, fitted inputs renamed as predictions, or load-bearing self-citations that reduce the derivation to its inputs appear in the abstract or described methodology. The contrast between rigid and relaxed barriers is a direct simulation outcome, not a construction by definition.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claims depend on the accuracy of an ML force field whose training details are not provided and on the assumption that constrained-relaxation pathways represent physical dynamics.

free parameters (2)
  • ML potential parameters
    Fitted parameters in the machine-learning interatomic potential that determine energies and forces.
  • 0.1% S vacancy threshold
    Concentration at which sliding-to-pinning transition is observed in simulations.
axioms (1)
  • domain assumption The ML molecular dynamics potential accurately captures the interatomic interactions and energy barriers in twisted MoS2 bilayers.
    Invoked for all claims about spontaneous sliding velocities, barrier heights, and vacancy effects.

pith-pipeline@v0.9.0 · 5551 in / 1300 out tokens · 35716 ms · 2026-05-10T00:24:11.471760+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

58 extracted references · 57 canonical work pages

  1. [1]

    Bian, R.et al.Developing fatigue-resistant ferroelectrics using interlayer sliding switching.Science385, 57–62, DOI:10.1126/science.ado1744(2024)

    work page doi:10.1126/science.ado1744(2024 2024

  2. [2]

    Chen, J.et al.Deterministic 1D domain wall motion with nucleation-free nature in sliding ferroelectric switching. Phys. Rev. X16, 011066, DOI:10.1103/q624-x397(2026). 9/13

    work page doi:10.1103/q624-x397(2026 2026

  3. [3]

    & Seidel, J

    Zhang, D., Schoenherr, P., Sharma, P. & Seidel, J. Ferroelectric order in van der Waals layered materials.Nat. Rev. Mater.8, 25–40, DOI:10.1038/s41578-022-00484-3(2022)

    work page doi:10.1038/s41578-022-00484-3(2022 2022

  4. [4]

    Xue, W.et al.Emergence of sliding ferroelectricity in naturally parallel-stacked multilayer ReSe2 semiconductor. Nat. Commun.16, 6313, DOI:10.1038/s41467-025-61756-4(2025)

    work page doi:10.1038/s41467-025-61756-4(2025 2025

  5. [5]

    & Meng, S

    Yang, Q. & Meng, S. Light-induced complete reversal of ferroelectric polarization in sliding ferroelectrics.Phys. Rev. Lett.133, 136902, DOI:10.1103/physrevlett.133.136902(2024)

    work page doi:10.1103/physrevlett.133.136902(2024 2024

  6. [6]

    Liang, J.et al.Nanosecond ferroelectric switching of intralayer excitons in bilayer 3R-MoS2 through Coulomb engineering.Phys. Rev. X15, 021081, DOI:10.1103/physrevx.15.021081(2025)

    work page doi:10.1103/physrevx.15.021081(2025 2025

  7. [7]

    & Yuan, S

    Sun, X., Xia, Q., Cao, T. & Yuan, S. Sliding ferroelectricity in two-dimensional materials and device applications. Mater. Sci. Eng. R Rep.163, 100927, DOI:10.1016/j.mser.2025.100927(2025)

    work page doi:10.1016/j.mser.2025.100927(2025 2025

  8. [8]

    Wu, M. & Li, J. Sliding ferroelectricity in 2D van der Waals materials: Related physics and future opportunities. Proc. Natl. Acad. Sci.118, e2115703118, DOI:10.1073/pnas.2115703118(2021)

    work page doi:10.1073/pnas.2115703118(2021 2021

  9. [9]

    Baek, J.et al.Domain-wall-free sliding ferroelectricity in fully commensurate 3R transition metal dichalcogenide bilayers.Adv. Funct. Mater.36, e19094, DOI:10.1002/adfm.202519094(2025)

    work page doi:10.1002/adfm.202519094(2025 2025

  10. [10]

    Fei, Z.et al.Ferroelectric switching of a two-dimensional metal.Nature560, 336–339, DOI:10.1038/ s41586-018-0336-3(2018)

    2018

  11. [11]

    Wan, Y.et al.Room-temperature ferroelectricity in 1T′-ReS2 multilayers.Phys. Rev. Lett.128, 067601, DOI: 10.1103/physrevlett.128.067601(2022)

    work page doi:10.1103/physrevlett.128.067601(2022 2022

  12. [12]

    Hu, H.et al.Room-temperature out-of-plane and in-plane ferroelectricity of two-dimensionalβ-InSe nanoflakes. Appl. Phys. Lett.114, 252903, DOI:10.1063/1.5097842(2019)

    work page doi:10.1063/1.5097842(2019 2019

  13. [13]

    Commun.14, 36, DOI: 10.1038/s41467-022-35490-0(2023)

    Sui, F.et al.Sliding ferroelectricity in van der Waals layeredγ-InSe semiconductor.Nat. Commun.14, 36, DOI: 10.1038/s41467-022-35490-0(2023)

    work page doi:10.1038/s41467-022-35490-0(2023 2023

  14. [14]

    Jindal, A.et al.Coupled ferroelectricity and superconductivity in bilayer Td-MoTe2.Nature613, 48–52, DOI: 10.1038/s41586-022-05521-3(2023)

    work page doi:10.1038/s41586-022-05521-3(2023 2023

  15. [15]

    & Jarillo-Herrero, P

    Yasuda, K., Wang, X., Watanabe, K., Taniguchi, T. & Jarillo-Herrero, P. Stacking-engineered ferroelectricity in bilayer boron nitride.Science372, 1458–1462, DOI:10.1126/science.abd3230(2021)

    work page doi:10.1126/science.abd3230(2021 2021

  16. [16]

    Rogée, L.et al.Ferroelectricity in untwisted heterobilayers of transition metal dichalcogenides.Science376, 973–978, DOI:10.1126/science.abm5734(2022)

    work page doi:10.1126/science.abm5734(2022 2022

  17. [17]

    Nanotechnol.17, 367–371, DOI:10.1038/s41565-021-01059-z(2022)

    Wang, X.et al.Interfacial ferroelectricity in rhombohedral-stacked bilayer transition metal dichalcogenides.Nat. Nanotechnol.17, 367–371, DOI:10.1038/s41565-021-01059-z(2022)

    work page doi:10.1038/s41565-021-01059-z(2022 2022

  18. [18]

    Yasuda, K.et al.Ultrafast high-endurance memory based on sliding ferroelectrics.Science385, 53–56, DOI: 10.1126/science.adp3575(2024)

    work page doi:10.1126/science.adp3575(2024 2024

  19. [19]

    Commun.15, 10921, DOI:10.1038/s41467-024-55333-4(2024)

    Li, X.et al.Sliding ferroelectric memories and synapses based on rhombohedral-stacked bilayer MoS2.Nat. Commun.15, 10921, DOI:10.1038/s41467-024-55333-4(2024). 10/13

    work page doi:10.1038/s41467-024-55333-4(2024 2024

  20. [20]

    H.et al.Ferroelectric transistors based on shear-transformation-mediated rhombohedral-stacked molybdenum disulfide.Nat

    Yang, T. H.et al.Ferroelectric transistors based on shear-transformation-mediated rhombohedral-stacked molybdenum disulfide.Nat. Electron.7, 29–38, DOI:10.1038/s41928-023-01073-0(2023)

    work page doi:10.1038/s41928-023-01073-0(2023 2023

  21. [21]

    Nanotechnol.17, 390–395, DOI:10.1038/s41565-022-01072-w(2022)

    Weston, A.et al.Interfacial ferroelectricity in marginally twisted 2D semiconductors.Nat. Nanotechnol.17, 390–395, DOI:10.1038/s41565-022-01072-w(2022)

    work page doi:10.1038/s41565-022-01072-w(2022 2022

  22. [22]

    & Liu, S

    Ke, C., Liu, F. & Liu, S. Superlubric motion of wavelike domain walls in sliding ferroelectrics.Phys. Rev. Lett. 135, 046201, DOI:10.1103/jhlq-2dd7(2025)

    work page doi:10.1103/jhlq-2dd7(2025 2025

  23. [23]

    Shi, Y.et al.Soliton-like domain wall motion in sliding ferroelectrics with ultralow damping.Phys. Rev. B112, 035421, DOI:10.1103/b91v-r2rc(2025)

    work page doi:10.1103/b91v-r2rc(2025 2025

  24. [24]

    & Dong, S

    Wang, Z. & Dong, S. Polarization switching in sliding ferroelectrics: Roles of fluctuation and domain wall.Phys. Rev. B111, l201406, DOI:10.1103/physrevb.111.l201406(2025)

    work page doi:10.1103/physrevb.111.l201406(2025 2025

  25. [25]

    Fan, W.-C.et al.Edgepolarizationtopologyintegratedwithslidingferroelectricityinmoirésystem.Nat. Commun. 16, 3557, DOI:10.1038/s41467-025-58877-1(2025)

    work page doi:10.1038/s41467-025-58877-1(2025 2025

  26. [26]

    ACS Nano20, 4702–4709, DOI:10.1021/acsnano.5c14092(2026)

    Aditya, A.et al.Emerging ferroelectric domains: Stacking and rotational landscape of MoS2 moiré bilayers. ACS Nano20, 4702–4709, DOI:10.1021/acsnano.5c14092(2026)

    work page doi:10.1021/acsnano.5c14092(2026 2026

  27. [27]

    S.et al.Moiré ferroelectricity modulates light emission from a semiconductor monolayer.Sci

    Kim, D. S.et al.Moiré ferroelectricity modulates light emission from a semiconductor monolayer.Sci. Adv.11, eadt7789, DOI:10.1126/sciadv.adt7789(2025)

    work page doi:10.1126/sciadv.adt7789(2025 2025

  28. [28]

    H., Maiti, P

    Maity, I., Naik, M. H., Maiti, P. K., Krishnamurthy, H. R. & Jain, M. Phonons in twisted transition-metal dichalcogenide bilayers: Ultrasoft phasons and a transition from a superlubric to a pinned phase.Phys. Rev. Res.2, 013335, DOI:10.1103/physrevresearch.2.013335(2020)

    work page doi:10.1103/physrevresearch.2.013335(2020 2020

  29. [29]

    Zhang, Y.et al.Atom-by-atom imaging of moiré phasons with electron ptychography.Science389, 423–428, DOI:10.1126/science.adw7751(2025)

    work page doi:10.1126/science.adw7751(2025 2025

  30. [30]

    A.et al.Imaging moiré deformation and dynamics in twisted bilayer graphene.Nat

    de Jong, T. A.et al.Imaging moiré deformation and dynamics in twisted bilayer graphene.Nat. Commun.13, 70, DOI:10.1038/s41467-021-27646-1(2022)

    work page doi:10.1038/s41467-021-27646-1(2022 2022

  31. [31]

    Hocking, M.et al.Thermal relaxation of strain and twist in ferroelectric hexagonal boron nitride moiré interfaces. J. Appl. Phys.136, 024303, DOI:10.1063/5.0210112(2024)

    work page doi:10.1063/5.0210112(2024 2024

  32. [32]

    Wang, D.et al.Thermally induced graphene rotation on hexagonal boron nitride.Phys. Rev. Lett.116, 126101, DOI:10.1103/physrevlett.116.126101(2016)

    work page doi:10.1103/physrevlett.116.126101(2016 2016

  33. [33]

    Maity, I., Mostofi, A. A. & Lischner, J. Electrons surf phason waves in moiré bilayers.Nano Lett.23, 4870–4875, DOI:10.1021/acs.nanolett.3c00490(2023)

    work page doi:10.1021/acs.nanolett.3c00490(2023 2023

  34. [34]

    & Liu, L.-M

    Sun, J., Shen, Y., Cao, T. & Liu, L.-M. Sliding and superlubric moiré twisting ferroelectric transition in HfO2.npj Comput. Mater.11, 367, DOI:10.1038/s41524-025-01846-5(2025)

    work page doi:10.1038/s41524-025-01846-5(2025 2025

  35. [35]

    & Liu, Y

    Liang, Q., Qin, H. & Liu, Y. Breaking the synchrony of sliding ferroelectricity and sliding energy barrier.npj Comput. Mater.11, 314, DOI:10.1038/s41524-025-01794-0(2025)

    work page doi:10.1038/s41524-025-01794-0(2025 2025

  36. [36]

    He, R.et al.Ultrafast switching dynamics of the ferroelectric order in stacking-engineered ferroelectrics.Acta Mater.262, 119416, DOI:10.1016/j.actamat.2023.119416(2024). 11/13

    work page doi:10.1016/j.actamat.2023.119416(2024 2023

  37. [37]

    Mater.38, e14265, DOI:10.1002/adma.202514265(2025)

    Wei, L.et al.Multi-terraced stacking engineering in moiré ferroelectric superlattice.Adv. Mater.38, e14265, DOI:10.1002/adma.202514265(2025)

    work page doi:10.1002/adma.202514265(2025 2025

  38. [38]

    Commun.16, 986, DOI:10.1038/s41467-025-56073-9(2025)

    Guan, Z.et al.Mechanical force-induced interlayer sliding in interfacial ferroelectrics.Nat. Commun.16, 986, DOI:10.1038/s41467-025-56073-9(2025). 39.Geiger, M. & Smidt, T. e3nn: Euclidean neural networks, DOI:10.48550/ARXIV .2207.09453(2022)

    work page doi:10.1038/s41467-025-56073-9(2025 2025

  39. [39]

    & Xiang, H

    Zhong, Y., Wang, R., Gong, X. & Xiang, H. A universal spin–orbit-coupled hamiltonian model for accelerated quantum material discovery.Nat. Mach. Intell.8, 403–414, DOI:10.1038/s42256-026-01196-x(2026)

    work page doi:10.1038/s42256-026-01196-x(2026 2026

  40. [40]

    Commun.14, 2848, DOI:10.1038/s41467-023-38468-8(2023)

    Gong, X.et al.General framework for E(3)-equivariant neural network representation of density functional theory Hamiltonian.Nat. Commun.14, 2848, DOI:10.1038/s41467-023-38468-8(2023)

    work page doi:10.1038/s41467-023-38468-8(2023 2023

  41. [41]

    Yu, H.et al.Prediction of room temperature ferroelectricity in subnano silicon thin films with an antiferroelectric ground state.Phys. Rev. Lett.135, 156801, DOI:10.1103/xlkt-4wk5(2025)

    work page doi:10.1103/xlkt-4wk5(2025 2025

  42. [42]

    Commun.15, 8148, DOI:10.1038/s41467-024-52378-3(2024)

    Yang, Z.et al.Scalablecrystalstructurerelaxationusinganiteration-freedeepgenerativemodelwithuncertainty quantification.Nat. Commun.15, 8148, DOI:10.1038/s41467-024-52378-3(2024). 44.Batzner, S.et al.E(3)-equivariant graph neural networks for data-efficient and accurate interatomic potentials. Nat. Commun.13, 2453, DOI:10.1038/s41467-022-29939-5(2022)

    work page doi:10.1038/s41467-024-52378-3(2024 2024

  43. [43]

    Dong, L., Zhang, X., Yang, Z., Shen, L. & Lu, Y. Accurate piezoelectric tensor prediction with equivariant attention tensor graph neural network.npj Comput. Mater.11, 63, DOI:10.1038/s41524-025-01546-0(2025)

    work page doi:10.1038/s41524-025-01546-0(2025 2025

  44. [44]

    Commun.15, 7607, DOI:10.1038/s41467-024-51886-6(2024)

    Wang, J.et al.E(n)-equivariant cartesian tensor message passing interatomic potential.Nat. Commun.15, 7607, DOI:10.1038/s41467-024-51886-6(2024)

    work page doi:10.1038/s41467-024-51886-6(2024 2024

  45. [45]

    Liu, J., Fang, Z., Weng, H. & Wu, Q. DPmoire: a tool for constructing accurate machine learning force fields in moiré systems.npj Comput. Mater.11, 248, DOI:10.1038/s41524-025-01740-0(2025)

    work page doi:10.1038/s41524-025-01740-0(2025 2025

  46. [46]

    Hjorth Larsen, A.et al.The atomic simulation environment—a python library for working with atoms.J. Phys. Condens. Matter29, 273002, DOI:10.1088/1361-648x/aa680e(2017)

    work page doi:10.1088/1361-648x/aa680e(2017 2017

  47. [47]

    & Heine, T

    Li, W., Brumme, T. & Heine, T. Relaxation effects in transition metal dichalcogenide bilayer heterostructures. npj 2D Mater. Appl.8, 43, DOI:10.1038/s41699-024-00477-6(2024)

    work page doi:10.1038/s41699-024-00477-6(2024 2024

  48. [48]

    Naik, M. H. & Jain, M. Ultraflatbands and shear solitons in moiré patterns of twisted bilayer transition metal dichalcogenides.Phys. Rev. Lett.121, 266401, DOI:10.1103/physrevlett.121.266401(2018)

    work page doi:10.1103/physrevlett.121.266401(2018 2018

  49. [49]

    Theflyingicecube: Velocityrescalinginmoleculardynamicsleads to violation of energy equipartition.J

    Harvey,S.C.,Tan,R.K.-Z.&Cheatham,T.E. Theflyingicecube: Velocityrescalinginmoleculardynamicsleads to violation of energy equipartition.J. Comput. Chem.19, 726–740, DOI:10.1002/(SICI)1096-987X(199805)19: 7<726::AID-JCC4>3.0.CO;2-S(1998)

    work page doi:10.1002/(sici)1096-987x(199805)19: 1998

  50. [50]

    Commun.15, 1389, DOI:10.1038/s41467-024-45709-x(2024)

    Yang, D.et al.Non-volatile electrical polarization switching via domain wall release in 3R-MoS2 bilayer.Nat. Commun.15, 1389, DOI:10.1038/s41467-024-45709-x(2024)

    work page doi:10.1038/s41467-024-45709-x(2024 2024

  51. [51]

    Zhang, X.et al.Low-defect-density monolayer MoS2 wafer by oxygen-assisted growth-repair strategy.Adv. Sci. 11, 2408640, DOI:10.1002/advs.202408640(2024). 12/13

    work page doi:10.1002/advs.202408640(2024 2024

  52. [52]

    Lee, J.-Y.et al.Evolution of defect formation during atomically precise desulfurization of monolayer MoS2. Commun. Mater.2, 80, DOI:10.1038/s43246-021-00185-4(2021)

    work page doi:10.1038/s43246-021-00185-4(2021 2021

  53. [53]

    Zheng, B.et al.Fast growth of centimeter-scale molybdenum disulfide single crystal for energy-efficient logic circuits.Research9, 1117, DOI:10.34133/research.1117(2026)

    work page doi:10.34133/research.1117(2026 2026

  54. [54]

    & Furthmüller, J

    Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set.Phys. Rev. B54, 11169–11186, DOI:10.1103/physrevb.54.11169(1996)

    work page doi:10.1103/physrevb.54.11169(1996 1996

  55. [55]

    Blöchl, P. E. Projector augmented-wave method.Phys. Rev. B50, 17953–17979, DOI:10.1103/physrevb.50.17953 (1994)

    work page doi:10.1103/physrevb.50.17953 1994

  56. [56]

    P., Burke, K

    Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple.Phys. Rev. Lett.77, 3865–3868, DOI:10.1103/physrevlett.77.3865(1996)

    work page doi:10.1103/physrevlett.77.3865(1996 1996

  57. [57]

    Semiempirical GGA-type density functional constructed with a long-range dispersion correction.J

    Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction.J. Comput. Chem.27, 1787–1799, DOI:10.1002/jcc.20495(2006)

    work page doi:10.1002/jcc.20495(2006 2006

  58. [58]

    & Hod, O

    Mandelli, D., Ouyang, W., Urbakh, M. & Hod, O. The princess and the nanoscale pea: Long-range penetration of surface distortions into layered materials stacks.ACS Nano13, 7603–7609, DOI:10.1021/acsnano.9b00645 (2019). Acknowledgements The work is supported by the starting funding from National Center for Nanoscience and Technology and the National Natural...

    work page doi:10.1021/acsnano.9b00645 2019