Evaluation of the H2 Evolving Activity of Benzenehexathiolate

as a viable strategy for improving catalytic performance and durability while maintaining a high degree of control over the structure and properties o...
0 downloads 6 Views 819KB Size
Subscriber access provided by UNIV OF MISSOURI ST LOUIS

Article

Evaluation of the H2 Evolving Activity of Benzenehexathiolate Coordination Frameworks and the Effect of Film Thickness on H2 Production Courtney A. Downes, Andrew J. Clough, Keying Chen, Joseph W. Yoo, and Smaranda C. Marinescu ACS Appl. Mater. Interfaces, Just Accepted Manuscript • DOI: 10.1021/acsami.7b15969 • Publication Date (Web): 18 Dec 2017 Downloaded from http://pubs.acs.org on December 19, 2017

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 free 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 accessible to all readers and 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.

ACS Applied Materials & Interfaces 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 10 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 Applied Materials & Interfaces

Evaluation of the H2 Evolving Activity of Benzenehexathiolate Coordination Frameworks and the Effect of Film Thickness on H2 Production Courtney A. Downes, Andrew J. Clough, Keying Chen, Joseph W. Yoo, and Smaranda C. Marinescu* Department of Chemistry, University of Southern California, Los Angeles, CA, 90089, USA KEYWORDS Metal dithiolenes, metal-organic frameworks, electrocatalysis, coordination polymers, solar energy conversion, hydrogen evolution

ABSTRACT: The design of earth-abundant catalysts for the electrochemical production of H2 from water is important for the realization of a sustainable energy future. Incorporation of molecular catalysts into extended frameworks has emerged as a viable strategy for improving catalytic performance and durability while maintaining a high degree of control over the structure and properties of the catalytic active site. Here, we investigate benzenehexathiolate (BHT) coordination frameworks as electrocatalysts for the hydrogen evolution reaction (HER) in pH 1.3 aqueous solutions. The electrocatalytic HER activity of BHT based coordination frameworks follows the order of CoBHT > NiBHT > FeBHT. CoBHT operates at an overpotential of 185 mV, the lowest observed overpotential of the reported metal dithiolene-based metal organic frameworks and coordination polymers to date. To further understand the properties that dictate electrocatalytic activity, the effect of film thickness on the HER performance of CoBHT, a parameter that has not been extensively explored for electrocatalytic coordination frameworks, was examined. As the thickness was increased to ~1 µm, charge and proton transfer through CoBHT was hindered, the number of electrochemically accessible active sites decreased, and the mechanical robustness of the modified electrode was diminished. The observed thickness-dependent HER activity of CoBHT highlights the importance of practical electrode construction and offers insight into how to optimize proton and electron transfer properties and active site densities within coordination frameworks without reducing the mechanical robustness of the immobilized catalysts.

INTRODUCTION To sustainably meet the rising energy demand associated with rapid global population and economic growth, continued development and deployment of renewable energy resources is paramount.1 The unparalleled abundance of solar energy, while attractive to meet the massive global energy demand, is hindered by its intermittent nature, resulting in a mismatch between energy supply and demand.2 Sustainable hydrogen production from water splitting (2H2O → 2H2 + O2) has emerged as a promising pathway for the storage and conversion of solar energy.2-3 Hydrogen is a valuable energy carrier that can be transformed into electricity using fuel cell technology or used in the production of industrially relevant chemicals such as ammonia and methanol.4-5 Hydrogen is currently generated from steam methane reforming whereby four molecules of H2 are produced per CH4; however, one molecule of CO2 is also produced.6 A sustainable energy future will require the carbon-neutral production of H2 and water splitting is a viable pathway to do so. Electrocatalysts that are costeffective and scalable must be developed for the hydrogen evolution (HER) and oxygen evolution reactions (OER) to facilitate water splitting at high efficiencies with minimal energy input. Earth-abundant homogeneous7-8 and heterogeneous9-11 catalysts have been extensively explored as alternatives to

platinum, the benchmark electrocatalyst for the HER. Inherent problems with homogeneous catalysts such as the lack of solubility and stability in aqueous media and the diffusion (to the electrode surface) dependent activation of the solubilized catalyst limit their practical feasibility. Similarly, heterogeneous catalysts have ill-defined active sites making mechanistic understanding, rational design, and optimization difficult. These short-comings can be overcome by combining the advantageous properties of homogeneous and heterogeneous catalysts. The heterogenization of molecular catalysts has been accomplished through a variety of methods12-13 such as surface adsorption,14-18 electropolymerization,19-21 covalent attachment,2224 and incorporation into extended frameworks.25-31 These systems can display the improved stability and robustness associated with heterogeneous catalysts while the welldefined catalytic active sites reminiscent of molecular systems can be extensively studied and modified to give insight into the mechanism of catalysis. Molecular catalysts incorporated into extended frameworks such as coordination polymers and metal-organic frameworks (MOFs) have been investigated as electrocatalysts for energy converting applications.32-34 The use of MOFs as intrinsic electrocatalysts, however, has been limited because of their traditionally insulating nature. The facilitation of efficient charge and electron transfer is an

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

important property for the promotion of high electrocatalytic activity. The ability to deposit thin film coordination frameworks directly onto electrode surfaces has improved their electrocatalytic performance because thin film MOFs exhibit higher conductivities than the bulk materials.35-36 However, most of these systems still exhibit insufficient charge transfer properties inhibiting high electrocatalytic activity. Additionally, the weak coordination bonds in the frameworks lead to structural instability in aqueous conditions. The development of conductive coordination frameworks has been a breakthrough in the application of MOFs for electrocatalytic energy conversion.37 The use of redox active ligands has facilitated charge transport through the frameworks and stronger metal-ligand bonds have improved stability in aqueous acidic and/or alkaline media. The few successful electrocatalytic MOFs displaying high activity and stability for the HER,28-30 oxygen reduction reaction (ORR),31 oxygen evolution reaction (OER),38-41 and CO2 reduction (CO2RR)42-43 all exhibit intrinsic electrical conductivity. The conductive nature of metal dithiolene based coordination frameworks makes these systems attractive candidates for electrocatalytic applications.36, 44-46 We have extensively explored the use of trinucleating and dinucleating dithiolene ligand scaffolds to synthesize 2D and 1D metal dithiolene coordination frameworks28 and polymers,25-26 which perform as robust electrocatalysts for the HER. To further investigate the variables that dictate the HER activity of dithiolene-based coordination frameworks, we examine here the effect of the coordinated metal by performing detailed electrochemical analysis of CoBHT, NiBHT, and FeBHT (BHT = benzenhexathiolate, Figure 1). FeBHT was synthesized through an interfacial reaction previously reported in the literature for the isolation of BHT-based coordination frameworks.28, 44 The synthesis and conductivity of NiBHT was reported by Nishihara and coworkers,44-45 however, the electrocatalytic HER activity was not studied. We have previously explored the electrocatalytic HER performance of CoBHT with the system displaying an overpotential of 340 mV and Faradaic efficiency of 97% in pH 1.3 aqueous solutions and a Tafel slope of 108 mV/dec (pH 2.6).28 A modified synthetic method was utilized here for the isolation of CoBHT and electrochemical analysis was performed using glassy carbon electrodes (GCE). Electrochemical techniques such as electrochemical impedance spectroscopy (EIS) and double-layer capacitance (Cdl) measurements, which provide insight into the charge transfer/HER kinetics and an estimation of the electrochemically active surface area, respectively, were not employed in the previous report on CoBHT. This analysis is used here to improve the fundamental understanding of the parameters that dictate the observed HER activity of CoBHT. Additionally, we explore the effect of film thickness on the HER activity of CoBHT. For layered materials, thicker films correspond to higher bulk catalyst loadings (higher number of available active sites) leading to improved HER performance.47-48 However, it has been observed that a plateau can be reached whereby increases in catalyst loading do not generate higher electrocatalytic current densities.18, 42, 47, 49 Poor diffusion of electrons and protons

Page 2 of 10

through thick films limits the number of electrochemically accessible active sites even as bulk catalyst loading increases, thus inhibiting the HER activity.47, 49-50 The longterm durability of thicker films can also be problematic. High mass loadings can lead to catalyst modified electrodes exhibiting limited mechanical robustness with cracking and delamination from the electrode surface occurring during electrochemical testing.

Figure 1. Structure of the metal benzenehexathiolate (BHT) coordination frameworks (CoBHT, NiBHT, and FeBHT) studied here.

Previous investigation into the HER performance of monolayer and bulk powders of cobalt dithiolene and mixed dithiolene-diamine coordination frameworks revealed increases in the overpotential to achieve 10 mA/cm2 and the Tafel slope as the method of catalyst deposition was changed from a monolayer (0.8 ± 0.1 nm thickness) to bulk powder (0.12 mg catalyst prepared as an ink).30 These differences are expected as thin film dithiolene- and diamine- based coordination frameworks are more conductive than the bulk materials because of a reduction in the number of grain boundaries, which inhibit charge transfer.35-36, 44-45 Further investigation into the effect of film thickness on the HER activity by increasing the amount of deposited catalyst from a monolayer to multilayer films was not performed. The influence of film thickness on the activity of intrinsic MOF electrocatalysts, although an important parameter for optimizing catalytic performance, has not been extensively explored. To the best of our knowledge, the most comprehensive study on the effect of MOF film thickness on catalytic activity was performed by Yaghi, Yang, and coworkers on a cobalt porphyrin based MOF, Al2(OH)2TCCP-Co (TCPP-H2 = 4,4’,4”,4”(porphyrin5,10,15,20-tetrayl)tetrabenzoate).42 The CO2 reduction performance increased with the number of atomic layer deposition (ALD) cycles before reaching maximum activity at 50 cycles (MOF film thickness of ~30-70 nm). At higher loadings, the performance decreased because of poor charge and mass transport properties. Extensive analysis of the influence of film thickness on the HER performance of coordination frameworks is necessary to facilitate the continued development of these systems as practical catalysts for energy converting devices. By varying the film thickness from ~20-1000 nm, we look to understand how efficiently charges and protons move through CoBHT. We can then determine the optimal

ACS Paragon Plus Environment

Page 3 of 10 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 Applied Materials & Interfaces

conditions to maximize the number of accessible active sites and HER activity without increasing the resistance of the system and limiting electron and proton transfer.

RESULTS AND DISCUSSION Electrochemical Analysis of FeBHT and NiBHT. To investigate the electrocatalytic hydrogen evolving activity of CoBHT, NiBHT, and FeBHT, we performed detailed electrochemical measurements in a three-electrode configuration with N2-saturated pH 1.3 H2SO4 electrolyte, where all of the potentials are referenced to the reversible hydrogen electrode (RHE). Multiple electrodes modified with the BHT-based frameworks were prepared and electrochemically tested with the results presented in the Supporting Information. We have recently shown that the deposition method used here, immersion of the electrode through the interfacially grown film, leads to small variances in the bulk catalyst loading even for films of the same thickness.51 These small variances in bulk loading and the morphological differences the deposition method presented here induces for each sample (how the sample adheres to the substrate) can influence the measured Cdl and Rct values. Thus, it is important to establish the trends observed for the electrocatalytic HER activity are valid for multiple modified electrodes. The EIS response for all three systems studied was best described by the two-time constant serial model (2TS), where the high frequency semi-circle is related to the surface porosity of the electrode or the contact between the electrode and the catalyst layer and the low frequency semi-circle is related to the charge transfer kinetics of HER (Figure S1).48, 52-54 NiBHT and FeBHT were readily deposited on GCE following the interfacial synthesis and electrochemically explored as hydrogen evolving catalysts. FeBHT exhibits overpotentials ranging from 473-541 mV (Figure 2a, Table 1, Figure S3, Table S1). Increases in Cdl resulted in reduced HER activity for FeBHT, which was not expected (Figure S3-S5, Table S1). The rise in Cdl was accompanied by an increase in the Rct and Tafel slope. Enhancements in the Cdl (1.4 to 2.6 mF/cm2) corresponded to increases in the Rct from 102.8 to 221.3 Ω (η = 470 mV) and Tafel slope from 112 to 143 mV/dec (Figure S4-S6, Table S1). These results indicate poor charge and proton transfer for the FeBHT samples with higher catalyst loading. The reproducibility of the polarization curves was also limited for FeBHT. Successive scans generated smaller overpotentials and improved charge transfer properties. This could be correlated to the observed removal of FeBHT from the electrode surface during electrochemical testing as we have shown that FeBHT is more active for the HER at lower catalyst loadings. For NiBHT, overpotentials as low as 331 mV were achieved at a Cdl value of 6.29 mF/cm2 (Figure 2, Table 1, Figure S8, Table S2). A broad reduction feature is observed at ~0 V for NiBHT and is assigned to reduction of the framework. This feature has been previously reported for NiBHT and the related nickel benzene-1,2,4,5-tetrathiolate coordination polymer (Figure S9).26, 44 Higher electrocatalytic activity was associated with larger Cdl values and reduced Rct (Table S2-S3). The value of Rct was dependent on the applied overpotential as at more negative potentials

the charge transfer resistance was reduced indicating faster HER kinetics (Figure S10-S11, Table S3). The intensity of the low frequency feature in the Bode plot (Figure S11) was diminished and shifted to higher frequencies at larger overpotentials. The Tafel slope did not change with Cdl, but the exchange current density (j0) increased with more electrochemically accessible active sites (Figure S12).49 Controlled potential electrolysis (CPE) of NiBHT at 0.72 V revealed a steady current response for the first hour followed by a gradual decline in generated current (Figure S13). The current drop-off was a response to the delamination of NiBHT from the electrode surface as the black insoluble catalyst was seen suspended in solution and minimal NiBHT was observed on the GCE surface after 7 h of electrolysis. ICP-OES studies revealed no solubilized nickel species confirming the decline in activity was due to catalyst delamination not dissolution. Analysis of the gaseous products by gas chromatography after 1 h of CPE confirmed H2 production with the Faradaic efficiency (FE) ranging from 65-75%. Most nickel dithiolene based catalysts operate at low to moderate efficiencies (FE ~ 6080%)26, 55-56 in comparison to cobalt dithiolene catalysts, which perform at high efficiencies (FE > 90%).14, 25, 28, 57 Xray photoelectron spectroscopy (XPS) following CPE revealed nickel and sulfur features analogous to those of the as-prepared NiBHT (Figure S14-S15).44

Figure 2. (a) Polarization curves for CoBHT (pink), NiBHT (cyan), and FeBHT (blue); scan rate: 100 mV/s, (b) The current density difference (Δj = ja – jc) at 0.15 V vs RHE plotted against the scan rate (markers); the Cdl values are estimated through linear fitting of the plots (dashed lines). All measurements were carried out in N2-saturated pH 1.3 H2SO4 solutions at room temperature.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 4 of 10

Table 1. Summary of the Electrocatalytic Hydrogen Evolution Properties of Metal Benzenehexathiolate Coordination Frameworks Sample

η @ 10 mA/cm2 (mV)

Cdl (mF/cm2)

Rct (Ω)

Tafel slope (mV/dec)

j0 (A/cm2)

CoBHT

185

27.1

64.4 (η = 270 mV)

88

10-4.9

NiBHT

331

6.29

218.5 (η = 370 mV)

67

10-7.8

FeBHT

473

1.43

102.8 (η = 470 mV)

119

10-6.7

Comparison of the HER Activity of CoBHT, NiBHT, and FeBHT. Polarization curves of CoBHT, NiBHT, and FeBHT are presented in Figure 2a highlighting the influence of the metal center on the HER activity of benzenehexathiolate-based coordination frameworks. CoBHT displayed the lowest overpotential, 185 mV, in comparison to NiBHT (331 mV) and FeBHT (405 mV). The increased activity of CoBHT is correlated to a higher electrochemically accessible surface area (Figure 2b) and lower charge transfer resistance at smaller overpotentials (Table 1, Figure 3a). NiBHT and FeBHT exhibit very large Rct values at -0.27 V vs RHE inhibiting high electrocatalytic HER at low overpotentials (Figure S16-S17, Table S4), however, at larger overpotentials, the Rct values decreased significantly (Figure 2a, Figure S16-S17, Table S4). For NiBHT, Rct decreased from 5,316 (η = 270 mV) to 218.5 Ω (η = 370 mV), while for FeBHT, a comparable Rct to CoBHT and NiBHT was not achieved until η = 470 mV. The Rct for FeBHT decreased from 1.8 × 104 (η = 270 mV) to 102.8 Ω (η = 470 mV). CoBHT is the superior HER catalyst stemming from the significantly enhanced charge transfer through the material at modest overpotentials and greater number of electrochemically accessible active sites. Higher electrocatalytic HER activity was achieved for CoBHT at all film thicknesses investigated (~20-1000 nm) in comparison to NiBHT (131 nm) even as the number of electrochemically accessible active sites decreased for CoBHT (Table 2). The larger j0 of CoBHT at all thicknesses compared to NiBHT highlights the higher intrinsic electrocatalytic HER activity of CoBHT (Table 1, Table 2). To further improve our understanding of how the metal center influences the electrocatalytic HER activity, the conductivity of the BHT-based coordination frameworks is currently under investigation in our laboratory to gain a complete understanding of the origin of the distinct electrochemical responses observed for CoBHT, NiBHT, and FeBHT. Tafel analysis revealed Tafel slopes of 88, 67, and 119 mV/dec for CoBHT, NiBHT, and FeBHT, respectively (Figure 3b, Table 1). As with previously reported metal dithiolene-based coordination frameworks and polymers, the Tafel slopes indicate the Volmer discharge reaction is the rate-limiting step.25-26, 28-30 The Tafel slope of FeBHT, 119 mV/dec, is much larger than that of CoBHT and NiBHT. The electronic structure of the analogous molecular complex, [Fe(bdt)2]2⁻, could provide insight into the origins of the larger Tafel slope and greater barrier for Hads. The Fe 3d orbitals are much higher in energy than the ligand (benzenedithiolate) orbitals resulting in a predominately metal centered singly occupied molecular orbital (SOMO). The

electronic structure is best described as a spin triplet FeII coordinated to two innocent benzenedithiolate ligands, [FeII(bdt2⁻)2]2⁻.58-59

Figure 3. (a) EIS spectra measured at -0.27 V, -0.37 V, and 0.47 V vs RHE presented as Nyquist plots (markers) with respective fits (solid line) for CoBHT (pink), NiBHT (cyan), and FeBHT (blue), respectively. (b) Tafel analysis of CoBHT (pink), NiBHT (cyan), and FeBHT (blue). All measurements were carried out in N2-saturated pH 1.3 H2SO4 solutions at room temperature.

Computational studies propose the molecular complexes, [Co(bdt)2]⁻ and [Ni(bdt)2]⁻ (bdt = benzene-1,2dithiolate), undergo relatively similar HER mechanisms with the thiolate moieties of the ligand scaffold the most likely site for initial protonation.55, 60 The non-innocent nature of these complexes arises because the SOMO of [Co(bdt)2]⁻ contains signiVicant mixed-metal ligand character, while the SOMO for [Ni(bdt)2]⁻ is mostly ligandbased.58-59 The lack of non-innocence in the benzenedithiolate ligands and metal dominant SOMO of [Fe(bdt)2]2⁻

ACS Paragon Plus Environment

Page 5 of 10 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 Applied Materials & Interfaces

could give rise to a different mechanism for HER than the computational studies propose for [Co(bdt)2]⁻ and [Ni(bdt)2]⁻. The increase in effective nuclear charge as the metal is changed from Fe to Co to Ni for the benzenedithiolate complexes modulates the metal and ligand character in the molecular orbitals, potentially resulting in different HER mechanisms and rate-limiting steps for the frameworks. More comprehensive theoretical studies are necessary to further understand the differences in the measured Tafel slopes.

Evaluation of the Thickness-Dependent HER Activity of CoBHT. Looking to understand proton and charge transfer in CoBHT to optimize the HER performance, we varied the amount of BHT used in the synthesis of CoBHT, which resulted in films of different thicknesses. The average film thickness ranged from 23 nm to 1 μm (Figure S19S22). It has been shown that hydrogen production can be dependent on the film thickness because of changes to the electron, charge, and proton transfer properties and the number of available and accessible active sites.47-50, 61-62 The effect of film thickness on the electrochemical hydrogen evolving activity of CoBHT was investigated on glassy carbon electrodes. To distinguish the electrochemical results of films of different thicknesses, the electrodes studied are denoted as CoBHT-X where X corresponds to the film thickness in nanometers. Scanning electron microscopy (SEM) images revealed morphological differences for the CoBHT films of different thicknesses (Figure S23-S27). Films within the thickness range of 20-250 nm exhibited a sheet like morphology expected for thin films. Puckering and wrinkling of the film was observed, which can help to expose more active sites for HER. The 1000 nm films displayed an extensively cracked and rough morphology with puckered and overlapping flakes of film. This type of morphology has the potential to expose more active sites and the interior of the thick film, which can induce higher electrocatalytic activity. However, the cracked morphology gives rise to extensive grain boundaries that can potentially inhibit electron transfer from the electrode surface to the exposed active sites, thus reducing the number of electroactive catalytic sites. Polarization curves of CoBHT modified GCE with 5 different film thicknesses, CoBHT-23, CoBHT-57, CoBHT-157, CoBHT-244, and CoBHT-1000 are presented in Figure 4a. As the film thickness was increased from 23 to 244 nm, the overpotential to achieve 10 mA/cm2 decreased from 246 mV to 185 mV (Figure 4a, Table 2). The reduced overpotential was accompanied by an increase in the Cdl from 2.1 to 27.1 mF/cm2. Within the thickness range 23 to 244 nm, an increase in the amount of deposited CoBHT resulted in the expected increase in the number of accessible active sites (Figure 4b, Table 2). As the thickness of the films increased to 1000 nm, reduced HER activity was observed and an overpotential of 213 mV was measured. Concomitantly, the number of electrochemically accessible active sites did not improve as the film thickness was increased from 244 to 1000 nm indicating not all the cobalt sites were electrocatalytically active. The substantial increase in mass loading for CoBHT-1000 did not result in a greater number of electrochemically accessible active sites as compared to films in the thick-

ness range of 20-250 nm, suggesting insufficient electron migration from the electrode surface through the 1000 nm thick film. The inability to electrochemically access the greater number of cobalt active sites associated with the large increase in film thickness from 244 to 1000 nm inhibits the electrocatalytic HER activity of CoBHT-1000.

Figure 4. (a) Polarization curves for CoBHT-23 (purple), CoBHT-57 (blue), CoBHT-157 (orange), CoBHT-244 (red), and CoBHT-1000 (green); scan rate: 100 mV/s, (b) The current density difference (Δj = ja – jc) at 0.15 V vs RHE plotted against the scan rate (markers); the Cdl values are estimated through linear fitting of the plots (dashed lines). All measurements were carried out in N2-saturated pH 1.3 H2SO4 solutions at room temperature.

The charge transfer resistance, extrapolated from EIS experiments, can reveal differences in the ability to move charge through CoBHT as the film thickness is modulated. EIS analysis at -0.27 V vs RHE for CoBHT-23, CoBHT-57, CoBHT-157, CoBHT-244, and CoBHT-1000 is presented in Figure 5a. The Rct decreases from 116 to 64 Ω as the film thickness increases from 23 to 244 nm. The improved charge transfer kinetics for CoBHT-244 facilitates electrocatalytic hydrogen evolution at smaller overpotentials. The intensity of the major feature of the Bode plot was reduced and shifted to lower frequencies for the most active electrode, CoBHT-244, confirming the superior HER kinetics (Figure S33). Further increases in film thickness to 1000 nm (CoBHT-1000) led to diminished charge transfer properties and HER kinetics as evidenced by the larger Rct (144.6 Ω) and the greater intensity of the major feature of the Bode plot (Figure S33). The larger charge transfer resistance as the film thickness reached 1000 nm indicates the charge transfer from the electrode to the CoBHT sur-

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 10

face and finally to the protons in solution is faster for thinner films.50

Table 2. Summary of the Electrochemical Properties of CoBHT with Different Film Thicknesses for the Hydrogen Evolution Reaction Sample

Average Film Thickness (nm)

η @ 10 mA/cm2 (mV)

Cdl (mF/cm2)

Rct (Ω) @ η = 270 mV

Tafel slope (mV/dec)

j0 (A/cm2)

CoBHT-1000

1000

213

15.9

144.6

62

10-6.3

CoBHT-244

244

185

27.1

64.4

88

10-4.9

CoBHT-157

157

192

22.8

78.3

92

10-4.8

CoBHT-57

57

225

3.10

100.8

76

10-5.6

CoBHT-23

23

246

2.11

116.2

77

10-5.6

Figure 5. (a) EIS spectra measured at -0.27 V vs RHE presented as Nyquist plots (markers) with respective fits (solid line) for CoBHT-23 (purple), CoBHT-57 (blue), CoBHT-157 (orange), CoBHT-244 (red), and CoBHT-1000 (green) (b) Comparison of the current density (mA/cm2) versus thickness (nm) at -0.20 V vs RHE (blue), -0.22 V vs RHE (green), -0.25 V vs RHE (pink), and -0.27 V vs RHE (purple). All measurements were carried out in N2-saturated pH 1.3 H2SO4 solutions at room temperature.

CoBHT-1000 exhibited an overpotential superior to CoBHT-23 and CoBHT-57, however, the Rct indicates sluggish charge transfer properties as the HER reaction rate is enhanced at larger overpotentials. Analyzing the electrocatalytic current density versus the film thickness at four

potentials (-0.20, -0.22, -0.25 and -0.27 V vs RHE) revealed the HER activity increased until a thickness of 244 nm (Figure 5b). At 1000 nm, smaller current densities were generated at all potentials compared to CoBHT-157 and CoBHT-244. For CoBHT-1000, the higher catalyst loadings associated with increased film thickness did not result in larger current densities and lower overpotentials as compared to films in the thickness range of ~20-250 nm because of poor charge transfer properties and fewer electrochemically accessible active sites (Table 2). At small overpotentials, CoBHT-1000 was catalytically competent and generated higher current densities than CoBHT-57. However, at larger overpotentials (η = 270 mV), the HER activity of CoBHT-57 (25.1 mA/cm2) surpassed the activity of CoBHT-1000 (23.1 mA/cm2) confirming the results of the EIS experiments performed at this potential (Table 2, Figure 5). Sluggish electron and charge transfer and poor proton permeation through the thick film limited the catalytic performance at larger overpotentials (Figure 5).50 Since variances in bulk catalyst loading, which can lead to deviations in the measured Cdl and Rct values even for films of the same thickness, are expected because of the deposition utilized here (immersion of the electrode through the film), multiple CoBHT-modified electrodes were evaluated to assess the validity of the observed thickness-dependent HER activity. CoBHT-1000 modified electrodes reproducibly exhibited insufficient charge transfer and proton diffusion properties at large overpotentials limiting the HER activity (Figure S34). The larger electrochemically accessible surface areas achieved for CoBHT-1000 modified electrodes in comparison to CoBHT-23 and CoBHT-57 (Table 2, Table S5) initially resulted in higher HER performance (lower overpotentials to achieve 10 mA/cm2), however, this high activity could not be sustained at larger overpotentials as poor charge and proton diffusion rendered CoBHT-1000 catalytically incompetent. Therefore, an optimal film thickness of ~250 nm has been identified for CoBHT whereby proton diffusion, charge transport, and the number of electrochemically accessible active sites are maximized. The Tafel slopes were extrapolated from similar potential windows and ranged from 62-92 mV/dec for CoBHT

ACS Paragon Plus Environment

Page 7 of 10 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 Applied Materials & Interfaces

modified GCE (Figure S35, Table 2). The Tafel slopes correspond to a rate-liming Volmer discharge step. The cobalt dithiolene-based coordination frameworks and polymers previously reported as HER catalysts all displayed large Tafel slopes suggesting the reaction was limited by insufficient adsorption of H+.25, 28, 30 Similar Tafel slopes and exchange current densities (Table 2) were measured for samples with relatively analogous film thicknesses (CoBHT-23/CoBHT-57 and CoBHT-157/CoBHT-244). Identifying if different barriers for Hads arise as the thickness is modulated from ~20-1000 nm is currently under investigation. The examination into the role of film thickness on the HER activity for CoBHT revealed an optimal thickness range (~250 nm) where charge transport, proton transfer, and the number of electrochemically accessible active sites were maximized. We have shown that the diminished charge and proton transfer properties and reduced number of electrochemically accessible active sites (even at higher bulk catalyst loadings) led to the poorer than expected HER performance of CoBHT-1000. In addition to the fundamental properties such as charge and proton transfer that dictate the electrocatalytic activity, practical device construction is an important factor to consider when determining the optimal catalyst loading and deposition method. The mechanical robustness of the deposited catalyst was greatly decreased as the film thickness was increased to 1000 nm. Cracking of the film and delamination from the electrode surface during electrochemical testing was common for CoBHT-1000, thereby limiting the practical feasibility of electrodes constructed with such high mass loadings. Therefore, it is important to consider the mechanical robustness of the catalyst-modified electrodes, in addition to electrocatalytic performance, when selecting an appropriate deposition method and catalyst thickness.

CONCLUSIONS In summary, cobalt, nickel, and iron coordination frameworks based on benzenehexathiolate were synthesized via an interfacial reaction and electrochemically investigated as hydrogen evolving catalysts in pH 1.3 aqueous solutions. CoBHT exhibited the lowest overpotential of 185 mV. The ability to operate at lower overpotentials for CoBHT was related to the improved charge transfer properties in comparison to NiBHT and FeBHT. Understanding the effect of film thickness on electrocatalytic activity, an insufficiently explored area in coordination framework electrocatalysis, is vital for determining the scalability and practicality of these systems as catalysts for energy conversion applications. Our exploration into the movement of charge and protons through the CoBHT films revealed thickness dependent HER activity. As the film thickness was increased, the overpotential initially decreased as higher accessible active sites were realized. An optimal thickness was achieved at 244 nm where charge transport, proton diffusion and transfer, and the number of electrochemically accessible active sites were optimized. As the film thickness was increased further to ~1000 nm, the charge transfer and proton permeation through the thick films was insufficient to promote high electrocatalytic activity espe-

cially when the HER kinetics were enhanced at larger overpotentials. Although the cracked and rough morphology of CoBHT1000 allowed for the exposure of more active sites, the introduction of significant grain boundaries diminished the charge transfer properties and hindered the electrochemical accessibility of many of the available active sites. Diffusion of protons and electrons through the thick film limited the overall HER performance canceling out any potential gains associated with increased mass loading. The identification of an optimized thickness (~250 nm) is valuable for informing device construction as the higher mass loading films (~1000 nm) display minimal mechanical robustness limiting their practical viability. The optimized CoBHT modified GCE, presented here, operate at overpotentials ranging from 185-246 mV representing, to the best of our knowledge, the highest performing HER activity of metal dithiolene based coordination frameworks. Understanding the ability to move charge and protons through these systems has allowed for the design and construction of heterogenized molecular catalysts displaying HER activity more reminiscent of traditional heterogeneous catalysts.

ASSOCIATED CONTENT Supporting Information. Cyclic voltammetry data, EIS circuit, EIS fit data, double-layer capacitance measurements, Bode plots, electrolysis results, and XPS spectra. This material is available free of charge via the Internet at http://pubs.acs.org.

AUTHOR INFORMATION Corresponding Author * [email protected]

Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.

Notes The authors declare no competing financial interest.

ACKNOWLEDGEMENT We are grateful to the University of Southern California (USC) and USC Women in Science and Engineering (WiSE) for funding. We thank Buddhinie Jayathilake and Prof. Sri R. Narayan for assistance with ZSimpWin software. XPS and SEM data were collected at the Center for Electron Microscopy and Microanalysis, USC. AFM measurements were collected using Prof. Mark E. Thompson’s instrument.

ABBREVIATIONS BHT, benzenehexathiolate; GCE, glassy carbon electrode; GR, graphite electrode; CV, cyclic voltammetry; Cdl, double-layer capacitance; FE, faradaic efficiency; SEM, scanning electron microscopy; EIS, electrochemical impedance spectroscopy; XPS, X-ray photoelectron spectroscopy; AFM, atomic force microscopy

REFERENCES (1) Chu, S.; Majumdar, A., Opportunities and Challenges for a Sustainable Energy Future. Nature 2012, 488 (7411), 294-303.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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) Gray, H. B., Powering the Planet with Solar Fuel. Nat. Chem. 2009, 1 (1), 7. (3) Turner, J. A., Sustainable Hydrogen Production. Science 2004, 305 (5686), 972-974. (4) Kim, D.; Sakimoto, K. K.; Hong, D.; Yang, P., Artificial Photosynthesis for Sustainable Fuel and Chemical Production. Angew. Chem. Int. Ed. 2015, 54 (11), 3259-3266. (5) 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 (6321), eaad4998. (6) LeValley, T. L.; Richard, A. R.; Fan, M., The Progress in Water Gas Shift and Steam Reforming Hydrogen Production Technologies – A Review. Int. J. Hydrog. Energy 2014, 39 (30), 16983-17000. (7) Thoi, V. S.; Sun, Y.; Long, J. R.; Chang, C. J., Complexes of EarthAbundant Metals for Catalytic Electrochemical Hydrogen Generation under Aqueous Conditions. Chem. Soc. Rev. 2013, 42, 2388-2400. (8) McKone, J. R.; Marinescu, S. C.; Brunschwig, B. S.; Winkler, J. R.; Gray, H. B., Earth-Abundant Hydrogen Evolution Electrocatalysts. Chem. Sci. 2014, 5 (3), 865-878. (9) Faber, M. S.; Jin, S., Earth-Abundant Inorganic Electrocatalysts and their Nanostructures for Energy Conversion Applications. Energy Environ. Sci. 2014, 7 (11), 3519-3542. (10) Roger, I.; Shipman, M. A.; Symes, M. D., Earth-Abundant Catalysts for Electrochemical and Photoelectrochemical Water Splitting. Nat. Rev. Chem. 2017, 1, 0003. (11) Anantharaj, S.; Ede, S. R.; Sakthikumar, K.; Karthick, K.; Mishra, S.; Kundu, S., Recent Trends and Perspectives in Electrochemical Water Splitting with an Emphasis on Sulfide, Selenide, and Phosphide Catalysts of Fe, Co, and Ni: A Review. ACS Catal. 2016, 6 (12), 8069-8097. (12) Bullock, R. M.; Das, A. K.; Appel, A. M., Surface Immobilization of Molecular Electrocatalysts for Energy Conversion. Chem. Eur. J. 2017, 23 (32), 7626-7641. (13) Queyriaux, N.; Kaeffer, N.; Morozan, A.; Chavarot-Kerlidou, M.; Artero, V., Molecular Cathode and Photocathode Materials for Hydrogen Evolution in Photoelectrochemical Devices. J Photochem. Photobiol. C: Photochem. Rev. 2015, 25, 90-105. (14) Eady, S. C.; MacInnes, M. M.; Lehnert, N., Immobilized Cobalt Bis(benzenedithiolate) Complexes: Exceptionally Active Heterogeneous Electrocatalysts for Dihydrogen Production from Mildly Acidic Aqueous Solutions. Inorg. Chem. 2017, 56 (19), 11654-11667. (15) Eady, S. C.; MacInnes, M. M.; Lehnert, N., A Smorgasbord of Carbon: Electrochemical Analysis of Cobalt– Bis(benzenedithiolate) Complex Adsorption and Electrocatalytic Activity on Diverse Graphitic Supports. ACS Appl. Mater. Interfaces 2016, 8 (36), 23624-23634. (16) Eady, S. C.; Peczonczyk, S. L.; Maldonado, S.; Lehnert, N., Facile Heterogenization of a Cobalt Catalyst via Graphene Adsorption: Robust and Versatile Dihydrogen Production Systems. Chem. Commun. 2014, 50 (59), 8065-8068. (17) Tran, P. D.; Le Goff, A.; Heidkamp, J.; Jousselme, B.; Guillet, N.; Palacin, S.; Dau, H.; Fontecave, M.; Artero, V., Noncovalent Modification of Carbon Nanotubes with Pyrene-Functionalized Nickel Complexes: Carbon Monoxide Tolerant Catalysts for Hydrogen Evolution and Uptake. Angew. Chem. Int. Ed. 2011, 50 (6), 1371-1374. (18) Zhang, W.; Haddad, A. Z.; Garabato, B. D.; Kozlowski, P. M.; Buchanan, R. M.; Grapperhaus, C. A., Translation of LigandCentered Hydrogen Evolution Reaction Activity and Mechanism of a Rhenium-Thiolate from Solution to Modified Electrodes: A Combined Experimental and Density Functional Theory Study. Inorg. Chem. 2017, 56 (4), 2177-2187. (19) Cabrera, C. R.; Abruña, H. D., Electrocatalysis of CO2 Reduction at Surface Modified Metallic and Semiconducting

Page 8 of 10

Electrodes. J Electroanal. Chem. Interfacial Electrochem. 1986, 209 (1), 101-107. (20) O'Toole, T. R.; Margerum, L. D.; Westmoreland, T. D.; Vining, W. J.; Murray, R. W.; Meyer, T. J., Electrocatalytic Reduction of CO2 at a Chemically Modified Electrode. J. Chem. Soc., Chem. Commun. 1985, (20), 1416-1417. (21) O'Toole, T. R.; Sullivan, B. P.; Bruce, M. R. M.; Margerum, L. D.; Murray, R. W.; Meyer, T. J., Electrocatalytic Reduction of CO2 by a Complex of Rhenium in Thin Polymeric Films. J Electroanal. Chem. Interfacial Electrochem. 1989, 259 (1), 217-239. (22) Andreiadis, E. S.; Jacques, P.-A.; Tran, P. D.; Leyris, A.; Chavarot-Kerlidou, M.; Jousselme, B.; Matheron, M.; Pécaut, J.; Palacin, S.; Fontecave, M.; Artero, V., Molecular Engineering of a Cobalt-Based Electrocatalytic Nanomaterial for H2 Evolution under Fully Aqueous Conditions. Nat. Chem. 2013, 5 (1), 48-53. (23) Le Goff, A.; Artero, V.; Jousselme, B.; Tran, P. D.; Guillet, N.; Metaye, R.; Fihri, A.; Palacin, S.; Fontecave, M., From Hydrogenases to Noble Metal-Free Catalytic Nanomaterials for H2 Production and Uptake. Science 2009, 326 (5958), 1384-1387. (24) Huan, T. N.; Jane, R. T.; Benayad, A.; Guetaz, L.; Tran, P. D.; Artero, V., Bio-Inspired Noble Metal-Free Nanomaterials Approaching Platinum Performances for H2 Evolution and Uptake. Energy Environ. Sci. 2016, 9 (3), 940-947. (25) Downes, C. A.; Marinescu, S. C., Efficient Electrochemical and Photoelectrochemical H2 Production from Water by a Cobalt Dithiolene One-Dimensional Metal-Organic Surface. J. Am. Chem. Soc. 2015, 137 (43), 13740-13743. (26) Downes, C. A.; Marinescu, S. C., One Dimensional Metal Dithiolene (M = Ni, Fe, Zn) Coordination Polymers for the Hydrogen Evolution Reaction. Dalton Trans. 2016, 45 (48), 19311-19321. (27) Downes, C. A.; Marinescu, S. C., Bioinspired Metal Selenolate Polymers with Tunable Mechanistic Pathways for Efficient H2 Evolution. ACS Catal. 2017, 7 (1), 848-854. (28) Clough, A. J.; Yoo, J. W.; Mecklenburg, M. H.; Marinescu, S. C., Two-Dimensional Metal-Organic Surfaces for Efficient Hydrogen Evolution from Water. J. Am. Chem. Soc. 2015, 137 (1), 118-121. (29) Dong, R.; Pfeffermann, M.; Liang, H.; Zheng, Z.; Zhu, X.; Zhang, J.; Feng, X., Large-Area, Free-Standing, Two-Dimensional Supramolecular Polymer Single-Layer Sheets for Highly Efficient Electrocatalytic Hydrogen Evolution. Angew. Chem. Int. Ed. 2015, 54 (41), 12058-12063. (30) Dong, R.; Zheng, Z.; Tranca, D. C.; Zhang, J.; Chandrasekhar, N.; Liu, S.; Zhuang, X.; Seifert, G.; Feng, X., Immobilizing Molecular Metal Dithiolene–Diamine Complexes on 2D Metal–Organic Frameworks for Electrocatalytic H2 Production. Chem. Eur. J. 2017, 23 (10), 2255-2260. (31) Miner, E. M.; Fukushima, T.; Sheberla, D.; Sun, L.; Surendranath, Y.; Dincă, M., Electrochemical Oxygen Reduction Catalysed by Ni3(hexaiminotriphenylene)2. Nat. Commun. 2016, 7, 10942. (32) Sakamoto, R.; Takada, K.; Pal, T.; Maeda, H.; Kambe, T.; Nishihara, H., Coordination Nanosheets (CONASHs): Strategies, Structures and Functions. Chem. Commun. 2017, 53 (43), 57815801. (33) Downes, C. A.; Marinescu, S. C., Electrocatalytic MetalOrganic Frameworks for Energy Applications. ChemSusChem 2017, 10 (22), 4374-4392. (34) Solomon, M. B.; Church, T. L.; D'Alessandro, D. M., Perspectives on Metal-Organic Frameworks with Intrinsic Electrocatalytic Activity. CrystEngComm 2017, 19 (29), 40494065. (35) Sheberla, D.; Sun, L.; Blood-Forsythe, M. A.; Er, S.; Wade, C. R.; Brozek, C. K.; Aspuru-Guzik, A.; Dincă, M., High Electrical Conductivity in Ni3(2,3,6,7,10,11-hexaiminotriphenylene)2, a Semiconducting Metal–Organic Graphene Analogue. J. Am. Chem. Soc. 2014, 136 (25), 8859-8862. (36) Clough, A. J.; Skelton, J. M.; Downes, C. A.; de la Rosa, A. A.; Yoo, J. W.; Walsh, A.; Melot, B. C.; Marinescu, S. C., Metallic

ACS Paragon Plus Environment

Page 9 of 10 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 Applied Materials & Interfaces

Conductivity in a Two-Dimensional Cobalt Dithiolene MetalOrganic Framework. J. Am. Chem. Soc. 2017, 139 (31), 1086310867. (37) Sun, L.; Campbell, M. G.; Dincă, M., Electrically Conductive Porous Metal–Organic Frameworks. Angew. Chem. Int. Ed. 2016, 55 (11), 3566-3579. (38) Lu, X.-F.; Liao, P.-Q.; Wang, J.-W.; Wu, J.-X.; Chen, X.-W.; He, C.-T.; Zhang, J.-P.; Li, G.-R.; Chen, X.-M., An Alkaline-Stable, Metal Hydroxide Mimicking Metal–Organic Framework for Efficient Electrocatalytic Oxygen Evolution. J. Am. Chem. Soc. 2016, 138 (27), 8336-8339. (39) Wang, L.; Wu, Y.; Cao, R.; Ren, L.; Chen, M.; Feng, X.; Zhou, J.; Wang, B., Fe/Ni Metal–Organic Frameworks and Their BinderFree Thin Films for Efficient Oxygen Evolution with Low Overpotential. ACS Appl. Mater. Interfaces 2016, 8 (26), 1673616743. (40) Zhao, S.; Wang, Y.; Dong, J.; He, C.-T.; Yin, H.; An, P.; Zhao, K.; Zhang, X.; Gao, C.; Zhang, L.; Lv, J.; Wang, J.; Zhang, J.; Khattak, A. M.; Khan, N. A.; Wei, Z.; Zhang, J.; Liu, S.; Zhao, H.; Tang, Z., Ultrathin Metal–Organic Framework Nanosheets for Electrocatalytic Oxygen Evolution. Nat. Energy 2016, 1, 16184. (41) Duan, J.; Chen, S.; Zhao, C., Ultrathin Metal-Organic Framework Array for Efficient Electrocatalytic Water Splitting. Nat. Commun. 2017, 8, 15341. (42) Kornienko, N.; Zhao, Y.; Kley, C. S.; Zhu, C.; Kim, D.; Lin, S.; Chang, C. J.; Yaghi, O. M.; Yang, P., Metal–Organic Frameworks for Electrocatalytic Reduction of Carbon Dioxide. J. Am. Chem. Soc. 2015, 137 (44), 14129-14135. (43) Lin, S.; Diercks, C. S.; Zhang, Y.-B.; Kornienko, N.; Nichols, E. M.; Zhao, Y.; Paris, A. R.; Kim, D.; Yang, P.; Yaghi, O. M.; Chang, C. J., Covalent Organic Frameworks Comprising Cobalt Porphyrins for Catalytic CO2 Reduction in Water. Science 2015, 349 (6253), 1208-1213. (44) Kambe, T.; Sakamoto, R.; Hoshiko, K.; Takada, K.; Miyachi, M.; Ryu, J. H.; Sasaki, S.; Kim, J.; Nakazato, K.; Takata, M.; Nishihara, H., π -Conjugated Nickel Bis(dithiolene) Complex Nanosheet. J. Am. Chem. Soc. 2013, 135 (7), 2462-2465. (45) Kambe, T.; Sakamoto, R.; Kusamoto, T.; Pal, T.; Fukui, N.; Hoshiko, K.; Shimojima, T.; Wang, Z.; Hirahara, T.; Ishizaka, K.; Hasegawa, S.; Liu, F.; Nishihara, H., Redox Control and High Conductivity of Nickel Bis(dithiolene) Complex π-Nanosheet: A Potential Organic Two-Dimensional Topological Insulator. J. Am. Chem. Soc. 2014, 136 (41), 14357-14360. (46) Huang, X.; Sheng, P.; Tu, Z.; Zhang, F.; Wang, J.; Geng, H.; Zou, Y.; Di, C.-a.; Yi, Y.; Sun, Y.; Xu, W.; Zhu, D., A Two-Dimensional π–d Conjugated Coordination Polymer with Extremely High Electrical Conductivity and Ambipolar Transport Behaviour. Nat. Commun. 2015, 6, 7408. (47) Ji, S.; Yang, Z.; Zhang, C.; Liu, Z.; Tjiu, W. W.; Phang, I. Y.; Zhang, Z.; Pan, J.; Liu, T., Exfoliated MoS2 Nanosheets as Efficient Catalysts for Electrochemical Hydrogen Evolution. Electrochim. Acta 2013, 109, 269-275. (48) Vrubel, H.; Moehl, T.; Gratzel, M.; Hu, X., Revealing and Accelerating Slow Electron Transport in Amorphous Molybdenum Sulphide Particles for Hydrogen Evolution Reaction. Chem. Commun. 2013, 49 (79), 8985-8987. (49) McAteer, D.; Gholamvand, Z.; McEvoy, N.; Harvey, A.; O’Malley, E.; Duesberg, G. S.; Coleman, J. N., Thickness Dependence and Percolation Scaling of Hydrogen Production Rate in MoS2

Nanosheet and Nanosheet–Carbon Nanotube Composite Catalytic Electrodes. ACS Nano 2016, 10 (1), 672-683. (50) Lin, Y.-K.; Chen, R.-S.; Chou, T.-C.; Lee, Y.-H.; Chen, Y.-F.; Chen, K.-H.; Chen, L.-C., Thickness-Dependent Binding Energy Shift in Few-Layer MoS2 Grown by Chemical Vapor Deposition. ACS Appl. Mater. Interfaces 2016, 8 (34), 22637-22646. (51) Downes, C. A.; Marinescu, S. C., Understanding Variability in the Hydrogen Evolution Activity of a Cobalt Anthracenetetrathiolate Coordination Polymer. ACS Catal. 2017, 7 (12), 8605-8612. (52) Navarro-Flores, E.; Chong, Z.; Omanovic, S., Characterization of Ni, NiMo, NiW and NiFe Electroactive Coatings as Electrocatalysts for Hydrogen Evolution in an Acidic Medium. J. Mol. Cat. A: Chem. 2005, 226 (2), 179-197. (53) Kucernak, A. R. J.; Naranammalpuram Sundaram, V. N., Nickel Phosphide: the Effect of Phosphorus Content on Hydrogen Evolution Activity and Corrosion Resistance in Acidic Medium. J. Mater. Chem. A 2014, 2 (41), 17435-17445. (54) Ma, L.; Ting, L. R. L.; Molinari, V.; Giordano, C.; Yeo, B. S., Efficient Hydrogen Evolution Reaction Catalyzed by Molybdenum Carbide and Molybdenum Nitride Nanocatalysts Synthesized via the Urea Glass Route. J. Mater. Chem. A 2015, 3 (16), 8361-8368. (55) Zarkadoulas, A.; Field, M. J.; Artero, V.; Mitsopoulou, C. A., Proton-Reduction Reaction Catalyzed by Homoleptic Nickel–bis1,2-dithiolate Complexes: Experimental and Theoretical Mechanistic Investigations. ChemCatChem 2017, 9 (12), 23082317. (56) Zarkadoulas, A.; Field, M. J.; Papatriantafyllopoulou, C.; Fize, J.; Artero, V.; Mitsopoulou, C. A., Experimental and Theoretical Insight into Electrocatalytic Hydrogen Evolution with Nickel Bis(aryldithiolene) Complexes as Catalysts. Inorg. Chem. 2016, 55 (2), 432-444. (57) McNamara, W. R.; Han, Z.; Alperin, P. J.; Brennessel, W. W.; Holland, P. L.; Eisenberg, R., A Cobalt-Dithiolene Complex for the Photocatalytic and Electrocatalytic Reduction of Protons. J. Am. Chem. Soc. 2011, 133 (39), 15368-15371. (58) Ray, K.; Begum, A.; Weyhermüller, T.; Piligkos, S.; van Slageren, J.; Neese, F.; Wieghardt, K., The Electronic Structure of the Isoelectronic, Square-Planar Complexes [FeII(L)2]2- and [CoIII(LBu)2]- (L2- and (LBu)2- = Benzene-1,2-dithiolates):  An Experimental and Density Functional Theoretical Study. J. Am. Chem. Soc. 2005, 127 (12), 4403-4415. (59) Ray, K.; Petrenko, T.; Wieghardt, K.; Neese, F., Joint Spectroscopic and Theoretical Investigations of Transition Metal Complexes Involving Non-Innocent Ligands. Dalton Trans. 2007, (16), 1552-1566. (60) Solis, B. H.; Hammes-Schiffer, S., Computational Study of Anomalous Reduction Potentials for Hydrogen Evolution Catalyzed by Cobalt Dithiolene Complexes. J. Am. Chem. Soc. 2012, 134 (37), 15253-15256. (61) Gholamvand, Z.; McAteer, D.; Harvey, A.; Backes, C.; Coleman, J. N., Electrochemical Applications of Two-Dimensional Nanosheets: The Effect of Nanosheet Length and Thickness. Chem. Mater. 2016, 28 (8), 2641-2651. (62) Yu, Y.; Huang, S.-Y.; Li, Y.; Steinmann, S. N.; Yang, W.; Cao, L., Layer-Dependent Electrocatalysis of MoS2 for Hydrogen Evolution. Nano Lett. 2014, 14 (2), 553-558.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 Paragon Plus Environment

Page 10 of 10