\textit{\textbf{First-principles}} description of pumped inelastic X-ray scattering: example of K-edge RIXS in graphite

本論文は、光励起された物質における時間分解共鳴非弾性X線散乱（RIXS）を予測するために、ベテ・サルピーター方程式と実時間時間依存密度汎関数理論を組み合わせた第一原理フレームワークを提示し、グラファイトにおける様々な遅延時間での角度依存K吸収端RIXSスペクトルのモデリングに成功することによってその正確性を実証するものである。

要約

複雑な機械、例えばグランドピアノがどのように動いているのかを理解しようとしている場面を想像してみてください。あなたは曲の演奏を聴くことができます（これは通常の分光法にあたります）。しかし、**共鳴非弾性X線散乱（RIXS）**は、光で作られたハンマーで特定の鍵盤を叩き、その際に内部の弦やハンマーがどのように振動して音が響いたのかを詳細に分析するようなものです。これは、その機械が何でできているかだけでなく、その部品がどのように動き、相互作用しているのかまでも教えてくれます。

この論文では、この「光のハンマー」による実験がどのようなものになるかを正確に予測する、極めて精密な新しいコンピュータプログラムを紹介しています。しかも、その材料が2番目のより速い光のパルス（ポンプ）によって揺さぶられている場合でもです。

以下に、日常的な比喩を用いた彼らの研究の解説をまとめます。

1. 問題点：予測不可能なものの予測

科学者たちは、長年、X線を用いて材料の「スナップショット」を撮ることができました。しかし、材料が「眠り」から覚めるようにレーザーで「ポンプ（励起）」されているときに、そのスナップショットが正確にどのような見た目になるかを予測することは非常に困難でした。

• 従来の方法： 以前のコンピュータモデルは、群衆を見て、全員が一人で静止して行動していると仮定しているようなものでした。彼らは、人々（電子）が実際にどのように手を取り合い、共に動いているか（「エキシトン効果」と呼ばれる現象）を見落としていました。
• 新しい方法： 著者らは、高速な3D映画シミュレーターのように機能する新しいフレームワークを構築しました。これは単に個々の動きを見るのではなく、人々が互いに引き合いながらどのように一緒に踊っているか（群衆のダンス）を観察するものです。

2. 手法：2ステップのダンス

研究者たちは、シミュレーションを作成するために2つの強力なツールを組み合わせました。

• ステップ1（「ポンプ」）： 彼らはRT-TDDFTと呼ばれるツールを使用して、レーザーが材料に当たったときに何が起こるかをシミュレートしました。トランポリンに懐中電灯の光を当てる様子を想像してください。このツールは、光が当たった直後にトランポリンがどのように跳ね、その上の人々がどのように体重を移動させるかを計算します。これにより、レーザーパルスの直後に電子がどこにあるかという「非平衡」なマップが得られます。
• ステップ2（「プローブ」）： 次に、彼らは**ベテ・サルピーター方程式（BSE）**を使用しました。これは、X線がその跳ねているトランポリンとどのように相互作用するかについての、非常に正確なルールブックのようなものです。これは、弾き飛ばされた電子と、その電子が残した「穴（空席）」との間の複雑なダンスを計算します。

これらを組み合わせることで、光が入ってくるあらゆる角度と、光が出ていくあらゆる角度に対する「エコー（散乱X線）」を予測することができます。

3. テストケース：グラファイト（鉛筆の芯）

彼らの手法が機能することを証明するために、彼らはグラファイト（鉛筆の芯の成分）でテストを行いました。

• なぜグラファイトなのか？ これは紙の束のようなものです。各シート内の原子は強い接着剤で固められていますが（強い接着剤のように）、シート同士は緩くつながっています（積み重ねられた緩い紙のように）。このため、グラファイトは非常に「異方性」が高く、つまり、横から見るか上から見るかによって挙動が大きく異なります。
• 結果： コンピュータ・シミュレーションは、グラファイトが奏でる2種類の明確な「音符」を予測することに成功しました。
  • $\pi$（パイ）の音符： シート間を移動する電子から発生します（緩い紙のように）。
  • $\sigma$（シグマ）の音符： シート内でタイトに動く電子から発生します（強い接着剤のように）。
    シミュレーションは、横から光を当てると主に「接着剤」の音が聞こえ、上から光を当てると「紙」の音が聞こえることを示しました。これは実世界の実験結果と完璧に一致しました。

4. 「ポンプされた」実験：テーブルを揺らす

この論文で最もエキサイティングな部分は、X線を当てる前にグラファイトをレーザーで「ポンプ」したときに何が起こるかです。

• 比喩： グラファイトを穏やかな池だと想像してください。レーザーポンプは、池に石を投げ込み、波紋を作るようなものです。X線は、その波紋によって水がどのように変化したかを確認するために送られるソナーの信号です。
• 発見： グラファイトが「ポンプ」されると、シミュレーションは「音符」がわずかに変化することを示しました。低エネルギー領域に新しい、かすかな音が現れ、既存の音のボリュームが変化しました。
• 教訓： コンピュータは、短いレーザーパルスであっても、材料の電子的な「気分」を変化させ、静止状態とは異なる一時的な状態を作り出すことを予測しました。シミュレーションは実験データと非常によく一致しており、これらの微妙な変化さえも捉えることができました。これは、手法が「時間分解（映画のフレームごとの動き）」研究に有効であることを証明しています。

まとめ

簡単に言えば、この論文は次のように述べています。「私たちは、材料がレーザーによって激しく揺さぶられているときでも、その材料がX線に対してどのように反応するかを正確に予測できる、非常に精度の高い新しいコンピュータモデルを構築しました。」

彼らはグラファイトを用いてテストを行い、コンピュータの「予測」が実世界の「実験」と完璧に一致し、材料の内部構造（タイトなシート vs 緩い層）が異なる角度や異なる時間において光にどのように反応するかを正しく特定したことを示しました。これにより、科学者はあらゆる仮説に対して高価な実験を行う必要なく、材料がリアルタイムでどのように振る舞うかを理解するための強力な新しいツールを手にすることになります。

技術概要

技術要約：励起された非弾性X線散乱の第一原理記述

問題と背景
共鳴非弾性X線散乱（RIXS）は、元素、軌道、および運動量分解能を備えた電子励起を探索するための強力な分光手法である。RIXSの理論的枠組みは、独立粒子近似から、エキシトン効果を捉えるためのベテ・サルピーター方程式（BSE）を用いた多体摂動論（MBPT）へと進歩してきたが、ポンプ励起された（pumped）時間分解（time-resolved）RIXSを記述することには依然として大きな隔たりが存在する。既存の理論モデルは、以下の事項を同時に考慮する能力に欠けることが多い：

1. 正確なスペクトル記述に不可欠な、コアおよび価電子励起の両方におけるエキシトン効果。
2. 光ポンプによって誘起される非平衡電荷キャリア分布。
3. 異方性材料において不可欠な、入射および放出フォトンの任意の偏光幾何学。

時間分解RIXS（tr-RIXS）実験は台頭しつつあるが、理論的なモデリングには、電子-正孔相関と非平衡ダイナミクスを統合する枠組みが必要である。

手法
著者らは、偏光および時間分解RIXSスペクトルを予測するために、**実時間時間依存密度汎関数理論（RT-TDDFT）とベテ・サルピーター方程式（BSE）*を組み合わせたab initio*（第一原理）フレームワークを提示している。その手法は以下の通りである：

1. 理論的定式化： このアプローチは、二重微分断面積（DDCS）に関するクラマース・ハイゼンベルク公式に基づいている。RIXSプロセスは、2段階のメカニズムとして記述される：

   • 吸収： 入射X線光子（$\omega_1$）がコアホール中間状態（$|n\rangle$）を生成する。
   • 放出： 価電子がコアホールを埋め、光子（$\omega_2$）を放出し、最終的な電子-正孔対（$|f\rangle$）を残す。
     DDCSは、BSEの固有ベクトルおよび固有値から導出された遷移演算子を用いて計算される。

2. 非平衡状態の統合： ポンプされたケースをモデル化するために、本フレームワークは、光学ポンプの下でのコーン・シャム系を伝播させるためにRT-TDDFT（具体的にはexcitingコード内）を採用している。これにより、時間依存の非平衡電子占有数が得られる。

3. ワークフロー：

   • 基底状態： DFT計算が出発点となる。
   • BSE計算： 以下の2つの独立したBSE計算が行われる：
     • コアレベル励起用（コア振動子強度 $t^{(1)}$ を得るため）。
     • 価電子/光学励起用（最終状態の経路 $t^{(2)}$ を得るため）。
   • 非平衡拡張： tr-RIXSについては、RT-TDDFTから得られた非平衡占有数がBSE計算における制約条件として使用される。
   • 偏光の一般化： 実装は、入射光（$\mathbf{e}_1$）と放出光（$\mathbf{e}_2$）の任意の独立した偏光ベクトルを扱えるように拡張されており、従来の $\mathbf{e}_1 = \mathbf{e}_2$ という制限を超えている。

4. 計算の詳細： 計算は、全電子フルポテンシャルパッケージであるexcitingを用いて行われた。研究では、van der Waals補正を伴うGGA-PBE汎関数を用いたグラファイトの炭素K端が利用された。BSEには、光学スペクトルのために占有バンド8および非占有バンド15を、X線スペクトルのために伝導帯40個を含めた。

主な貢献

• 統一されたフレームワーク： 本論文は、RT-TDDFTとBSEを組み合わせ、エキシトン的な多体効果とポンプ誘起の非平衡ダイナミクスの両方を捉える、完全にパラメータフリーなab initioワークフローを確立した。
• 偏光の一般化： 著者らは、入射光と散乱光の両方の任意の偏光方向をサポートするようにRIXSの実装を一般化し、複雑な散乱幾何学や混合偏光の調査を可能にした。
• excitingへの実装： この手法は、フルポテンシャルコードexciting内に実装されており、断面積計算にはpyBRIXSパッケージを利用している。

結果：グラファイトK端RIXS
このフレームワークは、面内（$\sigma$）成分と面外（$\pi$）成分の間で強い電子異方性を持つ半金属であるグラファイトに適用された。

• 平衡スペクトル：

  • 計算された吸収スペクトル（光学およびコアレベルの両方）は実験的特徴を再現しており、特に $\pi \to \sigma^*$ および $\sigma \to \pi^*$ 遷移を捉えるためにBSEが必要であることを裏付けた。
  • 異方性： 平衡状態のRIXSスペクトルは強い偏光依存性を示した。面内偏光（$\mathbf{e} \parallel x$）は高エネルギーの$\sigma$由来の励起（$>8$ eV）を強調し、面外偏光（$\mathbf{e} \parallel z$）は低エネルギーの$\pi$由来の励起を強調した。
  • 散乱幾何学： 実験セットアップと一致する入射角$30^\circ$のシミュレーションでは、$\pi$と$\sigma$の寄与の非線形な混合が見られ、干渉効果と低エネルギー損失領域での $\pi \to \pi^*$ 遷移の増強が示された。

• ポンプされた（非平衡）スペクトル：

  • シミュレーションは、266 nmのレーザー（持続時間45 fs）によってポンプされ、遅延時間150 fsの状態で行われた。
  • スペクトルの変化： ポンプされたスペクトルは平衡状態のケースと酷似しており、修正は緩やかであった。しかし、両方の偏光において、明確な追加の低エネルギー損失特徴（数eV）が現れた。これは、フェルミ準位付近の電子分布の変化に起因する。
  • 混合偏光： 混合偏光の設定では、非平衡スペクトルは平衡状態と比較して、スペクトル強度の再分配と特徴の広がりを示した。
  • 実験との比較： 理論結果は、実験のピーク位置と一般的な特徴を良好に再現した。微小な偏差（例：ピーク幅や特定のシフト）は、現在の静的な枠組みには含まれていない電子-フォノン結合や、実験における自己吸収効果に起因すると考えられる。

意義と主張
著者らは、自らの研究が、複雑な材料（金属から大きなギャップを持つ半導体まで）に対して、任意の散乱幾何学でRIXSスペクトルを計算するための、一般的かつ完全にパラメータフリーな経路を提供するものであると主張している。グラファイトの特性的な異方性と、光励起下でのその進化を再現することに成功したことで、本手法は以下の能力を実証した：

1. 明確に区別された$\pi$および$\sigma$由来の特徴を解明すること。
2. 電子構造と偏光の相互作用を捉えること。
3. 光励起によって誘起される時間依存現象をモデル化すること。

本論文は、このアプローチが、固体における非平衡現象を理解するための、tr-RIXS実験を解釈するための堅牢なツールを提供することを結論付けており、特に基本励起の方向依存性を強調している。著者らは、電子-フォノン結合や自己吸収効果の除外といった既知の制限に起因する実験との不一致について、控えめな姿勢を保っている。