逆格子空間

提供: testwiki
ナビゲーションに移動 検索に移動

本講座では、固体物理学を学ぶための必須知識である基本格子ベクトル、基本逆格子ベクトル、逆格子空間の3項目について説明する。

基本格子ベクトル

結晶の中にある原子の数は1023をゆうに超える膨大な数ある。物理を記述するため、この原子一つ一つを(x,y,z)のように座標で対応付けるのは計算上現実的ではない。

そこで固体物理学では、(理想的な)結晶中の原子が3次元方向に規則正しく周期的に並んでいることに着目し、単位胞(単位格子とも)という周期の最小単位のなかでの原子配置を記述する基本格子ベクトルだけを定義する。こうすることで、結晶中のあらゆる位置の原子の座標は、「基本格子ベクトル+単位胞の平行移動」ですべて表されることとなる。このことを、「結晶が並進対称性をもつ」と表現することがある。

数式で表すならば、単位胞を構成する基本格子ベクトル𝐚1,𝐚2,𝐚3(3次元なのでベクトルは3つ)を用いることで、任意の原子位置は

𝐑=n1𝐚1+n2𝐚2+n3𝐚3

と表現できる。n1,n2,n3は任意の定数である。

左から単純立方格子、面心立方格子、体心立方格子

以下に、固体物理ではポピュラーな3つの立方格子である単純立方格子(SC)、面心立方格子(FCC)、体心立方格子(BCC)(右図を見よ)それぞれの基本格子ベクトルを表す。一辺の長さ(格子定数)はすべてaとする。

SC格子:

𝐚1=a(1,0,0),𝐚2=a(0,1,0),𝐚3=a(0,0,1)

FCC格子:

𝐚1=a(12,12,0),𝐚2=a(0,12,12),𝐚3=a(12,0,12)

BCC格子:

𝐚1=a(1,0,0),𝐚2=a(0,1,0),𝐚3=a(12,12,12)

基本逆格子ベクトル

基本格子ベクトル𝐚1,𝐚2,𝐚3から

𝐚i𝐛j=2πδij

で定義されるベクトル𝐛1,𝐛2,𝐛3基本逆格子ベクトルと呼ぶ。ここで、δijはクロネッカーのデルタといい、i=jとなる場合にのみ1、それ以外では0となる関数である。これを明示的に表せば

𝐛1=2π𝐚2×𝐚3𝐚1𝐚2×𝐚3

𝐛2=2π𝐚3×𝐚1𝐚1𝐚2×𝐚3

𝐛3=2π𝐚1×𝐚2𝐚1𝐚2×𝐚3

である。この3式が定義式を満たすことを確認してみるとよい。

ここまで基本逆格子ベクトルを天下り的に与えたが、実際これが何を意味するのかを説明しよう。

まず𝐚i𝐛j=2πδijを読み解くと、「𝐛1は(1)𝐚2,𝐚3に垂直かつ(2)𝐚1との内積が2πとなるベクトル」などである。(1)はさておき、(2)はどういう意味なのか。これは波動論において波数がk=2πλで与えられることと結びつけるとよい。すなわち、基本格子ベクトルが実空間(距離空間)における位置を表すベクトルだとすれば、基本逆格子ベクトルは波数空間において基本格子ベクトルに対応するベクトルである。

これらにより実空間と波数空間が対応付けられ、自在に行き来できるようになる(はずである)。

簡単な外積の練習として、SC格子、FCC格子、BC格子それぞれの基本逆格子ベクトルを求めてみよ。仮に基本格子ベクトルを前項のように定義すれば以下のようになるはずである。

SC格子:

𝐛1=2πa(1,0,0),𝐛2=2πa(0,1,0),𝐛3=2πa(0,0,1)

FCC格子:

𝐛1=2πa(1,1,1),𝐛2=2πa(1,1,1),𝐛3=2πa(1,1,1)

BCC格子:

𝐛1=2πa(1,0,1),𝐛2=2πa(0,1,1),𝐛3=2πa(0,0,2)

逆格子空間

逆格子空間とは、基本逆格子ベクトルで張られる空間である。この空間上の一点を逆格子ベクトル𝐆という。明示的に書けば、

𝐆=n1𝐛1+n2𝐛2+n3𝐛3

と表される(n1,n2,n3は整数)。

とはいえ、これが一体何の役に立つのだろう。その説明のため、金属物質中の電子密度について考えよう。

理想的な結晶構造を持つ金属であれば、その原子配列は周期的であるため、その周りにある電子も各軸方向𝐚1,𝐚2,𝐚3それぞれに周期的な分布となるはずである。

電子密度n(𝐫)は並進操作𝐓=u1𝐚1+u2𝐚2+u3𝐚3に対して不変である:

n(𝐫+𝐓)=n(𝐫)

すなわち、結晶の軸がどのように平行移動しても分布自体は変換されないということである。

この電子密度が周期Lをもつとき、フーリエ級数に展開した形は、

n(𝐫)=p1,p2,p3>0neexp(2πp1La1)exp(2πp2La2)exp(2πp3La3)

と表される。ここで、p1,p2,p3は整数であり、neはフーリエ係数[注 1]。また、𝐫=(a1,a2,a3)の各方向に対して展開をしていることに注意。

実際、n(𝐫+𝐓)を計算して並進操作に対して不変であることを確かめてみるとよい。

先ほどのn(𝐫)の形は、以下のようにまとめられる。

n(𝐫)=𝐏neexp(i2πL𝐏𝐫)

つまり、並進操作𝐓に対してn(𝐫)を不変とするようなベクトル𝐏=(p1,p2,p3)が必要である。

ここで思い出してもらいたいのが、前々節で述べた基本逆格子ベクトルの定義𝐚i𝐛j=2πδijである。

𝐫は基本格子ベクトルの組み合わせであることから、同様に𝐏を逆格子ベクトルの組み合わせとすればよいと分かる:

𝐏=p1𝐛1+p2𝐛2+p3𝐛3

これは本節の初めに述べた逆格子ベクトル𝐆そのものである。

まとめると、電子密度は

n(𝐫)=𝐆nGexp(i2πL𝐆𝐫)

と書ける。これを読み解くと、逆格子ベクトルは電子密度と一対一対応していると言える。この一対一対応こそ、逆格子ベクトルと逆格子空間が役に立つ理由である。

正確に言えば、逆格子空間の各点(=逆格子ベクトル)は、電子密度のような結晶内の周期的な物理量と一対一対応しているということである。

この性質の有用性が分かる事例を次節に紹介する。

逆格子空間の利用例 - ブラッグの法則

前節の一対一対応の原理を用いたのが以下に解説するブラッグの法則である。ブラッグ回折やブラッグ反射などとも呼ばれる。

ブラッグの法則は、結晶中に入射した光(X線など)が散乱される際に、ある格子定数L(およびその整数倍)とその入射角の組が一意に定められるという法則である。

波数𝐤iが結晶格子内に侵入し、原子に散乱されて波数𝐤fとなり出ていく様子。

簡単のため右図のように、二次元方向のみ考える。また、弾性散乱のみと考えるものとする。

波数𝐤iが結晶内で散乱されて波数𝐤fの波となるとき、散乱波の振幅が入射波と同位相になって強められるためには、どういった波数の条件が満たされなければならないだろうか?

前節でも述べたように、結晶の周期に伴う電子の分布n(𝐫)はその対応する逆格子ベクトル𝐆を用いて

n(𝐫)=𝐆nGexp(i2πL𝐆𝐫)

と表されるのだった。ここで、式内の次元を考えれば、𝐆は1/(長さ)の次元、すなわち波数の次元を持つことが分かる。

入射波𝐤iがその電子(分布)に散乱され、散乱波𝐤fになるということは、

𝐤f=𝐤i+𝐆

が成り立つということである[注 2]この式がブラッグの法則そのものである

この法則が具体的にどう有用なのかまだ不透明なので、先ほどの図を使ってもっと具体的に考えよう。

この図において、横方向に𝐚1、上方向に𝐚2の単位ベクトルが向いているとする。 ここで、逆格子ベクトル𝐆は図中上方向に伸びる。その大きさは、

|𝐆|=|h𝐛1+k𝐛2+l𝐛3|

=k|𝐛2|

=k2πL

と表される。ここで、kは整数であり、強め合う二つの散乱波の反射する結晶の"レイヤー"の枚数に一致する(つまり、k=1なら間隔Lで隣り合う散乱波同士の干渉を意味する)。

この式から、結晶間の間隔L、すなわち格子定数が導かれる。

実験的に言えば、入射波と散乱波の差(=散乱ベクトル)を測定することで、その結晶の逆格子ベクトルが分かり、そこから結晶の格子定数が測定できるということである。

特に今は弾性散乱を考えているので|𝐤i|=|𝐤f|であり、すなわち散乱波の結晶面に対する角度を測定するだけで結晶の格子定数が計算可能ということである。

逆格子空間はこのように波(波数)を意味する(フーリエ)空間であり、波の散乱を考える際には非常に便利だということがお判りいただけただろうか。

注釈

  1. フーリエ係数自体は複素数だが、電子密度n(𝐫)は実関数としている。本来はp1p3それぞれに対してフーリエ係数が存在することからpp番目の項の和を実数とする仮定が必要である。しかし、本講座ではこの辺りの議論を省略してneと括っている。
  2. この説明は物理学的にはあまり正しくない。散乱振幅を考えるべきであるが、直感的な理解の助けとなるようこのような説明とした。