Wigner crystallization in Bernal bilayer graphene

技術要約：ベルナル・ビレイヤーグラフェンにおけるウィグナー結晶化

問題提起
本論文は、垂直方向の変位電場（$U$）の影響下における、ベルナル・ビレイヤーグラフェン（BBG）におけるウィグナー結晶（WC）状態の性質と安定性を調査している。ねじれ（ツイスト）および非ねじれのグラフェン多層構造における近年の実験では、豊かな相関相が明らかになっているが、BBG特有のバンド構造（具体的には$U$によって誘起される「メキシカンハット」型分散）と、ベリー曲率、およびトリゴナル・ワーピング（三回対称性の歪み）がWC相に与える相互作用は、依然として未解決の理論的課題である。著者らは、電子密度（$n$）、変位電場（$U$）、および温度（$T$）の観点からWCの相境界を決定し、従来のWC相とは異なる、軌道磁化やバレー秩序を含む基底状態の特性を明らかにすることを目的としている。

手法
著者らは、WC状態をモデル化するために、単一粒子調和振動子（HO）記述を採用している。この枠組みにおいて、三角格子状のWC内の各電子は、隣接する電子によるクーロン反発によって形成される放物線型のポテンシャル井戸に閉じ込められた粒子として扱われる。この閉じ込め強度は電子密度によって決定される。

研究は主に以下の2つの理論的領域で進行する：

1. 回転対称モデル（トリゴナル・ワーピングなし）： 著者らは、低エネルギー分散を$K$および$K'$点に対して回転対称として扱うため、トリゴナル・ワーピングを無視する。これにより、ベクトルポテンシャルとしてのベリー接続（$A(p)$）を組み込んだ、有効な単一粒子シュレディンガー方程式を運動量空間で導出する。これにより、エネルギースペクトル、波動関数、およびリンデマン基準の解析的な推定が可能となる。
2. トリゴナル・ワーピング・モデル： 著者らは、回転対称性を破り、「メキシカンハット」型分散を離散的な「ミニバレー」（中央のポケットと3つのサイドポケット）へと分裂させるトリゴナル・ワーピングの効果を導入する。この場合、リンデマン比を計算し、相図を描画するために、運動量空間グリッド上での二次元シュレディンガー方程式の数値対角化を行う。

WCの安定性は、電子の波動関数の二乗平均平方根幅が格子定数の臨界割合（$\eta_c \approx 0.23$）を超えたときに融解が起こるという、リンデマン基準を用いて評価される。また、交換積分を計算して秩序温度を推定し、軌道角運動量演算子の期待値を通じて磁化を評価する。

主要な貢献と結果

• 変位電場依存性と誘電スクリーニング： 本論文は、変位電場がゼロの場合、低エネルギーでのバンド質量が消失し、強力なバンド間誘電スクリーニングが生じるため、BBGにおけるWC状態は阻止されることを確立している。WC相は、バンドギャップを開き、バンド底を平坦化するために十分な変位電場が印加されたときにのみ出現する。
• 自発的軌道磁化（トリゴナル・ワーピングなし）： トリゴナル・ワーピングがない場合、著者らはベリー曲率がWC電子の軌道角運動量（$\ell$）の遷移を駆動することを実証している。変位電場$U$が増加するにつれて、波動関数が囲むベリーフラックスが増大する。このフラックスが臨界値（$\tilde{\Phi} > 1/2$）を超えると、基底状態は$\ell=0$から$\ell=1$（または$K'$バレーにおける$\ell=-1$）へと遷移する。この遷移は、WC状態における軌道磁化の不連続なジャンプとして現れる。
• 二重再入融解（トリゴナル・ワーピングあり）： トリゴナル・ワーピングを含めると、低$U$における電子密度関数として、相図は特有の「二重再入（doubly re-entrant）」挙動を示す。
  1. 極めて低密度では、電子は低エネルギーの「サイドポケット」に結晶化する。
  2. 密度が増加すると、電子は臨界密度$n_{side}^c$でフェルミ液体（FL）へと融解する。
  3. より高い密度$n^*$において、電子がより重い質量の「中央ポケット」を占有することがエネルギー的に有利となり、システムは再びWC状態へと再凍結する。
  4. 最終的に、さらに高い密度$n_{central}^c$において、この中央ポケットのWCは再び融解する。
• ミニバレー秩序： トリゴナル・ワーピングの存在は回転対称性を破るため、角運動量$\ell$は定義が不明確になる。代わりに、基底状態は非自明な空間パターンである「ミニバレー秩序」によって特徴付けられる。著者らは、最近の研究（文献[102]）に基づき、電子がWC格子上で交互のストライプパターンとして特定のミニバレーへと偏極することを予測している。
• 相境界： 著者らは、WCの相境界に関する数値的な推定値を提供している。BBGにおけるWC状態の最大密度は、変位電場が約$0.7 \text{ V/nm}$のとき、およそ$1.2 \times 10^{10} \text{ cm}^{-2}$と推定される。最大融解温度は、密度とともに減少するものの、およそ$10 \text{ K}$程度と推定され、再入する中央ポケット領域の最大温度は$\sim 1 \text{ K}$と低い。

意義と主張
本論文は、BBGにおけるWC状態の包括的な理論的相図を提供し、非自明な量子幾何学（ベリー曲率）と格子対称性の破れ（トリゴナル・ワーリング）がいかにしてウィグナー結晶の性質を根本的に変えるかを強調している。

• 軌道効果： 本研究は、スピンに基づく秩序とは明確に異なる、純粋にベリー曲率と変位電場によって駆動される自発的軌道磁化のメカニズムを特定した。
• 複雑な相挙動： 「二重再入」融解の予測は、異なる有効質量とエネルギーを持つ異なるバンドポケット（中央 vs サイド）の間の競合から生じる新しい特徴として提示されている。
• 実験的背景： 著者らは、BBGにおけるWCの予測密度範囲（$n \lesssim 10^{10} \text{ cm}^{-2}$）が、近年の実験で観察されている絶縁状態の密度よりも低いことに触れ、現在のBBFにおける実験的観察は、無秩序（disorder）に支配されているか、あるいは異なる領域に対応している可能性を示唆している。彼らは、再入挙動やミニバレー秩序は、バンド構造がより平坦でベリーフラックスがより大きい（$N\pi$）、層数（$N$）が多いロムボヘドラル積層グラフェンにおいて、より観察しやすい可能性があると示唆している。
• 補完性： 本論文は、スピン秩序に焦点を当てた随伴論文（文献[1]）を補完するものであると明示されており、本論文は軌道効果と、トリゴナル・ワーリングを含む詳細な相図に焦点を当てている。

著者らは、HO近似の限界について謙虚な姿勢を維持しており、この近似がより厳密な量子モンテカルロ計算と比較してWC相の安定性を過大評価している可能性があること、また、ここで用いられている「アインシュタイン・フォノン」の極限は、非自明なチャーン数を持つ相（量子異常ホール結晶など）を再現できないことを認めている。