ナノクラスタの幾何学と物性を支配するマジック数学的関係式

概要

触媒、量子ドット光ルミネッセンス、プラズモン共鳴などのサイズ依存的表面現象は、金属および半導体ナノクラスタの配位化学によって支配されます。本研究では、構造シェル数 n と主要なクラスタ属性を関連付ける簡潔な「マジックフォーミュラ」を用いて、これらの依存性を定量化します。面心立方格子（fcc）、体心立方格子（bcc）、単純立方格子（sc）、六方最密充填構造（hcp）、ダイヤモンド型立方格子、および正多面体クラスタの19種類の異なる幾何学を解析しました。純粋な空間座標から隣接行列と距離行列を構築し、ウィーナー指数、超ウィーナー指数、逆ウィーナー指数、シェゲド指数などのトポロジカル指数を計算し、それらの閉形式が n のマジック関数として現れることを示しました。本結果により、sc格子が頂点次数情報量の観点から最も単純であり、触媒設計に関連する表面原子分散の定量的尺度が提供されることが明らかになりました。

はじめに

ナノクラスタにおけるマジック数の概念は、van HardeveldとHartogの1969年の研究[1]に遡り、ナノサイエンス時代より前から存在します。以来、マジック数列は二次元多角形、三次元多面体、フラーレン、幅広い元素クラスタで観察され、1000件以上のWeb of Science記録がこの現象を引用しています[2-6]。クラスタのサイズと形状がその物理的・化学的特性（光学的、触媒的、電子的、磁気的）をどう決定するかを理解することは、現代ナノテクノロジーにとって不可欠です[7]。

マジック数は、クラスタが完全な原子シェルを完成したときに生じ、その結果としての配置は、最近接配位数、総原子数、表面対バルク比を支配する一連の特徴的な整数を示します。グラフ理論の観点では、原子は頂点、結合は辺となり、玉ねぎの層のようなネストされたシェルが明らかになります。本研究は、fcc、bcc、hcp格子に関する先行研究[1, 2]を、sc、ダイヤモンド型立方、正多面体の導入、およびトポロジカル指数と表面原子分散の探求によって拡張したものです。

表面原子は小さなクラスタの特性を支配します。露出表面原子分率FEとその幾何学依存性は、触媒活性、プラズモン応答、量子ドット光ルミネッセンスに影響を与えます[8-9]。本研究ではFEを定量化し、単位格子寸法から導かれる相対サイズパラメータd_relで正規化することで、クラスタ幾何学を比較しています。

トポロジカル指数は長らく分子構造と物性の相関に用いられてきました。ウィーナー指数から始まり、超ウィーナー指数、逆ウィーナー指数、シェゲド指数など、接続パターンを捉えるグラフベースの指標群が発展してきました[10-12]。本稿では、検討したすべてのクラスタタイプについてこれらの指数の厳密なマジックフォーミュラを導出し、n に関する多項式関係が根底にあることを明らかにしました。

手法

各クラスタ幾何学について、原子座標をアルゴリズム的に生成し、隣接行列 A と距離行列 D を構築します。2原子 i と j のユークリッド距離 r_ij が r_ij < r_c を満たす場合に結合しているとみなし、ここで r_c = 1.32·r_min、r_min は最短結合長です。隣接行列の要素は以下で定義されます：

$$ \mathbf{A}(i,j)=\begin{cases} 1 &\text{if } r_{ij}< r_c \text{ and } i\neq j\\ 0 &\text{otherwise} \end{cases} $$

頂点 i の配位数は A の第 i 列の和に等しい。距離行列 D は全頂点対間の最短経路長を記録します：

$$ \mathbf{D}(i,j)=\begin{cases} 0 & i=j\\ d_{ij} & i\neq j \end{cases} $$

A から D を計算する効率的なアルゴリズムが存在します[13]。これらの行列を用いて、確立されたプロトコル[14]に従い、ウィーナー指数 W(G)、超ウィーナー指数 WW(G)、逆ウィーナー指数 rW(G)、シェゲド指数 Sz(G) を評価します。

クラスタサイズはシェル数 n（シェル番号は0から n）で定量化します。総原子数は

$$ N_T(n)=\sum_{c_n=1}^{c_{nM}} N_{c_n}(n) $$

ここで N_{c_n}(n) は配位数 c_n を持つ原子数、c_{nM} はクラスタ内の最大配位数です。表面原子は c_n < c_{ns} を満たすため、

$$ N_S(n)=\sum_{c_n=1}^{c_{ns}} N_{c_n}(n) $$

総結合数は二重計数から導かれます：

$$ N_B(n)=\frac12\sum_{c_n=1}^{c_{nM}} c_n\,N_{c_n}(n) $$

これらの式と隣接行列データを組み合わせることで、正多面体について先に導出された多項式フォーミュラ[2]との厳密な照合が可能になります。

各幾何学について、n の最初の4つの整数値（交互の振る舞いを示す構造については適切なパリティ部分集合）を計算し、データを n の三次多項式にフィッティングします。シンボリック代数により厳密な有理係数が保証されます。N_T(n)、N_B(n)、配位分布のマジックフォーミュラは補足表にまとめられています。

露出表面原子分率は以下で定義されます：

$$ \text{FE}=\frac{N_S}{N_T}\times 100\% $$

相対クラスタサイズは以下で表されます：

$$ d_{rel}=b\,N_T^{1/3},\qquad b=d_{at}^{-1}\left(\frac{6V_u}{\pi n_u}\right)^{1/3} $$

ここで d_at は共有結合原子直径、V_u は単位格子体積、n_u は単位格子あたりの原子数です。定数 b は fcc/hcp で1.105、bcc で1.137、sc で1.488、ダイヤモンド型立方で1.517となります。このスケーリングにより、結晶構造に依存しないクラスタ間比較が可能になります。

結果と考察

表20（ここでは未掲載）は2018年までのナノクラスタ研究の主要なマイルストーンをまとめ、遷移金属クラスタ合成技術の急速な拡大を示しています。

FCCクラスタ

貴金属や主要な触媒元素を含む8種の遷移金属がfcc格子を採用しています。fccクラスタの実験的合成は立方体、八面体、立方八面体、二十面体に及び、表面エネルギーが低い(111)面が好まれます[7]。本研究のfcc菱面十二面体のマジックフォーミュラは明確な偶数/奇数列を示し、シェル番号を合わせると文献と一致します（表1参照）。導出された式は[2, 4]の報告と一致し、明示的な n 依存形式を提供することでそれらを拡張しています。

BCCクラスタ

7種の遷移金属がbcc格子で結晶化し、鉄はこの構造を持つ唯一の磁性元素です。鉄ナノキューブが合成されており、fcc鉄クラスタも報告されています[26]。5つの標準的なbcc形状のマジックフォーミュラを提示します（表2）。bcc構造はfcc対応物よりも複雑な配位パターンを示します。

HCPクラスタ

12種の遷移金属がhcp構造を持ちますが、酸化しやすさと化学的関心の薄さから、離散クラスタとして合成された例は少ないです。六方双錐形をとる金クラスタが観察されています[27]が、切頭六方双錐形はα-Fe₂O₃でのみ見られます[28]。hcp幾何学のマジックフォーミュラを提供します（表3）。

正多面体クラスタ

正多面体（立方体、正四面体、正八面体、正二十面体、正十二面体）は原型的なクラスタ幾何学として機能します。二十面体、十二面体、正四面体、体心正四面体クラスタのマジックフォーミュラを提示します（表4-7）。特筆すべきは、十二面体が表面配位数7というユニークな値を示し、連続するシェルで混合的な表面/バルク配位が生じることです。金ナノクラスタはすべての正多面体形状で合成されており[29, 30]、理論予測を裏付けています。

ダイヤモンド型立方、単純立方、十面体クラスタ

シリコン、ゲルマニウム、ダイヤモンド炭素はダイヤモンド型立方格子を採用し、水素終端シリコンクラスタが8-15 nmの立方体形状を達成しています[31, 32]。このようなSi-Hクラスタの組成を Si_{8n³+6n²−9n+5}H_12n−8 と予測します（n はシェル数）。scクラスタはポロニウムでのみ実現されますが、従来の一次元ハイパーキューブ研究[4]を発展させ、完全なマジックフォーミュラを詳細化しました（表8）。

マジックトポロジカルフォーミュラ

ToposPro（v5.3.2.2）を用い、結晶学オープンデータベースのcifファイルから頂点次数エントロピーを計算しました：

$$ I_{vd}=\sum_{i=1}^{v} a_i \log_2 a_i $$

表22によれば、sc格子は複雑性ゼロであるのに対し、fcc、hcp、bcc、ダイヤモンド型立方格子は I_vd が増大します。これらの指標は安定なトポロジカル指数の導出難易度と相関します。例えば、bcc立方体のシェゲド指数は高複雑性のため解けませんでした。

ウィーナー、超ウィーナー、逆ウィーナー、シェゲド指数のマジックフォーミュラを表23-25に掲載します。多項式次数は7-9に達し、複雑な格子の豊かな接続性を反映しています。特筆すべきは、scクラスタが最も単純な形式を与える一方、fcc構造が最も扱いやすい式を提供することです。

分散

図2（下記参照）は、露出表面原子分率FEを相対サイズパラメータd_relに対してプロットしたものです。2.2 nmの白金クラスタ（d_rel≈7.5）では、FE値は二十面体47.9%、立方八面体52.8%、十面体57.5%、八面体58.9%となります。八面体形状が与えられたサイズで最高の表面露出を示し、その優れた酸素還元反応（ORR）活性と相関します[38]。より大きな13 nmクラスタ（d_rel≈25-30）でも同様の傾向が見られますが、FE差は縮小します。しかし二十面体クラスタの歪み効果が触媒性能に影響を与える可能性があります[40, 41]。

結論

19種類のナノクラスタ幾何学をカタログ化し、原子数、結合数、配位分布、4つの主要トポロジカル指数について厳密なマジックフォーミュラを導出しました。これらの結果の多くは文献に新規であり、クラスタ構造と物性予測のギャップを埋めます。特に定量化された表面原子分散は、触媒挙動の予測に有用な指標となります。これらの包括的データセットが、特性を調整したナノクラスタ設計を行う研究者にとって貴重なリソースとなることを期待します。

略語

bcc:

体心立方格子

cif:

結晶学情報ファイル

fcc:

面心立方格子

FE:

露出表面原子分率

hcp:

六方最密充填構造

ORR:

酸素還元反応

rW(G):

逆ウィーナー指数

Sz(G):

シェゲド指数

W(G):

ウィーナー指数

WW(G):

超ウィーナー指数