Gaia
Neri
,
Iain M.
Aldous
,
James J.
Walsh
,
Laurence J.
Hardwick
and
Alexander J.
Cowan
*
Department of Chemistry, Stephenson Institute for Renewable Energy, The University of Liverpool, UK. E-mail: a.j.cowan@liverpool.ac.uk
First published on 24th November 2015
The development of selective electrocatalysts for CO2 reduction in water offers a sustainable route to carbon based fuels and feedstocks. However, molecular catalysts are typically studied in non-aqueous solvents, in part to avoid competitive H2 evolution. [Ni(cyclam)]2+ (1) is one of the few known electrocatalysts that operate in water and 30 years after its report its activity remains a rarely surpassed benchmark. Here we report that [Ni(cyclam-CO2H)]2+ (cyclam-CO2H = 1,4,8,11-tetraazacyclotetradecane-6-carboxylic acid (2)) shows greatly enhanced activity versus1 for CO production. At pHs < pKa of the pendant carboxylic acid a large increase in catalytic activity occurs. Remarkably, despite the high proton concentration (pH 2), 2 maintains selectivity for CO2 reduction and is believed to be unique in operating selectively in such acidic aqueous solutions.
The mechanism for the reduction of CO2 to CO by 1 has been extensively studied,5,10–13 and although the exact nature of the active species has yet to be unambiguously identified, several factors have been made clear. Firstly, [Ni(cyclam)]+ adsorbs on to some metal electrodes including Sn, Pb,14 and Hg,5 and adsorption onto the electrode is key for efficient CO2 reduction.111 has also been shown to act as a homogeneous CO2 reduction catalyst when used with a glassy carbon electrode (GCE)15 however the level of activity was significantly lower than can be achieved on Hg, which in part may be due to suppression of catalyst degradation pathways on Hg.16 Indeed Hg remains a common electrode for fundamental studies such as that presented here. At pH 5 adsorption initiates at potentials positive of the formal NiII/I couple in solution (−1.3 VNHE) and a monolayer is formed at ca. −1 VNHE. The adsorbed NiI complex is predicted to bind in an η1-C mode to CO2, prior to the transfer of a second electron to the catalyst centre.13 Computational studies13,17,18 indicate that the structure of the adsorbed complex is a trans-I conformation,19 with the axial amine hydrogens aiding CO2 binding. In solution 1 can adopt five different isomeric forms11 and for clarity structures of 1 and 2 are drawn in a planar geometry (Fig. 1). Numerous studies have proposed that the decline in activity upon cyclam modification is due to conformational changes limiting the availability of the N–H group suggesting their critical role in catalysis.5,13,17 Finally, 1 is most active at pH ∼4–5 indicating that dissolved CO2 and not HCO3− or CO32− is the preferred substrate. At pH values less than 4H2 evolution dominates and CO2 selectivity is lost.9
Only a limited number of other classes of molecular CO2 reduction electrocatalysts for use in water are known.6,20–25 Of relevance are recent studies on organic catalysts including mercaptopteridine,26 an iridium pincer catalyst24 and very recently a water soluble iron porphyrin catalyst, labelled WSCAT,27 which preliminary data suggests is an extremely active catalyst at pH 6.7, although at lower pH values only H2 was produced. The limited pH range appears to be typical of CO2 reduction catalysts with most being studied between pH 6 and 7. In addition to WSCAT, Savéant et al.28,29 have also extensively studied other iron porphyrins for use in DMF. In an important breakthrough, a large increase in electrocatalytic activity for CO2 reduction to CO by an iron porphyrin modified with phenolic groups in DMF + 2.0 M H2O was reported.22 The acidic phenol groups on the catalyst framework acted as both a local proton source and to aid CO2 binding, greatly accelerating the proton coupled reduction of CO2, and similar approaches have now been employed by several groups studying a range of transition metal electrocatalysts for use in non-aqueous solvents.30–32 While these studies show that the addition of acidic groups can greatly accelerate the rate of reduction of CO2 in DMF, they have not been applied to catalysts that are active in water. Here we demonstrate that the modification of 1 with a carboxylic acid leads to a step change in catalytic activity in water with a five-fold increase in the observed rate constant (kobs), the turnover frequency per adsorbed catalyst, for 2 compared to 1, at −0.99 VNHE, near the foot of the catalytic wave. We also note an extremely high kobs = 3.4 (±1.0) × 103 s−1 at −1.25 VNHE. Perhaps most remarkable is that catalyst 2 operates in acidic conditions whilst maintaining selectivity towards CO2.
Bulk electrolysis experiments (Fig. S7†) confirm that 2 remains selective towards CO production at pH 2 and that its activity exceeds that of the parent catalyst 1. 5.2 ± 0.3 C of charge is passed within 1 hour during the electrolysis of an unstirred solution (3 μM) of 2 at −0.99 VNHE with a very good selectivity towards CO production >4:1 CO:H2 (faradaic efficiency (FE), total = 81%, H2 = 15 ± 5%, CO = 66 ± 9%, with errors being the result of 3 experiments). This corresponds to an average bulk turnover number of 591 for CO in 1 hour. In contrast 1 passes only 2.0 ± 0.2C in 1 hour with a lower selectivity 0.2:1 CO:H2 (FE, total = 86%, H2 = 73 ± 16%, CO = 13 ± 10%), and a bulk turnover number of ca. 45 for CO production. No liquid phase products were detected by NMR. [Ni(cyclam)]2+ and its derivatives are known to form inactive species in the presence of CO,34 however activity can be maintained through constant CO2 purging and experiments with 2 over a 7.5 hour period show activity being maintained, Fig. S8.†
To understand the factors behind the enhanced activity of 2 at low pH we have examined the electrochemical response of 1 and 2 over a wide pH (6–2) range using rotating disc electrode (RDE) voltammetry, (Fig. 1c and d), differential pulse voltammetry (DPV, Fig. S2 and S3†) and CV measurements (Fig. S4–6†). RDE measurements are employed to study the catalysis under CO2 as they minimise the effects of substrate diffusion and product inhibition, simplifying the analysis of the electrochemical response. Between pHs 6 and 4 RDE measurements of 2 under CO2 show only a slight increase in plateau current density, Fig. 1d. Between pH 3 and 2 a dramatic change is noted with a new reductive feature (ca. −0.95 V) growing in under CO2 as the pH decreases, which is shown above to be due to catalytic CO2 reduction. This leads to a large decrease in the potential necessary for catalysis between pH 5 and 2 of ca. 240 mV versus the normal hydrogen electrode. In the RDE measurements we define the potential necessary for catalysis as being when the current density exceeds 2 mA cm−2.35 In contrast with 1 we only measure a very small shift (ca. 50 mV) in the potential necessary for catalysis between pH 6 and 2, which will be at least in part due to the increased level of H2 production at low pHs. By pH 2 there is minimal separation of the RDE curves of 1 in the presence and absence of CO2, Fig. 1c. This step change in behaviour of 2 but not 1 is indicative of a change in catalytic mechanism for 2 between pH 3 and 2. Furthermore whilst the variation in overpotential for CO2 reduction brought about by a change in pH (0.18 V) is equivalent for both catalysts only 2 shows a significant change in potential necessary for catalysis. The lack of a pH dependence for 1 is further explored in the ESI (Fig. S2 and S3†) where we demonstrate that the NiII/I couple under argon is independent of pH.
In order to assess if the change in current density under CO2 with pH is due to the protonation of the carboxylic acid of 2 we have measured the pKa of this group by Fourier-Transform Infrared (FTIR) spectroscopy in solution (Fig. S9†). The spectra were recorded in a 0.1 mm path length CaF2 IR cell. The initially synthesised catalyst is prepared in basic conditions and the carboxylate has νas(CO2−) at 1575 cm−1 and νs(CO2−) modes at 1375 cm−1 in line with literature reports for similar complexes.36 Titration of a 0.1 M solution of 2 in D2O (initial pD = 9.61) with DCl showed the clear emergence of the carboxylic acid form of 2 with νas(CO) at 1706 cm−1 in D2O, with pKa ∼ 2.6. Deuterated solvents are required to avoid the δ(HOH) mode of H2O masking the spectral window of interest. There is an excellent correlation between the relative concentration of the protonated carboxylic acid in solution and the current density for 2 measured under CO2 using RDE at −0.99 VNHE (Fig. 2b) and −1.1 VNHE (Fig. S10†). This clearly shows that the enhancement in catalytic activity towards CO2 reduction is due to the availability of the protonated carboxylic acid group of 2 at low pH values. In contrast a similar pH titration of catalyst 1 shows no clear changes in the spectral and pH region studied.
It has been shown that for 1 the active CO2 reduction catalyst is adsorbed onto Hg electrodes.11 It is therefore important to ascertain if the active form of 2 is also an adsorbed species. The current density under CO2 of 2 on a GCE is found to be significantly lower than that measured on a HMDE (Fig. S11 and S12†) suggesting that the active catalyst is indeed surface adsorbed 2. The surface concentration of 2 on the HMDE electrode has been measured using double-potential-step chronocoulometry11 (Tables S1–S3†) and is found to be 2.0 (±0.2) × 10−10 mol cm−2 at −0.99 VNHE at pH 5, similar to that previously reported for 1, 1.6 (±0.2) × 10−10 mol cm−2.34 At pH 2 the surface concentration of both 1 and 2 are changed by a minimal amount (2.2 (±0.2) × 10−10 mol cm−2 and 1.6 (±0.2) × 10−10 mol cm−2 respectively at −0.99 VNHE), indicating that the large increase in activity of 2 cannot be attributed to a change in the surface coverage of the catalyst with pH.
The kinetic behaviour of 1 and 2 at pH 2 can be obtained from the RDE measurements carried out at different rotation rates (Fig. S15†). We calculate the kinetic activity of the catalyst from voltammetric sweep measurements as it has recently been highlighted that turnover frequencies obtained from long-term bulk electrolysis measurements at high current densities can be complicated by a range of factors including (i) substrate diffusion, (ii) product inhibition and (iii) catalyst deactivation.37
Using the limiting current obtained from the intercepts of Koutecký–Levich plots (Fig. S14 and S15†) we can obtain an apparent first order rate constant (kobs, s−1), i.e. the turnover frequency (TOF) per adsorbed catalyst using eqn (1).38
icat = nFΓAkobs | (1) |
However the most significant feature of 2 is its selectivity towards CO2 even under acidic conditions. All previously reported derivatives of 1 have shown predominantly hydrogen production outside of a small pH window7,11 and we note that the majority of CO2 reduction catalysts are reported at pHs close to neutral (5–7),24,25,27 making the ability of 2 to operate at pHs as low as 2 unusual. The correlation between the current density under CO2 and the protonation state of the carboxylic acid group of 2 suggests that the protonation state of the catalyst is an important factor in the enhanced TOF, and hence selectivity towards CO2 of 2 in acidic solutions (Fig. 2b). It may be envisaged that protonation of the carboxylic acid group leads to 2 being more readily reduced to form the active NiI catalyst, however DPV studies indicate the NiII/I couple to be pH independent under argon, Fig. S2 and 3.† Alternatively previous studies have shown that the presence of a local proton source can accelerate CO2 reduction and it is viable that the acid group may also aid catalysis here.28,30–32 In the homogenous reduction of CO2 by 1 in acetonitrile a proton concentration dependent peak in the CV, similar to the feature observed by RDE (Fig. 1d) here at ca. −0.95 VNHE was reported.16 This peak was assigned to the reduction of a protonated CO2 adduct, with this proton dependent electron transfer becoming the rate limiting step in CO2 catalysis under certain conditions. It is feasible that the protonated carboxylic acid is acting as local proton source during the reduction of a CO2 adduct here. Such an interaction is geometrically feasible. The cobalt analogue of 2 has been reported for use in dye-sensitized solar cells with binding of the –CO2H group directly to the metal centre.36,39 However we do recognise that the empirical nature of the relationship in Fig. 2b does not provide direct evidence of the functional role of the carboxylic acid. We are currently also unable to discount the role of other potential mechanistic aspects including a possible change in structure of the adsorbed catalyst or change in the catalysts CO2 affinity and further mechanistic studies are currently underway.
The enhanced activity of 2 and the ability to maintain selectivity towards CO2 across a wide pH range are highly desirable traits. It is likely that robustness towards local pH fluctuations and the ability to be employed in electrolysers using proton exchange membranes will be advantageous for any practically applicable catalyst. However operating at pH 2 does have implications regarding the overpotential for catalysis. The potential necessary for catalysis with a current density of 2 mA cm−2 in Fig. 1 is ca. −0.9 VNHE, corresponding to an overpotential of ca. −0.65 V versus the apparent equilibrium potential for CO2 reduction to CO at pH 2. Whilst not dissimilar to other previous studies in aqueous solutions,9,27 it is significantly higher than typically required in solvents such as DMF, CH3CN and ionic liquids indicating that further improvements in molecular catalysts for use in aqueous solutions are still required.
Footnote |
† Electronic supplementary information (ESI) available: Including full experimental details, surface coverage measurements and supporting electrochemical measurements. See DOI: 10.1039/c5sc03225c |
This journal is © The Royal Society of Chemistry 2016 |