結晶格子ビューワー 戻る
材料科学 / 結晶学

結晶格子ビューワー

単純立方・BCC・FCC・HCP・ダイヤモンド・NaCl・グラフェン(2D)の結晶構造をアイソメトリックビューで回転アニメーション表示。配位数・充填率・径分布関数g(r)を計算します。

結晶設定

結晶構造
格子定数 a (Å)
Å
表示セル数1×1×1
回転速度
計算結果
配位数
充填率
最近接距離
原子数/UC
回転速度
1.0
格子
動径分布
理論・主要公式

充填率(原子充填率):

$$\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」を切り替えてみて。原子の詰まり方が大きく異なるでしょ?この詰まり方の違いが、鉄が硬くなったり柔らかくなったりする理由なんだ。
🙋
え、そうなんですか?「充填率」って表示されてますけど、これが高いとどうなるんですか?
🎓
充填率が高いほど、原子同士が密に詰まっているということだ。FCCやHCPは約74%で最密構造。だからアルミニウム(FCC)は加工しやすいんだ。逆に、ダイヤモンド構造を選んでみて?充填率は約34%と空隙が多いだけど、その隙間のない強い結合が世界一硬い理由なんだよ。
🙋
なるほど!でも、この「最近接距離」のスライダーを動かすと、原子の大きさが変わるみたいです。これは何を変えてるんですか?
🎓
いいところに気づいたね!それは原子の「半径」を変えているんだ。実務では、異なる元素を混ぜた合金を設計する時、この原子半径の違いが特に重要。無理に大きい原子を詰め込むと格子が歪んで、強度が上がったり逆に割れやすくなったりする。シミュレーターで「BCC」のまま半径を大きくしてみると、中心の原子が隣とぶつかりそうになるのがわかるよね?

物理モデルと主要な数式

結晶構造の基本は、原子が空間的に繰り返し配列した「単位格子」です。格子定数 $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$ が変化します。

よくある質問

画面上の結晶モデルをマウスでドラッグすると、任意の角度に回転させて観察できます。また、自動回転アニメーションは初期状態で有効になっており、視点を変えたい場合はドラッグ操作で一時的に停止し、再度ドラッグを止めると自動回転が再開します。
配位数は、選択した結晶構造において最近接原子の数を自動でカウントします。充填率は、単位格子内の原子総体積を単位格子の体積で割った値(n·(4/3)πr³ / a³)を、各構造の理論値に基づいて計算・表示します。原子半径rは格子定数aから自動決定されます。
グラフェンは2次元構造ですが、本ツールではアイソメトリックビューで立体的に表示されます。厚み方向は極薄く表現されますが、回転操作により蜂の巣格子の平面構造を斜めからも観察でき、結合の様子や六員環の配列を直感的に理解できます。
g(r)は、ある原子から距離rの位置に他の原子が存在する確率密度を示します。ピークが立つ位置が最近接・第二近接原子間距離に対応し、結晶の周期性や原子配列の規則性を定量的に評価できます。例えばFCCでは第一ピークがa/√2、第二ピークがaに現れます。

実世界での応用

金属材料設計:自動車のボディに使われる鋼板は、主にBCC構造の鉄(フェライト)です。強度と加工性のバランスを取るために、炭素などの元素を添加し、結晶粒の大きさ(シミュレーターの「表示範囲」に相当)を制御します。CAEでは、この微細構造をモデル化して衝突安全性をシミュレーションします。

半導体デバイス:スマートフォンの頭脳であるCPUは、シリコンのダイヤモンド構造が基礎です。この規則正しい格子に、わずかに不純物(ドーパント)を添加することで、電気を通しやすくしたり遮断したりする性質(半導体特性)を生み出しています。

電池材料開発:リチウムイオン電池の正極材料(例:コバルト酸リチウム)は、層状の結晶構造(HCPに近い)を持ちます。リチウムイオンがその層の間を出入りすることで充放电が行われます。充填率やイオンの通り道(格子間隙)の設計が、電池容量や寿命を決めます。

軽量高強度材料:航空機の機体に使われるアルミニウム合金は、FCC構造をベースとしています。FCCはすべり面が多く加工性が良いため、複雑な形状に成形できます。ここに銅やマグネシウムなどを添加して析出強化させ、軽量でありながら高い強度を実現しています。

よくある誤解と注意点

まず、「原子は硬いボール」という単純化モデルをそのまま信じすぎないでください。NovaSolverで表示される球体は、原子の「電子雲」の広がりを便宜的に半径で表したもの。実際の化学結合では、電子が共有されたり、軌道が混ざり合ったりするので、単純な幾何学的接触とは異なります。例えば、充填率が低いダイヤモンド構造が極めて硬いのは、この「硬いボール」モデルだけでは説明できず、強い共有結合の方向性が鍵です。

次に、パラメータを独立に変えられると誤解しがちです。ツールでは「最近接距離」スライダーを動かすと原子半径が変わり、格子定数も連動して変わります。しかし実材料では、格子定数は元素の種類でほぼ決まり、自由に変えられるものではありません。例えば、純鉄のBCC構造(フェライト)の格子定数は約0.286nm。ここに直径の大きいモリブデン原子を添加すると、無理やり格子を「引き伸ばす」ことで材料を強化します。ツールでBCCの半径を大きくして歪みを作る操作は、まさにこの固溶強化のイメージを理解するのに役立ちます。

最後に、「表示範囲」を単位格子だけに見ていると全体像を見失う点に注意。実在する材料は、この単位格子が数億〜数兆個も集まった「多結晶体」で、それぞれの結晶粒の向き(方位)がバラバラです。CAEで材料の異方性(方向による強度の違い)を評価する時は、この集合体としての振る舞いをシミュレーションします。ツールで表示範囲を広げて格子を繰り返し表示するのは、この「単結晶」から「多結晶」へのイメージを掴む第一歩です。

使い方ガイド

  1. 結晶構造型(立方晶・HCP・ダイヤモンド)をs-aで選択
  2. 格子定数aをv-aで入力(例:Al単体ではa=4.05Å、Fe-BCC型ではa=2.87Å)
  3. s-cellsNumで超格子セル数を指定(1~5範囲で計算負荷調整)し、3D構造を生成
  4. v-rotで回転角度を設定し、[100][110][111]など主要結晶方位を確認
  5. 配位数・充填率・径分布関数(RDF)が自動算出される

具体的な計算例

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Åの整数倍関係を確認できます。

実務での注意点