To understand the local ion transport in a BPMEA, a 1D isothermal continuum model was developed based on previous work by Weng et al.50 and Lees et al.51 The model domain includes an ionomer-immersed porous cathode catalyst layer (CL) and CEL of the BPM. The catalyst layer consists of nickel-nitrogen-carbon catalyst (NiNC-IMI) mixed with either cation- (Nafion, a sulfonated fluoropolymer52) or anion-exchange ionomer (Sustainion, an imidazolium functionalized styrene polymer31). We chose a steady-state model to capture the ionic behaviour at the initial stage of BPMEA operation. The model was fit to the experimental CO Faradaic efficiency data collected from NiNC-IMI catalyst layers with Sustainion ionomers by adjusting the electrochemical parameters such as the transfer coefficients and exchange current densities for CO2 reduction and the competitive HER (see Fig. 2a). To validate the model, the same kinetic parameters were then used to predict the CO Faradaic efficiency data collected with the Nafion ionomer. Details of the models (e.g., equations, boundary conditions, and parameters) are described in the Supplementary Information.
Since previous work has shown the role of anolyte ion concentration and crossover on cathodic performance and stability in a BPMEA CO2 electrolyzer,40 we used the model to more deeply investigate two strategies to control the ionic transport within the cathode. Specifically, (i) the use of ionomer in the CL for selective ionic transport and (ii) promoting an increased proton transference number in the CEL. Including ionomer in the catalyst layer is a common effective approach to modulate local ionic transport. Further, maintaining a high proton flux from the CEL to the cathode is also a prerequisite for a stable and efficient BPMEA system. Near-unity water dissociation efficiencies are desired as they limit ionic interactions between the anode and cathode environments that can impact anolyte pH and salt precipitation at the cathode. Therefore, these two strategies are perceived as practical approaches to managing local ion transport and thus reaction microenvironment in the electrolysis cell.
We first examined the role of ionomer choice for the CL in determining local ionic concentrations in the catalyst layer. As presented in Fig. 2b and c, incorporating anion-selective Sustainion ionomer in the CLs (Sus-CLs) leads to a counterintuitively lower pH and CO32− concentration than the Nafion ionomer in the CL (Naf-CLs). Notably, the Sus-CLs case provided a substantially higher (7 folds) HCO3− concentration near the CEL|CL interface, and over the entire catalyst layer an average HCO3− concentration that is more than twice that of the Naf-CLs case. This discernible divergence is a result of the different fixed charges of the two ionomers. Anion exchange ionomers promote the transport of generated (bi)carbonates near the gas-liquid interface towards the BPM, while the Nafion rejects this transport and promotes (bi)carbonate accumulation near the generation point. The positive fixed charge of the Sus-CLs case then provides ample HCO3− available for acidification and CO2 regeneration near the CEL|CL interface.
In addition to the anion transport, it is important to assess the cation transport (H+ and K+). As further suggested from the calculated K+ profiles shown in Fig. S3a, the positively charged quaternary ammonium groups in the ionomer at least partially exclude K+ transport from CEL to CL in the Sus-CLs and thus likely contribute to the observed reduced pH in the CLs by lowering the required amount of OH− ions needed to balance the positive charge. By contrast, the concentration profiles predicted for Naf-CL indicate a more uneven ionic distribution as compared to Sus-CL. The negatively charged sulfonic groups in Nafion lead to an excessive amount of K+ in the Naf-CL (see Fig. S3a), which then fosters a high content of OH− (or high pH shown in Fig. 2b) and CO32− to maintain the charge neutrality.
Overall, both ionomer cases show the ability for CO2 converted to (bi)carbonates to be regenerated into CO2 by the proton flux from the CEL of the BPM. Interestingly, the concentration of HCO3− at Naf-CL shown in Fig. 2c decreases from around 1.0 M at CL|GDL interface down to 0.047 M at the CEL|CL interface. Such a steep decrease in HCO3− concentration is a result of increased local pH inside Naf-CL. Due to the high local pH inside Naf-CL (Fig. 2b), however, the regenerated CO2 tends to diminish and form back to CO32− within the Naf-CL (Fig. 2c). As such, the uneven ionic distribution in Naf-CL might not be ideal for an efficient electrochemical conversion of CO2 because it could cause a low utilisation efficiency of the reverted CO2 for electrochemical conversion. Similarly, the local CO2 concentration near the Sus-CL|CEL interface is slightly higher than Naf-CL|CEL interface (Fig. 2d) due to a higher concentration of HCO3− which promotes CO2 regeneration. However, the CO2 concentration throughout the bulk of the CL is similar for the Sus-CL and Naf-CL because of the constant excess CO2 supply provided at the CL|GDL interface.
Next, we use the Sus-CL model to investigate the impact of proton transference number on the local ionic transport across the CLs. The ionic conduction across the CEL|CL interface relies on the transport of H+, K+, and (bi)carbonate ions, with H+ being the primary charge carrier. The H+ transference number quantifies the fraction of the current crossing the CEL in the absence of concentration gradients, which is a result of proton transport from the BPM. The total ionic current is the sum of the H+ transport, K+ crossover from the anode and anion crossover from the cathode. By sweeping the H+ transference number from 0.75 to 0.95, as presented in Fig. 2e-g and Fig. S3b, we observed a decrease in pH, K+ and CO32− concentrations but an increase in HCO3− concentrations in the CLs. This trend is as expected because an increase in H+ transference number indicates a suppressed current associated with K+ cations and promotes carbonate acidification within the CLs. Additionally, as shown in Fig. 2g, increasing the H+ transference number does not notably increase the local CO2 concentration but reduces the CO2 concentration within the CEL due to the minimised (bi)carbonate crossover from CL to CEL. The downside of an increased H+ flux across CEL|CL interface could be HER out competing CO2 reduction through direct H+ reduction and limited cations available to activate CO2 electrochemical reduction.
The role of ionomer in CL on BPMEA performance
To validate the predictions of the model, we compared experimentally NiNC-IMI CLs prepared with Sustainion and Nafion ionomers. A comprehensive study of the NiNC-IMI catalyst has been published elsewhere53,54, so this study will focus on the role of the ionic transport on the BPMEA performance by using the NiNC-IMI as a model CL. The products from the BPMEA cells are primarily CO and H2 with minor formate or formic acid (see nuclear magnetic resonance results in Fig. S4). Because the BPM is effective in minimizing formate crossover to the anolyte55, it is challenging to calculate the formate FEs accurately by analyzing anolyte compositions. Hence, this work focuses mainly on the CO and H2 products from the electrolysis.
Figure 3b-3d show that the ionomer in the CL defines cathode activity and selectivity in the BPMEA configuration. 15 wt% Sustainion ionomer in the CL can profoundly improve the catalyst CO FE and partial current densities compared to the 15 wt% Nafion. Due to the different equivalent weight of the ionomers (1100 g/mol for Nafion; ~ 225 g/mol for Sustainion31), 15 wt% Sustainion contains approximately five times the concentration of ionic groups as compared to the 15 wt% Nafion. As shown in Fig. 3b, the CO FE is 69.7 ± 0.5% at 50 mA cm− 2 and 45.4 ± 0.5% at 200 mA cm− 2 for CLs with 15% Sustainion ionomer, much higher than the CLs based with Nafion (53.5 ± 4.2% at 50 mA cm− 2 and 21.4 ± 1.0% at 200 mA cm− 2). The Naf-CLs reach a CO limiting current density of ~ 41 mA cm− 2, more than two-fold lower than the Sus-CLs (> 105 mA cm− 2).
The notable enhancement in CO2-to-CO upon FEs and partial current densities is related to the local mass transport within the CLs that either improves the local concentration of CO2 or limits the local concentration of protons. Fig. S5 demonstrates no significant suppression of HER by the Sustainion ionomer, implying a similar local proton availability and/or catalytic HER activity in both cases. Hence, we postulate that the observed improvement likely originates from the enhanced bicarbonate local transport within Sustainion-based CLs, as predicted from our model in Fig. 2c. The increase in local bicarbonate concentration in the CL provides an even CO2 distribution in the CL and thus promote CO2 electroreduction.
In the BPMEA reported here, there are two sources of CO2: CO2 that dissolves from the gas feed and CO2 that is regenerated from bicarbonate via acid-base chemistry. The modelled CO2 concentration profiles (Fig. 2c and f) show that the potent proton flux from CEL yields a peak in CO2 concentration at the CEL|CL interface. The increased CO2 concentration results in an increase in the predicted CO partial current density near the CEL|CL interface as shown in Fig. S6. This result is consistent with our previous bicarbonate electrolyser model that shows high CO formation rates at the CEL|CL interface.51 In this BPMEA model, a high rate of CO formation is also observed at the CL|GDL interface. This difference is attributed to the additional CO2 source from the gas phase, which is absent in the bicarbonate electrolyser. These modelling results confirm that both regenerated and dissolved CO2 contribute to CO production in a BPMEA.
Since bicarbonate anions exhibit more facile transport in the anion-exchange Sustainion ionomer than the cation-exchange Nafion ionomer18 (see Fig. 2b), we predict that the Sustainion ionomer should allow a facile transport of bicarbonate anions for CO2 regeneration, providing ample local CO2 for electrochemical conversion across the CL structure (see Fig. 3a). By contrast, in the Naf-CL, the regenerated CO2 tends to be spatially localized, causing a limited proportion of catalyst surface in the CLs to be accessible by the regenerated CO2. Therefore, the even distribution of bicarbonate enables the Sustainion-based CL with a larger catalyst surface area accessible by the regenerated CO2 than a Nafion-based CL56. Thus, while both systems enable reasonable CO2 utilization efficiencies, Sus-CL case exhibits both higher CO2 utilization efficiencies, CO FE, and partial current densities at comparable cell voltages (Fig. 3d).
Figure 3e shows that the CO2-to-CO utilisation efficiency increases with current density. Here, the CO2-to-CO utilisation efficiency is defined as the ratio of the CO2 reduced to CO versus the total CO2 consumed in the cell (i.e., due to crossover and reaction). The observed CO2-to-CO utilisation trend as a function of current density could result from an increase in the H+ transference number across the CEL of the BPM, because water dissociation dictates the overall ionic current while the co- and counter-ion crossover is mass-transport limited at increased current densities.57 When one mole of CO is produced, for example, there will be two moles of electrons consumed, two moles of H+/K+ transported from the CEL, and two moles of OH− co-produced. The OH− can then be converted into one-mole carbonate or two moles of bicarbonate, which require two moles of protons to remove the (bi)carbonate species. Therefore, the H+ transference number determines the availability of the H+ to revert the generated (bi)carbonate species back to CO2. As the H+ transference number is always below 1, it means that some CO2 converted to (bi)carbonates are never recovered and will either precipitate or crossover to the anode. Consequently, high H+ transference numbers across the BPM are beneficial to lowering the local pH and regenerating CO2 from (bi)carbonate, hence, maximizing CO2-to-CO utilisation efficiency.
A further increase of Sustainion ionomer loading to 30 wt% in the CL lowered the FEs and limiting current density for CO production. (Fig. 3b-3d) Such decline in performance is likely attributable to the blockage of the CL pores and reactive sites by the ionomer itself, resulting in an increased diffusion length for CO2 gas from gas channels to the catalyst surface or loss of active surface in the CLs, respectively. However, the CO2-to-CO utilisation efficiency is not significantly impacted by the increment of the ionomer loading, implying that in-situ formed (bi)carbonate ions can be reverted to CO2 easily within the CL matrix.
More importantly, as shown in Fig. 3f, the anolyte shows only a slight decrease in pH after the cell testing, with a total charge of 6300 C passed through the cell during the test. The pH of the anolyte could have been decreased because of both carbon crossover from the cathode to the anode and cation crossover from the anode to the cathode. The minimal change in the anolyte pH is due to the unique BPM function, which supplies protons at the cathode to convert (bi)carbonate anions to CO2 for CO production and minimises ion crossover rates across the cell. This feature allows BPMEA cells to be operated stably using a PGM-free anode based on nickel and with a high CO2-to-product utilisation efficiency, which could not be easily done using monopolar membrane-based electrolysis cells.
The role of the H+ transference number on BPMEA performance
The role of proton transference number is investigated experimentally by examining the effect of K+ cation crossover on CO2 reduction using 0.1 M and 1 M KOH as the anolyte. A concentrated KOH is expected to accelerate the crossover rate of K+ due to the large concentration gradient across the membrane (see Fig. S7a). Because cations are essential in activating CO2 reduction, the results in Fig. S7b show that 1 M KOH anolyte can slightly improve CO FE at current densities > 150 mA cm− 2. This finding is consistent with our previous report highlighting the cations' vital role in enhancing CO FEs over silver electrodes.40 Fig. S7c suggests that the BPMEA cell with a dilute anolyte shows a higher CO2-to-CO utilisation efficiency than the concentrated anolyte. This trend matches our model predictions of the proton transference number in Fig. 2d-f, because a reduced K+ crossover rate achieved by using a dilute anolyte leads to an increase in proton transference number and acidification of (bi)carbonate within the CL. Therefore, a trade-off exists between product selectivity and CO2 utilisation efficiency when choosing the anolyte for the BPMEA cell.
Local H+ enrichment at CEL/CL interface
In addition to the distribution of (bi)carbonate ions within the CL, we studied the local reaction environment at the CEL|CL interface by introducing a hydrophilic porous spacer (65 um thick) between the CL and CEL. (Fig. 4a) This configuration is fundamentally different from the previously reported method41 that includes concentrated salts within the spacer. In this study, we applied the spacer pre-soaked with ultrapure water before the cell assembly. The ion conduction within the spacer solely depends on the ionic fluxes of protons, (bi)carbonate anions, and K+ cations that come from the fixed charge in the CEL and anolyte.
As shown in Fig. 4b, the spacer significantly suppresses the HER down to below 15% and boosts the CO FEs up to 91% at 50 mA cm− 2 and 88% at 100 mA cm− 2, which is almost comparable to MEA cells based on monopolar membranes. The comparison of the CO partial current densities vs. cell potentials, as shown in Fig. S8, suggests cells with and without spacer achieve similar CO partial current densities under similar cell voltages. Therefore, the observed CO FE enhancement is mainly attributable to the suppression of HER in the CLs. This finding also indicates that the majority of HER in the absence of spacer occurs at the CL|CEL interfaces due to the direct contact of the CL with excess protons; the spacer increases the retention time for protons to reach the catalyst surface.
The drawback to including the spacer in the BPMEA cell is the large ohmic loss due to the slow ionic conduction across the spacer, as verified by the electrochemical impedance spectroscopic analyses in Fig. S9. The absence of abundant water at the cathode side of the BPMEA accelerates the dehydration of the spacer, which causes further reduction of the ionic conductivity of the spacer and eventually rapid cell voltage overshoot (See Fig. S10). Nonetheless, the results highlight the importance of the cell configuration in determining the CO2-to-CO selectivity of the BPMEA cells.
Long-term stability and ion crossover
Finally, we evaluated the long-term stability of the PGM-free BPMEA cell, particularly as it relates to the anolyte pH and the stable use of nickel anodes. Figure 5a demonstrates that the CO Faradaic efficiency rapidly dropped from 73–64% within the first 1.5 h, then decreased linearly to 30% after 150 h at a degradation rate of 0.36% per hour. Meanwhile, the HER increased along with the loss in CO FE. Such discernible selectivity loss likely originates from an increase in contact area between catalyst and H+ and possible deactivation of the catalyst in the acidic environment. Both could contribute to the rise in the rate of HER and suppression of the rate of CO2 reduction.
The cell voltage also increased at a degradation rate of 0.086% per hour. By comparing the ohmic resistances obtained from electrochemical impedance (Fig. S11), we found that the degradation of the cell potential should be related to the increase of the polarization resistance over either the cathode or anode rather than the deterioration of the BPM.
On the other hand, as shown in Fig. 5b, the CO2-to-CO utilisation efficiency increased rapidly in the first 1.5 h and was subsequently sustained at above 80–90% across the measurement. The rapid rise in the CO2 utilisation within the first 1.5 h implies that the rapid drop in CO FE is mainly a result of an increase of proton local availability for HER in the CL. The subsequent linear degradation is likely associated with the instability of the catalyst layer in a strongly acidic environment.
Without replenishing the anolyte across the entire test, we observed no significant change in the pH value of the anolyte after the test, as shown in Fig. 5c. The titration results shown in Fig. 5d and Fig. S12 revealed that about 0.006 mol K+ cations have migrated to the cathode side, which is equivalent to 0.21% of the ionic current during the conditioning. Due to the anolyte's high alkalinity, most of the CO2 crossover should be converted to carbonate ions in the anolyte. Figure 5d shows that the anolyte after the test is composed of 0.0735 M OH− and 0.011 M CO32−, implying that the (bi)carbonate crossover contributed to < 0.774% of the ionic current. The extremely low K+ and CO2 crossover rates were essential in sustaining the Ni-based anode stability without replenishing the anolyte during the 150-h test. We posit that a proton transference number of ~ 99% was maintained over the course of the experiment. While the 0.21% K+ crossover is likely needed to maintain CO2 reduction at the cathode, the 0.774% carbonate crossover is all that needs to be avoided to maintain anolyte pH indefinitely.
After 150 h conditioning within the BPMEA cell, as revealed by the XRD results in Fig. S13, the Ni anode remained as metallic Ni, while additional minor phases related to potassium nickelates58 formed due to the long-term oxidising treatment. Importantly, no nickel carbonate phase was detected in the conditioned anode, implying that the BPMEA configuration stabilised the Ni-based anode. Ex-situ surface analyses further confirm that both fresh and conditioned anodes are predominantly covered with Ni(OH)2 species at the surfaces (Fig. S14), and that their structures shown in Fig. S15 and S16 remain intact after 150 h conditioning. The retained structural and chemical stability of the Ni anode is a main result of the alkaline local environment achieved by the BPM, which supplies hydroxide ions and suppresses K+ and CO2 crossover.
We also noticed that the inlet pressure of the cell rose after an interval of 12–30 h. The rise of the inlet pressure is mainly a result of the build-up of precipitated (bi)carbonate salts at the entrance of the gas channel in the cell, though there was no discernible salt precipitation elsewhere in the gas channels or the back of the cathode (see Fig. S17). The salt precipitation occurs typically when the concentration of the local (bi)carbonate exceeds the solubility (e.g., 8.03 M for K2CO3 and 3.62 M for KHCO3)26,28 at the cathode structure and the crystal can grow continuously to reach the gas channels, similar to efflorescence process, due to its hygroscopic and porous nature. The inlet pressure rise (due to the salt precipitation) elucidates that there is a continuous supply of K+ from the cation crossover of about 0.21% to the cathode for the salt precipitation, which cannot be easily mitigated, especially at a high current density.
To circumvent this issue, we applied ~ 1 mL pure water pulse to wash off the precipitated salt from the gas channel during operation when an increase of inlet pressure was observed. Unlike the reported hourly water flush to remove the salts at the cathode for the anion-exchange membrane electrode assembly cell14,29, the BPMEA cell requires a much less frequent water pulse, thanks to its high CO2 utilisation and controlled K+ concentrations at the cathode. The water pulse operation also showed negligible impact on the cell potential and CO Faradaic efficiency, as shown in Fig. 5a.