Diffusion Length Mapping for Dye-sensitized Solar Cells

The diffusion length (L) of photogenerated carriers in the nanoporous electrode is a key parameter that summarizes the collection efficiency behavior in dye-sensitized solar cells (DSCs). At present, there are few techniques able to spatially resolve L over the active area of the device. Most of them require contact patterning and, hence, are intrinsically destructive. Here, we present the first electron diffusion length mapping system for DSCs based on steady state incident photon to collected electron (IPCE) conversion efficiency (η IPCE) analysis. The measurement is conducted by acquiring complete transmittance (T DSC) and η IPCE spectra from the photo electrode (PE) and counter electrode (CE) for each spatial point in a raster scan manner. L(x, y) is obtained by a least square fitting of the IPCE ratio spectrum (IPCE R = η IPCE-CE η IPCE-PE). An advanced feature is the ability to acquire η IPCE spectra using low-intensity probe illumination under weakly-absorbed background light (625 nm) with the device biased close to open circuit voltage. These homogeneous conditions permit the linearization of the free electron continuity equation and, hence, to obtain the collection efficiency expressions (η COL-PE and η COL-CE). The influence of the parameter's uncertainty has been quantified by a sensitivity study of L. The result has been validated by quantitatively comparing the average value of L map with the value estimated from electrochemical impedance spectroscopy (EIS).


Introduction
It has long been recognized that the spatial distribution of local parameters, such as optical absorption, electroluminescence, diffusion length, etc., can provide valuable information on the stability, performance and degradation prediction of solar cells [1,2].Moreover, spatial visualization of these properties is crucial to better address the up-scaling procedure that every photovoltaic technology has to tackle [3].Mapping techniques can be grouped into ex situ and in situ methods according, respectively, to the need for cell disassembly or not [4].We will focus our attention on in situ mapping techniques that are based on optical probes, such as Raman [5,6], transmittance-reflectance and photo-luminescence [7] microscopy and electro-optical techniques [8], such as light-beam-induced current (LBIC) [9][10][11] and electro-luminescence [12] microscopy.The good performances of DSCs require improving the charge collection efficiency over the entire semiconducting electrode (generally, titanium dioxide (TiO 2 )) to the contact [13,14]; in particular, the electronic property and the surface area of the photoanode determine the current output of the device [13,15].Accordingly, good light harvesting ability, charge transport and low recombination are key challenges for hybrid solar cells [16,17].The charge collection capability of the photoanode is quantified with the diffusion length of the injected charges; recent papers present an LBIC-based technique [18][19][20] to directly extract the absolute value of L; unfortunately, despite the simplicity of the analysis, these techniques require a substrate contact patterning that limits the application to customized devices.Here, we present a new in situ technique, spectrally-resolved analysis by transmittance and efficiency mapping (SATEM), to extract the (L) distribution by contemporary mapping the transmittance (T DSC ) and the incident photons to current conversion efficiency (η IPCE ).As such, SATEM provides 2D maps of both electro-optical and topographic information simultaneously.The L estimation is based on an accurate modeling of the incident photon to collected electron (IPCE) spectra, as already discussed by Barnes et al. [21,22] and Jennings et al. [23].Therefore, spectrally-resolved mapping of the IPCE ratio (IPCE R ) and T DSC contains more information on physical and electrochemical solar cell properties than classical LBIC techniques.We introduced SATEM analysis in a recent article [24] focused on the degradation mechanisms induced by the reverse bias condition on the large area DSC module (five cells of 3.6 cm 2 each); spatially-resolved transmittance and IPCE spectra were combined with the resonance Raman technique to investigate the effect of by-products on the regeneration efficiency of the dye and the free charge recombination losses inside the TiO 2 film.Moreover, the ability to tune both the background and monochromatic beam was used to highlight areas with trap-limited charge collection and the limitation of iodide diffusion.Differently, in the present work, we demonstrate that the combination of an SATEM analysis with an accurate modeling can provide the spatial distribution of L. In order to confirm the results, electrochemical impedance spectroscopy (EIS) has been chosen as the independent determination of L; in particular, the average value of L provided by SATEM analysis (the steady state technique) is in agreement with the value extracted from the fitting of the EIS spectra (the dynamic technique).

Theoretical Basis
The diffusion length, considered normal to the plane of the film, is the average distance an injected electron can travel through the cell before recombining with tri-iodide ions or oxidized dye species.If L is lower than the TiO 2 film thickness (L < d), only a small injected charge fraction will be collected, making a long diffusion length desirable [25].Assuming a first order recombination in free electron concentration, the electron diffusion length can be defined as L = √ D n τ n , where D n is the diffusion coefficient of electrons in the conduction band of the TiO 2 and τ n is the free electron lifetime.

Estimation of L from Electrochemical Impedance Spectroscopy Measurements
Following the suggestion of Jennings et al. [23], EIS analysis has been chosen to quantitatively validate the SATEM results.In particular, EIS experiments can be performed under the same homogeneous conditions of SATEM (open circuit voltage bias condition with weekly-absorbed background light), unlike other dynamic techniques (transient photocurrent and intensity-modulated photocurrent spectroscopy [26]) that operate under short circuit conditions.The diffusion length of DSCs is well defined from the impedance Z of the cell estimated by EIS [14,[27][28][29].The EIS technique is based on the superimposing of a small amplitude harmonic AC voltage modulation on the DC voltage of the cell for a set of frequencies and by measuring the resulting AC current.The final impedance spectrum can be fitted with a suitable electrical circuit model generally based on a transmission line that models the TiO 2 layer [30].The model used in our analysis (Figure 1) is discussed in detail in several works [30][31][32].The distributed resistance r t (total resistance R t = r t d) represents the local resistance to electron transport in the TiO 2 layer that is related to the carrier concentration (n c ) and diffusion coefficient (D n ).The distributed component r ct (total resistance R ct = r ct d ) models the local resistance to the charge transfer across the TiO 2 -electrolyte layer that is determined by n c and the free electron lifetime (τ n ).The diffusion length is then defined as: Since the measured R ct and R t are an averaged value of the local contributions over the TiO 2 thickness, the homogeneous conditions are indispensable to obtain a reliable estimation of the collection capability of the photoelectrode and, hence, of the absolute value of L. In particular, EIS analysis under illumination at open circuit allows characterizing r t and r ct under the flat quasi-Fermi level unlike what happens under dark and forward bias [17,26,33].The disadvantage of the impedance method is that it becomes difficult to measure the transport phenomena at high light intensities, particularly in the case of cells with long diffusion lengths, where r t becomes much smaller than the series and cathode impedance at open circuit voltages and fitting becomes unreliable.

Estimation of L from Incident Photon to Collected Electron Measurements
The η IPCE is composed of three components, such that: where j SC is the short-circuit current, q is the elementary charge, φ is the photon flux density, η LH is the efficiency of photon harvesting by dye molecules, η sep is the charge separation efficiency and η col is the electron collection from the TiO 2 to the external circuit [14].
In semitransparent DSCs, the ones discussed in this work, η IPCE can be measured from the photo electrode (PE) or counter electrode (CE) side.In the following, we differentiate the parameters obtained with PE or CE side illumination by the PE or CE subscript, respectively.As proposed by Halme et al. [34], a simple analytic expressions of η LH,PE and η LH,CE can be derived from the optical model shown in Figure 2.
where T TCO , T CE and T EL are the transmittances of the TCO-coated glass substrate of the photoelectrode, the counter electrode and the free-electrolyte, respectively.R PE is the photoelectrode film reflectance, and d is its thickness.α is the absorption coefficient of the electrolyte-filled dyed photoelectrode film: where α D is the absorption coefficient of the dyed photoelectrode film, α EL is the absorption coefficient of the bulk electrolyte solution, P is the porosity of the film (assumed to be 0.5) and is the average optical mean path length parameter that takes into account light scattering phenomena (a value of 1.5 is assumed here).Analytic expressions of the collection efficiencies are obtained by the resolution of the steady state free-electron continuity equation: where η reg is the sensitizer regeneration efficiency expressed as: where n c is the equilibrium concentration of free electrons in the dark, G(x) is the position-dependent photogeneration rate, D n is the diffusion coefficient of free electrons, k n and k g are rate constants for recombination of electrons with tri-iodide and oxidized sensitizer molecules, respectively, k r is the rate constant for the sensitizer regeneration reaction with iodide, γ n and γ g are the reaction orders with respect to free electron concentration and n eqm is the free electron concentration at thermal equilibrium in the dark.

Linearization of the Free Electron Continuity Equation
Classical expressions of η col,PE-CE are obtained from a linearization of the free electron continuity equation, where a first order recombination with tri-iodide (γ n = 1) and a negligible recombination with oxidized dye (k g = 0) are assumed [25,32,35].These assumptions are the main cause of the discrepancies between different methodological approaches that estimate the diffusion length parameter [22,23,34,36].For a position-independent background generation rate, G bg , and electron concentration (i.e., open-circuit condition), n bg , the free electron continuity equation for n bg can be written as I − : Then, the background generation rate can be written as: Assuming small perturbations of G(x) around G bg and n c (x) around n bg , Equation ( 6) can be linearized resulting in a familiar equation with a linear recombination term: where: From Equation ( 8), it follows that existing solutions to the electron continuity equation remain valid, even for γ n = 1 and non-negligible recombination with oxidized sensitizer molecules.This conclusion applies to the employed characterization techniques; therefore, if the provided measurements are performed under homogeneous conditions, the expressions of electron collection efficiency η col,PE-CE can be described by the standard diffusion model of electron generation, transport and recombination at the nanostructured photoelectrode introduced originally for DSCs by Sodergen et al. [25]:

.2. Incident Photon to Collected Electron Ratio
Considering that η sep can be assumed independent upon the light direction (PE or CE side) [21,34], it follows from Equations (3), ( 4), ( 9) and (10) that the IPCE R is given by: The expression can be improved by introducing the analytic expression T DSC , that is experimentally obtained from SATEM: where T TCO|EL|CE = T TCO T EL T CE represents the transmittance of a cell composed of the CE and electrolyte, and its spectrum is extracted and averaged over the perimeter points around the active area of the cell.IPCE R can be rearranged as follow: The spatially-resolved absorbance coefficient (α) is obtained by inverting Equation ( 12): On the contrary, T TCO , R PE and d are provided by independent measurements, but it is important to reiterate that: the spatial variation T TCO can be assumed negligible; the spatial variation of d is averaged by increasing the spot size of the optical probe during the spectral resolved mapping; and as will be shown in the sensitivity analysis, the influence of R PE over the estimation of L is demonstrated to be irrelevant.Equation ( 13) can be fitted to experimental IPCE R with L and d as the only free-fitting parameters.By allowing also d to vary around the experimental value, a better fitting to the spectra is obtained without the influence of the value of L [23].

Scanning Apparatus (Spectrally-Resolved Analysis by Transmittance and Efficiency Mapping)
The SATEM system (Figure 3) is built around an inverted microscope DMI 5000 from Leica (Wetzlar, Germany) that includes a charge-coupled device (CCD) camera, a motorized XY stage and two optical entrances; one of the latter is dedicated to the monochromatic light coming from a Xenon lamp (200 W Apex Source Model 66450) and a monochromator ( 18 m Cornerstone 130) from Newport, Irvine, CA, USA.The spectral resolution is 2 nm on a range from 300 nm to 1000 nm.The sample is completely illuminated from the back-side with a red light-emitting diode array (Oslon ILR-ON 625 nm 20 W from Osram, Munich, Germany) that allows one to fix an irradiation level up to the 1.5 sun in the visible range of 2 cm 2 .The absorption length calculated at 625 nm over the TiO 2 film is ca. 10 µm, approximately equal to the commonly-used layer thicknesses.The background light of the red led array is fine-tuned in order to ensure the same V oc for PE and CE illuminations.The transmitted optical signal can be registered either by a large area photodetector (PD 2 , Model FDS 1010 from Thorlabs, Newton, NJ, USA) close to the back side of the device or an integrating sphere, depending on the device layout and scattering phenomena.The short circuit currents of the device and PD 2 are discriminated in a phase-sensitive detection system composed of an optical chopper (Newport Model 75159) and two digital lock-in amplifiers (Eg&g 7265 from Signal Recovery, Oak Ridge, TN, USA).Frequency modulation in the sub-Hz range is necessary to avoid the underestimation of the IPCE value [37].An embedded transimpedance amplifier connected to PD 1 and a data acquisition board (Model 9205 from National Instruments, Austin, TX, USA) provides the incident optical power information.Thanks to the simultaneous acquisitions of the incident (PD 1 ) and transmitted (PD 2 ) optical powers and short circuit current of the device, spectra and maps of both IPCE and transmittance will be independent of possible fluctuations of the monochromatic light source.The system is completely automated under the LabVIEW programming environment, and the IPCE/transmittance spectra are acquired and elaborated in parallel.The advantages of the SATEM approach over previous spectroscopic and scanned probe studies are: (i) direct spatial correlation between optical (transmittance) and electro-optical (η IPCE ) information; (ii) the ability to extract wide spectra (350-1100 nm) with 2 nm of resolution; (iii) differential IPCE with localized white bias light (tunable up to 1.5 sun) to investigate electron collection and mass transport limitations [38][39][40]; (iv) sub-Hz optical chopping frequency to avoid the underestimation of the IPCE spectrum [41][42][43][44]; and (v) lamp power fluctuation immunity thanks to simultaneous acquisition of incident and transmitted optical signals and the short circuit current of the sample.

Results and Discussion
In order to focus attention on the potential of the technique rather than the specific device under study, a standard and stable DSC configuration has been chosen; in particular, a commercial iodine-based electrolyte (HSE) and Z907 sensitizer [45] from Dyesol (Queanbeyan, Australia) ensured an excellent stability that helped during extended measurements.

Sensitivity Analysis
Since the spatial estimation of L requires a nontrivial fitting procedure with several parameters involved, we conducted a sensitivity study in order to quantify the influence of parameter uncertainty (e x i ) on the L. In particular, x i stands for d, α, R PE , T DSC , T TCO , η IPCE,PE and η IPCE,CE .The SATEM system acquires T DSC , η IPCE,PE and η IPCE,CE for each spatial point; therefore, its contribution over the L uncertainty (e L ) is related only to the measurement errors.Film thickness d is equal for each spatial point and assumed equal to the average value of the thickness 1D profile obtained with the profilometer system; this is a reliable assumption because diffusion length mapping is time consuming, and hence, a large step and spot sizes are adopted (600 × 600 µm 2 ).R PE can be assumed constant due to the extremely low sensitivity value.S L T TCO reaches high values, but assuming an uniform deposition of the TCO film, its effect can be neglected.The absorption coefficient of the film (α) is evaluated by Equation ( 14) where T TCO|EL|CE is evaluated outside the active area where only the CE and electrolyte are present; several acquisitions with a large spot area have demonstrated that T TCO|EL|CE can be assumed constant.The relative error of L is expressed as e L = ∑ i S L x i e x i , where S L x i are the sensibility functions of L with respect to the involved parameters (x i ) and are defined as follows: where x i0 and L 0 are the nominal values.Due to the implicity of Equation ( 13) with respect to L, Equation ( 15) cannot be directly used.By considering an implicit form equation as: F(x, L) is a continuous and differentiable function around nominal values (x 0 , L 0 ).Therefore, by implementing the implicit function theorem, it is possible to express S L x i as: As shown in Figure 4, L has the opposite sensitivity trends with respect to η IPCE,PE and η IPCE,CE .The variation in the thickness of the film is neglected due to the reason that the thickness variation has been averaged in the spot light area.Between 470 nm and 580 nm, S L T TCO and S L T DSC are around +2.5 and −1.5; thus, an equal error of +e% in the evaluation of these parameters becomes respectively +2.5e% and −1.5e% of error on L. S L R PE is approximated to R PE , which means very small values over the visible spectra.For the accuracy of the absorption coefficient, α can be calculated spatially by the total transmittance of the cell and the transmittance of TCO/EL/CE.The most important indication coming from Figure 4 is that the sensitivities of L are confined to acceptable values over a range between 450 nm and 550 nm.

Estimation of L by Spectrally-Resolved Analysis by Transmittance and Efficiency Mapping
The SATEM technique and the fitting procedure described above have been applied to investigate the influence of spatial inhomogeneities (fabrication defects, deposition problems, etc.) present in a typical DSC cell (Figure 5a).To investigate such inhomogeneities, which produce a variation in the electro-optical parameter of the cell, we performed a preliminary mapping at fixed wavelength (e.g., 530 nm) with a spatial resolution of 50 µm.Figure 5b-d shows the corresponding IPCE PE , IPCE CE and absorbance maps.Interestingly, all of the maps show some uncorrelated trends; in particular, IPCE CE is characterized by a relatively higher efficiency on the center of the active area.Focusing on the main goal of the paper, a spectrally-resolved analysis was performed under homogeneous conditions on a reduced number of points (64) with a spot size increased to 650 × 650 µm 2 .As discussed in the sensitivity analysis, this trick is adopted to average spatial fluctuations of film thickness (d) and T TCO transmittance.
The spectral range between 450 nm and 650 nm with 2 nm of resolution was investigated.Underestimation of the IPCE values [41,42] was avoided by modulating the monochromatic beam at 0.33 Hz with an integration time of 10 s, resulting in a total acquisition period of 7.1 h for one side of the cell.This is responsible for the higher absolute values of η IPCE compared to the high resolution scan (Figure 5b,c).Figure 6a,b shows the IPCE PE , IPCE CE and absorbance spectra that show an evident variation of IPCE CE and absorbance.IPCE R spectrum fittings are performed between 450 nm and 550 nm in order to minimize the sensitivity of L. For clarity, only three spatial points have been chosen to highlight the fitting results (Figure 6c).The normalized root mean square (NRMS) error for each point is plotted in the inset showing an average value of 0.95 that confirms the goodness of fit.The resulting diffusion length map (Figure 6d) reports an increase up to 19 µm located in the center of the device and several points with values lower than the film thickness (d = 11 µm).The histogram reports an average value of 12.1 µm that confirms the low η IPCE recorded on the device.The estimated L(x, y) can be used to readily calculate η col spectra from Equations ( 9) and ( 10) leaving η sep the only unknown parameter in Equation ( 2).As already reported by Halme et al. [34] and Jennings et al. [23], a second fitting procedure with the experimental η IPCE for CE and PE illumination allows one to estimate the wavelength-dependent η sep .Figure 7a shows the experimental (dot) and fitted (line) of η IPCE for a single spatial point (Point "H" reported on Figure 6d); the estimated η sep is reported in Figure 7b already with the calculated η col and η LH .The low η sep and the short L values are the main causes of the poor performance of the device under investigation.Unfortunately, the strong wavelength dependence of η sep cannot be uniquely associated with the injection or regeneration problems of the sensitizer; anyway, two causes can be hypothesized: a mismatch between the excited-state energy levels of the dye and the electron acceptor states in the TiO 2 (injection process) [23] or/and a poor regeneration efficiency that results in a strong dependence of the regeneration process with light intensity or [46,47].

Estimation of L by Electrochemical Impedance Spectroscopy
The diffusion length value obtained by the IPCE R technique was compared with EIS data fitted with the equivalent circuit shown in Figure 1.Measurement was performed under the same homogeneous circumstances as for the IPCE test (red light illumination and open circuit condition).In order to improve the quality of the fitting, all of the capacitors were replaced by constant phase elements (CPE) with an exponent ranging from 0.85 to 0.99.Interestingly, the EIS spectrum reported in Figure 8 resembles that of a Gerischer impedance instead of a classic transmission line [14] that is necessary to evaluate R ct and R t .However, as already demonstrated [23,34], the transmission line model can still be adopted for L > d/2.In this study, a slightly higher value of L (15.1 µm) with respect to the SATEM analysis is obtained.This result confirms the potentiality of the SATEM technique to map the distribution of the diffusion length over the sample accurately.

Materials and Methods
A 0.25 cm 2 active area cell was fabricated with transparent nanoporous TiO 2 paste (DSL 18NR-T) from Dyesol screen-printed through a 43T mesh screen onto fluorine-doped tin oxide (FTO)-coated glass substrate (NSG TEC 8/from Pilkington Group Limited, St Helens, UK).After 15 min of drying at 80 • C, TiO 2 films were sintered at 525 • C for 30 min.The resulting thickness of the TiO 2 films was measured with a surface profilometer Dektak 150 from Veeco and was observed to be 6.7 µm.TiO 2 electrodes were dyed for 16 h in a 0.3 mM solution of Z907 from Dyesol.After dyeing, the PEs were rinsed with ethanol in order to remove the excess dye before the assembling process.Platinized CEs were made by screen-printing platinum paste through a 100T mesh screen onto FTO-coated glass substrate.The CEs were dried at 80 • C and then sintered at 525 • C for 30 min.The two electrodes were laminated with Bynel from Solaronix (Aubonne, Switzerland).After the hot-melting step, the distance between the two electrodes was measured to be about 45 µm.The used electrolyte was a commercial iodine-based electrolyte HSE from Dyesol.

Conclusions
SATEM is the first scanning apparatus that implements the steady state IPCE R technique under homogeneous conditions for an indirect estimation of the diffusion length in a DSC.The current study is believed to give insight into the system design and analytic procedure for reliable estimates of spatially-resolved L. The homogeneous conditions are the main peculiarity of the technique.In particular, IPCE PE , IPCE CE and absorbance spectra were acquired with uniform background generation (bias light at 625 nm) and uniform background electron density (open circuit voltage).We highlighted the most delicate parameters through a depth sensitivity analysis of L; in particular, the high value of the sensitivity with respect to the film thickness (S L d ) has been mitigated by increasing the optical spot size to 650 µm, resulting in a lower spatial resolution of the L distribution.This point can be easily improved by arranging an experimental d(x, y) with a profilometer measurement over the photoanode instead of the average value.Moreover, a better design of the optical line will give more room of improvement.Good agreement was obtained with the integral value of L estimated by an EIS analysis under the same operative conditions, confirming the remarkable setting of SATEM analysis.We believe that these results provide a solid base to extend the technique to different semi-transparent technologies, such as perovskite-based cells [48][49][50] and organic cells [51][52][53]).Unfortunately, a role in transport and recombination can be played by other contributions (for instance drift) that may have to be taken into account for precise collection efficiency modeling [54,55].However, it has been demonstrated that the electric field contribution in organic solar cells arises mainly at a high level of irradiation (1 sun) as a consequence of a charge accumulation in the device [56], and we suppose it to be negligible in our experimental conditions.In perovskite solar cells, a recent article [57] demonstrated the field free condition of the methyl-ammonium lead iodide (MAPI) layer under short and open circuit conditions; moreover, as already assumed in several works [58][59][60], we believe that a diffusion-based transport model could be a starting point for an investigation with the SATEM analysis.Further research is required in this direction to unravel the limiting mechanisms of the IPCE and optical-based analysis for the extraction of the charge collection capabilities of hybrid-organic solar cells.

Figure 1 .
Figure 1.Small signal equivalent circuit of dye-sensitized solar cells (DSCs) used to fit electrochemical impedance spectroscopy (EIS) spectra.The transmission line is outlined by the red box; r t and r ct are discussed in the main text.The model includes also: distributed capacitance of TiO 2 (c µ ), series resistor (R s ), charge-transfer resistance (R TCO ) and double-layer capacitance (C TCO ) at the exposed transparent conducting oxide (TCO)-electrolyte interface, resistor (R CO ) and capacitor (C CO ) at the TCO-TiO 2 interface, Warburg element (Z d ) for Nernst diffusion of I − 3 in the electrolyte, charge-transfer resistance (R pt ) and double-layer capacitance (C pt ) at the platinized counter-electrode.

Figure 2 .
Figure 2. Optical layer structure of a typical DSC.PE: photo electrode; and CE: counter electrode.

Figure 3 .
Figure 3. Illustration of the set up.

Figure 4 .
Figure 4. L diff sensitivity functions with respect to α, d, T TCO , R PE and T DSC .For S L diff α , S L diff d and S L diff R PE , we have set d = 11 µm, L = 21 µm and assumed a typical film reflectance spectra.

Figure 5 .
Figure 5. (a) Photo of the DSC with 0.25 cm 2 of active area.Spatial maps of IPCE PE , (c) IPCE CE and (d) absorbance acquired at λ = 530 nm.

Figure 6 .
Figure 6.(a,b) IPCE PE and IPCE CE and absorbance spectra of 64 spatial points of the active area; (c) IPCE R fitting results from 450 nm and 550 nm (points are experimental values) of three spatial points.The normalized root mean square (NRMS) error for each point is plotted in the inset, showing an average value of 0.95 that confirms the goodness of fit.(d) Spatial map of the estimated L; the normal distribution reported below shows an average value of 12.1 µm.

Figure 7 .
Figure 7. (a) Fitting of the η IPCE PE and η IPCE CE by varying the η sep spectrum for Point "H" of Figure 6d.(b) Estimated η sep and calculated light harvesting and collection efficiencies for the both the PE and the CE side.

Figure 8 .
Figure 8. Impedance spectra under the open circuit condition (0.51 V).Fitting was achieved using the model shown in Figure 1 with an extracted diffusion length of 15.1 µm.