Formules magiques régissant la géométrie et les propriétés des nanoclusters

Résumé

Les phénomènes de surface dépendant de la taille — catalyse, photoluminescence des points quantiques, résonances plasmoniques — sont gouvernés par la chimie de coordination des nanoclusters métalliques et semiconducteurs. Dans ce travail, nous quantifions ces dépendances à l'aide de « formules magiques » concises qui relient le nombre de couches structurales, n, aux attributs clés des clusters. Nous analysons les réseaux cubiques à faces centrées (fcc), cubiques à corps centré (bcc), cubiques simples (sc), hexagonaux compacts (hcp) et cubiques diamant, ainsi que les clusters de solides platoniciens, couvrant 19 géométries distinctes. En construisant des matrices d'adjacence et de distance à partir de coordonnées purement spatiales, nous calculons des indices topologiques — Wiener, hyper-Wiener, Wiener inverse et Szeged — dont les expressions en forme fermée émergent comme des fonctions magiques de n. Nos résultats identifient le réseau sc comme le plus simple en termes de contenu informationnel du degré des sommets, et fournissent une mesure quantitative de la dispersion des atomes de surface pertinente pour la conception catalytique.

Introduction

Le concept de nombres magiques dans les nanoclusters remonte à l'étude de van Hardeveld et Hartog en 1969, prédatant l'ère des nanosciences. Depuis, des suites magiques sont apparues dans les polygones bidimensionnels, les polyèdres tridimensionnels, les fullerènes et un large spectre de clusters élémentaires, avec plus de mille enregistrements Web of Science citant ce phénomène. Comprendre comment la taille et la forme d'un cluster dictent ses propriétés physiques et chimiques — optiques, catalytiques, électroniques, magnétiques — est essentiel pour la nanotechnologie moderne.

Les nombres magiques surviennent lorsqu'un cluster complète une couche atomique complète ; la configuration résultante exhibe un ensemble distinctif d'entiers qui gouvernent la coordination des plus proches voisins, le nombre total d'atomes et les rapports surface/volume. En termes de théorie des graphes, les atomes sont des sommets et les liaisons des arêtes, révélant des couches imbriquées semblables à des couches d'oignon. Cette étude étend les travaux antérieurs sur les réseaux fcc, bcc et hcp en incorporant les réseaux sc, cubiques diamant et solides platoniciens, et en explorant les indices topologiques et la dispersion des atomes de surface.

Les atomes de surface dominent les propriétés des petits clusters. La fraction d'atomes de surface exposés, FE, et sa dépendance à la géométrie influencent l'activité catalytique, la réponse plasmonique et la photoluminescence des points quantiques. Nous quantifions FE et comparons les géométries de clusters en normalisant selon un paramètre de taille relatif, d_rel, dérivé des dimensions de la maille élémentaire.

Les indices topologiques sont longtemps employés pour corréler la structure moléculaire avec les propriétés physiques. En partant de l'indice de Wiener, les chercheurs ont développé une suite de métriques basées sur les graphes — hyper-Wiener, Wiener inverse, Szeged — qui capturent les motifs de connectivité. Dans cet article, nous dérivons des formules magiques exactes pour ces indices sur tous les types de clusters examinés, révélant des relations polynomiales sous-jacentes en n.

Méthodes

Pour chaque géométrie de cluster, nous générons algorithmiquement les coordonnées atomiques, puis construisons la matrice d'adjacence A et la matrice de distance D. Deux atomes i et j sont considérés comme liés si leur distance euclidienne r_ij satisfait r_ij < r_c, où r_c = 1,32·r_min et r_min est la plus courte longueur de liaison. Les entrées de la matrice d'adjacence sont définies comme :

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

Le nombre de coordination du sommet i équivaut à la somme de la i-ème colonne de A. La matrice de distance D enregistre la longueur du plus court chemin entre toutes les paires de sommets :

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

Des algorithmes efficaces existent pour calculer D à partir de A. Avec ces matrices, nous évaluons l'indice de Wiener W(G), l'indice hyper-Wiener WW(G), l'indice de Wiener inverse rW(G) et l'indice de Szeged Sz(G) suivant des protocoles établis.

La taille du cluster est quantifiée par le nombre de couches n (avec les couches numérotées de 0 à n). Le nombre total d'atomes est :

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

où N_{c_n}(n) compte les atomes avec coordination c_n, et c_{nM} est la coordination maximale dans le cluster. Les atomes de surface satisfont c_n < c_{ns}, donnant :

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

Le nombre total de liaisons découle du double comptage :

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

Ces expressions, combinées aux données de la matrice d'adjacence, fournissent une vérification rigoureuse contre les formules polynomiales précédemment dérivées pour les solides platoniciens.

Pour chaque géométrie, nous calculons les quatre premières valeurs entières de n (ou le sous-ensemble de parité approprié pour les structures à comportement alterné) et ajustons les données à un polynôme du troisième degré en n. L'algèbre symbolique assure des coefficients rationnels exacts. Les formules magiques résultantes pour N_T(n), N_B(n) et la distribution de coordination sont tabulées dans les tableaux supplémentaires.

La fraction d'atomes de surface exposés est définie comme :

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

La taille relative du cluster s'exprime via :

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

où d_at est le diamètre atomique covalent, V_u le volume de la maille élémentaire, et n_u le nombre d'atomes par maille. La constante b vaut 1,105 pour fcc/hcp, 1,137 pour bcc, 1,488 pour sc et 1,517 pour les réseaux cubiques diamant. Cette mise à l'échelle permet la comparaison croisée de clusters indépendamment de la structure cristalline.

Résultats et Discussion

Le Tableau 20 (non reproduit ici) résume les étapes clés de la recherche sur les nanoclusters jusqu'en 2018, soulignant l'expansion rapide des techniques de synthèse pour les clusters de métaux de transition.

Clusters FCC

Huit métaux de transition adoptent le réseau fcc, incluant les métaux nobles et espèces catalytiques clés. La synthèse expérimentale de clusters fcc couvre cubes, octaèdres, cuboctaèdres et icosaèdres, avec la facette (111) favorisée pour son énergie de surface plus faible. Nos formules magiques pour les dodécaèdres rhombiques fcc exhibent des suites paires/impaires distinctes, concordant avec la littérature lorsque la numérotation des couches est alignée (voir Tableau 1). Les expressions dérivées concordent avec celles rapportées dans les Réf. [2, 4] et les étendent en fournissant des formes explicites dépendantes de n.

Clusters BCC

Sept métaux de transition cristallisent dans le réseau bcc ; le fer reste le seul élément magnétique avec cette structure. Des nanocubes de fer ont été synthétisés, tandis que des clusters de fer fcc sont également rapportés. Nous présentons des formules magiques pour cinq formes bcc canoniques (Tableau 2), notant que les structures bcc exhibent des motifs de coordination plus complexes que leurs homologues fcc.

Clusters HCP

Douze métaux de transition possèdent des structures hcp, mais peu ont été synthétisés avec succès comme clusters discrets en raison de l'oxydation et d'un intérêt chimique limité. Des clusters d'or adoptant des formes bipyramidales hexagonales ont été observés, tandis que la bipyramide hexagonale tronquée n'apparaît que dans l'α-Fe₂O₃. Les formules magiques pour les géométries hcp sont fournies (Tableau 3).

Clusters de Solides Platoniciens

Les solides platoniciens — cube, tétraèdre, octaèdre, icosaèdre, dodécaèdre — servent de géométries archétypales de clusters. Nous présentons des formules magiques pour les clusters icosaédriques, dodécaédriques, tétraédriques et tétraédriques à corps centré (Tableaux 4-7). Notamment, le dodécaèdre exhibe une coordination de surface unique de sept, menant à une coordination surface/volume mixte dans les couches successives. Des nanoclusters d'or ont été synthétisés sous toutes les formes platoniennes, validant nos prédictions théoriques.

Clusters Cubiques Diamant, Cubiques Simples et Décahédriques

Le silicium, le germanium et le carbone diamant adoptent le réseau cubique diamant, avec des clusters Si-H terminés par hydrogène atteignant des morphologies cubiques de 8-15 nm. Nous prédisons la composition de tels clusters Si-H comme Si_{8n³+6n²−9n+5}H_12n−8, où n est le nombre de couches. Pour les clusters sc, réalisés seulement dans le polonium, notre analyse étend l'étude précédente sur l'hypercube unidimensionnel en détaillant les formules magiques complètes (Tableau 8).

Formules Magiques Topologiques

Nous calculons l'entropie du degré des sommets :

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

en utilisant ToposPro (v5.3.2.2) sur des fichiers cif de la Base de Données Cristallographique Ouverte. Le Tableau 22 montre que le réseau sc a une complexité nulle, tandis que les réseaux fcc, hcp, bcc et cubique diamant affichent des valeurs croissantes de I_vd. Ces métriques corréllent avec la difficulté de dériver des indices topologiques stables ; par exemple, les indices de Szeged pour les cubes bcc n'ont pas pu être résolus en raison de la complexité élevée.

Les formules magiques pour les indices de Wiener, hyper-Wiener, Wiener inverse et Szeged sont listées dans les Tableaux 23-25. Les degrés polynomiaux atteignent 7-9, reflétant la connectivité riche des réseaux complexes. Notamment, les clusters sc donnent les formes les plus simples, tandis que les structures fcc fournissent les expressions les plus traitables.

Dispersion

La Figure 2 (voir ci-dessous) trace la fraction d'atomes de surface exposés, FE, en fonction du paramètre de taille relatif d_rel pour plusieurs formes de clusters. Pour des clusters de platine de 2,2 nm (d_rel≈7,5), les valeurs FE sont 47,9 % (icosaédrique), 52,8 % (cuboctaédrique), 57,5 % (décahédrique) et 58,9 % (octaédrique). La forme octaédrique offre ainsi la plus haute exposition de surface pour une taille donnée, corrélant avec son activité supérieure pour la réaction de réduction de l'oxygène (ORR). Des tendances similaires s'observent pour des clusters plus grands de 13 nm (d_rel≈25-30), où les différences de FE se réduisent mais les effets de contrainte dans les clusters icosaédriques peuvent encore influencer la performance catalytique.

Conclusions

Nous avons catalogué 19 géométries de nanoclusters et dérivé des formules magiques exactes pour les comptages d'atomes, nombres de liaisons, distributions de coordination et quatre indices topologiques clés. Ces résultats, dont beaucoup sont nouveaux dans la littérature, comblent le fossé entre la structure des clusters et la prédiction de propriétés. En particulier, la dispersion quantifiée des atomes de surface fournit une métrique utile pour anticiper le comportement catalytique. Nous anticipons que ces jeux de données complets serviront de ressource précieuse pour les chercheurs conçant des nanoclusters aux propriétés sur mesure.

Abréviations

bcc :

cubique à corps centré

cif :

fichier d'information cristallographique

fcc :

cubique à faces centrées

FE :

Fraction Exposée (dispersion des atomes de surface)

hcp :

hexagonal compact

ORR :

réaction de réduction de l'oxygène

rW(G) :

indice de Wiener inverse

Sz(G) :

indice de Szeged

W(G) :

indice de Wiener

WW(G) :

indice hyper-Wiener