Alma
Arévalo
,
Enrique
Juárez-Francisco
,
Diego A.
Roa
,
Marcos
Flores-Alamo
and
Juventino J.
García
*
Facultad de Química, Universidad Nacional Autónoma de México, Circuito Interior, Ciudad Universitaria, Mexico City 04510, Mexico. E-mail: juvent@unam.mx
First published on 5th August 2024
This report includes the preparation of a new set of well-defined Cu(I) catalytic precursors of the type [Cu(diphosphine)(PPh3)NO3] and [Cu(diphosphine)NO3], fully characterized by regular analytical methods, including single-crystal XRD (X-ray diffraction). The new compounds were assessed to activate CO2 in an electrocatalytic process to yield oxalate selectively and with a relatively low overpotential. Some mechanistic insights into this process are also provided; oxalate is a valuable product for further chemical applications.
Since electrochemical methods can use renewable energy as the driving force for the chemical conversion of substances, acting on a redox-active molecule lowers energy demands like heating. In an electrocatalytic system, the reducing equivalents come directly from the electrode surface as electrons, which are transferred either by diffusion of the catalyst to the electrode surface or by catalyst adsorption onto the electrode surface.3,4
To date, there are some reviews about the electroreduction of CO2 and hydrogen sources to obtain formic acid, hydrocarbons, or alcohols through homogeneous catalysis using different transition metal complexes and typically using Pd, Ir, Ru, and more recently, Ni, Co, or Fe compounds, and macrocyclic ligands or bulky phosphines.5,6
Despite the wealth of reports about CO2 reduction leading to CO and HCO2H already published, a relatively small amount of them give oxalate as the main product since this involves CO2 reduction along with a C–C bond formation. Seminal papers by Nonaka7 describe the production of a mixture of formic acid, oxalic acid, and CO using CuCl2 and PPh3 with a faradaic efficiency (FE) of 73%. Tanaka8 outlines the selective oxalate generation in electrochemical CO2 reduction using triangular metal–sulfide clusters of Ir (FE 60%) and Co (FE 80%). Isobe9 describes the formation of C2O42− in the electrochemical reduction of CO2 using [(Ir(η5-C5Me5))3(μ3-S)2](BPh4)2 as a catalytic precursor with a FE of 60%. Another report by Tanaka10 and coworkers shows the almost selective oxalate formation in electrochemical reduction of CO2 catalyzed by mono- and di-nuclear Ru(bpy)2 compounds and unsymmetrical chelating ligand complexes that allowed to elucidate the inner sphere mechanism for the formation of oxalate (FE 70%).
Other examples include Jäger11 with a macrocyclic [N42−]Ni: (Ni-Etn-(Me/COOEt)Etn), that is a selective homogeneous catalyst for the electrochemical reduction of CO2 to oxalate. Wong12 outlines the use of an iron complex, [FeII(dophen)(N-MeIm)2]ClO4 (N-MeIm = 1-methyl-imidazole), as a catalyst for CO2 reduction, obtaining a mixture of carbon monoxide, formate and oxalate with a FE of 11%. Bouwman13 describes a binuclear Cu(I) complex that is oxidized in air by CO2, giving a tetranuclear Cu(II) complex containing two bridging CO2-derived oxalate groups, which are precipitated with lithium perchlorate as lithium oxalate, with an efficiency of 96%. Udugala-Ganehenege14 reports the electrochemistry of a hemicyclic Ni(II) complex that showed activity for CO2 reduction to oxalate, which was detected spectroscopically, but the FE is not disclosed.
The production of oxalate is relevant since it is the main component of lithium batteries (as lithium oxalate15) or derived to oxalic acid to be used as a solid lubricant for the separation/recovery of rare earth elements, metal treatment, bleaching agents,16 textile treatment, leather tanning, marble polishing or as an intermediate of pharmaceuticals or even as medicaments,17 agrochemicals and in organic synthesis.18 Some of the metals used are expensive (Ir and Rh), although they produce oxalate in good yields; the cobalt cluster has a high FE, but so far, Fe and Ni are not suitable catalysts for this process; they have low FE. Regarding earth-abundant metals, copper catalysts give the best FE, have low toxicity to humans, and are not expensive.
Considering the above, we turned our attention to the electroreduction of CO2 using Cu(I) catalysts that selectively produce oxalate in good yields at low overpotential, which is reported herein, along with the synthesis and characterization of four new copper compounds, three of them containing the fragment R2P(CH2)2–PR2, where R = iPr (dippe), R = Ph (dppe), R = ethyl (depe); and a closely related complex with 1,1-bis(diisopropylphosphino)ferrocene (dipf).
The different phosphine ligands (P-donor) were selected in the catalyst design based on their electronic (strong σ-donor abilities) and steric (strong binding capabilities owing to the chelate effect) properties; similar synthetic reports are known using N- and S-donor ligands.19,20 Additionally, the dipf ligand was used to investigate if there is any improvement in reactivity due to the presence of two metals, a redox ligand that enhances the copper function of the catalyst. Finally, nitrate was chosen based on its great capacity as a labile ligand to generate a vacant site at the catalyst.
Compound (2) was obtained as a white solid, washed with hexane, then with toluene, and finally dried in vacuum for 4 hours; an isolated yield of 80% was obtained after workup. Suitable single crystals for XRD were obtained for complex (2) from a saturated THF–toluene solution stored at −30 °C for 48 h in a drybox. The corresponding ORTEP (Oak Ridge Thermal Ellipsoid Program) plot for (2) is depicted in Fig. 1.
Fig. 1 ORTEP plot (50% probability) for complex (2). Labels of carbon atoms are omitted for clarity. |
Complex (2) is a diamagnetic d10 copper(I) complex suitable for NMR (nuclear magnetic resonance) study. The room temperature 31P{1H} NMR spectra for (2) displayed two signals at 0.45 ppm (PPh3) and 12.1 ppm (dippe) (Table 1, entry 1); both signals are broad singlets; thus, the coupling constants cannot be determined. Due to the air sensitivity of the dippe ligand complex, (2) was synthesized under an argon atmosphere; nevertheless, after preparation, (2) turned out to be air-stable as a solid for several days in an uncontrolled atmosphere; however, solutions of (2) are not stable and change from colorless to blue, characteristic of copper(II) compounds.
Entry | Complex | PPh3δP (ppm) | Biphosphine δP (ppm) |
---|---|---|---|
1 | [Cu(dippe)(PPh3)NO3] (2) | 0.45 | 12.1 |
2 | [Cu(depe)(PPh3)NO3] (3) | 4.45 | −4.5 |
3 | [Cu(dppe)(PPh3)NO3] (4) | 2.87 | −6.25 |
4 | [Cu (dipf)NO3] (5) | — | 3.13 |
As expected, upon coordination with the copper center, the chemical shift corresponding to each P-donor ligand moves to a lower field (Scheme 2). The 13C{1H} spectra of (2) show three singlets at the aliphatic region corresponding with the chelate phosphine and four signals for the aromatic rings of the PPh3 moiety. The signal at δ = 137.81 ppm corresponds to the ipso carbon of the phosphine phenyl ring, as shown in Fig. S4.†
Scheme 2 δ P values for copper compounds vs. δP for free diphosphines dippe, dppe,23 depe,24 and dipf,25 (PPh3 excluded). |
The IR (ATR-neat) data for the nitrate moiety of the new copper compounds are listed in Table 2.
Entry | Complex | NO sym. | NO asym. |
---|---|---|---|
1 | [Cu(dippe)(PPh3)NO3] (2) | 1299 | 1399 |
2 | [Cu(depe)(PPh3)NO3] (3) | 1297 | 1434 |
3 | [Cu(dppe)(PPh3)NO3] (4) | 1274 | 1434 |
4 | [Cu (dipf)NO3] (5) | 1286 | 1431 |
Suitable single crystals for XRD were also obtained for complex [Cu(depe)(PPh3)NO3] (3) and [Cu(dipf)NO3] (5) from saturated THF–hexane solutions stored at −30 °C under an argon atmosphere. The corresponding ORTEP plot for (3) is depicted in Fig. 2, and the ORTEP plot for (5) is in Fig. 3.
Fig. 2 ORTEP plot (50% probability) for complex (3). Labels of carbon atoms are omitted for clarity. |
Fig. 3 ORTEP plot (50% probability) for complex (5). Labels of carbon atoms are omitted for clarity. |
The XRD studies for (2), (3) and (5) complexes show that the copper atom is in a tetracoordinate environment. This coordination is frequently found in Cu(I) complexes; however, this is the first report of the [Cu(dippe)(PPh3)NO3] (2), [Cu(depe)(PPh3)NO3] (3) and [Cu(dipf)NO3] (5) crystal structures. Table 3 shows selected bond lengths [Å], angles [°], and torsion angles [°]. The sum of the internal bond angles26 centered at copper is 656.77° for (2), 657.08° for (3), and 643.09° for (5), and the tau(4)-descriptor for 4-coordination (τ4) values27 of 0.78 and 0.81 are evidence that (2) and (3) are close to the trigonal pyramidal geometry (TRP), while (5) shows a perfect TRP geometry with τ4 of 0.85. The puckering parameters of the five membered ring Cu(1)–P(1)–C(1)–C(2)–P(2) in (2), with Q(2) = 0.454 Å and φ(2) = 250.5(2)°, establish the closest pucker descriptor being enveloped on C1; while the same analysis for (3) results in Q(2) = 0.456 Å and φ(2) = 278.9(3)° that corresponds to a closest pucker descriptor being twisted on C(1)–C(2).
2 | 3 | 5 | |||
---|---|---|---|---|---|
Bond | [Å] | Bond | [Å] | Bond | [Å] |
Cu(1)–O(1) | 2.146(2) | Cu(1)–O(1) | 2.090(3) | Cu(1)–O(1) | 2.1598(15) |
Cu(1)–P(1) | 2.2944(9) | Cu(1)–P(1) | 2.2759(16) | Cu(1)–P(1) | 2.2305(6) |
Cu(1)–P(2) | 2.2904(9) | Cu(1)–P(2) | 2.2593(16) | Cu(1)–P(2) | 2.2383(6) |
Cu(1)–P(3) | 2.2688(9) | Cu(1)–P(3) | 2.2323(13) | Cu(1)–O(2) | 2.2580(16) |
Angle | [°] | Angle | Angle | [°] | |
O(1)–Cu(1)–P(3) | 97.27(7) | O(1)–Cu(1)–P(3) | 99.13(10) | O(1)–Cu(1)–P(1) | 122.00(5) |
O(1)–Cu(1)–P(2) | 111.01(7) | O(1)–Cu(1)–P(2) | 109.20(11) | O(1)–Cu(1)–P(2) | 113.29(4) |
P(3)–Cu(1)–P(2) | 124.54(3) | P(3)–Cu(1)–P(2) | 129.28(6) | P(1)–Cu(1)–P(2) | 117.96(2) |
O(1)–Cu(1)–P(1) | 108.62(7) | O(1)–Cu(1)–P(1) | 112.16(11) | O(1)–Cu(1)–O(2) | 58.44(6) |
P(3)–Cu(1)–P(1) | 125.08(3) | P(3)–Cu(1)–P(1) | 116.47(6) | P(1)–Cu(1)–O(2) | 117.31(4) |
P(2)–Cu(1)–P(1) | 90.25(3) | P(2)–Cu(1)–P(1) | 90.84(5) | P(2)–Cu(1)–O(2) | 114.09(5) |
Torsion angles | [°] | Torsion angles | [°] | Torsion angles | [°] |
P(1)–Cu(1)–O(1)–N(1) | 60.1(2) | P(1)–Cu(1)–O(1)–N(1) | −61.1(3) | P(1)–Cu(1)–O(1)–N(1) | −106.60(11) |
P(2)–Cu(1)–O(1)–N(1) | −37.6(2) | P(2)–Cu(1)–O(1)–N(1) | 38.0(3) | P(2)–Cu(1)–O(1)–N(1) | 102.84(12) |
P(3)–Cu(1)–O(1)–N(1) | −169.0(2) | P(3)–Cu(1)–O(1)–N(1) | 175.3(3) | P(1)–Cu(1)–O(2)–N(1) | 114.57(12) |
O(2)–N(1)–O(1)–Cu(1) | −19.4(4) | O(2)–N(1)–O(1)–Cu(1) | 11.1(6) | P(2)–Cu(1)–O(2)–N(1) | −101.44(12) |
For compound (2), probably the isopropyl group has a higher steric effect with triphenylphosphine, resulting in longer Cu(1)–P(1) and Cu(1)–P(2) bond distances (approx. 2.294 Å) compared to the average value reported (2.273 Å) at the Cambridge Crystallographic Data Centre (CCDC-2023). For compound (3), the distance Cu(1)–P(1) is almost identical to the average value reported (2.273 Å) in CCDC-2023. In both compounds, the oxygen atom of the nitrate group is bonded in a κ1-O coordination mode. In compound (5), the presence of the ferrocenyl group bonded between P(1) and P(2) opens the angle P(1)–Cu(1)–P(2), and increases the distance between the P atoms to 3.83 Å, causing a decrease in the Cu(1)–P distances to a value of 2.23 Å, similar to the average Cu–P value reported in the CDCD-2023 for P as a monodentate ligand. Additionally, the steric effects of ferrocene avoid the coordination of triphenylphosphine and favor the nitrate group to act as a bidentate κ2-O ligand to copper. Note that there is no interaction between the Fe and Cu atoms in (5). Selected bond lengths [Å] and angles [°] for (2), (3) and (5) are shown in Tables S2, S4, and S6,† respectively.
This study was performed at different scanning rates, and no additional oxidation or reduction peaks were observed (Fig. S32 and S33†). A similar study was assessed for all the copper complexes used here; the reduction values for Cu(I)/Cu(0) are somewhat similar and are summarized in Table 4.
Compound | E Cu(I)/Cu(0)Ar (V vs. Fc) |
---|---|
[Cu(PPh3)2NO3] (1) | −2.44 |
[Cu(dippe)(PPh3)NO3] (2) | −3.01 |
[Cu(depe)(PPh3)NO3] (3) | −2.89 |
[Cu(dppe)(PPh3)NO3] (4) | −2.97 |
[Cu (dipf)NO3] (5) | −2.91 |
As seen in Table 5, all the tested proton sources favored the HER process; the more acidic substances showed the best HER activity (PTSA and PhCO2H), and two showed modest HER activity (AcOH and PhOH), as shown in Fig. S36–S66.† However, using a proton source for the CO2RR (electrochemical reduction of carbon dioxide), the catalytic current (icat) was diminished compared with the icat obtained without a proton source, as seen in Fig. 6.
Since the CO2 reduction proceeds under aprotic conditions (vide supra), one of the possible outcomes is the oxalate formation11,13 (the other one, disproportionation to CO and carbonate was not observed) by a two electron process:30
The following experiments were done to shed some light on the catalytic cycle. From the previously discussed experiments, it seems that the active species was a Cu(0) catalyst; thus, complex (2) reacted with sodium-amalgam34 to reduce Cu(I); the Cu(0) intermediate is a paramagnetic intermediate subjected to an EPR spectroscopic study (I = 3/2), where the corresponding spectra are shown in Fig. 8. The g = 2.0012 agrees with closely related Cu(0) compounds35 and confirms a one-electron reduction process; when CO2 is added to the EPR tube, the signal corresponding to copper (0) is not observed due to the oxidation to the Cu(I) complex (red line, Fig. 8).
Regarding the nitrate moiety, it displays the κ1-O for (2) and (3) coordination and κ2-O for compound (5); it is known that nitrate is a weak ligand to Cu(I) which is confirmed by the Cu(I)–ONO2 bond distances obtained from the XRD data of the compounds reported here (Table 6), and these bond distances are longer compared to other complexes reported elsewhere,36 for instance, a Cu–O bond distance on copper nitrate37,38 (1.973 Å, average). These data support that the nitrate group generated in the current report may be prone to dissociate allowing a CO2 molecule to be coordinated and reduced, allowing for an inner sphere mechanism.
Compound | Cu(I)–O (Å) |
---|---|
(2) | 2.1416 |
(3) | 2.09 |
(5) | 2.1498, 2.2580 |
During the CPC experiment described before, the released gases were bubbled in a suspension of 20 mg of Wilkinson catalyst in 5 mL of ethanol.39 The color of the suspension changed from brick-red to a dark yellow precipitate. The IR (ATR) spectra of the solid (Fig. 9) showed a signal at 1655 cm−1 assigned to an N–O bent fragment,40 and the elemental analysis agrees for [Rh(PPh3)3Cl(NO)] (C54H45ONP3ClRh): %C, 67.89, %H, 4.75, %N, 1.46. Found: %C, 67.78, %H, 5.23, %N, 1.47. Fig. S67.†
To investigate if triphenylphosphine has any effect on the system at hand, a CV was done by adding 0.5 equivalents of PPh3; the test showed that PPh3 slightly inhibits the cathodic current of the process (Fig. 10).
To verify the CO2˙− formation, a CPC experiment was done under the same conditions described before, but 10 eq. of TMEDA were added to the matrix cell; at the end of the test, a GC–MS of the solution was obtained, and product (d), as shown in Scheme 3, was detected (Fig. S68†).
It is well known that aliphatic amines generate radical cations.41,42 An initial electron transfer occurs to obtain species (a), which is in equilibrium with (b), then (a) loses a proton and forms the imidium salt (c); finally (b) or (c) can react with the CO2˙− formed during the experiment at the cathode to yield (d), and the protonated species (d) may trap back H+ from the media. The addition of CO2˙− to double bonds has been previously reported,43 and such additions are common in free radical chemistry, Fig. S68.†
Considering the findings of the current report and a recent publication,44 the following simplified catalytic cycle is proposed as shown in Fig. 11.
The original Cu(I) complex is reduced to form a Cu(0) intermediate, then a reaction with CO2 occurs to produce a CO2 adduct, allowing for one-electron reduction to form the CO2˙− radical that dimerizes to yield oxalate.
In the current report, a typical S shape wave was not observed; therefore, an analysis of the onset of the catalytic current can be used by applying the FOWA (foot-of-the-wave-analysis) standard to determine kapp (apparent rate constant) using the CV traces obtained. Thus, the appropriate process for this case begins with the obtention of the value of i/ip:
(1) |
i P: catalyst maximum cathodic current without substrate: 22.1053 microA. (Fig. 4).
i: current obtained from CV (Fig. 5).
n: number of transferred electrons from the electrode to the catalyst: 1.
R: ideal gases constant, 8.314 J K−1 mol−1.
T: working temperature (K), 298.15 K.
F: Faraday constant, 96485 C mol−1.
κ app: as described (vide infra).
ν: working sweep speed: 0.1 V s−1.
n′: catalyst mole number for interchange and product formation: 1.
Then, a graph with coordinates y = i/ip and x = eqn (1) is obtained (see Fig. 13); from the slope (m) the value of kapp is obtained according to eqn (2):
(2) |
To calculate the TOF value for complex (2), the following equation was used.49
(3) |
R = 8.314 J K−1 mol−1.
T = 298.15 K.
F = 96485 C mol−1.
η = 0.39 V.
The values obtained for the CO2RR cathodic process are from Fig. 5.
From Fig. 13, a slope value of 19.984 is obtained, then using that in eqn (2), kapp = 309.8 s−1. Finally, from eqn (3), the TOF value for complex (2) is 154.9 s−1. For the other Cu(I) complexes, (3), (4), and (5), since all of them exhibited the same behavior and the same E°cat potential value (−3.0 V vs. Fc, average) giving oxalate as product, the study was limited to complex (2) due to this compound had the higher icat value; see the ESI.†
Fig. 6 Inhibition of the CO2RR process in the presence of proton sources (red line cathodic peak, 274.21 μA, purple line cathodic peak with PTSA, 136.55 μA). |
Fig. 7 IR (ATR) of calcium oxalate obtained from the CPC experiment of complex (2).31 |
Fig. 9 (a) IR (ATR) of the Wilkinson catalyst. (b) IR-ATR after CPC experiment, NO bent υ = 1655 cm−1.40 |
Fig. 11 Proposed catalytic cycle for the CO2RR to oxalate.45 |
Fig. 12 CPC for the CO2 reduction using (2). Electron number interchanged (α = 2); nprod (3.435 × 10−4), F (Faraday constant). Total charge (QT, 149.9 C). |
Fig. 13 Graph where y = i/ip and x = eqn (1). |
The following describes a particular procedure, but it was observed as a general methodology for all used substrates in terms of the use of the same molar ratios between the substrate and dippe, dppe, depe and 1,1′-bis(di-isopropylphosphine)-ferrocene.
Compounds (3) and (4) were synthesized as described with depe and dppe, respectively. For (5), the dipf was added as a yellow suspension.
1H, δ = 1.20 ppm (s, –CH3), δ = 1.91 ppm (s, –CH2-bridge), δ = 2.15 ppm (s, –CH–isopropyl), δ = 7.37 ppm (m, Ho y Hm –PPh3), δ = 7.54 ppm (t, Hp –PPh3). IR (ATR-neat): 3055 cm−1, 2954 cm−1, 2868 cm−1 (w, CC –PPh3 ring), 1399 cm−1 (s, N–O sym), 1299 cm−1 (s, N–O asymm.), 700 cm−1 (s, –PPh3 ring). Anal. Calcd. for: C32H47O3NP3Cu: %C, 59.11, %H, 7.3, %N, 2.15. Found: %C, 56.9, %H, 7.5, %N: 2.3. Melting point: 202 °C (d).
Compound (3). [(depe)Cu(PPh3)NO3]. NMR (600 MHz): 31P{1H} (THF-d8): δ = −4.5 (d, P–diphosphine), δ = 4.4 ppm (m, –PPh3). 1H: δ = 0.97 ppm (m, –CH3), δ = 1.63 ppm (m, –CH2), δ = 2.54 ppm (m, –CH2 bridge), δ = 7.39 ppm (m, m-PPh3), δ = 7.44 ppm (m, o-PPh3), δ = 7.54 ppm (ws p-PPh3), 13C{1H}, δ = 9.28 ppm (s, –CH3), δ = 17.87 ppm (s, –CH2–), δ = 23.26 ppm (s, –CH2-bridge), δ = 130.03 ppm (m, m-PPh3), δ = 131.22 ppm (ws, p-PPh3), δ = 135.36 ppm (d, o-PPh3, J = 15.09 Hz). IR (ATR-neat): 3053–2875 cm−1 (–CH3 and –CH2–, str.), 1434 cm−1 (N–O asym), 1337 cm−1 (C–H felx.), 1297 cm−1 (N–O sym), 1090–1028 cm−1 (C–H arom.), 694 cm−1 (C–H arom.). Anal. Calcd. for C28H39O3NP3Cu: %C, 56.7, %H, 6.62, %N, 2.3. Found: %C, 55.5, %H, 6.7, %N, 2.41. Melting point: 257 °C (d).
Compound (4). NMR (600 MHz): 31P{1H} (THF-d8): δ = −6.22 ppm (s, P-bridge), δ = 2.91 ppm (m, –PPh3). 1H: δ = 2.42 ppm (m, –CH-bridge), δ = 7.22–751 ppm (m, aromatic phosphines). IR (ATR-neat): 3052 cm−1 (C–H str. arom.), 1434 cm−1 (N–O asym), 1274 cm−1 (N–O sym), 1096–1022 cm−1 (C–H arom.), 692 cm−1 (C–H arom.). Anal. Calcd. for C44H43O3NP3Cu: %C, 66.8, %H, 5.4, %N, 1.8. Found: %C, 66.8, %H, 5.1, %N, 2.6. Melting point: 256 °C (d).
Compound (5). NMR (600 MHz): 31P{1H} (THF-d8), δ = 3.13 ppm (s). 1H, δ = 1.30 ppm (m, –CH3), δ = 2.29 ppm (m, –CH–), δ = 4.44 (s, –CH, Fc), δ = 4.53 ppm (s, –CH, Fc). 13C{1H}, δ = 20.59 ppm (m, –CH3), δ = 21.44 ppm (–CH–), δ = 72.48 ppm (–CH–, Fc), δ = 75.19 ppm (–CH–, Fc). IR (ATR-neat): 2958–2866 cm−1 (–C–H alkyl), 1431 cm−1 (N–O asym), 1286 cm−1 (N–O sym), 1024 cm−1 (Fc), 820 cm−1 (Fc). Anal. Calcd. for C22H36O3NP2FeCu: %C,48.59, %H, 6.67, %N, 2.58. Found: %C, 50.22, %H, 6.81, %N, 3.11. Melting Point: 183 °C (d).
Footnote |
† Electronic supplementary information (ESI) available: X-ray determinations have been deposited at CCDC with the following numbers: 2314253–2314255. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4cy00759j |
This journal is © The Royal Society of Chemistry 2024 |