Redox Mediator-Assisted Electrocatalytic Hydrogen Evolution from

18 hours ago - The Zr6-based metal–organic framework, NU-1000, was successfully functionalized with candidate catalysts – MoSx units – via SIM ...
0 downloads 0 Views 5MB Size
Subscriber access provided by University of Sunderland

Article

Redox Mediator-Assisted Electrocatalytic Hydrogen Evolution from Water by a Molybdenum Sulfide-Functionalized Metal–Organic Framework Hyunho Noh, Chung-Wei Kung, Ken-ichi Otake, Aaron W. Peters, Zhanyong Li, Yijun Liao, Xinyi Gong, Omar K. Farha, and Joseph T. Hupp ACS Catal., Just Accepted Manuscript • DOI: 10.1021/acscatal.8b02921 • Publication Date (Web): 12 Sep 2018 Downloaded from http://pubs.acs.org on September 12, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

Redox Mediator-Assisted Electrocatalytic Hydrogen Evolution from Water by a Molybdenum Sulfide-Functionalized Metal–Organic Framework Hyunho Noh,† Chung-Wei Kung,† Ken-ichi Otake,† Aaron W. Peters,† Zhanyong Li,† Yijun Liao,† Xinyi Gong,† Omar K. Farha†‡* and Joseph T. Hupp†* †

Department of Chemistry, and ‡Department of Chemical and Biomolecular Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States ‡ Department of Chemistry, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia ABSTRACT: The Zr6-based metal–organic framework, NU-1000, was successfully functionalized with candidate catalysts – MoSx units – via SIM (solvothermal deposition in MOFs) of molybdenum(VI), followed by reaction with H2S gas. The structure of the material, named MoSx-SIM, was characterized spectroscopically and through a single crystal X-ray diffraction measurement. These measurements and others established that the catalyst is monometallic, with mixed oxygen and sulfur coordination, and that it forms from a MOF-node-supported molybdenum-based cluster featuring only oxy ligands. Notably, the formal potential for the MOFgrafted complex, like that for the metal-sulfur active-site of hydrogenase, is nearly coincident with the formal potential for the hydrogen couple. Its effective concentration within the mesoporous MOF is several hundred millimolar, and its porous-frameworkbased immobilization/heterogenization obviates the need for aqueous solubility as a condition for use as a well-defined catalyst. MoSx-SIM was evaluated as an electrocatalyst for evolution of molecular hydrogen from aqueous acid. Although the MoS xfunctionalized framework exhibits catalytic behavior, the highly insulating nature of the support inhibits high electrocatalytic performance. Introduction of an archetypal redox mediator (RM), methyl viologen (MV2+), resulted in more than 20-fold enhancement in its catalytic performance on a turnover frequency (TOF) basis, implying efficient RM-assisted electron transfer (ET) to otherwise electrochemically non-addressable MoSx moieties. Electrochemical kinetic studies with additional viologens as mediators reveal an unexpected square-root dependence of overall reaction rate on mediator concentration, as well as sensitivity to the strength of RMŸ+ as a reductant. These observations, together with observations of potential-dependent H/D isotope effects and potential-dependent pH effects, provide useful insights into the catalysis mechanism and help to explain how the MOF-affixed monometallic catalyst can effectively catalyze a two-electron reduction reaction, i.e. hydrogen evolution from acidified water.

KEYWORDS Hydrogen evolution reaction, electrocatalysis, redox mediator, metal–organic framework, molybdenum sulfide, site-isolated catalyst

INTRODUCTION Molecular hydrogen is used in enormous quantity both for ammonia synthesis1-2 and for syngas-based (i.e., FischerTropsch based) production of hydrocarbons.3-4 It also finds limited, but growing, commercial use in electricity-producing fuels cells, including cells that power vehicles.5-6 H2 is the ultimate clean-fuel as its direct combustion (or its indirect, but much more efficient, electrochemical consumption in fuel cells) produces only water as a chemical product. Most molecular hydrogen is derived from fossil fuels – for example, from methane via steam reforming or from coal via gasification.7 In principle, a much more sustainable, and carbon-neutral, route would be electrochemical or photoelectrochemical water splitting (the reverse of fuel-cell chemistry) with solar photons as the energy input.8-10 Either variant requires electrocatalysts, both for oxygen evolution and for the hydrogen evolution reaction (HER).10-14 While state-of-the-art electrocatalysts for hydrogen evolution generally contain precious and expensive metals, most often platinum,8,15-17 earth abundant and inexpensive molybdenum sulfide (MoSx)-based catalysts have, at least at the laboratory scale, exhibited prom-

isingly high electrocatalytic activity.18-25 Experimental21-22,26-28 and computational29-30 studies suggest that unsaturated Mo centers, present on MoSx edge sites, are primarily responsible for the catalysis. Hence, amorphous MoSx materials that present high surface densities of catalytically active edge sites often outperform their crystalline counterparts.18-20,31 Unfortunately, amorphization both complicates identification and impedes structural characterization of catalytic sites. To resolve structural ambiguity, attain site-uniformity, and thereby facilitate the experimental and computational interrogation of heterogeneous catalysts and their atomic-scale activity, chemists have begun exploring metal–organic frameworks (MOFs) as catalyst supports. The nodes of many MOFs consist of oxy-metal clusters that can be viewed as structurally well-defined analogues of conventional metal-oxide supports (e.g., alumina, zirconia, or ceria powders). The most useful of these clusters present readily accessible catalyst-grafting points, defined by reactive or displaceable –OH and/or –OH2 ligands.32-40 The porous, crystalline nature of MOFs facilitates: a) exact assessments of the proton topology of candidate grafting points,32-33,40-41 b) sterically uniform access to these points and to installed catalysts, and c) uniformity of Lewis and

ACS Paragon Plus Environment

ACS Catalysis 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Brønsted acidities and/or basicities, for given structural types (e.g., bridging oxo/hydroxo, or terminal hydroxo/aqua sites).42 The uniformity or near-uniformity and periodicity of supporting nodes, as well as installed catalysts, facilitates experimental reproducibility and greatly enhances the practical value of computational modeling. Indeed, these kinds of systems lend themselves well to predictive reactivity and to mechanistically instructive experimental assessment of structureactivity relationships.32-37 Notably, the three-dimensional nature of hierarchically porous MOFs typically yields large internal surface areas and large densities of reactant-accessible catalysts.43 At the same time, the organic spacers constituting linkers between nodes can be remarkably effective for inhibiting or altogether preventing sintering or agglomeration of installed catalysts.44-45 These properties can be used to advantage to obtain structurally characterizable and catalytically active metallic, metal-oxide or -chalcogenide clusters, as well as monometallic species.46-47 Our preferred approach to system assembly has been postsynthetic installation of catalysts, either through vapor-phase atomic layer deposition in MOFs (AIM),35-37,48 or through solution-phase solvothermal deposition in MOFs (SIM).3235,39,49 The resulting arrays have proven catalytically competent for gas-phase chemical and for condensed-phase chemical, photochemical, and electrochemical transformations. Herein, we present the functionalization of a Zr6-based MOF, NU-1000 (Figure 1), with MoSx through SIM (the material here onwards named MoSx-SIM, with the precise structure of the resting state of the catalyst active-site being established by X-ray methods). The MOF consists of eight-connected Zr6(µ3– nodes and 1,3,5,8-(pO)4(µ3–OH)4(H2O)4(OH)4 benzoate)pyrene linkers (TBAPy4-).40-41 The framework offers multiple catalyst-grafting points41 and high hydrothermal stability,50 making it a potentially attractive heterogeneous support for electrocatalytic HER. Notably, the MOF can be readily physically integrated, in thin-film form, with conductive glass or other electrode materials – either before or after catalyst installation. The highly insulating nature of this MOF,51 however, prevents the vast majority of installed candidate electrocatalysts from being probed electrochemically or contributing to desired electrochemical transformations.14,52-53

Figure 1. Non-functionalized metal–organic framework, NU1000 with corresponding structures of the inorganic nodes and the organic linkers. Three of the twelve reactive hydroxo/aqua groups on the node are highlighted.

To overcome this problem, and enable delivery of reducing equivalents from the electrode to locations throughout the

Page 2 of 11

catalyst-functionalized MOF, we have introduced and examined three molecular-scale redox mediators (RMs, Scheme 1). Consistent with the hierarchical porosity of NU-1000 (Figure 1),40 we find that dissolved RMs can readily diffuse to the MOF-affixed MoSx catalysts and substantially boost overall catalytic rates. As described and discussed below, we additionally find that by changing the concentrations of RMs, by varying their formal potential (Ef), and by altering the amount of substrate (here via changes in pH), substantial mechanistic insight, specifically for aqueous electrocatalytic generation of H2 by heterogenized MoSx species, can be obtained.

Scheme 1. Proposed Redox Mediator-Assisted Electrocatalytic System with Various Redox Mediators

RESULTS AND DISCUSSION Characterization of MoSx-SIM: The MoSx-SIM synthesis is detailed in the Methods section. N2 adsorption-desorption isotherms reveal small decreases in Brunauer–Emmett–Teller (BET) surface area (from 2100 to 1700 m2/g) upon installation of molybdenum oxide clusters, the material named MoOxSIM synthesized according to the reported procedure,34 and to 1400 m2/g following sulfidation to MoSx-SIM (Figure 2a). (The differences in BET surface area are considerably lessened if account is taken of increases in sample mass upon introduction of molybdenum and upon replacement of selected oxygen-containing ligands by their sulfur-containing analogues and areas are expressed in volumetric (m2/cm3) rather than gravimetric (m2/g) form; see Table S1.) The preservation of the type IVc feature in the isotherm of MoSx-SIM implies that the MOF’s mesoporosity is retained. Indeed, density functional theory (DFT)-calculated pore size distributions (based on experimental N2 isotherm measurements) indicate only a modest decrease in mesopore size, i.e. from 30 to 27 Å (Figure S1). MoSx installation decreases the total pore volume by about 20% – more-or-less what one would anticipate based on scaling of pore volume with the square of the channel diameter; see Table S1. Thermogravimetric analysis (TGA) profiles of MoOx-SIM and MoSx-SIM, assessed at 600 °C, revealed additional residual masses of 27 and 32 wt %, respectively, relative to NU-1000. The added masses are attributed to the deposited MoOx or MoSx species that remain within the decomposed framework (Figure S2). The 5 wt % difference between MoOx-SIM and MoSx-SIM suggests that after treatment with H2S, partial replacement of oxo, hydroxo, and/or aqua ligands by sulfur-containing ligands (e.g. sulfide, disul-

ACS Paragon Plus Environment

Page 3 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

Figure 2. a) N2 isotherms at 77 K of NU-1000, MoOx-SIM, and MoSx-SIM. Filled and unfilled data points represent adsorption and desorption, respectively. The persistent step near P/P0 = 0.25 points to retention of the characteristic mesoporosity of NU-1000. b) The MoSx-anchored Zr6 node of the representative structure of MoSx-SIM showing all the crystallographically distinct MoSx units (Mo1-03) and the framework viewed along the c) c-axis (the minor structure Mo3 is omitted for clarity) and the d) a-axis. ure S5b).21 The ligation of sulfur(2-) by the Mo center is also fide, and/or HS–) has occurred. Scanning electron microscopy (SEM) images show that no changes in MOF crystallite morevident from the extended X-ray absorption fine structure phology40 accompany MoSx installation, and energy dispersive (EXAFS) spectrum (Figure S6). Namely, the phaseuncorrected peaks at 1.4 and 1.9 Å are attributable to Mo–O spectroscopy (EDS) line scanning shows a uniform distribution of Mo after H2S treatment (Fig ure S3a). As the Mo Land Mo–S single scattering events, respectively. The relatively low intensity of the features at 2.2–3 Å, attributed to a combiedge and the S K-edge energies significantly overlap, inductively-coupled plasma optical emission spectroscopy (ICPnation of multiple scattering events (e.g. Mo–O–Zr, Mo–O– OES) was used to determine the average Mo and S loadings, Mo, and/or Mo–S–Mo), suggests the presence of mononuclear MoSx or perhaps the MoSx clusters composed of only a few with the results being 2.6 ± 0.2 Mo/Zr6 and 2.5 ± 0.2 S/Mo; Mo atoms within the framework, similar to what has been the Mo loading is in agreement with the previously reported observed for MoOx-SIM.34 value for MoOx-SIM.34 As evidenced by EDS line scans (Figure S3d), we find that H2S exposure of unfunctionalized NUTo further elucidate the structure, we turned to single-crystal 1000 results in no detectable sulfur incorporation. X-ray diffraction measurements. (See the SI for measurement and refinement details.) While almost all of the MoOx units in Diffuse-reflectance infrared Fourier transform (DRIFT) specMoOx-SIM reside, in monometallic form, within the ~10 Å tra reveal that the intensities of peaks associated with the micropore between the two nodes parallel to the crystalMOF-node-based terminal/bridging –OH stretches at 3674 cm1 lographic c-axis (here onwards named "c-pores"; ca. 90% of for both MoOx-SIM and MoSx-SIM are much weaker than MoOx is estimated to reside within the c-pores; see Figures that of NU-1000 (Figure S4).40 This observation suggests that S7-10 and Table S2),55 treatment with H2S results in migration Mo remains grafted to the oxy-Zr6 node even after treatment of roughly two-thirds of the installed Mo to oxy ligands diwith H2S. Mo 3d X-ray photoelectron spectroscopy (XPS, rected toward the 31 Å mesopore, i.e., one bridging oxo and Figure S5a) revealed a 2.4 eV negative shift in binding energy one bridging hydroxo per isolated molybdenum ion, as shown with H2S treatment, implying molybdenum reduction from in Figure 2. The Mo sites directed toward mesopores are disoxidation-state VI to IV, and replacement of at least some oxo ordered over three positions, two of which are crystallographor hydroxo ligands with electron-rich sulfur-containing ligands. ically independent (labeled Mo1 and Mo3 in Figure 2b and in We note the similar shift in binding energy observed upon the SI). In contrast, the fraction Mo ions remaining in the ctreating thin layers of MoO3 with H2S.54 The 3d5/2 binding pores are disordered between two equivalent positions (Mo2 energy of MoSx-SIM is ~1 eV more positive to that of the 20 in Figure 2d). Like those in the mesopores, however, the mobulk MoS2, which may be attributed to coordination of Mo to lybdenum units are sited as monometallic complexes, rather the highly electron-deficient Zr6 node. S 2p XPS shows that than multi-molybdenum clusters. nearly all sulfur in the structure is in oxidation state 2– (Fig-

ACS Paragon Plus Environment

ACS Catalysis 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Electrocatalytic HER: The monomeric MoIV(SH)2 units structurally resemble edge sites of MoS2, i.e. sites considered primarily responsible for molybdenum disulfide’s electrocatalytic activity.21-22,26-28 Thus, the electrocatalytic HER performance of MoSx-SIM was examined; for experimental details, see the Methods section. As illustrated in Figure S11, cyclic voltammograms (CVs) of all samples in aqueous H2SO4 at pH = 1.2 show at the scan reversal point of –520 mV vs. RHE (reversible hydrogen electrode),56 current densities (j) that are small and that are “more similar than different.” Brief bulk electrolysis measurements (5 min. at –420 mV vs. RHE) yielded no detectable H2 with MOF films of either MoOx-SIM or NU-1000 during a bulk electrolysis experiment. In contrast, MoSx-SIM produced small, but detectable, amounts of H2, corresponding to a turnover number (TON) of 1.8 ± 0.5 and a turnover frequency (TOF) of 0.36 ± 0.09 min-1, if all molybdenum atoms are included in the calculations. It is tempting to attribute the poor catalytic performance to the aforementioned insulating character of the supporting framework,51 and resulting electrochemical inaccessibility of all but a small fraction of MoSx sites, i.e. those sites located proximal to the underlying conductive-glass electrode (fluorine-doped tin oxide (FTO)). This explanation is consistent with the comparatively small sizes of the voltammetric peaks assignable in Figure S11 to MoOx-SIM and MoSx-SIM. Indeed, integration of the peaks implies electrochemical accessibility of only those molybdenum ions lying within about 20 nm of the underlying electrode, corresponding to only 0.3 % of the total amount of species (see Figure S12 for the film thickness measurement via SEM that was used to calculate the aforementioned values). The voltammograms point to very similar or identical formal potentials (Ef) for framework-supported MoOx and MoSx i.e. values of roughly −20 mV vs. RHE – more-or-less ideal if the ultimate goal is a reversible, heterogeneous, single-metal-atom electrocatalyst. With these baseline results in hand, we turned to methyl viologen (MV2+) as a candidate redox mediator; see Scheme 1. The Ef for the MV2+/Ÿ+ couple (–360 mV vs. RHE), suggested that electrochemically generated MVŸ+ would have sufficient reducing power to drive the putative HER catalyst, Mo(SH)2. The small size of the mediator relative to the pores of the MOF, together with the high aqueous solubility of the chloride salt of MV2+, suggested to us that it would effectively permeate NU-1000 and shuttle electrons to otherwise nonaddressable Mo(SH)2 units. The voltammograms in Figure 3 show that the viologen couple: a) is chemically reversible in aq. H2SO4, b) can readily reach the FTO electrode and accept electrons, but c) is ineffective on its own, or in combination with Nafion, NU-1000, or MoOx-SIM, in catalyzing hydrogen evolution. In striking contrast, in the presence of MoSx-SIM the mediator yields a classic electrocatalytic wave, i.e. a sigmoidal voltammogram lacking hysteresis and plateauing well above the peak current for reduction of the mediator in the absence of catalyst (Figure 3). The onset of the catalytic wave, relative to the Ef of the mediator, contains information about the rate of reaction of the mediator with the catalyst (node-attached MoS2) and/or the rate of the reaction of the catalyst with H+. The height of the wave is a measure of the maximum catalytic rate for a given set of conditions and is typically limited by the rate of one of the component reactions shown or implied in Scheme 1, but, under some condition, can be limited by the rate at which

Page 4 of 11

the mediator can diffuse between the electrode and catalyst. The overall shape of the wave can, to a first approximation, be understood in terms of the fraction of mediator present in the catalytically useful reduced form at the electrode/solution (porous MOF) interface at a given applied potential. Returning to Figure 3, the introduction of 10 mM MV2+ increases the catalytic current at –500 mV vs. RHE (i.e. near the top of the catalytic wave) nearly 60-fold, and decreases the apparent Tafel slope (slope of plot of applied potential vs. log (catalytic current)) from –200 mV/decade to –86 mV/decade (Figure S13a). Electrolysis for five minutes at a fixed potential of –420 mV vs. RHE yielded a TON of 38 ± 5 hydrogen molecules per MoS2 site (based on measurement of the amount of H2 generated), and a TOF of 8 ± 1 min-1; for all chronoamperometric curves recorded during the bulk electrolysis, see Figure S14).57

Figure 3. CVs of the bare FTO, Nafion-coated FTO (labeled "Nafion"), NU-1000, MoOx-SIM, and MoSx-SIM measured in pH 1.2 H2SO4 solution with 10 mM of MV2+. A CV of MoSxSIM without MV2+ is included for comparison.

Powder X-ray diffraction (PXRD) patterns of MoSx-SIM after H2S treatment, after drop casting, and after electrolysis all show that the framework has stayed intact (Figure S16). SEM images of the MoSx-SIM electrode before and after electrolysis show no morphological differences and the EDS line scans show that Mo remains uniformly distributed within the MOF (Figure S3b and c). Transmission electron microscopy (TEM) images of the sample after catalysis (Figure S17) offer no indication of nanoparticle formation (for example, via agglomeration of the Mo(SH)2 units). ICP-OES measurements of a digested sample of MoSx-SIM, following use as a catalyst for bulk electrolysis indicated loss of a small fraction of catalyst (0.6 ± 0.1 Mo/Zr6). Notably, bulk electrolysis of the solution into which Mo had leached generated no detectable H2. Mediator Formal Potential Effects on the Electrocatalysis: We reasoned that if the rate of electrocatalysis is limited by the rate of electron transfer (ET) from the mediator to the catalyst, the system response would be sensitive to the Ef of the RM2+/Ÿ+, as a more negative Ef would correspond to a greater thermodynamic driving force (ΔGRMET) for mediator-tocatalyst ET. We selected for screening diaquat dibromide (DQ2+; Ef = −289 mV vs. RHE) and triquat dibromide (TQ2+; Ef = −531 mV vs. RHE); their structures are shown in Scheme 1. (See Figure S18 for CVs of bare FTO, nafion-coated FTO, NU-1000, and MoOx-SIM in the presence of 10 mM DQ2+ or

ACS Paragon Plus Environment

Page 5 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

Figure 4. CVs of MoSx-SIM in pH 1.2 H2SO4 solution with either b) 10 or c) 2 mM of various redox mediators. d) Faradaic efficiency (FE) plotted against the concentration of redox mediators. e) log(TOF) vs. log([RM2+]) plots, of which the slopes indicate the reaction order.

TQ2+.) Shown in Figure 4a are slow-scan-rate voltammograms of MoSx-SIM (on FTO) without and with 10 mM of each of the three candidate mediators. (Similar plots are shown in Figure 4b for RM concentrations of 2 mM.) Although the measurements are constrained (especially for TQ2+) by scan reversal at –520 mV, it is clear that plateau currents for electrocatalytic waves increase in the order DQ2+ < MV2+ < TQ2+. This trend is consistent with the notion that the rate of electrocatalysis is governed, at least in part, by the driving force and the rate of reaction for ET from the mediator to the catalyst. Figure 4c presents Faradaic efficiencies (FE values) for H2 formation via electrolysis at a fixed potential of –420 mV vs. RHE. Ideally, the FE values should approach 100%, provided only negligible amounts of RMŸ+ escape from the MOF film. Particularly for DQ2+ at concentrations of 5 and 10 mM, however, the FE values fall well below 100%. They can be largely recovered, however, by changing the potential for electrolysis to – 349 mV (= Ef(DQ2+/Ÿ+) – 60 mV). We have been unable to identify cause for the diminished FE. Given the sensitivity to electrolysis potential and to mediator concentration, however, we speculate that at sufficiently high concentrations, DQŸ+ dimerizes58 and that the dimer is only marginally competent, if at all, for delivery of reducing equivalents to catalytic Mo(SH)2 sites. As shown in the SI, a further complication is the chemical instability of TQŸ+ in acidic solutions. As such,

further investigation of triquat as a redox mediator was not pursued. By varying the RM2+ concentration and measuring the TOF through bulk electrolysis – again at –420 mV vs. RHE – a reaction order of ~0.5 was observed for both MV2+ and DQ2+. (See Figure 4d; for CVs with various concentrations, see Figure S21.) The unexpected fractional reaction order rules out many mechanisms for mediated electrocatalysis. A potential mechanism that does yield a reaction order of 0.5 for RM2+ is summarized by Equations 1–6. The required conditions are that: 1) the first or second ET to MoSx participates in, or precedes, the rate-determining step (RDS),59 2) electrogenerated RMŸ+ exists in a steady-state, and 3) due to electrochemical reduction, the concentration of RM2+ within the pores of MoSx-SIM at any given time is significantly lower than the concentration in bulk solution. (For the mathematical derivations, see the SI.) Note that steps 3 and 4, in some form, conceivably are interchangeable (i.e. the initial reduction step occurs first, followed by proton uptake, and then reduction by a second electron); if so, the square-root dependence of overall catalytic rate on RM2+ concentration would still hold (note that this simultaneously introduces complex pH dependence, see the SI).

ACS Paragon Plus Environment

ACS Catalysis 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 11

Figure 5. a) CVs of MoSx-SIM with pH or pD = 1.2 electrolyte solutions with 10 mM MV2+, which was used to derive b) KIE between –300 to –520 mV vs. RHE. c) CVs of MoSx-SIM with various pH electrolytes with 10 mM MV2+. d) log(TOF) vs. –pH plot of which the slope indicates the reaction order relative to H+. k1 catalytic rate (Figure 5d). Namely, when the pH is below 1.5 (1) 2RM 2RM2+ + 2eand the potential of the working electrode is fixed at –300 mV k-1 vs. RHE, a zero-order dependence on [H+] is observed, sugk2 gesting that RM-assisted ET alone defines the RDS, in agree(2) RM + (Mo IV(SH)2)2+ RM2+ + (Mo III(SH)2)+ ment with the kinetic isotope results. When the pH is ink-2 k3 creased from 1.6 to 2 and 2.5, the reaction order in protons RM + (Mo III(SH)2)+ (3) RM2+ + MoII(SH)2 changes to ~1. At pH 1.6, the rates for ET and proton transfer k-3 are evidently equal, or nearly so. k4 II + IV + Extending the applied potential to –450 mV vs. RHE, and (4) Mo (SH)2 + H (Mo (Hydride)(SH)2) k-4 therefore generating greater amounts of MVŸ+ (and boosting k5 the rates of ET steps 2 and 3), we find that the dependence of (5) (Mo IV(Hydride)(SH)2)+ (Mo IV(S)(SH))+ + H2 the overall reaction rate on [H+] persists down to pH 1 (the k6 lowest pH examined). The relevant data points in Figure 5d (6) (Mo IV(S)(SH))+ + H+ (Mo IV(SH)2)2+ admittedly are somewhat scattered (R2 = 0.95); nevertheless, a reaction order in protons of close to 1 (i.e. ~0.77) is clearly As the aforementioned applied potential (−420 mV vs. RHE) evident. We conclude that proton transfer defines, or particiis 60 mV negative of the Ef(MV2+/Ÿ+), the applied potential for pates in, the overall reaction rate. (See the SI for the rate exthe DQ2+-assisted system was also set to 60 mV negative of pression in the two aforementioned scenarios.) Equation 4, the mediator’s Ef (i.e. −349 mV vs. RHE). (Note that at this the formation of a molybdenum-bound hydride from H+ and a potential, Faradaic efficiencies are reasonably close to 100%; reduced catalyst, is one possibility for the overall mixed rate see Figure 4c.) Figure 4d shows that while TOF values are control.60 lower at –349 mV than at –420 mV, the apparent reaction order is unaffected. This observation implies that the first or Foot-of-the-wave Analysis: The preceding discussion is fosecond ET, i.e., the reaction in Equation 2 or 3, still particicused on identifying rate-defining steps, but leaves unquantipates in, or precedes, the RDS.53 fied the magnitudes of rate constants and thus, intrinsic catalytic activity. The foot-of-the-wave analysis method (FOWA) pH Effects on the Electrocatalysis: Despite the observed developed by Savéant and co-workers61 has been widely apsquare-root dependence of the overall rate of electrocatalysis plied to both homogeneous and heterogeneous electrocatalytic on mediator concentration, it is conceivable that subsequent systems in order to understand the intrinsic catalytic perforproton-dependent steps function as the reaction bottleneck (see mance of a system, where the performance is typically quantithe SI for details). Comparison of electrocatalytic waves for fied in terms of a potential-independent second-order rate conMoSx-SIM in H2SO4/H2O versus D2SO4/D2O suggests particistant (k) and overpotential (η) dependent TOFs.62 To our pation of H+ in the overall reaction rate (Figure 5a), but only knowledge, the FOWA approach has yet to be applied to meat potentials more negative than about –300 mV vs. RHE. As diated electrocatalysis of hydrogen evolution in aqueous solushown in Figure 5b, the H/D kinetic isotope effect (KIE) vartions. One potential roadblock to using this analysis in an ies from 1 at –300 mV (implying lack of participation H+ in aqueous system is that cathodic peak current of the catalyst the RDS) to greater than 2 at –500 mV vs. RHE (now imply(ip0) in the absence of the substrate (i.e., H+ or H2O) cannot be ing H+ participation). Support for this conclusion is provided trivially determined.63 However, in the present system, at sufby CVs in Figure 5c, where the pH of the electrolyte is varied 2+ ficiently low pH and sufficiently low applied potential, profrom 1 to 2.5 while holding the concentration of MV in solutons are enlisted only after the RDS. Under these special contion constant. Significant differences in catalytic current, at ditions, the FOWA can be implemented by defining RMŸ+ and different pH values, are evident only at applied potentials reduced MoSx as the catalyst and substrate, respectively. To more negative than –300 mV vs. RHE. this end, we define ip0 to be the peak current for the reduction Bulk electrolysis experiments at fixed applied potentials of – of RM2+ in the absence of MoSx, i.e., when a non-catalytic 300 or −450 mV vs. RHE, in a series of solutions ranging in MoOx-SIM-coated FTO is used as the electrode (see the SI pH from 1 to 2.5, also show that the applied potential and the for detailed calculations). Values of the second-order rate consolution proton activity act together to determine the signifistant, kRM, reflective of the kinetics of ET (i.e. Equations 2 cance (or not) of proton uptake in defining the overall electro-

ACS Paragon Plus Environment

Page 7 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

Figure 6. The foot-of-the-wave analysis of MoSx-SIM with 10 mM of MV2+ or DQ2+ in the pH 1.2 H2SO4 electrolyte. a) !!

!

! ! !!

vs.

! 1 + exp 𝐸 − 𝐸!" plots, of which the slopes indicate the value of the second order rate constant (k), and b) !!/∙! !" log(TOF)-η plots. achieve with a homogeneous monometallic catalyst or with a and 3), determined from the FOWA, are shown in Figure 6a non-porous heterogeneous catalyst. Notably, the Ef for the and Table 1. The larger value for kRM for the MV2+/Ÿ+ is consistent with the aforementioned order of the reducing power. MOF-grafted catalyst, like that for the metal-sulfur active-site Further using kRM yields log(TOF) versus η plot (where η is of hydrogenase, is nearly coincident with the formal potential for the hydrogen couple. In electrode-supported thin-film form, the overpotential for HER, Figure 6b). Intriguingly, the ordering of intrinsic TOF values at η = 0 mV (denoted as TOF0) the functionalized MOF is weakly electrocatalytic for hydrogen evolution from aqueous acid. The initially marginal elecis the opposite of that for kRM. Namely, the TOF0 in the trocatalytic activity of the material denoted MoSx-SIM, is DQ2+/Ÿ+-mediated system is more than an order of magnitude greater than that in the MV2+/Ÿ+-mediated system (Table 1), largely due to the insulating nature of the electrode-supported thin-film MOF–which leaves the overwhelming majority of behavior that is also reflected in the less negative onset potential for mediated electrocatalysis (Figures 4a-b). At low or the potentially catalytic sites electrochemically isolated from zero overpotential for the HER itself, the fraction of mediator the underlying electrode. The introduction of viologen-type redox mediators in the solution phase is followed by their present in reduced form (i.e., reactive form, RMŸ+) is considerpermeation of the functionalized MOF and diffusion through ably less for methyl viologen than for diquat – enoughless to more than offset the greater driving force and larger kRM value the framework to an underlying electrode poised to reduce the mediators. The reduced mediators then deliver electrons in for reaction of MVŸ+ versus with DQŸ+ with the MOF-noderepetitive fashion to otherwise isolated catalytic sites, thereby supported molybdenum catalyst. engendering sizeable catalytic currents and readily detectable Table 1. Second-order rate constants (kRM) and intrinsic production of molecular hydrogen. Mechanistic studies show TOF at zero overpotential (TOF0) of methyl viologen and that overall control of the rate of catalysis can be defined by diquat-assisted electrocatalysis estimated through the footmediator-to-catalyst ET, solution-to-catalyst proton transfer, of-the-wave analysis or both. We suggest that RM-assisted electrocatalysis can be profitably extended to other reactions and other catalysts supk TOF RM 0 Ef vs. RHE RM ported on highly porous, but problematically electrically insu(mV) (M–1 s–1) (×10–3 min–1) lating MOFs, and thus can further contribute to atomically precise understandings of mechanisms and of structureDQ2+ –289 16 ± 1 24 ± 2 activity relationships for other electrocatalytic systems. 2+

MV

–360

28 ± 2

1.9 ± 0.2

METHODS CONCLUSIONS The water-stable Zr6-based MOF, NU-1000, is amenable to functionalization with molybdenum and sulfur, elements that yield a catalyst that is functional for hydrogen evolution from acidified water. Single-crystal X-ray diffraction and related measurements show the catalyst to be a monometallic complex of molybdenum grafted to the hexa-zirconium(IV) node of the MOF by a bridging oxo ligand and a bridging hydroxo ligand, with the molybdenum center also coordinating, in its resting state, a pair of non-bridging SH– ligands. Siting the singlemetal-atom catalyst on the MOF permits it to be presented to the infilling aqueous solution at a concentration of more than 100 mM, in a zone just a few microns thick, exclusively proximal to the electrode – a desirable but challenging scenario to

MoSx-SIM Synthesis: MoOx-SIM was synthesized according to the reported procedure.34,64 Approximately 50 mg of the freshly activated MoOx-SIM was mounted onto a custom-made stainless steel sample holder, further mounted onto a sample holder chamber in an ALD reactor (Savannah S100, Cambridge Nanotech Inc.). The sample chamber was set to 125 °C at all times. Prior to any H2S pulsing, the sample was left in a N2 flow at approximately 0.6 Torr for 30 min. Subsequently, the stop-valve was closed and H2S was pulsed for 0.015 s followed by a 120 s exposure time to allow diffusion of the H2S molecules throughout the MoOx-SIM crystallites. The stop-valve was then opened and the sample chamber was purged under vacuum for 120 s afterwards to remove excess H2S. This process was repeated for 150 cycles (or 3 cycles/mg) and the resulting sample was stored in an argon-filled glove box. Electrocatalytic HER: All MOF-based electrodes were synthesized through simple drop casting of the sample suspended in a Nafion/ethanol solution. Briefly, 0.3 mL of Nafion solution was added to

ACS Paragon Plus Environment

ACS Catalysis 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

2.7 mL of ethanol and the mixture was dried using MgSO4. After filtration of the MgSO4 powder, 0.5 mL of the solution was used to disperse 6 mg of a MOF sample. In the case of MoSx-SIM, this was conducted using a N2-connected standard Schlenk line and a 2–5 mL Biotage microwave vial to minimize the exposure to oxygen. The mixture was sonicated for 5–10 min to form a suspension, and 6 µL of the suspension was dropcasted onto a FTO conducting glass with an exposed area of 0.25 cm2. The electrochemically active exposed area was controlled with insulating polyamide tape. The electrodes were left to dry and were used as the working electrode in a standard H-cell setup with a septum seal for collecting H2 gas. For MoSx-SIM-based electrodes, the drop casting and electrochemical cell fabrication were all done in an argon-filled glove box. Prior to any electrochemical measurements, the electrolyte (with or without RMs; for the triquat synthesis and characterization, see the SI) was purged with N2 gas for 20 min. CVs were taken with a 25 mV/s scan rate with a platinum coil and a Ag/AgCl electrode (3M) as the counter and reference electrodes, respectively (see the SI for details). Bulk electrolysis experiments were performed for 5 min and 200 µL of the headspace were injected into a GC for analysis of headspace composition. TOF was estimated by determining the total Mo loading via a UV-visible spectrum of a decomposed MoSx-SIM (for details in GC or TOF calculation, see the SI).

AUTHOR INFORMATION Corresponding Author *[email protected] *[email protected]

Author Contributions All authors have given approval to the final version of the manuscript.

Notes The authors declare no financial interest.

ASSOCIATED CONTENT The Supporting Information is available free of charge via the Internet at ACS Publication Website. Detailed Material Synthesis and Characterization (PDF) Crystallographic information file of MoOx-SIM (cif) Crystallographic information file of MoSx-SIM (cif)

ACKNOWLEDGMENTS This work was supported by the ANSER Center, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Award Number DE-SC0001059. H. N. gratefully acknowledges support from the Ryan Fellowship program of the Northwestern University International Institute of Nanotechnology. C.-W. K. acknowledges support from the Postdoctoral Research Abroad Program (105-2917-I-564-046) sponsored by the Ministry of Science and Technology (Taiwan). We thank Prof. Michael R. Wasielewski and his lab for the use of a gas chromatograph. This work made use of the IMSERC, J. B. Cohen X-Ray Diffraction, EPIC, and KECK II facilities of Northwestern University’s NUANCE Center, which has received support from the Soft and Hybrid Nanotechnology Experimental (SHyNE) Resource (NSF ECCS-1542205); the MRSEC program (NSF DMR-1720139) at the Materials Research Center; the International Institute for Nanotechnology (IIN); the Keck Foundation; and the State of Illinois, through the IIN. Use of the Advanced Photon Source is supported by the U.S. Department of Energy, Office of Science, and Office of Basic Energy Sciences, under Contract DE-AC02-06CH11357. Materials Research Collaborative Access Team (MRCAT, Sectors 10ID)

Page 8 of 11

operations are supported by the Department of Energy and the MRCAT member institutions.

REFERENCES (1) Timur, K.; E., S. M.; Anatoliy, S.; Malte, B.; Robert, S., The Haber–Bosch Process Revisited: On the Real Structure and Stability of “Ammonia Iron” under Working Conditions, Angew. Chem. Int. Ed. 2013, 52, 12723-12726. (2) Erisman, J. W.; Sutton, M. A.; Galloway, J.; Klimont, Z.; Winiwarter, W., How a Century of Ammonia Synthesis Changed the World, Nat. Geosci. 2008, 1, 636-639. (3) Jahangiri, H.; Bennett, J.; Mahjoubi, P.; Wilson, K.; Gu, S., A Review of Advanced Catalyst Development for Fischer-Tropsch Synthesis of Hydrocarbons from Biomass Derived Syn-Gas, Catal. Sci. Technol. 2014, 4, 2210-2229. (4) Khodakov, A. Y.; Chu, W.; Fongarland, P., Advances in the Development of Novel Cobalt Fischer−Tropsch Catalysts for Synthesis of Long-Chain Hydrocarbons and Clean Fuels, Chem. Rev. 2007, 107, 1692-1744. (5) Sarantaridis, D.; Atkinson, A., Redox Cycling of Ni-Based Solid Oxide Fuel Cell Anodes: A Review, Fuel Cells 2007, 7, 246258. (6) Marković, N. M.; Schmidt, T. J.; Stamenković, V.; Ross, P. N., Oxygen Reduction Reaction on Pt and Pt Bimetallic Surfaces: A Selective Review, Fuel Cells 2001, 1, 105-116. (7) Ogden, J. M., Prospects for Building a Hydrogen Energy Infastructure, Annu. Rev. Energy Env. 1999, 24, 227-279. (8) Vesborg, P. C. K.; Seger, B.; Chorkendorff, I., Recent Development in Hydrogen Evolution Reaction Catalysts and Their Practical Implementation, J. Phys. Chem. Lett. 2015, 6, 951-957. (9) Zeng, M.; Li, Y., Recent Advances in Heterogeneous Electrocatalysts for the Hydrogen Evolution Reaction, J. Mater. Chem. A 2015, 3, 14942-14962. (10) Ding, Q.; Song, B.; Xu, P.; Jin, S., Efficient Electrocatalytic and Photoelectrochemical Hydrogen Generation Using MoS2 and Related Compounds, Chem 2016, 1, 699-726. (11) McKone, J. R.; Marinescu, S. C.; Brunschwig, B. S.; Winkler, J. R.; Gray, H. B., Earth-Abundant Hydrogen Evolution Electrocatalysts, Chem. Sci. 2014, 5, 865-878. (12) Wang, X.; Li, W.; Xiong, D.; Liu, L., Fast Fabrication of Self-Supported Porous Nickel Phosphide Foam for Efficient, Durable Oxygen Evolution and Overall Water Splitting, J. Mater. Chem. A 2016, 4, 5639-5646. (13) Yin, H.; Zhao, S.; Zhao, K.; Muqsit, A.; Tang, H.; Chang, L.; Zhao, H.; Gao, Y.; Tang, Z., Ultrathin Platinum Nanowires Grown on Single-Layered Nickel Hydroxide with High Hydrogen Evolution Activity, Nat. Commun. 2015, 6, 6430. (14) Ramohlola, K. E.; Masikini, M.; Mdluli, S. B.; Monama, G. R.; Hato, M. J.; Molapo, K. M.; Iwuoha, E. I.; Modibane, K. D., Electrocatalytic Hydrogen Evolution Reaction of Metal Organic Frameworks Decorated with Poly (3-Aminobenzoic Acid), Electrochim. Acta 2017, 246, 1174-1182. (15) Khaselev, O.; Turner, J. A., A Monolithic PhotovoltaicPhotoelectrochemical Device for Hydrogen Production Via Water Splitting, Science 1998, 280, 425-427. (16) Esposito, D. V.; Hunt, S. T.; Stottlemyer, A. L.; Dobson, K. D.; McCandless, B. E.; Birkmire, R. W.; Chen, J. G., Low-Cost Hydrogen-Evolution Catalysts Based on Monolayer Platinum on Tungsten Monocarbide Substrates, Angew. Chem. Int. Ed. 2010, 49, 9859-9862. (17) Cheng, N.; Stambula, S.; Wang, D.; Banis, M. N.; Liu, J.; Riese, A.; Xiao, B.; Li, R.; Sham, T.-K.; Liu, L.-M.; Botton, G. A.; Sun, X., Platinum Single-Atom and Cluster Catalysis of the Hydrogen Evolution Reaction, Nat. Commun. 2016, 7, 13638. (18) Li, Y.; Wang, H.; Xie, L.; Liang, Y.; Hong, G.; Dai, H., MoS2 Nanoparticles Grown on Graphene: An Advanced Catalyst for the Hydrogen Evolution Reaction, J. Am. Chem. Soc. 2011, 133, 7296-7299. (19) Benck, J. D.; Chen, Z.; Kuritzky, L. Y.; Forman, A. J.; Jaramillo, T. F., Amorphous Molybdenum Sulfide Catalysts for

ACS Paragon Plus Environment

Page 9 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

Electrochemical Hydrogen Production: Insights into the Origin of Their Catalytic Activity, ACS Catal. 2012, 2, 1916-1923. (20) Tran, P. D.; Tran, Thu V.; Orio, M.; Torelli, S.; Truong, Q. D.; Nayuki, K.; Sasaki, Y.; Chiam, Sing Y.; Yi, R.; Honma, I.; Barber, J.; Artero, V., Coordination Polymer Structure and Revisited Hydrogen Evolution Catalytic Mechanism for Amorphous Molybdenum sulfide, Nat. Mater. 2016, 15, 640. (21) Li, H.; Tsai, C.; Koh, A. L.; Cai, L.; Contryman, A. W.; Fragapane, A. H.; Zhao, J.; Han, H. S.; Manoharan, H. C.; AbildPedersen, F.; Nørskov, J. K.; Zheng, X., Activating and Optimizing MoS2 Basal Planes for Hydrogen Evolution through the Formation of Strained Sulphur Vacancies, Nat. Mater. 2015, 15, 48. (22) Tsai, C.; Li, H.; Park, S.; Park, J.; Han, H. S.; Nørskov, J. K.; Zheng, X.; Abild-Pedersen, F., Electrochemical Generation of Sulfur Vacancies in the Basal Plane of MoS2 for Hydrogen Evolution, Nat. Commun. 2017, 8, 1-8. (23) Benck, J. D.; Hellstern, T. R.; Kibsgaard, J.; Chakthranont, P.; Jaramillo, T. F., Catalyzing the Hydrogen Evolution Reaction (HER) with Molybdenum Sulfide Nanomaterials, ACS Catal. 2014, 4, 3957-3971. (24) Seh, Z. W.; Kibsgaard, J.; Dickens, C. F.; Chorkendorff, I.; Nørskov, J. K.; Jaramillo, T. F., Combining Theory and Experiment in Electrocatalysis: Insights into Materials Design, Science 2017, 355, 1-12. (25) McCrory, C. C. L.; Jung, S.; Ferrer, I. M.; Chatman, S. M.; Peters, J. C.; Jaramillo, T. F., Benchmarking Hydrogen Evolving Reaction and Oxygen Evolving Reaction Electrocatalysts for Solar Water Splitting Devices, J. Am. Chem. Soc. 2015, 137, 4347-4357. (26) Ahn, H. S.; Bard, A. J., Electrochemical Surface Interrogation of a MoS2 Hydrogen-Evolving Catalyst: In Situ Determination of the Surface Hydride Coverage and the Hydrogen Evolution Kinetics, J. Phys. Chem. Lett. 2016, 7, 2748-2752. (27) Tang, Q.; Jiang, D.-E., Mechanism of Hydrogen Evolution Reaction on 1t-MoS2 from First Principles, ACS Catal. 2016, 6, 49534961. (28) Deng, Y.; Ting, L. R. L.; Neo, P. H. L.; Zhang, Y.-J.; Peterson, A. A.; Yeo, B. S., Operando Raman Spectroscopy of Amorphous Molybdenum Sulfide (MoSx) During the Electrochemical Hydrogen Evolution Reaction: Identification of Sulfur Atoms as Catalytically Active Sites for H+ Reduction, ACS Catal. 2016, 6, 7790-7798. (29) Tsai, C.; Chan, K.; Abild-Pedersen, F.; Norskov, J. K., Active Edge Sites in MoSe2 and WSe2 Catalysts for the Hydrogen Evolution Reaction: A Density Functional Study, Phys. Chem. Chem. Phys. 2014, 16, 13156-13164. (30) Tsai, C.; Abild-Pedersen, F.; Nørskov, J. K., Tuning the Mos2 Edge-Site Activity for Hydrogen Evolution Via Support Interactions, Nano Lett. 2014, 14, 1381-1387. (31) Wang, H.; Lu, Z.; Kong, D.; Sun, J.; Hymel, T. M.; Cui, Y., Electrochemical Tuning of Mos2 Nanoparticles on ThreeDimensional Substrate for Efficient Hydrogen Evolution, ACS Nano 2014, 8, 4940-4947. (32) Ji, P.; Manna, K.; Lin, Z.; Feng, X.; Urban, A.; Song, Y.; Lin, W., Single-Site Cobalt Catalysts at New Zr12(µ3-O)8(µ3-OH)8(µ2OH)6 Metal–Organic Framework Nodes for Highly Active Hydrogenation of Nitroarenes, Nitriles, and Isocyanides, J. Am. Chem. Soc. 2017, 139, 7004-7011. (33) Ji, P.; Manna, K.; Lin, Z.; Urban, A.; Greene, F. X.; Lan, G.; Lin, W., Single-Site Cobalt Catalysts at New Zr8(µ2-O)8(µ2-OH)4 Metal-Organic Framework Nodes for Highly Active Hydrogenation of Alkenes, Imines, Carbonyls, and Heterocycles, J. Am. Chem. Soc. 2016, 138, 12234-12242. (34) Noh, H.; Cui, Y.; Peters, A. W.; Pahls, D. R.; Ortuño, M. A.; Vermeulen, N. A.; Cramer, C. J.; Gagliardi, L.; Hupp, J. T.; Farha, O. K., An Exceptionally Stable Metal–Organic Framework Supported Molybdenum(VI) Oxide Catalyst for Cyclohexene Epoxidation, J. Am. Chem. Soc. 2016, 138, 14720-14726. (35) Ahn, S.; Thornburg, N. E.; Li, Z.; Wang, T. C.; Gallington, L. C.; Chapman, K. W.; Notestein, J. M.; Hupp, J. T.; Farha, O. K., Stable Metal–Organic Framework-Supported Niobium Catalysts, Inorg. Chem. 2016, 55, 11954-11961.

(36) Li, Z.; Peters, A. W.; Platero-Prats, A. E.; Liu, J.; Kung, C. W.; Noh, H.; DeStefano, M. R.; Schweitzer, N. M.; Chapman, K. W.; Hupp, J. T.; Farha, O. K., Fine-Tuning the Activity of Metal-Organic Framework-Supported Cobalt Catalysts for the Oxidative Dehydrogenation of Propane, J. Am. Chem. Soc. 2017, 139, 1525115258. (37) Peters, A. W.; Li, Z.; Farha, O. K.; Hupp, J. T., Toward Inexpensive Photocatalytic Hydrogen Evolution: A Nickel Sulfide Catalyst Supported on a High-Stability Metal–Organic Framework, ACS Appl. Mater. Interfaces 2016, 8, 20675-20681. (38) Jiao, L.; Wang, Y.; Jiang, H.-L.; Xu, Q., Metal–Organic Frameworks as Platforms for Catalytic Applications, Adv. Mater. 2017, 1703663, 1-23. (39) Gutov, O. V.; Hevia, M. G.; Escudero-Adán, E. C.; Shafir, A., Metal–Organic Framework (MOF) Defects under Control: Insights into the Missing Linker Sites and Their Implication in the Reactivity of Zirconium-Based Frameworks, Inorg. Chem. 2015, 54, 8396-8400. (40) Mondloch, J. E.; Bury, W.; Fairen-Jimenez, D.; Kwon, S.; DeMarco, E. J.; Weston, M. H.; Sarjeant, A. A.; Nguyen, S. T.; Stair, P. C.; Snurr, R. Q.; Farha, O. K.; Hupp, J. T., Vapor-Phase Metalation by Atomic Layer Deposition in a Metal–Organic Framework, J. Am. Chem. Soc. 2013, 135, 10294-10297. (41) Planas, N.; Mondloch, J. E.; Tussupbayev, S.; Borycz, J.; Gagliardi, L.; Hupp, J. T.; Farha, O. K.; Cramer, C. J., Defining the Proton Topology of the Zr6-Based Metal–Organic Framework Nu1000, J. Phys. Chem. Lett. 2014, 5, 3716-3723. (42) Klet, R. C.; Liu, Y.; Wang, T. C.; Hupp, J. T.; Farha, O. K., Evaluation of Bronsted Acidity and Proton Topology in Zr- and HfBased Metal-Organic Frameworks Using Potentiometric Acid-Base Titration, J. Mater. Chem. A 2016, 4, 1479-1485. (43) Islamoglu, T.; Goswami, S.; Li, Z.; Howarth, A. J.; Farha, O. K.; Hupp, J. T., Postsynthetic Tuning of Metal–Organic Frameworks for Targeted Applications, Acc. Chem. Res. 2017, 50, 805-813. (44) Kim, I. S.; Li, Z.; Zheng, J.; Platero-Prats, A. E.; Mavrandonakis, A.; Pellizzeri, S.; Ferrandon, M.; Vjunov, A.; Gallington, L. C.; Webber, T. E.; Vermeulen, N. A.; Penn, R. L.; Getman, R. B.; Cramer, C. J.; Chapman, K. W.; Camaioni, D. M.; Fulton, J. L.; Lercher, J. A.; Farha, O. K.; Hupp, J. T.; Martinson, A. B. F., Sinter-Resistant Platinum Catalyst Supported by Metal–Organic Framework, Angew. Chem. Int. Ed. 2018, 57, 909-913. (45) Li, Z.; Schweitzer, N. M.; League, A. B.; Bernales, V.; Peters, A. W.; Getsoian, A. B.; Wang, T. C.; Miller, J. T.; Vjunov, A.; Fulton, J. L.; Lercher, J. A.; Cramer, C. J.; Gagliardi, L.; Hupp, J. T.; Farha, O. K., Sintering-Resistant Single-Site Nickel Catalyst Supported by Metal–Organic Framework, J. Am. Chem. Soc. 2016, 138, 1977-1982. (46) Zhao, S. L.; Wang, Y.; Dong, J. C.; He, C. T.; Yin, H. J.; An, P. F.; Zhao, K.; Zhang, X. F.; Gao, C.; Zhang, L. J.; Lv, J. W.; Wang, J. X.; Zhang, J. Q.; Khattak, A. M.; Khan, N. A.; Wei, Z. X.; Zhang, J.; Liu, S. Q.; Zhao, H. J.; Tang, Z. Y., Ultrathin MetalOrganic Framework Nanosheets for Electrocatalytic Oxygen Evolution, Nat. Energy 2016, 1, 1-10. (47) Zhao, S.; Yin, H.; Du, L.; He, L.; Zhao, K.; Chang, L.; Yin, G.; Zhao, H.; Liu, S.; Tang, Z., Carbonized Nanoscale Metal–Organic Frameworks as High Performance Electrocatalyst for Oxygen Reduction Reaction, ACS Nano 2014, 8, 12660-12668. (48) Li, Z.; Peters, A. W.; Bernales, V.; Ortuño, M. A.; Schweitzer, N. M.; DeStefano, M. R.; Gallington, L. C.; Platero-Prats, A. E.; Chapman, K. W.; Cramer, C. J.; Gagliardi, L.; Hupp, J. T.; Farha, O. K., Metal–Organic Framework Supported Cobalt Catalysts for the Oxidative Dehydrogenation of Propane at Low Temperature, ACS Cent. Sci. 2017, 3, 31-38. (49) Cui, Y.; Rimoldi, M.; Platero-Prats, A. E.; Chapman, K. W.; Hupp, J. T.; Farha, O. K., Stabilizing a Vanadium Oxide Catalyst by Supporting on a Metal-Organic Framework, ChemCatChem 2018, 10, 1772-1777. (50) Howarth, A. J.; Liu, Y. Y.; Li, P.; Li, Z. Y.; Wang, T. C.; Hupp, J.; Farha, O. K., Chemical, Thermal and Mechanical Stabilities of Metal-Organic Frameworks, Nat. Rev. Mater. 2016, 1, 1-15.

ACS Paragon Plus Environment

ACS Catalysis 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(51) Wang, T. C.; Hod, I.; Audu, C. O.; Vermeulen, N. A.; Nguyen, S. T.; Farha, O. K.; Hupp, J. T., Rendering High Surface Area, Mesoporous Metal–Organic Frameworks Electronically Conductive, ACS Appl. Mater. Interfaces 2017, 9, 12584-12591. (52) Kung, C.-W.; Mondloch, J. E.; Wang, T. C.; Bury, W.; Hoffeditz, W.; Klahr, B. M.; Klet, R. C.; Pellin, M. J.; Farha, O. K.; Hupp, J. T., Metal–Organic Framework Thin Films as Platforms for Atomic Layer Deposition of Cobalt Ions to Enable Electrocatalytic Water Oxidation, ACS Appl. Mater. Interfaces 2015, 7, 28223-28230. (53) Kung, C.-W.; Audu, C. O.; Peters, A. W.; Noh, H.; Farha, O. K.; Hupp, J. T., Copper Nanoparticles Installed in Metal–Organic Framework Thin Films Are Electrocatalytically Competent for CO2 Reduction, ACS Energy Lett. 2017, 2394-2401. (54) Spevack, P. A.; McIntyre, N. S., A Raman and XPS Investigation of Supported Molybdenum Oxide Thin Films. 2. Reactions with Hydrogen Sulfide, J. Phys. Chem. 1993, 97, 1103111036. (55) We acknowledge the difference in the structures computationally predicted in ref. 34 and that of experimentally determined via single crystal X-ray diffraction measurement presented in this work. (56) Scans were limited to –520 mV vs. RHE in order to avoid potential complications due to electrochemical reduction of MOF linkers or, for diaquat, over-reduction of the mediator to its diradical (i.e., charge = 0) form. (57) As evident in Figure S15, extending the electrolysis time to 30 min resulted in significant loss in electrocatalytic activity. The mechanism of catalyst deactivation has yet to be determined. Nevertheless, to avoid complications due to deactivation, including potentially faulty conclusions regarding electrocatalysis mechanism, electrolysis times thereafter were limited to 5 minutes. (58) Braterman, P. S.; Song, J. I., Spectroelectrochemistry of Aromatic Ligands and Their Derivatives. 1. Reduction Products of 4,4'-Bipyridine, 2,2'-Bipyridine, 2,2'-Bipyrimidine, and Some Quaternized Derivatives, J. Org. Chem. 1991, 56, 4678-4682.

Page 10 of 11

(59) Here, we are loosely using the term, rate-determining step. As the later discussion reveals, multiple chemical steps can lead to overall reaction bottleneck, depending on the applied potential. (60) We note that sulfhydryl and hydride moieties in both amorphous and crystalline MoSx have been detected. See ref. 20, 2628. (61) Costentin, C.; Drouet, S.; Robert, M.; Savéant, J.-M., Turnover Numbers, Turnover Frequencies, and Overpotential in Molecular Catalysis of Electrochemical Reactions. Cyclic Voltammetry and Preparative-Scale Electrolysis, J. Am. Chem. Soc. 2012, 134, 11235-11242. (62) Hod, I.; Sampson, M. D.; Deria, P.; Kubiak, C. P.; Farha, O. K.; Hupp, J. T., Fe-Porphyrin-Based Metal–Organic Framework Films as High-Surface Concentration, Heterogeneous Catalysts for Electrochemical Reduction of CO2, ACS Catal. 2015, 5, 6302-6309. (63) Elgrishi, N.; Chambers, M. B.; Fontecave, M., Turning It Off! Disfavouring Hydrogen Evolution to Enhance Selectivity for Co Production During Homogeneous CO2 Reduction by CobaltTerpyridine Complexes, Chem. Sci. 2015, 6, 2522-2531. (64) On the basis of several studies (most of which have yet to be published), we find that single-crystal X-ray structures for metalion-functionalized MOFs can be more easily obtained by installing metals slowly in the presence of an inert solvent such as dry heptane (i.e. SIM), than by installing them via AIM, a comparatively rapid, vapor-phase method. As a rule, high-porosity, solvent-evacuated MOFs are poor thermal conductors and AIM/SIM of metal ions is appreciably exothermic. As such, thermal disorder at the metalinstallation stage is much more likely to accompany AIM of a vacuum-filled MOF than slow SIM of a solvent-filled MOF – where even a poorly thermally conductive solvent can function as a heat sink, thereby moderating the temperature increases accompanying MOF modification. Additionally, although not explored here, we find SIM to be more easily scalable than AIM, at least for AIM with our commercial ALD reactors.

ACS Paragon Plus Environment

Page 11 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

ACS Paragon Plus Environment

11