Super Atomic Clusters: Design Rules and Potential for Building Blocks of Materials

Electronic structure of triangular M3 (M = B, Al, Ga): nonclassical three-center two electron π bond and σ delocalization

The small molecule clusters have received more and more attention due to their widespread applications in chemical insulators, explosives, semiconductors and the high energy density materials industry. The electron deficiency of group IIIA elements endows their clusters with interesting properties. In this work, the electronic structures of M₃ (M = B, Al, Ga) have been investigated by means of a complete active space self-consistent field (CASSCF) method. The nature of the chemical bond has been analyzed using the quantum theory of atoms in molecules (QTAIM) and electron localization function (ELF) analyses. The following conclusions can be drawn: in M₃ (M = B, Al, Ga) clusters, two π electrons are shared by three atoms forming a 3c-2e delocalization π bond. Going from B₃ to Al₃ to Ga₃, more and more electrons move from the bond pair to the outside of the M atom, which leads to a gradual enhancement of the delocalization of σ electrons. Aromaticity and the adaptive natural density partitioning (AdNDP) analyses reveal the existence of the 3c-2e π bond and delocalization of σ electrons in the studied systems.

Covalently Linked, Two-Dimensional Quantum Dot Assemblies

Using nanoscale building blocks to construct hierarchical materials is a radical new branch point in materials discovery that promises new structures and emergent functionality. Understanding the design principles that govern nanoparticle assembly is critical to moving this field forward. By exploiting mixed ligand environments to target patchy nanoparticle surfaces, we have demonstrated a novel method of colloidal quantum dot (QD) assembly that gives rise to 2D structures. The equilibration of solutions of spherical and quasispherical QDs, including CdS, CdSe, and InP, with 2,2'-bipyridine-5,5'-diacrylic acid resulted in the preferential formation of 2D assemblies over the course of days as determined by transmission electron microscopy analysis. Small-angle X-ray scattering confirms the existence of the QD assemblies in solution. The dependence of the assembly on linker properties (length and rigidity), linker concentration, and total concentration was investigated, together with the data point to a mechanism involving ligand redistribution to create a patchy surface that maximizes the steric repulsion of neighboring QDs. By operating in an underexchanged regime, the arising patchiness results in enthalpically preferred directions of cross-linking that can be accessed by thermal equilibration.

Cyclodextrin-Assisted Hierarchical Aggregation of Dawson-type Polyoxometalate in the Presence of {Re6Se8} Based Clusters

The association of metallic clusters (CLUS) and polyoxometalates (POM) into hierarchical architectures is achieved using γ-cyclodextrin (γ-CD) as a supramolecular connector. The new self-assembled systems, so-called CLUSPOM, are formed from Dawson-type polyoxometalate [P₂W₁₈O₆₂]^6- and electron-rich rhenium clusters. It is worth noting that a cluster-based cation [{Re₆Se₈}(H₂O)₆]²⁺ on one hand and a cluster-based anion on the other hand [{Re₆Se₈}(CN)₆]^4- can be associated with the anionic POM. In the absence of the supramolecular connector, a "CLUSPOM salt" was obtained from aqueous solution of the cationic cluster and the polyoxometalate. In this solid, the arrangement between the polymetallic building blocks is mainly governed by long-range Coulombic interactions. In the presence of γ-CD, the Dawson anion and the cationic cluster are assembled differently, forming a hierarchical supramolecular solid, K₂[{Re₆Se₈}(H₂O)₆]₂{[P₂W₁₈O₆₂]@2γ-CD}·42H₂O, where the organic macrocycle acts as a ditopic linker between the inorganic building blocks. In such an edifice, the short-range molecular recognition dominates the long-range Coulombic interactions leading to a specific three-dimensional organization. Interestingly, the assembling of anionic POM [P₂W₁₈O₆₂]^6- with the anionic rhenium cluster [{Re₆Se₈}(CN)₆]^4- is also achieved with γ-CD despite the repulsive forces between the nanosized anions. The resulting solid, K₁₀{[{Re₆Se₈}(CN)₆]@2γ-CD}[P₂W₁₈O₆₂]·33H₂O, is built from 1:2 inclusion complexes {[{Re₆Se₈}(CN)₆]@2γ-CD}^4- linked by a POM unit interacting with the exterior wall of the organic macrocycle. Multinuclear NMR and small-angle X-ray scattering investigations support supramolecular preorganization in aqueous solution prior to crystallization.

Modified Lennard-Jones potentials for nanoscale atoms

A classical 6-12 Lennard-Jones (LJ) equation has been widely used to model materials and is the potential of choice in studies when the focus is on fundamental issues. Here we report a systematic study comparing the pair interaction potentials within solid-state materials (i.e., [Co₆ Se₈ (PEt₃ )₆ ][C₆₀ ]₂ , [Cr₆ Te₈ (PEt₃ )₆ ][C₆₀ ]₂ , [Ni₉ Te₆ (PEt₃ )₈ ][C₆₀ ]) using density functional theory (DFT) calculations and LJ parametrization. Both classical (6-12 LJ) and modified LJ (mLJ) models were developed. In the mLJ approach, the exponents 6 and 12 are replaced by different integer number n and 2n, respectively, and an additional parameter (α) is introduced to describe intermolecular distance shift arising within the geometric centers' approach (instead of the shortest interatomic distance between particles). A general LJ approach reexamination reveals that in the case of nanoatoms, the attractive term decays with distance as the inverse fourth power, and the dominating at short distances repulsive term decays as the inverse eighth power. The modification of the LJ equation is even more prominent for interaction profiles, where intermolecular distance corresponds to separation between geometric centers of particles. In this approach, the attractive term decays with distance as the inverse 12th power, while the repulsive term decays rapidly (as the inverse 24th power). Thus, the mLJ models (e.g., 4-8 LJ) rather than the 6-12 classical ones seem to be a better choice for the description of binary interactions of nanoatoms. The developed mLJ models and electronic structure characteristics give an insight into the explanation of the unique physicochemical properties of superatomic-based solid-state materials.

Heavy chalcogenide-transition metal clusters as coordination polymer nodes

While metal-oxygen clusters are widely used as secondary building units in the construction of coordination polymers or metal-organic frameworks, multimetallic nodes with heavier chalcogenide atoms (S, Se, and Te) are comparatively untapped. The lower electronegativity of heavy chalcogenides means that transition metal clusters of these elements generally exhibit enhanced coupling, delocalization, and redox-flexibility. Leveraging these features in coordination polymers provides these materials with extraordinary properties in catalysis, conductivity, magnetism, and photoactivity. In this perspective, we summarize common transition metal heavy chalcogenide building blocks including polynuclear metal nodes with organothiolate/selenolate or anionic heavy chalcogenide atoms. Based on recent discoveries, we also outline potential challenges and opportunities for applications in this field.

Anionic Copper Clusters Reacting with NO: An Open-Shell Superatom Cu18. J Phys Chem Lett 2020;11:5807-5814. [PMID: 32597656 DOI: 10.1021/acs.jpclett.0c01643]

Gas-phase metal clusters have been a subject of research interest for allowing reliable strategies to explore the stability and reactivity of materials at reduced sizes with atomic precision. Here we have prepared well-resolved copper cluster anions Cu_n^- (n = 7-37) and systematically studied their reactivity with O₂, NO, and CO. We found remarkable stability of an open-shell cluster Cu₁₈^-, which is comparable with the closed-shell clusters Cu₁₇^- and Cu₁₉^- within the picture of an electronic shell model. Even without having a magic number of valence electrons, intriguingly, the unpaired electron on the singly occupied molecular orbital of Cu₁₈^- is mainly contributed by the central copper atom, while the other 18 delocalized valence electrons occupy the lower-energy superatomic orbitals of the cluster. The finding of such an open-shell superatom Cu₁₈^-, with an electron configuration of 1S²1P⁶1D¹⁰2S¹||1F⁰, is interesting in the sense that an elementary cluster of coinage metal atoms could still behave as a superatom mimicking coinage metals like silver or gold atoms with an empty f orbital. The superatomic stability of this Cu₁₈^- cluster is reinforced by the unique electrostatic interaction between the Cu^- core and Cu₁₇ shell, which provides new insights into the chemistry of metal clusters.

On the Interaction between Superatom Al12Be and DNA Nucleobases/Base Pairs: Bonding Nature and Potential Applications in O2 Activation and CO Oxidation

The interaction between quasi-chalcogen superatom Al₁₂Be and DNA nucleobases/base pairs has been explored by searching for the most stable Al₁₂Be-X (X = DNA bases and base pairs) complexes. Our results reveal that Al₁₂Be prefers to combine with guanine by two Al-O and Al-N bonds rather than the other DNA bases, no matter in free state or base pair. The formed Al-N and Al-O bonds between Al₁₂Be and DNA bases proved to be strong polar covalent bonds by the Wiberg bond index, nature bond orbitals, atoms in molecules theory, localized molecular orbitals, and electron localization functions analyses. More importantly, it is found that the formed global minimum of Al₁₂Be-G has the ability to activate an oxygen molecule into a peroxide dianion ¹O₂ ^2-, which can further catalyze the CO oxidation via the Eley-Rideal mechanism with a small energy barrier of 7.78 kcal/mol. We hope that this study could not only provide an in-depth understanding on the intermolecular interaction between metallic superatoms and DNA at the molecular level but also attract more interest in designing and synthesizing superatom-based heterogeneous catalysts with DNA/nucleobases as basic building blocks.

A super stable assembled P nanowire with variant structural and magnetic/electronic properties via transition metal adsorption

By means of first-principles calculations, we systematically investigated the structure, stability and magnetic and electronic properties of one-dimensional P nanowire (1D-P10 NW) assembled by Pn subunits (n = 2, 8) and transition metal doped 1D-P10 NW. Our calculations showed that the assembled 1D-P10 NW is super stable in thermodynamic, dynamic, thermal and chemical perspectives. Moreover, when the assembled 1D-P10 NW is decorated with transition metals (TM = Ti ∼ Zn, Zr ∼ Mo), structural transformation occurs (to sandwich or quasi-sandwich chains), and various magnetic and electronic characteristics are introduced to the nanowire. Particularly, the sandwich chains 1D-Mn2@P10 and 1D-V1@P5 are a ferromagnetic semiconductor and a ferromagnetic half-metal, respectively, and the magnetic anisotropy energies are both ∼0.3 meV per Mn/V atom. Our theoretical studies proposed a super stable 1D P nanowire and also offer a feasible approach to reach P5-TM-P5-TM chains with diverse magnetic and electronic properties, as well as ferromagnetic vdW-type 2D systems, which are promising in nanoelectronic devices and spintronics.

Probing the structural evolution and electronic properties of divalent metal Be2Mgn clusters from small to medium-size

Bimetallic clusters have aroused increased attention because of the ability to tune their own properties by changing size, shape, and doping. In present work, a structural search of the global minimum for divalent bimetal Be₂Mg_n (n = 1-20) clusters are performed by utilizing CALYPSO structural searching method with subsequent DFT optimization. We investigate the evolution of geometries, electronic properties, and nature of bonding from small to medium-sized clusters. It is found that the structural transition from hollow 3D structures to filled cage-like frameworks emerges at n = 10 for Be₂Mg_n clusters, which is obviously earlier than that of Mg_n clusters. The Be atoms prefer the surface sites in small cluster size, then one Be atom tend to embed itself inside the magnesium motif. At the number of Mg larger than eighteen, two Be atoms have been completely encapsulated by caged magnesium frameworks. In all Be₂Mg_n clusters, the partial charge transfer from Mg to Be takes place. An increase in the occupations of the Be-2p and Mg-3p orbitals reveals the increasing metallic behavior of Be₂Mg_n clusters. The analysis of stability shows that the cluster stability can be enhanced by Be atoms doping and the Be₂Mg₈ cluster possesses robust stability across the cluster size range of n = 1-20. There is s-p hybridization between the Be and Mg atoms leading to stronger Be-Mg bonds in Be₂Mg₈ cluster. This finding is supported by the multi-center bonds and Mayer bond order analysis.

Bonding properties of a superatom system with high-Z elements: insights from energy decomposition analysis

Superatoms with high-Z elements have novel physicochemical properties, and a comprehensive and thorough view of their bonding properties plays a crucial role in the design of superatoms. Now, energy decomposition analysis shows increasingly prominent performance for understanding inter- and intra-molecular interactions, so the bonding properties of typical superatoms with high-Z elements, Th@Au₁₄, have been investigated here. It is found that under different electron occupation types of the fragments, the electrostatic interaction energy, polarization, and exchange repulsion energy change significantly in their intramolecular interaction energy components, resulting in quantitative or even qualitative differences in their main interaction energy. This indicates that the bonding properties of fragments are related to their electronic structures, and even has extraordinary reference value for the future regulation and control of interactions in superatoms with high-Z elements, which has great significance for superatom development.

Superatoms with high-Z elements have novel physicochemical properties, and a comprehensive and thorough view of their bonding properties plays a crucial role in the design of superatoms.

Superatomic Gallium Clusters in Dendrimers: Unique Rigidity and Reactivity Depending on their Atomicity

Superatoms have been investigated due to their possible substitution for other elements. The solution-phase synthesis of superatoms has attracted attention to realize the availability of superatoms. However, the previous approach is basically limited to the formation of a single cluster. Here, superatoms are investigated and the number of valence electrons in these superatoms is changed by designing the number of gallium atoms present. Based on the dendrimer template method, clusters consisting of 3, 12, 13, and other numbers of atoms have been synthesized. The halogen-like superatomic nature of Ga₁₃ is structurally and electrochemically observed as completely different to the other clusters. The gallium clusters of 13 and 3 atoms, which can fill the 2P and 1P superatomic orbitals, respectively, exhibit different reactivities. The 3-atom gallium cluster is suggested as being reduced to Ga₃ H₂ ^- due to the lower shift of energy levels in the unoccupied orbitals. The results for these gallium clusters provide candidates for superatoms.

Iodization Threshold in Size-Dependent Reactions of Lead Clusters Pbn+ with Iodomethane

Utilizing a magnetron-sputtering (MagS) source in tandem with a multiple-ion laminar flow tube (MIFT) reactor and a customized triple quadrupole mass spectrometer (TQMS), we have prepared clean Pb_n⁺ (n = 1-13) clusters and measured their reactivity with iodomethane under high carrier gas pressures. Strong size dependences are found for the reactivity of these cationic Pb_n⁺ clusters with CH₃I. For the Pb_n⁺ with n ≤ 4, iodinated clusters Pb_nI⁺ were found to be the dominant products, in strong contrast to n > 4 where no such products were seen. Quantum chemical studies show that with an increasing number of Pb atoms, the Pb-Pb interatomic interactions become stronger compared with the Pb-I bonding in Pb_nI⁺ clusters. Furthermore, the reactions of Pb_1-4⁺ with CH₃I have fairly small transition state energy barriers, in contrast to those for Pb_n>4⁺ clusters which have magnitudes that will prevent reactions under the ambient conditions.

Discovery of a Superatom inside the Fullerene Cage

The stability of endohedral clusterfullerenes is generally understood in terms of the metal cluster shape, cage structure, and metal-cage interactions, with the electronic state of the internal cluster mostly neglected. Herein, theoretical calculations reveal that the (Ti₃C₃)⁶⁺ unit of recently synthesized Ti₃C₃@I_h(7)-C₈₀ exhibits a superatomic state with a perfect closed-shell 1S²1P⁶1D¹⁰ electronic configuration in accordance with the famous jellium model. This "trapped superatom" features considerable aromaticity and hyperconjugation interactions never reported for other clusterfullerenes. Besides the localized two-center two-electron (2c-2e) Ti-C/C-C bonds, it also has two 3c-2e Ti-C-Ti bonds. Furthermore, the ring strain of the cyclopropane-like C₃ core is effectively released upon the metal coordination. All these factors greatly stabilize the (Ti₃C₃)⁶⁺ cluster, showing the critical role of metal-to-cage charge transfer and cage encapsulation in enhancing the stability of this exotic metal cluster.

A machine learning based deep potential for seeking the low-lying candidates of Al clusters

A Machine-Learning based Deep Potential (DP) model for Al clusters is developed through training with an extended database including ab initio data of both bulk and several clusters in only 6 CPU/h. This DP model has good performance in accurately predicting the low-lying candidates of Al clusters in a broad size range. Based on our developed DP model, the low-lying structures of 101 different sized Al clusters are extensively searched, among which the lowest-energy candidates of 69 sized clusters are updated. Our calculations demonstrate that machine-learning is indeed powerful in generating potentials to describe the interaction of atoms in complex materials.

Iron-Cluster-Directed Synthesis of 2D/2D Fe-N-C/MXene Superlattice-like Heterostructure with Enhanced Oxygen Reduction Electrocatalysis

Metal clusters, an emerging category of materials with molecular metal dispersity, have been proposed for versatile applications. In this work, we show an unexpected function of metal clusters, which can contribute to preparing 2D/2D superlattice-like heterostructures. The key step is to use metal clusters to adjust the surface charge of 2D nanosheets and, consequently, match the charge negativities per surface area for different 2D nanosheets, which facilitate the electrical-driven assembly of these nanosheets into a superlattice-like heterostructure in aqueous solutions. Accordingly, iron-cluster-directed cationic Fe-N-C nanosheets (Zeta potential: +30.4 mV) have been assembled with anionic MXene (Zeta potential: -39.7 mV) to produce a superlattice-like heterostructure characteristic of a lateral size of around tens of nanometers, a surface area of 30 m² g^-1, and ultrathickness of several nanometers with repeated dimensions of 0.4 and 2.1 nm. Potential application of the synthesized Fe-N-C/MXene heterostructure has been demonstrated for electrocatalytic oxygen reduction reaction (ORR) that shows a positive onset potential of 0.92 V, four-electron transfer pathway, and strong durability of 20 h in alkaline electrolyte. This work suggests that metal clusters can assist the assembly of low-dimensional architectures for energy-related applications.

How metallic are noble-metal clusters? Static screening and polarizability in quantum-sized silver and gold nanoparticles

Metallicity of nanoparticles can be defined in different ways. One possibility is to look at the degree to which external fields are screened inside the object. This screening would be complete in a classical perfect metal where surface charges arrange on the classical -i.e., abrupt - surface such that no internal fields exist. However, it is obvious that this situation is modified for very small clusters: the surface charges are "smeared out" at the surface, and the screening might be less complete. In the present work we ask the question as to how close small noble-metal clusters are to a classical metal. We show that, indeed, the screening is almost complete (≈96%) already for as little as one atomic layer of the coinage metals, silver and gold alike. At the same time, we show that quantum effects, viz., electronic shell closings and the Friedel-like oscillations of the density, play a role, meaning that the clusters cannot be described solely using the concept of screening in a classical metal.

Tuning the Electronic Properties and Performance of Low-Temperature CO Oxidation of the Gold Cluster by Oriented External Electronic Field

Conventional electronic rules, including Jellium and Wade-Mingos rules and so on, have long been successfully dedicated to design superatoms. These rules, however, rely on altering the intrinsic properties, for example, the compositions or the number of valence electrons, of clusters, which is relatively complicated and inconvenient to manipulate, especially in experiments. Herein, by employing density functional theory calculations, the oriented external electric field (OEEF) was demonstrated to possess the capability of precisely and continuously regulating the electronic properties of clusters at will, representing a novel and noninvasive methodology in constructing stable superatoms because it hardly changes the geometries of clusters. More interestingly, the active sites formed by the charge redistribution upon the introduction of an OEEF could significantly promote the catalytic performance of the low-temperature CO oxidation over clusters. Considering the convenient source of the OEEF, the findings highlighted here may boost the potential applications of superatom-assembly nanomaterials in catalysis and materials science.

Identification of Magic Numbers in Homonuclear Clusters: The ε3 Stability Ranking Function

With the rise of cluster-assembled materials, an index that is able to rank and identify stable clusters or molecules is of great interest in materials sciences and engineering. In the present work, we applied a stability ranking function (ε³) in nanoclusters formed by simple metals (Na, Mg), main group elements (Al), or transition metals (Ti, Cu). The ε³ function parameters are molecular properties derived from the wave function. These parameters can be divided into kinetic and thermodynamic descriptors, in which the kinetic descriptors are the ionization potential and electronic excitation energy, while the atomization free Gibbs energy is the thermodynamic one. This simple ε³ function was able to identify the possible magic numbers of the studied clusters across the periodic table in a good agreement with previous experimental and theoretical works.

Ultra stable superatomic structure of doubly magic Ga13 and Ga13Li electrolyte

We report the extreme thermal stability of the superatomic electronic structure for 13-atom gallium clusters and the Ga₁₃Li electrolyte. Using previously-validated first-principles simulations, [K. G. Steenbergen and N. Gaston, Phys. Rev. B: Condens. Matter Mater. Phys., 2013, 88, 161402-161405] we show that the superatomic shell progression of doubly-magic Ga₁₃^- remains stable up to 1000 K, making this cluster an ideal candidate for high-temperature applications requiring an exceptionally stable electronic structure. Using the neutral and cationic clusters for comparison, we quantify the extent to which cluster stability (geometric and electronic) is modified through addition or subtraction of a single electron. Finally, combining 13-atom gallium with lithium, we illustrate that superatomic closed-shell Ga₁₃Li exhibits the same exceptionally high thermal stability as naked Ga₁₃^-. For technological use as a superatomic electrolyte, we demonstrate that Ga₁₃Li has a low affinity to water as well as a low Li⁺ binding energy.

Redefining the Mechanism of O2 Etching of Aln- Superatoms: An Early Barrier Controls Reactivity, Analogous to Surface Oxidation

New insights into aluminum anion cluster reactivity with O₂ were obtained through temperature-dependent kinetics measurements. Overall reactivity is controlled by a barrier at an avoided crossing where charge is transferred from the cluster to the O₂, mechanistically similar to what occurs as O₂ approaches a bulk Al surface. Contrary to prior interpretations, spin conservation does not inhibit the reaction of clusters with an odd number of Al atoms. In fact, the only spin constraint in these systems is on the reactivity of even clusters due to repulsive surfaces, not previously recognized. Although the superatom nature of Al₁₃^- is manifest in its high electron binding energy (EBE), the mechanism of its reactivity is not special; it reacts with O₂ as if it were a small piece of bulk Al. These experiments highlight the sensitivity of Al cluster reactivity with O₂ to temperature and EBE, uncovering routes to industrial scale use of aluminum superatom-based materials.

Singlet and triplet states of the sandwich-type Be2B6 and Be2B7+ clusters. A test for the electron counting rules of aromaticity

The studied complexes exhibit double aromaticity in their triplet states in line with the predictions of Hückel and Baird's rules.

A designer 32-electron superatomic CBe8H12 cluster: core–shell geometry, octacoordinate carbon, and cubic aromaticity

A 32-electron CBe₈H₁₂ cluster is designed with cubic octacoordinate carbon. It features core–shell geometry, two-fold superatomic bonding, and cubic aromaticity.

Boron-based ternary Rb6Be2B6 cluster featuring unique sandwich geometry and a naked hexagonal boron ring

Boron-based ternary Rb₆Be₂B₆ cluster features a naked hexagonal boron ring and unique “Big Mac” sandwich shape, being stabilized collectively by four-fold 2σ/6π/6σ/2σ aromaticity.

Footprints of Nanoscale DNA-Silver Cluster Chromophores via Activated-Electron Photodetachment Mass Spectrometry

DNA-templated silver clusters (AgC) are fluorescent probes and biosensors whose electronic spectra can be tuned by their DNA hosts. However, the underlying rules that relate DNA sequence and structure to DNA-AgC fluorescence and photophysics are largely empirical. Here, we employ 193 nm activated electron photodetachment (a-EPD) mass spectrometry as a hybrid MS³ approach to gain structural insight into these nanoscale chromophores. Two DNA-AgC systems are investigated with a 20 nt single-stranded DNA (ssDNA) and a 28 nt hybrid hairpin/single-stranded DNA (hpDNA). Both oligonucleotides template Ag₁₀ clusters, but the two complexes are distinct chromophores: the former has a violet absorption at 400 nm with no observable emission, while the latter has a blue-green absorption at 490 nm with strong green emission at 550 nm. Via identification of both apo and holo (AgC-containing) sequence ions generated upon a-EPD and mapping areas of sequence dropout, specific DNA regions that encapsulate the AgC are assigned and attributed to the coordination with the DNA nucleobases. These a-EPD footprints are distinct for the two complexes. The ssDNA contacts the cluster via four nucleobases (CCTT) in the central region of the strand, whereas the hpDNA coordinates the cluster via 13 nucleobases (TTCCCGCCTTTTG) in the double-stranded region of the hairpin. This difference is consistent with prior X-ray scattering spectra and suggests that the clusters can adapt to different DNA hosts. More importantly, the a-EPD footprints directly identify the nucleobases that are in direct contact with the AgC. As these contacting nucleobases can tune the electronic structures of the Ag core and protect the AgC from collisional quenching in solution, understanding the DNA-silver contacts within these complexes will facilitate future biosensor designs.

Transformation between Hexagonal Prism and Antiprism of the Singly and Doubly Cr-Doped Ge12 Clusters

Structural transformation is a unique characteristic of atomic clusters, but it turns out very different from cluster to cluster. This theoretical study proves that the isomeric transformation between hexagonal prism and hexagonal antiprism is found for the doubly doped Cr₂Ge₁₂ cluster but not for singly doped CrGe₁₂ cluster. We confirm that the ground state of CrGe₁₂ is the distorted hexagonal prism C_2h at the ³B_g triplet state instead of various shapes predicted in the previous studies. Upon comparison between the estimation at the B3P86/6-311+G(d) level of theory and the detachment energies measured by photoelectron spectroscopy, hexagonal antiprismatic shape is identified as the most stable isomer of the Cr₂Ge₁₂ cluster and it is easy to transform to the hexagonal prism-a less stable isomer by the rotation of the hexagonal rings. That is the first evidence for the structural transformation between a hexagonal prism and an antiprism of the germanium clusters, referring to the ability of Ge-based clusters in the formation of tubular geometry by doping Cr atoms. All the low-energy isomers of both Cr-doped germanium clusters have high magnetic moments. Interestingly, there is a tuning in magnetic properties of Cr₂Ge₁₂ from the ferromagnetism of the lowest-lying hexagonal antiprism to the ferrimagnetism of the higher-energy hexagonal prism. The stronger Cr-Cr bond and stronger interaction between the Cr₂ moiety and the antiprism cage are accounted for by the higher stability of the hexagonal antiprismatic isomer.

A C20-based 3D carbon allotrope with high thermal conductivity

Stimulated by the high thermal conductivity of diamond together with the light mass and rich resources of carbon, a great deal of effort has been devoted to the study of the thermal conductivity of carbon-based materials. In this work, we systematically study the thermal transport properties of a three dimensional (3D) C20 fullerene-assembled carbon allotrope, HSP3-C34, in which all carbon atoms are in sp3 hybridization. The stability of HSP3-C34 is confirmed and its thermal conductivity is obtained by using first principles calculations combined with solving the linearized phonon Boltzmann transport equation. At room temperature, the thermal conductivity of HSP3-C34 is 731 W m-1 K-1, which is larger than those of many 3D carbon allotropes, such as BCO-C16 (452 W m-1 K-1), 3D graphene (150 W m-1 K-1) and T-carbon (33 W m-1 K-1). A detailed analysis of its phonons reveals that three acoustic branches are the main heat carriers at room temperature, and the optical branches gradually become important with increasing temperature. A further study on the harmonic and anharmonic properties of HSP3-C34 uncovers that the main reasons for the high thermal conductivity are the weak anharmonicity and large group velocity resulting from the strong sp3 bonding. This study provides new insights on searching for carbon allotropes with high thermal conductivity.

Combining proton and silaborane-based superhalogen anions - an effective route to new superacids as verified via systematic DFT calculations

Based on systematic DFT calculations, silaborane-based superhalogen anions, which obey the Wade-Mingos rule, are shown to be capable of giving rise to superacids via their combination with protons. Compared to previous carborane-based systems, the acidities of the composites here are stronger in both the gas phase and solution phase. Thus, the potential of candidates based on silaborane could be greater than those based on carborane in the search for ultra-strong acidic systems. Within a given group, a higher superhalogen anion vertical electron detachment energy (VDE) generally leads to stronger acidity. This consistency arises from the dominant role of the VDE, as established through the decomposition of the gas-phase acidity into different contributions. Thus, constructing superacids from superhalogens is a rational route whose future should be positive. Besides the VDE, other effects, i.e., the deformation energy (DE) and bond dissociation energy (BDE), could also be crucial, especially in terms of the differences between the acidities of composites belonging to different groups. A comparison between the results in the gas phase and solution phase indicates that complete calculations of both gas-phase ΔG_acid and solution-phase pK_a values are necessary to obtain an unbiased description of the acidity. The solvation free energies of the participants in the deprotonation process, especially the conjugate acid, are responsible for the discrepancies between gas phase and solution phase behavior.

Polymeric tungsten carbide nanoclusters: structural evolution, ligand modulation, and assembled nanomaterials

Seeking novel superatoms with tunable electronic and magnetic properties has attracted much interest due to their potential application in cluster assembly nanomaterials. By employing density functional theory (DFT) calculations, the recently observed superatomic WC cluster was adopted as the basic unit to construct larger polymeric clusters, namely (WC)_n (n = 2-7), and their structural evolution was explored to understand the growth pattern of these superatomic clusters into nanoscale materials. An unusual odd-even pattern in structural evolution was disclosed, in which the (WC)₂ unit is considered as the basic building block. Moreover, W₄C₄ is found to possess a cubic structure, based on which the CO and PH₃ ligands were attached to examine their ligation effects on W₄C₄. Theoretical results show that the electronic properties of W₄C₄ can be dramatically altered during the ligation process. Intriguingly, the continuous attachment of CO and PH₃ ligands strongly increases and decreases the electron affinities (EA) and ionization potentials (IP) of the ligated W₄C₄ clusters, respectively, leading to the formation of superhalogen and superalkali species with high magnetic moments. The observed ligand induced strategy highlighted here could serve as an effective way to tune the electronic and magnetic properties of clusters resulting in the formation of novel superatoms. Finally, studies on the geometrical and electronic structures of the W₄C₄ cluster solid unveil its special 3-D cubic honeycomb geometry and metallic properties with predominant contribution from the 5d of W, which may have potential applications in electro-catalysis.

Greatly Enhanced Electron Affinities of Au2n Cl Clusters (n = 1-4): Effects of Chlorine Doping

Au2n Cl- (n = 1-4) clusters are investigated by both laser ablation mass spectrometry and theoretical calculations. It is interesting to find that the electron affinities of neutral Au2n Cl (n = 1-4) clusters are much larger than those of corresponding pure Au2n clusters. Among them, the electron affinity of Au2Cl is 4.02 eV, which can be defined as a very unique superhalogen that is quite different from classical ones of M n X m (M = metal, X = halogen, and n < m). Natural bond orbital and highest occupied molecular orbital analyses indicate that the extra electron is predominantly delocalized over the positively charged metal moiety in these anionic Au2n Cl- clusters, which is the main reason for the large electron affinities of the corresponding neutral species.

Size, Composition, and Support-Doping Effects on Oxygen Reduction Activity of Platinum-Alloy and on Non-platinum Metal-Decorated-Graphene Nanocatalysts

Recent investigations reported in the open literature concerning the functionalization of graphene as a support material for transition metal nanoparticle catalysts have examined isolated systems for their potential Oxygen Reduction Reaction (ORR) activity. In this work we present results which characterize the ability to use functionalized graphene (via dopants B, N) to upshift and downshift the adsorption energy of mono-atomic oxygen, O* (the ORR activity descriptor on ORR Volcano Plots), for various compositions of 4-atom, 7-atom, and 19-atom sub-nanometer binary alloy/intermetallic transition metal nanoparticle catalysts on graphene (TMNP-MDG). Our results show several important and interesting features: (1) that the combination of geometric and electronic effects makes development of simple linear mixing rules for size/composition difficult; (2) that the transition from 4- to 7- to 19-atom TMNP on MDG has pronounced effects on ORR activity for all compositions; (3) that the use of B and N as dopants to modulate the graphene-TMNP electronic structure interaction can cause shifts in the oxygen adsorption energy of 0.5 eV or more; (4) that it might be possible to make specific doped-graphene-Ni_xCu_y TMNP systems which fall close to the Volcano Peak for ORR. Our results point to systems which should be investigated experimentally and may improve the viability of future fuel cell or other ORR applications, and provide new paths for future investigations of more detail for TMNP-MDG screening.

Ternary 12-electron CBe3X3+ (X = H, Li, Na, Cu, Ag) clusters: planar tetracoordinate carbons and superalkali cations

Molecules with planar tetracoordinate carbons (ptCs) are exotic in chemical bonding, and they are normally designed according to the 18-electron rule. Here we report on the viability of ptC clusters with as few as 12 valence electrons, which represent the lower limit in terms of electron counting. Specifically, we have computationally designed a class of ternary 12-electron ptC clusters, CBe3X3+ (X = H, Li, Na, Cu, Ag), based on a rhombic CBe32- unit. Computer structural searches reveal that the ptC species are global minima, whose C center is coordinated in-plane by three Be atoms and a terminal X atom via robust C-Be/C-X bonding, either covalent or ionic. The other two X atoms are on the periphery and each bridge two Be atoms. Bonding analyses show that the ptC core is governed by delocalized 2π/6σ bonding, that is, double π/σ aromaticity, which collectively conforms to the 8-electron counting. Additional 4 electrons contribute to peripheral Be-X-Be and Be-Be σ bonding. The delocalized 2π/6σ frameworks appear to be universal for all ptC clusters, ranging from 18-electron down to 12-electron systems. In other words, the ptC species are dictated entirely by the 8-electron counting. Predicted vertical electron affinities of these ptC clusters range from 3.13 to 5.48 eV, indicative of superalkali or pseudoalkali cations.

Electronic shells of a tubular Au26 cluster: a cage-cage superatomic molecule based on spherical aromaticity

Gold clusters, which display a variety of unusual geometric structures due to their strong relativistic effects, have attracted much attention. Among them, Au₂₆ has a high-symmetry tubular structure (D_6d) with a large HOMO-LUMO energy gap, but its electronic stability still remains unclear. In this paper, the electronic nature of the Au₂₆ cluster is investigated using the density functional theory method. Depending on the super valence bond model, the tubular Au₂₆ cluster with 26 valence electrons could be viewed as a superatomic molecule composed of two open cages based on spherical aromaticity, and its molecule-like electronic shell closure is achieved via a super triple bond (σ, 2π) between the two cages. Based on this new cage-cage superatomic structural model, a series of similar tubular clusters are predicted from the Au₂₆ skeleton. The two capped Au atoms are replaced by Cu, Ag and In atoms, respectively, to form tubular D_6d Au₂₄Cu₂ and Au₂₄Ag₂ (26e) and Au₂₄In₂ (30e) clusters, where the super triple bonds also exist. Moreover, tubular D_5d Au₂₀In₂ (26e) is obtained by replacing hexatomic Au₆ rings in the bulk of Au₂₄In₂ with pentagonal Au₅ rings. Chemical bonding analysis reveals that there is a super quintuple bond (σ, 2π, 2δ) between two open (Au₁₀In) cages, in accordance with the 26e Li₂₀Mg₃ superatomic molecule composed of two icosahedral superatoms. Our study proposes the new cage-cage structural model of superatomic molecules based on spherical aromaticity, which extends the range of the super valence bonding pattern and gives inferences for further study of superatomic clusters.

Data Driven Determination in Growth of Silver from Clusters to Nanoparticles and Bulk

The thresholds among atomic clusters, nanoparticles, and the bulk state have been ambiguous. A potential solution is to determine cluster growth toward bulk, but this is challenging to determine with experiments and computation. Data science is proposed to predict atomic cluster growth and determine the cluster-nanoparticle-bulk thresholds using Ag clusters as a prototype element. Supervised machine learning reveals that Ag cluster growth has nonlinear models where nonlinear machine learning is found to accurately predict binding energy. Unsupervised machine learning discovers three groups (cluster, semiclusters, and nanoparticles) where linear regression is used to predict the binding energy in each group. Furthermore, machine learning reveals the linear relationship between binding energy and the surface-to-volume ratio of Ag nanoparticles. This allows for a binding energy estimation of large Ag nanoparticles and a revelation of how Ag nanoparticles grow toward the bulk. Thus, data science is proposed as a powerful tool for determining cluster growth and thresholds for clusters, nanoparticles, and bulk states.

Boron-based inorganic heterocyclic clusters: electronic structure, chemical bonding, aromaticity, and analogy to hydrocarbons

This Perspective article deals with recent computational and experimental findings in boron-based heterocyclic clusters, which focuses on binary B-O and B-S clusters, as well as relevant ternary B-X-H (X = O, S, N) species. Boron is electron-deficient and boron clusters do not form monocyclic rings or linear chains. Boron-based heterocyclic clusters are intuitively even more electron-deficient and feature exotic chemical bonding, which make use of O 2p, S 3p, or N 2p lone-pairs for π delocalization over heterocyclic rings, facilitating new cluster structures and new types of bonding. Rhombic, pentagonal, hexagonal, and polycyclic clusters are discussed herein. Rhombic species are stabilized by four-center four-electron (4c-4e) π bonding, that is, the o-bond. An o-bond cluster differs from a typical 4π antiaromatic system, because it has 4π electrons in an unusual bonding/nonbonding combination, which takes advantage of the empty 2p_z atomic orbitals from electron-deficient boron centers. A variety of examples (notably including boronyl boroxine) possess a hexagonal ring, as well as magic 6π electron-counting, making them new members of the inorganic benzene family. Pentagonal clusters bridge rhombic o-bond systems and inorganic benzenes, but they do not necessarily favor 6π electron-counting as in cyclopentadienide anion. In contrast, pentagonal 4π clusters are stable, leading to the concept of pentagonal o-bond. One electron can overturn the potential energy landscape of a system, enabling rhombic-to-hexagonal structural transition, which further reinforces the idea that 4π electron-counting is favorable for rhombic systems and 6π is magic for hexagonal rings. The bonding analogy between heterocyclic clusters and hydrocarbons goes beyond monocyclic species, which allows rational design of boron-based inorganic analogs of polycyclic aromatic hydrocarbons, including s-indacene as a puzzling aromatic/antiaromatic system. Selected linear B-O clusters are also briefly discussed, featuring dual 3c-4e π bonds, that is, ω-hyperbonds. Dual ω-hyperbonds, rhombic or pentagonal o-bond, and inorganic benzenes share a common chemical origin. The field of boron-based heterocyclic clusters is still in its infant stage, and much new chemistry remains to be discovered in forthcoming experimental and theoretical studies.