支配纳米团簇几何结构与性质的神奇数学关系

摘要

尺寸依赖的表面现象——催化、量子点光致发光和等离子体共振——由金属和半导体纳米团簇的配位化学决定。本研究利用简洁的“神奇公式”量化这些依赖关系，将结构壳层数 n 与关键团簇属性相关联。我们分析了面心立方（fcc）、体心立方（bcc）、简单立方（sc）、六方密堆积（hcp）和金刚石立方晶格，以及柏拉图固体团簇，共涵盖 19 种不同几何结构。通过从纯空间坐标构建邻接矩阵和距离矩阵，我们计算了拓扑指数——维纳指数、超维纳指数、逆维纳指数和塞格德指数——其闭式表达式作为 n 的神奇函数浮现。结果表明 sc 晶格在顶点度信息含量方面最简单，并提供了与催化剂设计相关的表面原子分散度定量度量。

引言

纳米团簇中神奇数的概念可追溯至 van Hardeveld 和 Hartog 1969 年的研究 [1]，早于纳米科学时代。此后，神奇序列出现在二维多边形、三维多面体、富勒烯及广泛的元素团簇中，Web of Science 收录的引用文献超过一千篇 [2-6]。理解团簇尺寸和形状如何决定其物理化学性质——光学、催化、电子、磁学——对现代纳米技术至关重要 [7]。

当团簇完成一个完整原子壳层时出现神奇数；由此产生的构型表现出一组独特的整数，支配最近邻配位数、总原子数和表面/体积比。用图论术语表述，原子为顶点，键为边，呈现出类似洋葱层的嵌套壳层。本研究在先前关于 fcc、bcc 和 hcp 晶格的工作 [1, 2] 基础上，引入了 sc、金刚石立方和柏拉图固体，并探讨了拓扑指数和表面原子分散度。

表面原子主导小团簇的性质。暴露表面原子分数 FE 及其对几何结构的依赖影响催化活性、等离子体响应和量子点光致发光 [8-9]。我们量化了 FE，并通过归一化到从晶胞尺寸导出的相对尺寸参数 d_rel 来比较团簇几何结构。

拓扑指数长期用于关联分子结构与物理性质。从维纳指数开始，研究者开发了一套基于图的指标——超维纳、逆维纳、塞格德——捕捉连接模式 [10-12]。本文推导了所有研究团簇类型这些指标的精确神奇公式，揭示了 n 的多项式潜在关系。

方法

对每种团簇几何结构，我们算法生成原子坐标，构建邻接矩阵 A 和距离矩阵 D。若两原子 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 的前四个整数值（或交替行为结构的适当奇偶子集），拟合 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 团簇

八种过渡金属采用 fcc 晶格，含贵金属及关键催化物种。实验合成的 fcc 团簇包括立方体、八面体、立方八面体和二十面体，其中 (111) 晶面因表面能较低而受青睐 [7]。我们针对 fcc 菱形十二面体推导的神奇公式呈现明显的奇/偶序列，对齐壳层编号后与文献一致（见表 1）。导出表达式与文献 [2, 4] 吻合，并通过提供显式 n 依赖形式加以扩展。

BCC 团簇

七种过渡金属结晶为 bcc 晶格；铁是唯一具有此结构的磁性元素。已合成铁纳米立方体，也有 fcc 铁团簇报道 [26]。我们给出五种典型 bcc 形状的神奇公式（表 2），指出 bcc 结构比 fcc 对应物展现更复杂的配位模式。

HCP 团簇

十二种过渡金属具有 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 团簇仅见于钋，我们的分析通过详述完整神奇公式（表 8），扩展了先前一维超立方体研究 [4]。

神奇拓扑公式

我们利用 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 种纳米团簇几何结构，推导了原子数、键数、配位分布和四个关键拓扑指数的精确神奇公式。这些结果（许多为文献首见）连接了团簇结构与性质预测。特别是量化的表面原子分散度为预测催化行为提供了有用指标。我们预期这些综合数据集将成为此类研究人员设计具有定制性质纳米团簇的宝贵资源。

缩写

bcc:

体心立方

cif:

晶体学信息文件

fcc:

面心立方

FE:

暴露分数（表面原子分散度）

hcp:

六方密堆积

ORR:

氧还原反应

rW(G):

逆维纳指数

Sz(G):

塞格德指数

W(G):

维纳指数

WW(G):

超维纳指数