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

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

基本信息

论文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₂双分子层中堆叠顺序和旋转角度的联合效应。从五种高对称性双分子层堆叠开始，对顶层施加1°到120°的扭转角，揭示了各种弛豫摩尔结构。结果表明初始堆叠显著影响给定扭转角下形成的摩尔畴配置。虽然所有五种堆叠顺序在无扭转时都是亚稳态的，但它们形成两个摩尔等价类：AA/AB和AA'/A'B/AB'。具体地，初始AA和AB堆叠在0±3°附近产生三角形铁电畴，而AA'、A'B和AB'堆叠在60±3°附近产生三角形铁电畴。

研究背景与动机

研究问题

该研究要解决的核心问题是：堆叠顺序如何影响扭转二维材料双分子层中摩尔超晶格的结构和铁电性质。

重要性

扭转电子学的发展：摩尔超晶格展现出超导性、关联驱动的Mott绝缘态和量子反常霍尔效应等奇异性质
应用前景：在量子信息、光电子学和应变工程中具有重要应用潜力
理论完善：现有研究主要关注扭转角的影响，而堆叠顺序的作用相对未被充分探索

现有方法局限性

大多数研究集中在特定堆叠配置（如3R或2H）
缺乏对不同初始堆叠顺序系统性比较
对摩尔等价类的理解不够深入

研究动机

通过系统性的分子动力学模拟，建立堆叠顺序与摩尔畴形态之间的明确关系，为扭转电子学器件设计提供理论指导。

核心贡献

发现摩尔等价类：证明五种高对称堆叠形成两个摩尔等价类（3R类：AA/AB；2H类：AA'/A'B/AB'）
建立旋转关系：揭示60°旋转将3R堆叠转换为2H堆叠的旋转对称性
确定铁电区域：精确定位铁电三角畴的角度范围（3R类：0±3°；2H类：60±3°）
系统性分析方法：开发了基于分子动力学的堆叠域分类和偶极矩计算方法

方法详解

任务定义

输入：五种高对称性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