Skip to content

完全磁気多極子基底に基づくX線円二色性・線二色性の総和則

X線磁気円二色性(XMCD)とX線磁気線二色性(XMLD)の総和則は、吸収スペクトルをエネルギーで積分した量を、基底状態における角運動量・スピン・電荷(軌道占有)の異方性に結び付ける理論である。完全磁気多極子基底(complete magnetic multipole basis)を用いると、従来は混在して見えにくかったXAS/XMCD/XMLD の積分量を「スピンレス(s=0)」と「スピンフル(s=1)」の多極子として同一の代数構造で整理され、反強磁性体や交替磁性におけるXMCD・XMLDの起源を同一の代数構造で扱えるようになる。

alt text

参考ドキュメント

  1. Y. Yamasaki, Y. Ishii, N. Sasabe, Sci. Technol. Adv. Mater. 26, 2513217 (2025). doi:10.1080/14686996.2025.2513217
  2. B. T. Thole et al., Phys. Rev. Lett. 68, 1943 (1992). doi:10.1103/PhysRevLett.68.1943
  3. J. Okabayashi, PF NEWS Vol.38 No.4 (2021) 磁気線二色性を用いた垂直磁気異方性の研究

0. 記号

  • 電磁波

    • 入射光子エネルギー ω
    • 偏光ベクトル ϵ
    • 右回り・左回り円偏光を ϵ=+, と書く
    • 線偏光の方向を e と書く(例:ez
  • 電子状態

    • 初期状態(基底状態) |g
    • 終状態(励起状態) |f
    • 電子数演算子や角運動量は基底状態平均 g||g を用いる
  • 角運動量

    • 軌道角運動量 L、スピン角運動量 S
    • 一電子の軌道量子数を l、全角運動量を j とする
    • L2,L3 端ではコア軌道が 2p であり lc=1jc=1/2,3/2 である
  • スペクトル積分

    • 吸収係数 μϵ(ω)
    • エッジ領域でのエネルギー積分を edgedω と書く

1. 研究の狙い

従来の Thole–Carra 総和則は、L2,3 の 円二色性の積分から軌道モーメント Lz と、スピンモーメントに磁気双極子項(anisotropic magnetic dipole)を加えた組 Sz+7Tz を抽出できることを示した。一方で近年、反強磁性体やアルターマグネットなど「巨視的磁化がゼロでも XMCD/XMLD が出る」系が増え、信号の起源を“磁化そのもの”から切り離して再分類する必要が生じた。

Y. Yamasaki ら(2025)は、吸収過程を スピン非反転(s=0) と スピン反転(s=1) に分解し、さらに j±(内殻スピン軌道分裂した2本のエッジ)を 和/差で組み替えることで、スピンレス多極子とスピンフル多極子を選択的に取り出す枠組みを与えた。 この整理により、XMCD/XMLD の積分が「どの(スピンレス/スピンフル)多極子演算子の期待値」を見ているかが明示化される。

2. X線吸収の出発点(E1遷移)

本章では、総和則の出発点である「吸収係数の定義」と「エッジ積分による演算子積への帰着」を、式の意味が追えるように段階的に記述する。ここで扱うのは電気双極子(E1)近似であり、空間的には一次(r に一次)の相互作用のみを残す。

2.1 フェルミの黄金律

入射X線が物質に吸収される過程は、基底状態 |g から励起状態 |f への遷移として表される。時間依存摂動論の一次(線形応答)では、単位時間あたりの遷移確率 Wgf はフェルミの黄金律に従う:

Wgf(ω)=2π|f|H^(ω)|g|2δ(EfEgω).

ここで

  • Eg,Ef はそれぞれ |g,|f の固有エネルギーである。
  • δ(EfEgω) はエネルギー保存を表し、「光子がもつエネルギー ω を吸収したときに到達できる終状態のみ」が寄与することを意味する。

吸収係数 μϵ(ω) は、光の偏光状態 ϵ のもとでの吸収確率(あるいは吸収断面積)に比例する量として定義されるので、状態和をとって

μϵ(ω)f|f|T^ϵ|g|2δ(EfEgω)

と書ける。ここで T^ϵ は光と電子の相互作用(遷移)を与える演算子である。

E1(電気双極子)近似の遷移演算子

電磁場と電子の結合は一般に最小結合(pp+eA)から導かれるが、X線吸収の総和則では、空間依存 eikr を長波長展開し、最初の非自明項(一次)として電気双極子項を取る。すると遷移演算子は

T^ϵ=ϵr

である。ここで

  • r は電子位置演算子であり、多電子系では和 r=iri と理解する。
  • ϵ は偏光ベクトルであり、円偏光・線偏光の違いは ϵ の成分の取り方として入る。

この形が重要なのは、T^ϵ がベクトル(階数1のテンソル)演算子であるため、後で角運動量代数により

11=012

という分解(等方成分・円二色性成分・線二色性成分)が自然に現れるからである。

2.2 スペクトル積分と閉包関係

総和則で扱うのは「あるエッジ領域(例:L2L3)にわたる吸収強度の積分」である。そこで

Iϵedgedωμϵ(ω)

を考える。μϵ(ω) の定義を代入して

Iϵedgedωf|f|T^ϵ|g|2δ(EfEgω).

ここで δ 関数の積分を明示する。一般に

dω δ(EfEgω)=1

であり、積分範囲がエッジ全体を十分に含むとき(寄与する |f を取りこぼさないとき)に、終状態ごとの寄与がそのまま足し合わさる。比例定数をまとめると

Iϵf|f|T^ϵ|g|2.

この絶対値二乗をブラケットで書き直すと

|f|T^ϵ|g|2=g|T^ϵ|ff|T^ϵ|g.

よって

Iϵfg|T^ϵ|ff|T^ϵ|g.

ここで閉包関係(完全系)

f|ff|=1^

を用いる。すると

Iϵg|T^ϵ(f|ff|)T^ϵ|g=g|T^ϵT^ϵ|g

が得られる。

この変形が意味すること

この式変形の物理的意味は次の通りである。

  • スペクトルのエネルギー分布(どの ω に強度があるか)を捨てて、エッジ全体で積分すると「終状態の詳細」への依存が消え、基底状態期待値だけで書ける量に落ちる。
  • その結果、実験で測れる積分強度 Iϵ は、基底状態の演算子平均 T^ϵT^ϵ へと帰着する。
  • 総和則は、T^ϵT^ϵ を角運動量代数で分解し、Lz,Sz や四極子期待値など「解釈可能な基底状態量」を右辺に出現させる手続きである。

次章へつながる形への書き換え

双極子演算子を球面テンソルで表すと(M=0,±1

T^ϵ=ϵr=M=11ϵMr^M,T^ϵ=M=11ϵMr^M.

したがって

T^ϵT^ϵ=M,MϵMϵMr^Mr^M,

よって積分強度は

IϵM,MϵMϵMr^Mr^M

となる。

ここで r^M は階数1テンソルであるから、演算子積 r^r^ は角運動量加算により

11=012

に分解される。すなわち、Iϵ の中には

  • K=0 成分:等方吸収(XAS)
  • K=1 成分:擬ベクトル成分(円二色性 XMCD)
  • K=2 成分:二階テンソル成分(線二色性 XMLD)

が必然的に現れる。以降の章では、この分解を 3j,6j,9j 記号を用いて厳密に行い、最終的に Lz,Sz、四極子、異方性磁気双極子項などへ結び付ける。

3. 球面テンソルによる偏光の取り扱い

3.1 球面基底

ベクトル演算子 r を球面成分 rqq=0,±1)で表す。

r0=z,r±1=12(x±iy)

同様に偏光ベクトルも球面成分 ϵq をもつ。すると

T^ϵ=ϵr=q=11ϵqrq

である。

3.2 円偏光・線偏光

  • 右回り円偏光(入射方向を +z とする)

    ϵ+: (ϵ+1,ϵ0,ϵ1)=(1,0,0)
  • 左回り円偏光

    ϵ: (ϵ+1,ϵ0,ϵ1)=(0,0,1)
  • 線偏光 ez

    (ϵ+1,ϵ0,ϵ1)=(0,1,0)
  • 線偏光 ex

    (ϵ+1,ϵ0,ϵ1)=(12,0,12)

これにより、XMCD(円偏光差)とXMLD(線偏光差)は偏光テンソルの組合せとして統一的に扱える。

4. 角運動量代数での積分強度の分解

4.1 積分強度

先の式より、積分強度は

Iϵedgedωμϵ(ω)T^ϵT^ϵ=q,qϵqϵqrqrq

である。

rq は階数1の球面テンソルであり、テンソル積 rr は角運動量加算で

11=012

に分解される。したがって

rqrq=K=02Q=KKC1q,1qKQ[rr]Q(K)

が成り立つ(C1q,1qKQ はクレブシュ=ゴルダン係数)。

4.2 偏光テンソル

偏光側も同様に

PQ(K)(ϵ)q,qϵqϵqC1q,1qKQ

と定義すると

IϵK,QPQ(K)(ϵ)[rr]Q(K)

となる。ここで

  • K=0 は等方XAS
  • K=1 は円二色性(擬ベクトル)
  • K=2 は線二色性(四極子)

に対応する。

5. 角運動量代数による PMMj± の展開

本章では、電気双極子(E1)遷移演算子 r^M を角運動量代数で整理し、コア準位の jj 結合(j±=lc±12)と価電子(ホール)側の基底を接続する。結論として、スペクトル積分に現れる一般係数 PMMj± は、3j6j9j 記号で表現でき、どの多極子成分がE1-XAS/XMCD/XMLDに現れるかが三角条件(選択則)として自動的に選別される。

5.1 コアと価電子の基底

コア(jj 結合)基底

コア準位はスピン軌道相互作用により j=lc±12 に分裂しているため、コア状態は

|lc,12;jmj,j{j+=lc+12, j=lc12}

で表す。

この jj 結合状態は、軌道とスピンの直積基底 |lcmc|12σσ=±12)にクレブシュ–ゴルダン係数で展開できる:

|lc,12;jmj=mc,σClcmc, 12σjmj|lcmc|12σ.

価電子(ホール)側基底

価電子(ホール)側は、まず軌道とスピンを分けた直積基底

|lvm,12σ|lvm|12σ

を用いる(必要なら後で jv を導入して再結合する)。

電気双極子遷移演算子

電気双極子(E1)演算子は軌道空間で階数1の球面テンソルであり、球面成分 M=0,±1 を用いて

r^M

と書く。スピン空間には単位演算子として作用するため、E1ではスピンは保存される:

12σ|1^|12σc=δσ,σc.

したがって、E1遷移の全行列要素は

lvm,12σ|r^M|lcmc,12σc=lvm|r^M|lcmcδσ,σc

である。

軌道部分のウィグナー–エッカルト分解(3j 記号)

軌道部分はウィグナー–エッカルト定理により

lvm|r^M|lcmc=(1)lvm(lv1lcmMmc)lvr^lc.

ここで lvr^lc は縮約行列要素であり、磁気量子数に依存しない。

コアの jj 結合基底への接続

求めたいのはコアの |lc,12;jmj から価電子 |lvm,12σ への行列要素である。コア側を直積基底に展開すると

lvm,12σ|r^M|lc,12;jmj=mc,σcClcmc, 12σcjmjlvm|r^M|lcmcδσ,σc.

δσ,σc により σc=σ が固定されるので

lvm,12σ|r^M|lc,12;jmj=mcClcmc, 12σjmj(1)lvm(lv1lcmMmc)lvr^lc.

この段階で、角運動量の幾何学は 3jClm,sσjm の積として現れている。

価電子側も jj 結合で扱う場合(6j 記号)

価電子側も jv を導入して |lv,12;jvmv を用いると、全行列要素は

nvlvjvmv|r^q|nclcjcmc=(1)jvmv(jv1jcmvqmc)nvlvjvr^nclcjc

となる。

この縮約行列要素は、スピン s=12 を用いて

nvlvjvr^nclcjc=(1)lv+12+jc+1(2jv+1)(2jc+1){lvjv12jclc1}nvlvr^nclc

6j 記号で表される。

さらに軌道縮約行列要素は、動径積分 nv|r|nc を用いて

nvlvr^nclc=nv|r|nc(1)lv(2lv+1)(2lc+1)(lv1lc000).

よって、エッジ全体の積分で動径積分が共通因子として括られると、残る係数は 3j6j によって完全に決まる。

5.2 3j 記号・6j 記号・9j 記号

3j 記号とクレブシュ–ゴルダン係数

クレブシュ–ゴルダン係数と 3j 記号の関係は

Cj1m1, j2m2JM=(1)j1j2+M2J+1(j1j2Jm1m2M).

3j 記号が非零となるための条件として

m1+m2M=0,|j1j2|Jj1+j2

が必要であり、これが選択則(角運動量の三角条件)を与える。

6j 記号のRacah公式

6j 記号は結合順序の変更(再結合)を表す係数であり、Racah公式を明示すると

{abcdef}=Δ(a,b,c)Δ(a,e,f)Δ(d,b,f)Δ(d,e,c)z(1)z(z+1)!(zabc)!(zaef)!(zdbf)!(zdec)!(a+b+d+ez)!(a+c+d+fz)!(b+c+e+fz)!.

ここで

Δ(x,y,z)=(x+yz)!(xy+z)!(x+y+z)!(x+y+z+1)!.

角運動量 a,b, は整数または半整数であり、階乗 x! は一般にはガンマ関数 x!Γ(x+1) として解釈するのが厳密である。実際に和 z に寄与するのは、分母の階乗の引数が非負となる z の範囲に限られるため、有限個の項のみが残る。

6j 記号が非零であるためには、(a,b,c)(a,e,f)(d,b,f)(d,e,c) の各三つ組がそれぞれ三角条件を満たす必要がある。

9j 記号と 6j 記号の和

9j 記号は三つの角運動量結合の再結合を表し、6j 記号の和として

{abcdefghi}=x(1)2x(2x+1){adghix}{behdxf}{cfixab}

と書ける(位相の取り方は定義に依存するため、使用する規約に合わせて統一する必要がある)。

PMMj±9j 記号が現れる理由

E1-XASのスペクトル積分では、双極子演算子の積

r^Mr^M

が現れる。r^M は階数1テンソルであるため、テンソル積は

11=012

に分解される。この分解を、コアの jj 結合(j±)と価電子側のスピン構造(12)を保ったまま行うと、角運動量の結合順序を複数回入れ替えることになる。結合順序の入れ替えは 6j 記号で与えられ、二段以上の再結合をまとめると 9j 記号に集約される。

したがって、一般係数 PMMj± の角部分は、最終的に 3j6j9j 記号の積(あるいは 9j を用いた簡潔な形)で表される。ここで、各記号に付随する三角条件が、E1遷移に現れる多極子成分の選別条件として働く。

6. XMCD総和則の導出

6.1 円二色性の定義

L2,3 端の吸収を μ±(ω) とし、円二色性(XMCD)を

Δμ(ω)=μ+(ω)μ(ω)

と定義する。

さらに、実験で用いられる積分量として

pL3dωΔμ(ω),qL2+L3dωΔμ(ω),rL2+L3dω(μ+(ω)+μ(ω))

を導入する。

6.2 K=1 成分の抽出

円偏光では偏光テンソル PQ(K) のうち K=1 が寄与し、Q=0 成分が z 軸方向の擬ベクトルに対応する。したがって

qO^0(1)

という形に帰着する。O^0(1) は価電子の角運動量・スピンに結び付く演算子である。

6.3 軌道総和則

Tholeらは、L2,3 の両端を合わせた円二色性積分が軌道角運動量 Lz に比例することを示した。定数因子を含めて、磁気モーメント mL

mL=μBLz

と定義すると、総和則は

mL=43μBnhqr

と書ける。ここで nhd 殻のホール数である。

この係数 4/3 は、2p3d における lc=1,lv=2 と、jc=1/2,3/2 の加重和から角運動量代数で決まる。

6.4 スピン総和則と磁気双極子項

Carraらは、L3L2+L3 の組合せでスピンに比例する和が得られることを示した。 citeturn0search2
スピン磁気モーメント mS

mS=2μBSz

と定義すると

mS+7mT=2μBnhp2qr

となる。ここで mT は磁気双極子項(anisotropic magnetic dipole)に対応する寄与であり、価電子のスピン密度分布の異方性を反映する。

6.5 tz の等価演算子

完全磁気多極子基底の枠組みでは、スピンフル磁気双極子の成分として tz が現れ、等価演算子として

tz=14(3[lz(ls)]+2l2sz)

の形で書かれることが示されている(記法は文献に従う)。 -L27

この項は、反強磁性体でのXMCDの起源(巨視的磁化がゼロでも局所的スピン密度が異方的であること)を記述する有力な指標であり、Mn3Sn、RuO2、MnTe などへの言及がある。

7. XMLD総和則の導出

7.1 線二色性の定義

線偏光方向を変えた吸収の差として

Δμlin(ω)=μex(ω)μez(ω)

などを定義する。紙面上は、面外と面内を区別して次の2種を考える。

  • 面外XMLD(zXMLD)

    S20j±(Ixj±+Iyj±2Izj±)

    のように、偏光が xy 面内か z 方向かの差をとる。

  • 面内XMLD(xyXMLD)

    S22j±(Ixj±Iyj±)

    のように、偏光が xy の差をとる。

7.2 K=2 成分の抽出

線偏光差は偏光テンソルの K=2 成分で記述され、価電子側の期待値は二階テンソル(四極子)に対応する。したがって

dωΔμlin(ω)Q^(2)

となる。

7.3 スピンレス四極子総和則

スピンレス多極子としての四極子は Q2m(0;0) であり、面外成分 m=0

O20=12(3z2r2)

に対応する。

等価演算子として

Qzz=12(lz213l2)

が与えられる。

したがって、面外XMLDの積分値は Qzz の期待値に比例し、d3z2r2dx2y2,dxy などの占有差(電荷分布異方性)を反映する。

7.4 スピンフル四極子総和則

スピンフル四極子は、磁気双極子項や磁気八極子項を含み、面外成分に対してスピン軌道相互作用 λ(ls) の異方性を示す成分が現れる。

等価演算子の一例として

Pzz=12(3lzszls)

が与えられる。

これは磁気結晶異方性エネルギーと関係する量であり、van der Laan による線二色性と異方性の関係付けとも整合する。

7.5 面内XMLDと面内異方性

面内XMLDでは、面内異方性を反映するスピンフル八極子由来の項が現れ、lxsxlysy のような面内スピン軌道相互作用の異方性を表す演算子が導入される。

さらに、NiO のネールベクトルの空間分布可視化に面内XMLDが用いられた例が言及されている。

8. 完全磁気多極子基底

完全磁気多極子基底とは、単一中心(原子サイト)上の電子自由度に対して、角運動量代数と対称操作(空間反転・時間反転)に基づいて、観測可能な演算子を体系的に分類する枠組みである。ここでの要点は、同じ階数 l の「多極子」に見えても、空間反転 I と時間反転 T に対する性質が異なる複数の族が存在し、それらを Q,M,T,G の4種として整理する点にある。

8.1 4種多極子

空間反転と時間反転の基本

空間反転(パリティ)Irr を与え、時間反転 TppLLSS を与える(磁性の起源となる角運動量は時間反転で符号が反転する)。

多極子の階数を l とすると、位置ベクトル(極性ベクトル)に由来する球面調和成分 Ylm(r^) は空間反転で (1)l の符号をもつ。一方、磁気的な量(電流や角運動量)を含む多極子は、パリティの性質が l だけでは決まらず、もう一段の符号規則をもつ。

4種の分類

完全磁気多極子基底では、多極子を以下の4種に分類する。

  • 電気多極子 Q(electric multipole)
  • 磁気多極子 M(magnetic multipole)
  • 磁気トロイダル多極子 T(magnetic toroidal multipole)
  • 電気トロイダル多極子 G(electric toroidal multipole)

これらは IT に対する変換性で区別される。階数 l の一般的な符号規則を

  • 空間反転の符号:I:(1)πl
  • 時間反転の符号:T:±

としてまとめると、基本的な整理は次のようになる。

種類空間反転 I時間反転 Tl=1 の代表例
Qlm(1)l+電気双極子(極性)
Mlm(1)l+1磁気双極子(磁化)
Tlm(1)lトロイダル双極子(アナポール)
Glm(1)l+1+電気トロイダル双極子

この表の意味は、例えば l=1 の場合に

  • 磁気双極子 M1m は空間反転で偶(+1)、時間反転で奇()である
  • 磁気トロイダル双極子 T1m は空間反転で奇(1)、時間反転で奇()である

といった区別ができる、ということである。反強磁性体やアルターマグネットでは、巨視的磁化(M1m)がゼロでも、時間反転で奇な別種の多極子成分が選択的に現れる場合があり、この4分類が解釈の骨格になる。

8.2 スピンレス多極子

スピンレス(s=0)とは、演算子がスピン空間に単位演算子として作用し、軌道自由度のみで構成されることを意味する。パウリ行列の球面テンソル表示で σ0 を単位行列とすると、スピンレス多極子演算子は

X^lm(0;0)=X^lm(orb)σ0

である。ここで X^lm(orb) は軌道部分の球面テンソル演算子であり、XQ,M,T,G のいずれかを表す。

軌道多極子の「殻内」表現

実際の総和則では、ある殻(例:d 殻)内の演算子が重要になる。殻の角運動量を l とし、その基底 |lm を用いると、軌道テンソル演算子は一般に

X^κq(orb)=m1,m2lm1|X^κq(orb)|lm2cm1cm2

の形に書ける。κ は軌道テンソルの階数、q はその成分である。E1 の総和則が最終的に殻内の期待値に落ちる理由は、エッジ積分により終状態が閉包され、基底状態の殻内密度行列 cm1cm2 の線形結合に帰着するからである。

低階数の等価演算子

殻内での物理的意味を把握しやすくするため、低階数の例を挙げる。d 殻(l=2)に限定すると、スピンレス電気四極子(Q2m)は「軌道占有の異方性」を表す。等価演算子の代表形として

Q20(orb)(3lz2l(l+1)),Q22(orb)(lx2ly2),Q2,2(orb)(lxly+lylx)

が用いられる。これらは、例えば d3z2r2(dx2y2,dxy) の占有差、あるいは面内異方性を直接反映する量になる。

8.3 スピンフル多極子

スピンフル(s=1)とは、演算子がスピン空間の一次テンソル(パウリ行列)を伴い、スピンと軌道の角運動量結合で構成されることを意味する。基本式は

X^lm(1;k)=m1,nCκm1,1nlmX^κm1(orb)σ1n,(κ=l+k)

である。

指標 (1;k) の意味

  • 上付きの 1 はスピン側のテンソル階数が 1(ベクトル)であることを示す。
  • k は、最終的な全テンソル階数 l を得るために必要な軌道階数 κ を指定するためのラベルである。

角運動量加算の三角条件より

|κ1|lκ+1

が必要であるから、与えられた l に対して許される κ

κ{l1, l, l+1}

となり得る(非負整数の範囲で)。したがって k=κl

k{1,0,+1}

が基本的に許され、同じ l でもスピンフル多極子が複数系統(最大3系統)存在し得る。この多重性が、従来の「スピンモーメント」だけでは整理し切れない寄与(例:異方性磁気双極子項)が自然に現れる数学的理由である。

スピン球面テンソル成分

σ1n はパウリ行列から作る球面テンソル成分であり

σ00=σ0,σ10=σz,σ1,±1=12(σx±iσy)

である。スピン角運動量演算子 S を用いるなら Sα=(/2)σα によって同値に書けるが、本稿では無次元表示として σ を用いる。

変換性の直観

スピンは

  • 空間反転で偶(軸性ベクトル)
  • 時間反転で奇

である。したがってスピンフル多極子は、軌道部分 X^κm(orb)I,T 性とスピンの I,T 性が合成されたものとして理解できる。結果として、同じ「見かけの階数 l」でも、時間反転に対して奇な成分と偶な成分が混在し得る。XMCD や XMLD が「どの成分を選択的に拾うか」を議論する際に、この点が決定的になる。

9. 従来総和則と完全多極子表現の関係

従来の Thole–Carra 総和則は、L2,3 端の XMCD スペクトル積分から

  • 軌道磁気モーメント mL
  • スピン磁気モーメント mS と異方性磁気双極子項 mT の組合せ

を抽出する形式として確立している。完全多極子表現は、この mT が付け足しの補正項ではなく、「スピンフル多極子の独立成分」として必然的に出現することを明確にし、巨視的磁化がゼロでも XMCD が生じ得る状況を同一の代数構造で扱えるようにする。

9.1 従来の XMCD 総和則が返す演算子

L2,3 端に対して、実験で用いられる積分量を

pL3dωΔμ(ω),qL2+L3dωΔμ(ω),rL2+L3dω(μ+(ω)+μ(ω))

とする(Δμ=μ+μ)。すると従来総和則は、係数を含めて

mL=μBLznhqr,mS+7mT=2μB(Sz+7Tz)nhp2qr

の形で書ける(nh はホール数)。ここで Tz は、スピン密度の異方性を表す演算子の期待値である。

従来の議論では、この Tz が「補正」あるいは「実験幾何で効く項」として扱われることが多いが、完全多極子表現では、Tz はスピンフル多極子の独立成分として体系の中に位置付く。

9.2 完全多極子表現での読み替え

E1 の積分は角運動量代数で

IϵKQPQ(K)(ϵ)O^Q(K)

と書ける。XMCD は K=1 成分(擬ベクトル)を選び出すため、総和則が本質的に返すのは「K=1 の(スピンレスとスピンフルを含む)有効多極子」の期待値である。

  • 軌道総和則は、スピンレス磁気双極子 M10(0;0)(実体としては Lz)を抽出するものとして理解できる。
  • スピン総和則は、スピンフル磁気双極子のうち、複数の k 成分が混合した組合せ(SzTz の線形結合)を抽出するものとして理解できる。

ここで重要なのは、同じ K=1 に属する「ベクトル量」であっても、スピンフルの中に複数の独立成分があるため、XMCD が

Sz0

だけを意味しない、という点である。例えば

Sz=0,Tz0

となり得る場合、従来の表現では「磁化はゼロなのに XMCD が出る」ように見えるが、完全多極子表現では「スピンフル多極子の異方性成分が非零である」として矛盾なく分類できる。

9.3 異方性磁気双極子項の位置付け

異方性磁気双極子項は、直観的には「スピン密度の角度分布が球対称からずれている度合い」を表す。代表的な形として、各電子 i に対して

T=i(si3(sir^i)r^i)

を定義すると、Tz はその z 成分である。ここで r^i は角度依存(Y2m に対応する二階の角度構造)を含むため、Tz

  • スピン(階数1)
  • 軌道角度依存(二階の異方性)

が結合した量である。これは 8.3 節の角運動量結合

(κ1)l

の言葉では、「軌道の異方性(κ=2)とスピン(1)が結合した l=1 成分」を含む、と解釈できる。したがって mT は、完全多極子基底におけるスピンフル多極子の一部として自然に現れる。

9.4 反強磁性体・アルターマグネットへの意味

巨視的磁化がゼロであっても、結晶対称性と局所電子状態の異方性により、時間反転で奇な多極子が許される場合がある。そのとき XMCD は

  • M1m(通常の磁化)だけでなく
  • スピンフル成分(Tz を含む)や、対称性で許された別種の時間反転奇成分

を通じて生じ得る。

従来総和則の式そのものは変えずに、抽出している演算子の「分類」を更新することが完全多極子表現の効用である。これにより、XMCD と XMLD を同一の角運動量代数の言語で扱い、信号の起源を「スピン・軌道モーメント」という二分法から「スピンレス/スピンフル多極子」という整理へ接続できる。

10. 反強磁性体・アルターマグネットにおけるXMCDの読み替え

反強磁性体(巨視的には磁化がゼロ、あるいはほぼゼロ)であっても、XMCDが検出される場合がある。従来の直観では、XMCDは円偏光に対する吸収差であるから「スピン偏極(Sz)が非零である」ことと結び付けられやすい。しかし、総和則が直接に返すのは Sz そのものではなく、スピン項に加えて異方性磁気双極子項(magnetic dipole term)を含む組合せである。したがって、XMCDの存在は Sz0 を必要条件としない。

10.1 総和則が返す量

L2,3 端のXMCD総和則は、実験で定義される積分量 p,q,r を用いて

mL=μBLz,mS+7mT=2μB(Sz+7Tz)

の形に帰着する(係数と nh を含む最終式は前節の定義に従う)。ここで Tz は、スピン密度の角度分布(殻内の異方性)を含む量であり、単なる「磁化」ではない。

完全磁気多極子基底の言語では、この Tz(あるいは等価演算子としての tz)は、スピンフル多極子の成分として体系的に現れる。すなわち、XMCDは「スピンのベクトル成分」だけでなく、「スピンと軌道異方性が結合したベクトル成分」を同じ K=1 チャネルで拾い得る。

10.2 tz の物理的意味と、Sz=0 でもXMCDが出る条件

異方性磁気双極子項は、直観的には「スピン密度が球対称からどれだけ歪んでいるか」を表す。代表的な定義の一つは

T=i(si3(sir^i)r^i),Tz=Tz^

である。ここで r^i が角度依存を担うため、Tz はスピンの向きだけでなく、軌道(電荷密度)分布の異方性と密接に結び付く。従って

Sz=0,Tz0

が起こり得る。このとき、総和則で抽出される (mS+7mT) は非零となり、XMCDが有限に観測され得る。重要なのは、これが「局所的な異方性(対称性と占有の偏り)」に由来するため、巨視的な強磁性(ネット磁化)と同一視できない点である。

10.3 反強磁性体でXMCDが観測されるときの解釈の分岐

反強磁性体でXMCDが観測された場合、少なくとも次の寄与が混在し得る。

  • 微小なスピンカンティングや不完全補償による Sz 成分
  • 外部磁場による誘起磁化(線形・非線形応答)
  • スピン密度の異方性に起因する Tztz)成分
  • 軌道角運動量の非零成分 Lz(交換相互作用・SOC・結晶場の協奏による)

完全多極子表現が与える利点は、上のうち「Tz 成分」を、補正項としてではなく、独立の対称的自由度として扱える点にある。したがって、観測されたXMCDを「強磁性的に見える信号」と短絡せず、どの演算子が許容され、どの演算子が優勢かを、結晶対称性と測定幾何に基づいて吟味できる。

10.4 交替磁性におけるXMCDの位置付け

交替磁性は、巨視的には反強磁性的(ネット磁化がゼロ)でありながら、結晶対称性のためにバンドがスピン分裂し得るコリニア磁性相として提案・整理されている(参考文献のPRX論文群を参照)。ここで本質的なのは、時間反転対称性が磁気秩序により破れていても、磁化ベクトル(M1m)が必ずしも巨視的に許されない場合がある点である。

この種の系に対してXMCDを読む際には、次のように整理すると混乱が減る。

  • XMCDは K=1 チャネルの応答であり、ベクトル量に結び付く。
  • そのベクトル量は、Sz だけに限られず、LzTz を含む。
  • 交替磁性では、対称性により「通常の磁化とは異なるベクトル成分」が許容され得るため、XMCDの非零は必ずしも「磁化がある」を意味しない。

従って、交替磁性や反強磁性体では、XMCDとXMLDを組で測定し、さらに外部磁場・入射角・方位角依存を併用して、K=1(XMCD)と K=2(XMLD)の情報を分離しながら、SzTz の相対寄与を議論するのが有効である。

11. 磁気結晶異方性とXMLD

XMLDは線偏光方向の違いに対する吸収差であり、E1×E1 の枠組みでは K=2(二階テンソル)成分に対応する。したがってXMLDの積分は、殻内の四極子(電荷分布異方性)や、スピンと結び付いた二階テンソル(スピン軌道相互作用の異方性)を選択的に反映し得る。

11.1 XMLDが拾う演算子の階数

E1遷移では

11=012

であり、線偏光差は K=2 成分を抽出する。従って、XMLDの積分は形式的に

dωΔμlin(ω)Q=22ΔPQ(2)O^Q(2)

の形に書ける。ここで ΔPQ(2) は偏光テンソルの差、O^Q(2) は殻内の四極子演算子(スピンレス/スピンフルの双方を含む)である。

11.2 スピンレス四極子と電荷分布異方性

スピンレス四極子は軌道占有の異方性を表し、代表的には

Q20(orb)3lz2l(l+1),Q22(orb)lx2ly2

などで表される。これらは d 殻に対して、例えば面外(3z2r2 成分)と面内成分の占有差、あるいは面内の二軸異方性を反映する。従って、XMLDの積分値は、電荷分布の異方性(結晶場・歪・界面による軌道占有の偏り)を議論するための直接の指標となる。

11.3 スピンフル四極子とスピン軌道相互作用の異方性

スピンフル四極子は、スピンと結び付いた二階テンソルであり、SOCの異方性や磁気異方性に感度を持ち得る。概念的には

  • スピンの向きが変わると、SOC項 λls の有効平均が変わる
  • その変化が、占有の異方性や遷移確率の異方性として現れ、線偏光差に寄与する

という流れである。等価演算子としての候補は複数あり得るが、議論では「どの K=2 成分が測定幾何で抽出されるか」を明示し、そこで許される演算子の対称性を合わせることが重要である。

11.4 磁気結晶異方性との接続の考え方

磁気結晶異方性エネルギーは、SOCと結晶場(および占有)の協奏で決まる。元素選択的に見ると、XASの二色性は「同じSOCを担う殻の占有・遷移確率の異方性」を直接に測っているため、磁気異方性の議論に接続できる。

この接続は、次の二段で理解すると整理しやすい。

  1. XMLD積分 四極子期待値 O^(2)
  2. O^(2) の変化 SOC有効平均の変化 異方性エネルギーの変化

国内解説としてPF NEWSでは、XMCD総和則による軌道モーメントの異方性と、XMLD総和則による四極子量(電荷分布異方性)の評価を対比し、界面垂直磁気異方性の議論へ接続している。この種の整理は、測定で得られる量が「スピン・軌道モーメント」だけではないことを明確にする上でも有効である。

12. 実験データ処理での整合条件

総和則は、エッジ領域のエネルギー積分により閉包関係を用いて基底状態期待値へ落とす構造をもつ。そのため、スペクトル処理の一貫性が崩れると、積分値 p,q,r が系統的にずれ、抽出した演算子の解釈が不安定になる。以下では、総和則の前提と整合するための確認事項を、式とともに丁寧に記す。

12.1 積分範囲

12.1.1 エッジ積分の意味

総和則で用いる積分は、本質的には

Iϵ=edgedωμϵ(ω)T^ϵT^ϵ

という閉包に対応している。ここでの「edge」は、対象とするコア準位からのE1励起が支配的で、かつ同じ終状態空間(例:2p3d3d 空孔状態)を十分に覆う範囲でなければならない。

12.1.2 実務上の積分端の選び方

L2,3 端では L3L2 がスピン軌道分裂で分かれているが、寿命幅や多体効果、連続吸収の立ち上がりにより重なり得る。積分範囲の設定では次が重要である。

  • L3 積分(p)は、L3 の立ち上がりから L2 の立ち上がり直前までを含むようにする
  • L2+L3 積分(q,r)は、L2 を含め、白線強度が十分収束する高エネルギー側まで含める
  • 同一試料の異なる偏光(+/、線偏光方向)の積分端は必ず同じにそろえる

式で書けば、実験定義は

p=L3dω(μ+μ),q=L2+L3dω(μ+μ),r=L2+L3dω(μ++μ)

であるから、積分区間がずれると p,q,r の比が系統的に変化し、抽出する LzSz+7Tz の見かけが変わる。

12.2 正規化

12.2.1 なぜ正規化が必須か

総和則の最終式は q/r(p2q)/r のように比で表されるため、強度スケールが揃っていないと定量ができない。そこで、測定モードに応じた正規化が必要である。

12.2.2 前縁・後縁による基準化

XASでは、前縁(pre-edge)領域で基線を取り、後縁(post-edge)領域で連続吸収を規格化する。例えば、前縁を直線 aω+b、後縁を直線(あるいは緩やかな多項式)cω+d として、

μ~(ω)=μ(ω)(aω+b)cω+d

のように規格化してから積分する。重要なのは、偏光状態ごとに別々の基線を当てるのではなく、同一の基線・規格化条件を適用することである(偏光に依存しない連続吸収は、本来共通因子として扱われるべきであるため)。

12.2.3 検出モード固有の非線形

蛍光収量(FY)では自己吸収、全電子収量(TEY)では表面感度・帯電や飽和効果が問題になり得る。これらは「強度が真の吸収に比例しない」状況を作るため、総和則の前提(遷移確率に比例する強度)を崩す。

確認事項としては次が有効である。

  • 同一試料で TEY と FY の線形性を比較し、エッジ近傍の歪みの有無を確認する
  • 厚膜・高吸収の条件では FY の自己吸収補正を検討する
  • 磁場反転やヘリシティ反転で差スペクトルの符号と大きさが再現するかを確認する

12.3 ホール数 nh

12.3.1 nh が必要になる理由

総和則は、コアから価電子殻(例えば d 殻)への遷移の総強度が「ホールの数」に比例するという事実を用いる。そのため、磁気モーメントの絶対値を求めるには nh が必要になる。

3d 系では

nh=10nd

であり、ndd 殻の占有数である。合金や強い混成がある場合、nd は整数にならず、定義(投影の取り方、基底の選び方)に依存する。

12.3.2 nh の決め方の整理

nh の推定には複数の方法がある。

  • 化学量論・形式価数に基づく見積もり(迅速だが混成が強いと不安定)
  • 第一原理計算の局所投影からの nd 推定(投影の定義依存が残る)
  • 既知の標準試料や文献値に合わせた相対較正(比較には強い)

強い混成や電荷移動がある系では、絶対値の議論よりも、同一手順で得た nh を用いて差分(組成・歪・温度・界面条件に対する変化)を見る方が再現性の高い議論になる。

13. 応用の射程

総和則を完全磁気多極子基底で整理することの応用は、「磁化」や「スピン偏極」だけでは分類できない弱い二色性信号を、対称性と演算子の期待値として読み替えられる点にある。以下に、研究上の展開先を具体的に記す。

13.1 隠れた多極子の同定

巨視的磁化が小さい、あるいはゼロに近い状態でも、スピンフル多極子(とくに異方性成分)が非零になり得る。二色性の積分量を、許容される多極子の表に照らして分類すると、観測信号を「どの多極子が担っているか」という問いに落とし込める。交換バイアスや界面反強磁性のように、界面近傍で対称性が変化する系では、電場や歪の印加が Tz や四極子期待値を変える経路を持ち、二色性の変化を多極子の変化として追跡できる。

13.2 反強磁性体におけるXMCD起源分類

反強磁性体でXMCDが得られたとき、信号の起源は一つではない。総和則で抽出される演算子が

Sz+7Tz

であることを前提に、次の切り分けが可能になる。

  • 外部磁場依存で線形に増える成分は誘起磁化(Sz)の寄与を示唆する
  • 幾何(入射角・方位角)に強く依存する成分は Tz(異方性)の寄与を示唆する
  • XMLDの同時測定で四極子量の変化が見える場合、殻内異方性と連動している可能性が高い

このように、XMCDを単独で「強磁性的」と解釈せず、K=1K=2 を同時に測り、対称性に沿って解釈することで、反強磁性体における信号を定量的に分類できる。

13.3 異方性スピン軌道相互作用の抽出

XMLDは K=2 成分を通じて、スピンレス四極子(占有の異方性)とスピンフル四極子(SOCの異方性に結び付く量)を選択的に反映し得る。面外XMLD・面内XMLDを分けて積分し、それぞれが抽出する Q20Q22、さらにスピンフル二階テンソル成分に対応する演算子の期待値を比較することで、「SOCがどの方向に強いか」という情報を元素選択的に引き出す道が開ける。

13.4 磁気異方性エネルギーへの接続

磁気異方性エネルギーはSOCと占有の異方性で決まるため、XMLDで得た四極子期待値の変化を手掛かりに、異方性の起源(結晶場由来か、界面由来か、歪由来か)を議論できる。さらにXMCDで軌道モーメント異方性が評価できる場合、XMLDで電荷分布異方性、XMCDで軌道応答を同時に扱うことで、異方性の電子論的理解が進む。

13.5 交替磁性物質群への展開

交替磁性では、ネット磁化がゼロであっても対称性によりスピン分裂したバンド構造が許され得る。このとき、二色性信号は「磁化の有無」ではなく「どの時間反転奇成分が許され、どの二階テンソルが共存するか」によって特徴付けられる。完全多極子基底による整理を導入すると、物質ごとに許される多極子の集合と、測定幾何が抽出する成分の対応関係を明示した上で、XMCD/XMLDの弱い信号を読み解ける。

まとめと展望

完全磁気多極子基底に基づく総和則は、XMCD・XMLDのエネルギー積分が抽出する演算子を、スピンレス/スピンフルの多極子として体系的に整理し、XMCDが Sz だけでなく異方性磁気双極子項(tz)を含み得ること、XMLDが四極子期待値(電荷分布異方性やSOC異方性)に結び付くことを、同一の角運動量代数の枠組みで示すものである。

今後は、反強磁性体・交替磁性・界面系のように信号が弱く起源が複合的な系に対して、XMCD(K=1)とXMLD(K=2)の積分量を、結晶対称性と測定幾何に沿って分解し、SzTz、さらに四極子期待値の相対寄与を同時に定量する方向が有望である。これにより、従来のスピン・軌道モーメント解析だけでは到達しにくかった、元素選択的な異方性指標や隠れた秩序の同定へと議論を拡張できる可能性が高い。

その他参考にしたsources

  • P. Carra et al., Phys. Rev. Lett. 70, 694 (1993). doi:10.1103/PhysRevLett.70.694
  • G. van der Laan, Phys. Rev. Lett. 82, 640 (1999). doi:10.1103/PhysRevLett.82.640
  • H. Kusunose, R. Oiwa, S. Hayami, J. Phys. Soc. Jpn. 89, 104704 (2020). doi:10.7566/JPSJ.89.104704
  • S. Hayami, H. Kusunose, J. Phys. Soc. Jpn. 93, 072001 (2024). doi:10.7566/JPSJ.93.072001
  • L. Šmejkal, J. Sinova, T. Jungwirth, Phys. Rev. X 12, 040501 (2022). doi:10.1103/PhysRevX.12.040501
  • L. Šmejkal et al., Phys. Rev. X 12, 031042 (2022). doi:10.1103/PhysRevX.12.031042