Electronic and atomic structure, and magnetism of transition-metal clusters

Structural Evolution Study of Titanium-Vanadium-Niobium Nanoparticles from Single to Multicomponent Systems

In this article, the revised basin-hopping with mirror-rotation sampling combined with density functional theory, which was proposed by our previous study, was used to study the structural property of Ti n (n = 3m, m = 1-7), V n (n = 3m, m = 1-7), Nb n (n = 3m, m = 1-7), and Ti n V n Nb n (n = 1-7) systems. We found that equiatomic Ti n V n Nb n (n = 1-7) systems do not change their lowest energy structures relative to the same size Ti n (n = 3m, m = 1-7), V n (n = 3m, m = 1-7), and Nb n (n = 3m, m = 1-7) systems, and this indicates that the nanoparticles composed of titanium, vanadium, or niobium elements may have similar energy morphologies when the atomic number is the same. Based on the low-energy structural similarity of titanium-vanadium-niobium systems between single and multicomponent, we used the element space position replacement (ESPR) method to reconstruct the low-energy structure of Ti n V n Nb n (n = 1-7) systems. For the Ti7V7Nb7 system, the average sampling step of 10 separate searches of the BH-MRS method is 1226 more than that of the ESPR method to find the lowest energy structure (six-ring layered structure). The electronic property calculation shows that using equiatomic vanadium and niobium elements to replace titanium element in the Ti n (n = 3m, m = 1-7) system may not change its stability, and the Ti n (n = 3m, m = 1-7) system has better electron trapping ability than V n (n = 3m, m = 1-7), Nb n (n = 3m, m = 1-7), and Ti n V n Nb n (n = 1-7) systems. Our method and results can be helpful for the design of nanostructures of transition metals with better catalytic properties.

Exploring nonlinear correlations among transition metal nanocluster properties using deep learning: a comparative analysis with LOO-CV method and cosine similarity

A novel approach is introduced for the rapid and accurate correlation analysis of nonlinear properties in Transition Metal (TM) clusters utilizing the Deep Leave-One-Out Cross-Validation technique. This investigation demonstrates that the Deep Neural Network (DNN)-based approach offers a more efficient predictive method for various properties of fourth-row TM nanoclusters compared to conventional Density Functional Theory methods, which are computationally intensive and time-consuming. The feature space, also known as descriptors, is established based on a broad spectrum of electronic and physical characteristics. Leveraging the similarities among these clusters, the DNN-based model is employed to explore the correlations among TM cluster properties. The proposed method, in conjunction with cosine similarity, achieves remarkable accuracy up to 10-9 for predicting total energy, lowest vibrational mode, binding energy, and HOMO-LUMO energy gap of TM2, TM3, and TM4nanoclusters. By analyzing correlation errors, the most closely coupled TM clusters are identified. Notably, Mn and Ni clusters exhibit the highest and lowest levels of energy coupling with other TMs, respectively. Generally, energy prediction for TM2, TM3, and TM4clusters exhibit similar trends, while an alternating behavior is observed for vibrational modes and binding energies. Furthermore, Ti, V, and Co demonstrate the highest binding energy correlations with TM2, TM3, and TM4sets, respectively. Regarding energy gap predictions, Ni exhibits the strongest correlation in the smallest TM2clusters, while Cr shows the highest dependence in TM3and TM4sets. Lastly, Zn displays the largest error in HOMO-LUMO energy gap across all sets, indicating distinctive independent energy gap characteristics.

Investigating Molecular Adsorption on Graphene-Supported Platinum Subnanoclusters: Insights from DFT + D3 Calculations

Platinum (Pt) subnanoclusters have become pivotal in nanocatalysis, yet their molecular adsorption mechanisms, particularly on supported versus unsupported systems, remain poorly understood. Our study employs detailed density functional theory (DFT) calculations with D3 corrections to investigate molecular adsorption on Pt subnanoclusters, focusing on CO, NO, N2, and O2 species. Gas-phase and graphene-supported scenarios are systematically characterized to elucidate adsorption mechanisms and catalytic potential. Gas-phase Pt n clusters are first analyzed to identify stable configurations and assess size-dependent reactivity. Transitioning to graphene-supported Pt n clusters, both periodic and nonperiodic models are employed to study interactions with graphene substrates. Strong adsorbate interactions predominantly occur at single top sites, revealing distinct adsorption geometries and stabilization effects for specific molecules on Pt6 clusters. Energy decomposition analysis highlights the paramount role of graphene substrates in enhancing stability and modulating cluster-adsorbate interactions. The interaction energy emerges as a critical criterion within the Sabatier principle, crucial for distinguishing between physisorption and chemisorption. Hybridization indices and charge density flow tendencies establish direct relationships with stabilization processes, underscoring graphene's influence in stabilizing highly reactive subnanoclusters. This comprehensive investigation provides critical insights into molecular adsorption mechanisms and the catalytic potential of graphene-supported Pt nanoclusters. Our findings contribute to a deeper understanding of nanocatalysis, emphasizing the essential role of substrates in optimizing catalytic performance for industrial applications.

Hierarchical structures and magnetism of Co clusters: a perspective from integration of deep learning and a hybrid differential evolution algorithm

Theoretically determining the lowest-energy structure of a cluster has been a persistent challenge due to the inherent difficulty in accurate description of its potential energy surface (PES) and the exponentially increasing number of local minima on the PES with the cluster size. In this work, density-functional theory (DFT) calculations of Co clusters were performed to construct a dataset for training deep neural networks to deduce a deep potential (DP) model with near-DFT accuracy while significantly reducing computational consumption comparable to classic empirical potentials. Leveraging the DP model, a high-efficiency hybrid differential evolution (HDE) algorithm was employed to search for the lowest-energy structures of CoN (N = 11-50) clusters. Our results revealed 38 of these clusters superior to those recorded in the Cambridge Cluster Database and identified diverse architectures of the clusters, evolving from layered structures for N = 11-27 to Marks decahedron-like structures for N = 28-42 and to icosahedron-like structures for N = 43-50. Subsequent analyses of the atomic arrangement, structural similarity, and growth pattern further verified their hierarchical structures. Meanwhile, several highly stable clusters, i.e., Co13, Co19, Co22, Co39, and Co43, were discovered by the energetic analyses. Furthermore, the magnetic stability of the clusters was verified, and a competition between the coordination number and bond length in affecting the magnetic moment was observed. Our study provides high-accuracy and high-efficiency prediction of the optimal structures of clusters and sheds light on the growth trend of Co clusters containing tens of atoms, contributing to advancing the global optimization algorithms for effective determination of cluster structures.

Exploring the Structural and Electronic Properties of Niobium Carbide Clusters: A Density Functional Theory Study

This paper systematically investigates the structure, stability, and electronic properties of niobium carbide clusters, NbmCn (m = 5, 6; n = 1-7), using density functional theory. Nb5C2 and Nb5C6 possess higher dissociation energies and second-order difference energies, indicating that they have higher thermodynamic stability. Moreover, ab initio molecular dynamics (AIMD) simulations are used to demonstrate the thermal stability of these structures. The analysis of the density of states indicates that the molecular orbitals of NbmCn (m = 5, 6; n = 1-7) are primarily contributed by niobium atoms, with carbon atoms having a smaller contribution. The composition of the frontier molecular orbitals reveals that niobium atoms contribute approximately 73.1% to 99.8% to NbmCn clusters, while carbon atoms contribute about 0.2% to 26.9%.

DNA Framework-Guided Self-Limiting Aggregation for Highly Luminescent Metal Cluster Nanoaggregates

The photoluminescent properties of atomically precise metal nanoclusters (MCs) have garnered significant attention in the fields of chemical sensing and biological imaging. However, the limited brightness of single-component nanoclusters hinders their practical applications, and the conventional ligand engineering approaches have proven insufficient in enhancing the emission efficiency of MCs. Here, we present a DNA framework-guided strategy to prepare highly luminescent metal cluster nanoaggregates. Our approach involves an amphiphilic DNA framework comprising a hydrophobic alkyl core and a rigid DNA framework shell, serving as a nucleation site and providing well-defined nanoconfinements for the self-limiting aggregation of MCs. Through this method, we successfully produced homogeneous MC nanoaggregates (10.1 ± 1.2 nm) with remarkable nanoscale precision. Notably, this strategy proves adaptable to various MCs, leading to a substantial enhancement in emission and quantum yield, up to 3011- and 87-fold, respectively. Furthermore, our investigation using total internal reflection fluorescence microscopy at the single-particle level uncovered a more uniform photon number distribution and higher photostability for MC nanoaggregates compared to template-free counterparts. This DNA-templating strategy introduces a conceptually innovative approach for studying the photoluminescent properties of aggregates with nanoscale precision and holds promise for constructing highly luminescent MC nanoparticles for diverse applications.

Cage doping of Ti, Zr, and Hf-based 13-atom nanoclusters: two sides of the same coin

Transition metal nanoclusters can exhibit unique and tunable properties which result not only from their chemical composition but also from their atomic packing and quantized electronic structures. Here, we introduce a promising family of bimetallic TM@Ti12, TM@Zr12, and TM@Hf12 nanoclusters with icosahedral geometry, where TM represents an atom from groups 3 to 12. Density functional theory calculations show that their stability can be explained with familiar concepts of metal cluster electronic and atomic shell structures. The magnetic properties of these quasispherical clusters are entirely consistent with superatom electronic shells and Hund's rules, and can be tuned by the choice of the TM dopant. The computed cluster atomization energies were analyzed in terms of the elements' cohesive energy, Ecoh, and contributions from geometric distortion, Edis, surface energy, Es, and ionic bonding, Ei. Some clusters have anomalous stability relative to Ecoh + Edis + Es + Ei. We attribute this to superatomic character associated with a favorable atomic and electronic shell structure. This raises the possibility of designing stable superatoms and materials with tailored electronic and magnetic properties.

Magnetism and electronic properties of ConMoP (n = 1 ~ 5) cluster: a DFT study

CONTEXT

Hydrogen has emerged as a promising clean energy carrier, underscoring the imperative need to comprehend its adsorption mechanisms. This study delves into the magnetic and electronic properties of Co-Mo-P clusters, aiming to unveil their catalytic potential in hydrogen production. Employing density functional theory (DFT), we optimized cluster configurations and scrutinized their magnetic behaviors. Our investigation unveiled 16 stable configurations of the ConMoP (n = 1 ~ 5) cluster, predominantly in steric forms. The magnetic attributes were primarily ascribed to the d orbitals of Co metal atoms, with Co3MoP exhibiting exceptional magnetic characteristics. Analysis of density of state diagrams revealed the prevalence of spin-up α-electrons in d orbitals, while spin-down β-electrons attenuated overall magnetic properties. Localized orbital (LOL) analysis highlighted stable covalent bonds within the clusters, affirming their catalytic potential. Orbital delocalization index (ODI) analysis revealed diverse spatial distribution ranges for orbitals across different configurations, suggesting a progressive attenuation of off-domain properties with increasing cluster size. Furthermore, infrared spectroscopy unveiled distinct vibrational peaks in various configurations, indicative of unique infrared activities. These findings contribute to a nuanced theoretical understanding of Co-Mo-P clusters and pave the path for future research aimed at augmenting their catalytic efficiency in hydrogen production. This study underscores the viability of Co-Mo-P clusters as alternatives to conventional Pt catalysts, offering insights into the design of novel materials for sustainable energy applications. Further research is warranted to explore the behavior of the Co-Mo-P system under diverse reaction conditions, fostering advancements in materials and energy science.

METHODS

In this study, we harnessed the ConMoP (n = 1 ~ 5) cluster as a simulation platform for probing the local structure of the material. Our aim was to scrutinize the magnetism, electronic characteristics influenced by the varying metal atoms within these clusters. A systematic exploration involved incrementing the number of metal atoms and expanding the cluster size to elucidate the corresponding property variations. Density functional theory (DFT) calculations were pivotal to our methodology, employing the B3LYP hybrid functional implemented in the Gaussian 16 software package. The ConMoP (n = 1 ~ 5) cluster underwent optimization calculations and vibrational analysis at the def2-tzvp quantization level, yielding optimized configurations with diverse spin multiplet degrees. To comprehensively characterize and visually represent the stability, electronic features, and catalytic attributes of these configurations, we employed a suite of computational tools. Specifically, quantum chemistry software GaussView and wave function analysis software Multiwfn played integral roles. Through the integrated use of these computational tools, we acquired valuable insights into the magnetism, electronic characteristics of the ConMoP (n = 1 ~ 5) cluster, shedding light on their dependency on distinct metal atoms.

Determination of Ground State Structures of Snx - (x=21-35) Clusters

In this work, an unbiased global search with a homemade genetic algorithm was performed to investigate the structural evolution and electronic properties of Snx - (x=21-35) clusters with density functional theory (DFT) calculations. All the ground-state structures for all these Snx - (x=21-35) clusters have been confirmed by the comparison of the experimental and simulated photoelectron spectra (PESs). It has been revealed that all Snx - (x=21-35) clusters are tricapped trigonal prism (TTP)-based structures consisting of two (for sizes x=21-28) or three (for x=29-35) TTP units, with the remaining atoms adsorbed on the surface or inserted between TTP units. The gradually decreasing HOMO-LUMO gaps indicate that these clusters are undergoing semiconductor-to-metal transformation. The average binding energies show that the structural stabilities of Snx - clusters are not as good as that of silicon and germanium clusters. It found that sizes x=23, 25, 29, 33 show high relative stability.

Reactions of N2O and CO on neutral Rh10On clusters: a density functional study

Density functional theory calculations were performed to identify product, reactant and intermediate dissociative/associative structures for the oxygen abstraction and addition reactions: Rh10On + CO → Rh10On-1 + CO2, n = 1-5 and Rh10On + N2O → Rh10On+1 + N2, n = 0-4 reactions. In the case of the oxygen abstraction reactions, the energetics of the reaction path were very similar in energy regardless of the number of oxygen atoms on the Rh10On cluster, whereas for the addition of oxygen to the Rh10On cluster, the reaction was found to become significantly less exothermic with each successive addition of oxygen.

New Multicomponent Optimization Scheme for Equiatomic Vanadium-Titanium Nanoparticle Study

We present a new multicomponent structure prediction method named MRS-Swap searching, which is inspired by space symmetry and swap of different atomic species. For the pure titanium cluster, a new ground-state structure of the Ti20 cluster with higher symmetry relative to a previous study was found by our method. Based on the structural analysis of Tin (n = 2k, k = 2-11), Vn (n = 2k, k = 2-11), and TinVm (n = m = 2-11) systems, we found that the lowest energy structures of these three systems are very similar, which indicates that equiatomic vanadium-titanium-mixed clusters do not change their ground-state structure relative to the same size pure vanadium and titanium cluster. According to the structure-activity relationship, we conclude that the yield strength (σ) of macro vanadium-titanium alloy is between pure titanium and pure vanadium metal, and this can be expressed through σ(Ti) > σ(TiV) > σ(V). The X-ray diffraction results show that the V2Nb, TiVNb, and Ti2Nb alloys also have the same BCC structure, which may be related to their microstructure. Our method and results can be helpful for future multicomponent alloy design.

Stability, electronic and catalytic properties of ConMoP(n = 1 ~ 5) clusters: A DFT study

CONTEXT

The investigation of the stability, electronic properties, and catalytic activity of clusters Co_nMoP holds significant applications and implications in catalyst design, materials science, energy conversion and storage, and environmental protection. The study aims to delve into the unique features of the clusters Co_nMoP(n = 1 ~ 5), aiming to drive advancements in these related fields. The results obtained from the analysis revealed the stable configurations of the ten clusters, primarily characterized by steric structures. Furthermore, the energy of the clusters was found to increase continuously during growth, as indicated by calculations of atomic fragmentation energy and atomic binding energy. The researchers conducted an analysis of the Natural Population Analysis(NPA) charge, which revealed that Co atoms acted as electron donors, while P and Mo atoms acted as electron acceptors within the clusters. Additionally, an examination of the electrostatic potential indicated that Co and Mo atoms displayed nucleophilic tendencies, while P atoms exhibited electrophilic characteristics. Moreover, the density of states curves, HOMO and LUMO orbitals, and Kooperman's theorem were applied to the clusters Co_nMoP(n = 1 ~ 5).Through this study, a deeper understanding of the properties and behavior of clusters Co_nMoP has been achieved, shedding light on their potential as catalysts. The findings contribute to the existing knowledge of these clusters and provide a basis for further research and exploration in this field.

METHODS

In this study, we employed the clusters Co_nMoP(n = 1 ~ 5) to simulate the local structure of the material, enabling us to investigate the stability, electronic properties, and catalytic properties influenced by the metal atoms. By systematically increasing the number of metal atoms and expanding the cluster size, we explored the variations in these properties. Density functional theory (DFT) calculations were performed using the B3LYP hybrid functional implemented in the Gaussian09 software package. The clusters Co_nMoP(n = 1 ~ 5) underwent optimization calculations and vibrational analysis at the def2-tzvp quantization level, resulting in optimized configurations with different spin multiplet degrees. For data characterization and graphical representation of the stability, electronic properties, and catalytic properties of the optimized configurations, we utilized a range of computational tools. Specifically, the quantum chemistry software GaussView, wave function analysis software Multiwfn were employed. Through the comprehensive utilization of these computational tools, we gained valuable insights into the stability, electronic properties, and catalytic properties of the clusters Co_nMoP(n = 1 ~ 5) and their dependence on different metal atoms.

Metal nanoclusters: from fundamental aspects to electronic properties and optical applications

Monolayer-protected noble metal clusters, also called nanoclusters, can be produced with the atomic precision and in large-scale quantity and are playing an increasingly important role in the field of nanoscience. To outline the origin and the perspectives of this new field, we overview the main results obtained on free metal clusters produced in gas phase including mainly electronic properties, the giant atom concept, the optical properties, briefly the role of the metal atom (alkali, divalent, noble metal) and finally the atomic structure of clusters. We also discuss the limitations of the free clusters. Then, we describe the field of monolayer-protected metal clusters, the main results, the new offered perspectives, the added complexity, and the role of the ligand beyond the superatom concept.

Molecular adsorption on coinage metal subnanoclusters: A DFT+D3 investigation

Gold and silver subnanoclusters with few atoms are prominent candidates for catalysis-related applications, primarily because of the large fraction of lower-coordinated atoms exposed and ready to interact with external chemical species. However, an in-depth energetic analysis is necessary to characterize the relevant terms within the molecular adsorption process that can frame the interactions within the Sabatier principle. Herein, we investigate the interaction between Ag_n and Au_n subnanoclusters (clu, n = 2-7) and N₂ , NO, CO, and O₂ molecules, using scalar-relativistic density functional theory calculations within van der Waals D3 corrections. The onefold top site is preferred for all chemisorption cases, with a predominance of linear (≈180°) and bent (≈120°) molecular geometries. A larger magnitude of adsorption energy is correlated with smaller distances between molecules and clusters and with the weakening of the adsorbates bond strength represented by the increase of the equilibrium distances and decrease of molecular stretching frequencies. From the energetic decomposition, the interaction energy term was established as an excellent descriptor to classify subnanoclusters in the adsorption/desorption process concomitant with the Sabatier principle. The limiting cases: (i) weak molecular adsorption on the subnanoclusters, which may compromise the reaction activation, where an interaction energy magnitude close to 0 eV is observed (e.g., physisorption in N₂ /Ag₆ ); and (ii) strong molecular interactions with the subnanoclusters, given the interaction energy magnitude is larger than at least one of the individual fragment binding energies (e.g., strong chemisorption in CO/Au₄ and NO/Au₄ ), conferring a decrease in the desorption rate and an increase in the possible poisoning rate. However, the intermediate cases are promising by involving interaction energy magnitudes between zero and fragment binding energies. Following the molecular closed-shell (open-shell) electronic configuration, we find a predominant electrostatic (covalent) nature of the physical interactions for N₂ ⋯clu and CO ⋯clu (O₂ ⋯clu and NO⋯clu), except in the physisorption case (N₂ /Ag₆ ) where dispersive interaction is dominant. Our results clarify questions about the molecular adsorption on subnanoclusters as a relevant mechanistic step present in nanocatalytic reactions.

Properties of Naked Silver Clusters with Up to 100 Atoms as Found with Embedded-Atom and Density-Functional Calculations

The structural and energetic properties of small silver clusters Agn with n = 2-100 atoms are reported. For n = 2-100 the embedded atom model for the calculation of the total energy of a given structure in combination with the basin-hopping search strategy for an unbiased structure optimization has been used to identify the energies and structures of the three energetically lowest-lying isomers. These optimized structures for n = 2-11 were subsequently studied further through density-functional-theory calculations. These calculations provide additional information on the electronic properties of the clusters that is lacking in the embedded-atom calculations. Thereby, also quantities related to the catalytic performance of the clusters are studied. The calculated properties in comparison to other available theoretical and experimental data show a good agreement. Previously unidentified magic (i.e., particularly stable) clusters have been found for n>80. In order to obtain a more detailed understanding of the structural properties of the clusters, various descriptors are used. Thereby, the silver clusters are compared to other noble metals and show some similarities to both copper and nickel systems, and also growth patterns have been identified. All vibrational frequencies of all the clusters have been calculated for the first time, and here we focus on the highest and lowest frequencies. Structural effects on the calculated frequencies were considered.

Emerging d-d orbital coupling between non-d-block main-group elements Mg and I at high pressure

d-d orbital coupling, which increases anisotropic and directional bonding, commonly occurs between d-block transition metals. Here, we report an unexpected d-d orbital coupling in the non-d-block main-group element compound Mg2I based on first-principles calculations. The unfilled d orbitals of Mg and I atoms under ambient conditions become part of the valence orbitals and couple with each other under high pressures, resulting in the formation of highly symmetric I-Mg-I covalent bonding in Mg2I, which forces the valence electrons of Mg atoms into the lattice voids to form interstitial quasi-atoms (ISQs). In turn, the ISQs highly interact with the crystal lattice, contributing to lattice stability. This study greatly enriches the fundamental understanding of chemical bonding between non-d-block main-group elements at high pressures.

Enhancing the Performance of Global Optimization of Platinum Cluster Structures by Transfer Learning in a Deep Neural Network

The global optimization of metal cluster structures is an important research field. The traditional deep neural network (T-DNN) global optimization method is a good way to find out the global minimum (GM) of metal cluster structures, but a large number of samples are required. We developed a new global optimization method which is the combination of the DNN and transfer learning (DNN-TL). The DNN-TL method transfers the DNN parameters of the small-sized cluster to the DNN of the large-sized cluster to greatly reduce the number of samples. For the global optimization of Pt₉ and Pt₁₃ clusters in this research, the T-DNN method requires about 3-10 times more samples than the DNN-TL method, and the DNN-TL method saves about 70-80% of time. We also found that the average amplitude of parameter changes in the T-DNN training is about 2 times larger than that in the DNN-TL training, which rationalizes the effectiveness of transfer learning. The average fitting errors of the DNN trained by the DNN-TL method can be even smaller than those by the T-DNN method because of the reliability of transfer learning. Finally, we successfully obtained the GM structures of Pt_n (n = 8-14) clusters by the DNN-TL method.

Structure and size-dependent vibrational and thermal properties of Ni clusters: A systematic ab initio approach

There is scarce information on the vibrational and thermal properties of small Ni clusters. Here, the outcomes of ab initio spin-polarized density functional theory calculations on the size and geometry effects upon the vibrational and thermal properties of Ni_n (n = 13 and 55) clusters, are discussed. For theses clusters a comparison is presented between the closed shell symmetric octahedral (O_h) and the icosahedral (I_h) geometries. The results indicate that the I_h isomers are lower in energy. Besides, ab initio molecular dynamics runs at T = 300K show that Ni₁₃ and Ni₅₅ clusters transform from their initial O_h geometries towards the corresponding I_h ones. For Ni₁₃, we also consider the lowest energy less symmetric layered 1-3-6-3 structure, and the cuboid, recently observed experimentally for Pt₁₃, which is competitive in energy but is unstable, as phonon analysis reveals. We calculate their vibrational density of states (νDOS) and heat capacity, and compare with the Ni FCC bulk counterpart. The characteristic features of the νDOS curves of these clusters are interpreted in terms of the clusters' sizes, the interatomic distance contractions, the bond order values as well as the internal pressure and strains of the clusters. We find that the softest possible frequency of the clusters is size and structure-dependent, being the smallest for the O_h ones. We identify mostly shear, tangential type displacements involving mainly surface atoms for the lowest frequency of the spectra of both I_h and O_h isomers. For the maximum frequencies of these clusters the central atom shows anti-phase movements against groups of nearest neighbor atoms. An excess of heat capacity at low temperatures with respect to the bulk is found, while at high temperatures a constant limiting value, close but lower to the Dulong and Petit value, is determined.

An Open-Shell Superatom Cluster Ta10- with Enhanced Stability by United d-d π Bonds and d-Orbital Superatomic States

We carried out a comprehensive study on the gas-phase reactions of Ta_n^- (n = 5-27) with nitrogen using a customized reflection time-of-flight mass spectrometer coupled with a velocity map imaging apparatus (Re-TOFMS-VMI). Among the studied tantalum clusters, Ta₁₀^- exhibits prominent mass abundance indicative of its unique inertness. DFT calculation results revealed a D_4d bipyramidal prolate structure of the most stable Ta₁₀^-, which was verified by photoelectron spectroscopy experiments. The calculations also unveiled that Ta₁₀^- has the largest HOMO-LUMO gap and second-order difference of binding energy among the studied clusters. This is associated with its well-organized superatomic orbitals, which consist of both 6s and 5d orbitals of tantalum atoms, allowing for splitting of superatomic 1D and 2P orbitals and an enlarged gap between the singly occupied molecular orbital (SOMO) and unoccupied β counterpart, which brings forth stabilization energy pertaining to Jahn-Teller distortion. Also, the SOMO exhibits a united d-d π orbital pattern that embraces the central Ta₈^- moiety.

Ionization potentials of MgN (N = 7-56) clusters formed by spontaneous collapse of magnesium foam in helium nanodroplets

The ionization potentials of magnesium clusters (Mg_N, N = 7-56) are determined by doping ultracold helium nanodroplets (He_M, M ≈ 52 000) with Mg atoms. Inspecting the particle size distributions resulting from non-resonant, short-wavelength, single-photon ionization gives evidence that beyond a certain ensemble size, the developing foam structure undergoes a spontaneous collapse on the way to the laser interaction region. As a result, hot Mg clusters form in the relaxation process. The spontaneous collapse manifests in a substantial change in the size distributions, when recording mass spectra at wavelengths shorter than 272 nm. Tracing individual Mg_N signals as a function of laser photon energy allows extraction of size-specific ionization potentials, which for small clusters show a good agreement with results obtained from density functional theory simulations. The further development is compared to calculations based on the liquid drop model. However, even when quantum effects are included, the simple scaling law is not able to reproduce the development of the ionization potentials. The results suggest that small neutral magnesium clusters behave as non-metallic. The comparison to electron affinities and band gaps obtained from photoemission experiments on Mg_N^- provides information on the charge state dependence of the non-metal-to-metal transition and properties like the Mulliken electron negativity.

Ab initio DFT simulation of electronic and magnetic properties of Tin+1 and FeTin clusters

We report a computational investigation of the electronic and magnetic properties of neutral Ti_n+1 and FeTi_n (n = 1-10) clusters using ab initio calculations based on density functional theory (DFT) within the generalized gradient approximation (GGA). The best structures for Ti_n+1 and FeTi_n clusters are planar for size n < 5, while from n = 5, they showed a compact three-dimensional cage structure. For the best structures of the FeTi_n clusters, the Fe atoms favor the peripheral position with the highest coordination with the neighboring Ti atoms. The evolution as a function of the size of the average binding energies (Eb/atom) and HOMO-LUMO gaps of Ti_n+1 and FeTi_n (n = 1-10) clusters are studied. The stability results show that the Ti_n+1 clusters have relatively higher stability than the FeTi_n cluster with the same size. In addition, the vertical ionization potentials and electron affinities, chemical hardness, and atomic magnetic moment of Ti_n+1 and FeTi_n (n = 1-10) clusters are also investigated.

Multi-atom cluster catalysts for efficient electrocatalysis

This review presents recent developments in the synthesis, modulation and characterization of multi-atom cluster catalysts for electrochemical energy applications.

Role of the OH-group in the adsorption properties of methanol, ethanol, and ethylene glycol on 15-atom 3d, 4d, and 5d transition-metal clusters

The adsorption of alcohols on transition-metal (TM) substrates has received the attention of many researchers due to the applications of alcohols in several technological fields. However, our atomic-level understanding is still far from satisfactory, in particular for the interaction of alcohols with finite-size TM clusters, where new effects can arise due to the presence of quantum-size effects. In this work, we report a theoretical investigation of the adsorption properties of methanol, ethanol, and ethylene glycol on 12 different 3d, 4d, and 5d TM₁₅ clusters based on density functional theory calculations within the semi-empirical D3 van der Waals corrections. From the correlation analysis of all the lowest- and high-energy configurations, we identified the adsorption modes of methanol, ethanol, and ethylene glycol on the TM₁₅ clusters, in which the OH group binds to the cationic TM sites via the O-TM and H-TM interactions. Due to the relatively weak alcohol-TM₁₅ interaction, the changes induced on the TM₁₅ clusters are small, except for Au₁₅ and Ru₁₅, where the bare cluster changes its structure to a nearby minimum in the potential energy surface. The adsorption energy for the alcohol/TM₁₅ systems is correlated to the combination of several parameters, in which the main contribution is connected with the O-TM interaction and the HOTM angles. Furthermore, the TM electronegativity is an important descriptor for the methanol and ethanol adsorption energies, while charge transfer is important for ethylene glycol.

Investigation of the Stability Mechanisms of Eight-Atom Binary Metal Clusters Using DFT Calculations and k-means Clustering Algorithm

Here, we report density functional theory calculations combined with the k-means clustering algorithm and the Spearman rank correlation analysis to investigate the stability mechanisms of eight-atom binary metal AB clusters, where A and B are Fe, Co, Ni, Cu, Ga, Al, and Zn (7 unary and 21 binary clusters). Based on the excess energy analysis, the six most stable binary clusters are NiAl, NiGa, CoAl, FeNi, NiZn, and FeAl, and except for FeNi, their highest energetic stabilities can be explained by the hybridization of the d- and sp-states, which is maximized at the 50% composition, i.e., A₄B₄. Based on the Spearman correlation analysis, the energetic stability of the binary clusters increases with an increase in the highest occupied molecule orbital-lowest unoccupied molecular orbital (HOMO-LUMO) energy separation, which can be considered as a global descriptor. Furthermore, reducing the total magnetic moment values increases the stability for binary clusters without the Fe, Co, and Ni species, while the binary FeB, CoB, and NiB clusters increase their energetic stability with a decrease in the cluster radius, respectively, i.e., an energetic preference for compact structures.

[Ni8(CNtBu)12][Cl]: A nickel isocyanide nanocluster with a folded nanosheet structure

The reaction of 1.75 equiv of ^tBuNC with Ni(1,5-COD)₂, followed by crystallization from benzene/pentane, resulted in the isolation of [Ni₈(CN^tBu)₁₂][Cl] (2) in low yields. Similarly, the reaction of Ni(1,5-COD)₂ with 0.6 equiv of [Ni(CN^tBu)₄], followed by addition of 0.08 equiv of I₂, resulted in the formation of [Ni₈(CN^tBu)₁₂][I] (3), which could be isolated in 52% yield after work-up. Both 2 and 3 adopt folded nanosheet structures in the solid state, characterized by two symmetry-related planar Ni₄ arrays, six terminally bound ^tBuNC ligands, and six ^tBuNC ligands that adopt bridging coordination modes. The metrical parameters of the six bridging ^tBuNC ligands suggest that they have been reduced to their [^tBuNC]^2- form. In contrast to the nanosheet structures observed for 2 and 3, gas phase Ni₈ is predicted to feature a compact bisdisphenoid ground state structure. The strikingly different structural outcomes reveal the profound structural changes that can occur upon addition of ligands to bare metal clusters. Ultimately, the characterization of 2 and 3 will enable more accurate structural predictions of ligand-protected nanoclusters in the future.

Stability Changes in Iridium Nanoclusters via Monoxide Adsorption: A DFT Study within the van der Waals Corrections

Small iridium nanoclusters are prominent subnanometric systems for catalysis-related applications, mainly because of a large surface-to-volume ratio, noncoalescence feature, and tunable properties, which are completely influenced by the number of atoms, geometry, and molecular interaction with the chemical environment. Herein, we investigate the interaction between Ir_n nanoclusters (n = 2-7) and polluting molecules, CO, NO, and SO, using van der Waals D3 corrected density functional theory calculations. Starting from a representative structural set, we determine the growth pattern of the lowest energy unprotected Ir_n nanoclusters, which is based on open structural motifs, and from the adsorption of a XO (X = C, N, and S) molecule, the preferred high-symmetric adsorption sites were determined, dominated by the onefold top site. For protected systems, 4XO/Ir₄ and 6XO/Ir₆, we found a reduction in the total magnetic moment, while the equilibrium bonds of the nanoclusters expanded (contracted) due to mCO and mNO (mSO) adsorption, with exceptions for systems with large structural distortions (4SO/Ir₄ and 6NO/Ir₆). Meanwhile, the C-O and N-O (S-O) bond strength decreases (increases) following an increase (decrease) in the C-O and N-O (S-O) distances upon adsorption. We show, through energetic analysis, that for the different chemical environments, relative stability changes occur from the most stable unprotected nanoclusters, planar square (Ir₄), and prism (Ir₆) to higher energy isomers. The change in the stability order between the two competing protected systems is feasible if the balance between the interaction energy (additive term) and distortion energies (nonadditive terms) compensates for the relative total energies of the unprotected configurations. For all systems, the interaction energy is the main reason responsible for stability alterations, except for 4SO/Ir₄, where the main contribution is from a small penalty due to Ir₄ distortions upon adsorption, and for 4NO/Ir₄, where the energetic effects from the adsorption do not overcome the difference between the binding energies of the unprotected nanoclusters. Finally, from energy decomposition and Hirshfeld charge analysis, we find a predominant covalent nature of the physical contributions in mOX···Ir_n interactions with a cationic core (Ir_n) and an anionic shell (XO coverage).

Energy Decomposition to Access the Stability Changes Induced by CO Adsorption on Transition-Metal 13-Atom Clusters

Our atomistic understanding of the physical-chemical parameters that drives the changes in the relative stability of clusters induced by adsorbed molecules is far from satisfactory. In this work, we employed density functional theory calculations to address this problem using CO adsorption on 13-atom transition-metal clusters, TM₁₃, namely, nCO/TM₁₃, where TM = Ru, Rh, Pd, and Ag, and n = 1-6. Unexpectedly, changes in the relative stability take place for all systems at a lower coverage, namely, at n = 3 (Ru₁₃), 4 (Rh₁₃, Ag₁₃), and 2 (Pd₁₃). To address the effects that lead to changes in the stability, we proposed an energy decomposition scheme for the binding energy of the nCO/TM₁₃ systems, which yields that the change in relative stability is dominated by the interaction energy and cluster distortion energy upon adsorption, where the interaction energy is higher for high-energy unprotected clusters. Furthermore, we characterized all adsorption parameters, which helps us to complement our atomistic understanding.

DFT insights into structural effects of Ni-Cu/CeO2 catalysts for CO selective reaction towards water-gas shift

The water-gas shift (WGS) reaction is a key step in hydrogen production, particularly to meet the high-purity H2 requirement of PEM fuel cells. The catalysts currently employed in large-scale WGS plants require a two-step process to overcome thermodynamic and kinetic limitations. Ni-Cu/CeO2 solids are promising catalysts for the one-step process required for small-scale applications, as the addition of Cu hinders undesired methanation reactions occurring on Ni/CeO2. In this work, we performed calculations on Ni4-xCux/CeO2(111) systems to evaluate the influence of cluster conformation on the selectivity towards water-gas shift. The structure and miscibility of CeO2-supported Ni4-xCux clusters were investigated and compared with those of gas-phase clusters to understand the effect of metal-support interactions. The adsorption of CO onto apical Ni and Cu atoms of Ni4-xCux/CeO2(111) systems was studied, and changes in the C-O bond strength were confirmed at the electronic level by investigating shifts in the 3σ and 1π orbitals. The selectivity towards WGS was evaluated using Brønsted-Evans-Polanyi relations for the C-O activation energy. Overall, a strengthening of the C-O bond and an increase in CO dissociation energy were verified on Cu-containing clusters, explaining the improvement in selectivity of Ni4-xCux/CeO2(111) systems.

Vanadium Cluster Neutrals Reacting with Water: Superatomic Features and Hydrogen Evolution in a Fishing Mode

Hydrogen evolution reaction (HER) is known as the heart of various energy storage and conversation systems of renewable energy sources. Here we observe the cluster reactions of a light transition metal, vanadium, with water in a gas-phase flow tube reactor. While HER products of V₁ and V₂ were not observed, the effective HER of water on neutral V_n (n ≥ 3) clusters reveals reasonable and size-dependent reactivity of the vanadium clusters. Superatomic features and reaction dynamics of V₁₀, V₁₃, and V₁₆ are highlighted. Among the three typical superatoms, V₁₀ and V₁₆ exhibit an abnormal superatomic orbital energy level order, 1S|2S|1P|1D..., where the energy-reduced 2S orbital helps to accommodate the geometric structure and hence reinforce the cluster stability. In comparison, V₁₃ bears a less symmetrical structure and reacts readily with water, allowing for recombination of a hydroxyl atom with an adsorbed hydrogen atom, akin to a fishing-mode HER process. The joint experimental and theoretical study on neutral V_n clusters clarifies the availability of superatom chemistry for transition metals and appeals further development of cluster theory based on electronic cloud/orbital analysis instead of simply counting the valence electrons. Also, we provide insights into the HER mechanism of metal clusters and propose a strategy to design new materials for portable fuel cells of hydrogen energy.

Gas-phase preparation and the stability of superatomic Nb11O15. Phys Chem Chem Phys 2021;23:15766-15773. [PMID: 34286767 DOI: 10.1039/d1cp02128a]

We report a study of the reactions of pure metal clusters Nbn- with dioxygen in the gas phase. It is found that the presence of low-concentration dioxygen reactants results in oxygen-addition products, whereas sufficient high-concentration dioxygen enables oxygen-etching reactions giving rise to molecular niobium oxides. Interestingly, in the presence of a suitable gas flow rate of an intermediate dioxygen concentration, a highly selective product Nb11O15- shows up in the mass spectra. Utilizing density functional theory (DFT) calculations, we have discussed the reactivities of Nbn- (3 ≤ n ≤ 14) clusters with oxygen, and unveiled the reasonable stability of Nb11O15- pertaining to its unique geometric structure with a D5h Nb@Nb10 core fully protected by 15 bridge-oxygen atoms. The oxygen-passivated Nb@Nb10O15- cluster exhibits a large HOMO-LUMO gap (1.46 eV) and effective multicenter bonds with remarkable superatom orbitals for all the 26 valence electrons of the Nb@Nb10 core corresponding to well-staggered energy levels. We illustrate the superatomic features in the Nb@Nb10 metallic core for which the adaptive natural density partitioning (AdNDP) analysis unveils thirteen 11c-2e bonds. Among them, one of the 11c-2e bonds accounts for the superatomic S orbital, three bonds correspond to superatomic P orbitals, another five display vivid D orbital characteristics, and the remaining four 11c-2e bonds are assigned to F orbital features. In addition, the net atomic charge of the center Nb atom is as high as -0.804 |e| rendering core-shell electrostatic interactions and the shielding effect of the Nb10O15 shell.

Carbon Monoxide Activation on Small Iron Magnetic Cluster Surfaces, FenCO, n = 1-20. A Theoretical Approach

The chemical activation of the carbon monoxide (CO) molecule on the surface of iron clusters Fe_n (n = 1-20) is studied in this work. By means of density functional theory (DFT) all-electron calculations, we have found that the adsorption of CO over the bare magnetic Fe_n (n = 1-20) clusters is thermochemically favorable. The Fe_n-CO interaction increases the C-O bond length, from 1.128 ± 0.014 Å, for isolated CO, up to 1.251 Å, for Fe₉CO. Also, the calculated wavenumbers associated with the stretching modes ν_CO are decreased, or red-shifted, as another indicator of the CO bond weakening, passing from 2099 ± 4 to 1438 cm^-1. Markedly, wavenumbers of vibrational modes ν_CO agree admirably well in comparison with experimental results reported for Fe_nCO (n = 1, 18-20), getting small errors below 2.6%. The C-O bond is enlarged on the Fe_nCO (n = 1-20) composed systems, as the CO molecule increases its bonding, charge transference, and coordination with the iron cluster. Therefore, small bare iron particles Fe_n (n = 1-20) can be proposed to promote the CO dissociation, especially Fe₉CO, which has been proven to obtain the most prominent activation of the strong C-O bond by means of the charge transference from the metal core.

Nanoalloys of Metals Which Do Not Form Bulk Alloys: The Case of Ag-Co

Ag and Co metals do not form macroscopic solid or liquid alloys. However, Ag_mCo_n clusters have been produced in dual-target dual-laser vaporization experiments, and the same occurs for other pairs of immiscible metals. We have performed density functional calculations to shed light on this phenomenon. The main result, obtained for clusters with sizes m + n not larger than 11, is that the cohesive energies justify that those clusters can be formed starting from free Ag and Co atoms, which is the case in the vaporization experiments. At the same time, mixing of Ag and Co is difficult even at the nanoscale. This is revealed by the application of several miscibility criteria. Those two features become, nevertheless, compatible in the clusters. Even if the cluster sizes considered are small, the emerging trend becomes clear: Co atoms form a core in the inner part, surrounded by a shell of Ag atoms on the surface. A consequence is that core-shell clusters can be formed from pairs of metals that do not form macroscopic alloys. The magnetic moments μ of the Ag_mCo_n clusters are due mainly to the Co atoms, and the presence of Ag induces a reduction in the magnitude of μ.

Theoretical study of interaction of Fe13O8@Zn48O48 cluster with dopamine: Magnetic and optical properties

In this study, we modelled the interaction of Fe₁₃O₈ and Fe₁₃O₈@Zn₄₈O₄₈ (core@shell) cluster with a biologically active dopamine molecule using density functional theory. First, the electronic, magnetic and optical properties of core@shell, Fe₁₃O₈@Zn₄₈O₄₈ cluster investigated and compared with isolated Fe₁₃O₈ and Zn₄₈O₄₈ clusters. Fe₁₃O₈@Zn₄₈O₄₈ cluster is found to be energetically stable. For Fe₁₃O₈ and Fe₁₃O₈@Zn₄₈O₄₈ clusters have the net magnetic moment 42 μ_B. The decrease in HOMO-LUMO gap of core@shell cluster as compared to that of isolated clusters reflects the higher reactivity. The results of the site dependent interaction of Fe₁₃O₈ and Fe₁₃O₈@Zn₄₈O₄₈ clusters with dopamine molecule are presented. The interaction strength is determined in terms of the cluster-dopamine complex binding energy and found to be enhanced for core@shell cluster than the Fe₁₃O₈. Furthermore, the calculated results predict that in presence of dopamine, the magnetic moment of Fe₁₃O₈ and Fe₁₃O₈@Zn₄₈O₄₈ cluster remains unaffected. The analysis of optical spectra of core@shell indicates the obvious red shift compared to Zn₄₈O₄₈ clusters. The optical spectra of Fe₁₃O₈@Zn₄₈O₄₈-dopamine shows the higher oscillator strength as compared to that of Fe₁₃O₈-dopamine complex. Fe₁₃O₈-dopamine complex gives rise to more quenched oscillator strengths as compared to that of bare iron oxide cluster. These results indicate interesting magneto-optical behaviour, which can be useful for biomedical applications.

Ab initio investigation of quantum size effects on the adsorption of CO2, CO, H2O, and H2 on transition-metal particles

Adsorption is a crucial preliminary step for the conversion of CO2 into higher-value chemicals, nonetheless, the atomistic understanding of how substrate particle size affects this step is still incomplete. In this study, we employed density functional theory to investigate the effects of particle size on the adsorption of model molecules involved in the CO2 transformations (CO2, CO, H2O and H2) on Con, Nin and Cun particles with different sizes (n = 13, 55, 147) and on the respective close-packed surfaces. We found significant size-dependence of the adsorption properties for physisorbed (linear) and chemisorbed (bent) CO2 on the substrates and distinct (symmetric or asymmetric) stretching of the C-O bonds, which can play a crucial role to understand the CO2 dissociation pathways. For CO and H2, some properties showed small oscillations, due to size effects that induced alternation of the adsorption site preference for different particle sizes; for H2O, the adsorption properties were almost independent of particle size. The presence of low-coordinated adsorption sites resulted in a trend for stronger adsorption and greater charge transfer for smaller clusters. Fixing the size-independent factors (e.g., type of metal), our results show that CO2 adsorption on transition-metal clusters is significantly affected by particle size, suggesting that substrate particle size could be a key factor to understand and control the catalytic transformations of CO2.

A small bimetallic Ag3Cu2 nanocluster with dual emissions within and against Kasha's rule

Single crystals of a small bimetallic Ag₃Cu₂ nanocluster protected by six ligands of 2,4-dimethylbenzene thiol are synthesized by a one-pot procedure of wet chemistry. This Ag₃Cu₂ nanocluster bears a trigonal bipyramid metallic core with two copper atoms located on both sides of a triangular Ag₃. Interestingly, the six Cu-Ag side edges of the trigonal bipyramid are fully protected by the six ligands giving rise to reinforced stability and high chemical purity. More interestingly, this Ag₃Cu₂ cluster shows strong dual fluorescence emissions in both ultraviolet visible (UV-vis) and near infrared (NIR) regions. Theoretical calculations reproduce the absorption and fluorescence spectra where the NIR emission at 824 nm is assigned to the S₁→ S₀ transition, while the simultaneous emission in the visible band is due to the radiation of highly excited states and is against Kasha's rule.

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.

Size and chemical order dependence of magnetic-ordering temperature and spin structure in Fe@Ni and Ni@Fe core-shell nanoparticles

The effect of particle size and chemical order on the temperature-dependent magnetic properties of Fe@Ni and Ni@Fe core-shell nanoparticles is studied in the framework of a classical spin Hamiltonian and Monte Carlo simulations. We found that the mean temperature-dependent magnetization and magnetic-ordering temperature are strongly affected by both the particle size (in size range of 4 to 16 nm) and core-shell chemical order. As a main result, we report the depression of the magnetic ordering-temperature with decreasing size of the elemental Fe and Ni nanoparticles. More specifically, in the case of Fe and Ni nanoparticles, the magnetic-ordering temperature is lowered by 40 (195 K) to 300 (175 K) compared to the bulk value for nanoparticle diameters ranging from 16 to 4 nm, respectively, consistent with previous theoretical data. We further provide a comprehensive insight into the magnetic properties of Fe@Ni and Ni@Fe nanoparticles, unveiling a rich and distinct magnetic-ordering temperature and spin structure that emphatically depends on the core/shell ratio.

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-small iridium nanoparticles as active catalysts for the selective and efficient reduction of nitroarenes

The ultra-small noble metal iridium nanoparticles (IrNPs) possessing super catalytic activity can be applied in the efficient and selective catalytic reduction of nitroarenes under mild reaction conditions for the first time.

A theoretical indicator of transition-metal nanoclusters applied in the carbon nanotube nucleation process: a DFT study

The center of gravity of the occupied d-states for nanoclusters is obtained as a good indicator to reveal the best candidates to the interaction with the carbon nanotubes.

CO, NO, and SO adsorption on Ni nanoclusters: a DFT investigation

Nickel nanoclusters are very promising for catalysis-related applications, especially involving chemical reactions with polluting molecules, such as carbon, nitrogen, and sulfur monoxides, which are directly or indirectly involved in serious environmental pollution problems. Therefore, it is of utmost importance to improve the understanding of the interaction between Ni nanoclusters and diatomic molecules, such as CO, NO, and SO, to provide insights into real subnano catalysts. Thus, here, we report an ab initio investigation based on density functional theory calculations within van der Waals D3 corrections to investigate the adsorption properties of CO, NO, and SO on Ni nanoclusters. From energetic and electronic criteria applied to Nin nanoclusters (n = 2-15), we selected Ni6 (octahedron) and Ni10 (triangular pyramid) nanoclusters as supports. According to our analyses, the molecular adsorption increases the stability of Ni nanoclusters, especially for Ni6 systems. The interaction intensity is larger for SO than for NO and CO in adsorbed systems, and the strong OS-Ni interaction is responsible for the well-known sulfur poisoning on transition-metal systems. The lowest energy adsorption sites are onefold for CO/Ni6, NO/Ni6, and CO/Ni10; twofold for NO/Ni10; and threefold for SO/Ni6 and CO/Ni10, where CO and NO molecules sustain linear and perpendicular geometries, while SO geometry changes to a bent configuration resulting from a sideways adsorption. The equilibrium bond lengths of the molecules expand upon adsorption, from 0.9% (NO/Ni6/10) to 11.3% (SO/Ni6/10), consequently, the internal molecular bond strengths decrease, since there is a reduction in the molecular stretching frequencies. This result occurs most strongly for SO followed by NO and CO systems, which was confirmed by an estimation of the energetic contribution of the distortion after the adsorption process. Thus, the strong S-Ni interaction, given by SO chemisorption on hollow sites with a sideways interaction, implies an energetic decrease and, consequently, a part of the energy gained from the SO-Ni interaction is from the SO and nanocluster distortions. Ultimately, using the energy decomposition analysis (from SAPT0) for XO/Ni6 systems, we improved the understanding of the CO and NO (SO) singlet (doublet) spin multiplicities' interaction with Ni6 nanoclusters.