Twistronic control of shift current in multilayer moiré system

技術要約：マルチレイヤー・モアレ系におけるシフト電流のツイストトロニクス制御

問題提起
非中心対称材料における光から電流への変換（バルク光起電力効果：BPVE）は、従来のp-n接合に代わる手法を提供し、ショックレー・クワイサー効率限界を克服できる可能性がある。BPVEのメカニズムの中でも、シフト電流（SC）は電子波関数の幾何学的特性に依存することから、特に重要である。遷移金属ダイカルコゲナイド（TMD）のような二次元（2D）材料は、このような用途への有望な候補であるが、ねじれたマルチレイヤー構成における大規模なモアレ超セルの非線形光学応答を正確にモデル化することは、計算上の課題となる。第一原理密度汎関数理論（DFT）計算は、比較的小さなねじれ角であっても数千の原子を含む大規模なモアレ超セルの計算において、計算コストが極めて高くなる。したがって、大規模な系のシフト電流に対するねじれ誘起制御を探索するために、これらの系の本質的な電子および幾何学的特徴を捉えることができる、効率的な計算フレームワークが必要とされている。

手法
著者らは、H-MoS₂系を調査するために、第一原理DFT計算からパラメータ化されたスレーター・ケスター（SK）タイトバインディング（TB）モデルを開発し、検証した。

• DFTベンチマーキング: 単層、二層、および三層のH-MoS₂に対して、Quantum ESPRESSOパッケージを用いてDFT計算を行った。TBパラメータのフィッティングを容易にするため、Mo 4dおよびS 3p軌道を用いて最大局在化ワニエ関数を構築した。
• タイトバインディングモデル: TBハミルトニアンはTBPLaSパッケージを用いて構築された。このモデルには、層内ホッピング（Mo–Mo, Mo–S, S–S）および層間相互作用が含まれる。ねじれた系については、AAスタッキング構成に対する指数関数的な距離依存性に従って層間SKパラメータをスケーリングし、モアレ格子内の原子間距離の変化を考慮するために滑らかなカットオフ関数を利用した。
• シフト電流計算: TBモデルに対してはWannierBerriパッケージを、DFTデータに対してはpostW90を用いて、シフト電流導電率テンソルを計算した。計算には150 × 150 × 1 k点のk点グリッドを用いた。先行研究において、MoS₂のシフト電流プロファイルへのスピン軌道相互作用の影響は無視できることが示されているため、スピン軌道相互作用は省略した。
• 対象範囲: 非ねじれの単層、二層、および三層、ならびに、モアレ格子ベクトルを定義する整数対から導出された特定のねじれ角（例：21.787°, 13.174°, 9.430°, 7.341°）を持つ、整合的なねじれ二層および三層を対象とした。

主な貢献と結果

1. モデルの検証: SK-TBモデルは、単層、二層、および三層のH-MoS₂について、DFTおよびWannier補間結果と比較した際、バンド端付近の電子バンド構造およびシフト電流導電率の主要なスペクトル特徴を正確に再現することが示された。このモデルは、追加のマルチレイヤー系へのフィッティングを行うことなく、太陽エネルギースペクトルの範囲における主要なピークの大きさおよびスペクトル位置を捉えることに成功した。
2. 層依存の増強: ねじれていないマルチレイヤー系では、シフト電流の大きさは層数におよそ比例することが判明した。これは、層間結合が電荷中心の位置を変化させ、結果として非線形応答の全体的な増強につながることを示唆している。
3. 対称性の破れとテンソルの活性化: ねじれ角の導入は、AAスタッキング構成の結晶対称性を破る。ねじれていない二層系では、対称性により特定の導電率テンソル成分（例：$\sigma_{yyy}^{sc}$）はゼロとなる。本研究は、ねじれによって、これら以前は禁止されていた成分が活性化されることを示した。ねじれ角が21.787°の場合、$\sigma_{yyy}^{sc}$成分は非ゼロとなり、ねじれ角とともに増大し、最終的には支配的な$\sigma_{xxx}^{sc}$成分と同程度の大きさになる。
4. ねじれ角によるチューナブル性: シフト距離（電荷中心の変位を示す実空間の尺度）の解析により、電子波関数のねじれによる修飾と非線形応答の増加との間の直接的な関連性が明らかになった。ねじれていない場合にゼロであった$y$方向のシフト距離は、検討した最大の角度において約1 Åに達するまで、ねじれ角とともに増加する。
5. 三層構成: 三層のねじれに関する解析も拡張した。中間層のみがねじれた構成では、$x$成分と$y$成分の両方への寄与が増加し、異方性が強化された。逆に、トップ層とボトム層が等しく反対の角度（$\pm\theta$）で回転している構成では、中央層の元の面内対称性が保持されるため、元の非ゼロのテンソル成分のみが残った。

意義と主張
本論文は、ねじれ角を、マルチレイヤーTMDベースの系におけるシフト電流生成をエンジニアリングするための効果的なパラメータとして確立している。主な意義は、適切にパラメータ化されたタイトバインディング法が、大規模なモアレ材料の非線形光学現象を探索するための実用的かつ計算効率の高い経路を提供し、従来の第一原理計算の限界を回避できることを示した点にある。

著者らは、モアレ系が2D材料における非線形光起電力現象を制御するための有望な戦略であることを、その結果に基づき主張している。具体的には、相対的な層の回転を通じて面内の光電流成分を活性化し、調整できる能力は、チューナブルな非線形応答を設計するための新しい自由度を提供する。本研究は、これらの予測に対する専用の実験的検証を求めており、導入された手法が、バルク光起電力用途のための、ねじれたファンデルワールス材料のより大きなモアレ超格子にも適用可能であることを示唆している。今後の課題として、ヘテロ二層系や相関モアレ系への展開が挙げられている。