X線回折(XRD)計算 戻る
材料科学 / 結晶解析

X線回折(XRD)計算ツール

結晶構造(FCC/BCC/SC/HCP/ダイヤモンド)とX線波長を選び、許容反射・d間隔・2θ・Scherrer幅をリアルタイム計算。回折パターンを可視化。

結晶構造パラメータ
結晶構造
格子定数 a (Å)3.61
X線源
X線源 / λ (Å)
Scherrer ブロード化
結晶子サイズ L (nm)30

基本式

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)

最強ピーク 2θ (°)
d-spacing (Å)
有効ピーク数
#hkld (Å)2θ (°)B (°)強度

X線回折(XRD)計算ツールとは

🧑‍🎓
「X線回折で結晶構造がわかる」って聞くけど、どうやって見分けるんですか?例えば鉄とアルミの回折パターンって、このツールで再現できますか?
🎓
ざっくり言うと、回折ピークが現れる角度の組み合わせが、結晶構造の「指紋」になるんだ。例えば、体心立方格子(BCC)の鉄と面心立方格子(FCC)のアルミでは、現れるピークの数と位置が全然違う。ツールの「結晶構造」をFCCとBCCで切り替えて、同じ格子定数で計算してみると、その違いが一目瞭然だよ。
🧑‍🎓
え、そうなんですか!「許容反射」って欄がありますけど、これが現れるピークの条件なんですね。でも、なぜFCCだと(100)面の反射は現れないんですか?
🎓
良いところに気づいたね。それは単位胞の中の原子の配置、つまり「構造因子」がゼロになるからだ。FCC構造では、h, k, lが全部奇数か全部偶数の時だけ回折強度が生まれる。だから(100)(1,0,0)は許されない。ツールで「結晶構造」をSC(単純立方)にすると、(100)が現れるのが確認できるよ。これが構造解析の基本だ。
🧑‍🎓
なるほど!下のグラフのピークの幅も、結晶子サイズLを変えると変わりますね。この「幅」から何がわかるんですか?
🎓
その通り。ピークがブロード(幅広)になるほど、結晶子(微結晶)のサイズが小さいことを意味する。これがScherrer(シェラー)式の応用だ。ツールで「結晶子サイズ L」のスライダーを、例えば100nmから10nmに小さくしてみて。回折角2θが大きい高角度のピークほど、幅が大きく広がるのが観察できるはずだ。実務では、ナノ材料の粒子径評価に頻繁に使われる手法さ。

物理モデルと主要な数式

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データの定量解析ができるようになる。ツールは、その結果を解釈するための「第一原理」を体感するための、最高の入り口なんだ。