Mechanistic Insights of Molecular Metal Polyselenides for Catalytic Hydrogen Generation

Molecular metal chalcogenides have attracted great attention as electrocatalysts for the hydrogen evolution reaction (HER). However, efficient utilisation of the active sites and catalytic performance modulation has been challenging. Here we explore the design of immobilized molecular molybdenum polyselenides [Mo 2 O 2 S 2 (Se 2 )(Se x )] 2– that exhibit efficient hydrogen evolution at low overpotential and stability over 1000 cycles. Density functional calculations provide evidence of a unimolecular mechanism in the HER process via the exploration of viable reaction pathways. The discussed findings are of a broad interest in the development of efficient molecular electrocatalytic materials.


Molecular metal chalcogenides have attracted great attention
as electrocatalysts for the hydrogen evolution reaction (HER).
However, efficient utilisation of the active sites and catalytic performance modulation has been challenging. Here we explore the design of immobilized molecular molybdenum polyselenides [Mo 2 O 2 S 2 (Se 2 )(Se x )] 2À that exhibit efficient hydrogen evolution at low overpotential and stability over 1000 cycles. Density functional calculations provide evidence of a unimolecular mechanism in the HER process via the exploration of viable reaction pathways.
The discussed findings are of a broad interest in the development of efficient molecular electrocatalytic materials.
Development of catalysts for the hydrogen evolution reaction (HER) is one of the key challenges in the realisation of a hydrogen economy based on electrocatalytic water splitting. The best performing HER catalyst is platinum, however its high cost and low natural abundance makes the development of noble-metal free catalysts an important area of research. 1 A variety of transition metal chalcogenides (ME x ) (M = Co, Fe, Ni, V, Nb, Ta, Mo, W; E = S, Se, Te) have been explored [2][3][4][5][6][7][8][9] in a quest for discovery of low production cost catalysts with high catalytic activity, and stability under aqueous corrosive conditions. Despite the extensive research conducted in this area, even the well explored chalcogenide, MoS 2 9-11 significantly underperforms as electrocatalyst for HER compared to Pt. This is due to the fact that most of the activity in MoS 2 comes from the edge sites which represent only a small fraction of the overall number of active sites in solid state materials like MoS 2 . Interestingly, molecular oxochalcogenides can incorporate by design a significant number of functional groups (which mimic very effectively the edge of solid-state chalcogenides) and have been shown to be at least as stable and effective as solid-state analogues. For example, the [Mo V 2 (S 2 ) 6 ] 2À cluster contains three disulphides per molybdenum, 12 while the hybrid family of LMo VI O(S 2 ) 2 (where L = bidentate ligand) type catalysts exhibited an indicative cooperative electronic effect between the inorganic and organic components. [13][14][15] When immobilised on an appropriate surface, they can effectively act as heterogeneous catalysts. Additionally, the design of the catalyst's key components (e.g. metal centres, type and number of chalcogens) can be pre-determined using a precise molecular chemistry approach. The all-inorganic molybdenum-and tungsten-based sulphide catalysts has shown overpotentials as low as À114 AE 3 mV at 10 mA cm À2 . 16 The presence of well-defined molecular structures rendered the computational modelling more informative and reinforced our understanding of the HER mechanism. Theoretical calculations demonstrated that the reduction occurs on the disulphide ligand, while the terminal oxo ligand stabilise the negative charge developed during the catalytic process. Interestingly, investigation of other molecular chalcoxide equivalents [17][18][19][20] for the HER and exploration of the edge based catalytic sites has not previously been attempted.
In this work, we investigated the effects of stoichiometric and structural control resulting from the appropriate combination of chalcogen elements and their ratio (e.g. O 2À /Se 2 2À ) within the coordination sphere of earth abundant transition metals and compare their mechanistic differences with their sulfur analogues. The dimeric nanocluster, [Mo 2 O 2 (m-Se) 2 (Se 2 )(Se 2 )] 2À 1, was synthesised following a facile method under mild conditions. The synthetic approach exploits the selenophilicity of molybdenum by exposing an alkaline molybdate solution to selenizing conditions. Alternatively, at elevated temperature and pressure regimes, the equilibrium shifts towards the formation of a mixture of species with the general formulae [Mo 2 O 2 (m-Se) 2 (Se 2 )(Se x )] 2À (where x = 2-4) 2.
During the reaction the species [Mo 2 O 2 (m-Se) 2 (Se 2 )(Se 2 )] 2À , [Mo 2 O 2 (m-Se) 2 (Se 2 )(Se 3 )] 2À and [Mo 2 O 2 (m-Se) 2 (Se 2 )(Se 4 )] 2À , seem to be in an equilibrium under harsh experimental conditions and found to co-crystallise within the same unit cell. Each method yielded orange-red crystals which were characterized by single crystal and powder X-ray structural analysis, FT-IR, 77 Se NMR and high-resolution electrospray ionization mass spectroscopy (ESI-MS). Subsequently, the [Mo 2 O 2 (m-Se) 2 (Se 2 ) 2 ] 2À 1 species was investigated as potential electrocatalyst to produce H 2 gas from acidic media where was found to exhibit excellent activity towards the hydrogen evolution reaction (HER). The oxoselenometalate compounds were synthesized under nitrogen atmosphere, following an optimised facile two-step process (see ESI †) utilizing low-cost materials, thus illustrating their potential for economically viable scalability.
The facile two-step synthetic procedure yields orange-red crystals and their single crystal X-ray diffraction (XRD) analysis shows that [(CH 3 ) 4 N] 2 [Mo 2 O 2 (m-Se) 2 (Se 2 ) 2 ] 1 crystallises in the orthorhombic system space group Pca2 1 with unit cell dimensions of a = 18.614(10), b = 11.422(6), c = 11.422(6) Å. The powder X-ray diffraction study conducted on the freshly synthesised material 1 (Fig. S5, ESI †) further confirm the purity of the synthesised material in the hexa-selenide form. Alternatively, the more challenging high temperature and pressure synthesis gives crystals containing a mixture of seleno-containing species with the formula that is by a diselenide (Se 2 2À ) and a tetraselenide (Se 4 2À ) group with average Se-Se bond distance of ca. 2.333 Å. The overall charge of the compound is compensated by two tetramethylammonium, (CH 3 ) 4 N + , cations per dimeric cluster.
In an effort to further characterize the isolated compounds, we employed high resolution electrospray ionization mass spectrometry (ESI-MS) [21][22][23] to determine the composition of the oxoselenometalate species in solution. The ESI-MS studies were performed in CH 3 CN. A series of singly charged distribution envelopes were observed, all in the region of 500-1000 m/z (Fig. S4, ESI †). More specifically, the distribution envelopes centred at 699.18 m/z value can be assigned to the intact compound [Mo 2 O 2 (m-Se) 2 (Se 2 ) 2 ] 2À . Additionally, 77 Se NMR spectrum of 1 recorded in DMF. The spectrum shows two distinct peaks corresponding to the two unique environments in the dimer located at shifts of 134 and 711 ppm relative to Me 2 Se (Fig. S6, ESI †).
To investigate the HER activity of the molecular catalyst we immobilised the catalyst on a glassy carbon electrode (GCE) support (see ESI †).
The amount of catalyst on the electrode surface was controlled by adjusting the loading concentration and volume. The performance of the catalyst improves as a function of the catalyst loading up to the value of 3.38 mmol cm À2 . Beyond this value the quality and mechanical stability of the produced electrocatalytic film degrades which leads either to similar or inferior performance due to detachment of the film from the electrode's surface. The results reported in Fig. 2 indicate that the cluster is highly active towards HER revealing a very small onset (40 mV) and low overpotential of À165 AE 3 mV (at j = 10 mA cm À2 ) as a function of the catalyst loading. Consequently, further increase of the applied overpotential results in a rapid rise of the cathodic current, which is  indicative of HER activity, in marked contrast to the bare GCE with carbon powder.
An important prerequisite of functional components used for large scale applications e.g. electrochemical water-splitting devices is the durability of the materials used. To assess this, we employed potential cycling as a simulation of the start-up/ shutdown conditions found in an electrolyser to provide us with information on the catalysts' long-term stability. Using catalyst 1 (3.38 mmol cm À2 ), we applied an accelerated scan rate of 100 mV s À1 between a potential range of 0.10 to À0.25 V (vs RHE) for 1000 cycles, Fig. 2b. At the end of the cycling, the catalyst showed negligible loss in cathodic current density, reaching an impressive density of j 4 100 mA cm À2 below 250 mV. These results attest to the electrochemical stability of [Mo 2 O 2 (m-Se) 2 (Se 2 ) 2 ] 2À 1 and its potential use in electrochemical HER devices. Finally, to confirm that the reaction taking place on the working electrode corresponds to the hydrogen evolution reaction (HER), the evolved gas was analysed using gas chromatography (Fig. S8, ESI †). As such, the faradaic efficiency of the molecular catalyst in 1 M H 2 SO 4 was determined to be 102.3% AE 3.0 with a TOF of 0.019 s À1 at À190 mV demonstrating the evolution of hydrogen at full faradaic efficiency.
To investigate the underlying processes and reaction pathways, the electronic structure of these molecules and their reduced analogues were studied on the basis of density functional theory (DFT). The geometries of both anions were optimised starting from the single crystal X-ray diffraction data (see ESI †). The electronic structure of these anions is similar and according to orbital analysis and charge distribution, both ligands -bidentate k 1,4 -Se 4 2À and side on Z 2 -Se 2 2À -formally contribute À2 charges each. Mo centres present +5 formal oxidation state with a Mo-Mo bond (2.855 Å) as corroborated by the bond valence sum calculations analyses and the crystal structure. Each perselenide moiety has a doubly filled p* 2 p* 2 group molecular orbital (MO) pair. Examination of the HOMO-LUMO gap provides a preliminary good estimation of the chemical potential of the species and indication of orbitals populated by the incoming electrons. 0 e À and 1 e À reduction of catalytic species are depicted in Fig. S9 (ESI †). In case of 0 e À , the energy gap is higher in [Mo 2 O 2 -(m-Se) 2 (Se 2 ) 2 ] 2À than [Mo 2 O 2 (m-Se) 2 (Se 2 )(Se 4 )] 2À and follows the same trend as sulphur analogues of these compounds. However, the HOMO-LUMO gap and LUMO energy are smaller in selenium analogues compared to their respective sulphur anions ( Fig. S9 and S10, ESI †). The first unoccupied MOs involve d metallic orbitals and a s-antibonding Se-Se orbital. Interestingly, in case of 1 e À reduced species, the energy needed to insert the second electron is higher for Se compounds than S analogues ( Fig. S9 and S10, ESI †). The half occupied orbital present Sigma Se-Se antibonding character. In order to investigate the electrochemical steps that are taking place it is important to identify the different pathways following the two-fold reduction of the molecular catalysts leading to protonation equilibria. Therefore, the stability of all possible isomeric forms was evaluated. Even though those unprotonated and short-lived species may not be detected and characterised in solution, evaluating their relevant stability is relevant as it can serve as a reference state for subsequent steps.
Considering the symmetric structure of the [Mo 2 O 2 (m-Se) 2 (Se 2 ) 2 ] 2À complex, we explored a series of isomers that can be formed during the two-variable proton/electron addition leading to the global panorama of possible HER equilibria (Fig. 3). The cleavage of perselenide bond upon 1 e À and 2 e À reduction was observed. It is expected based on molecular orbital analysis, because these anions feature s-antibonding Se-Se orbitals which might in principle be occupied by incoming electrons inducing the breaking of Se-Se bond. This behaviour was also observed in the sulphur analogues of these anions. The complete nomenclature of the obtained isomers and proton affinities are listed in the Fig. S12-S19 (ESI †). For compound 1, the most stable doubly reduced and single-protonated species is B1 since it stabilises the mitigation of charge developed in the newly formed diselenol species. Consequently, the second protonation step which derives from B1 is BB1, which exhibits a free energy value of À40.7 kcal mol À1 (see Fig. S12, ESI †). Taking into consideration the proximity of the neighbouring hydrogen atoms in this structure, that raised the possibility of the existence of a Tafel-type re-combination step. Further examination of this step for both anions revealed activation energies of 62 kcal mol À1 from double-protonated BB1 ( Fig. S21 and S26, ESI †) which are too high to be operative at room temperature. To this extend, departing from BA1 whereas one selenium and one oxygen atom are protonated, Taffel reaction activation energy is more favourable with an activation energy of 37.45 kcal mol À1 (Fig. S22, ESI †). In a similar manner, we aimed to examine the Heyrovsky step departing from single-and double-protonated [H 2 Mo 2 O 2 (m-Se) 2 (Se 2 )(Se x )] 2À (x = 2, 4) structures. This step is thermodynamically favourable across the range of molecular systems with calculated released energies within À46.12 to À4.51 kcal mol À1 (Fig. 4 and Fig. S20, S23-S25, ESI †). The activation energy for this step is in range of +2.54 to +18.62 kcal mol À1 for both anions and all reactant complexes taken into consideration. These E a values indicate that activation can happen at room temperature conditions. For compound 1, BB1 intermediate [H 2 Mo 2 O 2 (m-Se) 2 (Se 2 ) 2 ] 2À , the obtained transition state has a relative free energy value DG ‡ = +18.62 kcal mol À1 (Fig. S20, ESI †) which is comparable to the activation energy value obtained for compound 2 when the protons interact with the side on Z 2 -Se 2 2À terminal group (BB1, Fig. S24 (Fig. S25, ESI †). It is worth noting that single-protonated forms have significantly lower activation energies than the double-protonated ones. The activation barriers for compounds 1 (B1) and 2 are only 2.54 ( Fig. 4) and 5.18 kcal mol À1 (Fig. S23, ESI †), respectively. Thus, the model of the transition state for single protonated forms is the most probable candidate. For the non-reduced, single-and double-reduced species, their double-protonated forms are the most thermodynamically stable. Considering electrochemical routes, 2e À addition to double-protonated non-reduced species is non-favourable. In contrast, reduction reactions from single protonated and non-protonated forms are more feasible (Fig. 3). Thus, our results clearly point towards the Heyrovsky mechanism based on thermodynamic favourability and the high activation energies observed for the alternative Taffel step.
In conclusion, a more accessible and high-yield synthesis for the selenium analogue to the well-known [Mo 2 O 2 S 2 (S 2 ) 2 ] 2À structure makes it possible to explore the electronic and chemical properties of the molecule, which have previously remained unexplored. The performance of the molecular catalyst for the hydrogen evolution reaction was evaluated while we propose viable mechanistic pathways in relation to the interaction between catalysts and substrate while we identified intriguing differences with their sulfur analogues. The initial 1 e À reduction is easier compared to the sulphur containing analogue however the second reduction is more energetically demanding which is reflected to the observed slightly elevated Tafel slope. This indicates the significant electronic changes of the selenium containing structure relative to the sulphur. Notably, these results represent an attractive strategy for preparing next generation HER electrocatalysts. Understanding the electronic underpinnings of trends in activity as a function of the metal or chalcogen incorporated in the molecular species, can help inform the design of new catalysts and exploitation of cooperative and structural effects. Moreover, the development of facile low-cost synthesis of a molecular catalyst is highly desirable for the easier scale up in large scale applications.

Conflicts of interest
There are no conflicts to declare.