Dielectric function in WSe2

技術要約：WSe2単層ウィグナー結晶における誘電関数

問題提起
本論文は、2次元（2D）遷移金属ダイカルコゲナイド（TMD）単層、特に単層WSe2に焦点を当て、ウィグナー結晶状態の理論的な特性評価に取り組んでいる。低密度電子ガスにおけるウィグナー結晶化の存在は理論的に確立されており、近年TMDヘテロ構造において実験的に観測もされているが、これらの結晶状態の電子バンド構造および静的誘電応答を計算するための微視的な枠組みが必要とされている。課題は、連続的な並進対称性が失われた破れた対称性状態において、運動エネルギー、クーロン相互作用（特に2D材料に関連するケルディッシュ相互作用）、および交換効果の相互作用を正確にモデル化することにある。

手法
著者らは、完全にスピン偏極した2Dウィグナー結晶のバンド構造および静的誘電関数を計算するために、自己整合的なハートリー＝フォック（HF）数値手法を開発した。主要な手法構成要素は以下の通りである：

• ハートリー＝フォック定式化： 計算は、基底状態が完全にスピン偏極していると仮定した制限付きHF枠組み内で行われる。単一粒子軌道は、基本となる三角形状のウィグナー結晶格子の第1ブリルアンゾーンに属する波数$k$と、逆格子ベクトル$G$を用いた平面波基底$|k + G\rangle$によるブロッホ形式に展開される。
• 相互作用ポテンシャル： 電荷キャリア間の有効静電相互作用は、周囲の媒体による誘電遮蔽と2D材料の性質を考慮したケルディッシュ相互作用を用いてモデル化されている。
• 自己整合ループ： アルゴリズムは、周期的なガウス型電荷密度から初期化される。その後、密度高調波から導出されるハートリーポテンシャルと、占有固有ベクトルから導出されるフォック交換行列を反復的に更新し、収束するまでループを回す。密度高調波は、更新された固有ベクトルから再構成され、ループを閉じる。
• 数値実装：
  • ブリルアンゾーンのサンプリングには、対象となる運動量に対して既約$\Gamma-K-M$ウェッジを利用し、交換項は対称性による像から再構成された全ゾーンの総和によって評価される。
  • 逆空間における円形カットオフにより平面波基ースを選択し、三角格子が持つ$D_6$点群対称性を保持する。
  • クーロン/ケルディッシュポテンシャルの$q=0$における特異点は、逆空間における小さな円形セル内での解析的な平均化によって処理される。
• 誘電応答の計算： 静的誘電関数$\epsilon(q, 0)$は、絶縁体である折り畳みブリルアンゾーンの性質に適応させたリンドハード型の分極率公式を用いて計算される。これには、収束したHFバンド間の占有バンドと空バンド間の遷移に関する和が含まれ、収束した固有ベクトルと固有値を利用する。

主な貢献

1. 専門的なコードの開発： 著者らは、TMD単層に特化したケルディッシュ相互作用を持つ2Dウィグナー結晶の自己整合的HF方程式を解くことができる数値コードを提示している。
2. 微視的誘電応答： 本研究は、古典的なフォノンモデルや一様電子ガスの近似に頼ることなく、HFバンド構造および固有関数から直接、静的誘電関数を導出する方法を提供している。
3. 結晶秩序の特性評価： 結晶化の程度を定量化するために、以下の特定の指標を導入し利用している：
   • 第1スター秩序パラメータ（$M_1$）： 基本的な密度変調の振幅を測定する。
   • 密度コントラスト比（$C_n$）： 密度の極大値と極小値の間の実空間上の分離を定量化する。
   • 高対称点（$\Gamma, K, M$）におけるバンド構造の進化： 三角格子対称性のスペクトル的な署名として利用される。

結果
著者らは、弱結晶状態から強く局在化した状態への遷移を表す、さまざまな密度パラメータ（$r_s$）にわたる単層WSe2に対して彼らの手法を適用した：

• バンド構造の進化： $r_s$が増加する（密度が減少する）につれて、系は弱く変調された状態から堅牢なウィグナー結晶へと遷移する。高い$r_s$において、バンド構造は六角形ブリルアンゾーンの対称性と一致する明確な縮退パターンを示し、エネルギーバンドは$\Gamma-K-M$経路に沿って滑らかに進化する。
• 密度変調： 秩序パラメータ$M_1$は$r_s$とともに増加し（例：$r_s=4$で0.100から$r_s=50$で0.586へ）、これは電荷の局在化が強まっていることを示している。密度コントラスト比$C_n$も大幅に増加しており、格子サイトにおける明確な電荷密度ピークの形成を裏付けている。
• 誘電関数： 静的誘電関数$\epsilon_{HF}(q, 0)$は、$r_s$の増加に伴いピーク値が系統的に減少することを示す。これは電子の分極率の低下を反映している。電子がより局在化するにつれて、外部摂動に対する遮蔽効率が低下し、遮蔽された相互作用が裸のケルディッシュ相互作用に近づくためである。
• 収束性： 自己整合的な解は安定性を示し、最終的な密度プロファイルは、運動エネルギー、ハートリーエネルギー、および交換エネルギーのバランスによって決定される状態へと、初期のガウス近似から緩和される。

意義および主張
本論文は、低密度2D系における遮蔽、バンド構造の再構成、および相互作用効果を調査するための理論的ツールとして位置づけられている。著者らは、彼らの手法が、弱結晶解から局在化したウィグナー結晶状態への進化に関する一貫した微視的な記述を提供すると主張している。

本研究の意義は、原子層レベルの薄さを持つ半導体の最近の進展という文脈の中に置かれている。著者らは、計算されたスピン保存、凍結スピン/バレーセクターが、光ポンプによって偏極状態へと駆動され得るゲート制御可能なWSe2単層に対する、時間分解・ヘリシティ分解光分光法によってテスト可能であると示唆している。しかし、彼らは自らの結果が、熱力学的な液体・結晶転移の証明ではなく、あくまで「位相感受的な指標」として機能することを明示しており、決定的な転移の特定にはエネルギー比較や独立した相安定性解析が必要であることを注記している。本研究は、従来の低密度2Dウィグラー結晶理論と、TMDにおける電荷秩序の最近の光学研究を結びつけ、単層TMDの特定のケルディッシュ相互作用および光応答に適応した枠組みを提供している。