Vishakha Kaim1
Meenakshi Joshi
Matthias Stein
Sandeep Kaur-Ghumaan- 1Department of Chemistry, University of Delhi, Delhi, India
- 2Max-Planck-Institute for Dynamics of Complex Technical Systems, Molecular Simulations and Design Group, Magdeburg, Germany
Molecular hydrogen (H2) is one of the pillars of future non-fossil energy supply. In the quest for alternative, non-precious metal catalysts for hydrogen generation to replace platinum, biological systems such as the enzyme hydrogenase serve as a blueprint. By taking inspiration from the bio-system, mostly nickel- or iron-based catalysts were explored so far. Manganese is a known oxygen-reducing catalyst but has received much less attention for its ability to reduce protons in acidic media. Here, the synthesis, characterization, and reaction mechanisms of a series of four mono-nuclear Mn(I) complexes in terms of their catalytic performance are reported. The effect of the variation of equatorial and axial ligands in their first and second coordination spheres was assessed pertaining to their control of the turnover frequencies and overpotentials. All four complexes show reactivity and reduce protons in acidic media to release molecular hydrogen H2. Quantum chemical studies were able to assign and interpret spectral characterizations from UV–Vis and electrochemistry and rationalize the reaction mechanism. Two feasible reaction mechanisms of electrochemical (E) and protonation (C) steps were compared. Quantum chemical studies can assign peaks in the cyclic voltammetry to structural changes of the complex during the reaction. The first one-electron reduction is essential to generate an open ligand-based site for protonation. The distorted octahedral Mn complexes possess an inverted second one-electron redox potential which is a pre-requisite for a swift and facile release of molecular hydrogen. This series on manganese catalysts extends the range of elements of the periodic table which are able to catalyze the hydrogen evolution reaction and will be explored further.
Introduction
Molecular hydrogen (H2) is one of the most-promising energy carriers in the future due to its zero-carbon combustion (Camack et al., 2001; Penner, 2005). The design and development of novel molecular catalytic systems for the hydrogen evolution reaction (HER) have received significant attention during the last decade. The electrocatalytic reduction of a proton to molecular hydrogen (H2) might help in the industrial process development, leaving behind fossil fuels, and toward a hydrogen economy. (Schlögl, 2013) Most current industrial processes rely on platinum as a highly efficient catalyst for hydrogen generation. However, there is an urgent need to substitute precious platinum with more abundant first-row transition metal catalysts that take inspiration from biological systems (Artero and Fontecave, 2005; Cracknell et al., 2008).
In global industrial processes, hydrogen is used on a massive scale, for example, in ammonia production (Haber–Bosch) or as a highly selective reducing agent in the production of active pharmaceutical ingredients (APIs), to name only a few. So far, hydrogen production is based on fossil fuels and/or highly efficient (precious) metal catalysts (Lubitz and Tumas, 2007). The processes are relatively energy-inefficient and lead to emission of significant quantities of greenhouse gases. Evolution has led to a family of enzymes, namely, the hydrogenases that use hydrogen as an energy carrier for their metabolism. Those hydrogenases have developed their own bio-inorganic transition metal catalysts to produce or utilize hydrogen as an energy source. They do so at ambient temperature and normal pressure and make use of the non-precious metals iron and nickel. Their active sites have inspired the synthesis of many bio-mimetic and bio-inspired catalysts (Evans and Pickett, 2003).
A direct comparison of a platinum catalyst with a hydrogenase equivalent has shown that the bio-hydrogen production by the enzyme occurs at rates comparable to those of electrodeposited platinum and with less susceptibility to CO poisoning (Jones et al., 2002). Since the density of the catalytic sites for the enzyme is a lot lower than that for platinum, the enzyme’s active site must turn over faster than a corresponding platinum atom to attain the same overall current. The performance of an electrode-immobilized nickel catalyst was directly compared to the performance of the [NiFe]-hydrogenase enzyme. At acidic pH and in the presence of carbon monoxide, the enzyme outperformed the synthetic catalyst and pointed toward the utilization of bio-inspired complexes in fuel cells (Rodriguez-Maciá et al., 2015).
Drawing inspiration from the naturally occurring [FeFe]-, [NiFe]-, and hydrogenase enzymes, a plethora of complexes have been reported, for reviews see (Evans and Pickett, 2003; Bruschi et al., 2005; Kaur-Ghumaan and Stein, 2014; Simmons et al., 2014; Onoda and Hayashi, 2015; Schilter et al., 2016), most of which are homo- (FeFe) or hetero-bi-nuclear (NiFe, NiRu, and NiMn). The binuclear systems reported as electrocatalysts for the HER use mainly thiolates (with or without a pendant amine) and various electron-rich ligands (e.g., cyanides, phosphines, carbenes, π-accepting, and other non-innocent ligands) based on the structure known for the hydrogenase enzyme active sites (see Scheme 1). The binuclear systems reported hardly showed any activity under the conditions at which the enzyme functions and have limitations with respect to catalyst stability. The availability of an electron reservoir, balanced proton affinity, and a fine-tuning of the electronic and structural properties of a hydride were recently identified to be the most critical design aspects for bio-inspired HER catalysts (Sun et al., 2022).
SCHEME 1. Active site of the [FeFe] hydrogenase enzyme and bio-inspired mono-nuclear hydrogen-evolving complexes based on iron (Kaur-Ghumaan et al., 2010), nickel (Wilson et al., 2006; Brazzolotto et al., 2016), cobalt (Fang et al., 2014), or manganese (Kaim and Kaur-Ghumaan, 2021; Yang et al., 2021) and this work.
Mono-nuclear catalysts, however, are less explored but have the advantage of an easier synthesis and a higher density of catalytic sites (see above). Those reported contain, among others, Fe (Kaur-Ghumaan et al., 2010; Natarajan et al., 2017), Co (Wilson et al., 2006; Dempsey et al., 2009; Artero and Saveant, 2014; Fang et al., 2014), and Ni (Wilson et al., 2006; Brazzolotto et al., 2016) (see Scheme 1).
Some of the coordinating ligands have been implicated to be directly or indirectly involved in the HER as electron-donors (Bigi et al., 2010), be partially redox-active (Straistari et al., 2017), or function as proton relays (DuBois, 2014).
Manganese is an exquisite catalyst for oxygen reduction/evolution (Hirahara et al., 2014; Wiechen et al., 2014). Like iron, manganese is easily available and abundant and shows a benign environmental profile compared to platinum. Mn-complexes reported until now have been mostly studied as catalysts for electro- or photo-catalytic organic oxidations/reductions (C–H hydroxylation; P450-type activity), alkene epoxidation, and CO2 reduction in homogenous or heterogeneous phases; for reviews see (Fernandez et al., 2022; Sorais, 2022). The complexes designed for these catalytic transformations use porphyrin and corrole derivatives or strong-chelating non-heme pyridylamino ligands, π-acceptors (CO), and pyridine-, diimine-, or N-heterocyclic carbene-based ligands.
However, the performance of manganese complexes in catalyzing HER has not received much attention until now, despite its abundance in the earth’s crust. The isoelectronic properties of non-native MnI with the bio-inorganic FeII suggest that MnI-based complexes may also be suitable HER catalysts. However, there are very few Mn-based HER catalytic systems reported in the literature (Kaim et al., 2019; Kaim and Kaur-Ghumaan, 2019; Pan et al., 2019; Pan and Hu, 2020; Kaim and Kaur-Ghumaan, 2021; Yang et al., 2021) (see Scheme 1).
This prompted the design of a series of novel MnI-based catalysts as an attractive alternative for hydrogen production. We, here, report the synthesis of a series of distorted octahedral mono-nuclear manganese (I) complexes which all show catalytic activity toward hydrogen production in an acidic solution. Their electrocatalytic proton reduction behavior and HER using TFA as a proton source are demonstrated. Carbon monoxide ligands are also found in the biological system (see above), phosphine ligands have electronic effects very similar to those of small inorganic CO and CN− ligands, and sulfur coordination also occurs in the enzyme.
Upon variation of axial (PPh2(PhCOOH-p)) 1 and 3 vs. (PPh3) 2 and 4 and equatorial ligands containing heterocycles with only nitrogen 3 and 4 vs. one sulfur atoms 1 and 2, their individual effects on the electronic structure and catalytic performance can be investigated. The electrochemical features of these complexes in the absence and presence of TFA are analyzed using cyclic voltammetry.
Redox transitions and catalytic mechanisms can be rationalized by quantum chemical calculations. The calculations show that complexes 1–4 display an inverted redox potential, whereas the nature of the axial ligand has almost no effect on catalytic performance, and the bidentate equatorial heterocycle is partially redox-active and serves as an initial site of protonation, following a first one-electron reduction event. An intramolecular proton transfer toward the central metal ion occurs after the second reduction event. Finally, an end-on H2-Mn(I) complex is formed from which the release of molecular hydrogen is facile and thermodynamically feasible.
Materials and methods
Four mono-nuclear Mn-complexes were synthesized and structurally characterized using FTIR, 1H, and 31P NMR and MS spectrometry. For one complex, single crystals could be obtained with sufficient diffractions in order to be able to determine the crystal structure. The redox behavior of the four catalysts was investigated by cyclic voltammetry in the absence of acids to characterize the intrinsic redox processes for each catalyst. Finally, the catalytic behavior of the Mn-complexes toward the hydrogen evolution reaction (HER) was investigated in the presence of trifluoroacetic acid (TFA) in acetonitrile. The structural parameters, redox potentials, and catalytic behavior of all the four mono-nuclear Mn-complexes were obtained from quantum chemical calculations (using Density Functional Theory, DFT). These results are compared with the corresponding experimentally observed results in each section. The experimental and computational details are discussed in detail in the electronic supporting information (ESI).
Synthesis and characterization of mononuclear Mn complexes
Mononuclear Mn complexes fac-[(Mn(CO)3(κ2-S2NC7H4)(PPh2(PhCOOH-p))] 1, fac-[(Mn(CO)3(κ2-S2NC7H4)(PPh3)] 2 and fac-[Mn(CO)3(κ2-SN2C7H5)(PPh2(PhCOOH-p)] 3 were synthesized by the reaction of [(Mn(CO)3(µ-S2NC7H4))2] A and [(Mn(CO)3(µ-SN2C7H5)2)] B complexes with triphenylphosphine (TL1, PPh3) and 4-(diphenylphosphino)benzoic acid (TL2, PPh2(PhCOOH-p)) ligands in dichloromethane solvent for 72 h at room temperature as shown in Scheme 2. For comparison, the Mn analog fac-[(Mn(CO)3(κ2-SN2C7H5)(PPh3)] 4 was also synthesized according to the published protocol (Scheme 2) (Abedin et al., 2018). All the complexes 1–4 were characterized using different spectroscopic techniques such as FTIR, 1H, and 31P NMR spectroscopy and MS spectrometry (see Supplementary Figures S1–S7, ESI).
SCHEME 2. Scheme for synthesis of complexes 1–4.
Computational details
Most computational details are given in the ESI. The GGA BP86 functional (Perdew, 1986; Becke, 1988) and the hybrid B3LYP functional (Vosko et al., 1980; Lee et al., 1988; Becke, 1993; Stephens et al., 1994) with all-electron def2-TZVP basis sets (Weigend and Ahlrichs, 2005) and dispersion corrections (Grimme et al., 2010) were used. All structures were characterized as minima by the absence of imaginary frequencies. In order to account for solvation effects, the implicit COSMO solvation model was used (Klamt and Schüürmann, 1993; Schäfer et al., 2000).
Results and discussions
Crystallization and structural characterization of Mn-complex 1
For complex 1, single-crystals could be obtained, which were analyzed by single-crystal X-ray diffraction. The crystallographic parameters and selected bond lengths and bond angles of complex 1 are reported in Supplementary Tables S1 and S2 in the supporting information (see ESI). Since the coordination sphere of complexes 2–4 is very similar, no significant structural differences are to be expected. X-ray structure analysis confirms the equatorial chelation by 2-mercaptobenzothiazole through N1 and S1; axial 4-(diphenylphosphino)benzoic acid and three carbonyl groups completed a distorted octahedral geometry (see Figure 1A). The N and S atoms of the chelating ligand and two carbonyl groups occupy the equatorial positions. In the crystal, 1 dimerizes by forming two intermolecular hydrogen bonds between the benzoic acid entities (Figure 1B). (Braga et al., 2001)
FIGURE 1. (A) X‒ray crystal structure for fac-[(Mn(CO)3(κ2-S2NC7H4)(PPh2(PhCOOH-p))] 1. Hydrogen atoms have been omitted for clarity. (B) Hydrogen bonding in the crystal structure of 1.
Computational characterization of mono-nuclear Mn-complexes
Complexes 1–4 were structurally optimized using dispersion-corrected density functional theory (DFT). In the following section, only BP86 structural parameters are discussed throughout the manuscript, while B3LYP results are given in the ESI.
The structure of complex 2 is very similar to that of complex 1, except for the axial PPh3 ligand that replaces the PPh2(PhCOOH-p) ligand. Similarly, in complexes 3 and 4, Mn is coordinated with five ligands in a distorted octahedral geometry. The arrangement of axial PPh2Y and CO ligands in complexes 3 and 4 is the same as that in 1 and 2, respectively, except for the equatorial chelating ligand. Here, a bidentate 2-mercaptobenzimidazole ligand is present in an equatorial position in complexes 3 and 4 (see Figure 2). The first coordination sphere is identical in all four complexes, as shown in Figure 2.
FIGURE 2. Structures of complexes 1–4 with identical first coordination spheres and variation of atomic composition in the second shell.
Neither equatorial Mn‒N bond distances (2.06, 2.06, 2.06, and 2.06 Å, in 1–4) or Mn‒S bond distances (2.49, 2.49, 2.52, and 2.52 Å, in 1–4), nor axial Mn‒P bond distances (2.37, 2.36, 2.36, and 2.36 Å, in 1–4) are significantly influenced upon modifications of the equatorial and axial ligands. The mild distortion from octahedral symmetry can be seen in the angle < (COax-Mn-PR3) for 1–4 of 175°, 176°, 175°, and 175°; and the equatorial < (N‒Mn‒S) angle slightly deviating (<N‒Mn‒S = 64°-68°) from 90°. They, indeed, show that there are no large structural differences between complexes 1–4 (Supplementary Tables S3–S6).
Characterization of redox activity
Mn-complex reduction in the absence of acid
The intrinsic redox behavior of complexes 1–4 was investigated by performing cyclic voltammetric (CV) experiments in CH3CN using 0.1 M [N(n-Bu4)][PF6] as the supporting electrolyte at a scan rate of 0.1 Vs‒1. The CVs of the complexes for oxidation are shown in Supplementary Figure S8 (see ESI). Scan rate linear dependence studies show that all the redox processes are diffusion-controlled (see Supplementary Figure S8, ESI).
Complexes 1 and 3 (with PPh2(PhCOOH-p) ligand) each show two irreversible reduction peaks at
FIGURE 3. Cyclic voltammograms for reduction of complexes 1–4 in CH3CN at a scan rate of 0.1 V s−1. Assignment of redox transitions can be done by quantum chemical calculations (see text for details).
Complexes 2 and 4 (with PPh3 ligand) also exhibited two irreversible reduction peaks at
The structural interpretation of CV is difficult and can only be indirect. Quantum chemically calculated redox potentials allow the assignment of redox transitions to structural changes when calculated potentials and experiments agree to within a ∼100 mV. In order to assign redox CV peak positions, redox potentials for the sequential two-one electron reduction steps (Eqs 1, 2) plus one two-electron reduction step were calculated (Eq. (3)) (Supplementary Table S7, ESI).
The reduction potentials of complexes 1–4 are very close and occur at ‒1.63, ‒1.67, ‒1.62, and ‒1.65 V for the first and at ‒2.16, ‒2.18, ‒2.20, and ‒2.19 V for the second one-electron reduction, respectively. This shows that the electronic structure of all four complexes is very similar.
The two-electron reduction for complexes 1–4, however, is suggested to occur at potentials of ‒1.89, ‒1.92, ‒1.91, and ‒1.92 V and are thus at less negative potentials than the second one-electron reduction process (inverted redox potential).
Experimentally, the first one-electron reduction near ‒1.6 V for complexes 1–4 could not be resolved, which was also observed for mono-nuclear iron complexes (Natarajan et al., 2017).
The DFT calculated two-electron reduction potentials, as well as the second one-electron reduction potentials, are in excellent agreement with the experimentally observed peaks at ‒1.83, ‒1.93, ‒1.99, and ‒2.06 V and ‒2.19, ‒ 2.07, ‒2.30, and ‒2.36 V, respectively. This allows a definite assignment of observed reduction peaks in the CVs of complexes 1–4 corresponding to the formation of doubly reduced species (12‒ - 42‒), whereas the second CV peak is the transition 1-/12- to 4-/42-.
It can be suggested that the mono-reduced species 1-/2-/3-/4- are too short-lived to be resolved in the CV. Given the facile two-electron reduction process, it is likely that the intermediates 12-/22-/32-/42- are the first stable intermediates. This agrees with the suggested reaction mechanism (see below) and an EECC sequence of electrochemical (E) and chemical (C) events.
Computed structural changes during redox transitions
Quantum chemical calculations are able to elucidate structural changes in the electronic structure of transition metal complexes during redox transitions. Such detailed information is difficult to obtain experimentally. Upon the first one-electron reduction, the Mn−S bond in complexes 1‒–4‒ dissociates and opens a free coordination site at the Mn center and the ligand sulfur atoms. This breaking of the Mn−S bond can explain the irreversible character of the peaks in the CVs. The calculated Mn‒S distances are elongated (by 0.7–1.0 Å) to 3.24, 3.30, 3.44, and 3.51 Å in 1‒, 2‒, 3‒, and 4‒, respectively.
No further structural changes are observed in the reduced complexes 1‒ – 4‒. The calculated unpaired spin density is localized mainly on the Mn-atom (Figure 4 and Supplementary Figure S9) in the reduced species, which confirms that the first one-electron reduction event is primarily metal-centered. The close redox potentials for 1–4 show that minor changes in the second coordination sphere do not affect the reduction potential. From an analysis of the molecular orbitals HOMO and SOMO, it can be seen that the first one-electron reduction in 1–4 is mainly Mn-based and only partially involves the S-atom of the chelating ligand (Figure 4 and Supplementary Figure S9).
FIGURE 4. Display of (A) unpaired spin density distribution (at an isocontour value = 0.01e) and (B) highest occupied molecular orbitals of complexes 1 and 2 (at an isocontour value = 0.02e). The PhCOOH-p group in 1 affords a partial ligand-based second reduction event in 1-/12- whereas 2-/22- is metal-centered.
Addition of a second electron leads to the doubly reduced species 12‒, 22‒, 32‒, and 42‒ which are characterized by a further increase in Mn‒S bond lengths (by 0.2–0.4 Å) to 3.55, 3.66, 3.67, and 3.79 Å. Changes of other Mn‒ligand bond distances are negligible (see Supplementary Tables S3–S6, ESI). The structures of 12−–42− are distorted due to re-orientations of the axial ligands. The < (P‒Mn‒COaxial) angle deviates significantly from the ideal 180° (148°, 136°, 146°, and 127° in 12−–42− complexes). The second one-electron reduction is facile but less than the first one-electron step (by 49–57 kJ/mol).
The reduction 1-/12- and 3-/32- is partially ligand-based due to the delocalization of electrons over the PhCOOH-p group, giving stable resonating structures of the reduced ligand. For 2-/22- and 4-/42-, however, the second reduction is mainly metal-centered (Figure 4). The second reduction does not have any significant effect on the reduction potential of the four complexes.
Catalytic activity of mono-nuclear manganese complexes
The performance of complexes 1–4 as HER catalysts in the presence of trifluoroacetic acid (TFA) was explored in acetonitrile (pKa of TFA in CH3CN = 12.7, E0 in CH3CN = ‒ 0.89 V). Addition of TFA to 1-4 results in two irreversible waves with increasing peak intensity upon further addition of acid, thus indicating catalytic electrochemical proton reduction (Figure 5; Table 1 and Supplementary Figure S10, ESI). In the absence of a catalyst, negligible potential and currents were observed between ‒1.5 and ‒2.6 V for proton reduction in the presence of TFA in CH3CN at 0.1 Vs‒1 (Supplementary Figure S11, see ESI).
FIGURE 5. CVs for complex 1 (1 mM) in the absence (top curve, dashed black line) and presence of TFA (5–46 mM) at 0.1 Vs‒1.
TABLE 1. Electrochemical data for complexes 1–4 in acetonitrile.
Initial addition of acid shows reductions at or slightly more negative reduction potentials than the complex peaks and/or followed by the appearance of new peaks on the addition of ∼15–20 mM of acid. This suggests a possible reduction of the complexes followed by protonation with the addition of acid. The peaks shift toward more negative reduction potential with increasing acid concentration and level off at ∼189 mM (1), ∼123 mM (2), ∼100 mM (3), and ∼143 mM (4) of TFA. This shows that the increase in acid addition leads to a sustained reduction. The levelling off indicates that most of the reactive species (at the electrode/electrolyte interface) have already been oxidized or reduced.
The first peak can be assigned to be the catalytic wave with a half wave potential Ecat of ‒1.87 V (1), ‒1.88 (2), ‒1.88 (3), and ‒1.95 (4) (Figure 5; Table 1 and Supplementary Figure S10, ESI). The catalytic wave corresponding to the electrocatalytic reduction of protons to hydrogen was further confirmed by bulk electrolysis experiments in acetonitrile for each complex 1–4 (0.25 mM) between ‒2.0 to −2.25 V in CH3CN/0.1 M [N(n-Bu4)][PF6]/6 mM TFA (Supplementary Figure S12, ESI). From the controlled potential electrolysis (CPE) experiments for the complexes, TONs were calculated to be 0.19 (1), 0.44 (2), 0.58 (3), and 0.71 (4) (Fu et al., 2015a; b). Evaluation of CVs measured at different scan rates for different TFA concentrations indicates the dependence of icat on the scan rate (Supplementary Figure S13, ESI). From the linear plot of icat vs. [catalyst], it is evident that the reaction is first order in complex concentration for 1–4 at fixed acid concentrations (Supplementary Figure S14, ESI). The overpotential for complexes 1–4 was calculated using Evans’s method to be in the range of 0.96–1.06 V (using Ecat/2) (Felton et al., 2006; Helm et al., 2011).
The appearance of a second reduction peak in the CVs of 1–4 in the presence of TFA might originate from the formation of stable adducts between the acid (AH) and its conjugate base (A−). Homoconjugation complicates the determinations of the thermodynamics of H2 evolution since AH and AHA− may act as proton sources for hydrogen production. Homoconjugation in acetonitrile has been reported to occur at high concentrations of weak acids and influence the acid–base equilibrium above the association constant Kc. Trifluoroacetic acid (TFA) displays a high value of Kc = 103.9 and is therefore expected to undergo homoconjugation to a significant extent even at low concentrations. This process affects the pKa of the acid in CH3CN, leading to different E0 values for the acid (which can be calculated using Eq. 4), which can, thus, be challenging in determining accurate overpotentials (Appel and Helm, 2014).
Figure 6 shows the increase of catalytic current and TOF with acid concentration. Higher
FIGURE 6. (A)
The Mn complexes 1–4 showed activity with lower TOF/TON values than previously reported Mn complexes but improved
Mechanism of hydrogen evolution reaction
The electrochemical proton reduction can occur via ECEC or EECC (E = electrochemical and C = chemical step) mechanisms for all the four mononuclear Mn complexes. Changes in Gibbs free energies (ΔG) for each reaction step are reported in Schemes 3 and 4, while the results in the absence of solvent are reported in the ESI (Supplementary Tables S8–S11, ESI).
SCHEME 3. Proton reduction and hydrogen evolution mechanism of complexes 1 and 3.
SCHEME 4. Proton reduction and hydrogen evolution mechanism of complexes 2 and 4.
ECEC vs. EECC mechanistic pathways
After one-electron reduction (E), the Mn‒S bond dissociates to give 1‒ and 3‒ species (see above) and provides an accessible site for protonation at either Mn or S atoms (see Scheme 3).
In the ECEC pathway, protonation (C) at the sulfur atom is energetically favorable (ΔG of ‒1,148.1 and ‒1,153.6 kJ/mol) to form 1SH and 3SH species with S‒H bond distances of 1.45 and 1.44 Å, respectively. Direct protonation of Mn was not feasible, and the proton reverts to the sulfur again. Only after the second one-electron reduction step (E), an intramolecular proton transfer occurs spontaneously from the S atom to the Mn atom to form 1MnH− and 3MnH−. The Mn‒H bond lengths are 1.58 Å in 1MnH− and 3MnH−. The formation of 1MnH− and 3MnH− is energetically feasible with ΔG of ‒386.8 and ‒383.4 kJ/mol.
The second protonation step for 1MnH− and 3MnH− occurs at the manganese hydride to give end-on H2 complexes 1MnH2 and 3MnH2 (ΔG = ‒1,173.7 and ‒1,179.5 kJ/mol) with H‒H bond distances of 0.93 Å and 0.89 Å, respectively. The binding of H2 to manganese is end-on with different Mn‒H distances of 1.66, 2.05 Å in 1MnH2 and 1.68, 1.99 Å in 3MnH2. The second hydrogen atom only weakly interacts with the S (S-H = 1.76–1.99 Å).
The ECEC reaction mechanism is devoid of any significant structural changes, and the complex remains intact except for a Mn‒S bond cleavage in the first electron reduction step (see Supplementary Tables S3–S5). Finally, upon H2 release from 1MnH2 and 3MnH2, the Mn‒S bond re-forms to regenerate complexes 1 and 3. The release of molecular hydrogen is an exergonic process with ΔG values of ‒80.0 and ‒71.1 kJ/mol (Supplementary Scheme S1, see ESI).
In the EECC mechanism (Scheme 3), upon the addition of the second electron, the axial aromatic ligand (see above) significantly deviates from the ideal 180°, and the equatorial Mn‒S bond length increases by 0.2–0.3 Å in species 12‒ and 32‒. The remaining other Mn‒ligand bond distances remained almost unchanged. In 12‒ and 32‒ possible sites for protonation are available at both the Mn and S atoms, and the first protonation is found to occur at Mn to form 1MnH− and 3MnH− (a similar structure is observed in the ECEC mechanism, see above). The ΔG values of protonation to obtain 1MnH− and 3MnH− are ‒1,279.0 and ‒1,285.3 kJ/mol, respectively, and thus almost indistinguishable. The optimized structural parameters of all the reaction intermediates are listed in Supplementary Tables S3–S5 (ESI). The remaining steps of hydrogen generation and release in the EECC mechanism are the same (see Scheme 3) as discussed above for the ECEC mechanism.
Among the two possible hydrogen releasing mechanisms, the two-electron reduction of the catalyst followed by two protonation events, i.e., the EECC mechanism appears to be more plausible to occur in the presence of TFA (pKa of TFA in CH3CN = 12.7) since the mono-reduced species could not be observed in the CV. Furthermore, the di-anion more easily reduces the proton (by ∼131 kJ/mol) in the EECC mechanism than that of the mono-anion in the ECEC mechanism. Nevertheless, the use of a very strong acid (with pKa < 6–8) could favor the ECEC mechanism for the HER.
Complexes 2 and 4 follow the same EECC and ECEC mechanism as discussed above for complexes 1 and 3 and are thus not discussed in detail here. Here, we briefly discuss the EECC mechanism for complexes 2 and 4. The ECEC proton reduction mechanism and corresponding Gibbs free energies of complexes 2 and 4 are given in ESI, Supplementary Scheme S1. The optimized structural parameters of all the reaction intermediates are listed in Supplementary Tables S4–S6.
In the EECC mechanism (Scheme 4) of complexes 2 and 4, the protonation (C) occurs at the Mn atom of 22‒ and 42‒ (Mn‒H = 1.58 Å) to form 2MnH− and 4MnH− species. In the next step (C), the second proton binds to the terminal hydride to form 2MnH2 and 4MnH2 species with the H‒H distance of 0.97 Å and 0.95 Å, respectively. The Mn‒H distances are 1.65, 2.13 Å in 2MnH2 and 1.66, 2.14 Å in 4MnH2, which indicates an end-on metal-H2 complex. In the final step, the 2MnH2 and 4MnH2 complexes release molecular hydrogen with ΔG of ‒84.1 and ‒64.8 kJ/mol and regenerate the complexes 2 and 4.
It is worth mentioning that the presence of different ligands in complexes 1–4 causes only a minimal change in the structures as well as in the ΔG value (by 0.4–19 kJ/mol) of the electrochemical (E) and chemical (C) steps of the HER of these Mn-complexes.
Conclusion
Mononuclear Mn complexes 1–4 are able to catalyze the electrochemical proton reduction and the HER with TFA as a proton source. After careful characterization with NMR, FT-IR, UV–Vis spectroscopy, cyclic voltammetry, and quantum chemical calculations, we were able to assign and interpret the redox transitions and suggest mechanistic pathways for hydrogen generation.
The short-lived mono-reduced species were shown not to be part of the catalytic cycle. Only the two-electron reduced species (1–4)2‒ at −1.9 V and the second one-electron transition (1–4)1‒/(1–4)2‒ at −2.2 V can be observed. This implies a hydrogen evolution mechanism starting from the two-electron reduced intermediate along with two sequential protonation events at the Mn atom and then the formation of a loosely coordinated end-on H2 complex.
Variation of ligands in the second coordination sphere affects the redox potentials, reaction mechanism, and catalytic efficiency only to a minor degree.
The high cost of platinum catalysts used in hydrogen fuel cells is a limiting factor in the commercialization of fuel cell electric vehicles in addition to its negative environmental impact. The current global demand for platinum (approximately 200 t p. a.) is high compared to the available platinum reserves, and only about 17% of the used platinum is recycled (Duclos et al., 2017).
The quest for alternative, more cost-effective catalysts to maintain the efficiency of hydrogen fuel cells is ongoing.
With the introduction of manganese as an additional earth-abundant metal atom into the series of mono-nuclear hydrogen generating catalysts (in addition to Fe, Co, Ni, Ru), the synthetically accessible chemical space opens up. Its low cost and abundant availability, plus a benign environmental profile, make it an interesting candidate for the HER.
Approaches to control catalytic performance and catalyst stability by intra-molecular non-covalent (Natarajan et al., 2017) (Natarajan et al., 2022) and covalent (Pullen et al., 2019) interactions or the reaction mechanism by introducing additional sites of protonation (Wilson et al., 2006; Liu et al., 2013; Pandey et al., 2021) or modifications of the first coordination shell ligands will allow the design of more efficient catalyst candidates for sustainable hydrogen production.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material. CCDC entry 2189879 (1) contains the crystallographic data for this paper. These data can be obtained free of charge from the Cambridge Crystallographic Data Centre (www.ccdc.cam.ac.uk/data_request/cif).
Author contributions
VK performed the synthesis, crystallization, single-crystal structure determination, and CV experiments. MJ performed the computational work, structure, reaction mechanism, and redox potentials. VK and MJ wrote the first draft of the manuscript. MS and SK-G designed the study, supervised the work, and wrote the final version of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This work was supported by the Max Planck Society for the Advancement of Science (to MS and MJ), the Faculty Research Programme Grant-IoE from the Institution of Eminence (IoE), University of Delhi (Ref. No./IoE/2021/12/FRP to SK-G), University Grants Commission (UGC) and Council of Scientific and Industrial Research (CSIR), New Delhi, India (to VK). Gefördert durch die Deutsche Forschungsgemeinschaft (DFG) - TRR 63 “Integrierte chemische Prozesse in flüssigen Mehrphasensystemen” (Teilprojekt A4) - 56091768. Funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - TRR 63 “Integrated Chemical Processes in Liquid Multiphase Systems” (subproject A4)—56091768.
Acknowledgments
VK is grateful to University Grants Commission (UGC) for a Research fellowship and to Council of Scientific & Industrial Research (CSIR), New Delhi, India for a Senior Research fellowship. MJ and MS thank the Max Planck Society for the Advancement of Science for financial support.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2022.993085/full#supplementary-material
References
Abedin, T. S. M., Moni, M. R., Ghosh, S., Tocher, D. A., Hossain, G. M. G., Mobin, S. M., et al. (2018). Mn2(CO)6(μ-mbi)2 as a precursor for mono- and polynuclear complexes containing the 2-mercaptobenzimidazolate (mbi) ligand. Polyhedron 152, 164–171. doi:10.1016/j.poly.2018.06.039
Appel, A. M., and Helm, M. L. (2014). Determining the overpotential for a molecular electrocatalyst. ACS Catal. 4, 630–633. doi:10.1021/cs401013v
Artero, V., and Fontecave, M. (2005). Some general principles for designing electrocatalysts with hydrogenase activity. Coord. Chem. Rev. 249, 1518–1535. doi:10.1016/j.ccr.2005.01.014
Artero, V., and Saveant, J.-M. (2014). Toward the rational benchmarking of homogeneous H2-evolving catalysts. Energy Environ. Sci. 7, 3808–3814. doi:10.1039/c4ee01709a
Becke, A. D. (1988). Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A . Coll. Park. 38, 3098–3100. doi:10.1103/physreva.38.3098
Becke, A. D. (1993). Density‐functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 98, 5648–5652. doi:10.1063/1.464913
Bigi, J. P., Hanna, T. E., Harman, W. H., Chang, A., and Chang, C. J. (2010). Electrocatalytic reduction of protons to hydrogen by a water-compatible cobalt polypyridyl platform. Chem. Commun. 46, 958–960. doi:10.1039/b915846d
Braga, D., Maini, L., Grepioni, F., Elschenbroich, C., Paganelli, F., and Schiemann, O. (2001). Novel organometallic building blocks for crystal engineering. Synthesis and structural characterization of the dicarboxylic acid [Cr0(η6-C6H5COOH)2], of two polymorphs of its oxidation derivative [CrI(η6-C6H5COOH)2]+[PF6]-and of the zwitterionic form [CrI(η6-C6H5COOH)(η6-C6H5COO)]. Organometallics 20, 1875–1881. doi:10.1021/om0010722
Brazzolotto, D., Gennari, M., Queyriaux, N., Simmons, T. R., Pecaut, J., Demeshko, S., et al. (2016). Nickel-centred proton reduction catalysis in a model of [NiFe] hydrogenase. Nat. Chem. 8, 1054–1060. doi:10.1038/nchem.2575
Bruschi, M., Zampella, G., Fantucci, P., and De Gioia, L. (2005). DFT investigations of models related to the active site of NiFe and Fe hydrogenases. Coord. Chem. Rev. 249, 1620–1640. doi:10.1016/j.ccr.2004.12.007
R. Camack, M. Frey, and R. Robson (Editors) (2001). Hydrogen as a fuel (London: Taylor & Francis). doi:10.1201/9780203471043
Cracknell, J. A., Vincent, K. A., and Armstrong, F. A. (2008). Enzymes as working or inspirational electrocatalysts for fuel cells and electrolysis. Chem. Rev. 108, 2439–2461. doi:10.1021/cr0680639
Dempsey, J. L., Brunschwig, B. S., Winkler, J. R., and Gray, H. B. (2009). Hydrogen evolution catalyzed by cobaloximes. Acc. Chem. Res. 42, 1995–2004. doi:10.1021/ar900253e
Dubois, D. L. (2014). Development of molecular electrocatalysts for energy storage. Inorg. Chem. 53, 3935–3960. doi:10.1021/ic4026969
Duclos, L., Lupsea, M., Mandil, G., Svecova, L., Thivel, P.-X., and Laforest, V. (2017). Environmental assessment of proton exchange membrane fuel cell platinum catalyst recycling. J. Clean. Prod. 142, 2618–2628. doi:10.1016/j.jclepro.2016.10.197
Evans, D. J., and Pickett, C. J. (2003). Chemistry and the hydrogenases. Chem. Soc. Rev. 32, 268–275. doi:10.1039/b201317g
Fang, M., Wiedner, E. S., Dougherty, W. G., Kassel, W. S., Liu, T., Dubois, D. L., et al. (2014). Cobalt complexes containing pendant amines in the second coordination sphere as electrocatalysts for H2 production. Organometallics 33, 5820–5833. doi:10.1021/om5004607
Felton, G. a. N., Glass, R. S., Lichtenberger, D. L., and Evans, D. H. (2006). Iron-only hydrogenase mimics. Thermodynamic aspects of the use of electrochemistry to evaluate catalytic efficiency for hydrogen generation. Inorg. Chem. 45, 9181–9184. doi:10.1021/ic060984e
Fernandez, S., Bandomo, G. C. B., and Lloret-Fillol, J. (2022). “Manganese complexes for electro- and photocatalytic transformations,” in Manganese catalysis in organic synthesis. Editor Sortais, J.-B. (Weinheim: Wiley-VCH), 137–181. doi:10.1002/9783527826131.ch5
Fu, L.-Z., Zhou, L.-L., Tang, L.-Z., Zhang, Y.-X., and Zhan, S.-Z. (2015b). A molecular iron(III) electrocatalyst supported by amine-bis(phenolate) ligand for water reduction. Int. J. Hydrogen Energy 40, 8688–8694. doi:10.1016/j.ijhydene.2015.05.026
Fu, L.-Z., Zhou, L.-L., Tang, L.-Z., Zhang, Y.-X., and Zhan, S.-Z. (2015a). Electrochemical and photochemical-driven hydrogen evolution catalyzed by a dinuclear CoII–CoII complex. J. Power Sources 280, 453–458. doi:10.1016/j.jpowsour.2015.01.121
Grimme, S., Antony, J., Ehrlich, S., and Krieg, H. (2010). A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 132, 154104. doi:10.1063/1.3382344
Helm, M. L., Stewart, M. P., Bullock, R. M., Dubois, M. R., and Dubois, D. L. (2011). A synthetic nickel electrocatalyst with a turnover frequency above 100,000 s−1 for H2 production. Science 333, 863–866. doi:10.1126/science.120586
Hirahara, M., Shoji, A., and Yagi, M. (2014). Artificial manganese center models for photosynthetic oxygen evolution in photosystem II. Eur. J. Inorg. Chem. 2014, 595–606. doi:10.1002/ejic.201300683
Jones, A. K., Sillery, E., Albracht, S. P. J., and Armstrong, F. A. (2002). Direct comparison of the electrocatalytic oxidation of hydrogen by an enzyme and a platinum catalyst. Chem. Commun. 2002, 866–867. doi:10.1039/b201337a
Kaim, V., and Kaur-Ghumaan, S. (2019). Manganese complexes: Hydrogen generation and oxidation. Eur. J. Inorg. Chem. 2019, 5041–5051. doi:10.1002/ejic.201900988
Kaim, V., and Kaur-Ghumaan, S. (2021). Mononuclear Mn complexes featuring N, S-/N, N-donor and 1, 3, 5-triaza-7-phosphaadamantane ligands: Synthesis and electrocatalytic properties. New J. Chem. 45, 20272–20279. doi:10.1039/d1nj02104d
Kaim, V., Natarajan, M., and Kaur-Ghumaan, S. (2019). Dinuclear manganese carbonyl complexes: Electrocatalytic reduction of protons to dihydrogen. ChemistrySelect 4, 1789–1794. doi:10.1002/slct.201803754
Kaur-Ghumaan, S., Schwartz, L., Lomoth, R., Stein, M., and Ott, S. (2010). Catalytic hydrogen evolution from mononuclear iron(II) carbonyl complexes as minimal functional models of the [FeFe] hydrogenase active site. Angew. Chem. Int. Ed. 49, 8033–8036. doi:10.1002/anie.201002719
Kaur-Ghumaan, S., and Stein, M. (2014). [NiFe] hydrogenases: How close do structural and functional mimics approach the active site? Dalton Trans. 43, 9392–9405. doi:10.1039/c4dt00539b
Klamt, A., and Schüürmann, G. (1993). Cosmo: A new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient. J. Chem. Soc. Perkin Trans. 2, 799–805. doi:10.1039/p29930000799
Lee, C., Yang, W., and Parr, R. G. (1988). Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 37, 785–789. doi:10.1103/physrevb.37.785
Liu, T., Dubois, D. L., and Bullock, R. M. (2013). An iron complex with pendent amines as a molecular electrocatalyst for oxidation of hydrogen. Nat. Chem. 5, 228–233. doi:10.1038/nchem.1571
Lubitz, W., and Tumas, W. (2007). Hydrogen: An overview. Chem. Rev. 107, 3900–3903. doi:10.1021/cr050200z
Natarajan, M., Faujdar, H., Mobin, S. M., Stein, M., and Kaur-Ghumaan, S. (2017). A mononuclear iron carbonyl complex [Fe(μ-bdt)(CO)2(PTA)2] with bulky phosphine ligands: A model for the [FeFe] hydrogenase enzyme active site with an inverted redox potential. Dalton Trans. 46, 10050–10056. doi:10.1039/c7dt01994g
Natarajan, M., Kumar, N., Joshi, M., Stein, M., and Kaur-Ghumaan, S. (2022). Mechanism of diiron hydrogenase complexes controlled by nature of bridging dithiolate ligand. ChemistryOpen 11, e202100238. doi:10.1002/open.202100238
Onoda, A., and Hayashi, T. (2015). Artificial hydrogenase: Biomimetic approaches controlling active molecular catalysts. Curr. Opin. Chem. Biol. 25, 133–140. doi:10.1016/j.cbpa.2014.12.041
Pan, H.-J., and Hu, X. (2020). Biomimetic hydrogenation catalyzed by a manganese model of [Fe]-hydrogenase. Angew. Chem. Int. Ed. Engl. 59, 4972–4976. doi:10.1002/ange.201914377
Pan, H.-J., Huang, G., Wodrich, M. D., Tirani, F. F., Ataka, K., Shima, S., et al. (2019). A catalytically active [Mn]-hydrogenase incorporating a non-native metal cofactor. Nat. Chem. 11, 669–675. doi:10.1038/s41557-019-0266-1
Pandey, I. K., Agarwal, T., Mobin, S. M., Stein, M., and Kaur-Ghumaan, S. (2021). Switching site reactivity in hydrogenase model systems by introducing a pendant amine ligand. ACS Omega 6, 4192–4203. doi:10.1021/acsomega.0c04901
Penner, S. S. (2005). Steps toward the hydrogen economy. Energy (Amsterdam, Neth.) 31, 33–43. doi:10.1016/j.energy.2004.04.060
Perdew, J. P. (1986). Density-functional approximation for the correlation energy of the inhomogeneous electron gas. Phys. Rev. B 33, 8822–8824. doi:10.1103/physrevb.33.8822
Pullen, S., Maji, S., Stein, M., and Ott, S. (2019). Restricted rotation of an Fe(CO)2(PL3)-subunit in [FeFe]-hydrogenase active site mimics by intramolecular ligation. Dalton Trans. 48, 5933–5939. doi:10.1039/c8dt05148h
Rodriguez-Maciá, P., Dutta, A., Lubitz, W., Shaw, W. J., and Rüdiger, O. (2015). Direct comparison of the performance of a bio-inspired synthetic nickel catalyst and a [NiFe]-Hydrogenase, both covalently attached to electrodes. Angew. Chem. Int. Ed. 54, 12303–12307. doi:10.1002/anie.201502364
Schäfer, A., Klamt, A., Sattel, D., Lohrenz, J. C. W., and Eckert, F. (2000). COSMO Implementation in TURBOMOLE: Extension of an efficient quantum chemical code towards liquid systems. Phys. Chem. Chem. Phys. 2, 2187–2193. doi:10.1039/b000184h
Schilter, D., Camara, J. M., Huynh, M. T., Hammes-Schiffer, S., and Rauchfuss, T. B. (2016). Hydrogenase enzymes and their synthetic models: The role of metal hydrides. Chem. Rev. 116, 8693–8749. doi:10.1021/acs.chemrev.6b00180
Simmons, T. R., Berggren, G., Bacchi, M., Fontecave, M., and Artero, V. (2014). Mimicking hydrogenases: From biomimetics to artificial enzymes. Coord. Chem. Rev. 270, 127–150. doi:10.1016/j.ccr.2013.12.018
Stephens, P. J., Devlin, F. J., Chabalowski, C. F., and Frisch, M. J. (1994). Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields. J. Phys. Chem. 98, 11623–11627. doi:10.1021/j100096a001
Straistari, T., Fize, J., Shova, S., Réglier, M., Artero, V., and Orio, M. (2017). A thiosemicarbazone–nickel(II) complex as efficient electrocatalyst for hydrogen evolution. ChemCatChem 9, 2262–2268. doi:10.1002/cctc.201600967
Sun, L., Duboc, C., and Shen, K. (2022). Bioinspired molecular electrocatalysts for H2 production: Chemical strategies. ACS Catal. 12, 9159–9170. doi:10.1021/acscatal.2c02171
Vosko, S. H., Wilk, L., and Nusair, M. (1980). Accurate spin-dependent electron liquid correlation energies for local spin density calculations: A critical analysis. Can. J. Phys. 58, 1200–1211. doi:10.1139/p80-159
Weigend, F., and Ahlrichs, R. (2005). Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 7, 3297–3305. doi:10.1039/b508541a
Wiechen, M., Najafpour, M. M., Allakhverdiev, S. I., and Spiccia, L. (2014). Water oxidation catalysis by manganese oxides: Learning from evolution. Energy Environ. Sci. 7, 2203–2212. doi:10.1039/c4ee00681j
Wilson, A. D., Newell, R. H., McNevin, M. J., Muckerman, J. T., Rakowski Dubois, M., and Dubois, D. L. (2006). Hydrogen oxidation and production using nickel-based molecular catalysts with positioned proton relays. J. Am. Chem. Soc. 128, 358–366. doi:10.1021/ja056442y
Keywords: bioinorganic chemistry, manganese catalyst, hydrogen evolution, hydrogenase, redox activity, reaction mechanism, computational catalysis
Citation: Kaim V, Joshi M, Stein M and Kaur-Ghumaan S (2022) Mononuclear manganese complexes as hydrogen evolving catalysts. Front. Chem. 10:993085. doi: 10.3389/fchem.2022.993085
Received: 13 July 2022; Accepted: 05 August 2022;
Published: 07 October 2022.
Edited by:
José Cleiton Sousa dos Santos, University of International Integration of Afro-Brazilian Lusophony, BrazilReviewed by:
Nathan John DeYonker, University of Memphis, United StatesJuliana Serpa, University of International Integration of Afro-Brazilian Lusophony, Brazil
Copyright © 2022 Kaim, Joshi, Stein and Kaur-Ghumaan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Matthias Stein, matthias.stein@mpi-magdeburg.mpg.de; Sandeep Kaur-Ghumaan, skaur@chemistry.du.ac.in