Unraveling Hydrogen Adsorption on Transition Metal-Doped [Mo3S13]2- Clusters: Insights from Density Functional Theory Calculations

Molecular and dissociative hydrogen adsorption of transition metal (TM)-doped [Mo3S13]2- atomic clusters were investigated using density functional theory calculations. The introduced TM dopants form stable bonds with S atoms, preserving the geometric structure. The S-TM-S bridging bond emerges as the most stable configuration. The preferred adsorption sites were found to be influenced by various factors, such as the relative electronegativity, coordination number, and charge of the TM atom. Notably, the presence of these TM atoms remarkably improved the hydrogen adsorption activity. The dissociation of a single hydrogen molecule on TM[Mo3S13]2- clusters (TM = Sc, Cr, Mn, Fe, Co, and Ni) is thermodynamically and kinetically favorable compared to their bare counterparts. The extent of favorability monotonically depends on the TM impurity, with a maximum activation barrier energy ranging from 0.62 to 1.58 eV, lower than that of the bare cluster (1.69 eV). Findings provide insights for experimental research on hydrogen adsorption using TM-doped molybdenum sulfide nanoclusters, with potential applications in the field of hydrogen energy.

Investigation of Structures, Stabilities, and Electronic and Magnetic Properties of Niobium Carbon Clusters Nb7Cn (n = 1-7)

The geometrical structures, relative stabilities, and electronic and magnetic properties of niobium carbon clusters, Nb7Cn (n = 1-7), are investigated in this study. Density functional theory (DFT) calculations, coupled with the Saunders Kick global search, are conducted to explore the structural properties of Nb7Cn (n = 1-7). The results regarding the average binding energy, second-order difference energy, dissociation energy, HOMO-LUMO gap, and chemical hardness highlight the robust stability of Nb7C3. Analysis of the density of states suggests that the molecular orbitals of Nb7Cn primarily consist of orbitals from the transition metal Nb, with minimal involvement of C atoms. Spin density and natural population analysis reveal that the total magnetic moment of Nb7Cn predominantly resides on the Nb atoms. The contribution of Nb atoms to the total magnetic moment stems mainly from the 4d orbital, followed by the 5p, 5s, and 6s orbitals.

Superior Reactivity of Molybdenum-Sulfur Cluster Anions Mo5S2- and Mo5S3- toward Dinitrogen

Inspired by the fact that Mo is a key element in biological nitrogenase, a series of gas-phase MoxSy- cluster anions are prepared and their reactivity toward N2 is investigated by the combination of mass spectrometry, photoelectron imaging spectroscopy, and density functional theory calculations. The Mo5S2- and Mo5S3- cluster anions show remarkable reactivity compared with the anionic species reported previously. The spectroscopic results in conjunction with theoretical analysis reveal that a facile cleavage of N≡N bonds takes place on Mo5S2- and Mo5S3-. The large dissociative adsorption energy of N2 and the favorable entrance channel for initial N2 approaching are proposed as two decisive factors for the superior reactivity of Mo5S2- and Mo5S3-. Besides, the modulation of S ligands on the reactivity of metal centers with N2 is proposed. The highly reactive metal-sulfur species may be obtained by the coordination of two to three sulfur atoms to bare metal clusters so that an appropriate combination of electronic structures and charge distributions can be achieved.

The Oxygen Evolution Reaction at MoS2 Edge Sites: The Role of a Solvent Environment in DFT-Based Molecular Simulations

Density functional theory (DFT) calculations are employed to study the oxygen evolution reaction (OER) on the edges of stripes of monolayer molybdenum disulfide. Experimentally, this material has been shown to evolve oxygen, albeit with low efficiency. Previous DFT studies have traced this low catalytic performance to the unfavourable adsorption energies of some reaction intermediates on the MoS₂ edge sites. In this work, we study the effects of the aqueous liquid surrounding the active sites. A computational approach is used, where the solvent is modeled as a continuous medium providing a dielectric embedding of the catalyst and the reaction intermediates. A description at this level of theory can have a profound impact on the studied reactions: the calculated overpotential for the OER is lowered from 1.15 eV to 0.77 eV. It is shown that such variations in the reaction energetics are linked to the polar nature of the adsorbed intermediates, which leads to changes in the calculated electronic charge density when surrounded by water. These results underline the necessity to computationally account for solvation effects, especially in aqueous environments and when highly polar intermediates are present.

Making Sense of the Growth Behavior of Ultra-High Magnetic Gd2-Doped Silicon Clusters

The growth behavior, stability, electronic and magnetic properties of the Gd₂Si_n^- (n = 3-12) clusters are reported, which are investigated using density functional theory calculations combined with the Saunders 'Kick' and the Artificial Bee Colony algorithm. The lowest-lying structures of Gd₂Si_n^- (n = 3-12) are all exohedral structures with two Gd atoms face-capping the Si_n frameworks. Results show that the pentagonal bipyramid (PB) shape is the basic framework for the nascent growth process of the present clusters, and forming the PB structure begins with n = 5. The Gd₂Si₅^- is the potential magic cluster due to significantly higher average binding energies and second order difference energies, which can also be further verified by localized orbital locator and adaptive natural density partitioning methods. Moreover, the localized f-electron can be observed by natural atomic orbital analysis, implying that these electrons are not affected by the pure silicon atoms and scarcely participate in bonding. Hence, the implantation of these elements into a silicon substrate could present a potential alternative strategy for designing and synthesizing rare earth magnetic silicon-based materials.

Trimetallic clusters in the sumanene bowl for dinitrogen activation

It is of great importance to find catalysts for the nitrogen reduction reaction (NRR) with high stability and reactivity. A series of M₃ clusters (M = Ti, Zr, V, and Nb) supported on sumanene (C₂₁H₁₂) were designed as potential catalysts for the NRR by taking advantage of the high reactivity of trimetallic clusters and the unique geometric and electronic properties of sumanene, a bowl-like organic molecule. Detailed mechanisms of NN bond cleavage on C₂₁H₁₂-M₃ were investigated by DFT calculations and compared with those on bare M₃ clusters. M₃ in the sumanene bowl is very stable with large binding energies, which prohibits the cohesion of M₃ into M₆. In the bowl, M₃ has a (quasi-) equilateral triangle structure with lengthened M-M bonds, which is particularly beneficial to the N₂ transfer process on Ti₃ and V₃ clusters. The N-N bond can be dissociated by both M₃ and C₂₁H₁₂-M₃ clusters without the overall energy barriers. A blurring effect is found in which some geometric and electronic properties of different metal types become similar when M₃ is supported on the substrate. Our work demonstrates that sumanene is a suitable substrate to support M₃ in the activation of N₂ with enhanced stability and maintained a high level of reactivity compared to bare M₃.

Facile N≡N Bond Cleavage by Anionic Trimetallic Clusters V3-x Tax C4 - (x=0-3): A DFT Study

Activation of N₂ on anionic trimetallic V_3-x Ta_x C₄ ^- (x=0-3) clusters was theoretically studied employing density functional theory. For all studied clusters, initial adsorption of N₂ (end-on) on one of the metal atoms (denoted as Site 1) is transferred to an of end-on: side-on: side-on coordination on three metal atoms, prior to N₂ dissociation. The whole reaction is exothermic and has no global energy barriers, indicating that the dissociation of N₂ is facile under mild conditions. The reaction process can be divided into two processes: N₂ transfer (TRF) and N-N dissociation (DIS). For V-series clusters, which has a V atom on Site 1, the rate-determining step is DIS, while for Ta-series clusters with a Ta on Site 1, TRF may be the rate-determining step or has energy barriers similar to those of DIS. The overall energy barriers for heteronuclear V₂ TaC₄ ^- and VTa₂ C₄ ^- clusters are lower than those for homonuclear V₃ C₄ ^- and Ta₃ C₄ ^- , showing that the doping effect is beneficial for the activation and dissociation of N₂ . In particular, V-Ta₂ C₄ ^- has low energy barriers in both TRF and DIS, and it has the highest N₂ adsorption energy and a high reaction heat release. Therefore, a trimetallic heteronuclear V-series cluster, V-Ta₂ C₄ ^- , is suggested to have high reactivity to N₂ activation, and may serve as a prototype for designing related catalysts at a molecular level.

Non-Dissociative Activation of Chemisorbed Dinitrogen on One or Two Vanadium Atoms Supported by a Mo6 S8 Cluster

Adsorption of N₂ on Mo₆ S₈ ^q _V_x clusters (x=0, 1, 2; q=0, ±1) were systematically studied by density functional theory calculations with dispersion corrections. It was found that the N₂ can be chemisorbed and undergo non-dissociative activation on single or double metal atoms. The adsorption and activation are influenced by metal types (V or Mo), N₂ coordination modes and charge states of the clusters. Particularly, anionic Mo₆ S₈ ^- _V₂ clusters have remarkable ability to fix and activate N₂ . In Mo₆ S₈ ^- _V₂ , two V atoms prefer to adsorb on two adjacent S-Mo-S hollow sites, leading to the formation of a supported V…V unit. The N₂ is bridged side-on coordinated with these two V atoms with high adsorption energy and significant charge transfer. The bond order, bond length and vibration frequency of the adsorbed N₂ are close to those of a N-N single bond.