Surani Bin
Dolmanan‡
a,
Annette
Böhme‡
bc,
Ziting
Fan‡
def,
Alex J.
King
gh,
Aidan Q.
Fenwick
bi,
Albertus Denny
Handoko
a,
Wan Ru
Leow
j,
Adam Z.
Weber
gh,
Xinbin
Ma
de,
Edwin
Khoo
k,
Harry A.
Atwater
*bc and
Yanwei
Lum
*af
aInstitute of Materials Research and Engineering, Agency for Science, Technology and Research (A*STAR), Innovis, Singapore. E-mail: lumyw@nus.edu.sg
bLiquid Sunlight Alliance, California Institute of Technology, Pasadena, California, USA. E-mail: haa@caltech.edu
cDepartment of Applied Physics and Material Science, California Institute of Technology, Pasadena, California, USA
dJoint School of National University of Singapore and Tianjin University, International Campus of Tianjin University, Binhai New City, Fuzhou, China
eKey Laboratory for Green Chemical Technology of Ministry of Education, Collaborative Innovation Center of Chemical Science and Engineering, School of Chemical Engineering and Technology, Tianjin University, Tianjin, China
fDepartment of Chemical and Biomolecular Engineering, National University of Singapore, Singapore, Singapore
gDepartment of Chemical and Biomolecular Engineering, University of California Berkeley, Berkeley, California, USA
hLiquid Sunlight Alliance, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
iDepartment of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California, USA
jInstitute of Sustainability for Chemicals, Energy and Environment (ISCE2), Agency for Science, Technology and Research (A*STAR), Jurong Island, Singapore
kInstitute for Infocomm Research, Agency for Science, Technology, and Research (A*STAR), Connexis, Singapore
First published on 30th May 2023
Gas diffusion layers (GDL) have become a critical component in electrochemical CO2 reduction (CO2R) systems because they can enable high current densities needed for industrially relevant productivity. Besides this function, it is often assumed that the choice of catalyst and electrolyte play much more important roles than the GDL in influencing the observed product selectivity. Here, we show that tuning of the GDL pore size can be used to control the local microenvironment of the catalyst and hence, effect significant changes in catalytic outcomes. This concept is demonstrated using sputtered Ag films on hydrophobic PTFE substrates with 6 different pore sizes. Although Ag is known to be a predominantly CO generating catalyst, we find that smaller pore sizes favor the generation of formate up to a faradaic efficiency of 43%. Combined experimental and simulation results show that this is due to the influence of the pore size on CO2 mass transport, which alters the local pH at the electrode, resulting in reaction pathway switching between CO and formate. Our results highlight the importance of the local microenvironment as an experimental knob that can be rationally tuned for controlling product selectivity: a key consideration in the design of CO2R systems.
Furthermore, the reactor system has also been shown to be important. In the majority of early reports in the field, CO2R was conducted in H-type cells where CO2 is introduced into the system via continuous bubbling into the electrolyte.26,27 However, the low solubility of CO2 (33 mM) typically results in limiting current densities of only several tens of mA cm−2 due to significant mass transport limitations. To raise current densities towards industrially relevant productivity, catalyst particles are deposited onto gas diffusion layers (GDL), allowing CO2 mass transport limitations to be overcome28–34 and enabling current densities of >100 mA cm−2. This is due to the hydrophobic and porous nature of the GDL, resulting in the creation of thin layers of electrolyte over the catalyst particles. These thin layers of electrolyte have significantly lower CO2 transport diffusion lengths, thus facilitating rapid supply of reactants to the catalyst surface. Beyond this role, it is often assumed that the choice of catalyst material and electrolyte play more dominant roles compared to the GDL in controlling the observed product selectivity.
Although the mass transport of CO2 through the GDL should in principle be rapid, it is known that its effective diffusion coefficient is related to the porosity and average pore radius of the porous medium through the Bruggeman relationship.35 We therefore reasoned that tuning these parameters could be used to influence the mass transport of CO2, which directly impacts the catalyst microenvironment (local pH and CO2 reactant supply). This is because CO2 molecules can directly react with and hence neutralize electrochemically generated OH− to form bicarbonate and carbonate anions.27,36 The altered microenvironment could in turn result in a significant change in catalytic outcomes: an additional experimental knob to control CO2R selectivity beyond catalyst design and choice of electrolyte.
In this work, we demonstrate this concept using sputtered Ag films onto hydrophobic PTFE substrates with 6 different pore sizes as the GDL37,38 (Scheme 1). Even though Ag is well known to predominantly produce CO,39–43 we find that smaller PTFE pore sizes favor formate production up to a FE of 43%. Combined experimental and simulation results show that a decrease in GDL pore size slows down CO2 mass transport, leading to a higher local pH and hence reaction pathway switching from CO to formate. This pH trend was confirmed using a confocal microscopy setup44,45 equipped with a custom-built electrochemical cell and a pH sensitive fluorescent dye in the electrolyte. Our results highlight the importance of the properties of the GDL, which can significantly impact the catalyst local microenvironment and should be an important consideration for the design of CO2R systems.
Fig. 1 Multiphysics simulation results of varying GDL porosity on the (a) local pH and (b) CO2 concentration at various cathodic current densities of 100, 200, 300, 400 and 500 mA cm−2. A lower porosity is observed to result in a higher local pH and lower CO2 concentration. Detailed results can be found in the ESI (Fig. S1–S6†). |
Encouraged by these results, we began by sputtering 325 nm thick Ag films onto hydrophobic PTFE37,38 substrates with different pore sizes of 0.02, 0.1, 0.22, 0.45, 1.0 and 3.0 μm for use as gas diffusion electrodes (Fig. S7†). Each electrode will hence be termed as Ag/PTFE(X), where X is the pore size. Top-down scanning electron microscopy (SEM) images show the structure of these electrodes, with a web-like morphology of interconnected PTFE fibers coated conformally with Ag (Fig. 2). These SEM images reveal the 3D network of macro-scale pores that are inherently formed between the fibers, serving as pathways for reactant and product transport (Fig. S8–S10†). As would be expected, the PTFE substrates with larger pore sizes appear visibly more open and less dense.
Fig. 2 SEM images of hydrophobic PTFE substrates of various pore sizes coated with 325 nm of Ag using sputter deposition. The pore sizes are: (a) 0.02 μm, (b) 0.1 μm, (c) 0.22 μm, (d) 0.45 μm, (e) 1.0 μm and (f) 3.0 μm. Digital photographs (Fig. S7†) and more SEM images (Fig. S8–S10†) of the samples can be found in the ESI.† |
X-ray diffraction (XRD) characterization of the electrodes was performed (Fig. 3a), with Ag (111) observed as the dominant crystal facet and with no obvious differences between each of the Ag/PTFE with various GDL pore sizes. We also carried out cyclic voltammetry in a potential range where only non-faradaic processes occur to determine the double layer capacitance of each Ag/PTFE electrode (see methods section in ESI†). This gives an indication of the electrochemically active surface area (ECSA) since this value is directly proportional to the double layer capacitance.46 The results (Fig. S11 and S12†) show that despite pore size differences, the double layer capacitance and hence ECSA remains approximately within the same order of magnitude.
We then designed experiments to obtain a qualitative measure of the CO2 mass transport for the different pore size Ag/PTFE electrodes. Each electrode was assembled into a gas diffusion flow cell system (Fig. S13†), with a similar design to what was previously reported in the literature.30,37 15 ml of 1 M KOH was used as the electrolyte, which was continuously recirculated between the cathode chamber and an external centrifuge tube reservoir using a peristaltic pump. CO2 was flowed at a rate of 20 sccm, through a gas chamber in contact with the backside of Ag/PTFE. Without applying any current, we monitored the bulk pH of the electrolyte over a 120 min period by placing a pH probe into the external centrifuge tube reservoir. The results in Fig. 3b show that the bulk pH decreases significantly with time, as a result of the CO2 gas continuously diffusing from the backside of the Ag/PTFE and reacting with hydroxide in the electrolyte to form carbonate.36 We also observe that the bulk pH decreases more rapidly with increasing PTFE pore size. Importantly, this allows us to experimentally confirm that larger pore sizes do indeed facilitate faster CO2 mass transport.
Next, we sought to assess the influence of PTFE pore size on the product selectivity of the Ag/PTFE catalysts. Using the same flow cell system, we evaluated each Ag/PTFE under cathodic current densities of 100, 200 and 300 mA cm−2 in 1 M KHCO3 electrolyte and the FE data are shown in Fig. 4a–c. Based on the results, we observe that the formate FE appears to increase with decreasing pore size, from 24% for Ag/PTFE(3.0) up to a value of 43% for Ag/PTFE(0.02) at 200 mA cm−2. For better visualization, the formate FE is also presented as a contour plot (Fig. 4d), where the general trend of higher formate FE with smaller pore sizes is observed to hold true for all tested current densities. Also, the hydrogen FE tends to increase with larger pore size. These combined effects result in the CO FE initially increasing with pore size and then decreasing again, with a peak value of around 80% at 100 mA cm−2 for Ag/PTFE(1.0).
Fig. 4 Electrochemical CO2 reduction FE results with 1 M KHCO3 as the electrolyte. (a), (b) and (c) show the product FE data for Ag/PTFE as a function of GDL pore size under cathodic current densities of 100, 200 and 300 mA cm−2 respectively. (d) Is the corresponding color contour map of the HCOO− FE data for Ag/PTFE as a function of current density and GDL pore size. More data available in the ESI (Fig. S16†). |
Based on the bulk pH monitoring and simulation results, we hypothesized that this could be due to reduced CO2 mass transport at the smaller pore sizes, resulting in a higher local pH and, thus, switching selectivity towards formate. This selectivity switching was previously observed by Seifitokaldani et al., where CO2R was performed with Ag catalysts in KOH electrolyte47 with concentrations ranging from 0.1 M to 11 M. It was found that formate was produced with almost 60% FE in 11 M KOH, compared to only about 4% in 0.1 M KOH. Using DFT simulations, they concluded that this was due to the activation energy barrier for formate becoming lower than that compared to CO, in the absence of hydronium ions.
Hence, we employed a suite of experiments to further understand these initial observations and verify our working hypothesis. Firstly, we tested the Ag/PTFE catalysts in 2 M KHCO3, which has a stronger pH buffering ability as compared to 1 M KHCO3.48 In this case, we did not observe any significant differences in the formate FE as a function of GDL pore size (Fig. S14a and b†) at cathodic current densities of 100 and 200 mA cm−2. This suggests that the stronger buffer results in a similar local pH value for each of these cases, leading to a similar formate FE of around 14%. However, at the higher current density of 300 mA cm−2, the trend of higher formate FE with smaller pore size appears again (Fig. S14c†), with a FE of 19% for Ag/PTFE(3.0) as compared to a FE of 29% for Ag/PTFE(0.02). This results from the expected higher local pH rise with a larger current density and is therefore consistent with the notion that pH effects are indeed influencing the observed FE to formate.
To further investigate the effect of buffering, similar CO2R experiments were carried out with additional buffer conditions of 0.5 M and 1.5 M KHCO3 at 300 mA cm−2 for each pore size condition. The results for Ag/PTFE(0.22) are represented in Fig. 5a, where a trend of higher formate FE with lower buffer concentration is observed. This is because lower buffer concentrations result in a higher local pH,48 which then promotes the conversion of CO2 to formate. Formate FE for all Ag/PTFE samples under the different buffer conditions are shown as a colour contour map (Fig. 5b), where the trend of higher formate FE with a lower buffer concentration is observed to hold true for all GDL pore sizes.
Fig. 5 (a) HCOO− FE data for Ag/PTFE(0.22) as a function of KHCO3 concentration. (b) Color contour map of the HCOO− FE data for Ag/PTFE as a function of KHCO3 concentration and GDL pore size. More data for (a) and (b) are available in the ESI (Fig. S15–S17†). (c) Product FE data for the case where the CO2 feed stream was diluted with various amounts of N2. Ag/PTFE(0.45) was used as the cathode and 1 M KHCO3 was used as the electrolyte at an applied cathodic current density of 300 mA cm−2. (d) Graph showing the formate to CO ratio as a function of N2 gas dilution, based on the data shown in (c). |
We also carried out CO2R electrolysis experiments where the CO2 feed was diluted with N2. For a lower CO2 partial pressure, we expect the local pH to be higher due to fewer available CO2 to react with electrochemically formed OH−. For these experiments, Ag/PTFE(0.45) was used as the electrode and a constant current density of 300 mA cm−2 was applied. From the results (Fig. 5b and c), we observe that lower CO2 partial pressures do indeed result in a higher formate to CO ratio, consistent with our working hypothesis.
Furthermore, we also conducted in situ measurements using confocal microscopy with a pH sensitive fluorescent dye to provide experimental verification of the local pH trends as a function of GDL pore size and applied current density. Fluorescent confocal laser scanning microscopy enables imaging of the local pH in three spatial dimensions with a resolution of one micrometer under operating conditions.44,45 Such experiments were carried out using a custom-built electrochemical cell (Fig. 6a and S18†), consisting of a gas chamber for CO2 flow and an electrolyte chamber that is stacked above it. The electrolyte chamber is open at the top, which allows for a water immersion objective to be dipped into the electrolyte, in close proximity to the electrode surface.
Fig. 6 (a) Cross section of the custom-built electrochemical cell with water immersion objective for local pH measurements with APTS. (b) Local pH value averaged from zero to forty micrometers above the electrode surface as a function of current density for different Ag/PTFE GDLs. (c) Representative pH maps as a cross section through the plane perpendicular to the electrode surface for different Ag/PTFE samples at 100 mA cm−2. pH maps for other current densities can be found in Fig. S19.† |
We studied the local pH in the vicinity of the Ag/PTFE electrodes with the ratiometric fluorescent dye 8-aminopyrene-1,3,6-trisulfonic acid trisodium salt (APTS), which is dissolved in the electrolyte. APTS is a direct sensor of the local hydroxide activity and can be used to deduce the local pH (sensitive to values between 11.2 and 14). A more detailed analysis of the sensing mechanism of APTS can be found elsewhere.45 Measurements for all cases were performed using 1 M KHCO3 electrolyte with 200 μM APTS, which was constantly circulated through the electrochemical cell (more details in the ESI†).
The pH was mapped for a series of current densities between 20 mA cm−2 and 200 mA cm−2, in the plane perpendicular to the electrode surface, starting from a few micrometers below the electrode surface. The dimensions of each pH map are 120 μm in the x direction and 50 μm in the z direction (Fig. 6c and S19†). From the maps, a pH gradient can be clearly observed for all cases, with the pH being higher at points closer to the electrode surface. We averaged the pH in the area between the electrode surface and 40 micrometers above the electrode surface and plotted this as a function of current density (Fig. 6b). As expected, the pH increases as the current density increases since OH− is created as a by-product of CO2 reduction and hydrogen evolution.
Most importantly, for all current densities investigated, there is the clear trend that the local pH decreases with increasing pore size of the GDL. This is consistent with our preceding experimental and simulation results, that a larger GDL pore size can indeed better facilitate CO2 mass transport. This leads to more excess CO2 molecules that are available at the electrode surface to react with electrochemically generated OH−, leading to a lower pH value. These results are therefore strong experimental evidence for our hypothesis that tuning the GDL pore size can indeed directly impact the local pH. This then results in selectivity switching, leading to the observed increased selectivity towards conversion of CO2 to formate.
Finally, we sought to understand if this local microenvironment effect could also affect other catalysts for electrochemical CO2 reduction. Previous literature reports have indicated that an increased local pH can induce a higher selectivity towards multicarbon (C2+) products with Cu based catalysts. Hence to explore this effect, we prepared a series of samples by sputtering 325 nm of Cu onto PTFE substrates of different pore sizes (Fig. S20†). These catalysts were then tested at a constant cathodic current density of 200 mA cm−2 in 1 M KHCO3 electrolyte. The results (Fig. S21†) show that a smaller pore size does indeed lead to an increase in the FE towards C2+ products, as a consequence of the induced higher local pH.
However, the C2+ FE was observed to drop once the pore size becomes too small. This is because the local CO2 availability is expected to diminish at the smallest pore sizes, and these conditions are less favorable towards the formation of C2+ products. Our findings are consistent with the work by Strasser and co-workers, where they observed different “selectivity zones” within their Cu nanoparticle catalyst coatings on a gas diffusion layer.49 Zones closer to the gas diffusion layer experience higher local pH and increased CO2 availability, which enhances the C2+ selectivity. On the other hand, zones further away from the gas diffusion layer experience CO2 depletion, which reduces C2+ product formation. These observations are supported by another report in the literature,50 which showed that a lowered CO2 partial pressure suppresses the C2+ FE.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ta02558f |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2023 |