Brønsted Acid Scaling Relationships Enable Control Over Product

Jun 7, 2019 - Figure 1. Summary of representative macrocyclic cobalt ORR catalysts ..... to establish the catalytic rate law under both conditions (Fi...
0 downloads 0 Views 2MB Size
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.

Research Article Cite This: ACS Cent. Sci. 2019, 5, 1024−1034

http://pubs.acs.org/journal/acscii

Brønsted Acid Scaling Relationships Enable Control Over Product Selectivity from O2 Reduction with a Mononuclear Cobalt Porphyrin Catalyst Yu-Heng Wang,† Patrick E. Schneider,‡ Zachary K. Goldsmith,‡ Biswajit Mondal,† Sharon Hammes-Schiffer,*,‡ and Shannon S. Stahl*,† †

Department of Chemistry, University of WisconsinMadison, Madison, Wisconsin 53706, United States Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States

Downloaded via 91.243.90.106 on July 25, 2019 at 08:38:58 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.



S Supporting Information *

ABSTRACT: The selective reduction of O2, typically with the goal of forming H2O, represents a long-standing challenge in the field of catalysis. Macrocyclic transition-metal complexes, and cobalt porphyrins in particular, have been the focus of extensive study as catalysts for this reaction. Here, we show that the mononuclear Co-tetraarylporphyrin complex, Co(porOMe) (porOMe = meso-tetra(4-methoxyphenyl)porphyrin), catalyzes either 2e−/2H+ or 4e−/4H+ reduction of O2 with high selectivity simply by changing the identity of the Brønsted acid in dimethylformamide (DMF). The thermodynamic potentials for O2 reduction to H2O2 or H2O in DMF are determined and exhibit a Nernstian dependence on the acid pKa, while the CoIII/II redox potential is independent of the acid pKa. The reaction product, H2O or H2O2, is defined by the relationship between the thermodynamic potential for O2 reduction to H2O2 and the CoIII/II redox potential: selective H2O2 formation is observed when the CoIII/II potential is below the O2/H2O2 potential, while H2O formation is observed when the CoIII/II potential is above the O2/H2O2 potential. Mechanistic studies reveal that the reactions generating H2O2 and H2O exhibit different rate laws and catalyst resting states, and these differences are manifested as different slopes in linear free energy correlations between the log(rate) versus pKa and log(rate) versus effective overpotential for the reactions. This work shows how scaling relationships may be used to control product selectivity, and it provides a mechanistic basis for the pursuit of molecular catalysts that achieve low overpotential reduction of O2 to H2O.



potential of only 0.68 V (eqs 1 and 2).5 Nonetheless, generation of H2O2 is often kinetically favored because this reaction avoids the relatively high kinetic barrier associated with cleavage of the O−O bond.

INTRODUCTION The reduction of O2 to H2O, often called the oxygen reduction reaction (ORR), is among the most important chemical reactions on the planet due to its role in natural and artificial energy production.1−4 Extensive efforts have been directed toward understanding how this reaction is catalyzed by enzymes in nature and developing heterogeneous and molecular catalysts that enhance the rate and/or improve the thermodynamic efficiency of the reaction (i.e., lower the required overpotential). The four-electron, four-proton (4e−/ 4H+) ORR is a complex chemical transformation, and controlling the reaction selectivity presents a significant challenge. Generation of deleterious byproducts, such as superoxide or hydrogen peroxide, limits the energetic efficiency of the reaction and creates other complications, such as degradation of the carbon electrode in fuel cells. Formation of water is thermodynamically more favorable than formation of the partially reduced species. For example, the 4e−/4H+ reduction of O2 to H2O has a standard potential of 1.23 V, while the 2e−/2H+ reduction of O2 to H2O2 has a standard © 2019 American Chemical Society

O2 + 4H+ + 4e− → 2H 2O

E° = 1.23 V vs NHE

(1)

O2 + 2H+ + 2e− → H 2O2

E° = 0.68 V vs NHE

(2)

Molecular catalysts have played an important role in ORR studies because they are amenable to thorough mechanistic investigation, including systematic structure−activity analyses and characterization of reaction intermediates.6−10 Macrocyclic cobalt complexes are among the most widely studied classes of molecular ORR catalysts,11−15 with the first examples reported more than 50 years ago.16 Efforts to control the reaction selectivity have been an enduring focus of these studies. Mononuclear Co-macrocycles often favor the formation of Received: February 27, 2019 Published: June 7, 2019 1024

DOI: 10.1021/acscentsci.9b00194 ACS Cent. Sci. 2019, 5, 1024−1034

Research Article

ACS Central Science

Figure 1. Summary of representative macrocyclic cobalt ORR catalysts reported previously11−15 and their selectivities for O2 reduction.

H2O2 from O2 (eq 2).11−15 Although this reaction has been the focus of renewed attention due to growing interest in electrochemical hydrogen peroxide synthesis,17−19 historical efforts have been primarily focused on selective 4e−/4H+ reduction of O2 to H2O via strategic catalyst design (Figure 1). Prominent examples included cofacial porphyrins (a, b) with two redox-active Co centers available to deliver the four electrons needed for O2 reduction to water;20−22 4-pyridylsubstituted porphyrins with appended redox-active units, such as [Ru(NH3)5]2+ (c),23−25 that could provide additional electrons needed for O2 reduction at a single Co center; and “hangman” porphyrins (d) and related complexes capable of guiding proton delivery to the distal oxygen of a CoIII(OOH) intermediate to promote O−O bond cleavage.26,27 While some of the tailored catalysts in Figure 1a achieve high selectivity for the production of H2O (e.g., 99% with catalyst a),20 most generate mixtures of H2O and H2O2. Alongside these relatively sophisticated catalyst designs, certain mononuclear Co complexes have been identified that achieve high selectivity for H2O (up to 99% with the corrole catalyst e).28−31 Efforts to control ORR product selectivity with molecular catalysts include diverse catalyst classes, including Fe,32−34 Cu,35,36 and Mn37−39 complexes. New catalyst designs commonly feature multinuclear metal complexes11,13,20−22,35,36,40 or implement other approaches to control the relative rates of electron41−43 and/or proton transfer26,27,34,43,44 as a means to facilitate O−O cleavage and avoid H2O2 generation. Elucidation of the factors that contribute to product selectivity are complicated, however, by the manner in which these catalysts are studied. Molecular catalysts are often immobilized on an electrode prior to analysis, and the immobilization method (e.g., physisorption on edge-plane pyrolytic graphite or incorporation in a conductive “ink” containing Nafion polymer, etc.) can

influence the H2O/H2O2 product ratio, even with the same catalyst.45 These observations suggest that functional groups on the electrode surface or within the supporting matrix (e.g., sulfonic acid groups in Nafion), or intermolecular interactions between immobilized catalysts, influence the product selectivity and hinder interpretation of catalyst structure−selectivity data. Molecular ORR catalysts may be investigated under homogeneous conditions by using a chemical reductant. This approach facilitates mechanistic studies, and its use in the study of mononuclear macrocyclic Co complexes highlights the intrinsic preference for H2O2 formation with these catalysts (cf. Figure 1b).46−52 We recently studied a series of Co complexes bearing pseudo-macrocyclic N2O2-type ligands [e.g., bis(ketiminate), salen derivatives] that led to highly selective H2O2 production in methanol with decamethylferrocene (Fc*) as the reductant and acetic acid as the proton source.51,52 A parallel study showed that H2O is the sole product when phydroquinone was used as a combined source of electrons and protons.53,54 Mechanistic studies of the latter reaction showed that p-hydroquinone intercepts a CoII(HOOH) intermediate, facilitating O−O cleavage prior to release of H2O2. These results prompted us to consider whether similar H2O-selective O2 reduction could be achieved with independent sources of electrons and protons, similar to that required in the electrochemical ORR. Here, we show that the cobalt(II) porphyrin complex Co(porOMe) (1) (porOMe = meso tetra(4-methoxyphenyl)porphyrin)9,15 promotes highly selective formation of either H2O2 or H2O, depending on the identity of the Brønsted acid used in the reaction. The ORR thermodynamic potentials, EO2/H2O2 and EO2/H2O, exhibit Nernstian dependences on the acid pKa,55 while the CoIII/II potential of the catalyst is independent of the acid pKa. Selective formation of H2O is 1025

DOI: 10.1021/acscentsci.9b00194 ACS Cent. Sci. 2019, 5, 1024−1034

Research Article

ACS Central Science

similar to that described previously49,54,61−63 (see also Section VII in the Supporting Information). The product selectivity is high in all cases (≥93%, Table 1), but the product identity changes sharply, from H2O2 to H2O, between the pKa values of 2,6-(HO)2BA (3.6) and TFAH (4.9). Given the high product selectivity, the TOFs reported in Table 1 are calculated for the formation of a single product (H2O2 or H2O) and correspond to the initial rates of H2O2 or H2O formation (mM s−1) normalized to the catalyst concentration (mM). The results of these studies, summarized in Table 1, reveal a direct influence of the acid strength on the reaction rates and product selectivity. A plot of the logarithm of the TOF versus the pKa of the Brønsted acid reveals linear free energy correlations, with a deviation in slope corresponding to the change in product identity from H2O2 to H2O as the acid pKa increases (Figure 2). A similar correlation between acid pKa

observed when the CoIII/II potential is above the thermodynamic potential for O2 reduction to H2O2 (EO2/H2O2). Kinetic studies show that the logarithm of the catalytic turnover frequency [log(TOF)] scales linearly with the acid pKa but with different slopes for the formation of H2O and H2O2. These observations are complemented by kinetic, EPR spectroscopic, and voltammetric studies, as well as density functional theory (DFT) calculations, that illuminate the mechanistic basis for the formation of the different products. The ability to use scaling relationships to control product selectivity, as demonstrated herein, has broad implications for molecular electrocatalysis.



RESULTS AND DISCUSSION ORR Rates and H2O2/H2O Product Selectivity with Different Brønsted Acids. The reduction of O2 catalyzed by Co(porOMe) 1 was investigated in DMF with decamethylferrocene (Fc*) as the reductant and a series of different Brønsted acids, employing buffered conditions with an equimolar mixture of each acid and its conjugate base. Acid sources included the following (with their abbreviation and pKa in DMF in parentheses): [(DMF)H][ClO4] (DMF-H+, 1.6);56 ptoluenesulfonic acid (TsOH, 2.5);57 1-propanesulfonic acid (C 3 H 7 SO 3 H, 2.9); 56 2,6-dihydroxybenzoic acid (2,6(HO)2BA, 3.6);58,59 trifluoroacetic acid (TFAH, 4.9);59,60 oxalic acid, [(CO2H)2, 5.9];56 dichloroacetic acid (DCAH, 7.5);58,59 and maleic acid [CH2(CO2H)2, 7.9).58,59 The ORR rates were monitored by following the growth of the optical absorption band at 780 nm, corresponding to the conversion of Fc* to Fc*+ during the reaction (see Figure S5 in the Supporting Information).15,51,52 Upon completion of the reaction, the ORR selectivity was established by using iodometric titration and a TiIV(O)SO4 colorimetric assay to quantify the amount of H2O2 present in the reaction mixture,

Figure 2. Scaling relationship between log(TOF) and acid pKa (black; left axis), which exhibits a different correlation for the selective production of H2O2 or H2O (blue; right axis) catalyzed by 1 (see Table 1 for reaction conditions).

and ORR rates with a molecular catalyst was described recently by Pegis, Mayer, and co-workers,57 albeit without a change in product selectivity or slope of the correlation. Assessment of ORR Thermodynamics and Kinetics. The pKa of the Brønsted acid will influence the thermodynamic O2 reduction potential (i.e., the driving force for O2 reduction), and the thermodynamic potential for the ORR under nonaqueous conditions was determined as developed55 and implemented51,64,65 in several recent studies (see Section

Table 1. Selectivity of the ORR Catalyzed with 1 under Different Buffered Conditionsa entry

HA/A−

pKa (HA(DMF))

1 2 3 4 5 6 7 8

DMF-H+/DMF TsOH/NaOTsb C3H7SO3H/C3H7SO3Na 2,6-(HO)2BA/2,6-(NaO)2BAc TFAH/NaTFAd oxalic acid/lithium oxalate DCAH/NaDCAe maleic acid/sodium malonate

1.6 2.5 2.9 3.6 4.9 5.9 7.5 7.9

TOF (s−1) 7.7 5.5 2.6 2.0 9.0 7.1 5.3 4.5

× × × × × × × ×

10−2 10−2 10−2 10−2 10−3 10−3 10−3 10−3

selectivity H2O:H2O2 (%) 7:93 6:94 1:99 3:97 95:5 95:5 98:2 98:2

Reaction conditions: [acid] = [conjugate base] = 10 mM, [1] = 5 μM, [Fc*] = 1 mM, 1 atm O2, room temperature. bTsONa: sodium ptoluenesulfonate monohydrate. c2,6-(NaO)2BA: sodium 2,6-hydroxybenzoate. dNaTFA: sodium trifluoroacetate. eNaDCA: sodium dichloroacetate. a

1026

DOI: 10.1021/acscentsci.9b00194 ACS Cent. Sci. 2019, 5, 1024−1034

Research Article

ACS Central Science II of the Supporting Information for details). Briefly, the reference H+/H2 open-circuit potential (OCP) was determined under each of the buffered conditions (Figure S2),66 and the thermodynamic potentials for the O2/H2O2 and O2/ H2O redox couples in DMF were then estimated by adding 0.68 V (for O2/H2O2) and 1.23 V (for O2/H2O) to the H+/H2 OCP. Corrections associated with the free energy for transferring the H2O2 and H2O products of the ORR from water to DMF are negligible because the solvation free energy of H2O2 and H2O in DMF is similar to that in H2O.67,68 Application of this methodology provided the basis for a Pourbaix-like diagram correlating the thermodynamic potentials for H+/H2, O2/H2O2, and O2/H2O with the pKa of the Brønsted acids in organic media (filled circles, Figure 3). The

Figure 4. Linear free energy correlations between log(TOF) and the ηeff for the O2 reduction catalyzed by 1 under the conditions summarized in Table 1. A variation of this plot, with the ηeff values adjusted to account for the background concentration of H2O or H2O2, is provided in the Supporting Information (see Section VIII and Figure S11).

formed with weaker acids that contribute to lower ηeff. The inflection point occurs very close to the thermodynamic potential for the reduction of O2 to H2O2 (ηeff = 0.55 V). Thus, H2O2 is the observed product whenever the CoIII/II potential is below the O2/H2O2 potential. On the other hand, H2O is the observed product when the CoIII/II potential is above the O2/ H2O2 potential (i.e., where H2O2 production is thermodynamically unfavorable). Kinetic, EPR Spectroscopy, and Voltammetry Studies To Probe the ORR Mechanism under Strong and Weak Acid Conditions. The different slopes in the log(TOF) versus ηeff plot in Figure 4 implicate a change in the catalytic mechanism that depends on the strength of the acid present in the reaction mixture. In order to explore this hypothesis, mechanistic studies of the ORR catalyzed by 1 were performed with a representative strong and weak acid, [DMF-H][ClO4] (pKa = 1.6) and DCAH (pKa = 7.5). Initial rate data were collected to establish the catalytic rate law under both conditions (Figure 5, see also Section IX for full details). With [DMF-H][ClO4] as the acid, the catalytic rate exhibits a first-order dependence on [1], [H+], and [O2], but no dependence on [Fc*] (Figure 5a and eq 4). Analogous data with DCAH as the acid revealed that the catalytic rate exhibits a first-order dependence on [1], [H+], and [Fc*], but no dependence on [O2] (Figure 5b and eq 5).

Figure 3. Correlations of pKa(DMF) of acids with H+/H2 and O2/ H2O2 redox couples. The reduction potential of O2/H2O2 exhibits a Nernstian dependence on the pKa values of acids, while the E1/2(CoIII/II) of 1 is independent of the acid pKa. See Table 1 for the acid identities and their corresponding pKa values. The O2/H2O2 and O2/H2O redox potentials have been adjusted to account for the non-standard-state background concentrations of H2O and H2O2 used in this study. See Section IV of the Supporting Information for considerations.

potentials exhibit Nernstian trends, with slopes of 59 mV/pKa. Meanwhile, the CoIII/II redox potential of 1 (E1/2(CoIII/II)), determined under the same conditions, is unaffected by the identity of the Brønsted acid (open red circles, Figure 3). The latter observation suggests that the conjugate base of the acids does not coordinate to the Co center or otherwise alter the CoIII/II potential under the buffered conditions. For molecular electrocatalysts, an effective overpotential (ηeff) may be defined from the difference between the thermodynamic potential for O2 reduction under the reaction conditions (EORR) and the catalyst redox potential that initiates catalytic turnover, in this case E1/2(CoIII/II) (eq 3).69−71 The ηeff values, together with the catalytic ORR rates (i.e., TOFs) presented in Table 1, enable analysis of the linear free energy relationships (LFERs) between the log(TOF) and the ηeff values (Figure 4).72 ηeff = EORR − E1/2(Co

III/II

)

O2 /H 2O2 rate = −1/2 d[Fc*] /dt = k[Co][HA][O2 ]

(4)

O2 /H 2O rate = −1/4 d[Fc*] /dt = k[Co][HA][Fc*]

(5)

Electron paramagnetic resonance (EPR) spectroscopic analysis of 1 clearly showed the formation of a Co-superoxide [CoIII(O2•)] adduct under aerobic conditions (Figure 6a,b). The CoIII(O2•) species reacts rapidly upon the addition of a strong acid, DMF-H+, resulting in the disappearance of most of the EPR signal (Figure 6c). This behavior is consistent with an autoxidation pathway to generate CoIII species, as has been reported previously.73−76 In contrast, no reaction was observed

(3)

The LFER plot in Figure 4 shows that the ORR rate increases with increasing driving force (i.e., higher ηeff);33 however, the plot features an inflection point that correlates with the change in the product identity. H2O2 is formed with stronger acids that contribute to higher ηeff, while H2O is 1027

DOI: 10.1021/acscentsci.9b00194 ACS Cent. Sci. 2019, 5, 1024−1034

Research Article

ACS Central Science

Figure 5. Kinetic data for the reduction of O2 catalyzed by 1 in the presence of (a) DMF-H+ and (b) DCAH as the acid, obtained by monitoring the initial rates of Fc*+ formation. (a) For 2e−/2H+ O2 reduction to H2O2, the rate law exhibits first-order dependence on [Co], [HClO4], and [O2] but no dependence on [Fc*] (eq 4). (b) For 4e−/4H+ O2 reduction to H2O, the rate law exhibits first-order dependence on [Co], [DCAH], and [Fc*] but no dependence on [O2] (eq 5). See the Supporting Information for detailed experimental conditions.

Figure 6. X-band EPR spectra of CoII complex 1 in the absence and presence of O2 and, with the latter, in the presence of added acid. (a) EPR spectra of 1 (1 mM) in N2-saturated DMF. g values = 2.050, 2.015, 2.095; A = 15, 15, 20 G. (b) EPR spectra of 1 (1 mM) in O2-saturated DMF. g values = 2.050, 2.015, 2.095; A = 15, 15, 20 G. (c) EPR spectral evidence for reaction of CoIII(O2•) with HClO4. (d) EPR spectral evidence that CoIII(O2•) does not undergo protonation by DCAH. EPR parameters: microwave frequency = 9.46 GHz, microwave power = 10.4 mW, modulation frequency = 100 kHz, and modulation amplitude = 10 G. Temperature = 110 K. See the Supporting Information for additional details.

Figure 7. DPV of 1 under aerobic conditions (1 atm air) in the (a, b) absence and (c, d) presence of Fc* (the latter corresponding to catalytic conditions). Conditions: 25−100 μM 1, 5−10 mM acid, 4−6 mM Fc*, 0.1 M [NBu4][ClO4], 10 mL of DMF. See the Supporting Information for full experimental details.

between the CoIII(O2•) species and the weak acid, DCAH (Figure 6d). Differential pulse voltammetry (DPV) measurements provide complementary insights into the identity of the cobalt species present under the different conditions.77,78 Complex 1

exhibits an anodic peak at 0.28 V (vs Fc*+/0) in DMF under N2, and it is unaffected by the presence of DMF-H+ or DCAH (Figure S20). When a solution of 1 is exposed to 1 atm O2, a new cathodic peak is observed at lower potential (0.08 V) and is attributed to the CoIII(O2•) species detected by EPR 1028

DOI: 10.1021/acscentsci.9b00194 ACS Cent. Sci. 2019, 5, 1024−1034

Research Article

ACS Central Science

computed overpotentials show relatively good agreement with the experimental data (Table 3). On the other hand, the

spectroscopy (cf. Figure 6). Upon addition of [DMFH][ClO4], a single cathodic peak at substantially higher potential (0.37 V) is evident (Figure 7a), consistent with the conversion of CoIII(O2•) into a CoIII species derived from autoxidation of Co(II) in the presence of strong acid (postulated above, Figure 6c). When DCAH is used rather than [DMF-H][ClO4], the cathodic peak associated with the CoIII(O2•) species shifts to slightly higher potential (0.17 V, Figure 7b), an effect attributed to hydrogen bonding between the acid and the CoIII(O2•) species, similar to observations made previously with a N2O2-ligated CoIII(O2•) complex.52 DPV was also used to probe the nature of the resting-state Co species under catalytic conditions (i.e., 1 in the presence of O2, acid, and Fc*). An anodic peak was observed at 0.33 V with DMF-H+ as the acid source (Figure 7c). This result, together with the first-order dependence of the rate on [O2] and [H+] and zero-order dependence on [Fc*] (cf. Figure 5a), is rationalized by a CoII catalyst resting state that undergoes reaction with O2 and DMF-H+ in the turnover-limiting step(s) of the catalytic reaction. In contrast, a cathodic peak was detected at 0.36 V with DCAH as the acid source (Figure 7d), implicating a CoIII catalyst resting state under the weak acid reaction conditions. The zero-order dependence of the rate on [O2] (cf. Figure 5b) is consistent with a Co/O2 adduct as the resting state; however, the observed redox potential is considerably higher than that of the CoIII(O2•) species (cf. Figure 7b). Therefore, we tentatively assign the catalyst resting state to a CoIII(OOH) species, similar to that recently identified for the ORR catalyzed by a N2O2-ligated Co complex.52 The first-order dependence of the rate on [H+] and [Fc*] supports turnover-limiting electron−proton transfer (EPT) to this species, which could take place in a sequential or concerted process. Free Energy Profiles for the Co(porOMe)-Catalyzed ORR in the Presence of Weak and Strong Acid. The experimental data and mechanistic proposals presented above and summarized in Table 2 provided the basis for DFT

Table 3. Experimental and Computed Effective Overpotentials for the ORR Catalyzed by 1 To Produce H2O2 (cf. eq 3)a

strong acid (DMF-H+) weak acid (DCAH)

strong acid conditions

weak acid conditions

H2O2 rate ∝ [Co]1[HA]1[O2]1 CoII(porOMe)

H2O rate ∝ [Co]1[HA]1[Fc*]1 CoIII(OOH) species

ηeff(O2/H2O2) calculated (V)

+0.24 −0.16

+0.12b −0.28c

a

See Section XII in the Supporting Information for computational details. bO2 + 2H(DMF)2+ + 2[CoII]DMF → H2O2 + 2DMF + 2[CoIII]DMF2. cO2 + 4DCAH−DMF + 4[CoII]DMF → 2H2O + 4DCA− + 4[CoIII]DMF2.

relative free energies of intermediates along the reaction pathway (i.e., those observed experimentally and analyzed by DFT methods) will be influenced by the redox potential of the stoichiometric reductant (in this case, Fc*) that delivers electrons to the Co-based intermediates. The simplified free energy diagrams presented in Figure 8 seek to superimpose both of these considerations: (1) the overall reaction thermodynamics derived from the Co(porOMe) catalyst and DMF-H+ and DCAH as representative strong and weak acids and (2) the relative free energies of relevant reaction intermediates, insights into which were gained from DFT calculations (see Figure S25 in Section XII of the Supporting Information). The free energy diagrams in Figure 8 and the Pourbaix diagram in Figure 3 show how the formation of H2O2 becomes thermodynamically unfavorable when changing from a strong to a weak acid. The acid strength also influences the free energies of the reaction intermediates by changing the driving force for steps involving proton transfer (PT) and electron− proton transfer (EPT). DFT calculations indicate that the free energies associated with the progression from 1 → 1a → 1b → 1, leading to the formation of H2O2, are all downhill with a strong acid, and the corresponding steps leading to formation of H2O as the final product are even more thermodynamically favorable. The selective formation of H2 O 2 observed experimentally under strong acid conditions indicates that the reaction product is kinetically controlled and may be rationalized by a high kinetic barrier for O−O cleavage in the formation of water. When using a weak acid, formation of H2O 2 is thermodynamically unfavorable, and selective formation of H2O is observed. The required O−O cleavage step could proceed via intermediate 1b or 1c. The former pathway involves a relatively high-energy intermediate (1b), but it would allow the kinetically challenging O−O bond cleavage to occur as late as possible in the mechanism where it would also benefit from the increased driving force relative to EPT-induced O−O cleavage from the CoIII(OOH) species (1c). On the other hand, the latter pathway resembles wellestablished reactivity involving Fe III −OOH intermediates.41,84,85 Fundamental studies of other Co complexes86−90 have characterized relevant intermediates and provide a potential starting point for future studies to explore mechanistic questions concerning the O−O cleavage pathway. Implications for the ORR with Molecular Catalysts. The correlation between reaction rates and driving force in this system is best understood from the perspective of proton transfer, rather than electron transfer steps, as conveyed by the

Table 2. Summary of Experimental Studies of the ORR Mechanism product catalytic rate law proposed resting state

ηeff(O2/H2O2) experiment (V)

calculations to probe the relative free energies of intermediates under the strong and weak acid conditions. Free energy profiles were computed for the reactions with DMF-H+ and DCAH as the acid with Gaussian 09,79 using the BP86 exchange-correlation functional80,81 and the 6-31G**82 electronic basis set with additional diffuse basis functions on select oxygen atoms. The structures were optimized in the gas phase, followed by the calculation of solvation free energies in DMF using the SMD implicit solvation model.83 A complication associated with the use of DFT-based results to create free energy diagrams to rationalize all of the experimental observations presented herein is that CoII is the source of electrons used for O2 reduction in the catalytic cycle and provides the basis for the effective overpotential of the reaction (cf. eq 3). Thus, the overall thermodynamics of O2 reduction are defined with respect to the CoIII/II potential, and 1029

DOI: 10.1021/acscentsci.9b00194 ACS Cent. Sci. 2019, 5, 1024−1034

Research Article

ACS Central Science

Figure 8. Simplified free energy profiles for Co(porOMe)-catalyzed O2 reduction with a (A) strong and (B) weak acid, leading to the formation of H2O2 and H2O, respectively. These qualitative free energy profiles incorporate insights from both experimental and computational data, and they are not intended to convey precise quantitative information for reasons discussed in the text.

correlations in Figures 2−4. Specifically, the change in driving force (i.e., the effective overpotential) arises from the change in pKa of the Brønsted acid, which influences the thermodynamic potential for O2 reduction while not affecting the CoIII/II redox potential (cf. Figure 3). Thus, stronger Brønsted acids increase the difference between the ORR thermodynamic potential and the CoIII/II redox potential (i.e., the ηeff; see eq 3). The influence of the driving force on the reaction rates is evident from the rate laws for the formation of H2O2 and H2O in eqs 4 and 5, both of which feature a first-order dependence on acid concentration. Stronger acids will have a higher concentration of free H+, as defined by the Ka, and thereby lead to faster rates at higher driving force. The results presented herein demonstrate the ability to use scaling relationships to predict and control the product selectivity for H2O2 and H2O during the catalytic ORR. The thermodynamics of O2/H2O2 and O2/H2O exhibit a Nernstian dependence on the pKa of the Brønsted acid used in the reaction, while the catalyst CoIII/II potential is unaffected by the acid pKa (cf. Figure 3). These different correlations make it possible to access conditions that favor either H2O2 or H2O. With weak acids (high pKa), the CoIII/II potential is higher than the thermodynamic potential for O2 reduction to H2O2 (EO2/H2O2), and the ORR selectively generates H2O. With strong acids (low pKa), the CoIII/II potential is lower than the EO2/H2O2 potential, and the reaction selectively generates H2O2. This abrupt switch in product selectivity at the thermodynamic

EO2/H2O2 potential empirically aligns with the effective overpotential (ηeff) for the ORR, which is defined with respect to the E1/2(CoIII/II) (cf. eq 3 and Figure 4). It is perhaps surprising that this binary dependence of product selectivity on the relationship between the ηeff and EO2/H2O2 potential has not been demonstrated previously, but multiple considerations rationalize the lack of precedent. First, many ORR studies with macrocyclic Co catalysts have been conducted under strongly acidic, aqueous conditions with immobilized catalysts (cf. Figure 1a). The catalytic E1/2 values lie below the EO2/H2O2 potential in virtually all of these cases. The cofacial Co-porphyrin complex a in Figure 1a is a rare exception;91 however, the highly selective formation of H2O in this case has been attributed to the cofacial bimetallic structure, which is thought to promote O−O cleavage.11,22 The results presented herein, however, reveal that a binuclear catalyst structure is not required to achieve highly selective formation of H2O. Instead, the reaction simply needs to be conducted under conditions that are thermodynamically unfavorable with respect to the formation of H2O2. Otherwise, the H2O/H2O2 selectively will be dictated by the relative kinetic barriers leading to the two products; however, selective formation of H2O will necessarily feature a large overpotential in such cases. Another consideration is that many homogeneous ORR studies, including those with complexes h−l in Figure 1b, have been performed in organic solvents. A reliable methodology for defining the thermodynamic ORR potentials under these 1030

DOI: 10.1021/acscentsci.9b00194 ACS Cent. Sci. 2019, 5, 1024−1034

ACS Central Science



conditions was established only recently, however, and the ηeff for most ORR precedents in organic solvent was not determined.45 On the basis of the results described herein, it is likely that the highly selective formation of H2O2 with catalysts h−l in Figure 1b arises from catalyst CoIII/II potentials that fall below the EO2/H2O2 values. The results of this study also may be compared to a recent study by Nocera and co-workers.40 The authors analyzed ORR results with a series of Co- and Fe-based molecular catalysts and concluded that “high ORR selectivities for H2O is a result of large effective overpotentials for the reaction, achieved by the use of strong acids.” This conclusion contrasts the results presented herein, which achieve high selectivity for H2O by using a weak acid, resulting in a sufficiently low effective overpotential that formation of H2O2 is thermodynamically unfavorable. The majority of catalyst systems that exhibit high selectivity for H2O in the study by Nocera and co-workers (and the complementary studies by Mayer and co-workers33) operate with very high effective overpotentials (>1.2 V). It is possible that reactions with such a large driving force can access kinetically facile O−O cleavage pathways, either via direct O−O cleavage from M−OOH intermediates, similar to 1a, or via two-electron reduction of M−(HOOH) intermediates, analogous to 1b (cf. Figure 8). Safety Statement. No unexpected or unusually dangerous safety hazards were encountered.

Research Article

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acscentsci.9b00194.



Experimental details for data acquisition and additional kinetic data (PDF)

AUTHOR INFORMATION

Corresponding Authors

*E-mail: sharon.hammes-schiff[email protected]. *E-mail: [email protected]. ORCID

Zachary K. Goldsmith: 0000-0002-5556-4079 Biswajit Mondal: 0000-0002-2339-3076 Sharon Hammes-Schiffer: 0000-0002-3782-6995 Shannon S. Stahl: 0000-0002-9000-7665 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This research was supported as part of the Center for Molecular Electrocatalysis, an Energy Frontier Research Center, funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences. PES was supported by an NSF GRF (DGE-1752134). We acknowledge the NSF (CHE-0741901) and the Paul and Margaret Bender Fund for support of the EPR facility at UW-Madison.



CONCLUSIONS This study demonstrates that selective 4e−/4H+ or 2e−/2H+ reduction of O2 may be accomplished with the same monomeric cobalt porphyrin ORR catalyst, Co(porOMe) 1. The switch in selectivity is achieved by varying the pKa of the Brønsted acid used as the source of protons for the reaction. The acid pKa systematically modulates the thermodynamic ORR potentials for O2 → H2O2 and O2 → H2O in a Nernstian manner, but it does not influence the E1/2(CoIII/II) of 1. These differences in scaling relationships provide the basis for precise control over the product selectivity. The thermodynamically controlled formation of water arises when the CoIII/II potential is above the thermodynamic potential for the production of H2O2. In contrast, selective, kinetically favored formation of H2O2 is observed when the CoIII/II potential lies below the H2O2 thermodynamic potential. It will be important to extend the results of the present study to other catalyst systems, as the specific Brønsted acid scaling relationships identified here are not expected to apply universally. Changes in the identity and/or number of metal ions in the catalyst, the structure and electronic properties of the ancillary ligand, and specific components of the reaction system (solvent, Brønsted acid, conjugate base, etc.) could lead to changes in the reaction mechanism that will, in turn, influence the magnitude or nature of the scaling relationships. Nonetheless, nearly all proton-coupled redox process should be subject to analogous scaling relationships that could provide opportunities to manipulate product selectivity. The present study highlights the value of understanding the thermodynamic properties of such catalytic reactions as the basis for controlling product selectivity. This principle has important implications beyond O2 reduction, extending to CO2 and N2 reduction and alcohol oxidation in addition to a multitude of other protoncoupled redox processes.



REFERENCES

(1) Solomon, E. I.; Stahl, S. S. Introduction: Oxygen Reduction and Activation in Catalysis. Chem. Rev. 2018, 118, 2299−2301. (2) Babcock, G. T. How Oxygen Is Activated and Reduced in Respiration. Proc. Natl. Acad. Sci. U. S. A. 1999, 96, 12971−12973. (3) Katsounaros, I.; Cherevko, S.; Zeradjanin, A. R.; Mayrhofer, K. J. J. Oxygen Electrochemistry as a Cornerstone for Sustainable Energy Conversion. Angew. Chem., Int. Ed. 2014, 53, 102−121. (4) Shao, M.; Chang, Q.; Dodelet, J.-P.; Chenitz, R. Recent Advances in Electrocatalysts for Oxygen Reduction Reaction. Chem. Rev. 2016, 116, 3594−3657. (5) Bard, A. J.; Faulkner, L. R. Electrochemical Methods: Fundamentals and Applications; Wiley: Hoboken, NJ, 2001. (6) Dey, S.; Mondal, B.; Chatterjee, S.; Rana, A.; Amanullah, S.; Dey, A. Molecular Electrocatalysts for the Oxygen Reduction Reaction. Nat. Rev. Chem. 2017, 1, 0098. (7) Zhang, W.; Lai, W.; Cao, R. Energy-Related Small Molecule Activation Reactions: Oxygen Reduction and Hydrogen and Oxygen Evolution Reactions Catalyzed by Porphyrin- and Corrole-Based Systems. Chem. Rev. 2017, 117, 3717−3797. (8) Gewirth, A. A.; Varnell, J. A.; DiAscro, A. M. Nonprecious Metal Catalysts for Oxygen Reduction in Heterogeneous Aqueous Systems. Chem. Rev. 2018, 118, 2313−2339. (9) Pegis, M. L.; Wise, C. F.; Martin, D. J.; Mayer, J. M. Oxygen Reduction by Homogeneous Molecular Catalysts and Electrocatalysts. Chem. Rev. 2018, 118, 2340−2391. (10) Adam, S. M.; Wijeratne, G. B.; Rogler, P. J.; Diaz, D. E.; Quist, D. A.; Liu, J. J.; Karlin, K. D. Synthetic Fe/Cu Complexes: Toward Understanding Heme-Copper Oxidase Structure and Function. Chem. Rev. 2018, 118, 10840−11022. (11) Collman, J. P.; Wagenknecht, P. S.; Hutchison, J. E. Molecular Catalysts for Multielectron Redox Reactions of Small Molecules: The “Cofacial Metallodiporphyrin” Approach. Angew. Chem., Int. Ed. Engl. 1994, 33, 1537−1554.

1031

DOI: 10.1021/acscentsci.9b00194 ACS Cent. Sci. 2019, 5, 1024−1034

Research Article

ACS Central Science

(31) Shi, C.; Anson, F. C. (5,10,15,20-Tetramethylporphyrinato)cobalt(II): A Remarkably Active Catalyst for the Electroreduction of O2 to H2O. Inorg. Chem. 1998, 37, 1037−1043. (32) Costentin, C.; Dridi, H.; Savéant, J.-M. Molecular Catalysis of O2 Reduction by Iron Porphyrins in Water: Heterogeneous versus Homogeneous Pathways. J. Am. Chem. Soc. 2015, 137, 13535−13544. (33) Pegis, M. L.; McKeown, B. A.; Kumar, N.; Lang, K.; Wasylenko, D. J.; Zhang, X. P.; Raugei, S.; Mayer, J. M. Homogenous Electrocatalytic Oxygen Reduction Rates Correlate with Reaction Overpotential in Acidic Organic Solutions. ACS Cent. Sci. 2016, 2, 850−856. (34) Bhunia, S.; Rana, A.; Roy, P.; Martin, D. J.; Pegis, M. L.; Roy, B.; Dey, A. Rational Design of Mononuclear Iron Porphyrins for Facile and Selective 4e−/4H+ O2 Reduction: Activation of O−O Bond by 2nd Sphere Hydrogen Bonding. J. Am. Chem. Soc. 2018, 140, 9444−9457. (35) Thorseth, M. A.; Tornow, C. E.; Tse, E. C. M.; Gewirth, A. A. Cu Complexes that Catalyze the Oxygen Reduction Reaction. Coord. Chem. Rev. 2013, 257, 130−139. (36) Thiyagarajan, N.; Janmanchi, D.; Tsai, Y.-F.; Wanna, W. H.; Ramu, R.; Chan, S. I.; Zen, J.-M.; Yu, S. S.-F. A Carbon Electrode Functionalized by a Tricopper Cluster Complex: Overcoming Overpotential and Production of Hydrogen Peroxide in the Oxygen Reduction Reaction. Angew. Chem., Int. Ed. 2018, 57, 3612−3616. (37) Jung, J.; Liu, S.; Ohkubo, K.; Abu-Omar, M. M.; Fukuzumi, S. Catalytic Two-Electron Reduction of Dioxygen by Ferrocene Derivatives with Manganese(V) Corroles. Inorg. Chem. 2015, 54, 4285−4291. (38) Hooe, S. L.; Rheingold, A. L.; Machan, C. W. Electrocatalytic Reduction of Dioxygen to Hydrogen Peroxide by a Molecular Manganese Complex with a Bipyridine-Containing Schiff Base Ligand. J. Am. Chem. Soc. 2018, 140, 3232−3241. (39) Passard, G.; Dogutan, D. K.; Qiu, M.; Costentin, C.; Nocera, D. G. Oxygen Reduction Reaction Promoted by Manganese Porphyrins. ACS Catal. 2018, 8, 8671−8679. (40) Passard, G.; Ullman, A. M.; Brodsky, C. N.; Nocera, D. G. Oxygen Reduction Catalysis at a Dicobalt Center: The Relationship of Faradaic Efficiency to Overpotential. J. Am. Chem. Soc. 2016, 138, 2925−2928. (41) Collman, J. P.; Devaraj, N. K.; Decréau, R. A.; Yang, Y.; Yan, Y.L.; Ebina, W.; Eberspacher, T. A.; Chidsey, C. E. D. A Cytochrome c Oxidase Model Catalyzes Oxygen to Water Reduction Under RateLimiting Electron Flux. Science 2007, 315, 1565−1568. (42) Tse, E. C. M.; Barile, C. J.; Kirchschlager, N. A.; Li, Y.; Gewargis, J. P.; Zimmerman, S. C.; Hosseini, A.; Gewirth, A. A. Proton Transfer Dynamics Control the Mechanism of O2 Reduction by a Non-precious Metal Electrocatalyst. Nat. Mater. 2016, 15, 754− 760. (43) Gautam, R. P.; Lee, Y. T.; Herman, G. L.; Moreno, C. M.; Tse, E. C. M.; Barile, C. J. Controlling Proton and Electron Transfer Rates to Enhance the Activity of an Oxygen Reduction Electrocatalyst. Angew. Chem., Int. Ed. 2018, 57, 13480−13483. (44) Carver, C. T.; Matson, B. D.; Mayer, J. M. Electrocatalytic Oxygen Reduction by Iron Tetra-arylporphyrins Bearing Pendant Proton Relays. J. Am. Chem. Soc. 2012, 134, 5444−5447. (45) Rigsby, M. L.; Wasylenko, D. J.; Pegis, M. L.; Mayer, J. M. Medium Effects Are as Important as Catalyst Design for Selectivity in Electrocatalytic Oxygen Reduction by Iron−Porphyrin Complexes. J. Am. Chem. Soc. 2015, 137, 4296−4299. (46) Fukuzumi, S.; Mochizuki, S.; Tanaka, T. Efficient Reduction of Dioxygen with Ferrocene Derivatives, Catalyzed by Metalloporphyrins in the Presence of Perchloric acid. Inorg. Chem. 1989, 28, 2459− 2465. (47) Fukuzumi, S.; Okamoto, K.; Gros, C. P.; Guilard, R. Mechanism of Four-Electron Reduction of Dioxygen to Water by Ferrocene Derivatives in the Presence of Perchloric Acid in Benzonitrile, Catalyzed by Cofacial Dicobalt Porphyrins. J. Am. Chem. Soc. 2004, 126, 10441−10449.

(12) Anson, F. C.; Shi, C.; Steiger, B. Novel Multinuclear Catalysts for the Electroreduction of Dioxygen Directly to Water. Acc. Chem. Res. 1997, 30, 437−444. (13) Rosenthal, J.; Nocera, D. G. Role of Proton-Coupled Electron Transfer in O−O Bond Activation. Acc. Chem. Res. 2007, 40, 543− 553. (14) Mahammed, A.; Gross, Z. Metallocorroles as Electrocatalysts for the Oxygen Reduction Reaction (ORR). Isr. J. Chem. 2016, 56, 756−762. (15) Fukuzumi, S.; Lee, Y.-M.; Nam, W. Mechanisms of TwoElectron versus Four-Electron Reduction of Dioxygen Catalyzed by Earth-Abundant Metal Complexes. ChemCatChem 2018, 10, 9−28. (16) Jasinski, R. A New Fuel Cell Cathode Catalyst. Nature 1964, 201, 1212. (17) Fukuzumi, S.; Yamada, Y.; Karlin, K. D. Hydrogen Peroxide as a Sustainable Energy Carrier: Electrocatalytic Production of Hydrogen Peroxide and the Fuel Cell. Electrochim. Acta 2012, 82, 493−511. (18) Lu, Z.; Chen, G.; Siahrostami, S.; Chen, Z.; Liu, K.; Xie, J.; Liao, L.; Wu, T.; Lin, D.; Liu, Y.; Jaramillo, T. F.; Nørskov, J. K.; Cui, Y. High-efficiency Oxygen Reduction to Hydrogen Peroxide Catalyzed by Oxidized Carbon Materials. Nat. Catal. 2018, 1, 156− 162. (19) Kim, H. W.; Ross, M. B.; Kornienko, N.; Zhang, L.; Guo, J.; Yang, P.; McCloskey, B. D. Efficient Hydrogen Peroxide Generation Using Reduced Graphene Oxide-based Oxygen Reduction Electrocatalysts. Nat. Catal. 2018, 1, 282−290. (20) Collman, J. P.; Marrocco, M.; Denisevich, P.; Koval, C.; Anson, F. C. Potent Catalysis of the Electroreduction of Oxygen to Water by Dicobalt Porphyrin Dimers Adsorbed on Graphite Electrodes. J. Electroanal. Chem. Interfacial Electrochem. 1979, 101, 117−122. (21) Collman, J. P.; Denisevich, P.; Konai, Y.; Marrocco, M.; Koval, C.; Anson, F. C. Electrode Catalysis of the Four-electron Reduction of Oxygen to Water by Dicobalt Face-to-face Porphyrins. J. Am. Chem. Soc. 1980, 102, 6027−6036. (22) Durand, R. R.; Bencosme, C. S.; Collman, J. P.; Anson, F. C. Mechanistic Aspects of the Catalytic Reduction of Dioxygen by Cofacial Metalloporphyrins. J. Am. Chem. Soc. 1983, 105, 2710−2718. (23) Shi, C.; Anson, F. C. Electrocatalysis of the Reduction of Molecular Oxygen to Water by Tetraruthenated Cobalt Mesotetrakis(4-pyridyl)porphyrin Adsorbed on Graphite Electrodes. Inorg. Chem. 1992, 31, 5078−5083. (24) Shi, C.; Anson, F. C. Cobalt meso-Tetrakis(N-methyl-4pyridiniumyl)porphyrin Becomes a Catalyst for the Electroreduction of O2 by Four Electrons When [(NH3)5Os]n+ (n = 2, 3) Groups Are Coordinated to the Porphyrin Ring. Inorg. Chem. 1996, 35, 7928− 7931. (25) Steiger, B.; Anson, F. C. [5,10,15,20-Tetrakis(4-((pentaammineruthenio)-cyano)phenyl)porphyrinato]cobalt(II) Immobilized on Graphite Electrodes Catalyzes the Electroreduction of O2 to H2O, but the Corresponding 4-Cyano-2,6-dimethylphenyl Derivative Catalyzes the Reduction Only to H2O2. Inorg. Chem. 1997, 36, 4138−4140. (26) McGuire, R., Jr; Dogutan, D. K.; Teets, T. S.; Suntivich, J.; Shao-Horn, Y.; Nocera, D. G. Oxygen Reduction Reactivity of Cobalt(II) Hangman Porphyrins. Chem. Sci. 2010, 1, 411−414. (27) Dogutan, D. K.; Stoian, S. A.; McGuire, R.; Schwalbe, M.; Teets, T. S.; Nocera, D. G. Hangman Corroles: Efficient Synthesis and Oxygen Reaction Chemistry. J. Am. Chem. Soc. 2011, 133, 131− 140. (28) Schechter, A.; Stanevsky, M.; Mahammed, A.; Gross, Z. FourElectron Oxygen Reduction by Brominated Cobalt Corrole. Inorg. Chem. 2012, 51, 22−24. (29) Shi, C.; Steiger, B.; Yuasa, M.; Anson, F. C. Electroreduction of O2 to H2O at Unusually Positive Potentials Catalyzed by the Simplest of the Cobalt Porphyrins. Inorg. Chem. 1997, 36, 4294−4295. (30) Song, E.; Shi, C.; Anson, F. C. Comparison of the Behavior of Several Cobalt Porphyrins as Electrocatalysts for the Reduction of O2 at Graphite Electrodes. Langmuir 1998, 14, 4315−4321. 1032

DOI: 10.1021/acscentsci.9b00194 ACS Cent. Sci. 2019, 5, 1024−1034

Research Article

ACS Central Science (48) Honda, T.; Kojima, T.; Fukuzumi, S. Proton-Coupled ElectronTransfer Reduction of Dioxygen Catalyzed by a Saddle-Distorted Cobalt Phthalocyanine. J. Am. Chem. Soc. 2012, 134, 4196−4206. (49) Mase, K.; Ohkubo, K.; Fukuzumi, S. Efficient Two-Electron Reduction of Dioxygen to Hydrogen Peroxide with One-Electron Reductants with a Small Overpotential Catalyzed by a Cobalt Chlorin Complex. J. Am. Chem. Soc. 2013, 135, 2800−2808. (50) Mase, K.; Ohkubo, K.; Fukuzumi, S. Much Enhanced Catalytic Reactivity of Cobalt Chlorin Derivatives on Two-Electron Reduction of Dioxygen to Produce Hydrogen Peroxide. Inorg. Chem. 2015, 54, 1808−1815. (51) Wang, Y.-H.; Pegis, M. L.; Mayer, J. M.; Stahl, S. S. Molecular Cobalt Catalysts for O2 Reduction: Low-Overpotential Production of H2O2 and Comparison with Iron-Based Catalysts. J. Am. Chem. Soc. 2017, 139, 16458−16461. (52) Wang, Y.-H.; Goldsmith, Z. K.; Schneider, P. E.; Anson, C. W.; Gerken, J. B.; Ghosh, S.; Hammes-Schiffer, S.; Stahl, S. S. Kinetic and Mechanistic Characterization of Low-Overpotential, H2O2-Selective Reduction of O2 Catalyzed by N2O2-Ligated Cobalt Complexes. J. Am. Chem. Soc. 2018, 140, 10890−10899. (53) Anson, C. W.; Ghosh, S.; Hammes-Schiffer, S.; Stahl, S. S. Co(salophen)-Catalyzed Aerobic Oxidation of p-Hydroquinone: Mechanism and Implications for Aerobic Oxidation Catalysis. J. Am. Chem. Soc. 2016, 138, 4186−4193. (54) Anson, C. W.; Stahl, S. S. Cooperative Electrocatalytic O2 Reduction Involving Co(salophen) with p-Hydroquinone as an Electron−Proton Transfer Mediator. J. Am. Chem. Soc. 2017, 139, 18472−18475. (55) Pegis, M. L.; Roberts, J. A. S.; Wasylenko, D. J.; Mader, E. A.; Appel, A. M.; Mayer, J. M. Standard Reduction Potentials for Oxygen and Carbon Dioxide Couples in Acetonitrile and N,N-Dimethylformamide. Inorg. Chem. 2015, 54, 11883−11888. (56) The pKa values of [(DMF)H][ClO4] (1.6), C3H7SO3H (2.9), and oxalic acid (5.9) in DMF are derived experimentally from opencircuit potential measurements of EH+/H2 in DMF. See Section II of the Supporting Information for details. (57) Pegis, M. L.; Wise, C. F.; Koronkiewicz, B.; Mayer, J. M. Identifying and Breaking Scaling Relations in Molecular Catalysis of Electrochemical Reactions. J. Am. Chem. Soc. 2017, 139, 11000− 11003. (58) Kolthoff, I. M.; Chantooni, M. K.; Smagowski, H. Acid-base Strength in N,N-dimethylformamide. Anal. Chem. 1970, 42, 1622− 1628. (59) Maran, F.; Celadon, D.; Severin, M. G.; Vianello, E. Electrochemical Determination of the pKa of Weak Acids in N,Ndimethylformamide. J. Am. Chem. Soc. 1991, 113, 9320−9329. (60) Haindl, M. H.; Schmid, M. B.; Zeitler, K.; Gschwind, R. M. What Is Your Actual Catalyst? TMS Cleavage Rates of Diarylprolinol Silyl Ethers Studied by in situ NMR. RSC Adv. 2012, 2, 5941−5943. (61) Mair, R. D.; Graupner, A. J. Determination of Organic Peroxides by Iodine Liberation Procedures. Anal. Chem. 1964, 36, 194−204. (62) Fukuzumi, S.; Kuroda, S.; Tanaka, T. Flavin Analog-metal Ion Complexes Acting as Efficient Photocatalysts in the Oxidation of pMethylbenzyl Alcohol by Oxygen under Irradiation with Visible Light. J. Am. Chem. Soc. 1985, 107, 3020−3027. (63) Lee, Y.; Park, G. Y.; Lucas, H. R.; Vajda, P. L.; Kamaraj, K.; Vance, M. A.; Milligan, A. E.; Woertink, J. S.; Siegler, M. A.; Narducci Sarjeant, A. A.; Zakharov, L. N.; Rheingold, A. L.; Solomon, E. I.; Karlin, K. D. Copper(I)/O2 Chemistry with Imidazole Containing Tripodal Tetradentate Ligands Leading to μ-1,2-Peroxo−Dicopper(II) Species. Inorg. Chem. 2009, 48, 11297−11309. (64) Pegis, M. L.; McKeown, B. A.; Kumar, N.; Lang, K.; Wasylenko, D. J.; Zhang, X. P.; Raugei, S.; Mayer, J. M. Homogenous Electrocatalytic Oxygen Reduction Rates Correlate with Reaction Overpotential in Acidic Organic Solutions. ACS Cent. Sci. 2016, 2, 850−856. (65) A modified protocol was described and used in the following study: Passard, G.; Ullman, A. M.; Brodsky, C. N.; Nocera, D. G.

Oxygen Reduction Catalysis at a Dicobalt Center: The Relationship of Faradaic Efficiency to Overpotential. J. Am. Chem. Soc. 2016, 138, 2925−2928. (66) Roberts, J. A. S.; Bullock, R. M. Direct Determination of Equilibrium Potentials for Hydrogen Oxidation/Production by Open Circuit Potential Measurements in Acetonitrile. Inorg. Chem. 2013, 52, 3823−3835. (67) Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. Universal Solvation Model Based on Solute Electron Density and on a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions. J. Phys. Chem. B 2009, 113, 6378−6396. (68) Marenich, A. V.; Kelly, C. P.; Thompson, J. D.; Hawkins, G. D.; Chambers, C. C.; Giesen, D. J.; Winget, P.; Cramer, C. J.; Truhlar, D. G. Minnesota Solvation Database, version 2012; University of Minnesota: Minneapolis, 2012. (69) 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. (70) Appel, A. M.; Helm, M. L. Determining the Overpotential for a Molecular Electrocatalyst. ACS Catal. 2014, 4, 630−633. (71) Rountree, E. S.; McCarthy, B. D.; Eisenhart, T. T.; Dempsey, J. L. Evaluation of Homogeneous Electrocatalysts by Cyclic Voltammetry. Inorg. Chem. 2014, 53, 9983−10002. (72) The ORR studies were conducted with a background concentration of H2O2 or H2O to ensure a stable thermodynamic reference state (100 μM H2O2 and 100 mM H2O, respectively, under the conditions that led to selective formation of these products). For an LFER plot similar to Figure 3 that incorporates corrections to the ηeff values to account for these non-standard-state background concentrations of products, see Figure S11 in the Supporting Information. Additional details are provided in Section VIII of the Supporting Information. (73) Pignatello, J. J.; Jensen, F. R. Autoxidation of Transition-metal Complexes. Reaction of a 1:1 Cobalt-molecular Oxygen Complex with Acids to Yield Hydrogen Peroxide. Kinetics and Mechanism. J. Am. Chem. Soc. 1979, 101, 5929−5939. (74) Simandi, L. I.; Savage, C. R.; Schelly, Z. A.; Nemeth, S. Kinetics and Mechanism of the Oxygenation of Bis(dimethylglyoximato)cobalt(II) in Methanol. Inorg. Chem. 1982, 21, 2765−2769. (75) Nishinaga, A.; Kondo, T.; Matsuura, T. Oxygenation of Cobalt(II) Schiff base Complexes in Alcohols. Chem. Lett. 1985, 14, 905−908. (76) Pavlovic, D.; Asperger, S.; Dokuzovic, Z.; Jurisic, B.; Ahmeti, X.; Sertic, M.; Murati, I. Reactions of Cobalt(II) Protoporphyrin IX Dimethyl Ester, [COIIP], and [CoIIIP(Cl)] in Co-ordinating Aliphatic Alcohols. J. Chem. Soc., Dalton Trans. 1985, 1095−1101. (77) Brown, A. P.; Anson, F. C. Cyclic and Differential Pulse Voltammetric Behavior of Reactants Confined to the Electrode Surface. Anal. Chem. 1977, 49, 1589−1595. (78) Bard, A. J.; Faulkner, L. R. Electrochemical Methods: Fundamentals and Applications; Wiley: Hoboken, NJ, 2001. (79) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. 1033

DOI: 10.1021/acscentsci.9b00194 ACS Cent. Sci. 2019, 5, 1024−1034

Research Article

ACS Central Science D.; Farkas, O.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09; Gaussian, Inc.: Wallingford, CT, 2009. (80) Becke, A. D. Density-functional Exchange-energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A: At., Mol., Opt. Phys. 1988, 38, 3098−3100. (81) Brena, B.; Carniato, S.; Luo, Y. Functional and Basis Set Dependence of K-edge Shake-up Spectra of Molecules. J. Chem. Phys. 2005, 122, 184316. (82) Clark, T.; Chandrasekhar, J.; Spitznagel, G. W.; Schleyer, P. V. R. Efficient Diffuse Function-Augmented Basis Sets for Anion Calculations. III. The 3-21+G Basis Set for First-row Elements, Li− F. J. Comput. Chem. 1983, 4, 294−301. (83) Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. Universal Solvation Model Based on Solute Electron Density and on Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions. J. Phys. Chem. B 2009, 113, 6378−6396. (84) Chatterjee, S.; Sengupta, K.; Mondal, B.; Dey, S.; Dey, A. Factors Determining the Rate and Selectivity of 4e−/4H+ Electrocatalytic Reduction of Dioxygen by Iron Porphyrin Complexes. Acc. Chem. Res. 2017, 50, 1744−1753. (85) Huang, X.; Groves, J. T. Oxygen Activation and Radical Transformations in Heme Proteins and Metalloporphyrins. Chem. Rev. 2018, 118, 2491−2553. (86) For a sole example of a CoII/H2O2 adduct, see: Wallen, C. M.; Palatinus, L.; Bacsa, J.; Scarborough, C. C. Hydrogen Peroxide Coordination to Cobalt(II) Facilitated by Second-Sphere Hydrogen Bonding. Angew. Chem., Int. Ed. 2016, 55, 11902−11906. (87) Nguyen, A. I.; Hadt, R. G.; Solomon, E. I.; Tilley, T. D. Efficient C−H Bond Activations via O2 Cleavage by a Dianionic Cobalt(II) Complex. Chem. Sci. 2014, 5, 2874−2878. (88) Hong, S.; Pfaff, F. F.; Kwon, E.; Wang, Y.; Seo, M.-S.; Bill, E.; Ray, K.; Nam, W. Spectroscopic Capture and Reactivity of a Low-Spin Cobalt(IV)−Oxo Complex Stabilized by Binding Redox-Inactive Metal Ions. Angew. Chem., Int. Ed. 2014, 53, 10403−10407. (89) Wallen, C. M.; Bacsa, J.; Scarborough, C. C. Coordination of Hydrogen Peroxide with Late-Transition-Metal Sulfonamido Complexes. Inorg. Chem. 2018, 57, 4841−4848. (90) Nurdin, L.; Spasyuk, D. M.; Fairburn, L.; Piers, W. E.; Maron, L. Oxygen−Oxygen Bond Cleavage and Formation in Co(II)Mediated Stoichiometric O2 Reduction via the Potential Intermediacy of a Co(IV) Oxyl Radical. J. Am. Chem. Soc. 2018, 140, 16094−16105. (91) For the complexes in Figure 1a, values of E1/2 − EORR(H2O2) are as follows: a, +70 mV; b, −100 mV; c, −190 mV; d, −200 mV; e, −20 mV (pH = 4); f, −220 mV (Ar′ = p-C6H4Cl); g (Figure 1b), ECo(III/II) − EORR(H2O2) = −260 mV.

1034

DOI: 10.1021/acscentsci.9b00194 ACS Cent. Sci. 2019, 5, 1024−1034