An Optimization Path for Sb2(S,Se)3 Solar Cells to Achieve an Efficiency Exceeding 20%

Antimony selenosulfide, denoted as Sb2(S,Se)3, has garnered attention as an eco-friendly semiconductor candidate for thin-film photovoltaics due to its light-absorbing properties. The power conversion efficiency (PCE) of Sb2(S,Se)3 solar cells has recently increased to 10.75%, but significant challenges persist, particularly in the areas of open-circuit voltage (Voc) losses and fill factor (FF) losses. This study delves into the theoretical relationship between Voc and FF, revealing that, under conditions of low Voc and FF, internal resistance has a more pronounced effect on FF compared to non-radiative recombination. To address Voc and FF losses effectively, a phased optimization strategy was devised and implemented, paving the way for Sb2(S,Se)3 solar cells with PCEs exceeding 20%. By optimizing internal resistance, the FF loss was reduced from 10.79% to 2.80%, increasing the PCE to 12.57%. Subsequently, modifying the band level at the interface resulted in an 18.75% increase in Voc, pushing the PCE above 15%. Furthermore, minimizing interface recombination reduced Voc loss to 0.45 V and FF loss to 0.96%, enabling the PCE to surpass 20%. Finally, by augmenting the absorber layer thickness to 600 nm, we fully utilized the light absorption potential of Sb2(S,Se)3, achieving an unprecedented PCE of 26.77%. This study pinpoints the key factors affecting Voc and FF losses in Sb2(S,Se)3 solar cells and outlines an optimization pathway that markedly improves device efficiency, providing a valuable reference for further development of high-performance photovoltaic applications.


Introduction
With the increasing demand for renewable energy and a heightened focus on environmental sustainability, photovoltaic (PV) technology has emerged as a crucial method for converting solar energy into electricity.Thin-film solar cells, which utilize inorganic materials such as Cu(In,Ga)Se 2 (CIGS) and CdTe, have been commercialized, offering excellent power conversion efficiency (PCE) and especially lighter weights compared to silicon-based solar cells.However, the use of expensive rare earth elements (Ga and Te) and the toxicity of Cd present challenges to the broader application of these thin-film solar cells.Recently, researchers have shown significant interest in antimony selenosulfide semiconductor materials, Sb 2 (S , Se) 3 , due to their nontoxic nature, profuse existence in the Earth's crust, excellent ambient stability, large dielectric constant (ε r > 15), elevated relative absorption coefficient (α > 10 5 cm −1 ), and adjustable bandgap (E g ) ranging from 1.1 to 1.7 eV [1,2].Despite substantial efforts by researchers to enhance the PCE of Sb 2 (S,Se) 3 solar cells, the current champion PCE for Sb 2 (S,Se) 3 thin-film solar cells remains at 10.75% [3,4].Nevertheless, these devices still suffer from a substantial open-circuit voltage (V oc ) loss (defined as E g /q − V oc ) of approximately 0.8 V, and a fill factor (FF) loss (defined as the difference between the FF of ideal solar cells and the FF observed in practice) of approximately 15%.Moreover, there is still a considerable gap to bridge before reaching the Shockley-Queisser (S-Q) limit, which is based on the detailed balance principle and predicts a maximum efficiency of 33% for single-junction solar cells under standard illumination conditions [5].Therefore, continued research endeavors aimed at mitigating these V oc and FF losses in Sb 2 (S,Se) 3 solar cells are crucial to advancing towards this theoretical benchmark.
To achieve exceptional performance in Sb 2 (S,Se) 3 solar cells, we used a solar cell capacitance simulator (SCAPS) to analyze the limiting factors, including series resistance (R s ), shunt resistance (R sh ), bulk defects, interface recombination, and energy level arrangements between different functional layers.SCAPS is a powerful solar cell simulator that has been validated by numerous studies for accurately adjusting the parameters of each functional layer and interface, simulating photovoltaic device performance, optimizing device structures, and obtaining high-efficiency solar cells.For example, in our previous work [6], we achieved a PCE of over 15% by utilizing SCAPS for guiding interface engineering strategies and optimizing the design of PbS colloidal quantum dot solar cells.Similarly, Chen et al. [7] successfully prepared Sb 2 Se 3 solar cells with a certified efficiency of 6.5% using EDT-treated PbS colloidal quantum dots as hole transport layers (HTLs), guided by SCAPS simulation results.Bertens et al. [8] used SCAPS to predict the effect of doping concentration in the active layer on the performance of PbS colloidal quantum dot solar cells and achieved actual devices that matched the predicted performance.
This exhaustive investigation delves into the multifaceted factors that currently impede the enhancement of Sb 2 (S,Se) 3 solar cell efficiency, with the ultimate goal of identifying optimization strategies to attain a PCE exceeding 20%.By utilizing numerical methods, we constructed a theoretical model for the most efficient Sb 2 (S,Se) 3 solar cell structure reported thus far, featuring a fluorine-doped tin oxide (FTO)/CdS/Sb 2 (S,Se) 3 /Spiro-OMeTAD/Au configuration.Initially, through meticulous optimization of R s and R sh , we significantly reduced the FF loss from 10.79% to 2.80%, increasing the PCE from 10.72% to 12.18%.For the device exhibiting 12.18% efficiency, by refining the alignment of the conduction band between the electron transport layer (ETL) and the absorber layer, we achieved a notable increase in V oc , increasing it from 0.64 V to 0.76 V.This adjustment minimized the V oc loss and increased the PCE to 15.31%.Building upon this improved device, we further optimized the interface recombination, achieving a substantial 43% decrease in the V oc loss relative to that of the initial device, with FF losses decreasing below 1% and a substantial increase in the PCE to 23.75%.In the final phase of optimization, we targeted the absorber defects and thickness, resulting in an increase in the PCE to 26.77%, accompanied by an increase in the short-circuit current density (J sc ) to 30.93 mA•cm −2 .The simulation results revealed that under conditions of low V oc and FF, the internal resistance exerted a more substantial constraint than nonradiative recombination on the FF.Additionally, the spike-like energy band offset (EBO) increase at the FTO/ETL interface could impede the formation of a flat or slightly spiked EBO at the ETL/absorber interface, thus impeding further enhancements in the PCE.Importantly, mitigating interface recombination at the ETL/absorber interface has emerged as a paramount strategy for enhancing the overall performance of Sb 2 (S,Se) 3 solar cells, surpassing the importance of reducing recombination at the absorber/HTL interface.Remarkably, as the absorber layer thickness increased to 600 nm, the Sb 2 (S,Se) 3 solar cell achieved a noteworthy J sc exceeding 30 mA•cm −2 , consequently achieving a PCE of 26.77%.This study elucidates the critical factors influencing the performance of Sb 2 (S,Se) 3 solar cells, providing valuable guidance for further experimental work in this field.

Simulation Principles and Details
SCAPS has been meticulously crafted and generously offered by the University of Gent [16,17].It is specifically designed to simulate a wide array of solar cells, including crystalline silicon, thin films, and emerging photovoltaic technologies.The software incorporates sophisticated numerical algorithms to solve the relevant physical equations governing charge transport, recombination, and optical properties within the device.
The primary function of SCAPS is to solve the Poisson equation related to the electrostatic potential (ψ) and the continuity equation controlling the dynamics of free electrons and holes in one-dimensional heterojunctions [18].The specific Equations are as follows:

Simulation Principles and Details
SCAPS has been meticulously crafted and generously offered by the University of Gent [16,17].It is specifically designed to simulate a wide array of solar cells, including crystalline silicon, thin films, and emerging photovoltaic technologies.The software incorporates sophisticated numerical algorithms to solve the relevant physical equations governing charge transport, recombination, and optical properties within the device.
The primary function of SCAPS is to solve the Poisson equation related to the electrostatic potential (ψ) and the continuity equation controlling the dynamics of free electrons and holes in one-dimensional heterojunctions [18].The specific Equations are as follows: where ε represents the dielectric constant; q represents the electron charge; p and n symbolize the free hole and electron densities; N A and N D represent the densities of the acceptor and donor types, respectively; ρ n and ρ p indicate the densities of trapped electrons and holes, respectively; J n and J p denote the current densities of the electrons and holes, respectively; G designates the photogeneration rate; and, R n and R p are defined as the recombination rates of the electrons and holes, respectively.Additionally, µ n and µ p denote the electron and hole mobilities, respectively, and E Fn and E Fp represent the acceptor and donor quasi-Fermi levels, respectively.
To model the device structure via SCAPS, we compiled material parameters obtained from calculations and the relevant literature [3,4,9,[18][19][20][21][22][23][24][25], which are summarized in Table S1 (see the Supplementary Material).Additionally, we standardized the thermal velocities of electrons and holes to 10 7 cm•s −1 across all the layers.The specific bulk defects for the absorber layer are presented in Table S2, while Table S3 provides the detailed interface characteristics essential for the simulation.The absorption spectral data for FTO and CdS were obtained from the default file provided by SCAPS software (version 3.3.10;University of Gent, Belgium) support, and the absorption spectral data for the Sb 2 (S,Se) 3 material obtained from reference [4] were utilized in the simulation process.All the simulations were carried out at 300 K, neglecting radiative and Auger recombination, under the AM 1.5 G spectrum.

Device Validation and Optimization Route
Simulation 1 of the Sb 2 (S,Se) 3 solar cell was established utilizing the parameters outlined in Tables S1-S3 of the Supplementary Material.When R s = 3 Ω•cm 2 and R sh = 375 Ω•cm 2 , the simulation results strongly agree with the experimental findings reported in [4], particularly the current density-voltage (J-V) characteristic curves and external quantum efficiency (EQE) spectra, as evident in Figure 1b,c, respectively.To further confirm the accuracy of the model, a comparative analysis was conducted on how V oc correlates with variations in light intensity.According to the data presented in Figure 1d, the ideal factor (n ID ) values of the experimental and simulated devices are 1.35 and 1.39, respectively, revealing strong consistency between the measurements and simulations.Additionally, an n ID value between 1 and 2 indicates that interface recombination predominates in Sb 2 (S,Se) 3 solar cells [26].As depicted in Figure 1e, the band alignment diagram clearly illustrates that the Sb 2 (S,Se) 3 absorber layer is fully depleted, indicating a strong electric field throughout the entire layer.This robust electric field serves as a driving force, greatly enhancing the separation efficiency of photogenerated electron-hole pairs.
In our simulated device, the Sb 2 (S,Se) 3 absorber has an E g of 1.43 eV, which approaches the nearly ideal bandgap (E g = 1.34 eV) of the S-Q efficiency limit [5,27].However, the actual device performance only reaches a maximum efficiency of 10.75%, which is significantly below the S-Q efficiency limit.For instance, at an E g of 1.43 eV, the S-Q limit predicts a PCE of 32.54% with a V oc of 1.17 V and an FF of 89.53% [5].A strong correlation exists between FF and V oc , where the maximum achievable FF (FF m ) theoretically depends solely on V oc [3,28,29], assuming a negligibly small R s (R s ≈ 0) and an ideal R sh (R sh → ∞).

FF m =
ν oc − ln(ν oc + 0.72) The relationship between FF m and V oc is primarily influenced by n ID , which represents the nonradiative recombination rate.However, in actual devices, R s and R sh cannot be neglected.R s and R sh are considered in Equation ( 5), as shown in Equation (6) [30,31].
Applying Equations ( 5) and ( 6), the FF and FF m of Sb 2 (S,Se) 3 solar cells were determined on the basis of empirical V oc -FF plots, where the n ID value was varied from 1 to 2. According to Figure 1f, with R s = 0 Ω•cm 2 and R sh = 10 7 Ω•cm 2 , at a V oc of approximately 0.63 V (that is, from point 1 to point 2 on the graph), the calculated FF increased from 68.35% to 78.91%.On the other hand, by decreasing n ID from 1.39 to 1.0, the FF increased from 68.35% at 0.63 V to 79.94% at the S-Q limit, with a V oc of 1.17 V (that is, moving from point 1 to point 3).On the basis of these findings, after the respective optimizations, the FF constrained by the internal resistance (78.91%) was lower than that constrained by the nonradiative recombination (79.94%).Consequently, within the regimes of low V oc and FF, the internal resistance exerted a more significant limiting effect on the FF than nonradiative recombination did.These results highlight the rationale behind implementing a phased strategy employing staged optimization to reduce V oc loss and FF loss in Sb 2 (S,Se) 3 solar cells.This strategy involves initially optimizing internal resistances before addressing nonradiative recombination in Sb 2 (S,Se) 3 solar cells (that is, from point 1 to point 2 and then to point 4, as illustrated in Figure 1f).Additionally, point 5 in Figure 1f represents the high-efficiency Sb 2 (S,Se) 3 solar cell with the current lowest V oc loss [3], whereas points 1 ′ -5 ′ represent the results from five simulated devices obtained following our designed optimization route.

Optimization of the Series and Shunt Resistances
The series and shunt resistances are crucial factors influencing the performance of solar cells.Simulation 1 effectively replicates the experimental J-V curves, EQE spectra, and ideality factor (n ID ), aligning closely with published data [4].This validation highlights the predictive accuracy of subsequent simulations in predicting experimental outcomes.To explore the detailed effects of R s and R sh , R s was adjusted to fall within the range of 0 to 10 Ω•cm 2 , and the range specified for R sh extended from 10 2 to 10 7 Ω•cm 2 .The simulation calculation was carried out while the other parameters in Simulation 1 remained unchanged.As illustrated in Figure 2a-d, changes in R s have a more pronounced effect on the PCE and FF than variations in R sh do.When R s decreased to 0 Ω•cm 2 and R sh was maintained at 10 5 Ω•cm 2 , the FF notably increased from 68.35% to 79.55%.Conversely, the V oc and J sc values are least affected by changes in R s , which is consistent with previous reports [32,33].Moreover, Figure 2e demonstrates that n ID remains stable with variations in R s when R sh remains constant but shifts noticeably from 2.14 to 1.10 as R sh varies from 10 2 Ω•cm 2 to 10 5 Ω•cm 2 .Compared with Simulation 1, the J-V curve of Simulation 2 (R s decreased to 0 Ω•cm 2 and R sh increased to 10 5 Ω•cm 2 ) significantly improved the FF, whereas the PCE also increased from 10.72% to 12.57% (Figure 2f).This enhancement primarily results from a marked decrease in FF loss, decreasing from 10.79% to 2.80% as detailed in Table 1, emphasizing the critical role of R s and R sh optimization in enhancing solar cell performance.  a) V oc loss is equal to E g /q − V oc , where E g represents the bandgap of the active absorber and q is the elementary charge.(b) FF loss is equal to FF m − FF, where FF m denotes the maximum achievable FF across varying V oc values in a single solar cell.
The findings underscore the efficacy of a dual-pronged approach in enhancing solar cell efficiency and mitigating FF losses: reducing R s and augmenting R sh .An increase in R s stems from increased electrode and charge-transfer interface resistances, whereas a decrease in R sh is tied to reduced functional layer coverage and the proliferation of pinholes.Despite the inherent challenge of fully eradicating internal resistance, these results underscore the imperative for experimental endeavors to craft electrodes with exceptional conductivity and minimal sheet resistance, alongside strategies to diminish interface charge-transfer resistance.
The key determinants of this interfacial charge transfer resistance include the electrodesemiconductor contact resistance, the bulk resistance intrinsic to the semiconductor, barriers arising from interface band offsets, and recombination losses at the interface, exacerbated by impurities and defects [34].To increase the overall performance of Sb 2 (S,Se) 3 solar cells to unprecedented levels, concerted optimization endeavors targeting interface band alignment, interface recombination processes, and absorber defects are of paramount importance.

Optimization of the Energy Band Offset
Despite significant advancements in enhancing the efficiency of Sb 2 (S,Se) 3 solar cells through meticulous optimization of the internal resistance, resulting in a notable improvement in the FF, increasing the V oc remains a formidable challenge.To overcome this obstacle, it is imperative to meticulously optimize the internal energy band arrangements within the solar cell architecture.Optimal energy level alignment can facilitate the seamless transportation of photogenerated electrons toward the front contact and holes toward the back contact, emanating from the absorber layer, thereby reducing the interface-induced recombination losses in Sb 2 (S,Se) 3 solar cells.Establishing an ideal energy band offset (EBO) at the interfaces, which increases carrier transport properties and minimizes recombination losses, is of paramount importance [35].Notably, the formation of EBO, a crucial aspect of Sb 2 (S,Se) 3 solar cells, is closely linked to disparities in the electron affinity and energy gap across the material interface.Addressing these challenges through meticulous interface engineering is essential for mitigating recombination losses, enhancing V oc , and ultimately achieving higher overall efficiencies for these solar cells.The conduction band offset (CBO) existing at the junction of the ETL and absorber is fundamentally dictated by the disparity in electron affinities between these two layers (χ absorber and χ ETL , respectively), as follows: CBO = χ absorber − χ ETL (7) The interface between materials can exhibit three distinct band alignment configurations, each characterized by the CBO.A negative CBO induces a steep cliff-like EBO, whereas a positive CBO yields a more pronounced spike-like EBO.When the CBO is zero, a flat band alignment is observed [36].There is a general consensus among researchers that either a flat band alignment or a mild spike-like CBO is conducive to efficient charge transfer and the minimization of interface recombination losses [37].
Here, the band arrangements at the ETL/absorber interface were investigated by varying the χ ETL from 4.3 eV to 3.7 eV, with the χ absorber maintained at 3.81 eV.This range of χ ETL corresponded to a CBO ranging from −0.49 eV to 0.19 eV, as depicted in Figure S1a.
As the CBO shifted from negative to positive values, a distinct spike-like EBO emerged in the conduction band at the ETL/absorber interface, as expected.However, a remarkable observation was made at the FTO/ETL interface, where a considerably more prominent spike-like EBO manifested in the conduction band, as clearly shown in Figure S1b.This enhanced EBO brings additional obstacles to the transfer of electrons to the front electrode, ultimately hindering charge transfer with the increase in CBO.Therefore, when CBO decreases below −0.17 eV, the electrons can surmount the potential barrier at the FTO/ETL interface and reach the front electrode, thus avoiding computational convergence issues, as shown in Figure S1c.When χ ETL was tuned from 4.3 eV to 3.98 eV (corresponding to a CBO from −0.49 eV to −0.17 eV), remarkable linear enhancements were observed in the PCE, V oc , and FF, whereas the J sc remained relatively stable, as shown in Figure 3.This behavior originated from the reduction in cliff-like interfacial characteristics between the ETL and absorber, coupled with an enhanced built-in electric field at both the ETL/absorber and FTO/ETL interfaces (Figure S1d), resulting in reduced electron capture at the interface [38].Thus, the Sb 2 (S,Se) 3 solar cell reached its optimal performance at a CBO of −0.17 Certainly, the refinement of the CBO at the ETL/absorber interface significantly enhances the transport of photogenerated electrons.Nevertheless, the employment of an appropriate HTL is of paramount importance in facilitating the migration of photogenerated holes toward the back contact while simultaneously impeding the undesirable reverse flow of electrons.Thus, the valence band offset (VBO) at the interface between the absorber and HTL critically governs the transport dynamics of photogenerated holes toward the back electrode in the Sb2(S,Se)3 absorber layer.The VBO value can be calculated as follows: As shown in Figure 4a, the PCE demonstrates a consistent trend of an initial rise followed by a subsequent decline at all values of NA in the HTL, when the VBO varies from −0.4 eV to 0.2 eV.This VBO variation corresponds to a range of Ev_HTL from 4.84 eV to 5.44 eV, with Ev_absorber fixed at 5.24 eV.Notably, across NA values ranging from 10 16 cm −3 to 10 19 cm −3 , the peak PCE occurs at a VBO of −0.1 eV.In Figure 4b, the FF exhibits a similar trend to that of the PCE.However, Voc remains almost unchanged as the VBO increases to −0.1 eV, whereas Jsc initially remains constant and subsequently rapidly decreases after the VBO exceeds 0 eV (see Figure 4c,d).When the VBO exceeds 0 eV, the PCE, FF, and Jsc of the simulated device undergo rapid deterioration.This significant decline can be explained by the emergence of an energy barrier at the valence band of the absorber/HTL Certainly, the refinement of the CBO at the ETL/absorber interface significantly enhances the transport of photogenerated electrons.Nevertheless, the employment of an appropriate HTL is of paramount importance in facilitating the migration of photogenerated holes toward the back contact while simultaneously impeding the undesirable reverse flow of electrons.Thus, the valence band offset (VBO) at the interface between the absorber and HTL critically governs the transport dynamics of photogenerated holes toward the back electrode in the Sb 2 (S,Se) 3 absorber layer.The VBO value can be calculated as follows: As shown in Figure 4a, the PCE demonstrates a consistent trend of an initial rise followed by a subsequent decline at all values of N A in the HTL, when the VBO varies from −0.4 eV to 0.2 eV.This VBO variation corresponds to a range of E v_HTL from 4.84 eV to 5.44 eV, with E v_absorber fixed at 5.24 eV.Notably, across N A values ranging from 10 16 cm −3 to 10 19 cm −3 , the peak PCE occurs at a VBO of −0.1 eV.In Figure 4b, the FF exhibits a similar trend to that of the PCE.However, V oc remains almost unchanged as the VBO increases to −0.1 eV, whereas J sc initially remains constant and subsequently rapidly decreases after the VBO exceeds 0 eV (see Figure 4c,d).When the VBO exceeds 0 eV, the PCE, FF, and J sc of the simulated device undergo rapid deterioration.This significant decline can be explained by the emergence of an energy barrier at the valence band of the absorber/HTL interface.This barrier considerably impedes the seamless flow of photogenerated holes toward the back electrode, leading to their accumulation at the interface.This accumulation, in turn, intensifies the surface recombination processes, further reducing the performance metrics of the device.When VBO is −0.1 eV (E v_HTL = 5.14 eV) and N A ranges between 10 18 cm −3 and 10 19 cm −3 , the device exhibits optimal performance.As demonstrated by the simulation results, extreme positive or negative values of either the CBO or VBO adversely affect device performance.By optimizing the energy level alignment, we achieved an optimal CBO of −0.17 eV and VBO of −0.1 eV, which were applied in Simulation 3. Simulation 3 achieved a PCE of 15.31%, a Voc of 0.76 V, a Jsc of 24.86 mA•cm −2 , and an FF of 81.32%.Compared with Simulation 2, the device performance in Simulation 3 showed significant improvements, with a significant reduction of 15.19% in Voc loss and a noteworthy 21.80% increase in the PCE, as detailed in Table 1.These results highlight the critical role of precise control over energy level alignment in optimizing photovoltaic device performance.To achieve this optimization while keeping the absorber material's energy levels constant, the electron affinity of the ETL (CdS) needs to be reduced from 4.1 eV to 3.98 eV, and the electron affinity of the HTL (spiro-OMeTAD) needs to be reduced from 2.1 eV to 2.0 eV.In practice, it is necessary to optimize the conduction band energy levels of each functional layer material by adjusting doping elements or introducing specific functional groups to achieve better energy level matching, thereby constructing more efficient solar cell structures.

Optimization of Interface Recombination
Even if optimal band alignment is achieved in the device structure, the potential effect of nonradiative recombination due to interface defects resulting from lattice mismatches between functional layers must be considered.Sb2(S,Se)3 possesses an orthorhombic crystal structure, while CdS typically has a hexagonal structure and spiro-OMeTAD adopts a triclinic structure [39].The inherent structural incompatibility among the ETL, absorber, and HTL layers leads to the emergence of lattice mismatches at their junctions, fostering the development of interface defects that act as sites for nonradiative recombination.To simulate nonradiative recombination through the Shockley-Read-Hall (SRH) mechanism, the SRH recombination rate (RSRH) is considered.The RSRH consists of three main components: bulk recombination (Rbulk) within the absorber layer, interface re- As demonstrated by the simulation results, extreme positive or negative values of either the CBO or VBO adversely affect device performance.By optimizing the energy level alignment, we achieved an optimal CBO of −0.17 eV and VBO of −0.1 eV, which were applied in Simulation 3. Simulation 3 achieved a PCE of 15.31%, a V oc of 0.76 V, a J sc of 24.86 mA•cm −2 , and an FF of 81.32%.Compared with Simulation 2, the device performance in Simulation 3 showed significant improvements, with a significant reduction of 15.19% in V oc loss and a noteworthy 21.80% increase in the PCE, as detailed in Table 1.These results highlight the critical role of precise control over energy level alignment in optimizing photovoltaic device performance.To achieve this optimization while keeping the absorber material's energy levels constant, the electron affinity of the ETL (CdS) needs to be reduced from 4.1 eV to 3.98 eV, and the electron affinity of the HTL (spiro-OMeTAD) needs to be reduced from 2.1 eV to 2.0 eV.In practice, it is necessary to optimize the conduction band energy levels of each functional layer material by adjusting doping elements or introducing specific functional groups to achieve better energy level matching, thereby constructing more efficient solar cell structures.

Optimization of Interface Recombination
Even if optimal band alignment is achieved in the device structure, the potential effect of nonradiative recombination due to interface defects resulting from lattice mismatches between functional layers must be considered.Sb 2 (S,Se) 3 possesses an orthorhombic crystal structure, while CdS typically has a hexagonal structure and spiro-OMeTAD adopts a triclinic structure [39].The inherent structural incompatibility among the ETL, absorber, and HTL layers leads to the emergence of lattice mismatches at their junctions, fostering the development of interface defects that act as sites for nonradiative recombination.To simulate nonradiative recombination through the Shockley-Read-Hall (SRH) mechanism, the SRH recombination rate (R SRH ) is considered.The R SRH consists of three main components: bulk recombination (R bulk ) within the absorber layer, interface recombination at the ETL/absorber interface (R IF1 ), and interface recombination at the absorber/HTL interface (R IF2 ) [39][40][41], as shown in Equation ( 9): where σ n and σ p represent the capture cross-sections for electrons and holes, respectively; υ th denotes the thermal velocity; and n i and E i represent the intrinsic density and energy level of materials, respectively.Additionally, N t and E t are the density and energy level of the defect states, respectively, and S represents the surface recombination velocity (defined as S = σ n,p υ th N t ).
The influence of the surface recombination velocity, namely S 1 at the ETL/absorber interface and S 2 at the absorber/HTL interface, on the performance parameters of Sb 2 (S,Se) 3 solar cells is clearly illustrated in Figure 5a-d.When S 2 is held constant and S 1 varies between 10 6 and 10 10 cm•s −1 , the V oc remains below 0.80 V, with negligible changes in both the PCE and FF.However, when S 2 is less than 10 4 cm•s −1 and S 1 decreases from 10 6 to 10 2 cm•s −1 , the PCE gradually increases, ultimately exceeding 20%, and the V oc progressively increases to nearly 1.0 V.These findings underscore the more significant impact of S 1 than S 2 on the V oc and PCE of Sb 2 (S,Se) 3 solar cells.This result is likely attributed to the presence of a p-n junction formed between the p-type Sb 2 (S,Se) 3 absorber layer and the n-type CdS layer, which governs the behavior of holes as minority carriers within the Sb 2 (S,Se) 3 thin film.In contrast, a relatively weaker p + -p junction is formed with p-type Spiro-OMeTAD.Thus, the optimization of the ETL/absorber interface through meticulous interface engineering can likely lead to a more effective increase in the performance of Sb 2 (S,Se) 3 solar cells [42].where σn and σp represent the capture cross-sections for electrons and holes, respectively; υth denotes the thermal velocity; and ni and Ei represent the intrinsic density and energy level of materials, respectively.Additionally, Nt and Et are the density and energy level of the defect states, respectively, and S represents the surface recombination velocity (defined as S = σn,pυthNt).
The influence of the surface recombination velocity, namely S1 at the ETL/absorber interface and S2 at the absorber/HTL interface, on the performance parameters of Sb2(S,Se)3 solar cells is clearly illustrated in Figure 5a-d.When S2 is held constant and S1 varies between 10 6 and 10 10 cm•s −1 , the Voc remains below 0.80 V, with negligible changes in both the PCE and FF.However, when S2 is less than 10 4 cm•s −1 and S1 decreases from 10 6 to 10 2 cm•s −1 , the PCE gradually increases, ultimately exceeding 20%, and the Voc progressively increases to nearly 1.0 V.These findings underscore the more significant impact of S1 than S2 on the Voc and PCE of Sb2(S,Se)3 solar cells.This result is likely attributed to the presence of a p−n junction formed between the p-type Sb2(S,Se)3 absorber layer and the n-type CdS layer, which governs the behavior of holes as minority carriers within the Sb2(S,Se)3 thin film.In contrast, a relatively weaker p + -p junction is formed with p-type Spiro-OMeTAD.Thus, the optimization of the ETL/absorber interface through meticulous interface engineering can likely lead to a more effective increase in the performance of Sb2(S,Se)3 solar cells [42].Notably, Simulation 4 with optimized parameters of 10 2 cm•s −1 for S1 and 10 3 cm•s −1 for S2 attained an impressive PCE of 23.75%, with excellent characteristics, such as an nID of 1.08, a Voc loss of 0.45 V, and a negligible FF loss of 0.96%, as detailed in Table 1.As evident in Figure 5e, compared with simulation 3, optimizing the interface recombination velocities led to a notable increase in Voc from 0.76 V to 0.98 V, Jsc from 24.86 mA•cm −2 to 27.93 mA•cm −2 , and PCE from 15.31% to 23.75%.Figure 5f shows that this enhancement was attributed to the improved optical response of the device at wavelengths below 600 nm after reducing interface recombination.Reducing the trap density at the ETL/absorber Notably, Simulation 4 with optimized parameters of 10 2 cm•s −1 for S 1 and 10 3 cm•s −1 for S 2 attained an impressive PCE of 23.75%, with excellent characteristics, such as an n ID of 1.08, a V oc loss of 0.45 V, and a negligible FF loss of 0.96%, as detailed in Table 1.As evident in Figure 5e, compared with simulation 3, optimizing the interface recombination velocities led to a notable increase in V oc from 0.76 V to 0.98 V, J sc from 24.86 mA•cm −2 to 27.93 mA•cm −2 , and PCE from 15.31% to 23.75%.Figure 5f shows that this enhancement was attributed to the improved optical response of the device at wavelengths below 600 nm after reducing interface recombination.Reducing the trap density at the ETL/absorber interface from 10 14 cm −2 to 10 10 cm −2 or the capture cross-section from 10 −15 cm 2 to 10 −19 cm 2 is a significant challenge for practical experiments.Indeed, interface modification engineering can passivate surface defects and improve material smoothness and surface impurities, effectively reducing trap density and/or capturing cross-sections.For example, CdS thin films etched with N 2 H 4 exhibit higher transmittance, smoother surfaces, and fewer Cd oxychlorides [12].Additionally, adding insulating polymer layers can further reduce surface roughness and interface states, enhancing device performance [41].

Optimization of the Sb 2 (S,Se) 3 Absorber's Thickness and Bulk Defects
The thickness of the Sb 2 (S,Se) 3 absorber significantly affects various performance aspects of solar cells, including light absorption, carrier mobility, and recombination rates.Consequently, our subsequent simulations delved into the impact of varying the thickness and N t of the absorber to pinpoint the optimal configuration for enhancing device performance.When N t remained at 10 16 cm −3 and the thickness remained below 600 nm, as evident from Figure 6a-d, there was a slight increase in V oc and a slight decrease in FF with increasing thickness, whereas both the PCE and J sc significantly increased.The increase in absorber thickness led to greater photon absorption, resulting in a greater number of photogenerated carriers and thus improving the J sc and PCE.However, once the thickness exceeds 600 nm, the photon absorption approaches saturation and the PCE almost plateaus, as shown in Figure 6e.Furthermore, increasing the absorber thickness from 100 nm to 1200 nm led to a significant increase in the response of the EQE (see Figure 6f).Nevertheless, this increase in absorber thickness extended the transit distance for the charge carriers before collection and elevated the internal resistance of the device, consequently reducing the FF [42].At a thickness of 600 nm, the FF stabilized at nearly 86.01%.Consequently, with an Sb 2 (S,Se) 3 thickness of 600 nm and an initial N t of 6.28 × 10 12 cm −3 , Simulation 5 achieved a PCE of 26.77%, an FF of 87.09%, a V oc of 0.99 V, and a J sc of 30.93 mA•cm −2 , as summarized in Table 1.
Nanomaterials 2024, 14, x FOR PEER REVIEW 11 of 18 impurities, effectively reducing trap density and/or capturing cross-sections.For example, CdS thin films etched with N2H4 exhibit higher transmittance, smoother surfaces, and fewer Cd oxychlorides [12].Additionally, adding insulating polymer layers can further reduce surface roughness and interface states, enhancing device performance [41].

Optimization of the Sb2(S,Se)3 Absorber's Thickness and Bulk Defects
The thickness of the Sb2(S,Se)3 absorber significantly affects various performance aspects of solar cells, including light absorption, carrier mobility, and recombination rates.Consequently, our subsequent simulations delved into the impact of varying the thickness and Nt of the absorber to pinpoint the optimal configuration for enhancing device performance.When Nt remained at 10 16 cm −3 and the thickness remained below 600 nm, as evident from Figure 6a-d, there was a slight increase in Voc and a slight decrease in FF with increasing thickness, whereas both the PCE and Jsc significantly increased.The increase in absorber thickness led to greater photon absorption, resulting in a greater number of photogenerated carriers and thus improving the Jsc and PCE.However, once the thickness exceeds 600 nm, the photon absorption approaches saturation and the PCE almost plateaus, as shown in Figure 6e.Furthermore, increasing the absorber thickness from 100 nm to 1200 nm led to a significant increase in the response of the EQE (see Figure 6f).Nevertheless, this increase in absorber thickness extended the transit distance for the charge carriers before collection and elevated the internal resistance of the device, consequently reducing the FF [42].At a thickness of 600 nm, the FF stabilized at nearly 86.01%.Consequently, with an Sb2(S,Se)3 thickness of 600 nm and an initial Nt of 6.28 × 10 12 cm −3 , Simulation 5 achieved a PCE of 26.77%, an FF of 87.09%, a Voc of 0.99 V, and a Jsc of 30.93 mA•cm −2 , as summarized in Table 1.

Optimization of the Work Function at the Back Contact Layer
The variation in the work function of the back contact (φ BC ) directly impacts the energy band structure at the rear of solar cells, thereby impacting the collection of photogenerated carriers.We conducted a study to examine the influence of φ BC on the performance of Sb 2 (S,Se) 3 solar cells, with a specific focus on the difference ∆φ (referred to as E v_HTL − φ BC ).When ∆φ ranges from −0.44 eV to 0.46 eV, φ BC shifts within a range of 4.74 eV to 5.60 eV. Figure 7a reveals a consistent increasing trend in both the PCE and FF, reaching their maximum values at −0.14 eV.Beyond this ∆φ value, the PCE and FF remain stable.The downward curvature in the valence band and the characteristic S-shaped curve in the current-voltage characteristics when ∆φ falls below −0.14 eV, as depicted in Figure 7b,c, clarify this phenomenon.The varying energy levels across the HTL/back contact interface, causing the emergence of a Schottky barrier within the ∆φ range of −0.44 eV to −0.14 eV, are responsible for impeding the unhampered flow of holes from the HTL toward the back electrode, as shown in Figure 7b.Notably, when ∆φ approaches or exceeds −0.14 eV, the barrier gradually diminishes, allowing for improved hole transport and reduced resistance at the HTL/back contact interface, and all device parameters are saturated.In other words, when φ BC exceeds 5.04 eV, the back electrode and HTL form a good ohmic contact conducive to hole collection.The back contact employed in Simulation 5 is gold (φ BC = 5.1 eV), which falls within the optimized range, and the back contact of the Sb 2 (S,Se) 3 solar cell can be unchanged.

Optimization of the Work Function at the Back Contact Layer
The variation in the work function of the back contact (φBC) directly impacts the energy band structure at the rear of solar cells, thereby impacting the collection of photogenerated carriers.We conducted a study to examine the influence of φBC on the performance of Sb2(S,Se)3 solar cells, with a specific focus on the difference Δφ (referred to as Ev_HTL − φBC).When Δφ ranges from −0.44 eV to 0.46 eV, φBC shifts within a range of 4.74 eV to 5.60 eV. Figure 7a

Final Optimized Sb2(S,Se)3 Solar Cell Performance
After the above optimization analysis, the final values are as follow: the thickness of Sb2(S,Se)3 is 600 nm, and the electron affinities of the ETL (CdS) and HTL (spiro-OMeTAD) are optimized to 3.98 eV and 2.04 eV, respectively; the interface recombination velocities of ETL/absorber and absorber/HTL are optimized at 10 2 cm•s −1 and 10 3 cm•s −1 , respectively.Additionally, the back electrode with a work function greater than 5.04 eV (e.g., Au~5.1 eV) can form a better ohmic contact with HTL, and the series and shunt resistances of the device are optimized to 0 and 10 5 Ω•cm 2 , respectively.As shown in Table 1, the final optimized device (simulation 5) has achieved excellent performance parameters: Voc = 0.99 V, Jsc = 30.93mA•cm −2 , FF = 87.09%,and PCE = 26.77%.Compared with other simulation results presented in Table S4, the performance of our simulated Sb2(S,Se)3 solar cell is superior, due to the effective optimization path chosen based on Voc and FF losses to enhance the device structure.

Improvement in the PCE of the Sb2(S,Se)3 Solar Cells
Based on the preceding simulation outcomes, a two-phased approach is proposed to increase the efficiency of Sb2(S,Se)3 solar cells.Initially, an efficiency benchmark of up to 15% was achieved, and subsequently, the goal was to surpass the 20% efficiency threshold.

Final Optimized Sb 2 (S,Se) 3 Solar Cell Performance
After the above optimization analysis, the final values are as follow: the thickness of Sb 2 (S,Se) 3 is 600 nm, and the electron affinities of the ETL (CdS) and HTL (spiro-OMeTAD) are optimized to 3.98 eV and 2.04 eV, respectively; the interface recombination velocities of ETL/absorber and absorber/HTL are optimized at 10 2 cm•s −1 and 10 3 cm•s −1 , respectively.Additionally, the back electrode with a work function greater than 5.04 eV (e.g., Au~5.1 eV) can form a better ohmic contact with HTL, and the series and shunt resistances of the device are optimized to 0 and 10 5 Ω•cm 2 , respectively.As shown in Table 1, the final optimized device (simulation 5) has achieved excellent performance parameters: V oc = 0.99 V, J sc = 30.93mA•cm −2 , FF = 87.09%,and PCE = 26.77%.Compared with other simulation results presented in Table S4, the performance of our simulated Sb 2 (S,Se) 3 solar cell is superior, due to the effective optimization path chosen based on V oc and FF losses to enhance the device structure.

Improvement in the PCE of the Sb 2 (S,Se) 3 Solar Cells
Based on the preceding simulation outcomes, a two-phased approach is proposed to increase the efficiency of Sb 2 (S,Se) 3 solar cells.Initially, an efficiency benchmark of up to 15% was achieved, and subsequently, the goal was to surpass the 20% efficiency threshold.
Approaching 15% efficiency.Although the lowest V oc loss reported in Sb 2 (S,Se) 3 solar cells has an E g of 1.46 eV [3], the V oc value remains at approximately 0.69 V, which is notably far lower than the S-Q limit V oc of 1.19 V.This discrepancy is commonly observed in Sb 2 (S,Se) 3 solar cells with PCEs surpassing 10% [3,4,11,[13][14][15][43][44][45], highlighting a persistent challenge in their performance improvement [26].To overcome this challenge, feasible methods for enhancing carrier transport properties, minimizing recombination losses, and improving the internal resistance of Sb 2 (S,Se) 3 solar cells include the following: (1) Regulation of the absorption film to reduce the internal resistance.The Sb 2 (S,Se) 3 compound possesses a unique quasi-1D ribbon architecture characterized by [Sb 4 S(e) 6 ] n ribbons [46].This structural peculiarity engenders favorable grain morphology and anisotropic carrier mobility, most notably along the [hk1] direction.This specific transport pathway underscores the prospect of significantly elevated device performance, particularly in comparison with transport along the [hk0] direction.Different synthesis/fabrication strategies have been used to develop [hk1]-oriented Sb 2 (S,Se) 3 films.
A highly efficient method lies in rapid hydrothermal deposition, which kinetically promotes the growth of [221]-oriented crystals [15].Additionally, during synthesis processes such as close space sublimation (CSS) [47], rapid thermal evaporation (RTE) [45], and vapor transport deposition (VTD) [48], carefully balancing the temperature between the evaporation source and the substrate has been shown to be crucial for promoting [hk1] orientation and enhancing the crystallinity of Sb 2 (S,Se) 3 films.High-quality crystals, characterized by low defect density, extended carrier lifetimes, elevated electron and hole mobilities, and uniform crystal structures, play a pivotal role in reducing internal resistance and optimizing carrier transport in solar cells.(2) Band alignment absorber layer and charge-transport layers.To fully unlock the potential of Sb 2 (S,Se) 3 solar cells in terms of the PCE, the optimization of the carrier transport path has emerged as a paramount step forward in advancing their performance.This can be achieved through careful adjustment of CBO at the ETL/absorber interface and VBO at the absorber/HTL interfaces [49].The proper alignment of the energy bands between the absorber layer and charge-transport layers is essential for minimizing energy-level mismatches and reducing the interface energy barriers.
Optimizing the interface energy barrier involves using ETLs and HTLs with minimal defects or employing doping strategies to mitigate band mismatch.For example, nitrogen-containing functional groups in ethylenediamine (EDA) can coordinate with CdS, and shift the Fermi level of EDA-CdS toward the conduction band [50].Alternatively, in situ oxygen doping of CdS films results in a lower conduction band level compared to that of control CdS films [51].Various configurations of CdS-based ETLs such as SnO 2 /CdS [25,52], TiO 2 /CdS [53,54], ZnO/CdS [55], and Zn(O,S)/CdS [56] ensure well-matched band alignment, thereby facilitating efficient charge transport in Sb 2 (S,Se) 3 solar cells.
Exceeding 20% efficiency.Reducing SRH recombination is a critical aspect of improving the performance of Sb 2 (S,Se) 3 solar cells.By implementing appropriate passivation strategies and surface treatments, such as postselenization [57,58] and advanced interfacial engineering [59,60], it is feasible to significantly diminish the defect density within a material.This reduction not only mitigates the entrapment of charge carriers but also increases their lifetime.Achieving an efficiency of more than 20% in Sb 2 (S,Se) 3 solar cells requires suitable methods, including the following: (1) Regulation of the interface recombination.Annealing and interface treatments are essential for rectifying instability, smoothing out roughness, and enhancing conductivity within ETLs.This helps minimize the defect density and lower carrier recombination at the ETL/absorber junction.For example, oxygen-doped cadmium sulfide (CdS:O) facilitates the tailoring of (Sb 4 S(e) 6 ) n ribbons to the (221)-textured orientation [61], CdCl 2 -modified CdS films passivate surface defects [62], and KCl-treated CdS films promote absorption layer growth [10].Furthermore, enhancing the p-type characteristics of HTL materials through doping or developing low-cost and efficient alternatives is crucial for increasing hole mobility and improving material conductivity, such as high-mobility and work-function NiO x [63,64], excellent near-infrared-absorption PbS colloidal quantum dots [7,65], and high-hole-mobility and high-stability CuSCN [66].
These alternative approaches hold significant potential in increasing the performance of solar cells, primarily through the enhancement of hole mobility and the mitigation of interface recombination.(2) Balance of the relationship between the thickness of the absorber and bulk defects.Both theoretical calculations [67][68][69] and experimental measurements [3,70] have demonstrated the complex defect mechanisms in Sb 2 (S,Se) 3 .The thickness of the absorber layer influences both photon absorption and carrier transport [71].If the absorber layer is too thin, insufficient photon absorption occurs, and pinhole formation increases, thereby resulting in a low J sc .However, during the actual preparation process, overly thick layers increase photon absorption while also increasing the series resistance and carrier transport distance, thereby intensifying carrier recombination.The thickness of the absorber layer can be controlled by adjusting the synthesis/fabrication time and growth rate.Furthermore, the absorber bulk defects are reduced through empirical methods such as adjusting the Se/S atomic ratio [15,72], postselenization treatments [73], seeding materials [74,75], or additive engineering [4,76].

Conclusions
In summary, a numerical model of Sb 2 (S,Se) 3 solar cells was developed utilizing the SCAPS software (version 3.3.10;University of Gent, Belgium) package and validated against experimental data from the literature.This model was used to investigate the factors limiting V oc and FF at various stages of Sb 2 (S,Se) 3 solar cell development, and effective strategies were identified to increase the PCE of these cells.The simulation results indicated an improvement in the PCE from 10.72% to 12.57%, as well as a significant increase in the FF from 68.35% to 79.55% after addressing the limiting mechanisms related to R s and R sh .Further adjustments in the CBO at the ETL/absorber and VBO at the absorber/HTL interfaces led to an increase in the PCE from 12.57% to 15.31%, along with a notable 18.75% improvement in V oc .These results highlight the importance of reducing the internal resistance to enhance the FF and show the importance of optimal energy level arrangements for increasing V oc .
The subsequent optimization of interface recombination led to a notable increase in the overall performance of Sb 2 (S,Se) 3 solar cells, resulting in a remarkable PCE exceeding 20%.This was accompanied by a decrease in the V oc loss to 0.45 V, an increase in the J sc to 27.93 mA•cm −2 , and a reduction in the FF loss to 0.96%.These results highlight the crucial role of suppressing interface SRH recombination, particularly at the ETL/absorber interface, in improving the overall performance of devices.By striking a balance between the absorber thickness and the bulk defects, the J sc increased to over 30 mA•cm −2 at a thickness of 600 nm, consequently increasing the PCE to an impressive 26.77%.
This numerical analysis proposes an optimization roadmap aimed at reducing V oc and FF losses, thereby successfully enhancing the efficiency of Sb 2 (S,Se) 3 solar cells.While numerical simulations cannot completely replace experimental work, they undoubtedly provide significant guidance for future experimental efforts to improve the efficiency of Sb 2 (S,Se) 3 solar cells, helping to avoid blind experimental exploration and thus saving valuable time and resources.Furthermore, this study proposes a two-stage strategy based on simulation results to further boost the efficiency of Sb 2 (S,Se) 3 solar cells.We hope these findings can also be extended to other solar materials, such as Sb 2 Se 3 and Sb 2 S 3 , thereby broadening their applicability beyond current boundaries.

Supplementary Materials:
The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/nano14171433/s1, Figure S1 S1: Simulation parameters applied to different layers within the solar cell model; Table S2: Simulation pa-

Figure 1 .
Figure 1.(a) Schematic illustrations depicting the simulated device.(b) J−V curves, (c) EQE spectra, and (d) the ideal factor (nID) comparison between the experimental and simulated devices.(e) Energy band structure of the Sb2(S,Se)3 solar cell.(f) The FFm (solid lines) and FF (dotted lines) versus Voc for different nID values.Experimental data from the literature [4].

Figure 1 .
Figure 1.(a) Schematic illustrations depicting the simulated device.(b) J-V curves, (c) EQE spectra, and (d) the ideal factor (n ID ) comparison between the experimental and simulated devices.(e) Energy band structure of the Sb 2 (S,Se) 3 solar cell.(f) The FF m (solid lines) and FF (dotted lines) versus V oc for different n ID values.Experimental data from the literature [4].

Figure 2 .
Figure 2. Simulated results for (a) PCE, (b) FF, (c) V oc , (d) J sc , and (e) n ID with varying R s and R sh .(f) The J-V curves for Simulation 1 versus Simulation 2.

Figure 5 .
Figure 5. Simulated results for (a) PCE, (b) Voc, (c) Jsc, and (d) FF as a function of varied S1 and S2.(e) J-V curves and (f) EQE spectra for Simulation 3 and Simulation 4.

Figure 5 .
Figure 5. Simulated results for (a) PCE, (b) V oc , (c) J sc , and (d) FF as a function of varied S 1 and S 2 .(e) J-V curves and (f) EQE spectra for Simulation 3 and Simulation 4.

Figure 6 .
Figure 6.Simulated results for (a) FF, (b) PCE, (c) Voc, and (d) Jsc with varied Nt and thickness of the Sb2(S,Se)3 layer.(e) J-V curves and (f) EQE spectra with varied thickness of the Sb2(S,Se)3 layer.

Figure 6 .
Figure 6.Simulated results for (a) FF, (b) PCE, (c) V oc , and (d) J sc with varied N t and thickness of the Sb 2 (S,Se) 3 layer.(e) J-V curves and (f) EQE spectra with varied thickness of the Sb 2 (S,Se) 3 layer.
reveals a consistent increasing trend in both the PCE and FF, reaching their maximum values at −0.14 eV.Beyond this Δφ value, the PCE and FF remain stable.The downward curvature in the valence band and the characteristic S-shaped curve in the current−voltage characteristics when Δφ falls below −0.14 eV, as depicted in Figure7b,c, clarify this phenomenon.The varying energy levels across the HTL/back contact interface, causing the emergence of a Schottky barrier within the Δφ range of −0.44 eV to −0.14 eV, are responsible for impeding the unhampered flow of holes from the HTL toward the back electrode, as shown in Figure7b.Notably, when Δφ approaches or exceeds −0.14 eV, the barrier gradually diminishes, allowing for improved hole transport and reduced resistance at the HTL/back contact interface, and all device parameters are saturated.In other words, when φBC exceeds 5.04 eV, the back electrode and HTL form a good ohmic contact conducive to hole collection.The back contact employed in Simulation 5 is gold (φBC = 5.1 eV), which falls within the optimized range, and the back contact of the Sb2(S,Se)3 solar cell can be unchanged.

Figure 7 .
Figure 7. (a) The influence of Δφ on device performance, (b) the energy band diagrams, and (c) J-V curses with varied Δφ from −0.44 eV to 0.16 eV.

Figure 7 .
Figure 7. (a) The influence of ∆φ on device performance, (b) the energy band diagrams, and (c) J-V curses with varied ∆φ from −0.44 eV to 0.16 eV.
: (a) Energy band diagram of Sb 2 (S,Se) 3 solar cells under different χ ETL , (b) The conduction band diagram of the dashed circle in (a), (c) The conduction band diagram at CBO = −0.17eV, (d) The variation of built-in electric field at ETL/absorber and FTO/ETL interfaces with CBO varied from −0.49 eV to −0.17 eV; Table

Table 1 .
Photovoltaic parameters in five simulated devices.

Table 1 .
Photovoltaic parameters in five simulated devices.