結晶設定
充填率(原子充填率):
$$\eta = \frac{n \cdot \tfrac{4}{3}\pi r^3}{a^3}$$格子面間隔(ブラッグの法則に関連):$d_{hkl} = \dfrac{a}{\sqrt{h^2+k^2+l^2}}$(立方晶)
径分布関数:$g(r) = \dfrac{\rho(r)}{\rho_0}$(平均密度 $\rho_0$ との比)
$n$:単位格子内原子数、$r$:原子半径、$a$:格子定数(Å)、$(hkl)$:ミラー指数
単純立方・BCC・FCC・HCP・ダイヤモンド・NaCl・グラフェン(2D)の結晶構造をアイソメトリックビューで回転アニメーション表示。配位数・充填率・径分布関数g(r)を計算します。
充填率(原子充填率):
$$\eta = \frac{n \cdot \tfrac{4}{3}\pi r^3}{a^3}$$格子面間隔(ブラッグの法則に関連):$d_{hkl} = \dfrac{a}{\sqrt{h^2+k^2+l^2}}$(立方晶)
径分布関数:$g(r) = \dfrac{\rho(r)}{\rho_0}$(平均密度 $\rho_0$ との比)
$n$:単位格子内原子数、$r$:原子半径、$a$:格子定数(Å)、$(hkl)$:ミラー指数
結晶構造の基本は、原子が空間的に繰り返し配列した「単位格子」です。格子定数 $a$ を一辺の長さとする立方体を基本単位として、その中での原子の配置で構造が決まります。
$$ \text{充填率}= \frac{\text{単位格子内の原子の体積}}{\text{単位格子の体積}}= \frac{n \cdot \frac{4}{3}\pi r^3}{a^3}$$ここで、$n$は単位格子に含まれる実効的な原子数(例:BCCなら2)、$r$は原子半径、$a$は格子定数です。充填率は材料の密度や加工性に直結します。
最近接原子間距離 $d$ は、構造によって異なる式で表されます。これは原子同士がどのくらい近づいているかを示し、結合の強さに関わります。
$$ \text{BCC: }d = \frac{\sqrt{3}}{2}a \quad \text{FCC: }d = \frac{\sqrt{2}}{2}a $$BCCでは立方体の中心と頂点の原子が最近接(体対角線の半分)、FCCでは面心と頂点の原子が最近接(面の対角線の半分)です。シミュレーターで「最近接距離」を変えると、対応する格子定数 $a$ が変化します。
金属材料設計:自動車のボディに使われる鋼板は、主にBCC構造の鉄(フェライト)です。強度と加工性のバランスを取るために、炭素などの元素を添加し、結晶粒の大きさ(シミュレーターの「表示範囲」に相当)を制御します。CAEでは、この微細構造をモデル化して衝突安全性をシミュレーションします。
半導体デバイス:スマートフォンの頭脳であるCPUは、シリコンのダイヤモンド構造が基礎です。この規則正しい格子に、わずかに不純物(ドーパント)を添加することで、電気を通しやすくしたり遮断したりする性質(半導体特性)を生み出しています。
電池材料開発:リチウムイオン電池の正極材料(例:コバルト酸リチウム)は、層状の結晶構造(HCPに近い)を持ちます。リチウムイオンがその層の間を出入りすることで充放电が行われます。充填率やイオンの通り道(格子間隙)の設計が、電池容量や寿命を決めます。
軽量高強度材料:航空機の機体に使われるアルミニウム合金は、FCC構造をベースとしています。FCCはすべり面が多く加工性が良いため、複雑な形状に成形できます。ここに銅やマグネシウムなどを添加して析出強化させ、軽量でありながら高い強度を実現しています。
まず、「原子は硬いボール」という単純化モデルをそのまま信じすぎないでください。NovaSolverで表示される球体は、原子の「電子雲」の広がりを便宜的に半径で表したもの。実際の化学結合では、電子が共有されたり、軌道が混ざり合ったりするので、単純な幾何学的接触とは異なります。例えば、充填率が低いダイヤモンド構造が極めて硬いのは、この「硬いボール」モデルだけでは説明できず、強い共有結合の方向性が鍵です。
次に、パラメータを独立に変えられると誤解しがちです。ツールでは「最近接距離」スライダーを動かすと原子半径が変わり、格子定数も連動して変わります。しかし実材料では、格子定数は元素の種類でほぼ決まり、自由に変えられるものではありません。例えば、純鉄のBCC構造(フェライト)の格子定数は約0.286nm。ここに直径の大きいモリブデン原子を添加すると、無理やり格子を「引き伸ばす」ことで材料を強化します。ツールでBCCの半径を大きくして歪みを作る操作は、まさにこの固溶強化のイメージを理解するのに役立ちます。
最後に、「表示範囲」を単位格子だけに見ていると全体像を見失う点に注意。実在する材料は、この単位格子が数億〜数兆個も集まった「多結晶体」で、それぞれの結晶粒の向き(方位)がバラバラです。CAEで材料の異方性(方向による強度の違い)を評価する時は、この集合体としての振る舞いをシミュレーションします。ツールで表示範囲を広げて格子を繰り返し表示するのは、この「単結晶」から「多結晶」へのイメージを掴む第一歩です。
Cu(面心立方FCC、a=3.615Å)をセル数2×2×2で解析した場合:原子数32、配位数12、充填率74.05%。隣接原子間距離は3.615/√2≈2.556Åで計算され、RDFピークは2.55Å、3.615Å、4.43Åに出現。これらは{111}面間距離1.21Åの整数倍関係を確認できます。