Open Access Article
Mustafa Kareem
*ab,
Ethar Yahya Salihc,
Sampangi Rama Reddy B. Rd,
Praveen Priyaranjan Nayake,
Sridharan Sundharamf and
Sanjeev Kumarg
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), Bengaluru, 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 12th March 2026
Two-dimensional Ruddlesden–Popper (2DRP) halide perovskites have attracted considerable interest in photovoltaics due to their optoelectronic properties and superior environmental stability. Herein, a 2-(methylthio)ethylamine (MTEA+) bulky spacer was introduced as the organic interlayer with sulfur–sulfur bonding, which interacted electrostatically with the inorganic framework to induce the oriented (MTEA)2(MA)4Pb5I16 2DRP perovskite and lattice stabilization. We demonstrated, through numerical simulations, that the performance of 2DRP-based perovskite solar cells (PSCs) could be improved by adapting the perovskite film thickness, trap-state density, parasitic resistances, and temperature. Optical analyses revealed that the 2DRP perovskite dominated light-harvesting in the visible spectrum, while charge transport materials remained largely transparent. PSCs with 2DRP (n = 5) showed a maximum power conversion efficiency (PCE) of 20.79% with an outstanding open circuit voltage (VOC) of 1.49 V under standard AM1.5G illumination. Moreover, the proposed simulation predicted potentially improved thermal stability for the optimized device, theoretically retaining 95% of its initial performance. A simulated PCE of 20.79% represented a theoretical upper limit achievable only under optimized charge dynamic conditions.
Metal-halide perovskites have not only transformed photovoltaics but also emerged as a multifunctional material for a wide range of optoelectronics, including light-emitting diodes, photodetectors, and lasers, due to their high absorption coefficient, long charge-carrier diffusion length, tunable bandgap, facile solution processability, and defect tolerance.4 Perovskite solar cells (PSCs) have an important role in photovoltaic technologies because of their low cost, superior optoelectronic characteristics, and steadily increasing power conversion efficiency (PCE).5 The PCE of three-dimensional (3D) halide PSCs has improved from 3.8% to 27% in recent years, which is comparable with that of silicon photovoltaics.6 Nevertheless, the instability of 3D halide perovskites is a major hurdle that restricts the commercial progress of perovskite photovoltaics.7 Conversely, 2D-layered perovskites with atomically fine quantum well frameworks show higher atmospheric stability, in addition to the quantum confinement effect, tunable composition and structure homogeneity and have resulted in beneficial uses in optoelectronics.8
2D Ruddlesden–Popper (2DRP) perovskites have garnered significant interest for utilization in PSCs thanks to their enhanced moisture tolerance, improved thermal stability, low self-doping, and suppressed ion migration compared with their conventional 3D counterparts.9 2DRP perovskites can be identified by the fundamental chemical formula A2Bn−1MnY3n+1, where Y stands for a halide, M for a divalent metal, and A and B for monovalent cations. 2D perovskites are composed of one or more inorganic layers encased in bulky organic spacers. Strong hydrogen and Coulombic bonds tightly bind spacer ions, greatly enhancing structural stability.10 However, due to their high exciton binding energies and wide bandgap energies, 2DRP perovskites have shown lower efficiencies for photovoltaics than those reported for 3D devices, despite their enhanced stability.11 Recently, many bulky ammonium spacers, including allylammonium,12 2-thiophenemethylammonium,13 phenylethylamine,14 cyclohexylammonium,15 and polyethylenimine,16 have been studied to progressively improve the PCE of 2DRP PSCs. Ren and co-workers reported that the 2-(methylthio)ethylamine (MTEA+) spacer can operate as an efficient bulky organic cation in 2DRP perovskites. The MTEA+ spacer enables sulfur–sulfur (S–S) interactions, besides van der Waals forces, providing the perovskite lattice with improved charge-carrier dynamics and stabilization.10 Later, Yin and co-workers introduced the MTEA+ spacer into a tin-based 2DRP perovskite, which benefited from S–S non-covalent interactions. The authors demonstrated that these interactions induced the oriented crystallization of 2DRP perovskite, improved the film morphology, and reinforced the perovskite framework.17
However, most studies have focused on experimental evaluations of certain spacer cations, offering minimal theoretical understanding of the impact of device parameters and defect dynamics on photovoltaic performance. We wished to identify performance limits and dominant loss mechanisms for 2DRP PSCs. The reported high efficiency, therefore, represents an upper-bound scenario. We explored the optoelectronic characteristics of the lead-based 2DRP (MTEA)2(MA)4Pb5I16 (n = 5) perovskite with MTEA+ interlayer spacers. Using 1D solar cell capacitance simulator (SCAPS-1D) calculations, we characterized the effects of perovskite layer thickness, trap-state density (NTrap), and shallow acceptor density (NA). After optimization, the 2DRP device achieved a maximum PCE of 18% with a high VOC of 1.48 V. Then, we combined mobility-thickness-defect-density co-optimization to optimize 2DRP perovskite beyond 20% efficiency. Additionally, the MTEA+-based 2DRP PSCs retained 95% of their initial performance under 420 K in the SCAPS-1D simulation, indicating a theoretically enhanced thermal tolerance.
The standard AM1.5G spectral illumination with an intensity of 1000 W m−2 at 300 K was utilized as input solar radiation for the calculation. We used fluorine-doped tin oxide (FTO) as an illuminated side with a flat-band contact and a gold (Au) electrode on an opaque side with a work function of 5.1 eV. SCAPS-1D simulations closely aligned with experimental findings when critical conditions such as parasitic resistances and interfacial and bulk defects were considered. We set the series resistance (RS) and shunt resistance (RSH) at 4 Ohm cm2 and 1000 Ohm cm2, respectively. Two interfacial layers with a trap-state density (NTrap) of 1010 cm−2 were incorporated at the electron transport material (ETM)/perovskite and perovskite/hole transport material (HTM) interfaces (Table S1). The fundamental ETM, HTM, and absorber characteristics are listed in Table 1 and are derived from some earlier studies. The energetic distribution was considered as a single discrete level, rather than a Gaussian or exponential distribution, to simplify the recombination modeling while capturing the dominant trap-assisted recombination pathway.
| Parameters (units) | FTO | TiO2 | (MTEA)2(MA)4Pb5I16 | Spiro-OMeTAD |
|---|---|---|---|---|
| Thickness (nm) | 300 | 50 | 350 | 100 |
| Bandgap Eg (eV) | 3.5 | 3.2 | 1.8 | 3.01 |
| Electron affinity (eV) | 4.0 | 4.05 | 3.6 | 2.18 |
| Dielectric permittivity | 9 | 33 | 12 | 3 |
| Effective density of states in conduction band (cm−3) | 2.2 × 1018 | 1.0 × 1018 | 1.0 × 1019 | 2.2 × 1018 |
| Effective density of states in valence band (cm−3) | 1.8 × 1019 | 1.0 × 1019 | 1.0 × 1018 | 2.0 × 1017 |
| Electron mobility (cm2 V−1 s−1) | 20 | 18.5 | 1.3 × 10−2 | 7.9 × 10−3 |
| Hole mobility (cm2 V−1 s−1) | 10 | 1.0 × 10−4 | 1.7 × 10−2 | 7.9 × 10−3 |
| Electron thermal velocity (cm s−1) | 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 |
| Shallow donor density (cm−3) | 2.0 × 1019 | 1.0 × 1019 | — | — |
| Shallow acceptor density (cm−3) | — | — | 1.0 × 1015 | 1.2 × 1017 |
| Defect type | Neutral | Neutral | Neutral | Neutral |
| Trap-state density (cm−3) | 1.0 × 1015 | 1.0 × 1016 | 2.22 × 1015 | 1.0 × 1015 |
| Electron capture cross-section (cm2) | 1.0 × 10−16 | 1.0 × 10−16 | 1.0 × 10−16 | 1.0 × 10−16 |
| Hole capture cross-section (cm2) | 1.0 × 10−16 | 1.0 × 10−16 | 1.0 × 10−16 | 1.0 × 10−16 |
| Defect position above the valence band edge (eV) | 0.6 | 0.6 | 0.6 | 0.6 |
| Energetic distribution | Single | Single | Single | Single |
| References | 21 | 22–24 | 10 and 25 | 26–28 |
The perovskite film could have various defects, such as interstitials, vacancies, and grain boundaries, which would decrease the stability and efficiency of PSCs by acting as charge recombination sites. The NTrap within the 2DRP perovskite was changed from 1014 to 1018 cm−3, and the corresponding J–V plots are illustrated in Fig. 2a. All photovoltaic parameters were reduced progressively with increasing defects. As shown in Fig. 2b, a pronounced decline in JSC could be assigned to increased non-radiative recombination sites, which lower charge-carrier extraction. Additionally, a higher NTrap reduces carrier lifetime, shrinks the QFLS, and thereby limits VOC.32 Consequently, the PCE of devices decreased from 15.35% to 0.22% with defects varying from 1014 to 1018 cm−3 (Fig. 2c). The reduction of IPCE spectra across the visible range in Fig. 2d further confirmed the suppression in charge collection and JSC. Fig. 2e shows Nyquist plots derived from C–f measurements for PSCs at different NTrap values. The results demonstrated a narrowing in semicircles at elevated traps, implying a decrease in recombination resistance (Rrec). These findings were corroborated by recombination curves in Fig. 2f, which showed that a higher NTrap accelerated the recombination rate throughout the 2DRP absorber. According to the SRH equation,20 if NTrap is passivated, charge carrier lifetime increases and, hence, longer diffusion length and less recombination are reported. The best-performing device achieved a PCE of 15.35% at an NTrap of 1014 cm−3 and a charge-carrier diffusion length of 580 nm, which lies within the experimentally reported value.33,34 Moreover, research shows that 2D and quasi-2D perovskites usually have a lower NTrap than 3D perovskite frameworks. This is because of organic spacer cations and passivation effects, which support the idea of low trap density states.35
Increasing the perovskite thickness generally results in more photons being harvested, but thicker perovskites can be detrimental to carrier generation if carrier diffusion lengths are small. To evaluate how the PSC parameters are affected by the absorber thickness, the perovskite layer thickness was varied within a range of 200 nm to 600 nm. The obtained J–V characteristics are displayed in Fig. 3a. The device VOC was nearly constant, and JSC improved significantly from 11.52 to 16.54 mA cm−2 while increasing the 2DRP thickness from 200 to 600 nm (Fig. 3b). This current enhancement was due to an increasing photoresponse across the visible spectrum, enabling a high number of photo-induced carriers, as confirmed by IPCE curves in Fig. 3d. However, FF decreased moderately with an increase in thickness due to increasing trap-assisted recombination and the saturation current (Fig. 3c). A higher charge carrier recombination at a thick perovskite was further supported by decreasing Rrec (Fig. 3e) and increasing recombination rate (Fig. 3f). The simulated recombination Rrec had only a slight reduction (inset of Fig. 3e) because, under the fixed interface defect concentration and transport-layer properties utilized in our simulation, recombination was primarily limited by interface processes rather than bulk thickness effects. We selected the optimal thickness of 2DRP as 500 nm for further calculations. The PSC parameters were VOC = 1.47 V, JSC = 16.43 mA cm−2, FF = 67.55%, and PCE = 16.38%.
RSH and RS, which have a significant effect on solar cell performance, determine the shape and slope of J–V properties. In fact, RS originates from the electrical resistance of electrodes, as well as from the electrical dissipation that takes place in the perovskite, HTM, and ETM. Several carrier recombination pathways emerge as RSH, which is also greatly impacted by device architecture. Therefore, a low RSH could impact the collected photocurrent and cause a loss in photovoltage.36
Fig. 4a shows the J–V plots of 2DRP-based PSCs at different RS while maintaining RSH at 1000 Ohm cm2. For an RS > 8 Ohm cm2, JSC was downshifted slightly to 16.31 mA cm−2 (Fig. 4b). An increment in RS decreased the FF from 68.95% to 60.65% of PSCs. Consequently, the overall performance fell rapidly from 16.75% to only 14.55% due to ohmic losses that limited carrier extraction (Fig. 4c). Fig. 4d shows the J–V plots for various RSH ranging from 1000 to 5000 Ohm cm2. The JSC and VOC remained constant as RSH increased (Fig. 4e), but FF and efficiency gradually improved as RSH increased (Fig. 4f). A higher RSH reduces leakage pathways and minimizes parasitic recombination, leading to improved charge collection efficiency, as evidenced by the Nyquist plots in Fig. S1. At 5000 Ohm cm2, we obtained an FF = 73.52% and PSC = 18%. These trends in PSC behavior align well with the literature.37
To probe the behaviour of PSCs in an operating environment, we evaluated the performance of the PSC by changing the temperature from 300 K to 420 K. Fig. 4g depicts the J–V curves of simulated devices. Fig. 4h exhibits the tracking of normalized PCEs with temperature. VOC and PCE were slightly reduced with an increase in temperature, indicating a predication for high thermal stability. This reduction in PCE occurred because, as temperature rises, the reverse saturation current increases. Moreover, the increasing temperature can result in the thermal excitation of electrons, causing some vibrations and eventual instability, and improving charge-carrier recombination rates. Optimized devices maintained 95% of their initial efficiency at 420 K compared with 88% thermal stability for non-optimized devices (control), predicting high thermal stability. In SCAPS-1D, temperature influences carrier transport, recombination, and energy band alignment, but does not account for thermally induced degradation mechanisms such as phase decomposition, ion migration, or material volatilization. Therefore, the anticipated high thermal stability represented an upper performance limit under steady-state conditions rather than long-term operational durability.
For further optimization of 2DRP PSCs, we increased charge-carrier mobility from 10−5 to 0.1 cm2 V−1 s−1, which could be achieved experimentally by the vertical orientation of 2DRP films with heteroatom-containing spacers or employing spacers with stronger interlayer interactions. The improved mobility suppressed space-charge accumulation, thus lowering recombination rates and improving charge collection. Therefore, the simulated PSC realized a PCE of 20.79%, with a JSC of 17.17 mA cm−2, a VOC of 1.49 V, and an FF of 81.03% (Fig. S2). The substantial increase in VOC and FF indicates that improved vertical transport mitigated non-radiative losses. The inherent semiconducting characteristics of the 2DRP perovskite and relatively wide bandgap (1.8 eV), which sets a high theoretical VOC limit, allowed high VOC values in the optimized PSCs. Additionally, optimizing the defect density, band alignment, and charge mobility enabled QFLS to approach this thermodynamic limit with only a VOC deficit of 0.32 eV. Experimental results of charge transport in 2DRP perovskites have shown an effective carrier mobility in the range of 8.0 ± 0.6 cm2 V−1 s−1 for layered frameworks such as PEA2PbI4.38 Therefore, the carrier mobilities utilized in our theoretical calculations were consistent with optimized transport conditions and oriented films.
To analyze the influence of the transparent conducting oxide (TCO) and short-wavelength photons, a 300 nm-thick FTO film was included in the SCAPS-1D simulation. As shown in Fig. S3, inclusion of FTO led to an efficiency of 20.28% compared with 20.79% attained using the default TCO arrangement, with the minor decrease originating from reductions in JSC and FF, while VOC remained constant. We can state that the simulated 2DRP PSC was reliable even with applying TCO and that UV photons did not have a significant role in determining the overall performance.
To evaluate the effectiveness of our simulated PSC, we compared the PCE of our designed PSC with previously reported devices. Fig. 5 displays a comparative diagram for the state-of-the-art rigid PSCs reported in the literature. As exhibited in Fig. 5, the simulated PSC shows a competitive performance, positioning among the top-performing PSCs. This emphasizes how successfully MTEA+ cation spacer incorporation and material optimization improve device stability and efficiency.
![]() | ||
| Fig. 5 Summary of the performances for the state-of-the-art rigid PSCs with 2DRP absorber reported in the literature.13,17,34,39–48 | ||
Supplementary information (SI) is available. See DOI: https://doi.org/10.1039/d5ya00364d.
| This journal is © The Royal Society of Chemistry 2026 |