First-principles Floquet analysis from real-time propagation

技術要約：実時間伝播からの第一原理フロケ解析

問題提起
現実的な材料に対して第一原理シミュレーション内でフロケ理論を実用的に適用するには、重大な計算および解釈上の障壁が存在する。従来のアプローチでは、電子の自由度と光子の自由度を巨大なヒルベルト空間へと結合させた、拡張されたフロケ・ハミルトニアンを明示的に構築する必要がある。この構築は膨大な計算オーバーヘッドを招き、さらにフーリエ級数の強制的な打ち切り（最大フロケ成分数 $n_{max}$ の選択）を必要とする。加えて、標準的な第一原理手法である時間依存密度汎関数理論（TDDFT）は、フロケ固有状態ではなく、自然に実時間伝播された波動関数を出力する。さらに、実験的なシナリオでは、厳密な時間周期性を欠く有限パルスや過渡現象が伴うことが多く、これらは厳密なフロケの枠組みが直接的には適用できない条件であるが、ポンプ・プローブ実験においては「フロケ的な」電子構造が頻繁に観測されている。

手法
著者らは、拡張されたフロケ・ハミルトニアンの構築を回避し、伝播された波動関数から直接準エネルギーとフロケ状態を抽出する、実時間フロケ解析フレームワークを提案している。この手法は、主に以下の3つのコンポーネントで構成される：

1. 一周期進化演算子の再構成：
   静的なフロケ・ハミルトニアンを対角化する代わりに、本手法は伝播された状態間の重なりを用いて、一周期の時間進化演算子 $U(t+T, t)$ を再構成する。時刻 $t$ における正規直交基底 $\{|\psi_i(t)\rangle\}$ （例：コーン・シャム軌道）が与えられたとき、進化演算子の行列要素は次のように計算される：
   $$U_{ij}(t+T, t) = \langle \psi_i(t) | \psi_j(t+T) \rangle$$
   この行列を対角化することで、フロケ準エネルギー ($E_\alpha = -\arg(\lambda_\alpha)/T$) と、伝播された状態のサブスペース内における周期的なフロケ状態の成分 ($|\phi_\alpha(t)\rangle$) が得られる。このアプローチは、既存の時間伝播スキームに対して無視できる程度の計算コストしか追加しない。

2. 調和分解によるフロケバンドのアンフォールディング（展開）：
   （準エネルギーが駆動周波数 $\Omega$ モジュロで定義される）縮小ゾーン表現の曖昧さを解決するために、著者らはアンフォールディング手順を導入している。彼らは、フロケ状態のゼロ次高調波成分 ($n=0$) のノルムに基づいた重み $w_\alpha$ を定義する：
   $$w_\alpha = |u^{(0)}_\alpha|$$
   ここで、$|u^{(0)}_\alpha\rangle$ は周期的なフロケ状態 $|\phi_\alpha(t)\rangle$ をゼロ次高調波へ投影することで得られる。この重みは、あるフロケ状態が基礎となる平衡バンドの性質をどの程度強く保持しているかを定量化する。この展開手法により、物理的に重要な枝（大きな $w_\alpha$ を持つもの）を、高次のフォトンのレプリカ（小さな $w_\alpha$ を持つもの）から明確に区別することができる。

3. サイドバンドの対称性解析：
   再構成された実空間フロケ波動関数を用いることで、個々の光誘起サイドバンドの対称特性を直接解析することが可能になる。波動関数の空間分布を調べることにより、特定のフロケ・サイドバンド（$n=0, \pm 1, \dots$）の対称性（例：滑り鏡映対称性）を決定し、光学的な選択則を検証することができる。

主な結果
本フレームワークは、4つの異なる材料系に対する第一原理TDDFTシミュレーションを通じて検証された：

• グラフェン (2D): 円偏光（CPL）下において、本手法はディラック点におけるトポロジカルなギャップ開口を正常に特定した。アンフォールディングの手順は不可欠であることが示された。これを用いない場合、第一フロケブリルアンゾーン（FBZ）内のスペクトルは、折り畳まれたバンドが密集した判読不能な状態として現れる。重み付きのアンフォールディングは、ハイブリダイゼーション・ギャップとトポロジカル・ギャップを明確に明らかにした。
• 単層黒リンス (2D): 本手法は、光誘起サイドバンドの対称特性を解明する能力を示した。アームチェア配向の駆動下において、解析の結果、$n=0$ の価電子帯および伝導帯は平衡状態の滑り鏡映対称性（それぞれ奇および偶）を保持している一方で、$n=\pm 1$ のサイドバンドは逆の対称性を示すことが確認された。これにより、対称性によって許容されるハイブリダイゼーションによるバンド反発の直接的な説明が可能となった。
• 単層六方晶窒化ホウ素 (2D): 本アプローチは、厳密な時間周期性を超えた有限のCPLパルスにも適用された。パルス内の異なる時間窓（第2、第5、第8サイクル）で解析を行うことで、過渡的なフロケ特徴を捉えることに成功した。結果は、谷依存のバンド繰り込み（左円偏光に対して$-K$谷でより強い）が、パルスエンベロープに従って出現および消失することを示しており、過渡的な非平衡レジームにおける本手法の有用性を立証した。
• バルク硫化錫 (3D): 本手法は、完全に周期的な3Dバルク半導体へと拡張された。シミュレーションは、偏光依存のバンド繰り込みを明らかにした。具体的には、アームチェア配向の駆動は$\Gamma X$谷に影響を与え、ジグザグ配向の駆動は$\Gamma Y$谷に影響を与えており、これは従来の対称性に基づく解析と一致している。

意義と主張
本論文は、実時間シミュレーションとフロケ観測量を結びつけるための、実用的かつ計算効率の高いフレームワークを提供すると主張している。その主な貢献は以下の通りである：

1. 計算効率: 拡張されたフロケ・ハミルトニアンの構築と対角化という高コストなプロセスを回避し、標準的な時間伝播中に既に計算された状態の重みのみを利用する。
2. 解釈性: 高調波の重みに基づくアンフォールディング・スキームにより、折り畳まれたフォトンのレプリカの曖昧さを解決し、物理的なバンド端やハイブリダイゼーション・ギャップの識別を容易にする。
3. 汎用性: 本手法は2D単層材料から3Dバルク材料まで適用可能であり、厳密な時間周期性が存在しない有限パルス条件下にも拡張できる。
4. 分析の深化: 実空間の波動関数から、個々の光誘起サイドバンドの対称特性や選択則を直接抽出することを可能にする。

著者らは、フロケ理論が周期駆動された定常状態を構築するためのツールとしてだけでなく、現実的な材料における過渡的な準粒子ダイナミクスや非平衡電子構造のための強力な解析フレームワークとしても機能すると結論付けている。