A Molecular MoS2 Edge Site Mimic for Catalytic Hydrogen Generation

See allHide authors and affiliations

Science  10 Feb 2012:
Vol. 335, Issue 6069, pp. 698-702
DOI: 10.1126/science.1215868


Inorganic solids are an important class of catalysts that often derive their activity from sparse active sites that are structurally distinct from the inactive bulk. Rationally optimizing activity is therefore beholden to the challenges in studying these active sites in molecular detail. Here, we report a molecule that mimics the structure of the proposed triangular active edge site fragments of molybdenum disulfide (MoS2), a widely used industrial catalyst that has shown promise as a low-cost alternative to platinum for electrocatalytic hydrogen production. By leveraging the robust coordination environment of a pentapyridyl ligand, we synthesized and structurally characterized a well-defined MoIV-disulfide complex that, upon electrochemical reduction, can catalytically generate hydrogen from acidic organic media as well as from acidic water.

Molybdenite, or molybdenum disulfide (MoS2), the earliest form of molybdenum to be identified from ores, is one of the most widely used catalysts in industry today as the standard for hydrodesulfurization (HDS) of petroleum (1). In its nanoparticulate form, MoS2 has further demonstrated promise as an inexpensive alternative to platinum for the electrochemical and photochemical generation of hydrogen from water (26). As is the case with many inorganic solids, the catalytic activity of MoS2 is localized to rare surface sites, whereas the bulk material is relatively inert (79). High-resolution scanning tunneling microscopy studies and theoretical calculations performed on nanoparticulate MoS2 structures that form under sulfiding conditions implicate the formation of disulfide linkages or triangular MoS2 units along the fully sulfided catalytically active edges of the layered structure (1014). However, the precise molecular structures and modes of action of these sites remain elusive. Because of the bulk material’s layered structure, which favors the growth of plate-like crystals, a single crystal with a large edge dimension is extremely challenging to prepare (1). Here, we report the synthesis of a well-defined molecular analog of the proposed MoS2 edge structure, a side-on bound MoIV-disulfide complex. Electrochemical reduction of this molecule leads to the catalytic generation of hydrogen from acidic organic media as well as from acidic water, lending support to the proposed active site morphology in the more active heterogeneous catalyst.

Recent efforts from our laboratories have explored the chemistry of the PY5Me2 ligand [PY5Me2 = 2,6-bis(1,1-bis(2-pyridyl)ethyl)pyridine] and the ability of its MoIV-oxo complex [(PY5Me2)MoO]2+ to catalyze the production of hydrogen from neutral water (15). Along these lines, treatment of the MoII precursor [(PY5Me2)Mo(CF3SO3)]1+ with S8 at room temperature affords the MoIV-disulfide complex [(PY5Me2)MoS2]2+ (1) in 74% yield (Fig. 1) (16). This molecule represents a rare species in molybdenum coordination chemistry (1720) and a discrete molecular analog of the proposed triangular MoS2 edge sites of molybdenite. A single-crystal x-ray analysis of 1 revealed a distorted pseudo-octahedral geometry around the metal center, with the MoS2 plane twisting away from the plane of the axial pyridyl ligand by approximately 46°. The PY5Me2 ligand is slightly distorted around the central metal with a mean Mo-N distance of 2.19(2) Å and a mean N-Mo-N angle of 82(6)°. The mean Mo-S distance of 2.402(2) Å, S-S distance of 2.019(1) Å and S-Mo-S angle of 49.68(3)° are similar to those reported for Cp*2MoS2 (21) and CptBu2MoS2 (22), and together frame a molecular MoS2 isosceles triangle supported by the PY5Me2 ligand. The longer Mo-S distance compared with the 2.254(2) Å bond length for a terminal sulfido ligand bonded to MoIV (18) indicates single Mo-S bonds. Similarly, the S-S bond length is close to the distance of 2.060 Å observed between S atoms in S8 (23) as well as the distance between S atoms in a dinuclear MoV complex bearing a μ-S22– bridging ligand (24), consistent with a single bond between the two S atoms.

Fig. 1

Reaction of [(PY5Me2)Mo(CF3SO3)]1+ with sulfur to form [(PY5Me2)MoS2]2+ (1, top), and a model of the layered structure of MoS2 (bottom), highlighting the proposed surface reconstruction to give disulfide-terminated edge sites. The molecular structures depicted at the top are the results of single-crystal x-ray analyses. Green, yellow, light blue, red, blue, and gray spheres represent Mo, S, F, O, N, and C atoms, respectively; H atoms are omitted for clarity.

The cyclic voltammogram on a glassy carbon disc electrode (GCDE) of a 1-mM solution of 1 in acetonitrile with 0.1 M of (Bu4N)PF6 electrolyte shows a reversible set of redox peaks at E1/2 = 1.27 V versus the standard hydrogen electrode (SHE) (Fig. 2A), which we assigned to the [(PY5Me2)MoS2]2+/3+ couple. Scanning cathodically reveals two more reversible redox couples at E1/2 = –0.20 and –0.82 V versus SHE assigned to the [(PY5Me2)MoS2]2+/1+ and [(PY5Me2)MoS2]1+/0 couples. At still more negative potentials, a larger irreversible wave rising to a maximum at Ep = –1.49 V versus SHE is presumably due to an electrochemical-chemical (EC) process, perhaps involving the S22– ligand (fig. S2).

Fig. 2

Cyclic voltammograms on a glassy carbon disc electrode (A = 7.07 × 10−2 cm2) of (A) an acetonitrile solution containing 1 mM of [(PY5Me2)MoS2](CF3SO3)2 (1) and 0.1 M of (Bu4N)PF6 at a scan rate of 0.1 V/s; (B) acetonitrile solutions containing 0.5 mM of 1 and 0.1 M of (Bu4N)PF6 in the presence (blue line) and absence (red line) of 160 equiv. of acetic acid (79 mM), or 79 mM of the acid alone (gray line), at a scan rate of 0.1 V/s; (C) 1 M aqueous acetate buffer at pH 3 in the presence (red and blue lines) and absence (gray line) of 1.1 mM of 1 at a scan rate of 0.1 V/s. (D) Cyclic voltammograms on a mercury drop electrode (A ~ 1.16 × 10−2 cm2) of 0.05 M aqueous KHP buffer at pH 3 in the presence (red and blue lines) and absence (gray line) of 160 μM of 1 at a scan rate of 0.5 V/s. In (C) and (D), the red line indicates the initial scan and the blue line indicates the subsequent scan.

The accessibility of multiply reduced states for the MoIVS2 complex 1 prompted us to investigate its ability to reduce protons in organic solvents. The cyclic voltammogram (on a GCDE) of a 0.5-mM solution of 1 in acetonitrile containing 79 mM (160 equivalents) of acetic acid shows a sharp catalytic current after the MoIII/II couple with an onset potential of about –1.00 V versus SHE. The current reaches half its maximum value at –1.27 V and peaks at –1.50 V versus SHE. (Fig. 2B). The peak currents of the precatalytic waves as well as the catalytic wave vary linearly with the square root of the scan rate (fig. S4), indicating diffusion-limited redox processes, with the electrochemically active species freely diffusing in solution. The catalytic nature of the reductive wave at about –1.27 V versus SHE was verified by passing 3.00 C of charge during a bulk electrolysis. The experiment was conducted on a 55-μM solution of 1 in acetonitrile with 9 mM of acetic acid, using a glassy carbon plate electrode in a double compartment cell with the potential held at –1.09 V versus SHE for 2.75 hours. A control experiment performed under identical conditions, but without 1, showed the accumulation of just 0.679 C of charge (fig. S5). Subtracting the background contribution of direct acid reduction leads to 10.9 equiv. of charge generated per molybdenum center, which precludes a stoichiometric reaction with respect to molybdenum. Gas chromatography of the headspace gas during the bulk electrolysis confirmed the generation of H2.

Encouraged by this clear indication that 1 can operate as a molecular electrocatalyst in acetonitrile, we chose to explore its potential utility for the electrocatalytic reduction of protons in water, a much more attractive medium for the sustainable generation of hydrogen. Indeed, a sharp increase in current, indicative of catalytic water reduction, is evident in cyclic voltammograms (on a GCDE) of aqueous solutions of 1 at pH values of 3, 4, 5, and 6, with optimal activity observed at pH 3. We therefore conducted subsequent studies in aqueous media at pH 3 using potassium hydrogen phthalate (KHP) or sodium acetate aqueous buffers. The voltammogram (on a GCDE) of a 1.2-mM solution of 1 in aqueous acetate buffer at pH 3, acquired at a scan rate of 0.1 V/s, shows two irreversible reductions at peak potentials of −0.34 and −0.53 V versus SHE, respectively, before the onset of the catalytic current at a potential of about −0.65 V versus SHE (Fig. 2C). The voltammogram of the molybdenum-oxo compound [(PY5Me2)MoO]2+, taken under identical conditions, shows no catalytic activity within the solvent electrochemical window, confirming that the sulfide ligands are not hydrolyzed under acidic aqueous conditions. As observed in acetonitrile solution, the current response of the redox events of 1 on the glassy carbon electrode in acidic water shows a linear dependence on the square root of the scan rate, which is indicative of a diffusion-controlled process (figs. S6 and S7). A controlled potential electrolysis of a 92-μM solution of 1 in 3 M aqueous acetate buffer at pH 3 conducted using a glassy carbon rod electrode in a double-compartment cell at −0.75 V versus SHE for 1 hour led to the accumulation of 2.42 C of charge. In the absence of 1, an identical experiment led to the passage of only 0.015 C of charge (fig. S8). Subtraction of this background activity affords 14.3 equiv. of charge per molybdenum center, establishing the catalytic nature of the reduction. To determine whether deposition on the electrode occurs during an extended electrolysis, the electrode was carefully rinsed with distilled water and examined using molybdenum K-edge x-ray absorption spectroscopy, which indicated no molybdenum species present to the limit of x-ray fluorescence measurement (monolayer). Furthermore, the electronic adsorption spectra of the solution before and after the 1-hour electrolysis show that the (PY5Me2)Mo unit remains intact during that period (fig. S9).

Because of the background activity of direct water reduction at the glassy carbon electrode at high overpotentials, subsequent experiments were conducted on a mercury electrode, where no direct water reduction occurs below overpotentials of about −1 V (overpotential being the difference between the applied potential and the thermodynamic potential for water reduction at the same pH) (25). Similar to the data obtained on a glassy carbon electrode, a cyclic voltammogram on a mercury drop electrode recorded at a scan rate of 0.5 V/s shows two reductions, with a sharp wave corresponding to the first reduction at –0.16 V versus SHE, followed by a quasireversible couple at E1/2 = –0.38 V preceding the catalytic current (Fig. 2D). To probe the nature of the first reduction process, we investigated the pH dependence of the cathodic peak potential (Vc), summarized in the Pourbaix diagram depicted in Fig. 4 (assuming electrochemical reversibility for this redox couple). Here, Vc declines linearly with rising solution pH at a slope of –59.8(±0.8), which is within experimental error of the ideal value of –59 mV, corresponding to the association of one proton with the electron transfer (25). The observed pH dependence, taken together with the irreversibility of this peak, suggests a geometric distortion about the metal center consistent with the formation of a protonated reduced Mo-disulfide unit or cleavage of a Mo-S bond. Formation of a S-H bond is the first step in the proposed pathway for electrocatalytic proton reduction at MoS2 nanoparticle edge sites (10), as well as in dinuclear (CpMoμ-S)2S2CH2 complexes (26). Protonation of the bridging sulfido ligand is also proposed for H2 evolution by the nitrogenase enzyme (27). Furthermore, terminal S-H bonds at the edges of MoS2 HDS catalysts have been assigned, which has led to speculation that such species are the source of H2 for the hydrogenation of the desulfurized molecules (28). In contrast, the second reductive wave does not show a similar pH dependence, leading us to hypothesize that this process involves one electron reduction of the complex, which upon subsequent reduction and protonation forms the catalytically active species. Close contact between metal centers adsorbed on the electrode may also facilitate a multinuclear pathway, perhaps through the formation of S-H bonds that could homolytically or heterolytically cleave to yield H2.

Fig. 4

(A) Cyclic voltammograms of 0.3-mM solutions of [(PY5Me2)MoS2](CF3SO3)2 (1) in 0.05 M aqueous KHP buffers at pH values ranging from 3 to 7 on a Hg drop electrode at a scan rate of 1 V/s. (B) Pourbaix diagram showing the dependence of the peak potential of the first reductive wave on the solution pH: A linear fit gives a slope of −59.8 (± 0.8).

Analysis of the peak current of the cathodic and anodic scans of the second quasireversible redox couple as a function of scan rate (fig. S10) shows a linear relation, consistent with adsorption of the molecules on the electrode surface. Stabilization of the reduced molybdenum complexes through adsorption on the Hg electrode would be consistent with the positive shift of the redox events compared with the diffusion-limited redox processes observed with glassy carbon electrodes (25). For an idealized adsorption, which assumes reversible electron transfer with no interaction between adsorbed species, the number of electrons that give rise to the wave corresponds to 90.6/∆Ep,1/2, where ∆Ep,1/2 is the width of the wave at half the maximum current (25). At scan rates of 200 to 1000 mV/s, a mean width of 74(2) mV is observed for the anodic scan, which indicates a one-electron process. By observing the current response to variations in scan rate, we can quantify the coverage of metal sites on the electrode surface to be 3.4 to 3.6 × 10−10 mols/cm2. A rough calculation of the footprint of 1, obtained using the dimensions from the crystal structure, suggests monolayer coverage on the electrode [see supporting online material (SOM)]. Due to possible increases in surface area of the mercury when the pool is stirred during electrolysis, we use an estimated surface coverage of 10−10 mols/cm2 for subsequent hydrogen evolution rate calculations.

The cyclic voltammogram of a 130-μM solution of 1 at pH 3 shows a catalytic current with an overpotential of about −400 mV (−0.58 V versus SHE) (Fig. 3B). The catalytic performance of [(PY5Me2)MoO]2+ was once again evaluated at the same pH to establish that 1 did not form the molybdenum-oxo complex under aqueous, reducing conditions. As depicted in Fig. 3B, [(PY5Me2)MoO]2+ shows the onset of a catalytic current at an overpotential of −700 mV (−0.88 V versus SHE) under the same conditions. Thus, using the disulfide complex 1 gives an improvement of 300 mV in overpotential with respect to the analogous oxo complex for water reduction at pH 3. Controlled potential electrolysis experiments were conducted to assess the rate of hydrogen production at various overpotentials. Identical measurements were performed with and without the catalyst to subtract the background activity at each applied potential. Figure S11 shows the charge accumulated over 1-min intervals, with applied overpotential increasing from −428 mV until a saturation value of 6.36 C is reached at −828 mV (due to the voltage between the working and auxiliary electrodes exceeding the maximum voltage obtainable by the potentiostat at high current densities). Under these conditions, the pH change of the solution, as well as quantitative gas chromatographic analysis of the electrolysis-cell headspace, indicated that the catalyst performs at close to 100% Faradaic efficiency, where every electron is used for the generation of hydrogen (see figs. S12 and S13 and SOM). At an overpotential of −828 mV, the observed turnover frequency (TOF) reached a maximum of 280 moles of H2 per mole of catalyst per second.

Fig. 3

(A) Scan rate dependence of precatalytic waves for a 160-μM solution of [(PY5Me2)MoS2](CF3SO3)2 (1) in 0.05 M aqueous KHP buffer at pH 3, at scan rates ranging from 100 to 1000 mV/s in 100 mV/s increments on a mercury drop electrode (A ~ 11.6 × 10−3 cm2). (B) Cathodic scans of 130-μM solutions of 1 (orange line) or [(PY5Me2)MoO](CF3SO3)2 (green line) in 1 M aqueous acetate buffer at pH 3 and of the buffer solution alone (blue line), on a mercury pool electrode (A = 19.6 cm2) at a scan rate of 100 mV/s. (C) TOF versus overpotential for a 170-μM solution of 1 in seawater maintained at pH 3 with 1 M acetate buffer. The background solvent activity has been subtracted from the data. (D) Charge build up over time from electrolysis of a 66-μM solution of 1 (orange circles) in 3 M aqueous acetate buffer at pH 3, compared with the electrolyte solution alone (blue circles), with the cell operating at an overpotential of −780 mV. Overpotential = applied potential – E(pH 3).

We also evaluated the catalytic properties of 1 using a sample of California seawater, which was buffered at pH 3 by the addition of acetic acid and sodium acetate. As shown in Fig. 3C, the rates of hydrogen evolution for a given overpotential were similar to those observed in the studies conducted in distilled water, but with an even higher optimal TOF of 480 moles of H2 per mole of catalyst per second at an overpotential of −780 mV. This result highlights the robustness of the catalyst to the impurities found in seawater.

To assess the long-term stability of the catalyst, an extended electrolysis was conducted in a 3 M pH 3 acetate buffer at an overpotential of −780 mV (Fig. 3D). The total turnover number (TON) reached ~3.5 × 103 moles of H2 per mole of catalyst, showing that the catalyst is stable for long durations in aqueous media. This value is the lower bound for the catalyst lifetime using the total number of molecules in solution. Calculating the TON using the surface coverage of catalyst molecules on the electrode affords 1.9 × 107 moles of H2 per mole of catalyst. We report both numbers because we do not yet know the rate of exchange of the surface layer with the molecules in solution during the 23-hour period.

Comparing the electronic adsorption spectrum of a solution of the catalyst electrolyzed for 1 hour against a control shows that the (PY5Me2)Mo unit remains intact during this period (fig. S14). In addition, a 1H nuclear magnetic resonance study to quantify the amount of free ligand in solution after a 12-hour electrolysis indicated that less than 3% is lost from the initial catalyst loading, further supporting the stability of this molecular species under reductive, aqueous conditions (see SOM for details). However, adsorption on the mercury surface complicates the characterization of the catalytically active species, which we cannot as yet unambiguously identify. The analogous studies performed using glassy carbon electrodes in organic and aqueous media, where the catalyst is not adsorbed on the electrode surface, support a molecular species as a competent catalyst.

This first-generation molybdenum-disulfide complex exhibits considerably higher stability and comparable TOFs to hydrogenase enzymes (which show rates from 102 to 104 s−1) (29, 30), albeit with greater overpotentials. For comparison, water-tolerant cobalt-based molecular catalysts attached to glassy carbon electrodes can achieve TOFs of ~90 s–1 at an overpotential of −400 mV (31), and experimentally determined TOFs ranging from 2.3 × 10−4 to 2.2 × 10−3 s−1 at overpotentials ranging from −690 to −1040 mV have been reported for a cobalt complex on a mercury electrode (32); an analogous Co complex bearing the PY5Me2 ligand shows a TOF of ~0.3 s–1 at an overpotential of −900 mV on a mercury electrode (33). Furthermore, TOFs of ~0.07 s−1 at the equilibrium potential have been reported for carbon-supported [Mo3S4]4+ clusters (34), which are susceptible to desorption from the electrode. Other examples of molecular catalysts attached to electrodes display lower activity and lifetimes in aqueous solution (3537). A key advantage of 1 is the ability to form a layer of MoS2 units, analogous to constructing a sulfided edge of MoS2, where the dimensions of the layer can be defined by the size and shape of the electrode. Here, the electronic structure of the MoS2 units, and thereby perhaps the activity, stability, and required overpotential for proton reduction, may be adjusted through ligand modifications accessible through synthetic chemistry.

The ability to prepare, characterize, and evaluate molecular analogs of the active components of inorganic solids has broad implications for the design and optimization of functional metal sites, not the least of which is control over the density of these units. For example, recent electronic structure calculations conducted on nanoparticulate MoS2 indicate that only a quarter of the edge sites are used for hydrogen production (10). Increasing the number of active edge sites per unit volume by tailoring progressively smaller nanostructures or changing the electronics of the system to increase the enthalpy of hydrogen adsorption is a major challenge in inorganic materials and nanoscience. We present an alternative strategy using discrete molecular units, which in principle can be tailored to give a high density of catalytically active metal sites without the rest of the inactive bulk material.

Supporting Online Material

Materials and Methods

Figs. S1 to S14

Tables S1 and S2


References and Notes

  1. Acknowledgments: Initial synthetic and electrochemical studies were supported by the Department of Energy (DOE) Lawrence Berkeley National Laboratory Helios Solar Energy Research Center (51HE112B to C.J.C.). Further electrochemical studies are based on work performed by the Joint Center for Artificial Photosynthesis, a DOE Energy Innovation Hub, supported through the Office of Science of the U.S. Department of Energy under award DE-SC0004993 (J.R.L.). C.J.C. is an investigator with the Howard Hughes Medical Institute. We thank B. Lassalle, J. Yano, and V. Yachandra for obtaining the x-ray absorption spectroscopy spectra; J. Breen for fabrication of electrochemical cells; Y. Hwang for assistance with gas chromatography; A. T. Iavarone for obtaining the mass spectra; and the University of California President’s Postdoctoral Fellowship program for support of E.M. Structural parameters for 1 are available free of charge from the Cambridge Crystallographic Data Centre under identification number CCDC-798834. A patent has been filed by Lawrence Berkeley National Laboratory based on the results presented herein.
View Abstract

Stay Connected to Science

Navigate This Article