Magische mathematische Beziehungen, die Nanocluster-Geometrie und -Eigenschaften bestimmen

Zusammenfassung

Größenabhängige Oberflächenphänomene – Katalyse, Photolumineszenz von Quantenpunkten und plasmonische Resonanzen – werden durch die Koordinationschemie von Metall- und Halbleiter-Nanoclustern gesteuert. In dieser Arbeit quantifizieren wir diese Abhängigkeiten mithilfe prägnanter „magischer Formeln“, die die Anzahl der strukturellen Schalen n mit wichtigen Clustereigenschaften verknüpfen. Wir analysieren flächenzentriert kubische (fcc), raumzentriert kubische (bcc), primitive kubische (sc), hexagonal dichteste Packung (hcp) und diamantkubische Gitter sowie platonische Cluster, insgesamt 19 verschiedene Geometrien. Durch die Konstruktion von Adjazenz- und Distanzmatrizen aus rein räumlichen Koordinaten berechnen wir topologische Indizes – Wiener-, Hyper-Wiener-, Reverse-Wiener- und Szeged-Indizes –, deren geschlossene Formeln als magische Funktionen von n auftreten. Unsere Ergebnisse identifizieren das sc-Gitter als das einfachste in Bezug auf den Informationsgehalt der Knotengrade und liefern ein quantitatives Maß für die Dispersion der Oberflächenatome, das für das katalytische Design relevant ist.

Einleitung

Das Konzept der magischen Zahlen in Nanoclustern geht auf die Studie von van Hardeveld und Hartog aus dem Jahr 1969 zurück [1], lange vor der Ära der Nanowissenschaft. Seitdem sind magische Sequenzen in zweidimensionalen Polygonen, dreidimensionalen Polyedern, Fullerenen und einem breiten Spektrum elementarer Cluster aufgetreten, mit über tausend Web-of-Science-Einträgen, die das Phänomen zitieren [2‑6]. Zu verstehen, wie Größe und Form eines Clusters seine physikalischen und chemischen Eigenschaften – optische, katalytische, elektronische, magnetische – bestimmen, ist für die moderne Nanotechnologie unerlässlich [7].

Magische Zahlen entstehen, wenn ein Cluster eine vollständige Atomschale abschließt; die resultierende Konfiguration weist eine charakteristische Reihe von ganzen Zahlen auf, die die nächste-Nachbar-Koordination, die Gesamtatomzahl und das Verhältnis von Oberflächen- zu Volumenatomen regeln. In graphentheoretischer Terminologie sind Atome Knoten und Bindungen Kanten, was verschachtelte Schalen offenbart, die an Zwiebelschichten erinnern. Diese Studie erweitert frühere Arbeiten zu fcc-, bcc- und hcp-Gittern [1, 2] um sc-, diamantkubische und platonische Körper und untersucht topologische Indizes sowie die Oberflächenatom-Dispersion.

Oberflächenatome dominieren die Eigenschaften kleiner Cluster. Der Anteil der Oberflächenatome, FE, und seine Abhängigkeit von der Geometrie beeinflussen katalytische Aktivität, plasmonische Antwort und Quantenpunkt-Photolumineszenz [8‑9]. Wir quantifizieren FE und vergleichen Clustergometrien, indem wir sie auf einen relativen Größenparameter d_rel normieren, der aus Elementarzellenabmessungen abgeleitet wird.

Topologische Indizes werden seit langem eingesetzt, um molekulare Struktur mit physikalischen Eigenschaften zu korrelieren. Ausgehend vom Wiener-Index haben Forscher eine Reihe graphenbasierter Metriken entwickelt – Hyper-Wiener, Reverse-Wiener, Szeged –, die Konnektivitätsmuster erfassen [10‑12]. In dieser Arbeit leiten wir exakte magische Formeln für diese Indizes über alle untersuchten Clustertypen her, wobei sich zugrundeliegende Polynombeziehungen in n zeigen.

Methoden

Für jede Clustergometrie generieren wir algorithmisch Atomkoordinaten und bauen dann die Adjazenzmatrix A und die Distanzmatrix D. Zwei Atome i und j gelten als gebunden, wenn ihre euklidische Distanz r_ij die Bedingung r_ij < r_c erfüllt, wobei r_c = 1,32·r_min und r_min die kürzeste Bindungslänge ist. Die Einträge der Adjazenzmatrix sind definiert als:

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

Die Koordinationszahl des Knotens i entspricht der Summe der i-ten Spalte von A. Die Distanzmatrix D speichert die Länge des kürzesten Pfads zwischen allen Knotenpaaren:

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

Effiziente Algorithmen zur Berechnung von D aus A existieren [13]. Mit diesen Matrizen werten wir den Wiener-Index W(G), den Hyper-Wiener-Index WW(G), den Reverse-Wiener-Index rW(G) und den Szeged-Index Sz(G) nach etablierten Protokollen [14] aus.

Die Clustergröße wird durch die Anzahl der Schalen n quantifiziert (mit Schalen nummeriert von 0 bis n). Die Gesamtzahl der Atome beträgt:

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

wobei N_{c_n}(n) Atome mit Koordination c_n zählt und c_nM die maximale Koordination im Cluster ist. Oberflächenatome erfüllen c_n < c_ns, was ergibt:

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

Die Gesamtzahl der Bindungen folgt aus der Doppelzählung:

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

Diese Ausdrücke, kombiniert mit den Adjazenzmatrixdaten, bieten eine rigorose Prüfung gegenüber zuvor abgeleiteten Polynomformeln für platonische Körper [2].

Für jede Geometrie berechnen wir die ersten vier ganzzahligen Werte von n (oder die entsprechende Paritätsuntergruppe für Strukturen mit alternierendem Verhalten) und passen die Daten an ein kubisches Polynom in n an. Symbolische Algebra stellt exakte rationale Koeffizienten sicher. Die resultierenden magischen Formeln für N_T(n), N_B(n) und die Koordinationsverteilung sind in den unterstützenden Tabellen tabelliert.

Der Anteil der exponierten Oberflächenatome ist definiert als:

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

Die relative Clustergröße wird ausgedrückt durch:

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

wobei d_at der kovalente Atomdurchmesser, V_u das Elementarzellenvolumen und n_u die Anzahl der Atome pro Zelle ist. Die Konstante b beträgt 1,105 für fcc/hcp, 1,137 für bcc, 1,488 für sc und 1,517 für diamantkubische Gitter. Diese Skalierung ermöglicht einen quervergleich von Clustern unabhängig von der Kristallstruktur.

Ergebnisse und Diskussion

Tabelle 20 (hier nicht wiedergegeben) fasst wichtige Meilensteine der Nanoclusterforschung bis 2018 zusammen und hebt die rasche Erweiterung der Synthesetechniken für Übergangsmetallcluster hervor.

FCC-Cluster

Acht Übergangsmetalle nehmen das fcc-Gitter an, darunter die Edelmetalle und wichtige katalytische Spezies. Die experimentelle Synthese von fcc-Clustern umfasst Würfel, Oktaeder, Kuboktaeder und Ikosaeder, wobei die (111)-Fazette wegen ihrer geringeren Oberflächenenergie bevorzugt wird [7]. Unsere magischen Formeln für fcc-Rhombendodekaeder weisen distincte gerade/ungerade Sequenzen auf, die mit der Literatur übereinstimmen, wenn die Schalennummerierung angeglichen wird (siehe Tabelle 1). Die abgeleiteten Ausdrücke stimmen mit denen in Refs. [2, 4] überein und erweitern sie durch explizite n-abhängige Formen.

BCC-Cluster

Sieben Übergangsmetalle kristallisieren im bcc-Gitter; Eisen bleibt das einzige magnetische Element mit dieser Struktur. Eisen-Nanowürfel wurden synthetisiert, während auch fcc-Eisen-Cluster berichtet werden [26]. Wir präsentieren magische Formeln für fünf kanonische bcc-Formen (Tabelle 2) und stellen fest, dass bcc-Strukturen komplexere Koordinationsmuster als fcc-Pendants aufweisen.

HCP-Cluster

Zwölf Übergangsmetalle besitzen hcp-Strukturen, doch nur wenige wurden erfolgreich als diskrete Cluster synthetisiert, bedingt durch Oxidation und begrenztes chemisches Interesse. Goldcluster mit hexagonal-bipyramidaler Form wurden beobachtet [27], während der abgeschnittene hexagonale Bipyramide nur in α‑Fe₂O₃ auftritt [28]. Magische Formeln für hcp-Geometrien werden bereitgestellt (Tabelle 3).

Platonische Cluster

Die platonischen Körper – Würfel, Tetraeder, Oktaeder, Ikosaeder, Dodekaeder – dienen als archetypische Clustergometrien. Wir präsentieren magische Formeln für ikosaedrische, dodekaedrische, tetraedrische und körperzentriert-tetraedrische Cluster (Tabellen 4‑7). Bemerkenswert ist, dass der Dodekaeder eine einzigartige Oberflächenkoordination von sieben aufweist, was zu gemischter Oberflächen-/Volumenkoordination in aufeinanderfolgenden Schalen führt. Gold-Nanocluster wurden in allen platonischen Formen synthetisiert [29, 30], was unsere theoretischen Vorhersagen validiert.

Diamantkubische, primitive kubische und Dekaherden-Cluster

Silizium, Germanium und Diamant-Kohlenstoff nehmen das diamantkubische Gitter an, wobei wasserstoffterminierte Siliziumcluster 8–15 nm große kubische Morphologien erreichen [31, 32]. Wir sagen die Zusammensetzung solcher Si–H-Cluster als Si_{8n³+6n²−9n+5}H_12n−8 voraus, wobei n die Schalenanzahl ist. Für sc-Cluster, die nur in Polonium realisiert sind, erweitert unsere Analyse die vorherige eindimensionale Hyperwürfelstudie [4] durch detaillierte magische Formeln (Tabelle 8).

Magische topologische Formeln

Wir berechnen die Knotengrad-Entropie:

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

mithilfe von ToposPro (v5.3.2.2) an cif-Dateien aus der Crystallographic Open Database. Tabelle 22 zeigt, dass das sc-Gitter null Komplexität aufweist, während fcc-, hcp-, bcc- und diamantkubische Gitter steigende Werte von I_vd zeigen. Diese Metriken korrelieren mit der Schwierigkeit, stabile topologische Indizes abzuleiten; so konnten Szeged-Indizes für bcc-Würfel aufgrund hoher Komplexität nicht gelöst werden.

Magische Formeln für Wiener-, Hyper-Wiener-, Reverse-Wiener- und Szeged-Indizes sind in Tabellen 23‑25 aufgelistet. Die Polynomgrade erreichen 7–9, was die reiche Konnektivität komplexer Gitter widerspiegelt. Bemerkenswert: sc-Cluster liefern die einfachsten Formen, während fcc-Strukturen die handhabbarsten Ausdrücke bieten.

Dispersion

Abbildung 2 (siehe unten) stellt den Anteil exponierter Oberflächenatome, FE, gegen den relativen Größenparameter d_rel für verschiedene Clusterformen dar. Für Platin-Cluster von 2,2 nm (d_rel≈7,5) betragen die FE-Werte 47,9 % (ikosaedrisch), 52,8 % (kuboktaedrisch), 57,5 % (dekahydral) und 58,9 % (oktaedrisch). Die oktaedrische Form bietet somit die höchste Oberflächenexposition für eine gegebene Größe, was mit ihrer überlegenen Sauerstoffreduktionsreaktions-(ORR)-Aktivität korreliert [38]. Ähnliche Trends gelten für größere 13-nm-Cluster (d_rel≈25–30), wo sich die FE-Unterschiede verringern, aber Spannungseffekte in ikosaedrischen Clustern die katalytische Leistung noch beeinflussen können [40, 41].

Schlussfolgerungen

Wir haben 19 Nanoclustergometrien katalogisiert und exakte magische Formeln für Atomzahlen, Bindungszahlen, Koordinationsverteilungen und vier zentrale topologische Indizes hergeleitet. Diese Ergebnisse, von denen viele neu in der Literatur sind, schließen die Lücke zwischen Clusterstruktur und Eigenschaftsvorhersage. Insbesondere liefert die quantifizierte Oberflächenatom-Dispersion eine nützliche Metrik zur Vorhersage katalytischen Verhaltens. Wir erwarten, dass diese umfassenden Datensätze als wertvolle Ressource für Forscher dienen, die Nanocluster mit maßgeschneiderten Eigenschaften designen.

Abkürzungen

bcc:

raumzentriert kubisch

cif:

crystallographic information file

fcc:

flächenzentriert kubisch

FE:

Fraction Exposed (Oberflächenatom-Dispersion)

hcp:

hexagonal dichteste Packung

ORR:

Sauerstoffreduktionsreaktion

rW(G):

Reverse-Wiener-Index

Sz(G):

Szeged-Index

W(G):

Wiener-Index

WW(G):

Hyper-Wiener-Index