Mustafa
Kareem
*ab,
Ethar Yahya
Salih
c,
Malatesh
Akkur
d,
Satish Kumar
Samal
e,
Sridharan
Sundharam
f and
Sanjeev
Kumar
g
aCollege of Remote Sensing and Geophysics, Al-Karkh University of Science, Haifa St., Baghdad 10011, Iraq. E-mail: dr.mustafa@kus.edu.iq
bCollege of Science, University of Warith Al-Anbiyaa, 56001 Karbala, Iraq
cCollege of Energy and Environmental Sciences, Al-Karkh University of Science, Baghdad 10081, Iraq
dDepartment of Physics & Electronics, School of Sciences, JAIN (Deemed to be University), Bangalore, Karnataka, India
eDepartment of Electronics & Communication Engineering, Siksha ‘O’ Anusandhan (Deemed to be University), Bhubaneswar, Odisha-751030, India
fDepartment of Physics, Sathyabama Institute of Science and Technology, Chennai, Tamil Nadu, India
gDepartment of Physics, University Institute of Sciences, Chandigarh University, Mohali, Punjab, India
First published on 10th December 2025
We report high-performance organometallic perovskite solar cells (PSCs) utilizing a thermodynamically stable hematite (α-Fe2O3) as an electron transport layer (ETL). The incorporation of an α-Fe2O3 layer could provide suitable energy band alignment with the perovskite, as well as induce a deep valence band maximum, which promotes electron extraction from the conduction band of the perovskite and facilitates hole blocking. Using the solar cell capacitance simulator (SCAPS-1D) software for the optimized PSC under the typical AM 1.5G light spectrum, it was predicted to achieve a competitive power conversion efficiency (PCE) of 25.62% with a high short circuit current (JSC) of 23.58 mA cm−2, an open circuit voltage (VOC) of 1.286 V, and a fill factor (FF) of 84.39%. Moreover, high thermal stability of PSCs with exposure to a high temperature of 85 °C can be attained. Through a series of optimization processes, we conclude that a thinner α-Fe2O3 layer (10 nm) improves charge extraction and the transmittance of the visible light, while the decreased defect density significantly reduces recombination rates, thereby enhancing VOC and PCE. An optimum perovskite thickness of 800 nm was found to maximize light absorption. Additionally, controlled acceptor doping concentration (NA) enhanced carrier extraction and quasi-Fermi level splitting (QFLS), while high series resistance (RS) and low shunt resistance (RSH) were demonstrated to limit FF and efficiency.
Several n-type metal oxides, such as titania (TiO2), zinc, niobium, cerium, and tin, have been explored as ETLs in planar n–i–p PSCs.5,10 The ETL used most often is TiO2, and solar cells fabricated with this material have shown a PCE approaching 26%.11 However, in addition to its defect-rich surface, the TiO2 layer degrades PSCs when exposed to ultraviolet (UV) light, which impairs the operational stability of PSCs. Also, the fabrication of a TiO2 layer demands a high annealing temperature beyond 450 °C, which hampers its use in mass manufacturing.12 As a thermodynamically stable n-type iron oxide semiconductor, hematite (α-Fe2O3) has attracted much attention in optoelectronics owing to its low cost, abundance, non-toxicity, high chemical stability, and desirable energy band position.13 Fe2O3 has been utilized as a UV-stable photoelectrode in PSCs due to its UV photocatalytic resistance.14 Hu et al. applied α-Fe2O3 as the ETL in methylammonium lead triiodide (MAPbI3)-based PSCs and achieved an efficiency of 10.78% with enhanced durability as compared to a TiO2-based device. This enhancement is due to a high built-in potential (Vbi) across MAPbI3 in α-Fe2O3 PSCs, which results in improved charge transport and lower accumulation at the interface.15 Later, Qureshi et al. reported different solvents to improve the crystallinity and morphology of α-Fe2O3 ETLs. Among all, the ethanol-processed α-Fe2O3 layer leads to a champion efficiency of 13% with a suppressed hysteresis index of 0.04 in an n–i–p device.16 Recently, Bilal et al. demonstrated that dual interfacial modification of α-Fe2O3 enhances the efficiency and stability of planar PSCs. The 4-trifluorophenylammonium iodide dual modifier resulted in a performance improvement, yielding a champion efficiency of 15.63%.17
The present work reports for the first time, through systematic SCAPS-1D simulations, the potential of hematite-based ETLs for designing high-performance perovskite photovoltaics. Due to its deeper conduction band minimum than that of the perovskite, α-Fe2O3 can accelerate electron extraction from the perovskite absorber. We have tuned the layer thickness, trap-state density (NTrap), defect density, doping concentration of perovskite, resistive losses, and operating temperature. Numerically optimized α-Fe2O3-based PSCs let us achieve a PCE of 25.62%. Simultaneously, the optimized PSC exhibits high thermal stability, maintaining 93% of its initial performance at 85 °C and under AM 1.5G solar irradiation. The suggested approach to ETL optimization may represent a significant advancement for reliable and effective PSCs.
| Rrad = B(np − ni2) | (1) |
| RAuger = Cnn(np − ni2) + Cpp(np − ni2) | (2) |
![]() | (3) |
The Eg and electron affinity of the α-Fe2O3 film were fixed at 2.13 eV and 3.96 eV, respectively. α-Fe2O3 was prepared by the high pressure hydrothermal method, as reported by Khurram et al.21 The carrier-charge mobilities and dielectric constant were adopted from density functional theory-based studies,22,23 describing intrinsic properties of α-Fe2O3. A more accurate PSC can be developed by incorporating ETL/MAPbI3 and MAPbI3/HTL interfaces. Both interfaces are assumed to have a neutral defect density of 1010 cm−2 and a characteristic energy of 0.1 eV (see Table S1). In the initial stage, RS and RSH are considered to be 4 Ohm cm2 and 300 Ohm cm2, which are close to the experimentally obtained values. The key parameters are taken from previously published experimental and theoretical reports, as summarized in Table 1. The standard AM 1.5G illumination and a temperature of 300 K were utilized for the simulations.
| Parameters (units) | α-Fe2O3 | MAPbI3 | Spiro-OMeTAD | TiO2 | SnO2 | ZnO |
|---|---|---|---|---|---|---|
| Thickness [nm] | 40 | 550 | 100 | 10 | 10 | 10 |
| Bandgap Eg [eV] | 2.13 | 1.59 | 3.01 | 3.2 | 4.1 | 3.3 |
| Electron affinity [eV] | 3.96 | 3.9 | 2.18 | 4.05 | 4.39 | 4.17 |
| Dielectric permittivity | 26 | 25 | 3 | 33 | 7.25 | 3.7 |
| Effective density of states in the conduction band [cm−3] | 1.0 × 1017 | 2.0 × 1019 | 2.2 × 1018 | 1.0 × 1017 | 1.0 × 1017 | 1.0 × 1017 |
| Effective density of states in the valence band [cm−3] | 1.8 × 1018 | 1.0 × 1019 | 2.0 × 1018 | 1.8 × 1018 | 1.8 × 1018 | 1.8 × 1018 |
| Electron mobility [cm2 V−1 s−1] | 0.098 | 25 | 7.9 × 10−3 | 18.5 | 265 | 8.4 × 10−1 |
| Hole mobility [cm2 V−1 s−1] | 0.031 | 105 | 7.9 × 10−3 | 1.0 × 10−4 | 7.6 | 91.9 |
| Electron thermal velocity [cm s−1] | 1.0 × 107 | 1.0 × 107 | 1.0 × 107 | 1.0 × 107 | 1.0 × 107 | 1.0 × 107 |
| Hole thermal velocity [cm s−1] | 1.0 × 107 | 1.0 × 107 | 1.0 × 107 | 1.0 × 107 | 1.0 × 107 | 1.0 × 107 |
| Shallow donor density [cm−3] | 1.0 × 1018 | — | — | 1.0 × 1018 | 1.0 × 1018 | 1.0 × 1018 |
| Shallow acceptor density [cm−3] | — | 1.0 × 1016 | 1.2 × 1017 | — | — | — |
| Trap-state density [cm−3] | 1.0 × 1015 | 1.03 × 1016 | 1.0 × 1015 | 1.0 × 1015 | 1.0 × 1015 | 1.0 × 1015 |
| Electron capture cross-section [cm2] | 1.0 × 10−15 | 1.0 × 10−15 | 1.0 × 10−15 | 1.0 × 10−15 | 1.0 × 10−15 | 1.0 × 10−15 |
| Hole capture cross-section [cm2] | 1.0 × 10−15 | 1.0 × 10−15 | 1.0 × 10−15 | 1.0 × 10−15 | 1.0 × 10−15 | 1.0 × 10−15 |
| Defect position above the valence band edge [eV] | 0.6 | 0.6 | 0.6 | 0.6 | 0.6 | 0.6 |
| Energetic distribution | Single | Single | Single | Single | Single | Single |
| Ref. | 21–23 | 24–26 | 27–29 | 30–32 | 33–35 | 36 and 37 |
![]() | ||
| Fig. 1 Schematic illustration and calibration of the designed PSC structure. (a) Layered structure of the PSC. (b) Band diagram of the simulated device. (c) Electric field distribution. Inset: energy level alignment. (d) J–V plots of the fabricated and simulated cells under AM1.5 solar light.24 (e) IPCE spectra. (f) Nyquist plots extracted from AC impedance analysis. (g) Absorption coefficient spectra. (h) C–V and Mott–Schottky plots under dark conditions. (i) Recombination and generation rate profiles. | ||
Fig. 2a presents the J–V curves of PSCs with varying α-Fe2O3 ETL thicknesses. Fig. 2b shows that VOC and JSC reduce as the α-Fe2O3 thickness increases from 10 nm to 50 nm, indicating limited charge collection performance in thicker ETLs. A thinner α-Fe2O3 ETL could facilitate enhanced transmittance of the visible light, which contributes to more photon absorption of the MAPbI3 layer. The variations in FF and PCE are exhibited in Fig. 2c, implying that an optimum performance is obtained with thinner ETLs (10–20 nm); instead, thicker layers cause transfer barriers that deteriorate photovoltaic properties. The findings are further supported by the IPCE spectra in Fig. 2d, which demonstrate that PSCs with thinner α-Fe2O3 layers exhibit enhanced quantum efficiency across the visible spectrum. A thick ETL absorbs/scatters incoming light, thereby suppressing the photon flux reaching MAPbI3 and consequently lowering JSC. Fig. 2e shows that the device with an ETL with a thickness of 10 nm has a stronger electric field at the ETL/perovskite interface than the 50-nm ETL. This allows for separate carriers and minimizes recombination losses. Finally, Nyquist curves (Fig. 2f) show enlarged recombination resistance (Rrec) in the thinner α-Fe2O3, aligning well with enhanced charge extraction and mitigated recombination pathways. The 10-nm α-Fe2O3-based 2DRP device achieved a PCE of 14.47%.
Fig. S1 shows the J–V curves of 10-nm Fe2O3-based devices with varying defect densities in Fe2O3 ETLs. Using different defect concentrations in the 10-nm Fe2O3 layer induced ignorable changes in device performance due to the ultrathin ETLs, where charge-carrier transfer occurs by drift rather than bulk diffusion. Consequently, the light harvesting efficiency of PSCs is unchanged with increasing defects in ETLs (Fig. S2). Furthermore, the donor doping concentration (ND) in the 10-nm Fe2O3 layer was varied from 1014 cm−3 to 1018 cm−3, as shown in Fig. S3 and S4. Increasing the ND level led to a slight enhancement in efficiency (0.55%) owing to suppression in series resistance inside Fe2O3, implying that the ETLs perform as electron-selective rather than recombination sites.
The effect of perovskite NTrap on the PSC performance was systematically assessed while changing the values from 1015 cm−3 to 1017 cm−3, as shown in Fig. 3a. It is found that the device's JSC does not change while varying NTrap, while VOC gradually reduces to 0.91 V (Fig. 3b). High traps force photogenerated carriers to recombine before extraction, which narrows the QFLS.38 The IPCE responses (Fig. S5) show nearly identical patterns across different NTrap values, further corroborating the unchanged trend of JSC. On the other hand, both FF and PCE are markedly decreased with increasing defect concentration, as shown in Fig. 3c. The power lose in PSCs occurs because of trap-assisted non-radiative. This process is known as SRH recombination. According to the SRH model, a reduction in NTrap leads to an extension in the charge carrier's lifetime. This, in turn, results in higher diffusion lengths and a decrease in recombination. The Nyquist patterns shown in Fig. 3d prove this behavior by demonstrating smaller semicircle diameters at higher NTrap, indicating lower Rrec and higher carrier losses. As shown in Fig. 3e, when NTrap reduces, the recombination rate inside MAPbI3 and at the interfaces decreases. Additionally, the electric field alignment across the PSC (Fig. 3f) slightly improves by decreasing NTrap, providing better charge extraction. From the band diagram structure (Fig. 3g and h), the QFLS variations with increasing NTrap were calculated and are plotted in Fig. 3i. Increasing NTrap led to noticeable declines in QFLS magnitude due to increased nonradiative recombination losses. The same pattern is observed in the calculated VOC, indicating voltage losses.
The device performance using the FTO/α-Fe2O3/MAPbI3/spiro-OMeTAD/Au configuration is analyzed with perovskite thicknesses ranging from 400 nm to 800 nm in increments of 100, and the corresponding J–V plots are illustrated in Fig. 4a. According to Fig. 4b and c, the JSC and PCE parameters notably improve as the MAPbI3 thickness increases, while VOC remains almost constant, and FF shows a minor improvement. The enhancement in FF is due to improved carrier collection performance, which decreases resistive losses and mitigates recombination processes. The highest JSC (23.61 mA cm−2) and PCE (17.36%) are realized at a thickness value of 800 nm owing to improved photon absorption.39 The IPCE plots indicate that thicker perovskites show a higher spectral response for visible light, consistent with improved absorption and JSC (Fig. 4d). Collectively, the impedance analysis shown in Fig. 4e and the recombination rate profiles shown in Fig. 4f exhibit limited changes in Rrec and recombination dynamics, particularly near interfaces.
Fig. 5a shows the J–V performances of α-Fe2O3-based PSCs at different perovskite doping densities of shallow acceptor (NA), ranging from 1015 to 1019 cm−3. As displayed in Fig. 5b and c, VOC and FF showed an improved pattern with increasing NA concentration, leading to a higher PCE value of 21.17% at 1019 cm−3. It was found that an increase in NA concentration enhances device performance by increasing the Vbi and depletion width at the Fe2O3/MAPbI3 junction, thereby promoting faster charge-carrier separation.40 However, the JSC revealed a slight reduction due to the carrier scattering effect. The bulk recombination in perovskites is significantly minimized at high NA values due to suppressing trap-mediated SRH processes, as confirmed in Fig. 5d. The electric field distribution in Fig. 5e exhibits an intensified field at a high NA concentration, facilitating more efficient carrier collection. Also, high NA doping can improve band bending and result in an increased QFLS, as shown in the energy band diagram (Fig. 5f), which matches well with the observed VOC increase.
Individual RSH and RS resistances significantly affect the efficiency of the solar cells, as they control the slopes and shapes of the J–V features, which are produced by stacking various materials, metal contacts, the semiconductor–metal interface, and an inadequate PSC fabrication approach. To evaluate the impact of RS on the PSC, the value of RS is altered from 2.0 to 10 Ohm cm2, while the value of RSH is kept constant at 300 Ohm cm2, as depicted in Fig. 6a. It can be seen that VOC and JSC (Fig. 6b) are almost independent of variation in RS, while FF severely drops, leading to a reduction in PCE to 18.67%, with an increase in RS to 10 Ohm cm2 (Fig. 6c), which stems from resistive power losses and impeded charge extraction. After that, the value of RSH was varied from 500 Ohm cm2 to 4000 Ohm cm2 but the RS value was kept fixed at 2 Ohm cm2 and the changes in J–V curves are exhibited in Fig. 6d. The results indicate that the photovoltaic parameters tend to improve with the increase in RSH until 3000 Ohm cm2 (Fig. 6e and f), and after that, they slightly change with further increase in RSH. These findings corroborate well with those reported in ref. 41. According to the modified diode equation,41 a low RSH induces leakage current through recombination routes for photogenerated carriers, leading to decreased device performance, whereas higher RSH values suppress unwanted leakage currents, as supported by increasing Rrec (Fig. S6). Therefore, the device achieved an optimized PCE of 25.56% with a JSC of 23.60 mA cm−2, a VOC of 1.28 V, and an FF of 84.30% at an RS of 2 Ohm cm2 and an RSH of 4000 Ohm cm2.
To observe the influence of operating temperature on the solar cells, the value of temperature is changed from 25 to 85 °C, depicted in Fig. 7a. Increasing the temperature does not have an impact on JSC, indicating that light harvesting and carrier generation are not affected. It is found that VOC and FF are progressively reduced with increase in the temperature of PSCs, attributed to the induced reverse saturation current in the PSCs. Interestingly, an increase in the temperature from 25 to 85 °C resulted in a minor reduction in the overall PCE (Fig. 7b), where the PSC retained about 93% of its initial performance, demonstrating its high thermal stability. An increase in temperature creates deformation stress, resulting in a forward movement that increases NTrap at the interfaces. This influences the diffusion length, thereby limiting device performance. To address this hypothesis, we performed temperature-dependent interface defect-density simulations at low (1010 cm−2) and high (1015 cm−2) defect concentrations of Fe2O3/MAPbI3 and MAPbI3/spiro-OMeTAD interfaces. The simulated J–V plots of PSCs with different temperatures are shown in Fig. S7, and the corresponding efficiencies are shown in Fig. S8. Increasing the temperature results in a reduction in VOC for both low- and high-defective devices due to increased charge recombination. A PSC with 1015 cm−2 concentration reveals high drops in PCE under both temperature conditions. Therefore, interfacial defect passivation is important to achieve efficient and stable PSCs. Without modification, defective interfaces hinder carrier transfer, leading to energy loss and low power output.17 The decreasing interfacial defects are demonstrated to enhance the IPCE spectrum (Fig. S9), suggesting that the modification can increase charge carrier generation and minimize the recombination rate (Fig. S10).
Then, the performance of Fe2O3-based PSCs was compared with that of other typical semiconductors at the same optimized MAPbI3 and spiro-OMeTAD layers. Fig. 7c shows the J–V characteristics of different devices simulated at room temperature. All calculations were conducted with optimized series and shunt resistances and an ETL thickness of 10 nm. As depicted in Fig. 7d, all devices yielded close values of JSC because the perovskite layer mainly governs the light absorption process. A device with the Fe2O3 ETL achieved a VOC value higher than that of a tin oxide (SnO2) device, and the value closely matched that of titanium dioxide (TiO2) and zinc oxide (ZnO) devices. In addition, the Fe2O3-based PSC demonstrates strong FF (Fig. 7e), thus surpassing the performance of PSCs with ZnO and SnO2 ETLs due to better conduction band alignment and suppressed recombination losses. A comparison of our results with those of experimentally reported α-Fe2O3-based PSCs (Fig. 7f) demonstrates that this study attains a markedly superior PCE (25.62%), exceeding prior reports from 2017 to 2025 (Table 2), thereby illustrating the benefits of the optimized device design, even under thermal stress conditions.
| Structure | V OC (V) | J SC (mA cm−2) | FF (%) | PCE (%) | Year | Ref. |
|---|---|---|---|---|---|---|
| FTO/Ni–Fe2O3/MAPbI3/spiro-OMeTAD/Au | 0.92 | 22.35 | 69.10 | 14.20 | 2017 | 42 |
| FTO/Fe2O3/MAPbI3/spiro-OMeTAD/Au | 0.65 | 20.4 | 16.58 | 10.78 | 2017 | 15 |
| FTO/Fe2O3/MAPbI3/spiro-OMeTAD/Au | 1.55 | 11.27 | 33.00 | 5.70 | 2018 | 43 |
| FTO/Fe2O3/PCBM/MAPbI3/spiro-OMeTAD/Au | 0.91 | 23.2 | 68.00 | 14.2 | 2018 | 44 |
| FTO/Ti–Fe2O3/MAPbI3/spiro-OMeTAD/Ag | 1.10 | 21.56 | 75.00 | 17.85 | 2018 | 45 |
| FTO/Fe2O3/PbI2/Cs0.05(FA0.83MA0.17)0.95Pb(I0.83Br0.17)3/spiro-OMeTAD/Au | 1.06 | 23.48 | 77.10 | 19.32 | 2022 | 46 |
| ITO/Fe2O3/(FAPbI3)0.97(MAPbBr3)0.03/spiro-OMeTAD/Au | 0.99 | 24.79 | 58.31 | 14.33 | 2022 | 47 |
| ITO/Fe2O3/FA1−xMAxPb(I1−yBry)3/spiro-OMeTAD/Au | 0.98 | 23.45 | 54.74 | 12.61 | 2022 | 14 |
| FTO/Fe2O3/MAPbI3/spiro-OMeTAD/Au | 0.70 | 21.10 | 56.70 | 8.60 | 2023 | 48 |
| ITO/Fe2O3/Cs0.05(FA0.83MA0.17)0.95Pb(I0.83Br0.17)3/spiro-OMeTAD/WO3/Au | 1.05 | 19.21 | 64.48 | 13.01 | 2023 | 16 |
| FTO/TiO2/Fe2O3@SnO2/CsPbBr3/carbon | 1.60 | 7.88 | 80.85 | 10.23 | 2023 | 49 |
| ITO/Fe2O3/SnO2/MAPbI3/spiro-OMeTAD/Au | 1.14 | 21.29 | 75.13 | 18.24 | 2024 | 50 |
| ITO/Fe2O3/CF3PhAI/Cs0.05(FA0.83MA0.17)0.95Pb(I0.83Br0.17)3/CF3PhAI/spiro-OMeTAD/Au | 1.12 | 19.01 | 72.90 | 15.63 | 2025 | 17 |
| FTO/Fe 2 O 3 /MAPbI 3 /spiro-OMeTAD/Au | 1.28 | 23.58 | 84.39 | 25.62 | 2025 | This work |
The data will be available from the corresponding author on reasonable request.
| This journal is © the Owner Societies 2026 |