5種類の結晶構造の充填率・配位数・格子定数をリアルタイム計算。ミラー指数からブラッグ回折角まで一気に求めよう。
立方晶面間隔: $d_{hkl}= a/\sqrt{h^2+k^2+l^2}$
結晶構造の幾何学から、単位格子の一辺の長さ（格子定数 a）と原子半径 r の関係、および原子が占める体積の割合（充填率）を計算します。構造ごとに原子の配置が異なるため式が変わります。
$$a_{\text{SC}}= 2r, \quad a_{\text{BCC}}= \frac{4r}{\sqrt{3}}, \quad a_{\text{FCC}}= \frac{4r}{\sqrt{2}}$$SC（単純立方）: 立方体の辺上に原子が並ぶ。BCC: 体心と頂点の原子が接するのは体対角線方向。FCC: 面心と頂点の原子が接するのは面の対角線方向。
X線回折の条件はブラッグの法則で記述されます。立方晶では、ミラー指数 (hkl) で表される結晶面の間隔 d は格子定数 a から求められます。
$$2d_{hkl}\sin\theta = n\lambda, \quad d_{hkl}= \frac{a}{\sqrt{h^2 + k^2 + l^2}}$$$d_{hkl}$: 面間隔、$\theta$: ブラッグ角（回折角）、$n$: 回折次数（通常1）、$\lambda$: X線の波長。この条件を満たす角度 $\theta$ でのみX線が強め合い、回折ピークが観測されます。
材料開発・同定：未知の金属や合金のサンプルにX線回折実験を行い、得られた回折ピークの角度からブラッグの法則を用いて面間隔を逆算。既知の結晶構造のデータベースと照合することで、材料を同定したり、混ざり物（異相）を検出したりします。
半導体製造：シリコンウェハーはダイヤモンド構造をとります。成膜した薄膜の結晶性（単結晶か多結晶か）や格子ひずみを、X線回折プロファイルのピーク位置や幅から評価し、デバイス性能の管理に役立てています。
金属の強度設計：鉄（Fe）は温度によってBCC（常温、α鉄）とFCC（高温、γ鉄）で構造が変わります。この相変態を利用した熱処理（焼入れ・焼き戻し）により、鋼の強度や靭性を制御します。充填率の違いが変態時の体積変化に関わります。
電池材料研究：リチウムイオン電池の正極材料などは、リチウムの挿入・脱離に伴って結晶格子が膨張・収縮します。その格子定数の微小な変化を高精度X線回折で追跡し、電池の劣化メカニズムの解明やより安定な材料の開発に応用されています。
まず、このツールを使う上で気をつけてほしいのが「原子半径 r」と「格子定数 a」の関係だ。シミュレーター上で r を動かすと a も連動して変わるけど、これは原子が剛体球で、かつ隙間なく接している理想的な状態を仮定した計算だからね。実際の材料では、原子間の結合の性質（金属結合、共有結合など）によってこの関係は崩れる。例えば、実在する鉄（BCC）の格子定数から逆算した「原子半径」は、この幾何学的な式 $a = 4r / \sqrt{3}$ から求める値と完全には一致しないんだ。ツールはあくまで基本原理を理解するための「理想モデル」と心得よう。
次に、ブラッグ条件の計算で「回折角 θ」と表示される角度について。これは入射X線と結晶面のなす角（ブラッグ角）であって、よくある誤解は「検出器の位置（2θ）」と混同すること。実際のX線回折装置では、試料に対して入射角θ、検出器は2θの位置に移動して測定する。ツールでλ=0.154 nm (CuKα線)、FCCのa=0.3615 nm (Al)、hkl=(111)とするとθ≈19.3°と計算されるけど、実験プロトコルではこの角度を試料の傾き角として設定し、検出器は約38.6°の位置でピークを探すことになるんだ。
最後に、HCP（六方最密構造）の「充填率」がFCCとほぼ同じ約74%になる理由だ。3Dモデルを見ると原子の並びが全然違うのに、と思うよね。ここで理解すべきは、「最密充填」という状態が本質的に同じ空間充填効率を持つということ。FCCはABCABC…、HCPはABAB…という積層順序の違いだけなんだ。充填率や配位数（どちらも12）といった大局的な性質は同じになる。ツールで両方の構造を回転させて、最密面（FCCなら{111}面、HCPなら底面）がどう連なっているか観察してみると、この「積層の違い」が実感できるよ。
このツールで扱う計算は、CAEの中でも「マルチスケールシミュレーション」の入り口に直結する。原子レベルの結晶構造（ナノスケール）の情報は、その上のメゾスケール（結晶粒の集まり）やマクロスケール（部品全体）の材料挙動を決める基礎データになるんだ。例えば、分子動力学（MD）シミュレーションでは、ここで学ぶFCCやBCCを初期原子配置として設定し、変形や破壊の過程を計算する。ツールで確認した面間隔 $d_{hkl}$ は、MDのポテンシャル関数を設定する時の基準長さにも関連してくるぞ。
もう一つの大きな応用が残留応力測定だ。X線回折で得られるピーク位置は、結晶格子そのものが伸び縮みするとシフトする。この原理を利用して、溶接部や表面処理を施した部品の内部に残る応力を非破壊で評価できる。ツールで「格子定数 a」をわずかに増減させてみてくれ。例えば、aを0.360 nmから0.362 nmに0.5%増やすだけで、回折角θが明確にシフトするのが計算できるだろう？ これが実際の応力測定技術の根幹にある「格子ひずみの測定」そのものなんだ。
また、添加剤製造（3Dプリンティング）でもこの知識は生きる。金属粉末をレーザで溶融・凝固させるプロセスでは、急冷によって微細な結晶組織が生じる。この生成相がFCCかBCCか、あるいは非平衡相かは、製品の機械的特性を左右する。製造後の部品にX線回折を適用して結晶相を同定し、プロセス条件を最適化するのに、ブラッグの法則と結晶構造の知識が不可欠なんだ。
このツールに慣れたら、次のステップとして「逆格子空間」の概念を取り入れてみよう。ツールで計算している面間隔 $d_{hkl}$ の逆数 $1/d_{hkl}$ を考えると、実はこれが逆格子空間における点（逆格子点）までの距離になる。X線回折図形は、この逆格子空間の写像と解釈できるんだ。この視点を得ると、単なる角度計算だった回折条件が、幾何学的で直感的な操作に変わる。まずは、立方晶の $1/d_{hkl} = \sqrt{h^2+k^2+l^2} / a$ という式を、実空間の格子定数 a と逆空間の関係として捉え直す練習から始めると良い。
数学的には、ベクトルと線形代数の基礎が役に立つ。結晶面のミラー指数 (hkl) は、実はその面に垂直な法線ベクトルの成分表示と深く関係している。また、HCP構造で出てくる4指数表示 (hkil) の「i」が $i = -(h+k)$ となる理由を理解するには、六方晶系の基底ベクトルを理解する必要がある。ツールでHCPを表示しながら、底面の六角形の軸（a1, a2）とその和の方向を考えてみると、この対称性が見えてくるはずだ。
実務に近い次のトピックとしては、「回折ピークの幅」に注目することを勧める。ツールは完全な単結晶を想定したシャープな回折条件だけを計算するが、実際の材料では結晶粒が小さい（微結晶）場合や格子に歪みがある場合、ピークはブロード化する。このブロードニングの解析（シェラー式やウィリアムソン・ホールプロット）を通じて、材料の微細構造を定量的に評価する手法が広く使われている。まずは、完全な結晶と不完全な結晶の回折プロファイルの違いをイメージ図で比較してみると、学びが深まるよ。