[论文解读] Twister: Construction and structural relaxation of commensurate moir\'e superlattices

Saismit Naik, Mit H. Naik+2 more|arXiv (Cornell University)|Feb 15, 2021

Graphene research and applications参考文献 76被引用 39

一句话总结

Twister 是一个基于 Python 的开源软件包，采用共格晶格法，可为任意二维异质结或同质结体系构建一致的莫尔超晶格，通过经典力场实现精确的结构弛豫。该工具自动化识别最小面积超晶格，并与 LAMMPS 接口，实现高效弛豫，相比 DFT 方法显著降低计算成本，同时保持高精度。

ABSTRACT

Introduction of a twist between layers of two-dimensional materials which leads to the formation of a moir\'e pattern is an emerging pathway to tune the electronic, vibrational and optical properties. The fascinating properties of these systems is often linked to large-scale structural reconstruction of the moir\'e pattern. Hence, an essential first step in the theoretical study of these systems is the construction and structural relaxation of the atoms in the moir\'e superlattice. We present the Twister package, a collection of tools that constructs commensurate superlattices for any combination of 2D materials and also helps perform structural relaxations of the moir\'e superlattice. Twister constructs commensurate moir\'e superlattices using the coincidence lattice method and provides an interface to perform structural relaxations using classical forcefields.

研究动机与目标

• 开发一种通用工具，用于构建任意二维异质结与同质结体系中的一致莫尔超晶格。
• 通过基于经典力场的高效弛豫方法，解决莫尔超晶格结构弛豫的计算瓶颈问题。
• 自动化识别在用户定义的应变容差和扭转角范围内的最小面积超晶格。
• 通过与 LAMMPS 集成，实现从超晶格构建到原子弛豫的全流程自动化，降低非专家用户的使用门槛。
• 提供一个开源、并行化且可扩展的框架，用于莫尔体系的理论研究。

提出的方法

• 采用共格晶格法，通过寻找晶格矢量的整数组合，实现上下两层的近似完美对齐，从而识别一致的超晶格。
• 通过形变张量或均匀缩放对各层施加应变，当在合理尺度下无法实现完美共格时，减小超晶格尺寸。
• 在晶格矢量索引（n1, n2, m1, m2）的四维空间中进行系统搜索，全面探索所有可能的超晶格组合。
• 利用 NumPy 和 MPI4Py 实现高性能、多处理器并行计算，高效生成候选超晶格。
• 基于用户定义的失配阈值、超晶格矢量夹角和应变限制，实施筛选以识别最优超晶格。
• 在超胞内生成原子坐标，并与 LAMMPS 接口，实现基于经典分子动力学或能量最小化法的结构弛豫。

实验结果

研究问题

• RQ1 对于任意二维异质结体系（包括晶格常数不匹配的情况），其最小面积的一致超晶格是什么？
• RQ2 如何在不依赖昂贵的 DFT 计算的前提下，利用经典力场高效模拟莫尔超晶格的结构弛豫？
• RQ3 在保持物理精度的前提下，对二维层施加应变以获得更小、计算上可处理的超晶格，其适用范围有多大？
• RQ4 经典力场弛豫结果在多大程度上能复现 DFT 模拟中观察到的莫尔超晶格结构重构？
• RQ5 能否构建一个统一的自动化框架，实现从超晶格构建到结构弛豫的全流程简化，适用于多种二维材料？

主要发现

• Twister 已成功为五个测试体系构建一致超晶格：扭转双层 MoS2、MoS2/MoSe2、石墨烯/hBN、hBN/MoSe2 以及扭转双层黑磷。
• 该软件包以高精度识别出 MoS2/MoSe2 和扭转双层黑磷的共格扭转角，包括 16.1° 和 4.265° 等角度。
• 使用经典力场对扭转双层 MoS2、石墨烯和 MoS2/MoSe2 进行结构弛豫，所得重构结果与 DFT 结果高度一致，但计算成本仅为后者的几分之一。
• 采用 MPI 并行化后，中等规模体系的超晶格搜索时间从数小时缩短至数分钟，且在 8–16 个核心上表现出良好的性能扩展性。
• 该软件包实现了精确的原子数验证与边界条件处理，确保超胞中无原子遗漏或重复计数。
• 与 LAMMPS 的集成实现了从超晶格生成到弛豫的无缝衔接，显著降低了非专家用户的使用门槛。