X線Compton散乱による電池正極材料の第一原理解析：Li濃度非破壊定量とredox軌道同定#

PHASE logo

稼働中のリチウムイオン電池の内部でLiがどこにどれだけ存在するかを知ることは、性能・劣化・安全性のいずれにとっても本質的な情報です。しかし多くの分析手法はプローブの侵入深さが足りず、専用のテストセルでしか使えません。これに対し、100 keV級の高エネルギーX線を用いたCompton散乱分光は、密閉された市販電池をそのまま貫通してバルクの電子運動量分布を測ることができ、非破壊・in-operandoでのLi濃度計測を可能にします。その鍵となるのが、Compton profile $J(p_z)$ の形状パラメータ $S$ がLi濃度 $x$ に比例するという性質であり、実験値からLi濃度を読み取るにはあらかじめ第一原理計算で較正線を引いておくことが重要です。本解析では、第一原理計算ソフトウェアAdvance/PHASE (ver.5.0)に実装されたCompton profile計算機能を用いて、代表的な正極材料であるLi_xCoO₂について脱リチウム過程（ $x$ = 1, 2/3, 1/3, 0）のCompton profileを計算しました。得られた $S$ パラメータはLi濃度に対して高い直線性を示し、電池診断のための較正線を第一原理計算のみから構築できることを示します。さらに、Compton profileの差分 $\Delta J(p_z)$ を解析することで、充放電で出入りする電子がどの軌道に属するか（redox軌道）を運動量空間で直接可視化しました。

Keywords: 第一原理計算 (DFT), リチウムイオン電池, 正極材料, LiCoO₂, X線Compton散乱, 電子運動量密度, Compton profile, Sパラメータ, redox軌道, in-operando計測

Compton profileとは何を測っているのか#

入射X線のエネルギーが電子の束縛エネルギーより十分に大きいとき（衝撃近似）、非弾性散乱スペクトルの形は散乱ベクトル方向の電子運動量分布そのものになります。これをCompton profileと呼び、電子運動量密度 $\rho(\mathbf{p})$ の2重積分として定義されます [1]。

$J(p_z) = \iint \rho(\mathbf{p})\, dp_x\, dp_y , \qquad \rho(\mathbf{p}) = \sum_{n\mathbf{k}} f_{n\mathbf{k}} \left| \int \Psi_{n\mathbf{k}}(\mathbf{r})\, e^{-i\mathbf{p}\cdot\mathbf{r}} d\mathbf{r} \right|^2$

$J(p_z)$ は基底状態の一体波動関数を運動量表示で直接見る量であり、光電子分光や陽電子消滅のような行列要素を考慮する煩雑さがありません。バルク敏感で、温度・磁場・圧力といった試料環境も選ばず、金属容器の中の試料さえ測れます。この「素直さ」が、第一原理計算と実験を突き合わせる上での大きな利点です。

なお、運動量分布は原子核近傍の急峻な波束の構造に強く依存するため、通常の擬ポテンシャルを用いた計算では高運動量成分が欠落してしまいます。Advance/PHASEでは、この問題を解決する「PAW法に基づく全電子再構成」を適用しており、総和則 $\int J(p_z)\,dp_z = N$ を満たす高精度なプロファイルを算出しています [2]。

計算モデルと計算条件#

層状岩塩型のLi_xCoO₂（空間群 R $\bar{3}$ m、O3型積層）の六方晶単位胞（3化学式単位、Li₃Co₃O₆）を用いました。図1に示すように、CoO₂八面体シートとLi層が $c$ 軸方向に交互に積層した構造です。脱リチウムは、この単位胞に含まれる3枚のLi層を1枚ずつ取り除く「ステージ型」の規則配列でモデル化し、 $x$ = 1, 2/3, 1/3, 0 の4組成を作成しました。

この配列を採用した理由は単に最小のモデルだからではありません。全組成で格子ベクトルを固定した同一の六方晶セルを用いることで計算上は三方晶の点群操作が保たれるため、4組成すべてでk点メッシュ・対称性による全ブリルアンゾーン展開・運動量空間の分解能が完全に同一になります。

LixCoO2のモデル構造
図1. Li_xCoO₂（R $\bar{3}$ m, O3型）の計算モデル。 $c$ 軸を縦にとった側面図です。緑=Li、青=Co、赤=O、破線の円はLi空孔を表します。格子定数はLiCoO₂の実験値（ $a$ = 2.8161 Å, $c$ = 14.0536 Å）に固定し、酸素の内部座標のみを各組成で構造最適化しています。

格子定数を全組成で共通に固定したのは意図的な設定です。こうすることで、Compton profileの変化が格子変形ではなくLi脱離に伴う電子状態の変化（酸化還元）だけを反映するようになり、 $S$ パラメータとLi濃度の関係を純粋に電子的な効果として取り出せます。

表1. 計算条件の概要

項目	設定
擬ポテンシャル	PAW
交換相関汎関数	GGA-PBE
カットオフエネルギー（波動関数 / 電荷密度）	30 Ry / 270 Ry
k点メッシュ（構造最適化 / Compton profile）	6×6×2 / 12×12×12
単位胞	六方晶 3化学式単位
Compton profile	オンサイト異方性とテール補正: 有効
$p_z$ メッシュ	−10 〜 +10 a.u.、1001点（刻み 0.02 a.u.）
投影方向	[001]（層に垂直）、面内、マジック角（※）

※ 実験は多結晶（粉末）試料に対して行われるため、比較すべき理論量は球面平均したCompton profileです。本解析では、2次のルジャンドル多項式 $P_2(\cos\theta) = 0$ を満たす、すなわち $\cos^2\theta = 1/3$ となる「マジック角」方向を第3の投影軸に選びました。この方向では異方性の主要項（ $l=2$ ）がちょうど消えるため、単一方向の計算で粉末平均を良い精度で代表できます。

計算結果と考察#

数値精度：総和則による検証#

Compton profileの計算がどれだけ信頼できるかは、まず総和則 $\int J(p_z)\,dp_z = N$ （ $N$ は価電子数）で確認できます。LiCoO₂（ $N$ = 66）の場合、平面波部分とPAWオンサイト項の厳密な内訳の和は 65.987 電子となり、価電子数を 0.02 % の差で再現しました。出力窓 $|p_z| \le$ 10 a.u. 内での積分値は 65.908 電子で、残りは全電子波動関数が本来もつ 10 a.u. を超える高運動量テールです（これは誤差ではなく物理です）。

このように総和則が高い精度で満たされることは、平面波部分とオンサイト補正項が整合して評価できていることを意味します。総和則は形式的には自明ではなく、擬波動関数のノルムとPAW補強項の両方を正しく足し上げて初めて成立します。実際、同じ系を 4×4×1 の粗いk点メッシュで計算すると2 %程度の欠損が残り、密なメッシュで0.1 %を切ります。総和則は、Compton profile計算の収束を確認するための実用的な指標として使えます。

Li濃度に対するCompton profileの変化#

図2に、4組成の価電子Compton profileを示します。Liを挿入するほど $p_z \approx 0$ 近傍のピークが高くなります。しかもその増加分は $|p_z| \lesssim$ 1 a.u. の狭い領域に集中しており、 $|p_z| \gtrsim$ 2 a.u. の裾はほとんど変化しません。これは、Liの2s電子が正極ホスト（O 2p / Co 3d混成帯）に供与されると、その電子が実空間で広がった（＝運動量空間では狭い）状態を占めるためです。実験で観測される「Li濃度が高いほどピークが高い」という傾向 [3] を、そのまま再現しています。

LixCoO2の価電子Compton profile
図2. Li_xCoO₂ の価電子Compton profile（粉末平均相当、実験分解能0.1 a.u.相当のGauss広がりを付与）。右は低運動量領域の拡大図です。曲線はいずれも各組成の価電子数（ $x$ に応じて63〜66電子）に規格化されており、縦軸は単位胞あたりの値です。

Sパラメータ較正線 ― 本解析の主結果#

Compton profileの形状変化を1つの数値に凝縮するために、Suzukiらの定義 [3] に従って形状パラメータ $S$ を導入します。

$S = \dfrac{S_L}{S_H}, \qquad S_L = \int_{-d}^{d} J(p_z)\,dp_z, \qquad S_H = \int_{-10}^{-d} J(p_z)\,dp_z + \int_{d}^{10} J(p_z)\,dp_z$

実験のCompton profileには内殻電子の寄与が含まれます。一方、PAW計算が与えるのは価電子profileですので、比較のために自由原子の内殻profile（Li 1s、O 1s、Co [Ar]殻）を全電子原子計算から求めて加算しました。内殻はLi脱離では変化しないため、 $S$ の $x$ 依存性はすべて価電子部分から来ます。低運動量領域の境界は文献 [3] と同じ $d$ = 6 a.u. としています。

SパラメータのLi濃度依存性
図3. 全profile（価電子＋内殻）から求めた $S$ パラメータのLi濃度依存性と最小二乗直線。3つの投影方向の結果がほぼ重なっており、 $S$ が投影方向にほとんど依存しない（＝粉末試料でも使える）指標であることが分かります。

図3のとおり、 $S$ はLi濃度 $x$ に対して高い直線性を示し、最小二乗フィットから

$S = 1.1042\,x + 14.7224$

が得られました。決定係数は $R^2 > 0.999$ と非常に高い直線性を示し、傾き 1.1042 は、Li濃度が 0.1 変化すると $S$ が 0.1104 変化することを意味します。 $p_z \approx 0$ のピーク高さ自体は局所的な非線形性を示しますが、広い運動量窓で積分して実質的に電子数を数える $S$ パラメータを用いることで、その非線形性が希釈されて良好な直線関係が得られます。実験側の $S$ の統計誤差は積算光子数 $I_S$ に対して $\Delta S \propto I_S^{-1/2}$ でスケールしますので [3]、到達可能なLi濃度分解能はこの傾きと光子統計だけで決まります。傾きが大きいほど少ない積算で高い濃度分解能が得られることになります。

興味深いことに、得られた切片（14.72）はスピネルLi_xMn₂O₄について実験的に報告されている較正線 $S = 0.59x + 15.3$ [3] の切片と近い値を示しました。一方で、傾きに関してはLi_xCoO₂（1.10）がLi_xMn₂O₄（0.59）の約2倍となっています。これは、「Li挿入に伴い低運動量成分が増加する」という根本的な物理メカニズムは共通しているものの、その感度（傾き）は遷移金属の種類や結晶構造に強く依存することを示しています。このように材料ごとに異なる応答を示すからこそ、未知の新規バッテリー材料に対して、第一原理計算を用いてあらかじめ専用の較正線を高精度に引けることが、本手法の実用上の強みとなります。

各組成の計算結果は表2にまとめています。

表2. 各組成の計算結果まとめ

$x$	組成（単位胞）	価電子数	$S$（全profile, $d$=6 a.u.）	$\int J\,dp_z / N$ (%)
1	Li₃Co₃O₆	66	15.8265	99.87
2/3	Li₂Co₃O₆	65	15.4588	99.86
1/3	Li₁Co₃O₆	64	15.0906	99.86
0	Co₃O₆	63	14.7223	99.86

Redox軌道の同定 ― どの軌道の電子が出入りしているのか#

Compton profileの真価は、単なる濃度計であることを超えて電子構造そのものを読み出せる点にあります。異なるLi濃度のprofileの差 $\Delta J(p_z)$ をとると、内殻や充放電に関与しない価電子の寄与が消え、実際に出入りした電子の運動量分布だけが残ります [4]。

Compton profileの差分によるredox軌道の同定
図4. 左：移動電子1個あたりに規格化した差分Compton profile $\Delta J(p_z)$ 。参考として自由原子のO 2p、Co 3d、Li 2s軌道のprofile（いずれも1電子に規格化）を破線で示します。右： $\Delta J$ の累積積分。 $|p_z|$ をどこまで積分すると移動電子数（3個）に達するかを示しています。

3つの組成差（ $x$ = 1 → 0、2/3 → 0、1/3 → 0）は、移動電子1個あたりのピーク値が 0.47 / 0.48 / 0.48 a.u.⁻¹ と2 %以内の差で一致し、形も重なります。脱リチウム過程の全域を通じて同じ軌道が一貫して酸化還元に関与していることを示す結果です。その形は自由原子のCo 3d軌道よりも明らかに狭く、自由原子のO 2p軌道よりもさらにやや狭い分布になっています。運動量分布の広がりは実空間の局在度と逆数関係にあるため、これは「充放電で出入りする電子は、Coサイトに単独で強く局在した3d電子というより、強い共有結合によってO 2p成分を多く含む広がった混成軌道に属する」ことを意味します。この結果は、Li_xCoO₂の軟X線吸収分光で指摘されているO 2p性ホールの重要性 [5] や、スピネルLi_xMn₂O₄における共有結合性のO 2pレドックスの報告 [4] とよく整合しています。

なお、 $\Delta J$ は $|p_z|$ ≈ 1.6 a.u. でわずかに負に転じ、 $|p_z|$ ≈ 1.8 a.u. で極小となります。3組成すべてで符号・位置ともに再現し（ピーク値の 3 – 5 %、実験分解能0.1 a.u.相当の広がりを付与した場合）、積分すると移動電子1個あたり −0.005 個（ $p_z>0$ 側）に相当します。Li挿入によってこの運動量域の重みがわずかに減る、すなわち一部の電子が運動量空間で狭く（実空間で広く）なることを意味します。

図4右の累積積分は、移動した3電子のうち 95 % が $|p_z| \lt$ 1.5 a.u. の低運動量領域に収まっていることを示します。 $S$ パラメータが $d$ = 6 a.u. という広い窓を使いながらLi濃度に鋭敏に反応するのは、まさにこの「移動電子が低運動量側に集中する」性質によるものです。

層状構造に由来する異方性#

単結晶試料が得られる場合には、方向別のCompton profileから結合の異方性を直接読み取ることができます。図5は層に垂直な[001]方向と面内方向のprofileの比較です。

Compton profileの方向異方性
図5. $x$ = 1（LiCoO₂）と $x$ = 0（CoO₂）の方向別価電子Compton profile（左）と、[001]方向と面内方向の差（右、実線）。右図の点線は差を $p_z$ = 0 から積分して電子数に換算した累積量（右軸、 $\pm p_z$ 両側を数えた値）で、総和則によりいずれも大きな $p_z$ でゼロに戻ります。k点メッシュに起因する離散化を平滑化するため、図5では0.25 a.u.のGauss広がりを用いています。

LiCoO₂では層に垂直な[001]方向のprofileが面内方向より鋭いピークを持ちますが、裾野は逆にわずかに厚くなります。これは、Liが供与する層間の非局在電子（ $p_z$ 方向に狭い成分）と、CoO₂層内への2次元的閉じ込め（不確定性原理により $p_z$ 方向に広い成分）という2つの効果の重ね合わせです。脱リチウム後のCoO₂ではLi由来の鋭いピーク成分が大きく減少する一方、CoO₂層への閉じ込め効果による高運動量側の振動構造はむしろ大きくなります。

単結晶試料が用意できる場合には、このように方向別Compton profileから充放電に伴う結合性の変化 ― 層間の電子がどれだけ失われ、層内の閉じ込めがどう際立つか ― を直接読み取ることができます。

まとめ#

第一原理計算ソフトウェアAdvance/PHASEのCompton profile計算機能（PAW全電子再構成）を用いて層状正極材料Li_xCoO₂の脱リチウム過程を解析し、総和則を0.02 %の誤差で満たす高精度なプロファイルを得ました。算出された形状パラメータ $S$ はLi濃度に対して良好な直線関係を示すため、密閉された実電池のin-operando Compton散乱測定からLi濃度を定量するための「較正線」を第一原理計算のみから構築できることが実証されました。さらに、充放電に伴うプロファイルの差分 $\Delta J(p_z)$ を解析することで、実際に酸化還元に寄与するredox軌道がCo 3dよりもO 2p性の強い広がった状態であることを運動量空間で直接同定しました。本手法は、組成やドーピングが異なる多様な正極材料に対しても同様に適用可能であり、バルク透過性に優れたCompton散乱実験と第一原理計算を連携させることで、リチウムイオン電池の非破壊診断技術と新規材料設計の双方を促進するアプローチとなります。

※ 本解析はGGA-PBEの非スピン分極計算に基づいています。低スピンCo³⁺/Co⁴⁺を含むLi_xCoO₂ではこの扱いが妥当ですが、Mn系のように大きな磁気モーメントをもつ正極材料ではスピン分極を考慮する必要があります。本機能はスピン分解出力（磁気Compton profile）にも対応しています。また、Li空孔配列はステージ型の規則配列を仮定しています。 $S$ パラメータはブリルアンゾーン積分量であり空孔配列の詳細には鈍感ですが、単結晶の方向依存性を詳細に議論する場合には配列の影響を別途検討する必要があります。

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

お問い合わせ

参考文献#

M. J. Cooper, "Compton scattering and electron momentum determination", Rep. Prog. Phys. 48, 415 (1985).
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).
K. Suzuki, B. Barbiellini, Y. Orikasa, S. Kaprzyk, M. Itou, K. Yamamoto, Yung Jui Wang, H. Hafiz, Y. Uchimoto, A. Bansil, Y. Sakurai, and H. Sakurai, "Non-destructive measurement of in-operando lithium concentration in batteries via x-ray Compton scattering", J. Appl. Phys. 119, 025103 (2016).
K. Suzuki, B. Barbiellini, Y. Orikasa, N. Go, H. Sakurai, S. Kaprzyk, M. Itou, K. Yamamoto, Y. Uchimoto, Yung Jui Wang, H. Hafiz, A. Bansil, and Y. Sakurai, "Extracting the Redox Orbitals in Li Battery Materials with High-Resolution X-Ray Compton Scattering Spectroscopy", Phys. Rev. Lett. 114, 087401 (2015).
T. Mizokawa, Y. Wakisaka, T. Sudayama, C. Iwai, K. Miyoshi, J. Takeuchi, H. Wadati, D. G. Hawthorn, T. Z. Regier, and G. A. Sawatzky, "Role of Oxygen Holes in Li_xCoO₂ Revealed by Soft X-Ray Spectroscopy", Phys. Rev. Lett. 111, 056404 (2013).