Abstract

Since the discovery of graphene in 2004, the unique properties of two-dimensional materials have sparked intense research interest regarding their use as alternative materials in various photonic applications. Transition metal dichalcogenide monolayers have been proposed as transport layers in photovoltaic cells, but the promising characteristics of group IV–VI dichalcogenides are yet to be thoroughly investigated. This manuscript reports on monolayer Ge2Se2 (a group IV–VI dichalcogenide), its optoelectronic behavior, and its potential application in photovoltaics. When employed as a hole transport layer, the material fosters an astonishing device performance. We use ab initio modeling for the material prediction, while classical drift–diffusion drives the device simulations. Hybrid functionals calculate electronic and optical properties to maintain high accuracy. The structural stability has been verified using phonon spectra. The E–k dispersion reveals the investigated material’s key electronic properties. The calculations reveal a direct bandgap of 1.12 eV for monolayer Ge2Se2. We further extract critical optical parameters using the Kubo–Greenwood formalism and Kramers–Kronig relations. A significantly large absorption coefficient and a high dielectric constant inspired the design of a monolayer Ge2Se2-based solar cell, exhibiting a high open circuit voltage of Voc = 1.11 V, a fill factor of 87.66%, and more than 28% power conversion efficiency at room temperature. Our findings advocate monolayer Ge2Se2 for various optoelectronic devices, including next-generation solar cells. The hybrid quantum-to-macroscopic methodology presented here applies to broader classes of 2D and 3D materials and structures, showing a path to the computational design of future photovoltaic materials.

Introduction

Reducing fossil fuels and their harmful environmental impact requires improvements in green, sustainable energy sources. Among the various sources of green energy generation, solar energy has been identified as the most promising and expedient because it has the potential to address the current energy demand without making a hazardous impact on the environment . Henceforth, researchers have made continuous efforts to design efficient and robust PV devices and solar cells. The systematic study of various solar cells in the last few decades has led to many successful breakthroughs in terms of the stability, efficiency, and cost of PV technology. In the past few decades, perovskite solar cells (PSCs) have emerged as a groundbreaking technology in the field of renewable energy because of their remarkable efficiency and relatively low manufacturing cost . These solar cells are based on perovskite-structured compounds, which have demonstrated excellent light absorption, charge-carrier mobilities, and tunable bandgaps .

Despite the rapid advancements in PSC technology, some critical issues, such as long-term instability, poor device scalability, and the use of toxic compounds, have yet to be resolved . Recent advancements have seen the integration of two-dimensional (2D) materials with PSCs, opening new avenues for enhancing their performance and addressing current challenges with the PSCs . Integrating 2D materials in PSCs can improve their performance. The 2D materials can provide protective layers that work like a shield to the perovskite materials to protect them from environmental degradation caused by moisture and oxygen and provide better device scalability to the PSCs, which makes large-scale production more feasible . Also, the mechanical robustness and flexibility of the 2D materials will give extra room to design the PSCs for a wide range of applications areas . The high carrier mobility and enhanced optoelectronic characteristics can improve the device performance of conventional PSCs , and the use of non-toxic 2D materials can improve the performance of PSCs without health risks.

Many researchers recently proposed an ultrathin transport layer-based solar cell as the most expedient photo device for sunlight harvesting. The general scheme of a solar cell comprises various layers including the active/absorber layers, the electron transport layer (ETL) and the hole transport layer (HTL) . Both HTL and ETL play a crucial role in achieving a high performance of PV devices. The most common HTL material is spiro-OMeTAD, but it is very expensive . Furthermore, a large number of materials have been reported as ETLs/HTL materials, such as carbon , copper iodide , CuSbS2, NiO, CuSCN , ZnO, and CdS . These materials promote efficient extraction of charge carriers and suppress recombination between the perovskite film and the electrodes. But in all these cases, the transport layer thicknesses are taken as sufficiently high. During the last few decades, researchers claimed that ETLs/HTLs of two-dimensional materials could achieve high power conversion efficiency (PCE) and optimum device performance . The key features of the 2D materials, rendering them viable for PV applications were discussed in .

Because of the advantages of 2D materials, tremendous efforts have been made to invent highly stable and efficient solar cells based on 2D materials. A number of remarkable works have been recently reported using numerous 2D materials, including MXene’s, transition-metal dichalcogenides (TMDCs) and van der Waals structures . Among the various sets of 2D materials, TMDCs (with the general formula of MX2, where M is a transition metal and X is a chalcogen) attract the broad attention of the research community regarding their extensive applications in photodevices. This is mainly because of the exceptional optoelectronic behavior of the TMDCs, especially the layer-dependent bandgap and the high absorption of incident sunlight . Motivated by this, various TMDCs, including MoS2, WS2, WSe2, and TiS2, have been used in PSCs, either as HTL or as ETL. Yin et al. reported a PCE of approximately 17.3% with TiS2 nanosheets suspended in isopropanol alcohol in a PSC. 2D sheets of SnS2 have been proposed by Zhao et al. as a prospective ETL for PSCs. It achieved a maximum PCE of around 20%. Similarly, other PSCs using different 2D layers of TMDCs have been reported (summarized below in Table 3). Despite the pervasive use of TMDC-based 2D materials in solar cells, their device performance still leaves room for improvement.

Recently, a new class of dichalcogenides (group IV–VI metal chalcogenides) has been successfully synthesized . Its remarkable features, such as low cost, abundance, environmental friendliness, and many interesting physical properties, indicate a broad range of possible applications in sustainable electronics and photonics . Notably, a high absorption coefficient has been reported in for an SnS monolayer across the direct absorption edge at 1.3–1.5 eV, rendering it potentially applicable in solar cells and photodetectors. Similarly, Huang et al. showed that a band transition (from indirect to direct) occurs in SnSe2 when moving from bulk to monolayers, elucidating its application potential in optoelectronic devices. More group IV–VI dichalcogenides have been discussed in for different applications. Although some inceptive work has been done using monolayers of metal dichalcogenides, there are still plenty of opportunities to explore. In line with the same, we have studied the thermoelectric behavior of monolayer Ge2Se2 in our previous work .

This article aims to investigate the optoelectronic properties of the group IV–VI dichalcogenide monolayer Ge2Se2 using density functional theory (DFT) and to demonstrate its potential application in photovoltaic solar cells. The proposed Ge2Se2 monolayer exhibited excellent thermodynamic stability, higher carrier mobilities (due to the presence of valleys in the CB/VB), and good optical response (an interband transition in the visible region). Considering the proper design criteria and selecting an appropriate material for absorbers/ETL, we have designed and modeled the solar cell in a FTO–TiO2–CsSn0.5Ge0.5I3–Ge2Se2–Ag configuration, which outperforms, in terms of their key parameters, other reported works in the literature. The proposed solar cell configuration is environmentally friendly (absence of toxic materials), of low cost, and highly efficient. The computational methodology followed here is depicted in Figure 1, and we provide detailed corresponding instructions and references in the Methods section.

Figure 1: Optoelectronic property investigation of a group IV–VI dichalcogenide monolayer (Ge2Se2), and its application potential for hole transport layers in Ge2Se2/CsSnGeI3/TiO2-based solar cells. Our computational approach combines density functional theory and atomistic modeling with (macroscopic) device modeling to obtain an otimized device performance.

Results and Discussion

The computational approach can be better explained by separating it into two subsections, namely (a) material simulation using DFT and (b) device simulation using a macroscopic approach, which will be discussed in section Computational Methods.

Structural and stability analysis

The schematic atomistic model for monolayer Ge2Se2 is shown in Figure 2. For better understanding, a side view of monolayer Ge2Se2 with the corresponding lattice parameters is given in Figure 2a, and the top view of the supercell of the order 4 × 4 is shown in Figure 2b. From Figure 2a, we can observe that the unit cell of monolayer Ge2Se2 consists of four atoms (two for each Ge and Se). More precisely, we can observe two sub-layers corresponding to a crystallogen (Ge) and a chalcogen (Se), which are designate as top and bottom moieties. From Figure 2, we can also notice that the structure of monolayer Ge2Se2 is identical to the puckered phosphorene (simple orthorhombic structure) with space group symmetry of Pmn21 () in 2D space.

Figure 2: Geometry and stability of Ge2Se2 monolayers. (a) Side view with the corresponding structural parameters and (b) top view. The exhibited structure is a simple orthorhombic structure with space group symmetry of Pmn21 () in 2D space. (c) Brillouin zone, where the Γ-X-S-Y-Γ path is used to plot the phonon spectra, and (d) phonon dispersion, confirming the structure’s dynamical stability (Figure 2a,b were redrawn from , Figure 2d was re-plotted from ).

The geometry-optimized lattice parameters a and b are found to be 4.13 and 3.99 Å, respectively. The bond lengths between consecutive crystallogen and chalcogen (Ge–Se) are 2.59 and 2.67 Å in the horizontal and vertical directions, respectively. All optimized structural parameters excellently match with previously reported values . In addition to the lattice parameters for the monolayer Ge2Se2, we have optimized the vacuum slab in a normal direction to avoid interlayer interactions. The optimized vacuum level was found to be 23 Å. Furthermore, to test the stability of the crystal structure, we have computed the phonon band dispersion for monolayer Ge2Se2 within the first Brillouin zone (Figure 2c). The calculated phonon spectra along the high symmetry path Γ-X-S-Y-Γ in the first Brillouin zone are shown in Figure 2d.

The phonon spectrum for monolayer Ge2Se2 has twelve vibrational modes, of which three are acoustical (low-frequency) modes (transverse acoustic, longitudinal acoustic, and flexural acoustic). The remaining nine (high-frequency) modes correspond to optical modes. The flexural acoustic mode is an out-of-plane transverse acoustic mode, similar to those in other two-dimensional materials such as graphene, phosphorene, and stanene, demonstrating a quadratic nature near the Γ point . Figure 2d clearly depicts that all phonon bands throughout the spectrum have non-imaginary frequency, further confirming the investigated structure’s dynamic stability.

Electronic and transport properties

The electronic properties of materials play a crucial role in predicting the material behavior and transport parameters and, hence, its prospective applications. We investigated the electronic properties of monolayer Ge2Se2 using its electronic band structure and density of states, as shown in Figure 3. To maintain a high degree of accuracy, we calculated the band structure using the HSE06 functional, and the band dispersion along the high symmetry paths Γ-X-S-Y-Γ in the first Brillouin zone, as shown in Figure 3.

Figure 3: The electronic band structure of monolayer Ge2Se2, revealing a direct bandgap of 1.12 eV, and the corresponding density of states (The electronic band structure was re-plotted from ).

Figure 3 depicts that the monolayer Ge2Se2 is a direct-bandgap semiconductor with a bandgap of the order of 1.12 eV. Valence band maximum (VBM) and conduction band minimum (CBM) are located along the Γ-X path. The computed bandgap value and its dispersion nature are consistent with earlier reported works . The optimum bandgap value and its CBM/VBM positions make it suitable for photovoltaic applications and provide a guideline for selecting the absorber layers in the solar cells. Notably, multivalley and flat bands are present in the energy band structure. The formation of multiple valleys provides extra carrier pockets for transportation, further increasing carrier mobility. Also, the larger valley degeneracy will increase the density of states (DOS) effective mass without influencing the carrier mobility. We have also calculated the total DOS to understand and justify the band structure calculation. Figure 3 shows that the calculated DOS is consistent with the electronic band structure and the bandgap values, further validating the results.

Other properties that will be useful in the simulation of solar cells, such as the effective mass of charge carriers (electrons/holes), conduction/valence band density of states, electron/hole mobility, electron affinity, and work function can be derived from the initial band energy calculation.

We calculated the effective masses of electrons and holes as = 0.167me and = 0.1768me, respectively, which are very close to the values ( = 0.17me, and = 0.17328me) reported in the literature . With that, we estimated the effective DOS values in the CB and the VB to be 0.1732 × 1019 /cm3 and 0.1887 × 1019 /cm3, respectively. The electron and hole thermal velocities for monolayer Ge2Se2 were estimated as 9.43 × 105 m/s and 2.668 × 105 m/s, respectively. The predicted charge carrier mobilities were found to be 187.003 cm2/(V·s) and 300.451 cm2/(V·s) for holes and electrons, respectively. The high carrier mobility further suggests a potential application as a transport layer in solar cells as well as in the channel materials of FETs. It is important to note that although in the case of solar cells, the transmission of charge carriers takes place in the out-of-plane direction, which is usually considered to be low compared to the in-plane mobility (especially in the bulk materials), the charge carriers still move quickly across the thin HTL because of the ultrathin layer and the strong electric field gradient from the active layer to the electrode. Finally, higher in-plane mobility of the monolayers can compensate the potential limitations of out-of-pane mobility for its prospective uses in solar cells, V-FETs, TFETs, and memristors .

In the present calculation, we have calculated the average planar electrostatic potential along the plane as demonstrated in Figure 4. The figure depicts the vacuum energy level and the positions of HOMO and LUMO. The calculated work function and electron affinity for monolayer Ge2Se2 are about 3.313 and 3.873 eV, respectively, which are comparable with previously reported values of 3.39 eV and 3.9 eV .

Figure 4: Macroscopic potential distribution vs lattice vector in the normal direction for monolayer Ge2Se2. In a periodic system, the electron affinity and work function is defined as EA = EVac − ELUMO, and Eϕ = EVac − EF, where EF is the Fermi energy level, Eϕ is the work function, EA is the electron affinity, EVac is the vacuum energy level, and ELUMO is the LUMO energy level, respectively. The calculated work function and electron affinity are about 3.313 and 3.873 eV, respectively.

A high degree of agreement for all estimated electronic parameters with earlier reported values validates the accuracy of our computation.

Optical properties

As discussed in the previous section, monolayer Ge2Se2 has structural stability and exhibits an optimal electronic behavior with a direct bandgap of 1.12 eV. The DFT results indicate potential applications of the investigated structure in optoelectronic devices. The optical responses at varying photonic energies and the corresponding optical wavelength are summarized in Figure 5, depicting the dielectric function εr (real and imaginary parts) and the refractive index n for monolayer Ge2Se2. The real and imaginary components of the dielectric function correspond to dispersive and absorption effects of the material. We observe that the investigated structure exhibits strongly anisotropic optical properties for which the parametric values depend on the direction of the chosen plane.

Figure 5: Optical responses of monolayer Ge2Se2; real (solid) and imaginary (dashed) part of the dielectric constant and refractive index vs wavelength (left) and photon energy (right), for in-plane and out-of-plane excitation relative to the sheet. The teal-colored areas mark the regions of potential metallic behavior (Re(εr) < 0 or, equivalently, Re(n) < Im(n)), occurring between 180 and 517 nm (2.4 and 6.9 eV) for in-plane, and between 188 and 207 nm (6.0 and 6.6 eV) for out-of-plane excitation.

Dielectric function and absorption coefficient

The static dielectric constants (real part of εr at ω = 0), along the in-plane (XX) and out-of-plane (ZZ) direction were calculated as 6.37 and 2.4, respectively, suggesting that monolayer Ge2Se2 is highly polarized along the in-plane direction.

The imaginary part of the dielectric function shows peaks between 2 and 3 eV, along the XX direction, which indicate an interband transition in the visible region (Figure 5). Interestingly, we can observe tracking bands (the region where the energy difference between the conduction and valence bands is approximately constant ) around the high-symmetry point S and between X and Y in the first Brillouin zone. Around these points, we notice a small gradient and, thus, large DOS, which further leads to higher transition possibilities.

The absorption coefficient elucidates the rate at which light intensity diminishes while entering the material. It primarily depends upon the imaginary part of the refractive index k(λ) = Im and the wavelength of the incident light λ, as per αabs = 4πk(λ)/λ. The highest peak of the absorption coefficient lies in the visible region, which we attribute to the direct transition from VBM to CBM. The other peaks correspond to the transition from the top valence band to the conduction band in a band region of the Brillouin zone around the S point. The maximum absorption coefficient is 5 × 105 cm−1 along the XX direction. The highest peak for the ZZ direction is calculated as 5.8 × 105 cm−1, corresponding to the UV region. We attribute this maximum to the transition from VBM to CBM, while the other peaks designate a transition from VBM to CBM around a wide range of spectra at the S point.

Refractive index and extinction coefficient

The real part of the material’s refractive index defines the angle of refraction that the light undergoes while entering a medium. Higher refractive indices refer to lower refraction angles measured from the interface normal. The static refractive indices for monolayer Ge2Se2 in the direction parallel to the plane (XX) and normal to the plane (ZZ) are approximately 2.5 and 1.4, respectively, as shown in Figure 5. These refractive indices could enable the use of the investigated material as an inner layer coating between the absorber and substrate/electrodes in the design of solar cells.

Device Modeling

The high carrier mobility, optimum bandgap, and suitable optical characteristics of monolayer Ge2Se2 as discussed in the previous section, further motivated us to design a photovoltaic solar cell using monolayered Ge2Se2. The easy fabrication, structural stability, and proper band alignment with the absorber material make it a promising hole transport layer (HTL) in the proposed solar cell. The proposed solar cell is a low-cost, environmentally friendly, scalable Perovskite solar cell (PSC), with Ge2Se2 as a HTL. The schematic of the proposed solar cell is shown in Figure 6a.

Figure 6: Schematic of the proposed PSC using Ge2Se2 as HTL; (a) device setup consisting of stacked layers of FTO–TiO2–CsSn0.5Ge0.5I3–Ge2Se2–Ag; (b) band offset among different PSC layers, demonstrating the ease of charge-transfer from the active layers to the respective transport layers.

The proposed structure consists of a lead-free Perovskite layer as an absorber, titanium dioxide (TiO2) as ETL, and Ge2Se2 as HTL. A 100 nm thick layer of FTO is used to encapsulate the transport and absorber layers and treated as a conductor in this simulation. Regarding an absorber material corresponding to the proposed HTL, a minimum valence band offset between these two layers has been considered as a necessary condition for achieving optimum performance. In addition to this, environmental friendliness, stability, and cost-effectiveness were also considerations to choose CsSn0.5Ge0.5I3 as an absorber layer. After selecting the absorber layer, the ETL has been chosen based on a matching conduction band offset with the absorber layer and the high electron mobility. Among the various options, TiO2 offers minimum band offset and ease of synthesis. The front and back contacts are supposed to be ohmic and are made of silver (Ag) to maintain proper band alignment and reduce the band offset. It has been reported that tuning the CBM and VBM positions to keep a low band offset leads to superior device performance. Figure 6b exhibits the position of CBM/VBM in consecutive layers of the proposed structure. This further demonstrates the ease of charge transfers from the absorber to the respective transport layers.

The thickness and other simulation parameters for the active layer and the ETL were taken from the literature , while all the simulation parameters for Ge2Se2 have been derived from the DFT calculations as discussed in the previous section. All simulation parameters are summarized in Table 1

Table 1: Simulation parameters for the proposed PSC.

Device parameter FTO Ge2Se2 [our calculations] CsSn0.5Ge0.5I3 TiO2 Bandgap (Eg) (eV) 3.4 1.12 1.5 3.26 Electron affinity (eV) 4.5 3.31 3.9 4.2 Dielectric constant (εr) 9.1 6.37 28 10 Conduction band effective DOS (cm−3) 1.1 × 1019 0.1732 × 1019 1.0 × 1019 2.0 × 1017 Valence band effective DOS (cm−3) 1.1 × 1019 0.1887 × 1019 1.0 × 1019 6.0 × 1017 Electron thermal velocity (m·s−1) 1 × 107 9.47 × 105 1 × 107 1 × 107 Hole thermal velocity (m·s−1) 1 × 107 2.66 × 105 1 × 107 1 × 107 Electron mobility (cm2·V−1·s−1) 20 300.451 974 100 Hole mobility (cm2·V−1·s−1) 10 187.003 213 25 Performance estimation

Using the proposed model, we have estimated the performance of the solar cell in terms of open-circuit voltage (Voc), short-circuit current density (Jsc), fill factor (FF), and power conversion efficiency (η). The performance of solar cells largely depends on various parameters, including the thickness of transport layers, absorber layers, and absolute temperature. Therefore, to get the optimum performance of the proposed device, we have optimized the performance of the cell module by varying the layer thicknesses of HTL, ETL, and absorber. Also, we have tested the performance of the device under the influence of varying defects densities in the active layer and the interfaces. The optimization corresponding to different parameters is summarized in the following sections.

Effect of HTL thickness variation

While the absorber layer is the crucial component of a solar cell, the role of the HTL cannot be overlooked because the choice of HTL plays a significant role in ensuring efficient charge transport, reducing recombination losses, and overall device stability. Therefore, numerous efforts have been made in the past to investigate suitable HTL materials for efficient solar cells . The main function of the HTL is to capture the holes generated from the absorber layer and transfer them to the respective electrode. The thickness of the HTL also influences the device performance significantly. To optimize the HTL thickness for optimal device performance, we assumed a layer thickness of monolayer Ge2Se2, ranging from 1 to 10 nm. Figure 7a depicts the change of the performance parameters upon varying the HTL thickness.

Figure 7: PSC performance parameters as functions of (a) HTL thickness, (b) ETL thickness, (c) absorber thickness, and (d) absorber defects. The device exhibits optimal performance at layer thicknesses of ETL, absorber, and HTL of 20 nm, 2 µm, and 1 nm, respectively.

It can be observed from Figure 7a that the device parameters decrease with higher HTL thickness. This behavior is mainly due to the appearance of a higher series resistance at larger HTL thicknesses. Also, a thick HTL may increase the chances of charge carrier recombination. Nevertheless, the performance degradation is very slight because of the ultrathin HTL layer proposed here. It is worth mentioning that there is an incremental rise in the short-circuit current density with thicker HTL. This is due to less carrier losses owing to reduced recombination at the interface of HTL and absorber layer and better surface cleavage because of the thick HTL . Figure 7a indicates that 1 nm of HTL yields superior performance in terms of PCE, Voc, Jsc, and FF with values of 28.148%, ≈1.11 V, 28.70 mA·cm−2, and ≈87.77% FF, respectively.

Effect of ETL thickness variation

In conjunction with the HTL, the ETL is responsible for collecting electrons generated in the active absorber layer. It is also necessary to optimize the ETL thickness to get an optimal device behavior. The change of the device parameters with varying ETL thickness is shown in Figure 7b. A higher thickness of the ETL leads to a worse device performance. When the layer thickness increases, the series resistance in the cell will also increase, and the light transmittance through the ETL decreases. The latter affects photogeneration in the perovskite layer, reducing Jsc and, consequently, the cell efficiency. Also, one can notice that the fill factor initially increases and then decreases mainly due to the small reverse saturation current for lower values of the ETL thickness. Later on, at higher ETL thicknesses, the fill factor decreases with increase in series resistance, and the cell efficiency decreases . Since the ETL thickness is low, the impact on the performance degradation is insignificant. After examining all parameters, we found an optimized ETL thickness of 20 nm in our proposed device. The optimized values of the solar cell parameters at the optimal ETL thickness are 28.11%, ≈1.11 V, 28.68 mA·cm−2, and ≈87.81% for PCE, Voc, Jsc, and FF, respectively.

Effect of absorber thickness variation

The absorber layer thickness is the most important and critical device parameter. The role of the absorber layer is to absorb the photons and generate charge carriers required for the current flow. It is intuitive to have a thick absorber layer for a higher current density. However, because the absorber thickness further affects the recombination rate, optimizing it for optimal device performance is necessary. The impact of varying absorber thickness on the different device parameters is summarized in Figure 7c. The absorber layer thickness was varied from 1 to 10 µm to obtain the optimal performance of the proposed device. Figure 7c shows that the short-circuit current density (Jsc) increases with higher absorber thickness up to a thickness of 4 µm, before reaching a constant value. The maximum short-circuit current density (Jsc) was calculated as 28.86 mA·cm−2 corresponding to an absorber thickness of 4 µm. The current variation followed the Beer–Lambert law, which justifies that more absorption will occur when the thickness of the absorber layer increases. A thicker absorber layer allows more photons to get absorbed, which leads to the generation of more electron–hole pairs. These generated electron–hole pairs, collected by the respective transport layers, contribute to the higher current density .

It is worth noticing that there is a significant impact of the absorber thickness on almost all parameters. The PCE increases from 27.4% to 28.11%, and the short-circuit current Jsc value reaches up to 28.86 mA·cm−2. We observed that there is a significant reduction in FF and Voc at higher absorber thicknesses, primarily because of the increasing series resistance and higher losses . We set the optimal absorber thickness at 2 µm, above which the parameter variation was only very little.

Effect of defects in the absorber layer

Halide pervoskites are generally considered as defect-resilient; yet, deep level defects shows significant impact on the device performance . Therefore, it is important to analyze their impact on the cell performance for better device design. Figure 7d shows the variation of the cell parameters at varying absorber defect densities from 1014 cm−3 to 1018 cm−3. Figure 7d shows that higher defect densities lead to worse device performance. This is due to the decrement in the lifetime of the charge carriers because of the shorter diffusion length, which results from the higher defect densities . We can observed that the performance parameters decrease significantly above defect densities of 1015 cm−3.

Effect of interface defects between active and transport layers

The photovoltaic behavior of halide perovskites is significantly affected by the properties of the interfaces. This is because a majority of defects arises at the junction between two layers during the fabrication of the device, caused by mismatches in the lattice structures, which lead to dangling bonds. These defects at the junctions play a critical role in the performance of the device as they contribute to increased local recombination at the interface. Within this particular PSC, two critical interface layers exist, namely, the interface between the absorber and the ETL and the interface between the absorber and the HTL. We analyzed the device performance by varying the defect densities at both interfaces within the range of 106 cm−3 to 1015 cm−3. Because of the ultrathin HTL layer, the absorber/HTL interface is more sensitive, and the device performance is anticipated to degrade above interface defect densities of 107 cm−3. However, the ETL/absorber interface exhibits excellent performance up to interface defects of 1012 cm−3 (Table 2), after which the device performance starts to deteriorate (Figure 8).

Table 2: Performance of the proposed PSC under 1.5G solar illumination.

Device parameter Value Open circuit voltage (Voc) 1.11 V Short circuit current density (Jsc) 28.70 mA·cm−2 Fill factor (FF) 87.66% Power conversion efficiency (PCE) 28.10%

Figure 8: Variation of PSC performance parameters with interface defects between (a) ETL and active layer, and (b) active layer and HTL. Because of the ultrathin HTL layer, the interface at the absorber/HTL is more sensitive, and the device performance degrades for interface defect densities above 107 cm−3. However, the ETL/absorber interface exhibits an excellent performance up to the interface defects of 1012 cm−3, after which the device performance starts to deteriorate.

The computed current density vs voltage curve is shown in Figure 9a, and the energy band diagram corresponding to the different layers is demonstrated in Figure 9b. We have also compared the proposed device with the combination of different transport layers in Table 3. It demonstrates that the proposed PSC structure excels in terms of performance parameters compared to other reported structures.

Figure 9: (a) Current density vs voltage curve for the proposed PSC, demonstrating a short-circuit current density and an open-circuit voltage of 28.7 mA·cm−2 and 1.11 V, respectively. (b) Energy band diagram for the proposed solar cell, where EC stands for conduction band and EV stands for valence band. FN and FP stand for quasi-Fermi levels of electrons and holes, respectively, at T = 300 K.

Table 3: Comparison of the PSC performance with earlier reported works.

Device configuration Voc (V) Jsc (mA·cm−2) FF (%) PCE (%) Ag/TiO2/PSK/MoS2/ITO 0.93 26.24 83.0 29.43 Au/Cu2O/PSK/TiO2/FTO 0.96 15.8 59.0 8.93 Au/MoS2/PSK/WS2/FTO 0.96 27.3 87.62 22.17 Au/Spiro–OMeTAD/PSK/TiS2/FTO 0.95 20.05 66.90 12.75 Ag/Ge2Se2/PSK/TiO2/FTO [This work] 1.11 28.70 87.66 28.10 Computational Methods

The methodology presented here (illustrated in Figure 10) is more generally valid and will apply to broader classes of 2D and 3D materials and structures. Therefore, we believe this section might serve as a mini tutorial for the computational design of future photovoltaic materials.