[論文レビュー] 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ベースのオープンソースパッケージであり、任意の2次元ヘテロバイレイヤーまたはホモバイレイヤー系に対して、一致格子法を用いて共通モアレ超格子を構築する。このツールは最小面積超格子の自動特定を実行し、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.

研究の動機と目的

• 任意の2次元ヘテロバイレイヤーおよびホモバイレイヤー系における共通モアレ超格子の構築を目的とした汎用ツールの開発。
• モアレ超格子における構造的弛緩の計算的ボトルネックを解消し、高効率な古典的力場に基づく弛緩を可能にすること。
• ユーザーが定義した歪み許容範囲およびねつぎ角度範囲に応じた最小面積超格子の自動特定。
• LAMMPSを用いた超格子構築から原子的弛緩までのワークフローを簡素化し、非専門家が利用しやすくすること。
• 理論的研究のためのオープンソースで並列処理可能かつ拡張可能なフレームワークの提供。

提案手法

• 一致格子法を用い、格子ベクトルの整数線形結合によって上層と下層の完全な整合を達成する条件を特定することで、共通モアレ超格子を同定する。
• 完全な一致が現実的なスケールで達成できない場合、歪みをデフォーメーションテンソルまたは均一スケーリングにより層に適用することで、超格子サイズを小さくする。
• 格子ベクトルのインデックス（n1, n2, m1, m2）の4次元探索を実施し、すべての可能な超格子組み合わせを体系的に調査する。
• 複数プロセッサ上で候補超格子の高性能並列計算を実現するため、NumPy および MPI4Py を活用する。
• ユーザーが定義した不一致閾値、超格子ベクトル間の角度、歪み制限に基づくフィルタリングを実装し、最適な超格子を特定する。
• 超格子細胞内での原子座標の生成と、古典的分子動力学またはエネルギー最小化に基づく構造的弛緩のためのLAMMPSとのインターフェースを実装する。

実験結果

リサーチクエスチョン

• RQ1 格子定数の不一致を含む任意の2次元ヘテロバイレイヤー系において、最小面積の共通モアレ超格子は何か？
• RQ2 高価なDFT計算に依存せずに、古典的力場を用いた効率的なモアレ超格子の構造的弛緩シミュレーションはどのように実現できるか？
• RQ3 物理的正確性を維持しつつ、計算的に取り扱いやすい小さな超格子を実現するために、2次元層にどの程度の歪みを適用できるか？
• RQ4 古典的力場による弛緩が、DFTシミュレーションで観測された構造的再配置をどの程度正確に再現できるか？
• RQ5 多様な2次元材料を対象とした、超格子構築から構造的弛緩までのワークフローを統合的かつ自動化したフレームワークを構築できるか？

主な発見

• Twister は、ねつぎビレイヤーMoS2、MoS2/MoSe2、グラフェン/hBN、hBN/MoSe2、ねつぎビレイヤー黒リンの5つのテスト系において、共通モアレ超格子を正常に構築した。
• MoS2/MoSe2 およびねつぎビレイヤー黒リンにおいて、高精度に共通ねつぎ角度を特定し、16.1° や 4.265° といった角度を含む。
• 古典的力場を用いたねつぎビレイヤーMoS2、グラフェン、MoS2/MoSe2 の構造的弛緩は、DFT結果と整合性がありながらも、計算コストを大幅に削減した。
• MPI並列処理の導入により、中程度のスケールの系において、超格子探索の実行時間を数時間から数分に短縮し、8〜16コアで良好なスケーリング性能を示した。
• 原子数の正確な検証と境界条件の適切な処理を実現し、超格子内での原子の漏れや重複を回避した。
• LAMMPSとの統合により、超格子生成から弛緩までのシームレスな移行が可能になり、非専門家が利用する障壁を顕著に低減した。