Yaping
Pan‡
ab,
Han
Liu‡
ab,
Jiawei
Liu
ab,
Linrui
Wen
b,
Kejie
Lao
ab,
Shuirong
Li
d,
Xiaoliang
Fang
d,
Huakun
Wang
*c,
Hua Bing
Tao
*ab and
Nanfeng
Zheng
ab
aState Key Laboratory for Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China. E-mail: hbtao@xmu.edu.cn
bInnovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen 361005, China
cFujian Key Laboratory of Digital Simulations for Coastal Civil Engineering, Xiamen University, Xiamen 361005, China. E-mail: hkwang@xmu.edu.cn
dCollege of Energy, Xiamen University, Xiamen 361005, China
First published on 13th February 2024
A better understanding of current density distribution, i.e. the true electrochemical reaction behavior, in electrochemical energy devices has long been desired. However, it is considered very challenging to comprehensively observe the current density within real PEMWE devices at the micrometer scale as the components are densely assembled by metal plates. Thus, an ex situ method for micron-level current density detection is proposed in this study, where the highly active but unstable characteristic of RuO2 in the oxygen evolution reaction (OER) in acidic media was fully utilized. Compared to the average current density at the macrometer scale, the detection of current density at the micrometer scale can provide more valuable insights for the rational design of low iridium loading membrane electrodes. Our study showed that the border area of the PTL/CL interface exhibited the highest current density under any working condition. The untouched area exhibited the lowest current density because of poor conditions, which can be attributed to the combined effects of conductivity and mass transfer. Besides, an experimentally validated numerical model was also proposed, and the simulation results indicated that porosity was the most critical factor affecting the performance of PEMWE, and the optimized porosity was about 0.3. The new insights into the localized current density derived from this method open a new avenue for optimization of the performance and cost reduction of PEMWE.
Unfortunately, in practical PEMWE devices, the tight assembly of their components using dense metal plates makes it challenging to accurately observe and probe internal information. Several detection methods have been reported to directly study the current density distribution of PEMWE and PEMFCs. Approximately 100 literature reports on the methods for analyzing the current density distribution of PEM cells were surveyed in this study; however, only 24.3% of those have focused on techniques for investigating the current density distribution of PEMWE, which included electrochemistry, printed circuit boards (PCBs),12–14 sensing wires,15,16 and high-speed and microscale visualization systems (HMVSs),17–20 as illustrated in Fig. 1(a). However, each of these methods has certain limitations. For instance, methods based on electrochemistry are restricted to analyzing current density distribution on a macro-electrolytic cell scale. Moreover, while the resolution of PCBs and sensing wires is higher compared with electrochemistry, their device assembly is difficult and does not reflect or give true information on real working conditions. The utilization of HMVSs allows for the identification of current density distribution by observing the information of bubbles, and it has revealed that the primary origin of current density can be traced back to the interface between a porous transport layer (PTL) and catalyst layer (CL). While these methods offer information on the micron-scale current density distribution, their qualitative analysis is generally insufficient for a comprehensive and profound understanding of the current density distribution under working conditions, especially in the contact area between the PTL and CL. This limitation hinders the loading reduction and efficient utilization of Ir-based catalysts, as shown in Fig. 1(b). Hence, there is a crucial need to develop novel approaches for assessing smaller scales of current density distribution under actual operating conditions.
RuO2 has garnered considerable interest as a catalyst for the OER in acidic conditions owing to its exceptional activity. However, its inherent instability poses a significant drawback, given that Ru species tend to separate as volatile high-valence oxides during the OER process,21–26 as confirmed by RRDE27 and eQCM28 tests, as shown in Fig. S1.† To study the stability of RuO2 under real working conditions, we conducted tests with a single PEMWE cell with RuO2 as the anode catalyst. Aging tests of RuO2 as the anode catalyst were carried out at 1.8 V, as shown in Fig. S2,† and we accidentally discovered that the electrochemical performance of the aged single cell could be partially recovered when the Ti mesh was rotated 90 degrees, as shown in Fig. 2(a and b). This phenomenon suggested the non-uniform distribution of current density over the CL surface. At the same time, EIS diagnostics were used to obtain Rct and Rohm, where the former represents the charge transfer for electrochemical reactions of the OER and the latter represents the proton and electron transport resistance. A 20 mV perturbation amplitude was used to ensure high-quality results. After the aging test, the amorphization of the catalyst or the dissolution of high-valent compounds led to a sharp increase in both Rct and Rohm. Simultaneously, significant interface resistance occurred in the high-frequency region, denoted as Rc, impeding electron transmission and diminishing the electrochemical performance of PEMWE, as illustrated in Fig. 2(b), S3, and Fig. S4.†29 When the Ti mesh was rotated 90 degrees, Rct, Rohm and Rc were moderately reduced, consistent with the polarization curve trend.
Fig. S5† shows the morphology of the fresh and aged CL as measured by scanning electron microscopy (SEM), as shown in Fig. S5(a).† Clear differences could be observed between the two samples, with the aged CL exhibiting regular indentations from the Ti mesh. To distinguish between the surface areas of the CL, the region pressed by the Ti fiber was labeled as “(2),” while the remaining area was denoted as “(1)”. Atomic force microscopy (AFM) was employed to further investigate the differences in the CL composition. AFM images of the fresh and aged CL were obtained at various locations and analyzed for assessing the topography, current mapping, and stiffness. In the topography images, the fresh and aged CL at different locations had a similar roughness, as shown in Fig. S5(b).† The electronic current at the CL was influenced by the content of RuO2 and the ionomer, with a higher RuO2 content leading to a higher electronic current. At a bias voltage of 1 V, the electronic current on the CL was detected and appeared grainy, reflecting the multi-component composition of the CL, including RuO2 and the ionomer, as shown in Fig. 2(d–f).30,31 The area of (1) had the second highest conductivity, while the area of (2) had the lowest conductivity, indicating that the content of RuO2 followed the opposite order. Stiffness mapping revealed that region (2) was stiffer than region (1), as shown in Fig. 2(g–i). The AFM data further revealed the non-uniform distribution of the current density on the surface of the CL.
In this work, we developed an ex situ visualization method to study the current density distribution at the micrometer scale by exploiting the unstable property of RuO2 under acidic OER conditions, as shown in Fig. 3(e). Scanning electron microscopy (SEM) and energy dispersive X-ray spectroscopy (EDS) were employed to evaluate the Ru element content at various locations on the CL. The region with a lower content of the Ru element exhibited a higher distribution of current density, attributable to a faster electrochemical reaction and a greater precipitation of the Ru element. Since IrO2 is a stable OER catalyst, the Ir element was utilized as the calibration element to determine the Ru element content on the CL surface.32,33 As shown in Fig. S6,† aging tests of IrO2 + RuO2 as the anode catalyst were performed at 1.5, 1.8, and 2.1 V, respectively, until the current density reached 0.1 mA cm−2. Obviously, the cell voltage at 1.5 V took the longest time, lasting 16 h. Also, the endurance time was 2 h at 1.8 V and 2.1 V. The content of Ru element was determined at seven symmetrical positions at the interface between the Ti fiber and the aged CL using SEM and EDS, as illustrated in Fig. 3(d) and S9.† In order to better analyze the experimental results, the PTL/CL interface was classified into three distinct regions: the untouched area, border area, and touched area. Among these, the untouched area consisted of site1 and site7, the border area comprised site2 and site6, while the touched area encompassed site3, site4 and site5, as shown in Fig. 3(e).
The Ru element content on the CL surface was analyzed in the aging test at 1.5 and 1.8 V as depicted in Fig. 3(a and b). Site2 and site6 exhibited the lowest Ru element content and highest current density distribution for the CL, which were located at the border area of the PTL/CL interface, fully meeting the requirements for electron/proton conduction and water/gas two-phase flow transport.34,35 Site1 and site7 exhibited the highest Ru element content due to the poor conductivity of CL, which was consistent with conclusion drawn in other literature.36,37 Additionally, it was observed that the Ru content initially increased and then decreased from site3 to site5, with site4 exhibiting the highest content. The touched area (site3, site4, and site5) at the PTL/CL interface was compressed due to compression force, leading to challenges for oxygen removal and water delivery. Consequently, site4 exhibited a higher Ru element content. Similar trends were observed in the aging test at 2.1 V, except for differences in the Ru element content in the touched area, which exhibited a higher Ru element content compared to in the untouched area. The difference in Ru content in the touched area for the high aging voltage was attributed to the distinct supply of reactant water. Under high aging voltage, the OER reaction rate was extremely fast and a large number of oxygen bubbles accumulated, leading to catalyst utilization in the CL being relatively uneven due to the poor water supply.38–40
For further quantitative analysis of the current density distribution, a numerical model was developed (see Fig. S10†), and the model information is shown in Table S1.† Fig. S11† presents the predicted porosity and electric potential distribution at the CL surface. Compared with the untouched area of the PTL/CL interface, the pores within the CL in the touched area shrank or even closed entirely due to the compression force from the titanium fibers, making it challenging for water/gas to effectively be transported, as shown in Fig. S11(a).† Besides, due to the significantly lower electrical conductivity of the CL compared to the PTL, the potential dropped rapidly near the titanium fibers, with minimal electric potential in the untouched area, as shown in Fig. S11(b).† To comprehensively analyze the influence of porosity and conductivity, we introduce the concept of a reaction rate index (RRI, where RRI = current density × mass concentration of water) to reveal the reaction rate on the CL surface. Hence, the RRI along the horizontal direction (at a specific depth in the CL domain) was quite different, with the highest values observed in the border area, followed by the touched area, and the lowest values found in the untouched area, as shown in Fig. 4(b). Next, the influence of the porosity and electric conductivity of the CL on the RRI were further studied, and the simulation results are shown in Fig. 4(c and d). The simulation results indicated that, given the considerably lower conductivity of the CL compared to the PTL, variations in the conductivity of the CL had a negligible effect on the distribution of the electric field, leading to a nearly constant RRI value.36,37 While the porosity had great impact on the RRI, the value of the RRI increased and then decreased for porosities ranging from 0.1 to 0.7, with the highest value attained at 0.3, which resulted from the competition between mass transfer and the active catalyst sites. Simultaneously, the RRI value experienced a rapid increase within the initial few milliseconds, followed by a stabilization of the growth rate over the subsequent tens of milliseconds. This phenomenon arose from the significantly faster electron transport compared to mass transport (by five to six orders of magnitude).
Drawing from the experimental and simulation data, a schematic representation of the current density distribution at the PTL/CL interface is presented in Fig. 4(a). Since the voltage of the catalyst layer (U) remained constant in the chronoamperometry tests, the current density was determined by the impedance R, including Rct and Rlimit. Among these, Rct remained constant while Rlimit was composed of Rohm and Rmt, as shown in eqn (1).
(1) |
At the border area of the PTL/CL interface, the electric charge conduction and water/gas two-phase flow transport occur more easily, resulting in the lowest impedance and highest current density distribution. However, for the untouched area between the PTL and CL, due to the poor in-plane conduction of the CL compared to the PTL, the catalyst cannot be efficiently utilized. Besides, mechanically induced stress and strain by the PTL and continuous creep under hydrated environments results in micro-crack formation, which hinder utilization of the catalyst in the untouched area.35 Therefore, the current density distribution was lowest due to the higher Rohm. The touched area of the current density distribution thus varies depending on the operating conditions, such as the applied cell voltage and assembly compression force, as so requires simultaneous consideration of the water/gas transport and electron conduction, representing the interplay between Rohm and Rmt.
In summary, a straightforward and ex situ visible method to identify the current density distribution at the micrometer scale in the CL surface by utilizing the poor stability of RuO2 was proposed, and a numerical model was also developed, which was validated by the test results. The experimental and numerical results revealed that the distribution of the current density of the CL surface was uneven, indicating that the catalyst was not fully utilized. The border area of the PTL/CL interface exhibited the highest current density under any working conditions, due to meeting the requirements for electron/proton conduction and mass transport. The untouched area exhibited the lowest current density due to the poor conductivity of the CL. The untouched area exhibited varying current density under different operating conditions, which could be attributed to the combined effects of conductivity and mass transfer, such as the applied cell voltage and assembly compression force. Therefore, identifying the current density distribution at the micrometer scale provides valuable guidance for reducing precious metal loading, and designing PTL/CL interfaces more efficiently.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4cy00031e |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2024 |