基本式
Bragg の法則:$2d_{hkl}\sin\theta = \lambda$
面間隔(立方晶):$$d_{hkl}=\frac{a}{\sqrt{h^2+k^2+l^2}}$$
Scherrer式:$$B = \frac{K\lambda}{L\cos\theta}$$(K=0.9)
| # | hkl | d (Å) | 2θ (°) | B (°) | 強度 |
|---|
結晶構造(FCC/BCC/SC/HCP/ダイヤモンド)とX線波長を選び、許容反射・d間隔・2θ・Scherrer幅をリアルタイム計算。回折パターンを可視化。
Bragg の法則:$2d_{hkl}\sin\theta = \lambda$
面間隔(立方晶):$$d_{hkl}=\frac{a}{\sqrt{h^2+k^2+l^2}}$$
Scherrer式:$$B = \frac{K\lambda}{L\cos\theta}$$(K=0.9)
| # | hkl | d (Å) | 2θ (°) | B (°) | 強度 |
|---|
X線回折の基本はBragg(ブラッグ)の法則です。結晶内の平行な原子面で反射したX線が強め合う条件を表し、回折角θと面間隔dの関係を決めます。
$$2d_{hkl}\sin\theta = n\lambda$$ここで、$d_{hkl}$は格子面間隔(Å)、$\theta$はブラッグ角(度)、$n$は反射次数(通常1)、$\lambda$はX線の波長(Å)です。この式から、回折が観測される角度2θを計算できます。
立方晶系の結晶では、面間隔$d_{hkl}$はミラー指数(hkl)と格子定数aを用いて以下のように計算されます。
$$d_{hkl}=\frac{a}{\sqrt{h^2+k^2+l^2}}$$$a$は格子定数(Å)、$h, k, l$はミラー指数です。この式により、異なる(hkl)面に対応する回折角を求めることができます。六方晶(HCP)などでは、もう一つの格子定数cを含む少し複雑な式になります。
回折ピークの半値幅(広がり)Bから、結晶子の平均サイズLを見積もるのがScherrer(シェラー)の式です。
$$B = \frac{K\lambda}{L\cos\theta}$$$B$は回折ピークの半値幅(ラジアン)、$K$は形状因子(通常0.9)、$\lambda$はX線波長(Å)、$L$は結晶子サイズ(Å)、$\theta$はブラッグ角です。ピークが広がるほどLは小さく、結晶子が微細であることを示します。
新材料開発・同定:合成した未知の粉末試料のXRDパターンを測定し、データベースと照合したり、このツールのような計算結果と比較することで、物質の同定や結晶構造の確認を行います。例えば、新しい電池材料の純度評価に不可欠です。
ナノ材料の特性評価:Scherrer式を用いて、ナ粒子や薄膜の結晶子サイズを評価します。触媒材料では、粒子サイズが性能に直結するため、XRDによるサイズ測定は標準的な手法です。
残留応力測定:結晶に応力がかかると格子定数が変化し、回折ピークの角度がシフトします。このシフト量を精密に測定することで、材料内部の残留応力を非破壊で評価できます。溶接部の検査などで応用されます。
薄膜の構造解析:基板上に成長した薄膜の結晶方位や格子定数を評価します。半導体デバイス製造では、エピタキシャル成長した薄膜の品質管理にXRDが広く用いられています。
このツールを使い始めるとき、いくつかハマりやすいポイントがあるから気をつけてね。まず「格子定数aと結晶子サイズLを混同する」ケース。格子定数は原子同士の距離(単位胞の大きさ)で、Å(オングストローム)単位。一方、結晶子サイズLは、単結晶のように秩序だった領域がどこまで続いているかの「広がり」で、通常はnm(ナノメートル)単位だ。例えば、格子定数が約3.6Åのアルミニウムでも、ナノ粒子ならLは10nm(=100Å)程度と、全く異なるスケールなんだ。
次に「Scherrer式は万能ではない」という点。ツールでLを小さくするとピーク幅が広がるのは確かだけど、実測の幅は「結晶子サイズ」以外にも、装置の分解能や結晶の歪み(ひずみ)の影響を強く受ける。だから、実験データからLを逆算する時は、装置関数の補正や、複数のピークから歪みの影響を分離する「Williamson-Hallプロット」といった手法が必要になることが多いよ。
最後に、「計算される2θは理想値」であることを忘れないで。実際の測定では、試料の吸収や表面の粗さ、オフセットなどでピーク位置が数0.1°ずれることはザラ。同定作業では、ピークの「相対的な位置関係のパターン」を重視しよう。例えばFCC構造で、計算上(111)ピークと(200)ピークの2θの比がどうなるか、といった大局的な見方を身につけることが実務では役立つんだ。
XRDの計算が活躍するのは、何も材料科学だけじゃない。隣接する多くの工学分野で、結晶情報がキーになるんだ。半導体工学では、エピタキシャル成長した薄膜の結晶品質評価や、歪みシリコン中の応力測定に高分解能XRD(HR-XRD)が必須だ。ツールで計算するピーク位置が、基板と薄膜の格子不整合による「ピークのシフト」として観測され、そこから応力が定量化される。
触媒化学の分野でも重要だ。担持金属触媒では、ツールで想定するような単純な構造ではなく、非常に小さな結晶子(Lが数nm)が表面に分散している。得られるXRDパターンはブロードなピークしかなく、まさにScherrer式が威力を発揮して、活性サイトの大きさ(粒子径)を推定するのに使われる。例えば自動車の排ガス浄化触媒の性能は、この粒子径と強く相関しているんだ。
さらに残留応力測定や相変態の追跡にも応用される。例えば、焼入れした鋼はマルテンサイト相が生じるが、この相は体心正方晶(BCT)となり、通常のBCC鉄とはピーク位置が微妙に異なる。ツールで格子定数を少し変えて計算した結果と実測パターンを比較することで、変態の進行度合いをモニターできる。このように、計算ツールは「何が起こっているのか」を解釈するための強力な羅針盤になるんだ。
ツールの操作に慣れたら、次は背後にある理論を深掘りしてみよう。最初のステップは「逆空間の概念」に触れること。我々がツールで計算している回折角は、実は「逆格子ベクトル」の物理的表現なんだ。ブラッグの法則 $2d\sin\theta = \lambda$ は、逆空間では散乱ベクトル $\vec{k}-\vec{k_0}$ が逆格子点 $\vec{g}_{hkl}$ に一致する条件、つまり $\vec{k}-\vec{k_0} = \vec{g}_{hkl}$ と非常にシンプルに書ける。この視点を得ると、なぜ結晶構造によって反射の条件(構造因子)が変わるのかが、原子の配置をフーリエ変換するという操作として理解しやすくなるよ。
数学的には、フーリエ級数・変換の基礎を学ぶのが近道だ。結晶からの回折強度は、単位胞内の電子密度分布のフーリエ変換の自乗(つまり、構造因子の絶対値の2乗)で与えられる。ツールで「許容反射」と出ているのは、この構造因子がゼロでない(hkl)の組み合わせなんだ。例えばダイヤモンド構造で(200)が禁止される理由を、基底に2個の原子があることからフーリエ変換で説明できるようになれば一人前だ。
次の具体的なトピックとしては、「リートベルト解析」を調べてみることを勧める。これは、ツールでやっているような個々のピークの計算を超えて、測定された全体的な回折プロファイルを、格子定数、結晶子サイズ、微ひずみなど多数のパラメータで同時にフィッティングする手法だ。市販のソフトウェアの根幹をなす技術で、これを学べばXRDデータの定量解析ができるようになる。ツールは、その結果を解釈するための「第一原理」を体感するための、最高の入り口なんだ。