ハミルトニアンの選び方で固体を考える 第1部 第3章 §1.3.1 ★★★ 第1 Born 近似で散乱振幅がポテンシャルの Fourier 変換になる。章の残りはすべてこの1式の読み替え
前の章(第1部第2章)では、結晶を格子と基本構造(basis)の組で表し、格子の周期を持つ関数は逆格子ベクトル G だけを波数に持つ Fourier 級数で書けることを見た。ただしそこまでは幾何の約束であり、実在の結晶の格子定数や原子の並びをどうやって知るのかには答えていない。この章はその答えである。波長が原子間距離と同じ程度の波(X 線・電子線・中性子線)を結晶に当て、散乱された波の強さを角度ごとに測ると、強い散乱は散乱ベクトルが逆格子ベクトルに一致する方向にだけ現れ、その強さは単位胞の中の原子の並びで決まる。
この章で選ぶものは、入射粒子1つと、静止した結晶が作るポテンシャルだけである。電子の多体の問題も、原子の熱振動も、いったん全部捨てる。当てる道具は散乱理論の第1 Born 近似で、この近似のもとでは散乱振幅がポテンシャルの Fourier 変換そのものになる。そこから3つのものが出てくる。散乱が起きる方向を決める Laue 条件(Bragg の式・Ewald 球)、強度を決める構造因子と消滅則、そして探り針ごとに異なる原子形状因子である。最後に、捨てたもの(多重散乱・熱振動・結合による電子密度の歪み)を数え、それぞれを第V巻(近似の数理)の言葉で札にする。
選ぶ H
質量 m の粒子(中性子や電子)が、結晶の作るポテンシャル V(r) の中を飛ぶとする。この章で選ぶハミルトニアンは
H^=2mp^2+V(r^),V(r)=R∑j∑vj(r−R−dj)(1.3.1)
である。R は格子ベクトル、dj は単位胞の中の j 番目の原子の位置、vj はその原子1個が作るポテンシャルである。この形を書いた時点で、すでに大きく削っている。結晶をつくる原子は止まっていて(熱振動なし)、入射粒子は結晶の状態を変えず(弾性散乱だけ)、結晶全体のポテンシャルは原子1個ずつのポテンシャルの和で書ける、と仮定した。何を削ったかは 1.3.5 節で数え直す。
X 線は粒子でなく電磁波なので、式 (1.3.1) の形ではない。X 線は電子にだけ散乱され、その振幅は電子の個数密度 n(r) の Fourier 変換になる(1.3.4 節)。結論の形は中性子・電子と同じなので、この節では式 (1.3.1) で話を進める。
当てる道具:第1 Born 近似
エネルギー E=ℏ2k2/(2m) の平面波 eik⋅r を入射させる。定常状態の Schrödinger 方程式 (∇2+k2)ψ=(2m/ℏ2)Vψ は、外向きの球面波を表す Green 関数 −eik∣r−r′∣/(4π∣r−r′∣) を使って積分方程式
ψ(r)=eik⋅r−2πℏ2m∫∣r−r′∣eik∣r−r′∣V(r′)ψ(r′)d3r′(1.3.2)
に書き直せる(Lippmann–Schwinger 方程式)。検出器は結晶から遠くにあるので、r=∣r∣≫∣r′∣ としてk∣r−r′∣≈kr−k′⋅r′(k′=kr/r は検出器の方向を向いた、大きさ k の波数ベクトル)と展開すると、ψ≈eik⋅r+f(k′,k)eikr/r の形になる。ここで右辺の積分の中の ψ を入射波 eik⋅r′ で置き換えるのが第1 Born 近似である。散乱された波がもう一度散乱される寄与を捨てる近似なので、「1回散乱の近似」とも言う。結果は
f(k′,k)=−2πℏ2m∫V(r)e−iq⋅rd3r=−2πℏ2mV~(q),q:=k′−k(1.3.3)
となる。q を散乱ベクトルと呼ぶ。微分散乱断面積は dσ/dΩ=∣f∣2 で、検出器が数える粒子の数はこれに比例する。Fourier 変換の符号は記法の正本の流儀(V~(q)=∫Ve−iq⋅rd3r、2π は逆変換の側)に合っている。
式 (1.3.3) がこの章のすべての出発点である。散乱振幅は、ポテンシャルを散乱ベクトル q という1つの波数で Fourier 成分に分けたものに比例する。波長を変えたり結晶を回したりして q を動かしながら ∣f∣2 を測ることは、ポテンシャルの Fourier 成分の絶対値を1つずつ読み取ることに当たる。光学で Fraunhofer 回折の振幅が開口の Fourier 変換になるのと同じ構造で(光学と量子エレクトロニクス 第6章)、違うのは「開口」が3次元に周期的に並んだ原子だという点だけである。
格子の和と基本構造の和に分ける
式 (1.3.1) の V を式 (1.3.3) に入れ、積分変数を原子ごとに r=R+dj+s とずらすと、
V~(q)=L(q)R∑e−iq⋅R単位胞の和j∑v~j(q)e−iq⋅dj(1.3.4)
と、格子だけで決まる和 L(q) と、単位胞の中身だけで決まる和の積に分かれる。この分解は「結晶=格子+基本構造」という前の章の分け方を、そのまま散乱振幅に写したものである。L(q) は散乱の方向を決め(1.3.2 節)、単位胞の和は各方向の強さを決める(1.3.3 節)。
例題 1.3.1 原子 N 個の列の干渉関数 ★★★
オリジナル間隔 a で x 軸に並んだ同じ原子 N 個(x=0,a,…,(N−1)a)に、散乱ベクトル q=qex で散乱が起きるとする。
- (1) 格子の和 L(q)=∑n=0N−1e−iqna の絶対値の2乗が ∣L∣2=sin2(Nqa/2)/sin2(qa/2) になることを示せ。
- (2) ∣L∣2 が最大になる q とその高さ、最大の点から最初の零点までの q の幅を求めよ。
- (3) N=8 のとき、主極大と主極大の間の副極大の最大のものの高さは主極大の何割か。数値で答えよ。
解答
(1) 等比級数の和で L=(1−e−iNqa)/(1−e−iqa)。分子と分母からそれぞれ e−iNqa/2・e−iqa/2 をくくり出すとL=e−i(N−1)qa/2sin(Nqa/2)/sin(qa/2) となり、絶対値の2乗をとれば示す式になる。
(2) 分母が 0 になる qa=2πh(h は整数)で、分子も 0 になり、極限は ∣L∣2=N2 である。これが最大値で、それ以外の q では∣L∣2≤1/sin2(qa/2) が N によらない有限の値にとどまる。最大の点 q=2πh/a から q を増やすと、分子 sin(Nqa/2) が初めて 0 になるのはΔq=2π/(Na) のときである。
(3) 数値で ∣L∣2/N2 を qa/2π∈(0,1) で調べると、主極大のすぐ隣の副極大(qa/2π≈0.18 付近)が最大で、その高さは 0.0525、つまり主極大の約 5% である。N が大きくなると、この比は関数 sin2x/x2 の最大の副極大の値 0.0472(x≈1.430π)に近づく。
この結果の読み方を書いておく。主極大の高さは N2、幅は 1/N に比例するので、主極大の下の面積は N に比例し、原子1個あたりでは N によらない。結晶が大きくなるほど、散乱はちょうど q=2πh/a の点に鋭く集まる。これが次節の Laue 条件の1次元版である(図 1.3.1)。

図 1.3.1 原子 N=8 個の列の干渉関数 ∣L(q)∣2/N2(例題 1.3.1)。qa/2π が整数の所に高さ N2・幅 2π/(Na) の主極大が立ち、その間には N−2 個の小さな副極大がある。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com