Open Access Article
Simone
Bertini
a,
Motiar
Rahaman
a,
Abhijit
Dutta
a,
Philippe
Schollhammer
b,
Alexander V.
Rudnev
ac,
Fredric
Gloaguen
*b,
Peter
Broekmann
*a and
Martin
Albrecht
*a
aDepartment of Chemistry, Biochemistry &Pharmacy, Universität Bern, Freiestrasse 3, 3012 Bern, Switzerland. E-mail: martin.albrecht@dcb.unibe.ch; peter.broekmann@dcb.unibe.ch
bUMR 6521, CNRS, Université de Bretagne Occidentale, CS 93837, 29238 Brest, France. E-mail: frederic.gloaguen@univ-brest.fr
cA.N. Frumkin Institute of Physical Chemistry and Electrochemistry Russian Academy of Sciences, Leninskii pr. 31, 119071 Moscow, Russia
First published on 12th April 2021
Strategies for the conversion of CO2 to valuable products are paramount for reducing the environmental risks associated with high levels of this greenhouse gas and offer unique opportunities for transforming waste into useful products. While catalysts based on nickel as an Earth-abundant metal for the sustainable reduction of CO2 are known, the vast majority produce predominantly CO as a product. Here, efficient and selective CO2 reduction to formate as a synthetically valuable product has been accomplished with novel nickel complexes containing a tailored C,O-bidentate chelating mesoionic carbene ligand. These nickel(II) complexes are easily accessible and show excellent catalytic activity for electrochemical H+ reduction to H2 (from HOAc in MeCN), and CO2 reduction (from CO2-saturated MeOH/MeCN solution) with high faradaic efficiency to yield formate exclusively as an industrially and synthetically valuable product from CO2. The most active catalyst precursor features the 4,6-di-tert-butyl substituted phenolate triazolylidene ligand, tolerates different proton donors including water, and reaches an unprecedented faradaic efficiency of 83% for formate production, constituting the most active and selective Ni-based system known to date for converting CO2 into formate as an important commodity chemical.
While several complexes based on Earth-abundant Mn, Fe, and Ni metals have been known to catalyze the reduction of CO2,20–25 the vast majority of these catalysts produce CO as a predominant product.26,27 Only two Ni systems have been reported to yield formate,28 namely Sauvage's Ni(cyclam) system from over 30 years ago,29–32 and Fontecave's Ni(III) catalyst,33 though selectivity is a major issue in both systems due to significant formation of CO. Formate formation is highly desirable as it constitutes a pathway to convert waste to a valuable product for synthesis, hydrogen storage, formic acid fuel cells, and for other industrial uses.
Inspired by the work of Kirchner and others on manganese pincer complexes,34–36 which demonstrated a key relevance of the metal-hydride intermediate to promote CO2 insertion rather than direct CO2 bonding, we became interested in exploiting the potential of triazolylidenes37,38 as a specific subclass of N-heterocyclic carbene (NHC) ligands39–41 for imparting such reactivity. In order to increase the robustness of the M–CNHC bond and hence the reliability of the carbene as a spectator ligand to the nickel active site,42 oxygen chelating groups were introduced.43 Here we show that this approach provides a set of new, tunable, and highly active catalysts for the electrochemical reduction of CO2 to formate. Catalyst screening in half-cell measurements revealed that the most active system accomplishes unrivalled faradaic efficiency and outstanding selectivity towards formate, largely outperforming currently known catalyst systems based on nickel.
The 1H NMR spectra of the complexes reveal the presence of only one isomer. While complexes 3a,b and 4 feature the phenolate proton resonances in the expected aromatic region (Fig. S20†), the spectrum of complex 3c is distinct with a markedly upfield shifted phenolate ortho proton (δH = 4.8; Fig. S20†). This shift of more than 2 ppm suggests a trans arrangement of the two C,O-bidentate ligands, with the phenolate Hortho influenced by ring current anisotropy of the mesityl group of the other ligand. No such effect has been observed in the NMR spectrum of complex 3b featuring a phenyl-substituted triazolylidene, which suggests a cis configuration. These ligand arrangements were unambiguously confirmed by single crystal X-ray diffraction (Fig. 1). In the trans configuration of 3c, the mesityl ring is essentially orthogonal to the carbene–Ni plane and close to Hortho (Hortho⋯arene 2.75 Å), which accounts for the ring anisotropy observed by NMR spectroscopy. The cis configuration of the other complexes places the carbene wingtip groups R in close proximity, which results in a substantial distortion from square-planar to tetrahedral (cf. τ4 ≥ 0.15).44 The distortion is larger for bulky wingtip groups (Bu, Ph, in 3a,b, respectively) than with methyl substituents (4), and negligible in the trans complex 3c (τ4 < 0.01; Table S7†).
All four nickel complexes 3a–c, 4 display (quasi)reversible redox processes around +0.7 V and −2.0 V, tentatively attributed to NiII/NiIII and NiII/NiI transitions, respectively (Fig. 2, S21 and 22;† potentials vs. Ag/AgCl).33 Comparison of the redox potentials consistently indicates that triazolylidenes are stronger donor ligands than imidazolylidene,45 and that the wingtip group R directly affects the electron density on the nickel center with oxidation potentials increasing along the series R = Bu (E1/2 = 0.59) < R = Mes (E1/2 = 0.64) < R = Ph (E1/2 = 0.69; Table 1).
| Entry | Complex |
E
pc b |
E 1/2(NiII/III) b |
k
obs
[s−1] |
|---|---|---|---|---|
| a All values in V vs. Ag/AgCl; 1 mM MeCN solution of Ni complex with (Bu4N)PF6 as supporting electrolyte, 250 mV s−1 scan rate, Fc+/Fc as internal standard with E1/2 = 0.36 V vs. Ag/AgCl. b Cathodic peak potential Epc for NiII/NiI reduction and half-wave potential E1/2 (in parentheses ΔEp = Epa − Epc in mV for NiII/NiIII redox process) c Catalytic conditions; 1 mM complex in MeCN, AcOH (0.8 M), (Bu4N)PF6 as supporting electrolyte, kobs determined by foot of the wave analysis (see ESI for details†). | ||||
| 1 | 3a | –2.16 | 0.59 (130) | 440 |
| 2 | 3b | –2.09 | 0.69 (110) | 300 |
| 3 | 3c | –2.12 | 0.64 (120) | 10 |
| 4 | 4 | –1.92 (110) | 0.75 (120) | 200 |
The electrocatalytic performance of complexes 3a–c, 4 was first investigated in H+ reduction. A significant cathodic current was observed when acetic acid (AcOH) was present as proton donor in a MeCN solution of the Ni complex. Increasing the AcOH concentration from 5 to 400 eq. with respect to the Ni complex led to an enhanced current density, indicative of catalytic H+ reduction (Fig. 3). Extraction of kobs by foot-of-the-wave analysis (FOWA)46,47 reveals a direct influence of the steric and electronic properties of the ligand in promoting catalytic reduction (Table S8†). Specifically, the trans arrangement imposed by the very bulky mesityl wingtip groups is strongly deactivating, while all the cis-complexes are active. Moreover, the activity of the cis-complexes directly correlates with the ligand donor properties deduced from CV data: the alkyl-substituted triazolylidene induces more than a twice higher active than the analogous imidazolylidene (kobs = 440 s−1 for 3avs. 200 s−1 for 4).48
Prompted by the promising catalytic activities, complexes 3a–c, 4 were evaluated as catalyst precursors for electrochemical reduction of CO2. Initial experiments with a 1 mM MeCN solution of the nickel complex 3a under a CO2 atmosphere reveal a significant enhancement of the cathodic current upon complex reduction, indicative of CO2 transformation (Fig. 4). For complex 4 the reversibility of the NiII/NiI reduction was lost upon saturation with CO2 gas. All the complexes were active in the process, with only the trans isomer 3c induces low catalytic current, suggesting lower activity. The catalytic current enhances further when the reaction was performed in the presence of MeOH (40 eq. with respect to the Ni complex), suggesting a beneficial role of proton sources.49 Trifluoroethanol and phenol show similar effects, though the current increase is largest when using MeOH (Fig. 4). Blank measurements indicate no catalytic current with CO2-saturated solutions in the absence of the Ni complex, or when the complex is reduced in the presence of methanol yet without CO2 (Fig. S34†).
The robustness of the catalytically active species over time was investigated by chronoamperometry at −1.9 V vs. Ag/AgCl in MeOH/MeCN 1
:
50 v
:
v (Fig. S26–S29†). The observed catalytic current is constant over 2 h, suggesting no significant degradation during that time. More extended reaction reveals a gradual decrease of activity. Comparison of the different complexes reveals catalytic activity for CO2 reduction follows the same trends observed for H+ reduction, with highest rates for the nickel complex with the strongest donating triazolylidene ligand (3a, kobs = 280 s−1), which is essentially twice as fast as the corresponding imidazolylidene analogue 4 (kobs = 150 s−1; Table 2). Again, the cis ligand arrangement is essential for ensuing catalytic activity as the nickel complex 3c with the ligands in trans configuration is almost inactive (kobs = 10 s−1), and enhanced electron density at the nickel center increases catalytic activity (3a > 3b > 4).
| Entry | Complex | FEHCOO- (4 h) [%] | FEHCOO- (8 h) [%] | FEH2+CO (4 h) [%] |
k
obs
[s−1] |
|---|---|---|---|---|---|
a General conditions: 1 mM complex, at −1.9 V vs. Ag/AgCl, glassy carbon working electrode and Pt foil as counter electrode (see ESI for details†) in MeOH/MeCN 1 : 50 v : v with 0.1 M (Bu4N)PF6 as supporting electrolyte.
b Determined from foot of the wave data treatment (see ESI for details†).
|
|||||
| 1 | 3a | 54 | 68 | 3 | 280 |
| 2 | 3b | 43 | 47 | 4 | 220 |
| 3 | 3c | 10 | 10 | 2 | 10 |
| 4 | 4 | 25 | 25 | 3 | 150 |
Product identification focused first on gas-phase analysis of volatiles products as most hitherto known Ni catalysts for CO2 reduction produce CO.10,11 Remarkably, only traces of H2 were detected in the gas phase, and CO quantities were below the detection limit (less than 2% faradaic efficiency after 4 h, Table 2). Analysis of the solution phase by ion-exchange chromatography (IC), HPLC and NMR measurements identified formate as the principal product of CO2 reduction with complexes 3a–c, 4. Complex 3a does not display only the highest activity but also imparts the highest selectivity. Optimization of reaction times affords faradaic efficiencies up to 70% for formate formation, one of the highest known so far for homogeneous Ni-based electrocatalysts.10,11
Ligand tailoring has been used to further improve the catalytic activity of these nickel complexes, in particular through introduction of electron-donating substituents on the phenolate. Substitution of the 4,6-positions moreover avoids undesired radical reactions of the phenolate, which may infer from one-electron-reduction upon catalyst activation. To this end the 4,6-Me2-phenol-triazolium salt 5a was nickelated via the procedure established for the synthesis of 3a, while the tBu analogue 5b required nBuLi as a stronger base to accomplish nickel complexation (Fig. 5a). The new O,C-bidentate chelated trz Ni(II) complexes 6a and 6b showed NMR characteristic reminiscent of those of 3a,b, indicating the formation of the cis-isomers exclusively. This structural assignment was further confirmed by single crystal X-ray diffraction studies of complexes 6a,b (Fig. 5b). The distortion from square planar geometry in these complexes is slightly less than in the parent complex 3a with τ4 values of 0.15 and 0.18 for 6a and 6b, respectively (cf. 0.21 for 3a, Table S7†).
The electrochemical properties of complexes 6a,b are comparable to those of complexes 3a–c, featuring a quasi-reversible oxidation at 0.49 and 0.45 V vs. Ag/AgCl, respectively, and an irreversible reduction below −2.0 V (Fig. S21 and 22†). The incorporation of donating groups on the phenolate lowers the redox potentials by about 90 (R = Me, 6a) and 150 mV (R = tBu, 6b relative to the parent phenolate complex 3a (Table 3). These potential shifts indicate enhanced electron density at the metal center and efficient electronic tailoring of the nickel center by ligand modifications.
| Entry | Complex |
E
pc b [V] |
FEHCOO- (8 h) [%] | FEH2+CO (4 h) [%] |
k
obs
[s−1] |
|---|---|---|---|---|---|
a General conditions: 1 mM complex, at −1.9 V vs. Ag/AgCl, glassy carbon working electrode and Pt foil as counter electrode (see ESI for details†) in MeOH/MeCN 1 : 50 v : v with 0.1 M (Bu4N)PF6 as supporting electrolyte.
b Cathodic peak potential Epc for NiII/NiI reduction in V vs. Ag/AgCl; 1 mM MeCN solution of Ni complex with (Bu4N)PF6 as supporting electrolyte, 250 mV s−1 scan rate, Fc+/Fc as internal standard with E1/2 = 0.36 V vs. Ag/AgCl.
c Determined from foot of the wave data treatment (see ESI for details†).
|
|||||
| 1 | 6a | –2.23 | 74 | 4 | 300 |
| 2 | 6b | –2.31 | 83 | 3 | 370 |
In these cis triazolylidene Ni complexes, lower reduction potentials are correlated with higher electrocatalytic CO2 reduction activity (cf.Table 2) Indeed, the stronger donor ligands in complexes 6a,b and the associated low reduction potential improved the catalytic performance. Foot of the wave data analysis revealed a higher kobs for CO2 reduction for 6a than for the parent complex 3a (300 vs. 280 s−1; Table 3), and more significantly for 6b with the lowest reduction potential (kobs = 370 s−1). Moreover, the robustness of the catalytically active species is demonstrated by a constant catalytic current over the first 2 h (Fig. 6 and S31†). Formate was identified as exclusive product of CO2 reduction also with complexes 6a,b with only traces of H2 and CO as side products. Complex 6b revealed the highest FE of 83% for formate formation and represents a new benchmark for Ni in this transformation. Moreover, these results demonstrate that ligand tailoring on a molecular level through incorporation of electron-donating substituents provides an efficient strategy for catalyst optimization for this CO2 reduction and leads in this case to a marked increase of the FE to unprecedented 83%.
While this work is focused on catalyst optimization for the CO2 reduction half-reaction, obviously also the availability of protons is elementary as the efficiency of the CO2 to formate transformation is associated with a coupled electron/proton transfer. Variation of the proton source has been accomplished by varying the additive to the aprotic solvent (MeCN) from MeOH to different proton donors including PhOH, iPrOH, CF3CH2OH (TFE), and H2O (Fig. 7 and S32†). All these proton sources are mediating CO2 reduction, yet with variable efficiency. Notably, water as the cleanest proton source is tolerated and affords considerable yields of formic acid. Best performance was achieved with TFE, reaching FEs up to 82% within 2 h. These data indicate that formate is exclusively formed at the cathode and not anodically, e.g. by partial MeOH oxidation, and that MeOH is not required to achieve high formate yields. Moreover, the broad tolerance of a variety of proton sources provides ample opportunities for optimizing also the oxidation half-reaction for designing an efficient full cell electrolyzer.
![]() | ||
| Fig. 7 Comparison of the faradaic efficiencies (FE) for formate in a CO2-saturated solution containing 1 mM complex 6b with 0.5 M of different proton donors (in MeCN containing 0.1 M (Bu4N)PF6 as supporting electrolyte). Electrocatalytic CO2 reduction was performed at −1.9 V vs. Ag/AgCl for 2 h using glassy carbon working electrode; the formate yield was quantified by ion exchange chromatography (Fig. S32†). | ||
Total FE values below 100% suggest parasitic side reactions and catalyst deactivation, which might be triggered by several factors, including acidification of the reaction medium due to the generated formate, (cf. H+ reduction with high HOAc concentrations above). This limitation should be easily mitigated, for example, by using a flow reactor. Catalyst deactivation during CO2 electroreduction was evidenced by a combination of high-resolution scanning electron microscopy (HR-SEM), energy dispersive X-ray spectroscopy (EDX) and X-ray diffraction (XRD) analyses with an electrode containing complex 6b after 2 h of operation at −1.9 V vs. Ag/AgCl. Post-electrolytic HR-SEM and XRD analyses identified Ni oxide nanoparticles on the working electrode surface, while EDX suggests formation of films that are too thin for detection (Fig. S37†). Specifically, XRD revealed trace 2θ intensities that are characteristic of NiO and Ni2O3 on the electrode surface (Fig. S38†). Nickel oxide formation may be rationalized by partial reduction of the complexed nickel to Ni0 under the cathodic potential applied during electrolysis, which induces complex decomposition and formation of Ni nanoparticles. These electrochemically generated Ni nanoparticles were probably transformed to Ni oxides (NiO/Ni2O3) only when the electrode was exposed to air after the experiment. Such partial complex degradation thus provides a plausible rationale for the gradual loss of activity during extended reaction times and may account for the uncompensated charge from chronoamperometric experiments as some cathodic charge is involved in the in situ reduction of Nin+ to Ni0 during the CO2 reduction process. It must be noted, however, that Ni nanoparticles are unable to catalyze CO2 reduction and only induce H2 formation (HER), independent of size and shape of these nanoparticles.50 Therefore the highly selective reduction of CO2 to formate as observed here is confidently attributed to the distinct catalytic activity of the molecularly defined NHC Ni complexes.
The production of formate as valuable product from CO2 reduction is very rare for nickel-based catalysts, as most Ni systems produce CO.9 The best complex of the series, complex 6b shows high faradaic efficiency and exquisite formate vs. CO selectivity. The two other known Ni systems generating formate are less efficient in comparison and produce considerable amounts of CO as by-product (Table S11†).10,11 The triazolylidene system presented here offers unique potential, due to the intrinsically high selectivity towards formate production, and because the catalytic rate can be further optimized due to the correlation of rational ligand design and reduction potentials with catalytic activity. Moreover, a variety of proton sources are tolerated, including water. The largest drawback is probably the high overpotential required in comparison to the Kubiak-Sauvage electrocatalyst (ΔE = 0.5 V).11
The poor catalytic activity of the trans isomer points to a key role of the complex geometry. Either, the distortion from square-planar in the cis complexes may facilitate the formation of a penta-coordinate nickel(I) intermediate, or more likely, the cis-arrangement of the two oxygen donors produces a proton acceptor pocket, especially after one-electron reduction to either a formal nickel(I) complex or a phenoxide radical anion.51,52 Proton chelation by the two oxygen units is therefore suggested to form the reduced complex A (Scheme 2), presumably through a proton-coupled electron transfer (PCET). This intermediate facilitates the formation of a nickel(III) hydride intermediate B for CO2 insertion and generation of the formate complex C, which is then readily reduced to release the formate product. Formation of the nickel hydride intermediate is presumed to be key for the product selectivity,34–36,53–55 as most known nickel catalysts for CO2 reduction bind CO2 directly and therefore produce CO rather than formate.20,29–32,54–59 Alternatively, the proton scavenging may localize the negative charge on the oxygen, which facilitates the formation of the η1−OCO Ni adduct B′ as another potentially critical intermediate for formate formation.60 Attempts to characterize the putative nickel(I) intermediate by spectroelectrochemistry were not successful, as a solution of complex 3a or 4 did not produce any EPR signal after electrolysis at −2.0 V vs. Ag/AgCl, presumably because the nickel(I) species is not sufficiently stable in the absence of substrates. Nonetheless, this mechanistic proposal provides a rationale for the strong divergence of catalytic activity of the cis vs. trans complexes.
:
1). Single crystals suitable for X-ray diffraction analysis were obtained by slow diffusion of pentane into a CH2Cl2 solution of the complex.
:
2) gave complex 4 as a yellow crystalline powder (90 mg, 70%). 1H NMR (400 MHz, CD2Cl2): δ 7.33–7.12 (m, 3H, 2HPhO + 1Himi), 7.03 (t, J = 8.5 Hz, 1H, HPhO), 6.92 (s, 1H, Himi), 6.61 (t, J = 8.5 Hz, 1H, HPhO), 3.19 (s, 3H, CH3–N). 13C{1H} NMR (101 MHz, CD2Cl2): δ 159.22 (C–O), 157.62 (Cimi–Ni), 127.77 (CPhO–N), 127.62 (CPhO–H), 124.06 (Cimi–H), 122.42 (CPhO–H), 118.23 (Cimi–H), 117.57 (CPhO–H), 113.64 (CPhO–H), 37.01 (CH3–N). HR-MS (ESI): calcd for C20H17N4NaNiO2 [M + Na]+m/z = 426.0597 (found 426.0603). Anal. calcd for C20H17N4NiO2 (404.08): C, 59.45; H, 4.24; N, 13.87. Found: C, 59.81; H, 3.95; N, 14.01.
:
1) to afford complex 6a as a bright yellow crystalline solid (100 mg, 71%).1H NMR (400 MHz, CD2Cl2): δ 7.27 (s, 1H, HPh), 6.82 (s, 1H, HPh), 3.88 (s, 3H, N–CH3), 2.18 (s, 3H, CH3), 2.16 (s, 3H, CH3), 1.91 (s, b, 1H, CH–Pr), 1.81 (s, b, 1H, CH–Pr), 1.56 (s, b, 1H,CH–Et), 1.43 (s, b, 1H, CH–Et), 1.23–1.14 (m, 2H, CH2–CH3), 0.66 (t, J = 4 Hz, CH3). 13C{1H} NMR (101 MHz, CD2Cl2): δ 156.33 (C–O), 146.28 (Ctrz–Ph), 142.68 (Ctrz–Ni), 131.26 (CPh–H), 130.67 (CPh–H), 126.21 (CPhO–N), 121.10 (CPh–H), 117.13 (CPh–H), 36.47 (CH3–N), 32.25 (CH2–Pr), 24.99 (CH3), 22.98 (CH2–Et), 20.74 (CH3), 17.15 (CH2–CH3), 13.72 (CH3). HR-MS (ESI): calcd for C30H40N6NiO2 [M + H]+m/z = 575.2639 (found 575.2638). Anal. calcd for C30H24N6NiO2 (575.38): C, 62.62; H, 7.01; N, 14.61. Found: C, 62.33; H, 7.12; N, 14.54.
:
5) to obtain complex 6b as a dark yellow powder (280 mg, 62%). 1H NMR (400 MHz, CD2Cl2): δ 7.44 (d, J = 2.1 Hz, 1H, HPhO), 7.30 (d, J = 2.1 Hz, 1H, HPhO), 4.23 (s, 3H, CH3), 2.34 (t, J = 6.0 Hz, 2H, CH2–C3H7), 1.67 (s, b, 1H, CH2–C2H5), 1.47 (s, b, 1H, CH2–C2H5), 1.37–1.29 (m, 10H, (CH3)3–C + CH2–CH3), 1.28–1.24 (m, 10H, (CH3)3–C + CH2–CH3) 0.76 (t, 3H, J = 6.0 Hz, CH3). 13C{1H} NMR (100 MHz, CD2Cl2): δ 154.24 (C–O), 145.50 (CPh), 141.85 (CPh), 138.89 (CTrz–Ni), 128.53 (CPh–H), 126.22 (CPh–N), 122.42 (CTrz–Bu), 121.58 (CPh–H), 38.89 (C(CH)3), 36.19 (C(CH)3), 34.73 (CH3–N), 31.70 (CH3), 29.23 (CH3), 28.52 (CH2–C3H7), 23.49 (CH2–C2H5), 23.25 (CH2–CH3), 13.77 (–CH3). HR-MS (ESI): calcd for C42H64N6NiO2 [M + H]+m/z = 743.4517 (found 743.4495). Anal. calcd for C42H64N6NiO2 (743.70): C, 67.83; H, 8.67; N, 11.30. Found: C, 67.76; H, 8.88; N, 11.35.
Footnote |
| † Electronic supplementary information (ESI) available: Synthesis of all new triazolium salts, NMR spectra and cyclic voltammograms (CVs) of new compounds, NMR and CV experiments in acidic media, details on X-ray diffraction and on electrocatalytic set-up and product identification (PDF), crystallographic data for complexes 3a–c, 4, 6a–b (CIF). CCDC 2004181–2004184, 2050371 and 2050372. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d1gc00388g |
| This journal is © The Royal Society of Chemistry 2021 |