Analytical Model for Current–Voltage Characteristics in Perovskite Solar Cells Incorporating Bulk and Surface Recombination

The effects of surface recombination on the steady-state carrier profiles and photocurrent in perovskite solar cells are investigated in this paper. The continuity equations for both holes and electrons are solved considering carrier drift and diffusion under the exponential carrier generation profile in the perovskite layer and considering both bulk and interface carrier recombination. An analytical expression for the solar-induced photocurrent is derived. The rate of carrier recombination at the interfaces has a very significant effect on the carrier profile, photocurrent, and, hence, on the charge collection efficiency. The external current density is calculated considering the dark current and nominal solar spectrum-induced photocurrent. The proposed model is fitted and verified with published experimental results from various publications. The fittings of the model with experimental results provide information about the interface and bulk charge carrier transport parameters.


Introduction
The halide perovskite solar cells (PSCs) have received huge attention in recent years [1,2] due to their higher power conversion efficiency (PCE) and cheaper manufacturing cost.The PCE of PSCs has exceeded 25% in single-junction cells, which is comparable to highly researched silicon solar cells [3,4].The perovskite materials possess several intrinsic qualities, such as a high absorption coefficient over the entire sun spectra and good charge carrier transport properties.
The general form of the halide perovskite materials is ABX3 (where A = CH3NH3 (abbreviated as MA+), HC(NH 2 ) 2 (abbreviated as FA+), and Cs; B = Pb 2+ or Sn 2+ ; and X is a halogen atom such as Cl − , I − , or Br − ).The halide perovskites, such as MAPbI 3 , MAPbI 3−x Cl x , and (FAPbI 3 ) x (MAPbBr 3 ) 1−x , have been examined for photovoltaic solar cell applications.The bandgap of these perovskites varies within 1.5-1.6 eV [5,6], which is appropriate for photovoltaic solar cells.Both crystalline and polycrystalline forms of these perovskites can be used.
The perovskite photovoltaic solar cells fall within the thin-film photovoltaic structure category.The general structure is transparent electrode/HTL/bulk perovoskite/ETL/ electrode type (where HTL and ETL stand for hole transport layer and electron transport layer, respectively).The top contact is made of transparent conductive oxide, and the bottom contact can be a metal.The bulk perovskite layer is an undoped (intrinsic type or slightly n or p type) layer, and its thickness is less than 1 µm [7].The ETL and HTL are much thinner than the bulk perovskite layer.The structure is equivalent to ETL-i-HTL or HTL-i-ETL.The conventional structure is the ETL-i-HTL type, which provides higher efficiency.On the other hand, the HTL-i-ETL structure provides better stability and lower cost, with lower efficiency.Though the ETL and HTL are not conventional n and p layers, the carrier transport properties of these layers are like n and p layers.This means that the electron transport is much more favorable than the hole transport in ETL, and the opposite is true in HTL.The energy band of the ETL and HTL must align with the perovskite layer so that photogenerated electrons and holes can easily pass towards the corresponding electrodes.The ETL and HTL also play a very crucial role in the interface properties of the perovskite layer and, thus, the control carrier profile across the perovskite layer and the photocurrent as well.It is very crucial to investigate the quantitative effects of the interfaces on the current-voltage (J-V) characteristics of these photovoltaic cells.
There are a limited number of papers found in the literature that describe a physicsbased model for calculating the photocurrent in perovskite solar cells [8][9][10].These publications did not consider the combined effects of several main charge carrier transport mechanisms, such as carrier diffusion, bulk, and surface recombination, for calculating the photocurrent.Many previous publications [11][12][13][14] experimentally demonstrated that the PCE highly depends on the transport and contact properties of the ETL and/or HTL layers by changing the composition of these ETL and HTL layers.However, as per my knowledge, only one paper in the literature [15] proposed a physics-based analytical model for studying the current-voltage characteristics in perovskite solar cells incorporating a few important charge carrier transport mechanisms, such as carrier diffusion, bulk, and surface recombination.However, the mathematical derivation of the previous model was based on a few unrealistic assumptions, such as the same surface recombination velocity for both electrons and holes at both interfaces and the same carrier transport properties for both holes and electrons.
In this paper, the previous model [15] has been extended by incorporating unequal values of surface recombination velocity for electrons and holes.The effects of surface recombination on steady-state electron and hole profiles are examined by solving the carrier continuity equation across the perovskite layer and by varying the surface recombination velocities at both interfaces.The theoretical model is compared with published experimental results, and hence, the physical mechanisms for the variation in the charge collection efficiency with the change in ETL and HTL are investigated.The charge carrier transport properties are also extracted.

Mathematical Model
The sun excitation creates an exponential distribution of electron-hole pairs (EHPs) across the bulk perovskite layer, and thus, the carrier generation rate G is the following [16]: and where E ph = hc/λ is the incident photon energy, h is the Plank constant, c is the speed of light, λ is the wavelength of the incident photon, x is the distance from the radiationreceiving electrode, R is the refection loss from the front surface, G 0 (λ) is the EHP carrier generation rate at x = 0, α(λ) is the absorption coefficient of the perovskite material, and I 0 is the intensity of the sun light.The sunlight generated electrons and holes drift under the built-in electric field.The carriers may also move by the diffusion process.Some of the carriers are lost during their travel towards the charge-collecting electrodes.The monomolecular recombination is the dominant carrier loss mechanism because of the low bimolecular recombination coefficient in perovskites and relatively low light intensity [17].Since our aim is to develop a compact analytical model, we make a few assumptions, such as [18] (1) the electric field F is approximately uniform [19] across the relatively thin perovskite layer under normal sun exposure, (2) the drift mobility (µ) and carrier lifetime (τ) in the active layer for each carrier type (electrons and holes) are uniform, and (3) the intrinsic carrier concentration for the relatively large bandgap (~1.5 eV) perovskites is negligible compared to the photogenerated carrier concentration.
The excess carrier continuity equations for holes and electrons in the bulk of the ETL-i-HTL structure under steady-state condition are the following: where δ and D are the photogenerated excess carrier concentration and diffusion coefficient, respectively, while the subscripts n and p represent electrons and holes, respectively.As mentioned previously, the ETL and HTL are not conventional n and p layers of homogeneous material.The electron concentration at the top surface (x = 0) of the perovskite layer may not be equal to the electron concentration of the ETL; rather, it should be determined by solving the recombination current density at the surface.The same argument applies to the hole concentration at the bottom interface (x = W).Therefore, all the boundary conditions for the carrier concentrations are determined by the surface recombination current.The boundary values of hole and electron current densities are the following [9,20,21]: where q is the elementary charge, W is the thickness of the perovskite layer, J is the conduction current density, and S is the surface recombination velocity.The first subscripts n and p on S represent electrons and holes, while the second subscripts t and b represent the top and bottom interfaces, respectively.Considering the boundary conditions in Equation ( 5) and solving Equations ( 3) and ( 4), the excess carrier concentrations are the following: where The photocurrent density for the incident photon of wavelength λ can be conveniently determined using the Shockley-Ramo theorem [22,23], which is given by the following: The entire sun spectra contribute to the photocurrent and, thus, the total photocurrent density: Most of the external voltage drops can occur across the perovskite layer because the very thin ETL and HTL are more conductive than the perovskite layer.Thus, the expression of the electric field F is the following: where V bi is the built-in potential, R s is the effective series area resistance, and V and J are the external load voltage and current density, respectively.The external current density from a solar cell is given by the following [24,25]: where J d (V) and R p are the dark current density and the shunt area resistance, respectively.The expression of the dark current density is given by the following [15,26]: where V t = k b T/q, T is the absolute temperature, and k b is the Boltzmann constant.The expression of dark current (Equation ( 12)) has a voltage-independent constant term J c , which is considered as a fitting parameter.

Results and Discussion
The Equations ( 8) to (12) are simultaneously solved to find the net current density J at a particular external voltage V.The air mass (AM) 1.5 global spectra are used as the irradiation to the solar cells [27].The mobility and lifetime of carriers vary widely in perovskite materials depending on the material compositions and deposition conditions [28].The electron and hole mobilities in crystalline perovskites vary from 10 to more than 100 cm 2 /Vs, whereas in polycrystalline perovskites, these are ~1-10 and ~5-20 cm 2 /Vs, respectively [29,30].The carrier lifetime usually varies in the range 0.1-100 µs in perovskite films [7].

Effects of Surface Recombination
Figures 1 and 2 show the photogenerated electron concentration profiles across a 0.5 µm-thick perovskite layer in an ETL-i-HTL type structure as a function of surface recombination velocities at both ends.The following parameters are assumed [15]: V bi = 1.2 V, µ p = 10 cm 2 V −1 s −1 , µ n = 5 cm 2 V −1 s −1 , τ n = τ p = 10 µs, and F = 1 V/µm.The absorption depth (i.e., 1/α) for the entire sun spectra is smaller than 0.2 µm.As shown in Figures 1 and 2, the carrier profile is more dependent on the surface recombination velocity of (Figure 2a) the carrier type that moves towards the bottom electrode (holes in this case).Again, the surface recombination velocity at the top contact for this carrier type critically controls the carrier profile (Figure 1a) and, thus, the overall photocurrent as well.In general, the photogenerated carrier concentration decreases with an increase in surface recombination velocity.The surface recombination velocity increases with increasing trapped carrier concentration at the interface between the perovskite and ETL or HTL layer.The interface properties between the ETL and perovskite have a more significant effect on the carrier profile.
Figures 1 and 2 show the photogenerated electron concentration profiles across a 0.5 µm-thick perovskite layer in an ETL-i-HTL type structure as a function of surface recombination velocities at both ends.The following parameters are assumed [15]: Vbi = 1.2 V, µp = 10 cm 2 V −1 s −1 , µn = 5 cm 2 V −1 s −1 , τn = τp = 10 µs, and F = 1 V/µm.The absorption depth (i.e., 1/α) for the entire sun spectra is smaller than 0.2 µm.As shown in Figures 1 and 2, the carrier profile is more dependent on the surface recombination velocity of (Figure 2a) the carrier type that moves towards the bottom electrode (holes in this case).Again, the surface recombination velocity at the top contact for this carrier type critically controls the carrier profile (Figure 1a) and, thus, the overall photocurrent as well.In general, the photogenerated carrier concentration decreases with an increase in surface recombination velocity.The surface recombination velocity increases with increasing trapped carrier concentration at the interface between the perovskite and ETL or HTL layer.The interface properties between the ETL and perovskite have a more significant effect on the carrier profile.Among four transport parameters (i.e, drift mobilities and carrier lifetimes of holes and electrons), the drift mobility for the carriers that drift towards the back electrode has the most significant effect on the overall charge collection efficiency.Considering that holes are drifting towards the bottom electrode, the effect of hole drift mobility on the hole conduction current density at V = 0 is shown in Figure 3.The surface recombination velocities are S pb = S pt = S nt = S nb = 10 4 cms −1 .The effects of other transport parameters of the bulk (i.e., electron drift mobility, electron lifetime, and hole lifetime) are almost negligible for this very thin (less than 1 µm) perovskite layer.Note that the effects of the hole and electron transport properties will be interchanged in the HTL-i-ETL type structure.Among four transport parameters (i.e, drift mobilities and carrier lifetimes of holes and electrons), the drift mobility for the carriers that drift towards the back electrode has the most significant effect on the overall charge collection efficiency.Considering that holes are drifting towards the bottom electrode, the effect of hole drift mobility on the hole conduction current density at V = 0 is shown in Figure 3.The surface recombination velocities are Spb = Spt = Snt = Snb = 10 4 cms −1 .The effects of other transport parameters of the bulk (i.e., electron drift mobility, electron lifetime, and hole lifetime) are almost negligible for this very thin (less than 1 µm) perovskite layer.Note that the effects of the hole and electron transport properties will be interchanged in the HTL-i-ETL type structure.As mentioned earlier, the surface recombination velocity at the top contact for the carriers that drift towards the back electrode critically controls the carrier profile and, hence, the overall photocurrent.Therefore, it is instructive to show this effect on the overall conversion efficiency.Figure 4 shows the effect of the surface recombination velocity As mentioned earlier, the surface recombination velocity at the top contact for the carriers that drift towards the back electrode critically controls the carrier profile and, hence, the overall photocurrent.Therefore, it is instructive to show this effect on the overall conversion efficiency.Figure 4 shows the effect of the surface recombination velocity at the top surface (S pt ) on the overall PCE of the same cell as in Figures 1 and 2. The other surface recombination velocities are S pb = S nt = S nb = 2 × 10 4 cms −1 .The carrier transport parameters in Figure 4 are the same as in Figure 1, with R s = 3 Ωcm 2 , R p = 10 4 Ωcm 2 , and J c = 3 × 10 −9 mAcm −2 .As shown in Figure 4, the surface recombination velocity at the top surface (S pt ) has a strong effect on the PCE.As mentioned earlier, the surface recombination velocity at the top contact for the carriers that drift towards the back electrode critically controls the carrier profile and, hence, the overall photocurrent.Therefore, it is instructive to show this effect on the overall conversion efficiency.Figure 4 shows the effect of the surface recombination velocity at the top surface (Spt) on the overall PCE of the same cell as in Figures 1 and 2. The other surface recombination velocities are Spb = Snt = Snb = 2 × 10 4 cms −1 .The carrier transport parameters in Figure 4 are the same as in Figure 1, with Rs = 3 Ωcm 2 , Rp = 10 4 Ωcm 2 , and Jc = 3 × 10 −9 mAcm −2 .As shown in Figure 4, the surface recombination velocity at the top surface (Spt) has a strong effect on the PCE.

Model Fitting with Measurement Data
The theoretical model fitting with the experimental results of a Cs 0.05 (FA 0.83 MA 0.17 ) 0.95 Pb(I 0.83 Br 0.17 ) 3 -based solar cell on the FTO substrate is shown in Figure 5.The material composition of the structure is FTO/TiO 2 (TiO 2 /APTES/PASCA-Br)/Cs 0.05 (FA 0.83 MA 0.17 ) 0.95 Pb(I 0.83 Br 0.17 ) 3 /PCBM/Ag [14], which is the ETL-i-HTL type.The perovskite layer thickness is 0.4 µm.The symbols and solid lines represent the experimental data and the model fittings to the experimental results, respectively.The three structures have different interfacial materials at the ETL/perovskite interface.The transport properties of electrons and holes are not equal in perovskite materials.Since these values are not precisely known, they are assumed equal in this experimental fitting in order to reduce the fitting parameters.The mathematical model considers separate values for electron and hole transport parameters, and thus, one can easily change these values in the calculation if these values are known.The transport parameters are assumed as µ p = 10 cm 2 V −1 s −1 , µ n = 5 cm 2 V −1 s −1 , and τ n = τ p = 1 µs, which are typical values for perovskites.Since the charge collection efficiency mostly depends on the carriers that move towards the bottom electrode (holes in this case), the assumed values are accurate for this carrier type.The transport parameters for other types of carriers have little effect on the overall performance, and these mentioned values may not be very accurate.Again, to reduce the fitting parameters, the surface recombination velocity for electrons is assumed to be the same as that for holes.The other fitting parameters are given in Table 1.The change in the ETL layer changes the interface properties of the ETL and perovskite layer and so does the surface recombination velocity, S pt .The reduction in S pt and R s improves the PCE.
overall performance, and these mentioned values may not be very accurate.Again, to reduce the fitting parameters, the surface recombination velocity for electrons is assumed to be the same as that for holes.The other fitting parameters are given in Table 1.The change in the ETL layer changes the interface properties of the ETL and perovskite layer and so does the surface recombination velocity, Spt.The reduction in Spt and Rs improves the PCE.The theoretical fitting with the experimental J-V characteristics of a (FAPbI3)0.95(MAP-bBr3)0.05-basedsolar cell is shown in Figure 6.The structure is the ETL-i-HTL type, and its material composition is FTO/SnO2/(FAPbI3)0.95(MAPbBr3)0.05/SpiroOMeTAD/Au.The symbols and solid lines represent the experimental data (extracted from Figure 16 of ref. [31]) and the model fittings to the experimental results, respectively.The structure mentioned above is the "control" structure, and the other structure has a 2D passivation of (BA)2PbI4 between the perovskite and HTL layer by the SIG processing method.The 2D passivation improves the interface properties and reduces the surface recombination velocity, as evident from Table 1, which leads to a higher open circuit voltage and better PCE.    1.

PCE (%)
Cs 0.05 (FA 0.83 MA 0.17 ) 0.95 Pb(I 0.83 Br 0.17   The theoretical fitting with the experimental J-V characteristics of a (FAPbI 3 ) 0.95 (MAPbBr 3 ) 0.05 -based solar cell is shown in Figure 6.The structure is the ETL-i-HTL type, and its material composition is FTO/SnO 2 /(FAPbI 3 ) 0.95 (MAPbBr 3 ) 0.05 /SpiroOMeTAD/Au.The symbols and solid lines represent the experimental data (extracted from Figure 16 of ref. [31]) and the model fittings to the experimental results, respectively.The structure mentioned above is the "control" structure, and the other structure has a 2D passivation of (BA) 2 PbI 4 between the perovskite and HTL layer by the SIG processing method.The 2D passivation improves the interface properties and reduces the surface recombination velocity, as evident from Table 1, which leads to a higher open circuit voltage and better PCE.

Conclusions
A physics-based mathematical model for describing the current-voltage characteristics in perovskite solar cells has been developed by incorporating carrier drift, diffusion, and recombination in the bulk and interfaces.The individual effects of surface recombination at both ends on steady-state electron and hole profiles across the perovskite layer are studied.The theoretical calculations are compared with published experimental data and show excellent agreement.This theoretical model can serve as a tool for extracting charge carrier transport and cell parameters.The mathematical model can also be applied to other solar cell technologies if the cell structure is similar to PSC (for example, the structure of a bulk heterojunction organic solar cell is similar to PSC).
Funding: The author acknowledges the financial support from NSERC through its Discovery Grant program.

Figure 1 .Figure 1 .
Figure 1.The profile of photogenerated electron concentration across the perovskite layer in an ETLi-HTL type structure for various recombination velocities at both ends.

Figure 1 .Figure 2 .
Figure 1.The profile of photogenerated electron concentration across the perovskite layer in an ETLi-HTL type structure for various recombination velocities at both ends.

Figure 2 .
Figure 2. The profile of photogenerated hole concentration across the perovskite layer in an ETL-i-HTL type structure for various recombination velocities at both ends.

Figure 3 .
Figure 3. Hole conduction current densities versus distance from the radiation-receiving electrode for different hole mobilities.Hole lifetime, τ p = 10 µs.

Figure 4 .
Figure 4. PCE versus S pt of a MAPBI 3 -based ETL-i-HTL type solar cell.

Figure 5 .
Figure 5. Current−voltage characteristics of a Cs 0.05 (FA 0.83 MA 0.17 ) 0.95 Pb(I 0.83 Br 0.17) 3 -based solar cell on FTO substrate.Symbols represent extracted experimental data from Figure 4 of [14].The legends represent compositions of the ETL layer.
Figure 5. Current−voltage characteristics of a Cs 0.05 (FA 0.83 MA 0.17 ) 0.95 Pb(I 0.83 Br 0.17) 3 -based solar cell on FTO substrate.Symbols represent extracted experimental data from Figure 4 of [14].The legends represent compositions of the ETL layer.

Table 1 .
The extracted parameters for Figures5 and 6are summarized in Table

Table 1 .
The extracted parameters for Figures5 and 6are summarized in Table1.