結晶系の選択
結晶構造
ミラー指数 (hkl)
計算結果
d面間隔
—
ブラッグ角 θ (Cu Kα)
—
2θ
—
配位数
—
充填率
—
原子/単位格子
—
理論式
立方晶系 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線回折（XRD）は材料の結晶構造・格子定数・残留応力を非破壊で測定する標準手法。CAEモデルの材料パラメータ（弾性係数・熱膨張係数）も結晶構造から予測されます。
X線粉末回折パターン（Cu Kα, λ=1.54Å）
結晶構造と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) で指定される結晶面の面間隔 d は、格子定数 a を用いて以下のように計算されます。この値は、その面の原子の密度や、X線回折の条件に直接関わります。
$$d_{hkl}= \frac{a}{\sqrt{h^2+k^2+l^2}}$$
$d_{hkl}$: (hkl)面の面間隔 [Å]
$a$: 立方晶の格子定数 [Å] (シミュレーターのパラメータ)
$h, k, l$: ミラー指数（整数）
X線回折のピークが現れる角度（ブラッグ角θ）は、ブラッグの法則によって与えられます。ここでは、X線源として広く使われる銅(Cu)のKα線の波長を用いています。
$$n\lambda = 2d\sin\theta \quad (\lambda_{\text{CuK}\alpha}=1.54\,\text{Å})$$
$n$: 回折の次数（通常は1）
$\lambda$: X線の波長 [Å]
$d$: 面間隔 [Å]
$\theta$: ブラッグ角 [度] (グラフの横軸は通常2θで表示されます)
この式から、面間隔dが求まれば、回折ピークの角度θが計算できます。
実世界での応用
新規材料の構造同定：研究室で新しく合成された金属間化合物やセラミックスの粉末X線回折(XRD)パターンを取得し、このシミュレーターのように既知の結晶構造（BCC, FCC, HCPなど）と格子定数を仮定して計算されたピーク位置と照合します。これにより、その材料がどのような結晶構造を持っているかを決定します。
合金の相分析：鉄鋼材料にはフェライト（BCC）やオーステナイト（FCC）など、温度や組成によって異なる結晶構造の相が存在します。XRDパターンのピークを解析することで、材料中にどの相が、どの程度の体積分率で存在するかを定量的に評価できます。
残留応力測定：加工や溶接を受けた部品では、結晶格子自体が歪み、面間隔dが変化します。この変化はXRDピークのわずかな角度シフトとして現れます。ピークシフト量を精密に測定することで、材料表面に残留している応力の大きさと方向を非破壊で評価できます。
CAE材料パラメータの推定：弾性係数や熱膨張係数などの材料物性は、結晶構造や原子間結合の強さに強く依存します。第一原理計算などのCAE手法では、このツールで扱うような原子配列（結晶構造と格子定数）を出発点として、これらの巨視的な材料パラメータを予測することがあります。
よくある誤解と注意点
まず、このツールで「ミラー指数は面の向きだけを表す」と思いがちだけど、実は面間隔も決めるってことを押さえよう。例えばFCC構造で(200)面と(100)面は平行だよね？でも、このツールで見ると(200)面の方が(100)面よりずっと密に原子が並んでいる（面間隔が半分）のがわかる。X線回折では(200)面として別のピークが出るんだ。ここを混同すると、ピークの指数付けを間違えるぞ。
次に、「充填率」は原子の大きさを変えても変化しないって点。ツールで原子半径をスライダーで動かすと、原子が重なって見えるけど、充填率の計算式はあくまで「剛体球が最密充填したとき」の理論値だから、半径を変えても表示された値は変わらない。実際の材料では原子半径の概念自体が曖昧だから、これは理想的なモデルでの話だよ。
最後に、実務でXRDパターンを見るときの落とし穴。このシミュレーターのグラフは理想的な完全結晶のパターンだ。実際の材料には格子欠陥や微細結晶、残留応力があるから、ピークがブロードニング（幅広がり）したり、シフトしたりする。例えば工具鋼の焼入れ後は、マルテンサイト相の発生でピークが大きく広がる。ツールのシャープなピークだけを鵜呑みにしないようにね。
関連する工学分野
このツールで計算される「面間隔d」と「ブラッグ角θ」は、CAEを含む様々な工学分野の基礎データとして使われているんだ。
まずは半導体デバイス工学。シリコンウェハーの結晶方位はデバイス特性に直結する。ツールで(100)面と(111)面の原子配置の違いを見てみると、(111)面の方が原子密度が高い（面間隔が狭い）のがわかるよね。この密度の違いがエッチング速度や酸化速度に影響し、微細加工プロセス設計の根拠になるんだ。
次に金属疲労・破壊解析。結晶粒のすべり面は、最も面間隔の広い密面（例えばFCCなら{111}面）で起こりやすい。ツールでいろんな面の間隔を比較すれば、なぜ特定の面で変形が進むのかイメージできる。有限要素法(FEM)で結晶塑性モデルを使う時は、このような結晶学的な情報が材料パラメータとして入力されるぞ。
あとは電池材料開発も面白い応用例だ。リチウムイオン電池の正極材料（例えば層状岩塩型）では、リチウムイオンが結晶の特定の面間（例えば(003)面間）を移動する。この拡散経路と面間隔は相関が強く、充放電特性を左右する。XRDで格子定数の変化を追うことで、充電中の材料の状態をモニターできるんだ。
発展的な学習のために
このツールに慣れたら、次は「なぜそのピークが現れる（or 消える）のか？」を数式レベルで追いかけてみよう。例えばBCC構造で(100)面の回折ピークがなぜ消えるのか？それは構造因子の計算に立ち戻れば理解できる。BCCの場合は、単位格子の中心と角の原子からの散乱波の位相差が半波長ずれて打ち消し合う条件になるんだ。数式で書くと、ミラー指数の和(h+k+l)が奇数の時に散乱強度がゼロになる。この理屈がわかれば、FCCやNaCl型でのピーク出現則も自分で導けるようになる。
数学的には、逆格子空間の概念を学ぶと全てがつながる。実空間の面間隔$$d_{hkl}$$は、逆格子空間におけるベクトル$$\vec{g}_{hkl}$$の長さの逆数（$$d = 2\pi / |\vec{g}|$$）に対応する。このツールで格子定数aを大きくすると、実空間ではユニットセルが大きくなるけど、逆空間では逆格子ベクトルが短くなり、面間隔dが増加する。この双対的な関係を理解すると、電子線回折(SAED)パターンなどの解析も格段に楽になるぞ。
推奨する次のステップは、「格子定数が一定でない」世界に進むこと。実際の合金は固溶体なので、組成変化に伴って格子定数が連続的に変化する（ヴェガードの法則）。また、複相材料の場合は、ツールで異なる構造・格子定数を設定した2つのシミュレーションパターンを重ね合わせて、現実のXRDパターンを解釈する練習をしてみると、一気に実践力が上がるよ。