Skip to main content

ORIGINAL RESEARCH article

Front. Chem., 07 October 2022
Sec. Catalytic Reactions and Chemistry
This article is part of the Research Topic Hot Topic: Chemical Reactions & Catalysis for a Sustainable Future View all 6 articles

Mononuclear manganese complexes as hydrogen evolving catalysts

  • 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
www.frontiersin.org

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 14 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)32-S2NC7H4)(PPh2(PhCOOH-p))] 1, fac-[(Mn(CO)32-S2NC7H4)(PPh3)] 2 and fac-[Mn(CO)32-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 14 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
www.frontiersin.org

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 24 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
www.frontiersin.org

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 14 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
www.frontiersin.org

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 14) or Mn‒S bond distances (2.49, 2.49, 2.52, and 2.52 Å, in 14), nor axial Mn‒P bond distances (2.37, 2.36, 2.36, and 2.36 Å, in 14) 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 14 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 14 (Supplementary Tables S3–S6).

Characterization of redox activity

Mn-complex reduction in the absence of acid

The intrinsic redox behavior of complexes 14 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 Epc1 = ‒1.83 (1), ‒1.99 (3) V and E2pc = ‒2.19 (1), ‒2.30 (3) V (Figure 3).

FIGURE 3
www.frontiersin.org

FIGURE 3. Cyclic voltammograms for reduction of complexes 14 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 Epc1 = ‒1.93, ‒2.06 and E2pc = ‒2.07, ‒2.36 (Figure 3) vs. Fc/Fc+, respectively. This clearly indicates only minor control of redox transitions by equatorial and axial coordination. The irreversibility of the redox transitions clearly indicates large and irreversible structural re-arrangements.

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).

1+e 1(1)
1+e 12(2)
1+2e 12(3)

The reduction potentials of complexes 14 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 14, 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 14 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 14 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
www.frontiersin.org

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 14 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
www.frontiersin.org

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
www.frontiersin.org

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 acidbase 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).

EAH/H2, A0=EH +/H202.3RT/F×pKa(1)

Figure 6 shows the increase of catalytic current and TOF with acid concentration. Higher icatip and TOF values for complex 4 are slightly larger than the values for complexes 13, indicating faster turnover of complex 4.

FIGURE 6
www.frontiersin.org

FIGURE 6. (A) icatip (first catalytic peak) vs [acid]; (B) dependence of TOF (s−1) on [acid] for complexes 1 (■), 2 (●), 3 (▼), and 4 (▲) (1.5 mM) in 0.1 M [N(n−Bu4)][PF6]/CH3CN.

The Mn complexes 14 showed activity with lower TOF/TON values than previously reported Mn complexes but improved icatip values at comparable overpotentials (Kaim et al., 2019; Kaim and Kaur-Ghumaan, 2019; Pan et al., 2019; Pan and Hu, 2020; Yang et al., 2021). The mononuclear Mn complex with heterocyclic (N,S; S,S) and phosphaadamantane ligands (Kaim and Kaur-Ghumaan, 2021) showed activity in both organic and partially aqueous media. The Mn complexes reported in this article open the possibility of further exploring electrocatalysis in aqueous media and also the path of immobilizing the catalyst onto electrode surfaces via the -COOH substituted phosphine ligand.

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
www.frontiersin.org

SCHEME 3. Proton reduction and hydrogen evolution mechanism of complexes 1 and 3.

SCHEME 4
www.frontiersin.org

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 14 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 14 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 (14)2 at −1.9 V and the second one-electron transition (14)1/(14)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

CrossRef Full Text | Google Scholar

Appel, A. M., and Helm, M. L. (2014). Determining the overpotential for a molecular electrocatalyst. ACS Catal. 4, 630–633. doi:10.1021/cs401013v

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Becke, A. D. (1993). Density‐functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 98, 5648–5652. doi:10.1063/1.464913

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

R. Camack, M. Frey, and R. Robson (Editors) (2001). Hydrogen as a fuel (London: Taylor & Francis). doi:10.1201/9780203471043

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Dubois, D. L. (2014). Development of molecular electrocatalysts for energy storage. Inorg. Chem. 53, 3935–3960. doi:10.1021/ic4026969

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Evans, D. J., and Pickett, C. J. (2003). Chemistry and the hydrogenases. Chem. Soc. Rev. 32, 268–275. doi:10.1039/b201317g

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Lubitz, W., and Tumas, W. (2007). Hydrogen: An overview. Chem. Rev. 107, 3900–3903. doi:10.1021/cr050200z

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Penner, S. S. (2005). Steps toward the hydrogen economy. Energy (Amsterdam, Neth.) 31, 33–43. doi:10.1016/j.energy.2004.04.060

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Schlögl, R. (Editor) (2013). Chemical energy storage. Berlin: De Gruyter. doi:10.1515/9783110266320

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Sorais, J.-B. (2022). Manganese catalysis in organic synthesis. Weinheim: Wiley VCH.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., He, S., Wu, Q., Zhang, P., Chen, L., Huang, F., et al. (2021). A bio-inspired mononuclear manganese catalyst for high-rate electrochemical hydrogen production. Dalton Trans. 50, 4783–4788. doi:10.1039/d1dt00672j

PubMed Abstract | CrossRef Full Text | Google Scholar

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, Brazil

Reviewed by:

Nathan John DeYonker, University of Memphis, United States
Juliana 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, bWF0dGhpYXMuc3RlaW5AbXBpLW1hZ2RlYnVyZy5tcGcuZGU=; Sandeep Kaur-Ghumaan, c2thdXJAY2hlbWlzdHJ5LmR1LmFjLmlu

Disclaimer: 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.