《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

共鳴 X 線回折実験と多重極子秩序状態の観測
埼玉大学

研究機構 / 大学院理工学研究科

道村 真司

序論
局在 f 電子系において、 多重極子自由度は超伝導や四極子近藤効果など特異な物性を引き出す未解明部分の多い

重要な秩序変数である。しかし、従来の多重極子研究では、多重極子モーメントの秩序変数や波数ベクトル等を実験で直接的に決定することができず、マクロ物性と理論的解釈に強く依存している状況であった。しかし、2000 年以降の最近では、多重極子秩序状態の観測に対して共鳴 X 線散乱実験の有効性が認知されてきた。その結果、推論ではなく実験的な直接的証拠として多重極子秩序状態とその物性への影響を論じることが出来るようになり、これまでの理論的考
察の検証等、実験と理論の両輪が効率的に機能し始めた。
また、共鳴 X 線散乱の実験環境についても、黎明期は 1K 以下の極低温や高磁場下での実験環境がなく限られた物質に限った実験であったが、2010 年以降、実験環境の整備によりそれらの問題が解決されつつあり、対象物質の幅も広がってきている。
共鳴 X 線散乱は、f 電子系のみならず d 電子系においても幅広く活用される。放射光を用いた共鳴 X 線散乱 (回折)実験は、実験室で扱う X 線回折実験に比べて入射エネルギーや偏光条件等のパラメーターが増えることにより一見複雑であるように感じる。たしかに、共鳴 X 線散乱の理解には、X 線の偏光や分散補正項、蛍光など多くの知識が必要であるが、一つ一つを理解し整理することにより自分で実験をすることも容易である。また、その知識は X 線による磁気構造解析 (共鳴 X 線磁気散乱実験) や価数評価 (X 線吸収分光実験) などに必要な基礎知識である。共鳴 X 線散乱を理解することにより、多くの X 線分光実験に対しても理解が容易になるに違いない。
このテキストでは、
第一章：X 線のエネルギーと偏光電磁場
第二章：古典的な X 線の散乱過程と吸収(電子の散乱・回折現象の簡易モデルによる Thomson 散乱の導出 及び Rayleigh 散乱, Compton 散乱, 発光の解説)
第三章：量子力学な X 線の散乱過程(電子と電磁場との相互作用ハミルトニアンによる Thomson 散乱, 非共鳴磁気散乱, 共鳴散乱の導出)
第四章：非共鳴 X 線回折実験の実際と解析
第五章：共鳴 X 線回折実験の構成となっている。
第一章, 第二章では、これまで学部で学んだ X 線回折の復習とともに、回折現象において X 線の偏光が意味をもつことや吸収や発光等が共鳴散乱とどのような関係があるかを意識するきっかけを記したつもりである。
第三章では、古典的な解釈では定量的に表現できない共鳴散乱を二次摂動を取り入れることにより、表式化し Thomson
散乱, 非共鳴磁気散乱, 共鳴散乱の３つの項に分離できることまでを示した。
第四章以降では、実際に解析で用いる表式の紹介と計算例を示した。共鳴 X 線回折実験には、回折実験の基礎的な知識が重要である。このテキストでは第一章, 第二章に重点をおいており、基礎的な知識をもって頂きたい。第四章以降に関しては、広島大学の松村准教授が丁寧にまとめられている [9]。このテキストと比較しながら、松村准教授のテキストを参照して頂きたい。

39

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

1

X 線のエネルギーと偏光

1.1

電磁波とエネルギー

X 線は電場 E と磁場 H の互いに垂直な電磁波である。簡単のため、電場を単純な調和振動 (振動数 ω, 振幅 E0 ) とし、Maxwell 方程式より磁場を電場を用いて表わすと
E
H= (0, E0 sin(ωt), 0)√
ε0
E0 sin(ωt))= (0, 0,
µ0

(1.1)
(1.2)

となり、真空の誘電率 ε0 と真空の透磁率 µ0 を用いた関係により、電場か磁場のどちらかが決まれば、他方も決まる。
ここでは、電場の方向を y 軸としている。
また、空間上の電磁波を考える上では、式 (1.1) のような振動数 ω や波長 λ よりも式 (1.3) のような伝播ベクトルk(|k|=2π/λ=ω/c) と位置 r を用いる場合が多い。

E(r) = ε̂E0 exp(ik · r)

(1.3)

電場方向の単位ベクトル ε̂ は偏光の単位ベクトルであり、k は逆格子空間を考える上でも非常に有用な量である。
いずれにしても、電場と磁場のエネルギーは等しく、以下のように表記できる。

UE
UH

1
ε0 E022
1
1= µ0 ⟨H 2 ⟩ = µ0 H02 = ε0 E02 = U E2
2= ε0 ⟨E 2 ⟩ =電磁波としてのエネルギーはポインティングベクトル S を用いて表わされ、√
1 ε0 2⟨S⟩ = ⟨E × H⟩ =E
2 µ0 0

(1.4)
(1.5)

(1.6)

である。ポインティングベクトルは、単位時間あたりに単位面積を通過する電磁場のエネルギー (エネルギー密度) である。
図. 1.2.1 に電磁波の波長とエネルギースケールを示した。X 線は 10−2 Å ∼102 Å の波長領域の電磁波であり、エネル
ギースケールにすると 数 eV ∼ 120 keV である。実験室系 X 線回折実験で利用される X 線管球は表. 1[1, 2] のように

5 keV ∼ 80 keV である。また、電磁波のエネルギー E の波長 λ と振動数 ω との関係は
E=

hc
= h̄ω
λ

(1.7)

ここで、h, c, h̄ = h/2π はプランク定数, 光の速度, ディラック定数である。
余談ではあるが、中性子や電子等の粒子の波のエネルギーはE = (h̄k)2 /2m(1.8)

ここで、m は中性子や電子など の粒子の質量である。
実際に物理定数を代入すると、波長 [Å] と次の関係が得られる。
電磁波 ： 12.4/λ[keV ](1.9)

中性子 ： 150.4/λ2 [eV ](1.10)

2

電子 ： 81.8/λ [meV ]40(1.11)

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

図 1.2.1 電磁波の波長とエネルギースケール

1.2

X 線の偏光

共鳴 X 線回折ではこの偏光状態が非常に大切である。偏光状態は、直線偏光と円偏光, 楕円偏光の三つに大別され、
水平偏光、垂直偏光、45 度偏光、右回り偏光、左回り偏光が代表的である。(前者３つが直線偏光、後者２つが円偏光。)
これらの５つの偏光状態は２つの電場の重ね合わせで記述できる。図. 1.2.2 には、偏光状態の模式図を示す。z 軸方向に伝播ベクトル k をとる。電場は k に直交するので、x, y 軸の電場成分 (Ex , Ey ) は

Ex

Ey

=

Ex0 cos(ωt)=

E0 cos(ϕ)cos(ωt)=

Ey0 cos(ωt + δ)=E0 sin(ϕ)cos(ωt + δ)

(1.12)と表わすことができる。E0 は電場の絶対値、ϕ, δ は (Ex ,Ey ) の振幅比と位相差のパラメーターである。
表 1 X 線管球とエネルギー線源

Kα1 (keV)

Kα2 (keV)

Kβ1 (keV)CrFe

5.41472
6.403845.405509
6.390845.94671
7.05798

CoCu

6.93032
8.047786.91530
8.027837.649438.90529Mo
W

17.47934
59.31824

17.3743
57.981719.6083
67.244341(1.13)

《講義ノート》

1.2.1

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

直線偏光

直線偏光は電場ベクトルの方向が一定である状態である。そのため、Ex , Ey の位相差 δ は 0 または ±π である。δ=0とすると、ϕ は電場ベクトルの水平面からの角度となる。したがって、

δ = 0◦ , ϕ = 0◦

:

水平偏光

(1.14)

δ = 0◦ , ϕ = 45◦: 45◦ 偏光

(1.15): 垂直偏光(1.16)◦

δ = 0 , ϕ = 90

◦

となる。

1.2.2

円偏光

円偏光は電場ベクトルの大きさが一定、かつ xy 面内で回転している状態である。そのため、ϕ=45◦ , δ=±π/2 である。回転方向は δ の符号に依存し、したがって、

δ = +90◦ , ϕ = 45◦: 右周り円偏光(1.17)

δ = −90◦ , ϕ = 45◦:

左周り円偏光(1.18)

となる。

1.2.3

楕円偏光

式 (1.12) で表わされる偏光状態の内、直線偏光でも円偏光でもない両者が混じった偏光状態は楕円偏光となる。この時、(ϕ, δ) の条件は (ϕ, δ)=(0, 任意), (±π, 任意), (±45◦ , ±π), (∓45◦ , ±π) を除く全てである。当然、ϕ は長軸の水平面からの角度、δ は短軸と長軸の比となる。

1.2.4

無偏光

これまで紹介した偏光状態の電場ベクトルは、ある時間では特定の大きさと方向を向いており、それらの時間変化は周期的性をもっていた。この状態を「完全偏光」と呼ぶ。一方で、全く周期性のない偏光状態を「無偏光」と呼ぶ。無偏光状態は、電場ベクトルの大きさと方向が時間的空間的に揺らいでいると考えれば良いだろう。現実には、波動である光のなかには必ず周期的成分が存在し、完全な無偏光状態はない。逆に、揺らぎを 0 にすることも難しく、完全な偏光状態もない。そのため、全ての光は偏光状態と無偏光状態の混ざった偏光状態となっている。

42

《講義ノート》

図 1.2.2

2

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

代表的な偏光状態とイメージ図 [3]。電磁波のイメージ；時間を止めたときの電場ベクトル。偏光楕円；電磁波の進行方向から見たときの電場ベクトルの時間変化の軌跡。(Ex0 , Ey0 )；電場ベクトルを直交する x 軸成分と y 軸成分に分離し、時間を止めたときの電場のベクトル。

X 線の古典的な散乱過程と吸収共鳴 X 線回折実験の解析には量子力学的な記述を理解しなければならない。しかし、量子力学的な記述でイメージを掴むことは難しい。そのため、多くの大学での X 線回折に関する講義は、古典的な記述での理解から始める。実際、X
線管球を線源とした実験室系での X 線回折実験を理解するには古典的な記述で十分である。しかし、限られた講義時
間では、X 線回折における偏光や吸収、発光等について深く踏み込むことはあまりない。共鳴 X 線回折実験では、この偏光と吸収の理解が重要である。本章で、「実験室系での X 線回折実験でも偏光が関係していること」と「散乱, 吸収, 発光についての理解」を古典的な記述でイメージしてほしい。
散乱, 吸収, 発光について概要を述べる。原子に入射した X 線は、原子内の電子と相互作用し、散乱, 吸収, 発光の 3つに大別された現象を生じる。散乱は、入射 X 線の電磁場としての性質を強く反映する。特に、弾性散乱であるトムソ
ン散乱は電磁波と電子との相互作用として説明できる。一方、非弾性散乱であるコンプトン散乱は光子と電子との衝突(相互作用) として理解され、X 線の粒子性を反映している。弾性/非弾性に関わらず、散乱波は X 線と相互作用した一つの電子からの放射である。そのため、古典的に密度分布をもった複数の電子を構成要素にもつ原子による散乱波は、
電子分布中の微少体積からの散乱の重ね合わせとみなせる。
残る吸収と発光の 2 つは、コンプトン散乱と同様に X 線の粒子性を強く反映する。原子と X 線との相互作用で生じる吸収は、光電吸収と呼ばれ、原子内の電子と光子の非弾性的衝突であり、光子は電子にエネルギーを与えることで消失する。光電効果により飛び出した光電子のエネルギーは、入射光子のエネルギーに比べて結合エネルギー分のロスをもつため、非弾性的現象である。光電子の飛び出しにより、球面波状に広がる波を散乱するため、非弾性散乱と捉えることもできる。発光は主に光電効果で生じた内殻の空孔への外殻電子の遷移に伴う光子の生成である。遷移経路により蛍光や燐光と区別されるが、これら発光は吸収に伴う現象であり非弾性的現象である。吸収と発光は原子内に束縛された電子を考慮した現象であるため、１つの自由電子と１つの光子の現象ではない。

43

《講義ノート》

図 2.1.1

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

電気双極子モーメントによる X 線の散乱波 [4]。y 軸方向から入射された電場により原点位置の電子が z 軸方向に双極子振動し、散乱波を放射している。ただし、放射波は x<0 のみを表示している。

2.1

散乱

2.1.1

散乱と回折の概要

X 線回折実験にあたり、区別の曖昧な表現として、散乱と回折の違いがある。これらの使い分けは、明確にされていない場合が多い。散乱と回折の違いを考えるために、電子、原子、結晶を散乱体として説明する。まずは、X 線を電子に入射した場合、原子に入射された X 線は原子内の電子等の散乱体に吸収され、その散乱体を振動させた結果として X線が放射される。入射された X 線が、そのまま散乱 X 線として放射されることはないので、勘違いしないようにしま
しょう。この放射された X 線は、特定の方向で強く観測されることはなく、基本的にあらゆる方向で観測される。基本
的に全方位に放射されるこの現象は「散乱」といえる。もちろん、あらゆるといえども規則性はある。例えば、散乱体が電子一つであり弾性散乱を生じる場合を Thomson 散乱と呼ぶが、Thomson 散乱では、図. 2.2.1[4] のような電子一つを中心に放射されるドーナツ状の散乱波になる。さらに、原子核と電子の電磁相互作用による束縛状態である原子を散乱体とする弾性散乱の散乱波は、原子内に存在する複数の電子の電子分布を古典描像として数密度関数で表し、電子分布中の微少体積要素での Thomson 散乱波の重ね合わせと考える。電子の数密度関数は連続的に変化をするため、電子と同様に原子でも干渉の結果としての散乱振幅が特定の方向で強まることはない。この散乱振幅は、原子散乱因子または原子形状因子と呼ばれる。
一方、原子が周期的に配列した結晶では、個々の原子からの散乱波が干渉し合い、散乱波の振幅は位相の揃うある特定の方向で強くなる。この干渉により特定の方向の散乱 X 線が強くなる現象が「回折」である。この回折現象の起きる条件は Bragg 条件としてよく知られている。波長 λ をもつ入射 X 線が周期間隔 (面間隔 d) をもつ結晶に Bragg 条件(λ=2dsinθ) を満たす θ で入射すると、散乱角 2θ 方向で位相が揃った回折 X 線を得られる。この周期性を取り入れた際の散乱振幅は、構造因子または結晶構造因子と呼ばれている。
物性実験 (特に物質合成) を生業とする多くの研究者が、「結晶構造を調べる」ために X 線回折実験を行う。その際、当然我々は回折 X 線を観測するが、観測した回折 X 線は、結晶を構成する原子による散乱と周期性を含んでいる。

44

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

図 2.1.2 電子に X 線を入射した際の散乱模式図

2.1.2

散乱断面積

「散乱」を議論する上で、重要な量は「微分散乱断面積」である。微分散乱断面積は、粒子 1 つがもつ散乱効率である。図. 2.2.2 に散乱体 (電子) に X 線を入射した際の散乱模式図を示す。ここで、入射X 線と散乱 X 線の電場を E in ,√E rad とする。このとき、入射 X 線のエネルギー密度は 21

ε0 2
2µ0 Ein であり、電場の振幅の二乗 (Ein ) に比例することよ2
り、光子の数密度 (単位面積当りの光子数) は Ein/h̄ω に比例する。したがって、単位面積、単位時間当りの光子数であ2
る入射 X 線強度のフラックス Φin は、光子の速度は c と電場を用いて、cEin/h̄ω に比例する。

Φin ∝

c 2
E
h̄ω in(2.1)

一方、実験での散乱 X 線強度は試料から距離 R 離れた検出器で観測される。検出器の立体角を ∆Ω とすると、検出器の面積は R2 ∆Ω である。よって、検出器で観測される散乱 X 線の X 線強度 Irad はIrad ∝c 2
2
R ∆ΩErad
h̄ω

(2.2)

と表わすことができる。先程述べたように、微分散乱断面積は、粒子 1 つがもつ散乱効率である。散乱 X 線強度を立
体角 ∆Ω で考えたので、散乱体に実効的に照射された入射 X 線強度も立体角当りで考えると ∆ΩΦin となる。よって、散乱効率である微分散乱断面積 (dσ/dΩ) は以下のように表わせる。

dσ
dΩ

=

Irad
|E 2rad |R2=∆ΩΦin|E 2in |(2.3)

また、微分散乱断面積を立体角で積分した σ は全断面積と呼ばれ、dσ/dΩ と同様に、入射 X 線がどの程度の割合で散乱されているかの尺度となる。

45

《講義ノート》

2.1.3

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

電磁波のスカラーポテンシャルとベクトルポテンシャル

Maxwell 方程式より電磁波のスカラーポテンシャル ψ とベクトルポテンシャル A は∫
ρ(ξ, η, ζ, t − rc )1
dv
ψ(x, y, z, t) =
(2.4)
4πε0
r
∫
J (ξ, η, ζ, t − rc )µ0
A(x, y, z, t) =
dv
(2.5)
4π
r
√
である。ここで、ρ は電荷密度、J は電流密度である。また、r= (x − ξ)2 + (y − η)2 + (z − ζ)2 , dv= dξdηdζ であり、
(x, y, z, t) が ψ や A の座標や時刻であるのに対して、(ξ, η, ζ, t) は積分点の座標や時刻となる。(ξ, η, ζ, t) を用いて、時刻 t でのポテンシャル ψ や A を求める場合には、積分点での時刻 (t − rc ) に遡って積分する。 rc は、積分点から ψ やA の座標までの光の移動所要時間である。
磁場 H と電場 E の Ψ との関係は、Maxwell 方程式より、
HE= rotA(2.6)∂= −gradψ − A∂t∫= c2 divAdt(2.7)

であるため、電磁波のベクトルポテンシャル (あるいはスカラーポテンシャル) は、H と E を求めるための重要な量で
ある。

2.1.4

X 線による 1 個の電子の振動 (電気双極子)

原子中の電子に X 線を入射した時、電子は振動数 ω の入射 X 線により振動し、散乱波を放射する。また、原子中の電子は固有の振動数 ω0 をもち、電磁波の放射による減衰力が存在する。簡単のため、図. 2.2.5 のように、z 軸方向に調和
振動的に直線偏光した X 線を原子周りの電子 1 つに入射した場合を考えるとその運動方程式は、電子の質量 me として

me

dz
d2 z
+ 2me ρ + me ω02 z
dt2
dt= eEin cos(ωt)(2.8)となり、入射 X 線の偏光と原子中の一つの電子の振動は同軸方向の調和振動となる。電子の調和振動の運動は、古典的に原子核 (電荷:+e) 周りの電子 (電荷:−e) が振動しており、z(t) は原子核と電子の相対距離と考えるとよい。第二項が電磁波の放射による減衰力であり、第三項は原子中の電子の固有振動 (非減衰) である。ただし、ρ=e2 ω 2 /12πε0 me c3 。
また、入射 X 線による強制振動は右辺である。
このとき、原子核と電子は時間的に変化する電気双極子p(t) = (0, 0, ez(t))(2.9)とみなせる。
したがって、z 軸方向に調和振動的に直線偏光した X 線を原子周りの電子 1 つに入射すると、その電子は式 (2.9) のような z 軸を主軸とした電子双極子的な振動となる。この時間依存性をもつ電気双極子モデルは、H.Hertz により立てられ、Hertz の双極子と呼ばれる。古典力学において、入射した X 線の電子による散乱波の正体は、この電子双極子的な振動運動により放射される電磁波である。

2.1.5

電気双極子が放射する電磁波

では、この電子双極子的な振動運動による電磁波は、どのような電磁波となるか。式 (2.6) に示したように、H と E
は電磁波のベクトルポテンシャル (あるいはスカラーポテンシャル) により求めることができる。

46

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

電気双極子振動では、電流密度 J は電荷密度 ρ と電荷の相対移動速度 v を用いてJ = ρv(2.10)となる。したがって、式 (2.4) と式 (2.9) より、電気双極子のベクトルポテンシャルはA(x, y, z, t) =
=
=
=
=

∫
J (ξ, η, ζ, t − rc )
µ0
dv
4π
r
∫
r
µ0 v(ξ, η, ζ, t − rc )ρ(ξ, η, ζ, t − )dv
4π
r
c
µ0 e
r
v(ξ, η, ζ, t − )4π rcµ0 e ∂rz(ξ, η, ζ, t − )4π r ∂tcrµ0 ∂p(ξ, η, ζ, t − )4πr ∂tc(2.11)であり、H と E は式 (2.6) に式 (2.11) を代入し、

H

=

E

=

µ0
∂
r
rot( p(ξ, η, ζ, t − ))
4πr
∂t
c
r
1 ∂2
r
µ0
{grad div p(ξ, η, ζ, t − ) − 2 2 p(ξ, η, ζ, t − )}4πε0 rcc ∂tc(2.12)(2.13)、極座標表示で示すと

r
1 sinθ ∂ 2
pz (ξ, η, ζ, t − ))2
4πc r ∂tc
µ0 sinθ ∂ 2rpz (ξ, η, ζ, t − ), 0)E(r, θ, ϕ) = (0,4π r ∂t2c

H(r, θ, ϕ) = (0, 0,(2.14)(2.15)

つまり、電子の運動方程式 (式 (2.8)) の解を用いた p(t) を求めればよい。
式 (2.8) の解を用いるとpx (t) =ex(t) = 0py (t) =ey(t) = 0(2.16)(2.17)
2pz (t) =ez(t) =

e Ein
1 cos(ωt + δ)2
me [(ω0 − ω 2 )2 + 4ρ2 ω 2 ] 2(2.18)≡ Ccos(ωt + δ)
(2.19)となるので、

HϕEθ

1 sinθ ∂ 2rpz (ξ, η, ζ, t − )24πc r ∂tc
1 sinθ
2(−Cω cos(ωt + δ))=
4πc r
e2 sinθ ω 2 Ein cos(ωt + δ)
= −4πcme r [(ω02 − ω 2 )2 + 4ρ2 ω 2 ] 12

=

= −

sinθ
ω 2 Ein cos(ωt + δ)e2
4πε0 c2 me r me [(ω02 − ω 2 )2 + 4ρ2 ω 2 ] 21

である。

47

(2.20)(2.21)

《講義ノート》

図 2.1.3

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

原子中の電子による X 線の散乱 (a)z 軸方向に調和振動的に直線偏光した X 線 (伝播ベクトル k//y 軸, 偏光ベクトル ε) を原点にある原子周りの電子へ入射している。電子は z 軸方向に電気双極子 p(t) を作り、散乱 X 線 (k’, ε’) を放射する。

したがって、原子中の電子一つに散乱された z 軸方向の調和振動的な直線偏光は、式 (2.20), (2.21) の様な電磁波となる。ただし、

C

=

tanδ

=

ρ =

e2 Ein1me [(ω02 − ω 2 )2 + 4ρ2 ω 2 ] 22ρω− 2(ω0 − ω 2
π e2 1e2 ω 2=3 ε0 me c λ12πε0 me c3(2.22)(2.23)(2.24)である。
図. 2.2.1 は式 (2.21) で表わされる電場の模式図である。電場は sinθ に依存するため、xy 面内にドーナツ型上に放射され、z 軸方向には放射されない。

2.1.6

Thomson(トムソン) 散乱

通常、我々が結晶構造解析や共鳴 X 線回折実験を行う場合、十数 keV から数 keV の X 線を用いる。式 (2.21) における原子内の電子の固有振動数 ω0 は、このエネルギー領域の振動数 ω に比べて無視できるほど小さく、ω ≫ ω0 である。ここで、吸収端から十分に離れたエネルギー (非共鳴エネルギー) での散乱を考える。吸収端については、後述の第 2.2 節を参照されたい。非共鳴エネルギーでは、電子軌道への電子の束縛が小さく ω0 → 0 とみなせる。これは、原子核による電子の束縛が小さい自由電子状態とも言える。

48

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

ω ≫ ω0 , ω0 → 0 であるから、式 (2.20), (2.21) はHϕEθ

e2 sinθEin cos(ωt + δ)
4πcme rsinθe2
Ein cos(ωt + δ)
= −4πε0 c2 me rsinθEin cos(ωt + δ)
= −re
r
= −

(2.25)(2.26)とみなせ、エネルギーに依存しない散乱であることがわかる。ここで re は、電子の静止エネルギー (me c2 ) = 半径

re の球面内に分布した電荷のエネルギー (e2 /4πε0 re ) から仮定される古典電子半径である。また、入射 X 線の電場 ( =Ein cosωt) に対して、散乱波は E θ ∝ −Ein cos(ωt + δ) より、散乱波は入射波と同じ振動数 (=エネルギー) である弾性
散乱だということが明らかである。さらに、式 (2.23) より tanδ = +∞ であるため、散乱波の位相は入射波に比べて-πずれている。
入射 X 線強度 Iin と散乱 X 線強度 Irad は

1
2

√

ε0 2
E
µ0 in

Iin

= ⟨S in ⟩ =Irad= ⟨S rad ⟩ = ⟨E × H⟩ = (= Iin re2

(2.27)
1
2

√

ε0 2
e2
sin2 θ
Ein )(−)2 22µ0
4πε0 c mersin2 θ
r2(2.28)

より、
dσは散乱断面積 dΩdσdΩ=Irad r 2= re2 sin2 θ
Iin(2.29)となる。
このような条件下 (ω ≫ ω0 , ω0 → 0) での電子と電磁波との散乱を「Thomson(トムソン) 散乱」と呼ぶ。実験室系のX 線回折実験で観測される X 線回折の散乱因子のほとんどが Thomson 散乱である。
吸収端近傍でのエネルギー (非共鳴エネルギー) における散乱では、エネルギーに依存しない Thomson 散乱項に加えてエネルギーに依存する分散補正項が加わる。分散補正項に関しては、後述の第 2.4.3 節にて触れる。
余談ではあるが、Rayleigh(レイリー) 散乱は、被散乱体と散乱体が X 線 (数十 eV ∼ 数 eV) と電子である Thomson 散
乱と違い、被散乱体と散乱体は可視光線 (3 eV ∼ 1.5 eV) と空気中の気体分子である。そのため、散乱条件は (ω0 ≫ ω)となる。この条件下での散乱波は、E θ ∝ λ−2 Ein cos(ωt + δ) であり、Thomson 散乱と同じ弾性散乱である。しかし、電場の振幅は入射波長の λ−2 に比例しているため、散乱強度は λ−4 のエネルギー依存性をもち、その位相差は 0 である。青色の波長 λ青 は赤色 λ赤 の約 0.5 倍 (λ赤 ∼ 2λ青 ) である。そのため、青色の散乱強度は赤色に比べて 16 倍程度大きい。我々は、昼間は散乱された青色の光を青空として観測し、夕方は散乱されにくい (青色が散乱され少なくなった)
赤色の光を夕日として観測している。
ついでに、Mie(ミー) 散乱は (ω0 ∼ ω) の条件下で生じる散乱であり、被散乱体と散乱体は可視光線 (3 eV ∼ 1.5 eV)と空気中の液体分子 (水滴等) である。散乱振幅の共鳴的なエネルギーに対応し、幅広いエネルギーで同程度の散乱強度をもつ。位相差は約 π/2 であり、ω0 = ω の時に π/2 となる。
図. 2.1.4 に各散乱の双極子放射の電場振幅の係数 C と入射エネルギーの関係を示す。各散乱の特性が良くわかる。

2.1.7

Compton(コンプトン) 散乱

Compton(コンプトン) 散乱は、Thomson 散乱や Rayleigh 散乱と違い、散乱波と入射波の振動数 (エネルギー) の異なる散乱である。非弾性散乱と呼ばれるエネルギー保存がされない散乱過程の一種である。Thomson 散乱や Rayleigh
散乱は、電磁波と散乱体との相互作用であり、散乱体に散乱された電磁波はそのエネルギーのほぼ全てを電子の振動に49

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

図 2.1.4 双極子放射の電場振幅の係数 C と入射エネルギーの関係 (a)Thomson(トムソン) 散乱 (b)Rayleigh(レイリー) 散乱と Mie(ミー) 散乱

変換していた。一方、Compton 散乱は電磁波の波としての性質に加えて、粒子としての性質を併せもった散乱過程である。図. 2.1.5(a) と (b) に弾性散乱と非弾性散乱の模式図を示す。弾性散乱は電子を弾きだすことはなく電子はその場で振動することにより散乱波を放射する。一方、非弾性散乱は、電子が弾き出された際に散乱波を放射するが、電磁波のエネルギーの一部は弾き出された電子の運動に利用されるため、散乱波のエネルギーは入射波に比べて小さくなる。
実際の Compton 散乱は、一定エネルギー以上の X 線を当てなければ生じない。このことから、光は量子化した粒子(光子) である仮説が立てられ、Arthur Holly Compton が実験的にも証明している。
また、Compton 散乱とは逆に、高エネルギーの電子が低エネルギーの光子 (マイクロ波や赤外線) に散乱される逆Compton 散乱も存在する。散乱体である光子はよりエネルギーの高いＸ線やγ線となる。

2.2

吸収

2.2.1

原子内のエネルギー準位

吸収過程を考える前に、原子内のエネルギー準位を確認しておこう。原子内の電子は、各電子軌道のエネルギー準位に束縛されている。スピン-軌道相互作用によるスピン角運動量と軌道角運動量の合成角運動量の量子数を j とすると、各電子軌道のエネルギー準位は、j と軌道角運動量の量子数 l, 主量子数 n に依存する殻で整理すると模式的には
図. 2.2.1 のようになる。p 軌道や d 軌道は、スピン-軌道相互作用によって 2 つのエネルギー準位に分裂している。分
裂後の各軌道を、j の値を用いて、2p1/2 , 2p3/2 と表わした。

2.2.2

光電吸収と吸収端図. 2.2.2,. 2.2.3 に吸収過程の模式図を示す。図. 2.2.2 では、原子に吸収されたＸ線 (光子) のエネルギーは電子に供給され、内殻電子が原子から飛び出している。その結果、原子がイオン化される吸収過程は光電吸収として知られる、原子から飛び出した電子は光電子と呼ばれ、光電子分光実験では光電子のエネルギーから内殻電子の状態を解析している。また、光子のエネルギーの大きさによっては、内殻電子は連続エネルギー帯内の空軌道へ遷移する (図. 2.2.3 左図
から中央図)。空軌道への遷移は吸収効率の変わり目でもあり、空軌道への遷移エネルギーは吸収端エネルギーと呼ばれる。。例えば、Al の L 殻の 2s 軌道の電子が M 殻の 3p 軌道へ励起する吸収を L1 吸収端と呼び、その励起エネルギーを L1 吸収端エネルギーと呼ぶ。同様に、Ce の K 殻の 1s 軌道の電子が O 殻の 5p 軌道へ励起する吸収は K 吸収端であり、その励起エネルギーは K 吸収端エネルギーと呼ばれる。図. 2.2.4 に Na の K 吸収端近傍の透過率のエネルギー依存性を示す。透過率は透過 X 線強度を入射 X 線強度で割った値であり、Na の厚さは 2µm で計算している。Na では、K吸収端である 1.0721keV 以下のエネルギーの X 線を入射しても L 殻以上のエネルギー準位からの電子励起しかないが、50

《講義ノート》

図 2.1.5

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）(a) 弾性散乱；(例；Thomson 散乱：入射 X 線の電場により振動させられた電子が散乱 X 線を放射する。入射 X 線のエネルギーと散乱
X 線のエネルギーは等しい。)(b) 非弾性散乱 (例；Compton 散乱：運動量 h̄k[エネルギー h̄ck] を持った光子が、電子に衝突し、運動量 h̄k′ となって散乱される。電子は衝突時に運動量 h̄(k − k′ ) をもち、散乱後の光子はエネルギー h̄c(k − k′ ) をロスする。) (c) 吸収-発光の放射過程 (例；蛍光 X 線：入射 X 線のエネルギーを吸収した内殻電子が外殻に遷移する。空孔となった内殻へ外殻電子が遷移
し、遷移と尾ともに蛍光 X 線を散乱する。)

51

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

図 2.2.1 各軌道電子のエネルギー準位と蛍光 X 線放射の遷移。矢印は蛍光 X 線放射の遷移。

K 吸収端以上のエネルギーの X 線を入射すると K 殻の電子も励起することが可能になり急激に透過率が下がる。吸収端エネルギー以下では内殻電子の他軌道への遷移が少なくなり、エネルギーと共に透過率は緩やかに減少する。

2.2.3

蛍光 X 線と燐光 X 線の放射吸収により内殻電子が高エネルギー準位の空軌道や原子外へ励起された結果、励起後の内殻には空孔が生じる。空孔が生じた内殻には、外殻から電子が遷移し再び内殻が満たされる。電子がエネルギーの高い最外殻から低い内殻へ遷移するため、開放された外殻電子と内殻電子の結合エネルギーの差とほぼ同じエネルギーをもつ光子が放射される。

(図. 2.2.3) これが蛍光 X 線である。電子のもつエネルギーは遷移と共に熱エネルギーとしても放射されるため、蛍光 X
線は吸収端エネルギーに比べてわずかに低い。フッ素の蛍光 X 線エネルギー (Kα1 [2p1/2 →1s1/2 ]) は 676.8 eV であるが、K 吸収端エネルギーは 696.7 eV と 10 eV 低い。
図. 2.2.5 には、F 原子に X 線を入射し K 殻電子が M 殻へ遷移した場合を仮定し、蛍光と燐光に関する過程を模式的に示した。高エネルギー準位の空軌道に励起された電子は最外殻 (励起状態で最も低いエネルギー準位) まで落ちてくる(図. 2.2.5(a)-(c))。この時、光子の放射はなく熱エネルギーの放射のみの無放射過程である。電子遷移では、遷移前後の電子スピンの向きは同じ (スピン選択律 : スピン角運動量の変化 ∆ sz =0) であるので、図. 2.2.5(c) の状態から (d) の K
殻へ蛍光 X 線を放射しながら遷移する。蛍光 X 線の放射過程はこのスピン選択律に従っており、図. 2.2.5(c) のように遷移する電子と遷移先の電子の向きが対 (逆) であるため、励起一重項状態と呼ばれる。図. 2.2.5 の場合は 2p3/2 →1s1/2なので、Kα2 である。電子遷移は内殻から K 殻へもあり得るので、2p1/2 (L1) →1s1/2 (K) も取り得る。その場合は Kα1と呼ばれる蛍光 X 線が放射される。しかし一方で、励起した状態のスピンの向きが変化してもエネルギー差が少ない
場合やエネルギーを得する場合は、図. 2.2.5(c) の状態から (e) の状態のようなスピンの状態変化を取り得る。これは項間交差と呼ばれる。(e) の状態から K 殻へ遷移する際にも X 線を放射するが、その際の X 線は蛍光と区別して「燐光」と呼ばれる。(e) の状態から K 殻への遷移はスピン多重度が異なる２つの状態間の遷移であるため、(e) の状態は三重項状態とみなせる。励起三重項状態となる経緯については、スピン軌道相互作用により励起された際にスピン選択律が破られる場合や励起一重項状態と励起三重項状態のエネルギー準位が近いために項間交差を生じる場合があるが、いずれにしても励起三重項状態から遷移する際にはパウリの排他率を満たすために遷移する電子スピンの向きを交換する過程を踏むために遷移確率は低く燐光として放射されるエネルギーは蛍光よりも低い。また、燐光は (e) の状態を経由するため、蛍光よりも時間的に遅れて放射される。
まとめると蛍光：スピン多重度が同じ 2 つの状態間での遷移に伴う放射52

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

図 2.2.2 X 線の吸収と光電子の放出。光電子分光の Na 原子に X 線を照射すると、外殻電子が弾き飛ばされ Na+ イオンとなる。

図 2.2.3 X 線の吸収と蛍光 X 線の放出。

53

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

図 2.2.4 Na の透過率のエネルギー依存性 (K 吸収端近傍) の計算値 [5]。Na の厚さは 0.2µm で計算している。

燐光：スピン多重度が異なる 2 つの状態間での遷移に伴う放射
である。

2.2.4

Auger(オージェ) 電子の放出

遷移する際の開放エネルギーは, 蛍光 X 線や燐光 X 線のように X 線として放出されるだけではない。図. 2.2.6 のよ
うに、L 殻から K 殻へ遷移した開放エネルギーが別の外殻 (ここでは M 殻) の電子を原子外に放出する場合もある。このように二次的に放出された電子を Auger(オージェ) 電子と呼ぶ。
以上のように、吸収により、光電子や Auger 電子の放出、蛍光 X 線と燐光 X 線の放射が生じる。

2.3

遷移選択則と E1,E2 遷移

電子が他軌道へ遷移する場合、始状態 (|i⟩, 遷移前) と終状態 (|f ⟩, 遷移後) を結びつける行列要素が必要である。波
動関数の異方性を最も簡略化した近似は双極子近似であり、遷移行列要素は電気双極子モーメント R=-el(r)(e は電荷,

l(r) は正電荷と負電荷の距離) である。
電気双極子モーメントによる遷移確率は、始状態と終状態の波動関数で遷移行列要素で挟んだ際の空間積分の二乗であるから、始状態と終状態の波動関数を ψi , ψf 、空間座標を τ とすると、∫⟨f |R|i⟩ = ψf∗ · R · ψi dτ

(2.30)の二乗に比例する。奇関数の空間積分は 0 であるため、奇関数の電気双極子モーメントを遷移行列要素にもつ E1 遷移の終始の波動関数 (ψi , ψf ) は (偶関数, 奇関数) あるいは (奇関数, 偶関数) でなければならない。ここでは詳しく述べな
いが、この偶奇性 (パリティ) による制約の他に、Wigner―Eckart の定理等から遷移前後での軌道角運動量の変化 ∆Lが ±1 であることが導かれる。スピン角運動量 s も含めた条件 (選択則) は遷移前後での軌道角運動量の変化∆l = ±1

遷移前後での全角運動量の変化∆j = 0, ±1

遷移前後でのスピン角運動量の変化

54∆sz = 0

《講義ノート》

図 2.2.5

図 2.2.6

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

F 原子に X 線を入射し K 殻電子が M 殻へ遷移した場合を仮定した場合の蛍光 X 線 (c)-(d) と燐光 X 線 (e)-(f) の放射過程。ただし、
蛍光 X 線は Kα2 の場合。

(a)(b) 入射 X 線により K 殻電子が M 殻へ励起される。(c)M 殻に励起された電子ではなく、別の内殻 (L 殻) の電子が K 殻へ遷移し、その遷移エネルギーが外殻の M 殻電子へ。(d) エネルギーを得た外殻電子が Auger 電子として放出される。

55

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）となる。ただし、上記は原則であり、細かな禁則条件も存在するので注意が必要である。電気双極子モーメントを遷移行列要素とする遷移が電気双極子 (E1) 遷移である。高次の多重極子モーメントを遷移行列要素とする遷移も存在し、電気四極子モーメントを遷移行列要素とする遷移が電気双四極子 (E2) 遷移と呼ばれる。一般に、高次になるほど遷移行
列要素は小さくなる。E2 遷移の選択則は、

∆l

= ±2

∆j

= 0, ±1, ±2

∆sz

= 0

となる。こちらも細かな禁則条件に注意が必要である。

4f 電子系の希土類イオンでは、最外殻の 4f 軌道の内側に閉殻の 5d 軌道が存在する。そのため、希土類イオンの L
吸収端の E1 遷移と E2 遷移は表. 2 のようになる。E1 遷移と E2 遷移のエネルギーは数 eV 離れているが、私は E1 遷移と E2 遷移を別々に記載している文献に心当たりはない。基本的に、E2 遷移の遷移確率は、E1 遷移のものよりも小さい。そのため、吸収端エネルギーは E1 遷移を表記しているものと思われる。L 吸収端エネルギーは表. 2 のように

LIII ,LII ,LI で記載されている。
表 2 希土類イオンの L 吸収端エネルギーと電気双極子 (E1) 遷移, 電気四極子 (E2) 遷移

E1 遷移
E2 遷移

LIII
2p3/2 ↔(5d3/2 , 5d5/2 )2p3/2 ↔(4f5/2 , 4f7/2 )LII
2p1/2 ↔5d3/2
2p1/2 ↔4f5/2

LI
2s1/2 ↔(5p1/2 , 5p3/2 )2s1/2 ↔(5d3/2 , 5d5/2 )La
Ce

5.48275.72345.8906
6.1642

6.2663
6.5488Pr
Nd
Pm

5.96436.20796.45936.44046.72157.01286.83487.12607.4279SmEu

6.7162
6.97697.31187.61717.73688.0520Gd
Tb
Dy

7.24287.51407.79017.93038.25168.58068.37568.70809.0458HoEr

8.07118.35798.91789.26439.39429.7513Tm
Yb
Lu

8.64808.94369.24419.61699.9782
10.3486

10.1157
10.4864
10.8704
単位 [keV]

2.4

1 つの原子による散乱と原子散乱因子

2.4.1

Thomson 散乱の偏光因子

完全な直線偏光状態の X 線を原子中の一つの電子に入射すると、Thomson 散乱波は式 (2.26) のような直線偏光となること示した。まずは、この Thomson 散乱の微分散乱断面積と偏光について整理する。直線偏光の電場は式 (1.3) のように表せるので、Thomson 散乱の前後で位相が π ずれることと弾性散乱であることに注意すると、Thomson 散乱における入射電場 Ein と散乱位置からの距離 R での散乱電場 Erad の比は以下のように書ける。

56

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

Erad (r, t)
Einexp(ikR)sinθR
exp(ikR) ˆ′= −re|ε · ε̂|R= −re(2.31)(2.32)

ここで k は入射 X 線の波数ベクトル、ε̂, εˆ′ は散乱前後の単位偏光ベクトルである。
したがって、微分散乱断面積は式 (2.3) をもちいてdσdΩ= re2 |εˆ′ · ε̂|2(2.33)となる。この式 (2.33) は、入射偏光と散乱偏光のなす角 θ で表わした式 (2.29) と同じことを示している。式 (2.33) で表わすことにより、Thomson 散乱の散乱断面積は偏光に依存しない項 (r2 ) と偏光状態の項 |εˆ′ · ε̂|2 に分離できることeがわかる。偏光に依存しない re2 は電子による散乱波の基本単位 (１つの電子による散乱強度) である。また、re は電子の Thomson 散乱長と呼ばれる「１つの電子による散乱振幅」である。式 (2.31) で表わされるように、re は偏光による散乱振幅の減衰 (偏光因子) を除いた散乱振幅となる。

2.4.2

1 つの原子による散乱 ; Thomson 散乱と原子散乱因子

さて、式 (2.31) では、入射電場と散乱電場の比として表わしたため、散乱波の波数ベクトル k′ が現れていなかった。
これは、電子一つの散乱を考えていたためである。(図. 2.4.1(a))
ここからは、電荷密度 ρ(r) で分布した複数個 (Z 個) の電子で構成された単一原子による X 線散乱を考える。図. 2.4.1(b)にその模式図を示す。入射 X 線は、電荷分布中の点 A や点 B で等の各所で散乱される。したがって、散乱波はこの電荷分布中の各所の体積要素からの散乱波の重ね合わせとなる。また、X 線の電子による散乱は、主に弾性散乱である

Thomson 散乱である。そのため、入射 X 線の波数ベクトル k と散乱 X 線の波数ベクトル k′ の大きさは変化しない(|k′ |=|k|)。次に、電荷分布中の点 A と点 B が距離 r 離れているとして、それらの散乱波の位相差 ∆δ(r) を考える。光路差を考えると、位相差 ∆δ(r) は∆δ(r) = (k′ − k) · r ≡ Q · r(2.34)となる。
電荷分布を考慮した原子による散乱で k と k′ を必要とするのは干渉を考慮し、この位相差が現れるためなのである。
逆を言えば、式 (2.31) で散乱波の波数ベクトル k′ を必要としないのは電荷一つからの散乱では重ね合わせの概念がな
く、位相差を考える必要がないからということである。さらに、k と k′ よりも Q を重要視するのは、r 離れた点の位相差を決めるのは、入射 X 線や散乱波の個別の波数ベクトルと r の関係ではなく、常に Q と r の関係だからである。
ここで定義した Q は、散乱ベクトルと呼ばれ、図. 2.4.1(b) から想像できるように、弾性散乱では Q=|Q|=2|k|sinθ=(4π/λ)sinθ
である。
最後に、原子による散乱長を考える。ただし、ここでの散乱長は偏光による影響を除く。各体積要素からの散乱波の重ね合わせ (干渉) を考慮しなければ、位置 r における微少体積要素 dr の散乱振幅は Thomson 散乱長 −re に電子数

(ρ(r)dr) を掛けるだけで良い。干渉を考慮すると、位相因子として exp(−iQ · r) をさらに加えるので、原子の全散乱
長は∫−re ρ(r)exp(−iQ · r)dr ≡ −re f0 (Q)(2.35)となる。
このように、原子の全散乱長は電子一つからの Thomson 散乱長 (−re ) とその係数に分解できる。係数部分の f0 (Q)は、原子散乱因子あるいは原子形状因子と呼ばれ、電荷密度の異なる個々の原子に固有の値である。基準となる電子一つからの Thomson 散乱長 (−re ) は、波数ベクトル k′ が意味をもたないため Q 依存性をもたないが、原子一つからの全散乱長 (−re f0 (Q)) は図. 2.4.2(a) のような Q 依存性を示す。International Tables for Crystallography Vol. C には、

57

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

図 2.4.1 (a) 電子一つによる散乱と (b) 原子一つによる散乱, (c) 結晶による回折原子の電子密度が球対称を持っているとみなせる場合の原子散乱因子が (sinθ) / λ(=Q/4π) の関数として掲載されている。
式 (2.35) と後述の結晶構造因子は式の形は似通っているが、原子散乱因子は連続的な電荷分布での干渉であり、Q 依
存性にピーク構造をもたない (図. 2.4.2(a))。一方、連続的な分布ではなく、原子の周期的な配列による結晶構造因子のQ 依存性は、図. 2.4.2(c) のようにピーク構造をもつ。図. 2.4.2(c) の横軸を 2θ に変換し縦軸を結晶構造因子の 2 乗にとると、見慣れた X 線回折パターンに見える。ただし、当然ではあるが、実験で得られる回折強度は単純な結晶構造
因子の 2 乗に比例せず、結晶構造因子に偏光因子やローレンツ因子等の補正項が加わっている。
偏光因子は、第 2.4.1 節で示したように εˆ′ · ε̂ であり、偏光因子を取り入れた原子散乱因子は−re f0 (Q) (εˆ′ · ε̂)(2.36)

である。

2.4.3

原子散乱因子の分散補正項

吸収を考える上でも触れたが、原子内の電子は量子力学的な離散的エネルギー準位をもち、原子核に束縛されている。

K 殻電子のような内殻電子ほど強く束縛されるため、入射 X 線による電子振動が抑制された結果、吸収端近傍では原子散乱因子も小さくなる。すなわち、原子核による束縛を考慮した原子散乱因子では、吸収端近傍で式 (2.35)(b) の f0 にはないエネルギー (h̄ω) 依存性が現れるということである。その原子散乱因子の減少分は、分散補正項 f (h̄ω)+if (h̄ω)と呼ばれる複素数の形を取る。

58

《講義ノート》

図 2.4.2

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

散乱因子の Q 依存性。(a)Ce の原子散乱因子 (f0 (Q))；原子一つからの全散乱長は −re f0 (Q) である。散乱体である電子は密度 ρ(r) をもち、連続している。[6] (c) 結晶構造因子の模式図；原子散乱因子と異なり、散乱体は距離 R の間隔で周期的に配列している。そのため、Q 依存性は周期構造のフーリエ変換で不連続になる。

電子のうける束縛を考慮すると、式 (2.21) は式 (2.26) のように ω0 → 0 とみなすことはできない。そのため、式 (2.21) はEθ

e2
ω 2 Ein cos(ωt + δ)sinθ
2
4πε0 c me r me [(ω02 − ω 2 )2 + 4ρ2 ω 2 ] 21

=

−

=

−re

=
=
=
≃

ω2
1
[(ω02 − ω 2 )2 + 4ρ2 ω 2 ] 2
2sinθ
Ein cos(ωt + δ)
r

ω
sinθ
Ein cos(ωt + δ)1[(ω02 − ω 2 )2 + 4ρ2 ω 2 ] 2 r
ω2
sinθ
−re 2
Ein cos(ωt + δ)|(ω0 − ω 2 ) + 2iρω| r)(sinθ
ω 2 + 2iρω
Ein cos(ωt + δ)
−re 1 + 2 0 2(ω − ω0 ) + 2iρω
r
)
(
sinθ
ω02
−re 1 + 2
Ein cos(ωt + δ)2
(ω − ω0 ) − 2iρωr−re(2.37)となり、式 (2.31)-式 (2.35) と同様に偏光因子も考慮した原子の全散乱長は()∫
ω02−re 1 + 2
ρ(r)exp(−iQ · r)dr(ω − ω02 ) + 2iρω)∫(2ω02 ρωω 2 (ω 2 − ω 2 )= −re 1 + 2 0 2 2 0 2 2 + i 2
ρ(r)exp(−iQ · r)dr(ω − ω0 ) + 4ρ ω(ω − ω02 )2 + 4ρ2 ω 2≡ −re (f0 (Q) + f (h̄ω) + if (h̄ω))(2.38)となる。
また、原子核による束縛を考慮したのが分散補正項であるということは、分散補正項の散乱体は束縛された内殻電子である。束縛された内殻電子は、電子密度分布全体ではなく、原点に局在した状態であるため、原子の電子密度分布全体からの散乱である f0 とは異なり散乱角依存性はほとんどなく、エネルギー依存性だけとみなせる。そのため、Q 依
存性に加えて、h̄ω 依存性も考慮した原子散乱因子 f 0 は

f 0 (Q, h̄ω) = f0 (Q) + f (h̄ω) + if (h̄ω)(2.39)となり、f (h̄ω)+if (h̄ω) のように h̄ω の関数として表記される。
例として、図. 2.4.3 に Yb におけるエネルギー依存性を示す。各原子の分散補正項は http://skuld.bmsc.washington.edu/scatter/から得た。実験データの解析に分散補正項を用いる場合には、実験で得た吸収スペクトルを用いて分散補正項を得る。
これは分散補正項の虚部 f ” が吸収断面積に比例することを利用しており、実部 f ′ は虚部 f ” との Kramers-Kronig の関係式を用いて求めている。

59

《講義ノート》

図 2.4.3

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

Yb における分散補正項のエネルギー依存性。f ” は実験から得た吸収スペクトルから得ており、f ′ は f ” の Kramers-Kronig 変換で求めている。吸収端近傍の分散補正項は複雑であり、一般には測定試料毎に実験から求められる。sasaki table は東工大の佐々木先生が公開されていた分散補正項のエネルギー依存性である。出典元は International Tables for Crystallography Vol. C と思われる。

60

《講義ノート》

2.5

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

回折

最後に我々が観測する結晶からの散乱因子を考える。結晶は、図. 2.4.1(c) のように単位格子が規則的に並んだ状態である。
まず、単位格子内の j 番目の原子からの散乱因子を fj0 とすると、単位格子全体からの散乱因子 Fmol は、∑Fmol (Q, h̄ω) =fj0 (Q, h̄ω)exp(iQ · r j )(2.40)jとなる。
結晶からの散乱因子 Fcry は、単位格子を基本とした格子ベクトル R を用いて、Fcry (Q, h̄ω) = Fmol (Q, h̄ω)∑exp(iQ · Rn )(2.41)nである。ここで、格子ベクトル R は格子の基本ベクトル (a1 , a2 , a3 ) と整数 (n1 , n2 , n3 ) をもちいて

R n = n1 a 1 + n2 a 2 + n3 a 3(2.42)

で表わされる。
ここで、逆格子の単位ベクトル (a∗1 , a∗2 , a∗3 ) と格子点 (h, k, l) を用いた逆格子ベクトル G を考えると、

G = ha∗1 + ka∗2 + la∗3(2.43)

G · Rn = 2π(hn1 + kn2 + ln3 ) = 2π × 整数

(2.44)である。

であり、Q=G のときに、式 (2.41) の格子ベクトルに関する exp(iQ·Rn ) の位相が揃うことがわかる。また、exp(iQ·Rn) は G · Rn = 2π × 整数 より ±1 になる。
したがって、結晶構造因子の絶対値は|Fcry (Q, h̄ω)|= | ± Fmol (Q, h̄ω)|∑= |fj0 (Q, h̄ω)exp(iG · r j )|
j

|

=

∑

fj0 (Q, h̄ω)exp(2πi(hxj + kyj + lzj ))|(2.45)

jと表記される。
実際、逆格子ベクトル G(h, k, l) における弾性散乱では、

k′ − kG =|k|= 2π/λ(2.46)

より、

k′
k

′2

G2

= (G + k)=

G2 + 2G · k + k2

= 2G · k= 2kcos(
= 2

2
π
− θ)
2

2
2π
sin(θ)
λ61(2.47)

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

ここで、格子面間隔 d と逆格子ベクトル G の関係は, G=2nπ/d (n; 整数) であるから

2nπ
2π= 2 sin(θ)
dλnλ = 2dsinθ

(2.48)であり、Q=G は Bragg の条件を満たす。
原子散乱因子や結晶構造因子には、式 (2.33) で現れる偏光因子は省略されている。偏光因子も含めて、結晶からの微分散乱断面積を求めると、dσdΩ= | − re Fcry (Q, h̄ω)εˆ′ · ε̂|2
= re2 |Fcry (Q, h̄ω)εˆ′ · ε̂|2∑
= re2 |(εˆ′ · ε̂
)fj0 (Q, h̄ω)exp(iG · r j )||2
j

となる。

62

(2.49)

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

量子力学的な記述における電子による X 線の散乱

3

前章では、X 線の電子による散乱は振動する電場によって電子を振動させることであり、散乱波は電子の振動による電磁波の放射として扱った。これは電磁波と電子の相互作用を古典的に扱った模型であり、Thomson 散乱と呼ばれる散乱であった。
この章では、X 線の電子による散乱を電子と電磁場との相互作用ハミルトニアンから考え、Thomson 散乱だけでなく非共鳴 X 線磁気散乱と共鳴線散乱の３つの散乱振幅を求める。共鳴 X 線実験では、散乱ベクトルと散乱前後の偏光をパラメーターとして得た散乱強度から秩序変数 (多重極子モーメント) を議論する。したがって、最終的に散乱振幅を散乱ベクトルと散乱前後の偏光, 多重極子演算子の３つで整理する。

3.1

電子と電磁場との相互作用ハミルトニアン

一つの原子に X 線を入射した際に、原子に束縛されている電子が電磁場を受ける状況下では、電磁場を含める系全体のハミルトニアン H は以下のように表わされる。

H=)2 ∑ Ze2 ∑
∑ 1 (
e
e2
eh̄ ∑
pi − A(r i ) −+−si · (∇ × A(r i ))2mcri|ri − rj | mc i
i
i
i,j
(
)
(
) ∑
e
1
eh̄ ∑
†
si · E(r i ) × pi − A(r i ) +h̄ωk aku aku +− 2 2
2m c i
c
2

(3.1)

k,u

第１項は電磁場中における電子の運動エネルギー、第２項は原子殻からのポテンシャル、第３項は電子間のクーロン相
互作用、第４項は電子スピンと磁場との相互作用、第５項は相対論的効果によるスピン軌道相互作用、第 6 項は電磁場のエネルギーである。a†ku は波数ベクトル k, 偏光状態 u の光子の生成演算子であり、aku は消滅演算子である。また、

ωk は入射 X 線の振動数である。
ベクトルポテンシャル A を ∇ · A=0 となるようにとると、A と運動量 p(=-ih̄∇) との内積関係は
p · A = −ih̄∇ · A = 0

(3.2)

)2
e
e2
e
pi − A(r i ) = p2i − 2 A(r i ) · pi + 2 A(r i )2
c
c
c

(3.3)

である。
したがって、式 (3.1) の第１項は(となり、第１項が電子系のハミルトニアン、第 2,3 項が電子と電磁場との相互作用のハミルトニアンと分離できる。
また、Maxwell 方程式から定義される電場と ψ, A の関係は

E

=

−∇ψ −

1 ∂
A
c ∂t

(3.4)

より式 (3.1) の第５項はeh̄ ∑esi · E(r i ) × (pi − A(r i ))2m2 c2 i
c
{
}
∑
e ∑
∂A(r i )eh̄
si · (−∇ψ(r i ) × pi ) + 2
si · (× A(r i ))2m2 c2
c i
∂t
i

−
≈

(3.5)

となる。Aの一次項は式 (3.1) の第１項の一次項と比較して h̄ω/mc2 が係るだけ小さく、無視している。式 (3.1) を電子系のハミルトニアン H0 , 電磁場のハミルトニアン HR , 電磁場と電子の相互作用のハミルトニアン H ′ の三つに分

63

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

類すると以下のように表わされる。

H

=

H0

=

HR

=

H0 + HR + H ′∑ p2∑ Ze2 ∑
∑
e2
i
−
+
+
ξ(r i )h̄2 si · li
2m
r
|r
−
r
|
i
i
j
i
i
i,j
i
∑
k,u

H′

=

≡

1
h̄ωk (a†ku aku + )
2

(3.6)
(3.7)
(3.8)

e2 ∑
e2 h̄ ∑
∂A(r i )× A(r i ))A(r i )2 −si · (
22mc i
2m2 c4 i
∂t
e ∑
eh̄ ∑
−
A(r i ) · pi −si · (∇ × A(r i ))
mc i
mc i

(3.9)

H1′ + H2′ + H3′ + H4′

ただし、ξ(r i ) = (e/2m2 c2 r)(dψ(r i )/dr) である。これらのハミルトニアンの内、散乱に寄与するハミルトニアンは電磁場と電子の相互作用の H ′ になる。以降、H ′ の 4 つの項を H1′ , H2′ , H3′ , H4′ と呼称する。
また、H ′ は入射電磁場のベクトルポテンシャル A の関数である。古典的な記述でも示したように、電場や磁場は Aを求めることにより得られる。電磁波の場合、時間に依存しないベクトルポテンシャル A の演算子表記は√∑ 2πh̄c2
A =(εku aku eik·r + ε∗ku a†ku e−ik·r )V ωk(3.10)k,uとなる。εku は光子の偏光ベクトル, V は系全体の体積である。

3.2

遷移確率と散乱断面積

量子力学的な散乱過程は、時間に依存する摂動論により記述される。入射 X 線と電子の系全体の始状態 |i⟩ から終状
態 |f ⟩ への遷移は、入射 X 線 (電磁場) と電子の相互作用ハミルトニアン H ′ により与えられる。二次の摂動論の Fermiの黄金則では、中間状態を |n⟩ として、|i⟩ から |f ⟩ への遷移確率 W は以下のように表わされる。

W=∑ ⟨f |H ′ |n⟩⟨n|H ′ |i⟩
2π⟨f |H ′ |i⟩ +h̄Ei − Enn2

δ(Ei − Ef )(3.11)

W は単位時間当りの遷移確率なので、終状態の状態密度を ρ(Ef ) とすると、(単位時間当りの散乱 X 線の光子数)=W ρ(Ef ) である。したがって、散乱 X 線強度 Irad は検出器の立体角 ∆Ω に入る単位時間当りの散乱 X 線の光子数であるから、

Irad = W ρ(Ef )∆Ω, ρ(Ef ) =V ωk2′(2π)3 h̄c3

(3.12)となる。散乱断面積は入射 X 線のフラックス Φin を用いて以下のように表わせる。

dσ
dΩ

=

Irad
W ρ(Ef )
d2 σ
=
≡
∆ΩΦin
Φin
dΩdE

(3.13)

遷移エネルギーも関係してくるため、散乱断面積当りのみならずエネルギー当りも考慮しなければならないため、微分
散乱断面積を d2 σ/dΩdE と表記した。実験では Φin , ρ(Ef ) は一定であるので、遷移確率は微分散乱断面積に比例する。
式 (3.11) に戻り、遷移確率の内容を考える。式 (3.11) の第１項が一次の摂動項、第２項が二次の摂動項となっている。まずは、A に対して一次である H3′ を例に一次摂動, 二次摂動を考える。A の表式は式 (3.10) であるので、A の一次である H3′ の一次摂動は「|i⟩ から |f ⟩ への遷移では、電子状態 |a⟩ は変わらず、aku あるいは a†ku による光子の消滅あるいは生成が生じる」ことを示す。すなわち、これは光電吸収 (内殻電子の入射 X 線の光子エネルギーの吸収によ64

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

る高エネルギー準位の外殻への遷移と光子の消滅) である。ここで、|a⟩ を電子の始状態とし、X 線と電子の系全体の始状態 |i⟩ と区別する。) 一方、H3′ の二次摂動は電子の中間状態状態 |b⟩ を経ている。弾性散乱の場合、電子の始状態と終状態は同じ |a⟩ であるので単純に解釈すると、⟨n|H3′ |i⟩ では入射 X 線の吸収による電子状態の中間状態への遷移と入射光子の消失, ⟨f |H3′ |n⟩ では入射光子と同じエネルギーの散乱光子の生成と電子状態の中間状態から始状態 |a⟩ への還元である。しかし、中間状態への遷移は、入射 X 線光子が消滅する場合のみならず散乱 X 線光子が生成された場合もあり得る。前者の場合は中間状態から終状態への遷移で散乱 X 線光子を生成し、後者の場合は入射 X 線光子が消滅する。そのため、中間状態を経由には 2 つの経路があり得る。
何れにしても、一次の A のでは、吸収過程である一次摂動 (光子の消失のみ) は散乱を与えず、二次摂動は中間状態を経た吸収と散乱 (光子の消失と生成) を与える。
次に A に対する二次の H1′ の一次摂動を考える。二次の A の一次摂動は「|i⟩ から |f ⟩ への一次遷移でも、電子状態|a⟩ は変わらず、aku a†ku による散乱 (光子の消滅と生成) が生じる」ことを示す。これは、電子のエネルギー準位を基底状態に保ったまま散乱する Thomson 散乱である。

したがって、二次の A は一次摂動で散乱を与えることに比べて、一次の A は二次摂動でなければ散乱を与えず、二次摂動でのみ散乱を与える。また、H2′ や H4′ はスピンsも含み、前者は一次摂動で与えられる電子スピンによる磁気散乱項、後者は二次摂動で与えられる磁気散乱項であることが判る。

H1′ , H2′ , H3′ , H4′ を全て含んだ H ′ の微分散乱断面積 (|a, ku⟩ → |a′ , k′ u′ ⟩) は以下のようになる。
d2 σ
dΩdE a,ku→a′ ,k′ u′

(
=

e2
mc2

)2

1 ∑
+
m
′
b,j,j

⟨a′ |
{

∑

e−iκ·r j |a⟩(ε′ · ε) − ijh̄ωk ′ ∑⟨a |sj e−iκ·r j |a⟩ · (ε′ × ε)mc2
j′⟨a′ |(ε′ · pj ′ − ih̄sj ′ · k′ × ε′ )e−ik ·r j ′ |b⟩⟨b|(ε · pj + ih̄sj · k × ε)eik·r j |a⟩
Ea − Eb + h̄ωk + iΓb /2′⟨a′ |(ε · pj + ih̄sj · k × ε)eik·r j |b⟩⟨b|(ε′ · pj ′ − ih̄sj ′ · k′ × ε′ )e−ik ·r j′ |a⟩+Ea − Eb − h̄ωk′ − iΓb /2}2

×δ(Ea − Ea′ + h̄ωk − h̄ωk′ )(3.14)

ここで、散乱ベクトルを κ = k′ − k と表わした。|a′ ⟩, |a′ ⟩ は電子の始状態と終状態。弾性散乱の場合は散乱の前後で電子状態に変化はないので |a⟩=|a′ ⟩。|b′ ⟩ は電子の中間状態。二次摂動項である {} 内の第 1 項と第 2 項は、それぞれ入射光子が先に消滅する場合と散乱光子が先に生成される場合を示し、始状態と中間状態の電子エネルギーを Ea , Eb とす
ると、それぞれの中間状態のエネルギー En は前者が En =Eb 、後者が En =Eb + h̄ωk + h̄ωk′ 。始状態のエネルギーを Ei= Ea + h̄ωk であり、二次摂動項の分母は始状態と中間状態の差になる。Γb は、中間状態の寿命に対応するエネルギー
幅として導入している。この分母は、二次摂動が生じるのは”基本的に”入射 X 線の光子エネルギー h̄ωk が |Ea − Eb | を中心とした Γ の幅の中である時のみであることを示す。このように、h̄ωk と遷移エネルギーが一致したときにのみ現れる散乱を共鳴散乱、h̄ωk =|Ea − Eb | を共鳴条件と呼ぶ。また、共鳴エネルギーとは遷移エネルギー (|Ea − Eb |) を指す。
一方、共鳴条件以外でも現れる一次摂動ような散乱を非共鳴散乱とよぶ。勘違いしないで頂きたいのは、非共鳴散乱は共鳴条件下でも現れる。ただし、軌道磁気モーメントによる非共鳴磁気散乱項は二次摂動で表わされ、式 (3.14) の第３項の中に紛れ込んでいる。そのため、二次摂動項=共鳴散乱項とならないことに注意しなければならない。
したがって、式 (3.14) の || 内は、結局第１項は H1′ の一次摂動である非共鳴散乱 (トムソン散乱)、第２項は H2′ の一次摂動である電子スピンによる非共鳴磁気散乱, 第３項 ({} 内) は H3′ と H4′ の二次摂動である共鳴磁気散乱 (+軌道
磁気モーメントによる非共鳴磁気散乱) である。
次節では第３項を共鳴散乱項と非共鳴散乱項に分離し、散乱振幅を共鳴散乱項と非共鳴散乱項で整理する。

65

《講義ノート》

3.3

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

散乱振幅の共鳴散乱項と非共鳴散乱項式 (3.14) に電子の運動密度を表す演算子

J (k) =∑

(pj − ih̄k × sj )eik·r j(3.15)

jを導入すると、散乱振幅 F は以下のように表わされる。

F

= −




∑
e
h̄ωk ′ ∑′−iκ·r j′⟨a |sj e−iκ·r j |a⟩ · (ε′ × ε)
⟨a
|
e
|a⟩(ε
·
ε)
−
i
2
mc2 
mc
j
j
2

}
1 ∑ ⟨a′ |ε′ · J † (k′ )|b⟩⟨b|ε · J (k)|a⟩ ⟨a′ |ε · J (k)|b⟩⟨b|ε′ · J † (k′ )|a⟩
+
+
)
(
m
Ea − Eb + h̄ωk + iΓb /2
Ea − Eb − h̄ωk′ − iΓb /2(3.16)b

この散乱振幅をクロネッカーのデルタ δαβ と 3 階反対称テンソル εαβγ を用いると以下のように表記できる。


2 ∑
∑
∑
e
h̄ωk
sj e−iκ·r j |a⟩εαβγ
ε′α εβ ⟨a′ |F = − 2e−iκ·r j |a⟩δαβ − i 2 ⟨a′ |

mc
mc
jγ
j
α,β
}
1 ∑ ⟨a′ |J †α (k′ )|b⟩⟨b|J β (k)|a⟩ ⟨a′ |J β (k)|b⟩⟨b|J †α (k′ )|a⟩+(
(3.17)+)
m
Ea − Eb + h̄ωk + iΓb /2
Ea − Eb − h̄ωk′ − iΓb /2b

ここで、αβ は電子遷移の自由度に対応する。特に二次摂動項では E1 遷移に対しては ∆l = ±1 なので、立法調和関√√√√数の α, β=x, y, z に対応する。E2 遷移でに対しては α, β=(3z 2 − r2 )/2, 3(x2 − y 2 ), 3yz, 3zx, 3xy である。 γ は電子スピンの自由度で ±1 に対応すると思われる。
さて、式 (3.14) の後に述べたように、軌道磁気モーメントによる非共鳴磁気散乱項は二次摂動にあるため、これではまだ共鳴成分と非共鳴成分の分離ができていない。すなわち式 (3.17) の第３項にある軌道磁気モーメントによる非共鳴磁気散乱項を抽出する必要がある。非共鳴磁気散乱項が存在するということは、Eb − Ea ≫ h̄ωk あるいは Eb − Ea ≪ h̄ωk
であっても 0 でない成分が存在するということであり、これを抽出すれば良い。
二次摂動項の分母部分に着目し、Eb − Ea ≫ Γb かつ Eb − Ea について Γb を無視すると以下の変形することができる。ここで、遷移エネルギーを Eb − Ea = h̄ωr と表記すると

1
Ea − Eb + h̄ωk + iΓb /2

1
Ea − Eb − h̄ωk − iΓb /2

1
1
1
−
+
h̄ωk − h̄ωr + iΓb /2 h̄ωk
h̄ωk
h̄ωr
1
1
≈ (
)
+
h̄ωk h̄ωk − h̄ωr + iΓb /2 h̄ωk
1
1
1
=
+
−
−(h̄ωk − h̄ωr + iΓb /2) h̄ωk
h̄ωk
h̄ωr
1
1
≈ (
)
−
h̄ωk h̄ωk − h̄ωr + iΓb /2 h̄ωk=

いまさらであるが、弾性散乱を考えているので、h̄ωk = h̄ωk′ とした。２式とも左辺より非共鳴項 (右辺第 2 項) を抽出することになる。上式を用いて式 (3.17) をもう一度表記すると


∑
∑
h̄ωk
e2 ∑ ′
εα εβ ⟨a′ |sj e−iκ·r j |a⟩εαβγ
F = − 2e−iκ·r j |a⟩δαβ − i 2 ⟨a′ |

mc
mc
α,β

jγ

j

1 1
⟨a′ |J †α (k′ )J β (k) − J β (k)J †α (k′ )|a⟩+m h̄ωk)}()( ′ † ′⟨a |J α (k )|b⟩⟨b|J β (k)|a⟩ + ⟨a′ |J β (k)|b⟩⟨b|J †α (k′ )|a⟩1 ∑ h̄ωr+mh̄ωkh̄ωk − h̄ωr + iΓb /2
b66(3.18)

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

抽出した第３項に関しては、|b⟩ についての closure relation が適用でき、∑

n ⟨f |B|n⟩⟨n|A|i⟩ = ⟨f |BA|i⟩ となる。

式 (??) の第２項 (電子スピンによる非共鳴磁気散乱項) と第３項 (軌道磁気モーメントによる非共鳴磁気散乱項) をまとめると、散乱振幅は以下のように表わせる。


2 ∑
∑
e
F = − 2ε′α εβ ⟨a′ |e−iκ·r j |a⟩δαβ

mc
j

α,β

h̄ωk ′ ∑ −iκ·r j (pj × κ)γ αβγ
⟨a |
e
{i
A
+ sjγ B αβγ }|a⟩
2
mc2
h̄κ
j
)}
)( ′ † ′(∑1⟨a |J α (k )|b⟩⟨b|J β (k)|a⟩ + ⟨a′ |J β (k)|b⟩⟨b|J †α (k′ )|a⟩
h̄ωr
+
m
h̄ωk
h̄ωk − h̄ωr + iΓb /2−i(3.19)

b

ただし、

Aαβγ

= −2(1 − k̂ · kˆ′ )εαβγ

B αβγ= εαβγ − εαδγ k̂δ′ k̂β′ + εβδγ k̂δ′ k̂α′11 αβδ ′1ε (k̂δ k̂γ + k̂δ k̂γ′ ) + (k̂ × kˆ′ )α δβγ + (k̂ × kˆ′ )β δαγ−222k̂, kˆ′ は入射 X 線と散乱 X 線の波数ベクトルの単位ベクトルである。
最後に、δαβ や εαβγ ではなく、偏光ベクトルを用いて式 (??) を表わすと
F= (Fcharge + Fmag ) + Freso=

Fcharge

=

Fmag

=

Freso

=

Fnonreso + Freso∑e2e−iκ·r j |a⟩(ε′ · ε)
− 2 ⟨a′ |mc
j(3.20)∑(pj × κ)
h̄ωk
e2
′
(−i
)⟨a
|
· B + sj · D}|a⟩e−iκ·r j {i
mc2
mc2
h̄κ2
j
)
(
)( ′ † ′⟨a |J α (k )|b⟩⟨b|J β (k)|a⟩ + ⟨a′ |J β (k)|b⟩⟨b|J †α (k′ )|a⟩e2 ∑ ′
1 ∑ h̄ωr
− 2
εα εβ
mc
m
h̄ωk
h̄ωk − h̄ωr + iΓb /2−α,β(3.21)(3.22)

b

ただし、

B

=

κ2
− 2 (ε′ × ε)kD={(ε′ × ε) − (ε · kˆ′ )(ε′ × kˆ′ ) + (ε′ · k̂)(ε × k̂) − (kˆ′ × ε′ ) × (k̂ × ε)}となる。

Freso は共鳴散乱項であり、入射エネルギー h̄ωk が共鳴エネルギー h̄ωr 近傍の時 (|h̄ωk − h̄ωr | < iΓb /2) のみ値をもち得る。一方、Fcharge と Fmag は非共鳴散乱項であり、h̄ωk が h̄ωr 近傍であろうとなかろうと値をもち得る。散乱振幅は、このように非共鳴散乱項 Fnonreso と共鳴散乱項 Freso に分離できる。
さて、これらの散乱振幅は電子の運動密度演算子 J (k) と電子位置 r), 散乱ベクトル (κ), 散乱前後の偏光 (ε, ε′ ) で書かれている。しかし、我々が秩序変数として議論する演算子は、運動密度演算子ではなく磁気双極子モーメントなどの多重極子演算子 (自由度) である。次節では、これら散乱項を多重極子演算子と波数ベクトル, 散乱前後の偏光ベクトルで整理する。

67

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

4

非共鳴散乱振幅における多重極子演算子と波数ベクトル, 偏光因子

4.1

Thomson 散乱：電荷 (単極子モーメント) からの非共鳴散乱式 (3.21) は、第 3.2 節から導いたように A2 の一次摂動項である。また、古典的な散乱原理から導かれる Thomson
散乱項と等しく、以下にそれを示す。
電荷を表わす電気単極子モーメント q は以下のように電荷密度の実空間積分であり、電荷密度 ρ のフーリエ変換の関係よりq(κ) =∑

e−iκ·r j

j

=

∫

1
4πε0

(4.1)

ρ(r)e−iκ·r dv(4.2)

したがって、式 (3.20) は、|a⟩ を電荷分布の電子系の波動関数として

−

∑
e2
⟨a|
e−iκ·r j |a⟩ (ε′ · ε) =
2
mc
j

−

e2
⟨a|q(κ)|a⟩ (ε′ · ε)mc2e2
⟨a|q(κ)|a⟩ (ε′ · ε)mc2e2
= − 2 ⟨q⟩ (ε′ · ε)
mc
= −

あるいは

−

∑
e2
⟨a|
e−iκ·r j |a⟩ (ε′ · ε) =
2
mc
j

−

e2
mc2

∫

ψ∗

2

∑

∫

e−iκ·r j ψ dv (ε′ · ε)je1ρe−iκ·r dv (ε′ · ε)
mc2 4πε0e2 1− 2f0 (κ) (ε′ · ε)
mc 4πε0
−

=
=

(4.3)

(4.4)

である。
式 (4.4) より、古典的な散乱原理から導かれる「偏光因子を考慮した Thomson 散乱項 (式 (2.36))」と等しいことがわかる。
同時に、式 (4.3) と式 (4.4) より∫
f0 (κ) =ρe−iκ·r dv= 4πε0 ⟨a|q(κ)|a⟩(4.5)

であるから、κ=0 のとき f0 は電荷 (全電気単極子モーメント) であり、原子散乱因子 f0 (κ)(=4πε0 ⟨a|q(κ)|a⟩) は単極子モーメントの形状因子といえる。
最終的に、式 (3.21) はFcharge= −re 4πε0 ⟨q⟩ (ε′ · ε)

(4.6)

= (4πε0 × Thomson 散乱長) × (電気単極子モーメントの期待値 [電荷の大きさ]) × (偏光因子)
または= −re f0 (κ) (ε′ · ε)

(4.7)

= (Thomson 散乱長) × (電気単極子モーメントの形状因子 [原子散乱因子]) × (偏光因子)と表記できる実験では、逆格子点 (h, k, l) で実験するので回折角 2θ は与えられ、原子の種類より f0 (κ) あるいは ⟨q⟩ も与えられるので、Thomson 散乱の実験的パラメーターは偏光因子 (=ε′ · ε) のみとなる。

68

《講義ノート》

4.2

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

磁気モーメント (スピン磁気モーメント, 軌道磁気モーメント) からの非共鳴散乱：非共鳴磁気散乱式 (3.21) は、s · (∂ A(r i )× A(r i )) の 1 次摂動項であるスピン磁気モーメントに関係する非共鳴散乱項と A · p の 2∂t

次摂動項に紛れ込んでいる軌道磁気モーメントに関係する非共鳴散乱項である。
電荷と同様にスピン磁気モーメントを表わす双極子モーメント QS は以下のようにスピン磁気モーメント密度 mS の実空間積分であり、mS のフーリエ変換の関係より

QS (κ) =∑e−iκ·r j sj

j

= −

1
2µB

∫

(4.8)

mS (r)e−iκ·r dv

(4.9)

である。
電気単極子モーメントと原子散乱因子 (単極子モーメントの形状因子) との関係とくらべると、∫
fmagS (κ)(スピン磁気双極子モーメントの形状因子) =mS (r)e−iκ·r dv= −2µB ⟨a|QS (κ)|a⟩(4.10)といえる。
式 (3.21) の第二項において、D は電子状態とは関係がないので、式 (4.10) を用いて以下のように変形でき、スピン
磁気双極子モーメントの形状因子の項であることがわかる。

FmagS

= −

∑
e2
h̄ωk
′
(−i
)⟨a
|
e−iκ·r j sj · D|a⟩
mc2
mc2
j

h̄ωk
e2
(−i 2 )(−2µB )⟨a|QS (κ)|a⟩ · D
2
mc
mc
h̄ωk
e2
= − 2 (−i 2 )(−2µB )⟨QS ⟩ · D
mc
mc

= −

(4.11)

軌道磁気モーメントを表わす双極子モーメント QL も以下のように軌道磁気モーメント密度 mL の実空間積分であり、mL のフーリエ変換の関係より

QL (κ) =

−

1
µB

∫

mL (r)e−iκ·r dv

であり、軌道磁気双極子モーメントの形状因子は以下のようになる。

∫

fmagL (κ)(軌道磁気双極子モーメントの形状因子) =(4.12)mL (r)e−iκ·r dv= −µB ⟨a|QL (κ)|a⟩(4.13)

したがって、磁気形状因子は以下のようになり、κ=0 の場合がイオンの全磁気モーメントになる。

−µB (⟨a|QL (κ)|a⟩ + 2⟨a|QL (κ)|a⟩fmagL (κ) + fmagS (κ) =(4.14)

式 (3.21) の第 1 項は、軌道磁気双極子モーメントの形状因子と関係がないように見えるが、実は以下に示すように

κ に垂直な軌道双極子モーメントの形状因子成分の項である。

Q⊥L (κ) =∑e−iκ·r jj((pj × κ)ih̄κ2=κ̂ × (QL × κ̂)=QL − (QL · κ̂)κ̂

69)(4.15)

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

より

= −

Fmag

h̄ωk
e2
(−i 2 )⟨a|(Q⊥L · B + QS · D)|a⟩2
mc
mc(4.16)と変形できる。B, D は、式 (3.23), 式 (3.23) のままである。
また、Fmag は、Q⊥L ではなく QL を用いて、

Fmag

= −

e2
h̄ωk
(−i 2 )(⟨QL ⟩ · B ⊥ + ⟨QS ⟩ · D)
mc2
mc

(4.17)と変形できる。ここで、κ に垂直な B を B ⊥ としてB⊥= B − (B · κ̂)κ̂ = κ̂ × (B × κ̂)= 2(1 − k̂ · k̂′ ) {κ̂ × (ε′ × ε) × κ̂}}Q⊥L · B={κ̂ × (QL × κ̂)} · B={QL − (QL · κ̂)κ̂} · B=

QL · {B − (B · κ̂)κ̂}=

Q · B ⊥L

(4.18)(4.19)を用いた。
実験では、逆格子点 (h, k, l) で実験するので回折角 2θ は与えられ、原子の種類より f0 (κ) あるいは ⟨q⟩ も与えられる。そのため、非共鳴磁気散乱の実験的パラメーターは、スピン磁気モーメントと軌道磁気モーメントの大きさ, 偏光/
波数ベクトル因子 (=B, D) のみとなる。

4.3

非共鳴磁気散乱・共鳴散乱の回折実験配置

式 (4.17) に示したように、非共鳴磁気散乱項は磁気モーメントと X 線の偏光ベクトル、波数ベクトルを変数とする。
したがって、実験における非共鳴磁気散乱強度と X 線の偏光ベクトル及び波数ベクトルの条件を把握できれば、磁気モーメントを観測できるということである。
後述するように、共鳴散乱項も非共鳴磁気散乱項と同様に、多重極子モーメントと X 線の偏光ベクトル、波数ベク
トルを変数とする。そのため、実験方法としては共鳴散乱も非共鳴磁気散乱も変わらない。
その実験配置を図 (4.3.1) に示す。図右手側 (上流) より放射光 X 線が入射してくる。

入射 X 線
通常、入射されてくる放射光 X 線の偏光は水平偏光であり、下流の回折計が水平振りである場合には π 偏光 (偏光ベクトル ∥ 散乱面) である。

入射 X 線の偏光制御と偏光角度

π 偏光である入射 X 線の偏光を自在に制御するのが、移相子 (Phase Retarder) である。移相子により π 偏光を直線
偏光のみならず楕円偏光や円偏光へ変換することが可能である。
SPring-8 BL22XU では、ダイヤモンド結晶 2 枚を用いた移相子システムを構築している。ここでは詳しくは述べない
が、結晶を 2 枚用いることにより入射 X 線の角度分散による偏光度低下を抑制しており、入射 X 線の直線偏光度 99.9%に対して偏光制御後の線偏光度 97% 以上を維持している。
非共鳴磁気散乱・共鳴散乱の回折実験では直線偏光を用いており、σ 偏光に対する直線偏光角度を η と定義した。

70

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

図 4.3.1 非共鳴磁気散乱・共鳴散乱の回折実験配置

試料での回折
試料での回折は、Bragg 則 (λ = 2dsinθ) に従い、実験室系の回折と基本的に変わりはない。異なる点はアジマス角

Ψ と呼ばれる散乱ベクトル G 周りの回転角である。Ψ を回しても、散乱ベクトル G は変化しないため Bragg 則は満たされたままである。しかし、結晶内の磁気モーメントや電気四極子モーメントなどの双極子以上の多重極子モーメントは回転するため、多重極子モーメントと偏光/波数ベクトルの関係が変化し、非共鳴磁気散乱では式 (4.17) に従って散乱強度依存性が現れる。この強度依存性がアジマス角依存性と呼ばれる。
一方で、式 (4.7) に従う Thomson 散乱では、多重極子モーメントは単極子モーメントでありスカラー量である。そのため、アジマス角を回転させても単極子モーメントと偏光ベクトルの関係は変化せず、アジマス角依存性も現れない。
また、磁場印加方向は常に散乱面に垂直であり、アジマス角を回転させると磁場印加方向も変化する。

回折 X 線の偏光角度
直線偏光を用いた非共鳴磁気散乱・共鳴散乱の回折 X 線の偏光は直線偏光のままであるが、その偏光角度は変化し得る。回折後の偏光角度を ξ と定義した。

回折 X 線の偏光解析
回折 X 線の偏光は、高純度単結晶をアナライザー結晶とした偏光解析装置により解析する。原理的には Thomson 散
乱を利用しており、実験のエネルギーにおいてアナライザー結晶の Bragg 角 2θA が 90◦ となるような格子面間隔をもつアナライザー結晶を選択すると、アナライザー結晶の散乱面に垂直な偏光のみを取り出すことができる。

ϕA は試料の散乱面からの角度であり、ϕA =90◦ のとき、回折 X 線の π ′ (ξ=90◦ ) 偏光成分のみを取り出せる。また、
ϕA =0◦ のとき、回折 X 線の σ ′ (ξ=0◦ ) 偏光成分のみを取り出せる。

71

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

文章にすると、すこし混乱しそうになるが、ϕA =90◦ のとき、回折 X 線の π ′ 偏光成分はアナライザーにとって σ 偏
光成分である。Thomson 散乱の偏光因子 ε′ · ε では、π 偏光成分に掛かる sin2θA が σ 偏光成分には掛からない。そのため、2θA =90◦ では π 偏光成分が 0 となり、σ 偏光成分のみを取り出すことが出来る。ただし、2θA が 90◦ から 5◦ ず
れると π 偏光成分は 2% 程度混入してしまうので、一つのアナライザー結晶の反射指数では 500eV 程度しか利用できるエネルギー範囲がない。
偏光解析装置についてまとめると、ϕA のとき回折 X 線の ξ=ϕA の偏光成分のみを取り出せるということである。

4.4

非共鳴磁気散乱の回折実験シミュレーション式 (4.17) と図 (4.3.1) を用いて、非共鳴磁気散乱の回折実験での実験結果を計算する。
モデルとしては、以下の仮定で行う。

• 300 反射が禁制反射である結晶• 磁気秩序の波数ベクトル kmag =(1,0,0)
計算する磁気反射指数 300 反射(散乱ベクトル G=(3, 0, 0))• ⟨QS ⟩(x, y, z) = (0, 0, Sz )• ⟨QL ⟩(x, y, z) = (0, Ly , 0)
• Ψ=0◦ の時、実験室系の Z 軸は結晶系の c 軸と一致

4.4.1

非共鳴磁気散乱のアジマス角依存性 (π と σ)

第 4.3 節のアジマス角 Ψ の説明で触れたように、磁気モーメントや電気四極子モーメントなどの双極子以上の多重極子モーメントによる回折強度はアジマス角依存性をもち得る。実際に先に挙げたモデルでの確認する。式 (4.17) に従って、散乱振幅を計算する場合、結晶系 xyz 座標 (abc 軸) で表現し易い多重極子モーメントと実験室系 XYZ 座標で表現し易い波数ベクトルや偏光ベクトルを同じ座標系でまとめる必要がある。偏光状態を考える上では、実験室系 XYZ 座
標でまとめるのが安易であろう。
そこで、実験室系 XYZ で表現された波数ベクトルや偏光ベクトルは下記のようになる。

k = k(sinθ, cosθ, 0),k′ = k(−sinθ, cosθ, 0)

εσ = (0, 0, 1),ε′σ = (0, 0, 1)

επ = (−cosθ, sinθ, 0),ε′π = (−cosθ, −sinθ, 0)(4.20)

一方、結晶系 xyz 座標 (abc 軸) で表現されたスピン/軌道磁気モーメントは⟨QS ⟩ (x, y, z) = (0, 0, Sz )(4.21)⟨QL ⟩ (x, y, z) = (0, Ly , 0)(4.22)である。

72

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

Ψ=0◦ と 90◦Ψ=0◦ では、スピン/軌道磁気モーメントを XYZ 軸で表現するのは容易であり、⟨QS ⟩ (X, Y, Z)Ψ=0◦ = (0, 0, Sz )(4.23)⟨QL ⟩ (X, Y, Z)Ψ=0◦ = (0, −Ly , 0)(4.24)⟨QS ⟩ (X, Y, Z)Ψ=90◦ = (0, Sz , 0)(4.25)⟨QL ⟩ (X, Y, Z)Ψ=90◦ = (0, 0, Ly )(4.26)である。

Ψ=90◦ でも同様に、である。

σσ ′ , σπ ′ , πσ ′ , ππ ′ 散乱過程における散乱振幅を 2×2 行列で表現するとe2h̄ωk(−i 2 )(⟨QL ⟩ · B ⊥ + ⟨QS ⟩ · D)
mc2
mc (
)
′
′
Fσσ
Fπσ
e2
h̄ωk
′
′
= − 2 (−i 2 )
′
′
mc
mc
Fσπ
Fππ
′
′



2


S
sin2θ
−4L
cosθ(1−cos2θ)2zy− e (−i h̄ωk ) 



2
mc2

2
 mc

4L
cosθ(1−cos2θ)
S
sin2θ
y
z


Ψ=0◦





=


0
−Sz cosθ(1 − cos2θ)h̄ωk 
e2


− mc

2 (−i mc2 )Scosθ(1−cos2θ)
0
z


Ψ=90◦




Fmag = −(4.27)(4.28)となる。

πσ ′ 散乱過程に注目すると、Ψ=0◦ のときは軌道磁気モーメントからの信号を捉えているが、Ψ=90◦ になるとスピン
磁気モーメントからの信号を捉えるようになることがわかる。
このように、アジマス角が変われば、同じ πσ ′ 散乱過程でも全く違った信号に変わり、それは他の σσ ′ , σπ ′ , ππ ′ 散乱過程でも同じである。

任意の Ψ
今回のモデルの場合、任意の Ψ でのスピン/軌道磁気モーメントを XYZ 軸で表現するのは容易であり、⟨QS ⟩ (X, Y, Z)Ψ = (0, Sz sinΨ, Sz cosΨ)(4.29)⟨QL ⟩ (X, Y, Z)Ψ = (0, −Ly cosΨ, Ly sinΨ)(4.30)

である。

[3,0,0] のような単純な散乱ベクトル以外での XYZ 軸への変換は回転行列を利用する。
下記に一般的な xyz 軸から XYZ 軸への変換手順を示す。

”xyz 軸における散乱ベクトル G”と”Ψ=0◦ の時に Z 軸と一致する結晶軸を H0 (x, y, z)”をG(x, y, z) = (H, K, L)(4.31)

H0 (x, y, z) = (H0 dx, H0 dy, H0 dz)(4.32)

73

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）とする。
ここで、xyz 軸の基底ベクトル êx0 , êy0 , êz0 はêz0 (x, y, z) = (0, 0, 1)(4.33)

êx0 (x, y, z) = (1, 0, 0)(4.34)êy0 (x, y, z) = (0, 1, 0)(4.35)、Ψ=0◦ の時 XYZ 軸 (X, Y, Z)Ψ=0◦ の基底ベクトル êX0 , êY0 , êZ0 は、

êZ0 (x, y, z) = Ĥ0 (x, y, z)(4.36)êX0 (x, y, z) = −Ĝ(x, y, z)(4.37)êY0 (x, y, z) = êZ0 (x, y, z) × êX0 (x, y, z)(4.38)であるので、xyz 軸基底から XYZ 軸基底 (Ψ=0◦ ) への変換は
 
 
X
x
 
 
(êX0 , êY0 , êZ0 )  Y = (êx0 , êy0 , êz0 ) y 

(4.39)

Z Ψ=0◦
z
 


  
 
X
êX0
x
êX0
x
 


  
 
= êY0  (êx0 , êy0 , êz0 ) y  = êY0  y 
Y 
Z

Ψ=0◦

êZ0

z

êZ0

(4.40)zとなる。
さらに、ここから X̄Ψ=0◦ 軸周りの Ψ 軸回転を考えると
 

 
X
1
0
0
X
 

 
 Y  = 0 cosΨ sinΨ  Y 

Z



 
x


 
= 0 cosΨ sinΨ êY0  y 
0 −sinΨ cosΨZ Ψ=0◦
0 −sinΨ cosΨ
êZ0
z
Ψ

 
êX0
x

 
=  cosΨêY0 + sinΨêZ0  y 
1

0

0



−sinΨêY0 + cosΨêZ0
z
 
 
RΨ11 RΨ12 RΨ13
x
x

 
 
= RΨ21 RΨ22 RΨ23  y  = RΨ y 
RΨ31 RΨ32 RΨ33
z
z

êX0

(4.41)(4.42)

(4.43)

K
L
H,
RΨ12 = −,RΨ13 = −
|G|
|G|
|G|
(
)
H0 dx
H0 dz K
H0 dy LsinΨ + −+cosΨ
RΨ21 =|H 0 ||H 0 | |G||H 0 | |G|
(
)
H0 dy
H0 dz HH0 dx L
RΨ22 =sinΨ + −+cosΨ
|H 0 ||H 0 | |G||H 0 | |G|
(
)
H0 dz
H0 dx K
H0 dy H
RΨ23 =sinΨ + −+cosΨ
|H 0 ||H 0 | |G||H 0 | |G|()
H0 dy LH0 dz K
H0 dx
cosΨ +−sinΨ
RΨ31 =|H 0 ||H 0 | |G||H 0 | |G|
(
)
H0 dy
H0 dz HH0 dx L
RΨ32 =cosΨ +−sinΨ
|H 0 ||H 0 | |G||H 0 | |G|
(
)
H0 dz
H0 dx K
H0 dy H
RΨ33 =cosΨ +−sinΨ
|H 0 ||H 0 | |G||H 0 | |G|

RΨ11 = −(4.44)

74

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

1.8

1.8

1.6 ( i ) S z = 2, L y = 8, 2q =30°

1.6 ( ii ) S z = 2, L y = 2, 2q =30°

1.4

1.4

| Fmag || Fmag |

2

1.2

2

1.2

I π π'

1
0.8
0.6

0.6

I π σ'

0.4

0.4

0.2

0.2

0
-180

0
-180

-120

-60

0

I π π'

1
0.8

60

120

180

I π σ'-120

Y (degree)

-60

0

60

120

180

Y (degree)

図 4.4.1 非共鳴磁気散乱の (300) 反射でのアジマス依存性。Ψ=0◦ において Z ∥ c。

(i)Sz =2, Ly =8, 2θ=30◦ , (ii)Sz =2, Ly =2, 2θ=30◦

式 (4.29), (4.30) を用いた各散乱過程における散乱振幅はh̄ωke2(−i 2 )(⟨QL ⟩ · B ⊥ + ⟨QS ⟩ · D)
mc2
mc (
)
′
′
Fσσ
Fπσ
e2
h̄ωk
′
′
= − 2 (−i 2 )
′
′
mc
mc
Fσπ
Fππ
′
′

Fmag = −′Fσσ′ = Sz sin2θcosΨ′2
Fσπ′ = 4Ly cosθ(1 − cos2θ) cosΨ + Sz cosθ(1 − cos2θ)sinΨ′2Fπσ′ = −{4Ly cosθ(1 − cos2θ) cosΨ + Sz cosθ(1 − cos2θ)sinΨ}′Fππ′ = Sz sin2θcosΨ(4.45)
となる。
′
′
′
′
′
′
2
非共鳴磁気散乱においては、Fσσ′ = Fππ ′ , Fσπ ′ = −Fπσ ′ であり、実験では図 (4.4.1) のように回折強度 Iππ ′ (=∥Fππ ′ ∥ )
′
′
′
と Iπσ
′ あるいは Iσσ ′ と Iσπ ′ のアジマス角依存性が得られる。

ここでは、(i)Sz =2, Ly =8, 2θ=30◦ , (ii)Sz =2, Ly =2, 2θ=30◦ でのアジマス依存性を示した。
′′′′
Sz と Ly の値の違いで Iππ′ と Iπσ ′ の強度比が現れ、Iπσ ′ では最大値をとる Ψ の位置が変化する。Iππ ′ は Sz しか観測しないため (i) と (ii) で違いは現れない。これらの差から ⟨QS ⟩ と ⟨QL ⟩ を推定できる。
′′
ただし、実験では Iππ′ と Iπσ ′ の双方に同じ大きさの Scale 因子が掛かるため、⟨QS ⟩ と ⟨QL ⟩ の絶対値を見積もることは難しい。

4.4.2

非共鳴磁気散乱の偏光依存性

前節でアジマス角依存性により磁気モーメントを推定できること説明した。本節では、偏光依存性もアジマス角依存性と同じ情報をもつことを示す。
式 (4.6) で示したように、散乱振幅は磁気モーメントと偏光因子の内積である。後述するように、共鳴散乱での散乱
振幅も秩序変数である多重極子モーメントと偏光因子の内積外積である。アジマス角依存性は、偏光因子を固定し、試料に固定された多重極子モーメントをアジマス角周りで回して実験パラメータとしていた。逆に多重極子モーメントを固定して偏光因子を実験パラメータとする偏光依存性は、アジマス角依存性と比べて多重極子モーメントと偏光因子のどちらを実験パラメータとして動かしているかの違いしかなく、アジマス角依存性ともつ情報は同じである。
アジマス依存性と同様に偏光依存性を第 (4.4) 節のモデルで計算すると75

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

1.8

1.8

1.6 ( i ) S z = 2, L y = 8, 2q =30°
1.4 ξ = 0º

I π π'

30º 60º 90º

Y

=0°

Y

ξ = 0º30º 60º 90º
1.2

I π π'

2| Fmag |

2

1
0.8

1
0.8

0.6

0.6

0.4

0.4

0.2

0.2

0
-180

-120

-60

0

60

h (degree)

=0°

1.4

1.2

| Fmag |1.6 ( ii ) S z = 2, L y = 2, 2q =30°120

I π σ'

0
-180

180

-120

-60

0

60

h (degree)

I 120
π σ'

180

図 4.4.2 非共鳴磁気散乱の (300) 反射, Ψ=0◦ (Z ∥ c) での入射偏光依存性。

(i)Sz =2, Ly =8, 2θ=30◦ , (ii)Sz =2, Ly =2, 2θ=30◦

1.8

1.8

1.6 ( i ) S z = 2, L y = 8, 2q =30°

| Fmag |

2

1.2

30º
90º
60º

1
0.8

1

0.6

0.4

0.4

0.2

0.2
-120

-60

0

60

120

ξ = 0º
30º

0.8

0.6

0
-180

Y

=30°

1.4

ξ = 0º2

1.2

1.6 ( ii ) S z = 2, L y = 2, 2q =30°

=30°
| Fmag |

1.4

Y

0
-180

180

60º 90º-120h (degree)

-60

0

60

120

180

h (degree)

図 4.4.3 非共鳴磁気散乱の (300) 反射, Ψ=30◦ での入射偏光依存性。Ψ=0◦ において Z ∥ c。

(i)Sz =2, Ly =8, 2θ=30◦ , (ii)Sz =2, Ly =2, 2θ=30◦e2h̄ωk(−i 2 )(⟨QL ⟩ · B ⊥ + ⟨QS ⟩ · D)
mc2
mc
e2
h̄ωk
= − 2 (−i 2 )mc
mc(
4Ly (1 − cos2θ)2 [cosθsin(ξ − η)cosΨ + sin2θ{cos(ξ + η) − cos(ξ − η)}sinΨ]Fmag = −+Sz [cosθ(1 − cos2θ)sin(ξ − η)sinΨ + sin2θcos(ξ − η)cosΨ] )(4.46)(4.47)(4.48)

Ψ=0◦ の場合の結果を図 (4.4.2) のように示す。
ここでは、図 (4.4.1) と比較できるように、(i)Sz =2, Ly =8, 2θ=30◦ , (ii)Sz =2, Ly =2, 2θ=30◦ での入射偏光依存性を示した。回折 X 線の回折偏光成分として、ξ=0◦ , 30◦ , 60◦ , 90◦ を計算しているが、全て最大値が同じであり、アジマ
ス依存性ほど (i) と (ii) の差が現れない。これはモデルにおける磁気モーメントの方向が単純すぎるためである。しかし、(i) では ξ=0◦ 成分が η=30◦ 偏光入射で最大値となる。つまり η=30◦ 偏光入射は回折により ξ=0◦ 偏光になる。それに対して、(ii) では η=10◦ 偏光入射が ξ=0◦ 偏光になる。

Ψ=0◦ では単純すぎたので、Ψ=30◦ の場合の結果を図 (4.4.2) に示す。
Ψ=0◦ に比べて、回折偏光成分に明らかな差が現れる。このように偏光依存性には特徴が現れにくいアジマス角度と現れ易いアジマス角度が存在する。実験では可能な限り異なるアジマス角での偏光依存性測定をしておく必要がある。
しかし、アジマス角を回すことにより磁場方向が変わってしまうことを考えると、異なる散乱ベクトル (反射指数 HKL)で測定することが有効だろう。

76

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

偏光依存性のメリット・デメリット
偏光依存性を測定するメリットをいくつか挙げる。
メリット– アジマス依存性と異なり、偏光依存性は印加磁場方向が変化しない。

– 試料を回転させる必要がないため、X 線の試料への照射位置が動かない。
– アジマス依存性測定中に生じる散乱ベクトルの微調整の必要がない。
デメリット– 入射偏光制御装置の立ち上げが必要– エネルギーを大きく (1keV 以上) 動かす場合、入射偏光制御装置に用いるダイヤモンドを変更し制御装置の再立ち上げが必要
偏光依存性を測定する最大のメリットは磁場中低温での磁気モーメントの推定が可能になることである。コンプレッ
サー型の冷凍機及び低磁場の小型の電磁石であれば試料と共に回転させることが可能であるが、3K 以下に到達可能なヘリウムフロー型の冷凍機や超伝導マグネットは試料を回転させることは出来ない。そのため、アジマス依存性測定ではマグネット中の試料のみを回転させていた。そのため、磁場中でのアジマス依存性測定を測定すると磁場方向は常に変化し、磁場中での多重極子モーメントの推定は非常に困難であった。アジマスの代わりに偏光を回すことにより、試料と磁場の関係が固定され安定した状態で磁場中の多重極子モーメントを観測可能である。
試料を動かさないことに付随するメリットは他にも存在する。アジマス角依存性の測定では、Ψ 回転ステージの上に試料をセットするが、1mm2 以下の試料を狂いなく回転中心におくことや Ψ 回転軸と散乱ベクトルを 0.1◦ のズレもな
く設置することは不可能である。そのため、実験ではアジマス角を動かす度に X 線の照射位置がわずかに変化し、物理とは非本質的な強度の変化が存在した。また、散乱ベクトルのズレを調整するため、ωscan や χscan 等を常に行う必要があった。ダイヤモンド移相子による偏光制御は、X 線がダイヤを透過する際の π 偏光成分と σ 偏光成分の位相のズレを利用する。X 線の光路は変わることなく、偏光のみ変化するため試料への照射位置も変わらず非本質的な強度の変動もなく、散乱ベクトルの調整に伴う時間も省略できる。
一方、偏光依存性のデメリットは入射偏光制御装置の立ち上げの手間である。入射偏光制御装置の立ち上げ時間は 4-6
時間掛かる。また、アナライザーと同様に一つのセッティングで利用できるエネルギー範囲は狭く、同一元素の LII 吸
収端での実験後は LIII 吸収端で再度調整する必要がある。しかし、偏光制御のメリットは制御装置の立ち上げの手間を大きく上回る。

77《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

5

共鳴散乱振幅における多重極子演算子と波数ベクトル, 偏光因子5.1

共鳴散乱：多重極子モーメントからの共鳴散乱

第 3.2 節から導いたように、式 (3.22) は、A の二次摂動項である。非共鳴散乱では、式 (3.20), (3.21) より各々電気
単極子モーメント (電荷) と磁気双極子モーメント (磁気モーメント) からの非共鳴散乱 (式 (4.6), (4.17)) を導いたが、式 (3.22) からは、電気単極子, 磁気双極子, 電気四極子, 磁気八極子モーメントからの共鳴散乱が導かれる。その際、第2.3 節で述べたように共鳴過程には E1 遷移と E2 遷移があるため、各々の遷移で異なる散乱項が導かれる。以下に簡略
的に各散乱項の導出を示す。
式 (3.22) において, 式 (3.15) の運動量密度演算子を∑(pj − ih̄k × sj )(1 + ik · r j )

(5.1)

∑
∑ ′
i
i ′
′
mωr3 ⟨a | j ′ ε · r j ′ (1 − 2 k · r j ′ )|b⟩⟨b| j ε · r j (1 + 2 k · r j )|a⟩ωkh̄ωk − h̄ωr + iΓb /2
b
∑
∑
)
′
′
′
h̄2 ωr ⟨a | j ′ ε · ((lj ′ + 2sj ′ ) × k )|b⟩⟨b| j ε · ((lj + 2sj ) × k)|a⟩(5.2)

J (k) =jと近似すると、

Freso

e2 ∑
= − 2
mc
+

(

h̄ωk − h̄ωr + iΓb /2

4m ωk

(5.3)

さらに、下記のような原子の電気双極子と磁気双極子及び電気四極子の演算子を導入すると

Rα

=

∑

rjα

(5.4)

j

Mα

=

Qαβ

=

1 ∑
(ljα + 2sjα )2m j∑
rjα rjβ

(5.5)
(5.6)

j

Freso

=


∑ ⟨a′ |Rα − i Qαγ kγ′ |b⟩⟨b|Rβ + i Qβδ kδ′ |a⟩e2 ∑ mωr3 ∑ ′
2
2
− 2
εα εβ
mc
ωk
h̄ωk − h̄ωr + iΓb /2
α,β
γ,δ
b

∑ h̄2 ωr ∑∑ ⟨a′ |Mα × kγ′ |b⟩⟨b|Mβ × kδ |a⟩+ε′α εβ2 ωkh̄ωk − h̄ωr + iΓb /2
b

α,β

(5.7)

(5.8)

γ,δ

となる。
式 (5.7) の第 1 項は電気的な遷移であり, ⟨a′ |Rα |b⟩⟨b|Rβ |a⟩ の電気双極子 (E1) 遷移と ⟨a′ |Qαγ |b⟩⟨b|Qβδ |a⟩ の電気四極子 (E2) 遷移による散乱項である。これらの電気双極子や電気四極子は遷移を媒介するモーメントであり、秩序変数の電気双極子や電気四極子ではないので、混同しないこと。
式 (5.7) の第 2 項は磁気的な遷移 (M1 遷移) であるが、電気的な E1 遷移に比べて 5 桁程度小さく無視できる。

次節からは、E1, E2 遷移による散乱項を秩序変数 (電気単極子, 磁気双極子, 電気四極子, 磁気八極子モーメントの大きさ) と偏光/波数ベクトルで表わす。

78

《講義ノート》

5.2

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

E1,E2 遷移

秩序変数と偏光/波数ベクトルを用いた E1, E2 遷移の散乱の記述はいくつかあるが、最近では長尾さんと五十嵐さん
により記述された以下の定式を用いている。

(E1)
Freso

=

2
∑

(ν)

αE1 (ωk )

ν=0
(E2)
Freso

=

4
∑

2ν+1
∑

PE1,µ (ε, ε′ )⟨zµ(ν) ⟩,(5.9)

PE2,µ (ε, ε′ , k, k′ )⟨zµ(ν) ⟩(5.10)

(ν)

µ=1
(ν)

αE2 (ωk )

ν=0

2ν+1
∑

(ν)

µ=1
(ν)

(ν)

この定式は、E1 遷移と E2 遷移の干渉などのエネルギー依存性に対応しており, αE1 (ωk ) と αE2 (ωk ) がエネルギー依存をもつ係数である。ν は多重極子モーメントの rank であり、表 (3) の通りである。E1 遷移では四極子モーメント以下、E2 遷移では十六極子以下であるため、それぞれ ν ≤ 2, ν ≤ 4 である。アジマス依存性や偏光依存性はエネルギー
(ν)

(ν)

を固定して測定する。そのため、これらの依存性を測定する際には、αE1 (ωk ) と αE2 (ωk ) は一定値として扱って良く、

scale 因子と変わりない。
(ν)
続いて、⟨zmu ⟩ は多重極子演算子の期待値 (各多重極子モーメントの大きさ) である。µ は各多重極子モーメントにお(ν)

ける独立成分であり、⟨zmu ⟩ とともに表 (3) に定義している。
最後に、PE1,µ (ε, ε′ ) と PE2,µ (ε, ε′ , k, k′ ) は、⟨zmu ⟩ に対応した偏光/波数ベクトル因子である。具体的な表式は表 (4)
(ν)

(ν)

(ν)

に示している。例えば、PE1,µ (ε, ε′ ) における µ=1,2,3 は Jx , Jy , Jz に対応し、以下のようになる。
(1)−i′ (ε′ × ε)µ=1
 ′εy εz − ε′z εy= ε′z εx − ε′x εz 
ε′x εy − ε′y εx 1

PE1,1 (ε, ε′ ) =(1)= ε′y εz − ε′z εy(1)
PE1,2 (ε, ε′ )

=

′

′

−i (ε × ε)µ=2= ε′z εx − ε′x εz(1)
PE1,3 (ε, ε′ )′(5.11)(5.12)(5.13)(5.14)(5.15)=′−i (ε × ε)µ=3(5.16)=

ε′x εy − ε′y εx(5.17)

79

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

表 3 多重極子演算子の rank と点群 Oh の規約表現で分類した多極子演算子の慣用表記。また、Jx Jy Jz = Jx Jy Jz + Jy Jz Jx + Jz Jx Jy 。

Nortionrank(ν)01

規約表現

Γ1g (A1g )Γ4u (T1u )

1
1
2
3

2

Γ3g (Eg )Γ5g (T2g )

1
2
3
4
5

3

Γ2u (A2u )Γ4u (T1u )Γ5u (T2u )

1
2
3
4
5
6
7

4

Γ1g (A1g )Γ3g (Eg )Γ4g (T1g )Γ5g (T2g )等価演算子

立方調和関数

1
Jx

1
x

Jy
Jz

y
z

O20 = 21 {3Jz2 − J(J + 1)}√
O22 = 23 (Jx2 − Jy2 )/√Oyz = 23 (Jy Jz + Jz Jy )12
2
2 (3z − r )√32− y2 )
2 (x
√

ν

1

(0)
z1
(1)
z1
(1)
z2
(1)
z3
(2)
z1
(2)
z2
(2)
z3
(2)
z4
(2)
z5
(3)
z1
(3)
z2
(3)
z3
(3)
z4
(3)
z5
(3)
z6
(3)
z7
(4)
z1
(4)

2

z2

3

z3

4
5
6
7
8
9

(4)

(4)
z4
(4)
z5
(4)
z6
(4)
z7
(4)
z8
(4)
z9

3yz
√
3zx
√
3xy
√
15xyz
1
2
2
2 x(5x − 3r )1
2
2
2 y(5y − 3r )√Ozx = 23 (Jz Jx + Jx Jz )√Oxy = 23 (Jx Jy + Jy Jx )√Txyz = 615 Jx Jy Jz
Txα = 21 (2Jx3 − Jx Jy2 − Jz2 Jx )Tyα = 21 (2Jy3 − Jy Jz2 − Jx2 Jy )Tzα = 21 (2Jz3 − Jz Jx2 − Jy2 Jz )√Txβ = 615 (Jx Jy2 − Jz2 Jx )√Tyβ = 615 (Jy Jz2 − Jx2 Jy )√Tzβ = 615 (Jz Jx2 − Jy2 Jz )H40 = 45

√

7 4
4
4
3 [Jx + Jx + Jx

− 35 J(J + 1){J(J + 1) − 13 }]√H24 = − 45 [ 67 (Jx2 − Jy2 )Jz2
− {J(J + 1) − 56 }(Jx2 − Jy2 )]√H44 = 485 [35Jz4 − 30J(J + 1)Jz2+ 3J(J + 1)J(J + 1) − 2+ 25Jz2 − 7Jx4 + Jy4 − Jx2 Jy2 ]√αHx = 835 (Jy3 Jz − Jy Jz3 )√Hyα = 835 (Jz3 Jx − Jz Jx3 )√Hzα = 835 (Jx3 Jy − Jx Jy3 )√Hxβ = 85 [2Jx2 Jy Jz − (Jy3 Jz + Jy Jz3 )]√Hyβ = 85 [2Jy2 Jz Jx − (Jz3 Jx + Jz Jx3 )]
Hzβ =

√

5
2
3
3
8 [2Jz Jx Jy − (Jx Jy + Jx Jy )]

80

1
2
2
2 z(5z − 3r )√15x(y 2 − z 2 )√2
15y(z 2 − x2 )
√2
15
2
2
2 z(x − y )
√
5 21
3 4
4
4
4
12 (x + y + z − 5 r )
√
15
4
4
4
2
2 2
12 {7(2z − x − y ) − 6(3z − r )r }√152
− y 2 )(r2 − 7z 2 )
4 (x
√
35
yz(y 2 − z 2 )√235zx(z 2 − x2 )√2
35
xy(x2 − y 2 )√25yz(7x2 − r2 )√2
5zx(7y 2 − r2 )√2
52
22 xy(7z − r )《講義ノート》

5.3

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

散乱の結晶構造因子
(E1)

(E2)

ここまで扱った Freso , Freso は、いわば原子散乱因子である。
従って、例えばサイト n が 3 つの多重極子秩序状態の結晶構造因子 F は、結晶周期 e−iκ·r n を掛けたものとなる。
そこで、O22 型の四極子秩序が存在し z 軸方向に磁気秩序が誘起されたとすると、その E1 遷移を考えた F はF= αE1,n=1 (ωk )PE1,3 (ε, ε′ )⟨z3 ⟩e−iκ·r 1 + αE1,n=1 (ωk )PE1,2 (ε, ε′ )⟨z2 ⟩e−iκ·r 1
(1)
(1)
(1)
(2)
(2)
(2)
+α
(ω )P
(ε, ε′ )⟨z ⟩e−iκ·r 2 + α(ω )P(ε, ε′ )⟨z ⟩e−iκ·r 2
(1)

(1)

(1)

(2)

(2)

(2)

(5.18)(5.19)kk
3
2
E1,n=2E1,3
E1,n=2E1,2(1)(1)(1) −iκ·r 3
(2)
(2)
(2)
′
′
+αE1,n=3 (ωk )PE1,3 (ε, ε )⟨z3 ⟩e+ αE1,n=3 (ωk )PE1,2 (ε, ε )⟨z2 ⟩e−iκ·r 2(5.20)

ここで、四極子秩序が単純な反強的であり、n=1,3 を A 副格子, n=2 が A 副格子とは反転した B 副格子とする。同様に磁気秩序も単純な反強的であるとする。また、結晶周期の位相を ϕ として (n, ϕ) = (1, 0), (1, π), (2, 2π) とすると、

F= αE1,n=A (ωk )PE1,3 (ε, ε′ )⟨z3 ⟩(+1) + αE1,n=A (ωk )PE1,2 (ε, ε′ )⟨z2 ⟩(+1)
(1)

=

(1)

(1)

(2)

(2)

(2)

(5.21)

(1)
(1)
(1)
(2)
(2)
(2)
+αE1,n=B (ωk )PE1,3 (ε, ε′ )(−⟨z3 ⟩)(−1) + αE1,n=B (ωk )PE1,2 (ε, ε′ )(−⟨z2 ⟩)(−1)
(1)
(1)
(1)
(2)
(2)
(2)
+αE1,n=A (ωk )PE1,3 (ε, ε′ )⟨z3 ⟩(+1) + αE1,n=A (ωk )PE1,2 (ε, ε′ )⟨z2 ⟩(+1)
(1)
(1)
(1)
(2)
(2)
(2)
2{αE1,n=A (ωk )PE1,3 (ε, ε′ )⟨z3 ⟩ + αE1,n=A (ωk )PE1,2 (ε, ε′ )⟨z2 ⟩}
(2)
(2)
(2)
(1)
(1)
(1)
+{αE1,n=B (ωk )PE1,3 (ε, ε′ )⟨z3 ⟩ + αE1,n=B (ωk )PE1,2 (ε, ε′ )⟨z2 ⟩}(2)(5.22)(5.23)(5.24)(5.25)(2)

副格子が異なってもほとんどエネルギー依存性が変わらないのであれば、αE1,n=A (ωk ) = αE1,n=B (ωk ) として、

F

= 3αE1 (ωk )PE1,3 (ε, ε′ )⟨z3 ⟩ + 3αE1 (ωk )PE1,2 (ε, ε′ )⟨z2 ⟩
(1)

(2)

(1)

(1)

(2)

(2)

(2)

(2)

(1)

(5.26)
(1)

(2)

(2)

実験的に αE1 (ωk ) を決定することは難しい。αE1 (ωk ) を決定できない場合は、3αE1 (ωk )⟨z3 ⟩ と 3αE1 (ωk )⟨z2 ⟩ を
係数 A13 と A22 としてF= A13 PE1,3 (ε, ε′ ) + A22 PE1,2 (ε, ε′ )
(2)

(1)

(5.27)(2)となるため、共鳴 X 線回折実験で αE1 (ωk ) を決定できない場合は、秩序変数の絶対値を求めることは出来ない。しかし、共鳴 X 線回折という一つのプローブで複数の秩序変数の温度変化や磁場変化などの相対変化を詳細に追うことが可能である。そのため、複数の秩序変数が絡み合う多重極子秩序変数の研究には非常に有用である。
第 4.4 節と同様に、実験室系 XYZ 座標と結晶系 xyz 座標に注意し、アジマス角依存性や偏光依存性を計算出来る。
実験で得たアジマス角依存性や偏光依存性の結果を説明できるモデルを探し出すことで、秩序変数の絞り込みを行うことが可能となる。共鳴 X 線回折実験は、中性子回折実験で磁気構造を調べることと非常に似ている。共鳴 X 線回折実験で得られる結果は、多重極子秩序の構造モデルを絞り込むことが可能であって、決定できるわけではない。また、秩序変数として複数の多重極子モーメントの秩序構造情報を得られることはメリットではあるが、反面、実験結果を再現するためのパラメータが多くなり間違ったパラメータで実験結果を再現できてしまう可能性がある。したがって当然ではあるが、共鳴 X 線回折実験において解析を行う上では、磁気双極子に関しては中性子回折実験を参考にすべきであるし、その他の物性を熟慮してモデルを立てなければならない。

81

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

表 4 多重極子モーメントの独立成分に対応した幾何学 (偏光/波数ベクトル) 因子。
(0)

PE1,µ

=

(ε′ · ε);(µ = 1)

′

(1)
PE1,µ
(2)
PE1,µ

=
=

−i(ε × ε)µKµ (ε′ , ε);(µ = 1, 2, 3);(µ = 1, 2, 3, 4, 5)

PE2,µ

(0)

=

;(µ = 1)

(1)
PE2,µ

=

3
√
{(k′ · k)(ε′ · ε) + (k′ · ε)(ε′ · k)}
4 5−i 4√310 {(ε′ · ε)(k′ × k) + (k′ · k)(ε′ × ε)+(k′ · ε)(ε′ × k) + (ε′ · k)(k′ × ε)}µ− 2√314 {(ε′ · ε)Kµ (k, k′ ) + (k′ · k)Kµ (ε, ε′ )+Kµ (k′ × k, ε × ε′ )}
∑3
1
′
′
′
′
i 4√
j=1 {(k × k)j Kj+2 (ε , ε) + (ε × ε)j Kj+2 (k , k)2+(k′ × ε)j Kj+2 (ε′ , k) + (ε′ × k)j Kj+2 (k′ , ε)}

(2)
PE2,µ

=

PE2,µ=1

(3)

=

(3)

=

PE2,µ=j+1(3)
PE2,µ=j+4

=

√

(4)
PE2,µ=1(4)
PE2,µ=2=PE2,µ=3

(4)

=

(4)

=

PE2,µ=j+3(4)
PE2,µ=j+6==;(µ = 1, 2, 3, 4, 5)

√

5
′
′
′
′
2 {(k × k)j εj · εj + (ε × ε)j kj · kj
(1)+(k′ × ε)j ε′j · kj + (ε′ × k)j kj′ · εj } + 12 PE2,j
√
∑
3
3
′
′
′
′ ′′
′′
′
i 16
j ′ j ′′ εjj j (εj ′ εj − εj ′′ εj )2 {(k × k)j∑′′′+(ε × ε)j j ′ j ′′ εjj ′ j ′′ (kj ′ kj ′ − kj ′′ kj ′′ )∑+(k′ × ε)j j ′ j ′′ εjj ′ j ′′ (ε′j ′ kj ′ − ε′j ′′ kj ′′ )∑+(ε′ × k)j j ′ j ′′ εjj ′ j ′′ (kj′ ′ εj ′ − kj′ ′′ εj ′′ )}i 38;(µ = 1, 2, 3)

√ (0)
2 15
{ (kx′ kx ε′x εx + ky′ ky ε′y εy + kz′ kz ε′z εz ) − 5PE2,1 }√ 15 442(kx′ kx ε′x εx + ky′ ky ε′y εy − 2kz′ kz ε′z εz )4 √− 12 67 {(k′ × k)K2 (ε′ , ε) + (ε′ × ε)K2 (k′ , k)+(k′ × ε)K2 (ε′ , k) + (ε′ × k)K2 (k′ , ε)}
√
3 14
′
′
′
′
4 √(kx kx εx εx + ky ky εy εy )2
− 32 21{(k′ × k)K1 (ε′ , ε) + (ε′ × ε)K1 (k′ , k)+(k′ × ε)K1 (ε′ , k) + (ε′ × k)K1 (k′ , ε)}∑2√{Kj+2 (k′ , k) j ′ j ′′ εjj ′ j ′′ (ε′j ′ εj ′ − ε′j ′′ εj ′′ )8 6∑
+Kj+2 (ε′ , ε)j j ′ j ′′ εjj ′ j ′′ (kj′ ′ kj ′ − kj′ ′′ kj ′′ )∑+Kj+2 (k′ , ε)j j ′ j ′′ εjj ′ j ′′ (ε′j ′ kj ′ − ε′j ′′ kj ′′ )∑+Kj+2 (ε′ , k)j j ′ j ′′ εjj ′ j ′′ (kj′ ′ εj ′ − kj′ ′′ εj ′′ )}√3 [{7k ′ kj − 3(k′ , k)}Kj+2 (ε′ , ε)j
4 42;(j = 1, 2, 3);(j = 1, 2, 3);(j = 1, 2, 3)+{7ε′j εj − 3(ε′ , ε)}Kj+2 (k′ , k)=

K1 (a, b)
K2 (a, b)==

K3 (a, b)
K4 (a, b)
K5 (a, b)===

εjj ′ j ′′+{7ε′j kj − 3(ε′ , k)}Kj+2 (k′ , ε)+{7kj′ εj − 3(k′ , ε)}Kj+2 (ε′ , k)];(j = 1, 2, 3)

1(3az bz − a · b)
2√
3(ax bx − ay by )√2
3
(ay bz + az by )√2
3
(az bx + ax bz )√2
3
2 (ax by + ay bx )=+1

; (jj’j” = xyz, yzx, zxy)

=
=

-1
0

; (jj’j” = xzy, zyx, yxz)
; (jj’j” = otherwise)82

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

まとめ
共鳴 X 線散乱は直感的に取り扱うことが難しく、共鳴 X 線散乱だけを理解しようとしても非共鳴と共鳴の区別をもつことも怪しくなってくる。また、実験室系での X 線回折実験では散乱に伴う「X 線の偏光」を意識することが少な
く、偏光状態の観測が重要になる共鳴 X 線回折実験を理解し難い。
そこで、第一章と第二章では、X 線を照射された電子の運動を古典的に取り扱った場合の X 線の散乱や吸収を解説
し、非共鳴の散乱や散乱に伴う偏光状態の直感的なイメージを示した。これは、比較的イメージし易い古典的な描像と比較しながら、量子力学的な記述で求められる共鳴 X 線散乱を理解することにより、共鳴 X 線散乱の描像を直感的に得やすくなると考えたためである。
第三章では、古典的な解釈では定量的に表現できない共鳴散乱を二次摂動を取り入れることにより、表式化し Thomson
散乱, 非共鳴磁気散乱, 共鳴散乱の３つの項に分離できることまでを示した。これで、共鳴散乱と非共鳴散乱の区別がイメージできれば良い。実験では、「散乱 X 線の偏光ベクトルが入射 X 線に対してどのように変化したかという偏光ベクトルの変化」さらに、「その偏光ベクトルと散乱ベクトルの関係」という幾何学的な関係を得る。したがって、第三章で表式化した各々の散乱項において、偏光ベクトルと散乱ベクトルが秩序変数に対してどのような関係をもつかが実験結果を解析するために必要となってくる。
第四章では、実験で得られるあるいは与えるパラメータである偏光ベクトルと散乱ベクトルが、秩序変数とどのような関係であるかを具体的な表式で紹介した。実際に実験で得られるアジマス依存性や偏光依存性を一度計算するとよいだろう。また、今回は触れなかったが、これまでは E1 遷移による共鳴散乱と E2 遷移による共鳴散乱の干渉を実験デー
タの解析に取り入れることが困難であり、それぞれ単独に扱われてきたが、近年の理論家による定式化により、我々実験家にも解析できるようになるなど、この数年で解析方法が進歩し続けている。実験法に関しても、移相子を用いた入射偏光の制御と 1K 冷凍機, 超伝導マグネットの組み合わせにより、極低温 · 強磁場での研究が可能になった。X 線検出器も SDD(Silicon Drift Detector) の普及により蛍光 X 線の分離やエネルギー分解能の上昇, 検出下限の向上などの発展を続けている。
第四章以降に関しては、広島大学の松村准教授がまとめられている [9]。松村氏のテキストは非常に丁寧に書かれており、私は共鳴 X 線回折の基本の多くをここから学んだ。最近、テキストを更新され、私が盛り込みたかった内容も触れられている。勝手ながら、読者にはそちらを参照して頂きたい。松村氏のテキストを読むにあたり、もっていた方が良い基礎知識やイメージをこのテキストの第一章や第二章に書いたつもりである。イメージし易い古典的な X 線回折の知識と量子力学的な記述の比較は、量子力学的な記述である共鳴散乱の理解の一助になるはずである。

83

《講義ノート》

物性研究・電子版 Vol. 6, No. 1, 061203 （2016年11月・2017年2月合併号）

謝辞
最後に、テキストを編集して頂いた野原実教授 (岡山大学) や籐秀樹教授 (神戸大学) 及び播磨尚朝教授 (領域代表, 神戸大学) には、この様な執筆の機会を頂き、大変感謝しております。また、本稿の内容が書けるようになるまで多くの方々にお世話になっております。特に稲見俊哉主任研究員 (量子科学技術研究開発機構)、松村武准教授 (広島大学) には深く感謝致します。

参考文献[1] Values are from J. A. Bearden, ”X-Ray Wavelengths”, Review of Modern Physics, (January 1967) pp. 86-99[2] http://www.med.harvard.edu/jpnm/physics/refs/xrayemis.html[3] 大橋 治彦, 平野馨一 編, 放射光ビームライン光学技術入門, 日本放射光学会, (April 2009) pp.236[4] 雨宮慶幸 他, X 線物理学の基礎, 講談社, (May 2012) pp.8
[5] http://henke.lbl.gov/optical constants/filter2.html[6] International Tables for Crystallography[7] 桜井敏雄, X 線結晶解析の手引き, 裳華房, (May 1994)[8] 橋爪弘雄, 岩住俊明, 放射光 X 線磁気分光と散乱, IPC, (April 2007)[9] 松村武, 非共鳴及び共鳴 X 線回折を利用した磁気および多極子秩序状態の観測
http://home.hiroshima-u.ac.jp/tmatsu/Matsumura/Research.html

84