Compton散乱プロファイルの第一原理計算：運動量空間に見る化学結合・Fermi面・スピン

X線回折が実空間の電子密度 $\rho(\mathbf{r})$ を与えるのに対し、高エネルギーX線のCompton散乱は運動量空間の電子密度 $\rho(\mathbf{p})$ を直接測定できる、数少ない実験手法です。運動量密度は波動関数そのものの広がりを反映するため、共有結合の方向性、金属のFermi面形状、強磁性体のスピン磁気モーメントといった「価電子の素性」に非常に敏感です。加えて、100 keVを超える硬X線を用いるためバルク敏感かつ非破壊であり、近年では動作中（in operando）のリチウムイオン電池内部のLi濃度分布や、正極材料の酸化還元軌道の観測にも応用されています。ただし実験で得られるのは3次元の $\rho(\mathbf{p})$ そのものではなく、散乱ベクトル方向へ2重積分したComptonプロファイル $J(p_z)$ という1次元量です。したがって、観測されたプロファイルから三次元的な電子状態の知見を正確に抽出するためには、第一原理計算との定量的・直接的な比較検証が重要となります。本解析事例では、第一原理計算ソフトウェア Advance/PHASE (ver.5.0) に実装されたCompton profile計算機能（PAW全電子再構成に基づく実装）を用い、性質の異なる3つの物質——Si（共有結合性半導体）、Li（単純金属）、Fe（強磁性遷移金属）——について方向別Comptonプロファイルを計算しました。プロファイル計算における総和則の有効性、k点メッシュに起因するエイリアシングと実験分解能の適切な関係付け、およびFeにおけるスピン磁気モーメントの高い再現性についても定量的に示します。

Keywords: 第一原理計算 (DFT), Compton散乱, 電子運動量密度, Comptonプロファイル, PAW法, 全電子再構成, Fermi面, 磁気Comptonプロファイル, スピン磁気モーメント, リチウムイオン電池

Compton散乱と電子運動量密度

図1はCompton散乱測定の模式図です。入射X線（$\omega_1,\mathbf{k}_1$）が電子によって非弾性散乱される際、散乱X線（$\omega_2,\mathbf{k}_2$）のエネルギーは、静止電子を仮定した際のComptonシフトに伴う基準エネルギーを中心に広がりを持ちます。散乱X線のエネルギー移行が電子の束縛エネルギーより十分大きいとき（衝突近似, impulse approximation）、この基準エネルギーからのずれは、散乱ベクトル $\mathbf{K}=\mathbf{k}_1-\mathbf{k}_2$ に沿った電子の運動量成分 $p_z$ によるDopplerシフトとして観測されます。したがってCompton線の形状は、運動量密度を $\mathbf{K}$ に垂直な面内で積分した

$$\displaystyle J(p_z) = \iint \rho(\mathbf{p})\, dp_x\, dp_y , \qquad \rho(\mathbf{p}) = \sum_{n\mathbf{k}} f_{n\mathbf{k}} \left| \Psi_{n\mathbf{k}}(\mathbf{p}) \right|^2$$

を与えます [1]。ここで $\Psi_{n\mathbf{k}}(\mathbf{p})$ は波動関数の運動量表示です。定義から直ちに総和則

$$\displaystyle \int_{-\infty}^{\infty} J(p_z)\, dp_z = N \quad (N: \text{電子数})$$

が成り立ち、これは計算の正しさを検証する強力な指標になります。また $J(p_z)$ のFourier変換

$$\displaystyle B(r) = \int J(p_z)\, \cos(p_z r)\, dp_z$$

は相反形状因子（reciprocal form factor）と呼ばれ、プロファイル方向に沿った波動関数の自己相関関数に対応します。すなわち、運動量空間の測定量から実空間の結合の情報を取り出す橋渡しとなります。

図1. Compton散乱測定の模式図。(a) 入射・散乱X線と散乱ベクトル $\mathbf{K}$。$\mathbf{K}$ の向きが運動量軸 $p_z$ を定義します。(b) 散乱X線のスペクトル。静止電子なら鋭い線（破線）となるところ、電子の運動量分布によりDoppler広がりを持ちます。(c) 運動量密度 $\rho(\mathbf{p})$。Fermi面（濃い領域の境界）とUmklapp成分が現れます。(d) $\rho(\mathbf{p})$ を $p_x,p_y$ について積分した Comptonプロファイル $J(p_z)$。Fermi運動量 $p_\text{F}$ に折れ曲がり（break）が現れます。

第一原理計算におけるComptonプロファイルの算出と全電子再構成

標準的な擬ポテンシャル法を用いた第一原理計算では、価電子の波動関数は原子核近傍で節を持たない擬波動関数 $\tilde\psi$ に置き換えられています。実空間の電子密度や全エネルギーの算出にはこれで十分ですが、運動量密度は原子核近傍の急峻な波束の構造に強く支配されるため、擬波動関数のまま計算すると高運動量成分がほぼ完全に失われてしまいます。そこで、PAW（Projector Augmented-Wave）法などに基づく現代の計算手法では、以下の変換式

$$\displaystyle \Psi^{\rm AE}(\mathbf{p}) = c_{n\mathbf{k}}(\mathbf{G}) + \sum_{a}\sum_{i} e^{-i\mathbf{p}\cdot\mathbf{R}_a}\, i^{l}\, Y_{lm}(\hat{\mathbf{p}})\, \Delta\phi_{i}(|\mathbf{p}|)\, \langle \tilde p_i | \tilde\psi_{n\mathbf{k}} \rangle$$

を用いて、全電子（all-electron, AE）波動関数の運動量振幅を再構成してから $\rho(\mathbf{p})$ を組み立てます（$\Delta\phi_i$ はAE部分波と擬部分波の差のBessel変換です）[2]。これにより、擬ポテンシャルベースの計算効率を保ちながら、全電子法に匹敵する精度でComptonプロファイルを評価することが可能になります。さらに定量性を高めるためには、サイト内成分の異方性（$(l,m)\neq(l',m')$ 交差項）の考慮や、高運動量側における解析的なテイル補正による総和則の厳密化などが重要となります。

また、計算結果を実験値と適切に比較・解析するためには、いくつかの物理的な条件設定が求められます。代表的な留意点は以下の通りです。

• 射影方向の指定： 実験における散乱ベクトル $\mathbf{K}$ の向きに合わせて、結晶軸に対する射影方向（[100], [110], [111]など）を設定します。
• 実験分解能の考慮（Gauss広がり）： 離散的なk点メッシュに起因するエイリアシング（非物理的な微細振動）を取り除き、かつ実験装置の運動量分解能（半値全幅: FWHM）に合わせるため、計算されたプロファイルに適切なGauss関数による畳み込み（ブロードニング）を行います。
• スピン分解（磁気Compton散乱）： 遷移金属などの強磁性体を対象とする場合、スピン上向き・下向きのプロファイルを個別に算出することで、軌道磁気モーメントを含まない純粋なスピン磁気モーメントの分布（磁気Comptonプロファイル）を評価できます。

計算モデルと計算条件

計算対象の3系はいずれも実験格子定数を用いた完全結晶（プリミティブセル）です。Compton散乱の実験は室温で行われることが多く、運動量密度の比較では実験格子定数を用いるのが標準的なためです。交換相関汎関数はGGA-PBE、擬ポテンシャルはPAWを用いました。主な計算条件を表1に示します。

表1. 計算条件

項目	Si	Li	Fe
結晶構造	ダイヤモンド構造	bcc	bcc
格子定数	5.431 Å	3.491 Å	2.8665 Å
セル内原子数	2	1	1
擬ポテンシャル	PAW（GGA-PBE）
価電子数 / セル	8	1	8
カットオフ（波動関数 / 電荷密度）	25 / 225 Ry	25 / 225 Ry	36 / 324 Ry
k点メッシュ（Monkhorst-Pack）	16×16×16	48×48×48	16×16×16
既約k点数	408	2756	140
スピン分極	なし	なし	あり（強磁性）
射影方向	[100], [110], [111]
Gauss広がり FWHM	0.15 a.u.	0.15 a.u.	0.40 a.u.

※ Feの広がり幅0.40 a.u.は、放射光を用いた磁気Compton散乱実験の運動量分解能（0.4 a.u. 程度）に合わせています [3]。

計算結果と考察

1. Si：共有結合が生む方向異方性

図2にSiの価電子Comptonプロファイルを示します。総和則は $\int J\,dp_z = 7.9932$（目標値8.0000）で、誤差は0.08 %に収まっています。プロファイルは $p_z \lesssim 1$ a.u. に主ピークを持ち、方向依存性は $J(0)$ で数 % 程度です。

図2. Si（ダイヤモンド構造）の価電子Comptonプロファイル。(a) 3方向の $J(p_z)$。太線はGauss広がり（FWHM 0.15 a.u.）を与えたもの、細線は生ヒストグラム。(b) 方向間の差（異方性）。異方性を取ることで等方的な内殻寄与が完全に相殺されるため、実験と直接比較する上で最も信頼できる指標となります。

Comptonプロファイルの絶対値には、実験では内殻電子（Siでは1s, 2s, 2p の10電子/原子）の等方的な寄与が加わります。本計算はPAWの価電子（3s, 3p の4電子/原子）のみを扱うため、絶対値の比較には自由原子の内殻プロファイルを加算する必要があります。一方、方向間の差（異方性）では等方的な内殻寄与が完全に相殺するため、図2(b)の異方性は計算と実験を直接比較できる量になります。Siの異方性は共有結合が $\langle111\rangle$ 方向を向いていることの直接的な現れであり、実験的にもReedとEisenbergerによる測定以来、繰り返し議論されてきた対象です [4]。

この異方性はk点メッシュに対してよく収束しています。$J_{[110]}-J_{[100]}$ の振幅は 12³、16³、20³ のメッシュに対してそれぞれ 0.358、0.355、0.358 a.u. であり、1 %以内の差で一致します。ただし異方性の振幅は運動量分解能（広がり幅）に強く依存する点には注意が必要です。同じ計算結果でも、FWHMを0.15 → 0.30 → 0.40 a.u.と変えると振幅は 0.355 → 0.223 → 0.148 a.u.（原子あたり 0.177 → 0.112 → 0.074）と変化します。実験と比較する際には、必ず実験の運動量分解能に合わせた広がりを計算側にも与える必要があります。

異方性の物理的意味は、$J(p_z)$ をFourier変換した相反形状因子 $B(r)$ を見るとより明瞭になります（図3）。$B(r)$ はその方向に沿った波動関数の自己相関であり、原子間ベクトルの位置に構造を持ちます。

図3. Siの相反形状因子 $B(r)=\int J(p_z)\cos(p_z r)\,dp_z$。(a) 3方向の $B(r)$。全体形状は原子内の自己相関に支配されるため、方向依存性は一見小さく見えます。(b) 方向間の差 $\Delta B(r)$。破線は各方向に沿った最短の原子間ベクトルの長さで、[111]方向は Si-Si 結合そのもの（$\sqrt{3}a/4=4.44$ bohr）、[110]方向は $a/\sqrt{2}=7.26$ bohr、[100]方向は $a=10.26$ bohr です。$\Delta B$ の構造がこの距離スケールに分布することから、Comptonプロファイルの方向異方性が結合の方向性に由来することが分かります。なお $B(r)$ はk点メッシュ由来の高周波エイリアシング（後述）の影響を受けにくく、生のヒストグラムから直接計算できます。

k点メッシュとエイリアシング：実験分解能との対応

Comptonプロファイルの計算では、離散的なk点から得られる運動量 $\mathbf{p}=\mathbf{k}+\mathbf{G}$ を $p_z$ のヒストグラムとして積算します。この過程には、2つの離散的なスケールが存在します。1つは射影方向に沿った運動量サンプリング間隔 $\delta \simeq 2\pi/(aN_k)$（Siの16³メッシュで0.038 a.u.）であり、もう1つは出力ヒストグラムのビン幅 $\Delta = 2\,\text{pzmax}/(\text{npz}-1)$（本計算では0.05 a.u.）です。生のプロファイルに現れる櫛状の構造（図4a）は、これら2つの離散スケールによるうなり（ビート）に起因するため、$\delta$ そのものの周期とは一致しません。したがって、これは物理的な構造ではなく、エイリアシングと呼ばれる離散化の産物です。

図4. Siにおけるk点メッシュ依存性（$p_z\parallel[100]$）。(a) 生ヒストグラム。櫛状構造の周期・振幅はメッシュ密度に対して単調には変化せず、$\delta$ とビン幅 $\Delta$ の関係によって決まります（実測の主周期は 8³: 0.144、16³: 0.163、20³: 0.136 a.u.。12³ の場合は $\delta=0.051$ が $\Delta=0.05$ とほぼ一致するためうなりが生じず、振幅が一桁小さくなります）。(b) Gauss広がり（FWHM 0.15 a.u.）適用後。8³〜20³ のすべての結果が $J(0)$ の1 %以内の差で一致します。

計算実務における効果的な対応策は、エイリアシングの周期（本計算では0.15 a.u.程度）と同等以上のFWHMを持つGauss関数で畳み込み（ブロードニング）を行うことです。図4(b)に示す通り、この処理を適用することで、8³から20³の異なるk点メッシュから得られたプロファイルが、$J(0)$ において1 %以内の差異（20³を基準として 8³で 1.0 %、12³で 0.3 %、16³で 0.5 %）に良好に収束します。さらに、この広がり幅は実際の放射光Compton散乱実験における運動量分解能（0.1〜0.4 a.u.）と同程度であるため、計算結果を実験値と直接比較する処理としても妥当性を持ちます。（Advance/PHASEでは broadening_fwhm を指定することで、広がりを適用したデータが自動的に出力されます。）結論として、生のヒストグラムに現れる微細構造を物理的に有意なものとして解釈することは避けるべきです。

2. Li：自由電子的なプロファイルとFermi運動量

Liは最も単純な金属であり、価電子1個がほぼ自由電子的に振る舞います。図5(a)に示すように、計算されたプロファイルは自由電子ガス（FEG）モデル $J(p_z)=\frac{3}{4p_\text{F}}\left[1-(p_z/p_\text{F})^2\right]$ と非常によく一致します。

図5. bcc Liの価電子Comptonプロファイル（48×48×48 メッシュ）。(a) 3方向の $J(p_z)$ と自由電子ガスモデル（破線、同じ広がりを与えたもの）。細線は生ヒストグラムで、k点メッシュ由来のエイリアシングを示しています。(b) 方向間の差。構造は $p_z \lesssim p_\text{F}=0.591$ a.u. の領域に広がっており、Fermi面の非球形性と高次運動量成分の方向依存性を反映しています。振幅は $J(0)$ の2〜4 %程度と小さく、k点メッシュに対する収束は遅くなります。

表2. Liの $J(0)$ と自由電子ガスとの比較（同じ広がりを与えた値、FWHM 0.15 a.u.）

方向	[100]	[110]	[111]	自由電子ガス
\(J(0)\) (a.u.)	1.222	1.201	1.254	1.255

Fermi運動量は電子密度から $p_\text{F}=(3\pi^2 n)^{1/3}=0.591$ a.u. と求まります。図5(b)の異方性は $p_z \lesssim p_\text{F}$ の領域全体に広がっており（$p_z=0$ と $p_z\simeq0.3$ a.u. に極値を持ちます）、$p_\text{F}$ の近傍だけに局在しているわけではありません。$J(p_z)$ が $p_z$ 一定面上の運動量密度をすべて足し上げた量である以上、これは自然な結果です。異方性の起源は、Fermi面がbcc Brillouin帯境界（$N$点）方向にわずかに膨らんでいること（非球形性）と、波動関数の高次運動量成分（Umklapp成分）の方向依存性の両方にあります。高分解能Compton散乱はこの非球形性を実際に検出できる感度を持ち、Sakuraiらの測定によって定量的に確認されています [5]。

ただし、Liの異方性は $J(0)$ の2〜4 %程度（差曲線の最大値で $J_{[110]}-J_{[100]}$ 2.4 %、$J_{[111]}-J_{[100]}$ 2.6 %、$J_{[111]}-J_{[110]}$ 4.4 %）しかない微小量であり、k点メッシュに対する収束が遅い点には注意が必要です。同じ広がり幅（FWHM 0.15 a.u.）に揃えて比較しても、32³ → 48³ で異方性の振幅は方向対によって1〜3割変化します。Si（12³〜20³ で1 %以内）とは対照的です。実験・計算の双方でLiの異方性測定が難しいとされる理由がここにあります。

なお、Liについては $p_\text{F}$ におけるプロファイルの折れ曲がりの大きさ（Fermi面のbreak、$Z_\text{F}$）に、理論と実験の顕著な食い違いが知られています。電子相関を取り込んだ電子ガス理論が0.5〜0.75程度を与えるのに対し（相関を含まない独立粒子近似では占有数が0または1となるため、$Z_\text{F}\simeq1$ です）、高分解能Compton散乱の実験からは著しく小さい値が報告されています [6]。これは電子相関効果の直接的な現れと解釈されています。DFT計算のプロファイルを基準線として実験との差分を取ることで、相関効果そのものを抽出するという解析が可能になるわけで、第一原理計算がこの分野で果たす役割の典型例といえます。

3. Fe：磁気Comptonプロファイルとスピン磁気モーメント

強磁性体に円偏光X線を入射すると、スピン上向き・下向きの寄与の差

$$J_{\rm mag}(p_z) = J_{\uparrow}(p_z) - J_{\downarrow}(p_z)$$

だけを取り出すことができます。これが磁気Comptonプロファイル（MCP）です。MCPの最大の特長は、軌道磁気モーメントを含まず、スピン磁気モーメントのみに感度を持つことです。したがってその積分値は

$$\displaystyle \int J_{\rm mag}(p_z)\, dp_z = \mu_{\rm spin}$$

となり、磁化測定（スピン＋軌道）やXMCD（総和則に仮定が必要）とは独立に、スピンモーメントを決定できます。

bcc Feの磁気Comptonプロファイルを図6に示します。

図6. bcc Feのスピン分解Comptonプロファイル。(a) $p_z\parallel[100]$ におけるスピン上向き・下向きのプロファイルと総和。(b) 磁気Comptonプロファイル $J_\uparrow-J_\downarrow$。[110]・[111] では $p_z=0$ でへこみ $p_z\simeq1.3$ a.u. にピークを持ちますが、[100] では $p_z=0$ が最大で $p_z\simeq0.7$ a.u. に浅い谷が現れます。この方向依存性そのものが、局在した3d電子と広がった4s電子のスピン分極の重なり方が方向によって異なることを示しています。細線は生ヒストグラム。

表3. bcc Feのスピン磁気モーメント

求め方	値 (\(\mu_\text{B}\))
SCF（実空間のスピン密度の積分）	2.182
磁気Comptonプロファイルの積分 \(\int (J_\uparrow-J_\downarrow) dp_z\)	2.169
実験の目安：飽和磁化 2.22 \(\mu_\text{B}\) から軌道成分（約0.1 \(\mu_\text{B}\)）を除いたスピン成分	約2.1

表3補足. 磁気Comptonプロファイルの方向依存性（広がり後）

方向	\(J_\text{mag}(0)\)	ピーク位置	ピーク値	\(\int J_\text{mag}\,dp_z\) (\(\mu_\text{B}\))
[100]	0.396	\(p_z=0\)	0.396	2.169
[110]	0.300	1.32 a.u.	0.403	2.171
[111]	0.238	1.28 a.u.	0.443	2.171

表3に示すように、bcc Feのスピンモーメント（積分値）は方向によらず 2.169〜2.171 $\mu_\text{B}$ で一致しており、方向ごとに独立に計算された3つのプロファイルが同じ全スピンに帰着することも、実装の整合性を裏づけています。

運動量空間で積分して得られたスピンモーメント（2.169 $\mu_\text{B}$）は、実空間のスピン密度から得られる自己無撞着な値（2.182 $\mu_\text{B}$）と0.6 %以内の差で一致しました。この一致は、PAW全電子再構成が運動量空間で正しく規格化されていること（表2）の直接的な帰結であり、Comptonプロファイル計算の内部整合性を確認する良い指標になります。実験のスピン成分（約2.1 $\mu_\text{B}$）との差はGGA-PBEが遷移金属の磁気モーメントをやや過大評価する傾向によるもので、汎関数に起因する既知の系統誤差の範囲内です。

MCPの形状にも重要な情報があります。磁気Comptonプロファイルは、正にスピン分極した局在3d電子（実空間で局在＝運動量空間で広がる）と、負にスピン分極した遍歴4s電子（運動量空間では $p_z$ の小さい領域に集中）の重ね合わせです。本計算では [110]・[111] で $p_z\simeq1.3$ a.u. のピークと $p_z=0$ のへこみという典型的な3d的形状が現れる一方、[100] では $p_z=0$ が最大となり、$p_z\simeq0.7$ a.u. に浅い谷が生じました（表3補足に数値を示します）。この方向依存性は、3d軌道のスピン分極が運動量空間で異方的に分布していることの直接的な現れです。MCPの形状を実験と比較することで、スピンモーメントを軌道ごとに分解できる点がMCPの威力です。実際に、電池材料であるリチウムマンガンスピネルなどにおいて、スピン（磁性）を持つ遷移金属軌道の同定にMCPが応用できることが示唆されています [7]。

まとめ

第一原理計算ソフトウェアAdvance/PHASEのCompton profile計算機能（PAW全電子再構成）を用いた解析により、性質の異なる3つの物質の電子状態を運動量空間から定量的に評価しました。Siでは共有結合の方向性に由来する異方性、Liではフェルミ面の微小な非球形性、そしてFeでは局在3d電子と遍歴4s電子の寄与が反映された純粋なスピン磁気モーメントを正確に抽出できることを確認しました。また、離散的なk点サンプリングに伴うエイリアシング（非物理的な振動）は、実験装置の分解能と同等のGauss広がりを与えることで適切に処理できるという実務上重要な知見も示しました。バルク敏感かつ非破壊であるCompton散乱実験と第一原理計算を組み合わせることで、化学結合やスピン状態の深い理解が可能となり、蓄電池のin operando観測をはじめとする最先端のデバイス・材料解析に対する強力なアプローチとなります。

本解析の詳細や、研究への適用可能性に関するご相談はこちら

お問い合わせ

参考文献

1. M. J. Cooper, "Compton scattering and electron momentum determination", Rep. Prog. Phys. 48, 415 (1985).
2. I. Makkonen, M. Hakala, and M. J. Puska, "Calculation of valence electron momentum densities using the projector augmented-wave method", J. Phys. Chem. Solids 66, 1128 (2005).
3. Th. Tschentscher, J. E. McCarthy, V. Honkimäki, and P. Suortti, "High-energy magnetic Compton scattering experiments at ESRF", J. Synchrotron Rad. 5, 940 (1998).
4. W. A. Reed and P. Eisenberger, "Gamma-ray Compton profiles of diamond, silicon, and germanium", Phys. Rev. B 6, 4596 (1972).
5. Y. Sakurai, Y. Tanaka, A. Bansil, S. Kaprzyk, A. T. Stewart, Y. Nagashima, T. Hyodo, S. Nanao, H. Kawata, and N. Shiotani, "High-resolution Compton scattering study of Li: Asphericity of the Fermi surface and electron correlation effects", Phys. Rev. Lett. 74, 2252 (1995).
6. W. Schülke, G. Stutz, F. Wohlert, and A. Kaprolat, "Electron momentum-space densities of Li metal: A high-resolution Compton-scattering study", Phys. Rev. B 54, 14381 (1996).
7. H. Hafiz, K. Suzuki, B. Barbiellini, N. Tsuji, N. Yabuuchi, K. Yamamoto, Y. Orikasa, Y. Uchimoto, Y. Sakurai, H. Sakurai, A. Bansil, and V. Viswanathan, "Tomographic reconstruction of oxygen orbitals in lithium-rich battery materials", Nature 594, 213 (2021).

関連ページ

• 第一原理計算ソフトウェア Advance/PHASE
• 解析分野：ナノ・バイオ
• 産業分野：材料・化学