Investigation of niobium clusters: Bare and CO-adsorption

Computational modeling of the Nb 4 -CO chemisorption process

CONTEXT

The transition metal niobium (Nb) has attracted considerable attention from the scientific community due to its intriguing electronic properties and applications in catalysts suitable for chemical reactions. Thus, this work investigates the adsorption of the atmospheric polluting gas carbon monoxide (CO) by the niobium cluster (Nb4 ), to describe the reactive nature of Nb4 . This entire study was carried out by applying the Coupled-Cluster method and Density Functional Theory (through the HSE06 functional) and the def2-QZVP plus Def2-TZVP/C auxiliary basis set functions. The results of electronic structure calculations and IR vibrational spectra suggest that both Nb4 and the Nb4 -CO clusters can be considered stable. Furthermore, the obtained results also indicate that there is a chemisorption of carbon monoxide by the Nb4 niobium cluster. This feature can serve as motivation for future theoretical-experimental studies, as it suggests that the Nb4 cluster may have possible technological applications in automotive catalytic processes.

METHODS

Initial three-dimensional structures were constructed. Complete optimization of the geometry was performed in coupled cluster and density functional theory methods. From the optimization configuration, it was possible to investigate the stability, chemisorption process, binding energies, charge analysis, molecular orbital energies, and IR vibrational spectra of the systems.

Weak Interactions Initiate C-H and C-C Bond Dissociation of Ethane on Nbn + Clusters

The conversion of ethane into value-added chemicals under ambient conditions has attracted much attention but the mechanisms remain not fully understood. Here we report a study on the reaction of ethane with thermalized Nb_n ⁺ clusters based on a multiple-ion laminar flow tube reactor combined with a triple quadrupole mass spectrometer (MIFT-TQMS). It is found that ethane reacts with Nb_n ⁺ clusters to form both products of dehydrogenation and methane-removal (odd-carbon products). Combined with density functional theory (DFT) calculations, we studied the reaction mechanisms of the C-C bond activation and C-H bond cleavage on the Nb_n ⁺ clusters. It is unveiled that hydrogen atom transfer (HAT) initiates the reaction process, giving rise to the formation of Nb-C bonds and an elongated C-C distance in the HNb_n ⁺ CH₂ CH₃ motif. Subsequent reactions allow for C-C bond activation and a competitive HAT process which is associated with CH₄ removal or H₂ release, resulting in the production of the observed carbides.

Effects of All-Electron Basis Sets and the Scalar Relativistic Corrections in the Structure and Electronic Properties of Niobium Clusters

In this paper, an augmented all-electron double-ζ basis set is used in calculations of the structure and electronic properties of small niobium clusters. The B3PW91 and M06 DFT functionals with and without second order Douglas-Kroll-Hess (DKH) scalar relativistic corrections are also utilized. Furthermore, an additional d Gaussian type function is introduced in the standard basis sets in order to improve the description of the clusters orbitals in the valence band. Our findings show that the extra d function is important to yield accurate results of electronic properties and, in addition, the DKH corrections can be relevant when the all-electron basis sets are used in the calculations. Our best results are obtained with the M06 functional together with the DKH second order corrections and with the extra d function added to the all-electron basis set.

Evolution of the structural, energetic, and electronic properties of the 3d, 4d, and 5d transition-metal clusters (30 TMn systems for n = 2–15): a density functional theory investigation

Subnanometric transition-metal (TM) clusters have attracted great attention due to their unexpected physical and chemical properties, leastwise compared to their bulk counterparts.

Structures, spectra, and energies of niobium clusters from Nb13 to Nb20

A comprehensive theoretical investigation on structures and properties of niobium clusters in the range from 13 to 20 atoms, in three different charged states, is performed by using the BPW91 and M06 functionals and the cc-pVDZ-PP basis set. These species are predicted to prefer low spin ground state, i.e., singlet (for even electron) and doublet (for odd electron) systems. In terms of growth mechanism, a compact structure with one Nb encapsulated by a cage formed from five and six triangles is found to be favored over an icosahedral evolution. Unlike many 3d metals, whose volumes are much smaller, 13 and 19 Nb atoms clusters do not exist as icosahedra and double-icosahedra. A distinct case is Nb(15) as it bears a slightly distorted bcc structure. For some systems, several lower lying isomers are computed to be so close in energy that DFT computations cannot clearly establish their ground electronic states. The existence of structural isomers with comparable energy content is established for Nb(n) species with n = 13, 18, 19, and 20 in both neutral and charged states. The vibrational (IR) spectra are also calculated. While the spectra of smaller systems are strongly dependent on addition or removal of an electron from the neutral, the spectra of the larger size clusters are mostly independent of the charged state. The neutrals and their corresponding ions usually have a quite similar IR pattern. Electron affinities (EA), ionization energies (IE), average binding energies, dissociation energies, and frontier orbital energy gaps are evaluated. The computed EAs and IEs are generally in fair agreement with experiment. The Nb(15) system is observed to be stable and it can form a highly symmetric structure in all charged states with both open and closed electron shells.

Structural, electronic and magnetic effects of Al-doped niobium clusters: a density functional theory study

The application of the ab initio stochastic search procedure with Saunders "kick" method has been carried out for the elucidation of global minimum structures of a series of Al-doped clusters, Nb(n)Al (1 ≤ n ≤ 10). We have studied the structural characters, growth behaviors, electronic and magnetic properties of Nb(n)Al by the density functional theory calculations. Unlike the previous literature reported on Al-doped systems where ground state structures undergo a structural transition from the Al-capped frame to Al-encapsulated structure, we found that Al atom always occupies the surface of Nb(n)Al clusters and structural transition does not take place until n = 10. Note that the fragmentation proceeds preferably by the ejection of an aluminum atom other than niobium atom. According to the natural population analysis, charges always transfer from aluminum to niobium atoms. Furthermore, the magnetic moments of the Nb(n)Al clusters are mainly located on the 4d orbital of niobium atoms, and aluminum atom possesses very small magnetic moments.

Chemistry of (and on) transition metal clusters: a Fourier transform ion cyclotron resonance study of the reaction of niobium cluster cations with nitric oxide

The reactions of niobium cluster cations, Nb(+)(n) (n = 2-19), with nitric oxide have been investigated using Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR). The overall reaction rate constants are found to be in reasonable agreement with collision rates calculated using the surface charge capture model. The dominant reaction for small clusters (n <9) involves reaction-induced fragmentation resulting in the loss of either NbO or NbN. By contrast, the main reaction observed for the larger clusters (n> 11) is sequential NO chemisorption. Clusters n = 9, 10 exhibit both extremes of behaviour and are the only clusters upon which there is evidence of NO decomposition with N(2) loss observed whenever multiple NO molecules are co-adsorbed. The rate constants for each process have been determined as a function of cluster size.

Tuning ferroelectricity of niobium clusters by oxygen chemisorption

The chemisorption of atomic oxygen on Nb(N) (N=2-16) has been investigated by the scalar relativistic all-electron density-functional calculations with emphasis on its effect on the ferroelectricity of Nb(N). We have shown that the binding of O atom to niobium clusters is site- and size-dependent, for which the bridge-site doping is preferred in the smaller size range from N=2 to 8 and the threefold hollow site one for the larger Nb(N) with 9 < or = N < or = 16. Though the geometrical structures of Nb(N) are modified slightly when doped with an oxygen atom, their ferroelectric properties vary considerably, depending on cluster size and the O adsorption sites, which is mainly caused by the charge transfer between the oxygen atom and niobium clusters. The addition of oxygen can enhance the ferroelectricity of Nb(N) with moderate and essentially zero moments while induce only small dipole change for those Nb(N) with large moments. Thus, the big fluctuation of the Nb(N)'s dipole moments with size is greatly suppressed by the O doping, indicating that one might expect to tune the size-dependent ferroelectricity of Nb(N) by the chemical decoration.

Structural isomers and reactivity for Rh6 and Rh6+

The structure, energetics, and interconversion of isomers of Rh(6) and Rh(6)(+) are studied by using density functional theory with Gaussian basis sets, using guess structures derived from basin-hopping simulations, and obtained by using the Sutton-Chen potential. A large range of spin multiplicities is considered for each isomer. Our calculations suggest two low-lying structures as possible structural isomers: a square bipyramid and a trigonal prism. The reactivity of these two candidate structural isomers with respect to adsorption of nitric oxide is studied via location of reaction transition states and calculation of reaction barriers. Similarities and differences with surface reaction studies are highlighted. These data provide powerful evidence that structural isomerism, and not different spin states, is responsible for the observed biexponential reaction kinetics.

Coexistence of ferroelectricity and ferromagnetism in tantalum clusters

The atomic and electronic structures of Ta(N) (N=2-23) clusters have been determined in the framework of pseudopotential density-functional calculations, based upon an unbiased global search with guided simulated annealing to an empirical potential. It is found that the ground-state structures of Ta(N) are very similar to those of Nb(N), showing no preference for the icosahedral growth. Also, a size- and structure-dependent ferroelectricity is found in these tantalum clusters. More importantly, it is found that the ferroelectricity and ferromagnetism can coexist in the homogeneous transition-metal cluster, offering a possibility to obtain a new type of "multiferroic" materials composed of the clusters. Finally, the far-infrared spectroscopy is suggested to be an efficient tool to distinguish the ferroelectric clusters.

Relaxation dynamics and structural isomerism in Nb10 and Nb10+

The structure, energetics, and interconversion of isomers of Nb10 and Nb10+ are studied using density functional theory with Gaussian basis sets, using guess structures derived from basin-hopping simulations with the Finnis-Sinclair [Philos. Mag. A 50, 45 (1984)] potential. These results are used as input to a master equation approach to model the relaxation of these clusters. Ionization potentials are calculated for all relevant minima, as are the infrared spectra. On the basis of these data, and known experimental results, plausible explanations are given for the biexponential reaction kinetics observed for Nb10 and Nb10+ with respect to small molecule adsorbates. In principle, this approach could be extended to investigate any midsized transition metal cluster that exhibits structural isomerism.

Structure and stability of binary transition-metal clusters (NbCo)n (n⩽5): A relativistic density-functional study

Equilibrium geometries and electronic properties of binary transition-metal clusters, (NbCo)n (n < or = 5), have been investigated by means of the relativistic density-functional approach. The metal-metal bonding and stability aspects of these clusters have been analyzed on the basis of calculations. Present results show that these clusters exhibit rich structural varieties on the potential-energy surfaces. The most stable structures have a compact conformation in relatively high symmetry, in which the Nb atoms prefer to form an inner core and Co atoms are capped to the facets of the core. Such building features in clustering of the Nb/Co system are related to the order of bond strength: Nb-Nb>Nb-Co>Co-Co. As the binary cluster size increases, the Nb-Co bond may become stronger than the Nb-Nb bond in the inner niobium core, which results in a remarkable increment of the Nb-Nb bond length. Amongst these binary transition-metal clusters, the singlet (NbCo)4 in T(d) symmetry has a striking high stability due to the presence of the spherical aromaticity and electronic shell closure. The size dependence of the bond length and stability of the cluster has been explored.

Cluster Size Selectivity in the Product Distribution of Ethene Dehydrogenation on Niobium Clusters

Ethene reactions with niobium atoms and clusters containing up to 25 constituent atoms have been studied in a fast-flow metal cluster reactor. The clusters react with ethene at about the gas-kinetic collision rate, indicating a barrierless association process as the cluster removal step. Exceptions are Nb8 and Nb10, for which a significantly diminished rate is observed, reflecting some cluster size selectivity. Analysis of the experimental primary product masses indicates dehydrogenation of ethene for all clusters save Nb10, yielding either Nb(n)C2H2 or Nb(n)C2. Over the range Nb-Nb6, the extent of dehydrogenation increases with cluster size, then decreases for larger clusters. For many clusters, secondary and tertiary product masses are also observed, showing varying degrees of dehydrogenation corresponding to net addition of C2H4, C2H2, or C2. With Nb atoms and several small clusters, formal addition of at least six ethene molecules is observed, suggesting a polymerization process may be active. Kinetic analysis of the Nb atom and several Nb(n) cluster reactions with ethene shows that the process is consistent with sequential addition of ethene units at rates corresponding approximately to the gas-kinetic collision frequency for several consecutive reacting ethene molecules. Some variation in the rate of ethene pick up is found, which likely reflects small energy barriers or steric constraints associated with individual mechanistic steps. Density functional calculations of structures of Nb clusters up to Nb(6), and the reaction products Nb(n)C2H2 and Nb(n)C2 (n = 1...6) are presented. Investigation of the thermochemistry for the dehydrogenation of ethene to form molecular hydrogen, for the Nb atom and clusters up to Nb6, demonstrates that the exergonicity of the formation of Nb(n)C2 species increases with cluster size over this range, which supports the proposal that the extent of dehydrogenation is determined primarily by thermodynamic constraints. Analysis of the structural variations present in the cluster species studied shows an increase in C-H bond lengths with cluster size that closely correlates with the increased thermodynamic drive to full dehydrogenation. This correlation strongly suggests that all steps in the reaction are barrierless, and that weakening of the C-H bonds is directly reflected in the thermodynamics of the overall dehydrogenation process. It is also demonstrated that reaction exergonicity in the initial partial dehydrogenation step must be carried through as excess internal energy into the second dehydrogenation step.

Origin of spontaneous electric dipoles in homonuclear niobium clusters

Surprisingly large, spontaneous electric dipole moments recently observed in homonuclear niobium clusters below 100 K are explained using first-principles electronic structure calculations. The calculated moments for Nb(n) (n < or =15) generally agree with the experimental data. A strong correlation is found between the geometrical asymmetry of the cluster and electric dipole: its magnitude is proportional to the spread in the principal moments of inertia and its direction aligns with the axis of the largest principal moment. Charge deformation densities reveal directional, partially covalent bonds that stabilize structural asymmetry. Classical simulations of the deflection of a cluster in a molecular beam reveal that the electronic dipole may persist at higher temperatures, but is masked by the rotational dynamics of the cluster.