2025-11-19T09:04:14.166965

Emerging Ferroelectric Domains: Stacking and Rotational Landscape of MoS2 Moire Bilayers

Aditya, Irie, Dasgupta et al.

The structures and properties of moire patterns in twisted bilayers of two-dimensional (2D) materials are known to depend sensitively on twist angle, yet their dependence on stacking order remains comparatively underexplored. In this study, we use molecular dynamics simulations to systematically investigate the combined effects of stacking order and rotation in MoS2 bilayers. Beginning from five well-established high-symmetry bilayer stackings, we apply twist angles between 1 and 120 to the top layer, revealing a variety of relaxed moire structures. Our results show that the initial stacking significantly influences the moire domain configurations that emerge at a given twist angle. While all five stacking orders are metastable without twist, they form two moire-equivalent classes- AA/AB and AA',A'B,AB', i.e., for a given twist angle, structures within each class relax to the same moire configuration. Specifically, initial AA and AB stackings give rise to triangular ferroelectric domains near 0+/-3, while AA', A'B, and AB' stackings produce triangular ferroelectric domains near 60+/-3. At precisely 60 and 120 twists, the bilayers relax to into pure high-symmetry stackings, highlighting the rotational relationships between these configurations and explaining the shift of 60 in the ferroelectric rotational range. These findings demonstrate the critical role of stacking order in governing the rich moire landscapes accessible in twistronic systems.

academic

新興フェロ電気ドメイン：MoS₂モアレ二層膜の積層とロテーション景観

基本情報

論文ID: 2510.10831
タイトル: Emerging Ferroelectric Domains: Stacking and Rotational Landscape of MoS2 Moire Bilayers
著者: Anikeya Aditya、Ayu Irie、Nabankur Dasgupta、Rajiv K. Kalia、Aiichiro Nakano、Priya Vashishta
分類: cond-mat.mtrl-sci（凝縮系物理学-材料科学）
機関: 南カリフォルニア大学先進計算・シミュレーション協働研究所、熊本大学物理学科
論文リンク: https://arxiv.org/abs/2510.10831

概要

本研究は分子動力学シミュレーションを用いて、MoS₂二層膜における積層順序と回転角度の複合効果を系統的に調査した。5つの高対称性二層膜積層から開始し、上層に1°から120°の捻り角を適用して、様々な緩和モアレ構造を明らかにした。結果は、初期積層が与えられた捻り角で形成されるモアレドメイン配置に大きく影響することを示した。5つの積層順序すべてが無捻りの状態では準安定的であるが、2つのモアレ等価類を形成する：AA/ABおよびAA'/A'B/AB'。具体的には、初期AA積層とAB積層は0±3°付近で三角形フェロ電気ドメインを生成し、AA'、A'B、およびAB'積層は60±3°付近で三角形フェロ電気ドメインを生成する。

研究背景と動機

研究課題

本研究が解決しようとする中核的な問題は：積層順序が捻られた二次元材料二層膜のモアレ超格子の構造とフェロ電気特性にどのように影響するかである。

重要性

捻り電子工学の発展：モアレ超格子は超伝導性、相関駆動型モット絶縁体状態、量子異常ホール効果などの異常な性質を示す
応用前景：量子情報、光電子工学、ひずみ工学における重要な応用可能性
理論の完善：既存研究は主に捻り角の影響に焦点を当てており、積層順序の役割は相対的に十分に探索されていない

既存方法の限界

ほとんどの研究は特定の積層配置（3Rまたは2Hなど）に集中している
異なる初期積層順序の系統的比較が不足している
モアレ等価類の理解が不十分である

研究動機

分子動力学シミュレーションの系統的実施により、積層順序とモアレドメイン形態の間に明確な関係を確立し、捻り電子工学デバイス設計に理論的指針を提供する。

核心的貢献

モアレ等価類の発見：5つの高対称積層が2つのモアレ等価類を形成することを証明（3R類：AA/AB；2H類：AA'/A'B/AB'）
回転関係の確立：60°回転が3R積層を2H積層に変換する回転対称性を明らかにした
フェロ電気領域の特定：フェロ電気三角ドメインの角度範囲を正確に特定（3R類：0±3°；2H類：60±3°）
系統的分析方法：分子動力学に基づく積層ドメイン分類と双極子モーメント計算方法を開発した

方法の詳細

タスク定義

入力：5つの高対称性MoS₂二層膜積層（AA、AB、AA'、A'B、AB'）出力：異なる捻り角での緩和モアレ構造と積層ドメイン分布 制約条件：捻り角範囲1°-120°、ステップサイズ1°

モデルアーキテクチャ

1. 分子動力学シミュレーションフレームワーク

力場の選択：
- Stillinger-Weber (SW) 力場：層内Mo-S相互作用
- Kolmogorov-Crespi (KC) 力場：層間相互作用
幾何学的設定：直径4000 Åの円形二層膜薄片
境界条件：エッジ50 Å領域の原子を固定してスリップを防止

2. 積層ドメイン分類アルゴリズム

# 積層認識アルゴリズムの中核的考え方
for each Mo atom in top layer:
    identify 3 closest S neighbors
    create cylindrical region around Mo + 3S atoms
    assign signature based on atoms below in cylinder
    classify stacking type by unique signature