コンテンツにスキップ

化合物の結晶成長の第一原理シミュレーション：LaB6結晶の平衡形状#

PHASE logo

電子顕微鏡などの高性能な電子銃カソード（陰極材料）として広く用いられる六ホウ化ランタン（LaB6）は、低い仕事関数と優れた熱電子放出特性を持ちます。ナノ粒子や単結晶の特性は「最表面にどの結晶面が露出しているか」に強く依存するため、合成・結晶成長の環境を制御し、望ましい形状（ファセット）を持つ結晶を作り出すことが重要です。本解析では、第一原理計算ソフトウェアAdvance/PHASEを用いて精緻に算出したLaB6(001)面の表面エネルギーと、文献から得られた他の結晶面のデータを統合し、さらにPythonライブラリPymatgenに実装されているウルフの定理に基づく形状予測機能を組み合わせることで、材料を合成する際の化学環境（La-rich環境およびB-rich環境）の違いが、LaB6結晶の平衡形状にどのような変化をもたらすのかをシミュレーションしました。

Keywords: 第一原理計算, DFTシミュレーション, 表面エネルギー, 結晶成長, 平衡形状, ウルフの定理, LaB6

理論的背景と計算手順#

ナノ粒子や微結晶の熱力学的な安定形状は、全体の表面自由エネルギーを最小化するような多面体として記述されます。これを幾何学的に決定するのが「ウルフの定理（Wulff construction）」です。

表面エネルギーと化学ポテンシャル#

多成分系化合物における特定の結晶面の表面エネルギー は、周囲の環境（各構成元素の化学ポテンシャル ）に依存して変化します。一般に、表面エネルギーは以下の式で定義されます。

ここで、 は表面積、 はスラブモデルの全エネルギー、 は各元素の原子数です。化合物の組成比（化学量論）から外れた表面（非化学量論的表面）では、特定の元素が豊富な環境（Rich限界）か乏しい環境（Poor限界）かによって、最も安定な表面構造や表面エネルギーが変動します。

計算の手順#

  1. バルク計算: 結晶の安定な格子定数と各元素の化学ポテンシャルの基準値を決定します。
  2. 表面スラブ計算: 主要な結晶面に対して、十分な厚みを持つスラブモデルを作成し、原子位置の最適化を伴うDFT計算を行います。
  3. 高精度化処理: 数値的なノイズを排除するため、層数依存性を考慮した線形フィッティング法などを用い、収束した表面エネルギー値を抽出します [1]。
  4. 平衡形状の構築: 得られた各面の表面エネルギーをWulffの定理に基づき統合し、3D形状として可視化します。

LaB6(001)面における計算結果と文献との比較#

LaB6の表面特性については、Uijttewaalらによる先行研究 [2] で詳細な議論がなされています。論文内の図9（(001)面の詳細な相図）および図13（各面の最安定面のまとめ図）では、(001)面における「La終端」と「B終端」の安定性が極限環境で逆転する微細な挙動が示唆されていますが、この複雑な相転移の境界付近の挙動については、さらに詳細な理論的検証を行う余地があります。

本解析（および関連する基礎解析 [3]）において、Advance/PHASEを用いて実施した独自の高精度計算結果を、文献 [2] の各データ点と比較したものを表2に示します。その結果、Advance/PHASEによるバルク形成エネルギーや各終端面の表面エネルギーは、文献の図9および図13にプロットされた値と良好に一致することが確認されました。

表2. LaB6(001)面およびバルク物性の比較（Advance/PHASE vs. 文献値）

比較項目 文献値（Uijttewaalら [2]） Advance/PHASE 計算値 [3]
バルク形成エネルギー () -3.96 eV -3.92 eV
(001)La終端 La-rich 表面エネルギー 0.80 J/m2 (図13) 0.814 J/m2
(001)La終端 B-rich 表面エネルギー 2.634 J/m2 (図13) 2.639 J/m2
(001)B終端 La-rich 表面エネルギー 圧倒的に高い (図9) 圧倒的に高い
(001)B終端 B-rich 表面エネルギー 2.40 J/m2 (図9) 2.438 J/m2

※文献[2]の図から読み取った値（eV/nm2）を J/m2 に換算して引用

LaB6表面エネルギーの全体像#

今回の検証結果に基づき、LaB6(001)面についてはAdvance/PHASEによる計算値を、その他の主要な結晶面については文献 [2] の値を採用して、La-rich限界およびB-rich限界における表面エネルギーを整理しました（図1）。

LaB6 Surface Energies Histogram

図1. LaB6の主要な結晶面の表面エネルギー比較。(001)面（*印）はAdvance/PHASEによる計算値、その他の面は文献 [2] より採用。青はLa-rich環境、オレンジ色はB-rich環境を示します。

La-rich環境では(001)面が突出して安定ですが、B-rich環境に移行すると表面終端の切り替わりと結晶面間の相対安定性の変化により、(111)面が最安定面となります。

LaB6結晶成長時の平衡形状#

図1の表面エネルギーの異方性データを用い、Pymatgen [4] によって構築されたLaB6結晶の平衡形状（Wulff形状）を図2に示します。合成環境の違いが、結晶の外形に大きな変化をもたらすことがわかります。

Wulff shapes of LaB6 in La-rich and B-rich limits

図2. LaB6結晶の熱力学的平衡形状。左：La-rich環境下（(001)面が支配的）、右：B-rich環境下（(111)面が支配的）。

  • La-rich環境（図2左）： (001)面の表面エネルギーが圧倒的に低いため、結晶は(001)面に包まれ、立方体（Cube）に近い形となります。これは実験的に得られるLaB6単結晶の典型的な形状 [5] を再現しています。
  • B-rich環境（図2右）： 最安定面が(111)面にシフトするため、結晶形状は八面体（Octahedron）に近い多面体へと変貌します。角の部分には(100)面や(112)面がわずかに露出する複雑なファセットが形成されます。

まとめ#

本事例では、第一原理計算ソフトウェアAdvance/PHASEを用い、文献データと統合することで、LaB6結晶の表面安定性と形状予測を行いました。文献データの詳細な解析と自社計算結果の比較検証を通じ、(001)面における精緻な物理現象を正確に捉えることができました。表面エネルギーの異方性を評価することで、合成条件のコントロールによって特定の形状・ファセットを持つナノ粒子を設計するための確かな理論的指針が得られます。

本解析の詳細や、研究への適用可能性に関するご相談はこちら

お問い合わせ

参考文献#

  1. V. Fiorentini and M. Methfessel, "Extracting convergent surface energies from slab calculations", J. Phys.: Condens. Matter 8, 6525 (1996).
  2. M. A. Uijttewaal, G. A. de Wijs, and R. A. de Groot, "Ab initio and work function and surface energy anisotropy of LaB6", J. Phys. Chem. B 110, 18459 (2006).
  3. 電子銃陰極材料LaB6の表面エネルギー計算：表面終端と化学ポテンシャル
  4. S. P. Ong, W. D. Richards, A. Jain, G. Hautier, M. Kocher, S. Cholia, D. Gunter, V. Chevrier, K. A. Persson, and G. Ceder, "Python Materials Genomics (pymatgen) : A Robust, Open-Source Python Library for Materials Analysis", Comp. Mater. Sci. 68, 314 (2013).
  5. H. Zhang, J. Tang, Q. Zhang, G. Zhao, G. Yang, J. Zhang, O. Zhou, and L.-C Qin, "Field Emission of Electrons from Single LaB6 Nanowires", Adv. Mater. 18, 87 (2006).

関連ページ#