結晶構造・ミラー指数ビジュアライザー 戻る
材料科学シミュレーター

結晶構造・格子定数・ミラー指数ビジュアライザー

SC・BCC・FCC・HCP・NaCl型結晶のユニットセルを可視化し、ミラー指数(hkl)からd面間隔・ブラッグ角・X線回折パターンをリアルタイム計算。

結晶系の選択
結晶構造
格子定数 a
Å
ミラー指数 (hkl)
計算結果
d面間隔
ブラッグ角 θ (Cu Kα)
配位数
充填率
原子/単位格子
結晶
X線粉末回折パターン(Cu Kα, λ=1.54Å)
理論・主要公式

立方晶系 d面間隔:

$$d_{hkl}= \frac{a}{\sqrt{h^2+k^2+l^2}}$$

ブラッグの法則:

$$n\lambda = 2d\sin\theta \quad (\lambda_{\text{CuK}\alpha}=1.54\,\text{Å})$$

結晶構造とX線回折とは

🙋
このシミュレーターで「ミラー指数」を変えると、3Dモデルの赤い面が動きますけど、あれは何を表しているんですか?
🎓
あの赤い面は、君が入力したミラー指数 (hkl) で表される「原子が並んでいる面」だよ。例えば(100)なら立方体の各面に平行な面だし、(110)なら立方体を斜めに切る面になる。上のスライダーでh,k,lを1や2に変えてみると、面の傾きがどう変わるか、すぐに体感できるぞ。
🙋
なるほど!で、下の「X線回折パターン」のグラフにピークが出てきますが、あれは(100)面とかで反射した結果なんですか?
🎓
その通り!X線が結晶に入射すると、その(100)面や(110)面で反射する。ブラッグの法則を満たす角度θでだけ強め合ってピークが現れるんだ。シミュレーターでは、右側で選んだ結晶構造(例えばBCC)と格子定数aから、全ての可能な面の間隔dを計算し、ピーク位置をプロットしている。構造をFCCに変えてみると、現れるピークの組み合わせが変わるのがわかるよ。
🙋
現場では、この逆のことをやるんですよね?実験でピークを測って、結晶構造や格子定数を決める。
🎓
鋭いね!実務ではまさにその逆プロセスが重要だ。例えば新しく合成した合金のXRDパターンを取って、このツールのように「FCCでa=3.6Åだと、ピークはここに来るはず」とシミュレーション結果と照合する。これが材料同定の基礎だ。ツールの格子定数aを少しずつ動かしてみると、ピーク位置がどうシフトするか確認できる。残留応力で格子が歪むと、まさにこうしてピークが移動するんだ。

よくある質問

ブラッグの法則を満たすために、入力した(hkl)と格子定数a、X線波長λの組み合わせでsinθが1以下になる必要があります。また、構造因子がゼロとなる禁制反射(例:BCCの(100)面)ではピークは出現しません。画面のエラーメッセージをご確認ください。
異なります。本ツールで表示されるθはブラッグ角であり、実験のX線回折パターンで横軸に取られる2θ(入射X線と回折X線のなす角)の半分です。ピーク位置を実験と比較する際は、表示値を2倍にしてご利用ください。
NaCl型構造では、ナトリウムイオン(Na⁺)と塩化物イオン(Cl⁻)を異なる色(例:青と緑)で表示しています。これにより、両イオンが交互に配置した岩塩構造を視覚的に確認でき、ミラー指数で指定した面がどのイオン種を切断するかも把握できます。
本ツールでは、原子の表示サイズは格子定数に依存せず、各元素のイオン半径または共有結合半径に基づく固定値で描画しています。これは結晶構造の種類や原子の相対的な大きさをわかりやすくするためです。格子定数の変更は面間隔dや回折角度の計算にのみ反映されます。

実世界での応用

新規材料の構造同定:研究室で新しく合成された金属間化合物やセラミックスの粉末X線回折(XRD)パターンを取得し、このシミュレーターのように既知の結晶構造(BCC, FCC, HCPなど)と格子定数を仮定して計算されたピーク位置と照合します。これにより、その材料がどのような結晶構造を持っているかを決定します。

合金の相分析:鉄鋼材料にはフェライト(BCC)やオーステナイト(FCC)など、温度や組成によって異なる結晶構造の相が存在します。XRDパターンのピークを解析することで、材料中にどの相が、どの程度の体積分率で存在するかを定量的に評価できます。

残留応力測定:加工や溶接を受けた部品では、結晶格子自体が歪み、面間隔dが変化します。この変化はXRDピークのわずかな角度シフトとして現れます。ピークシフト量を精密に測定することで、材料表面に残留している応力の大きさと方向を非破壊で評価できます。

CAE材料パラメータの推定:弾性係数や熱膨張係数などの材料物性は、結晶構造や原子間結合の強さに強く依存します。第一原理計算などのCAE手法では、このツールで扱うような原子配列(結晶構造と格子定数)を出発点として、これらの巨視的な材料パラメータを予測することがあります。

よくある誤解と注意点

まず、このツールで「ミラー指数は面の向きだけを表す」と思いがちだけど、実は面間隔も決めるということを押さえよう。例えばFCC構造で(200)面と(100)面は平行だよね?でも、このツールで見ると(200)面の方が(100)面よりずっと密に原子が並んでいる(面間隔が半分)のがわかる。X線回折では(200)面として別のピークが出るんだ。ここを混同すると、ピークの指数付けを間違えるぞ。

次に、「充填率」は原子の大きさを変えても変化しないって点。ツールで原子半径をスライダーで動かすと、原子が重なって見えるけど、充填率の計算式はあくまで「剛体球が最密充填したとき」の理論値だから、半径を変えても表示された値は変わらない。実際の材料では原子半径の概念自体が曖昧だから、これは理想的なモデルでの話だよ。

最後に、実務でXRDパターンを見るときの落とし穴。このシミュレーターのグラフは理想的な完全結晶のパターンだ。実際の材料には格子欠陥や微細結晶、残留応力があるから、ピークがブロードニング(幅広がり)したり、シフトしたりする。例えば工具鋼の焼入れ後は、マルテンサイト相の発生でピークが大きく広がる。ツールのシャープなピークだけを鵜呑みにしないようにね。

使い方ガイド

  1. val-aで結晶系(SC/BCC/FCC)を選択し、sl-aで格子定数a(Å単位、例:Fe BCC時a=2.87)を入力
  2. inp-h、inp-k、inp-lにミラー指数を入力(例:(110)面の場合h=1,k=1,l=0)
  3. シミュレータが自動計算:d面間隔=a/√(h²+k²+l²)、ブラッグの法則2dsinθ=λ(Cu Kα λ=1.5406Å)からθと2θを算出
  4. 結晶構造の配位数・充填率・単位格子当たりの原子数も同時表示

具体的な計算例

Cu FCC結晶(a=3.615Å)の(200)面をX線回折解析:ミラー指数h=2,k=0,l=0を入力するとd=(3.615/√4)=1.808Å、Cu Kα線(λ=1.5406Å)のブラッグ角θ=56.3°、2θ=112.6°と計算されます。FCC構造は単位格子当たり4原子、充填率74%、配位数12を示します。

実務での注意点