Abstract
In crystalline semiconductors, absorption onset sharpness is characterized by temperature-dependent Urbach energies. These energies quantify the static, structural disorder causing localized exponential-tail states, and dynamic disorder from electron-phonon scattering. Applicability of this exponential-tail model to disordered solids has been long debated. Nonetheless, exponential fittings are routinely applied to sub-gap absorption analysis of organic semiconductors. Herein, we elucidate the sub-gap spectral line-shapes of organic semiconductors and their blends by temperature-dependent quantum efficiency measurements. We find that sub-gap absorption due to singlet excitons is universally dominated by thermal broadening at low photon energies and the associated Urbach energy equals the thermal energy, regardless of static disorder. This is consistent with absorptions obtained from a convolution of Gaussian density of excitonic states weighted by Boltzmann-like thermally activated optical transitions. A simple model is presented that explains absorption line-shapes of disordered systems, and we also provide a strategy to determine the excitonic disorder energy. Our findings elaborate the meaning of the Urbach energy in molecular solids and relate the photo-physics to static disorder, crucial for optimizing organic solar cells for which we present a revisited radiative open-circuit voltage limit.
Similar content being viewed by others
Introduction
Recently, research on organic solar cells has seen significant progress through the development of non-fullerene electron acceptors (NFAs), in particular delivering increases in single-junction efficiencies from 13.1 to 18.7 % between 2017 to 20211,2. This has revived ambitions to ultimately achieve industrial scale low-cost photovoltaics with low embodied manufacturing energy. While chemists have been productively synthesizing new materials, understanding of the opto-electronic properties lags behind â particularly in relation as to why these new NFA blends are so effective at photogeneration with low voltage losses. A particular area of intense interest is light absorption and how it is related to molecular dipole moments and their energetic disorder. Of course, this is also a very relevant question for the fundamental solid-state physics of molecular and disordered semiconductors. Furthermore, energetic disorder critically determines dominant radiative loss mechanisms limiting the open-circuit voltage in optoelectronic applications.
In general, semiconductors tend to absorb light at photon energies below the bandgap (sub-gap absorption) depending on the energetic disorder. In the case of inorganic semiconductors, the absorption coefficient \((\alpha )\) often displays an exponential tail below the bands. These so-called Urbach tails increase their broadening with temperature \(T\)3. The sub-gap \(\alpha\) generally follows the expression
where \(E\) is the photon energy, \({E}_{{\rm{on}}}\) is the energy onset of the tail, and \({E}_{{\rm{U}}}(T)\) is the Urbach energy quantifying the total energetic disorder of the system. Depending on the semiconductor, \({E}_{{\rm{U}}}\) generally varies between 10 to 100âmeV at room temperature4,5,6,7. In banded semiconductors, it has been suggested that \({E}_{{\rm{U}}}\left(T\right)={E}_{{\rm{U}},{\rm{D}}}\left(T\right)+{E}_{{\rm{U}},{\rm{S}}}\), where \({E}_{{\rm{U}},{\rm{D}}}\left(T\right)\) is a temperature-dependent dynamical disorder term related to the thermal occupation of phonon states4,8, while \({E}_{{\rm{U}},{\rm{S}}}\) is the width of the assumed exponential distribution of sub-gap states induced by static disorder9. However, a unifying theory describing the density of states and their absorption leading to Urbach tails for materials of different chemical bonding and morphology is still lacking10,11,12.
In non-excitonic amorphous semiconductors, the definition of a clear bandgap edge is often difficult due to large static disorder inducing sub-gap broadening. Organic semiconductors, which are excitonic and partially amorphous, display even more complex sub-gap features including intermolecular hybrid charge transfer (CT) states in technologically-relevant blends of electron donors (D) and acceptors (A), excitonic features13 and trap states14. CT states typically give rise to light absorption with Gaussian sub-gap spectral line-shapes. This has been attributed to intermolecular D:A transitions described by non-adiabatic Marcus theory15 or its extensions16,17,18,19, suggesting that the static disorder is governed by a Gaussian distribution of CT states. Whether a similar description applies for sub-gap absorption by intramolecular excitons is unclear. With the rise of NFA semiconductors with small energetic offset relative to donors20, the lack of CT state sub-gap spectral features in donor-acceptor blends has led to the increased usage of Urbach energies to understand disorder and sub-gap photo-physics. As a rule of thumb, it has been assumed that low Urbach energies are indicative of lower static disorder and hence expected to result in higher performance such as reduced charge recombination and higher charge carrier mobilities in a photovoltaic device21,22,23. Due to the complex and often convoluted spectral features in sub-gap light absorption of organic semiconductors, however, the Urbach energy carries significant ambiguity. Moreover, a long-standing debate24,25,26 on the distribution of density of states (DOS) defining the static disorder in organic semiconductors has been revived: does the DOS follow a Gaussian or an exponential distribution? With no doubt, the origin of the static disorder is of great importance for the classification and future development of organic semiconductors and their electro-optical properties. However, it has remained unclear how this important Figure-of-merit relates to the Urbach energy.
In this work, we show that the exciton sub-gap absorption in organic semiconductors, at photon energies well below the gap, is generally characterized by Urbach tails with characteristic energies equivalent to the thermal energy \({kT}\) as demonstrated by temperature-dependent external quantum efficiency (EQE) measurements (k is the Boltzmann constant). However, these Urbach tails are often convoluted with Gaussian line-shapes induced by trap states and/or CT states, resulting in erroneous, energy-dependent, Urbach energies larger than \({kT}.\) A simple model, combining the Gaussian distributed DOS with Boltzmann-like thermally activated optical transitions, is shown to reproduce the absorption profiles and to give an estimate for the Gaussian static disorder. Based on this model, the exciton sub-gap absorption is found to be composed of two regimes: near the onset, the sub-gap absorption is dominated by Gaussian static disorder resulting in strongly energy dependent Urbach energies, while thermal broadening dominates at energies well below the gap, where the Urbach energy approaches \({kT}\). As such, using an (energy-independent) Urbach energy as a probe for the static disorder in organic solar cells is meaningless. Finally, the effect of exciton static disorder and thermal broadening on the expected radiative voltage losses and power conversion efficiency (PCE) in organic solar cells based on low offset D:A blends is clarified.
Results
The spectral line-shape of \(\alpha\) and the absorptance \(A\) in the sub-gap tail are generally related via a modified Beer-Lambert law, \(A=\tilde{f}\alpha d\), where d is the thickness of the active layer and \(\tilde{f}\) is an energy-dependent correction factor accounting for optical interference27,28. For optically thin films with layer thicknesses of 100 to 150ânm, \(\tilde{f}\) is often assumed to be close to 1 (negligible optical interference effects). Moreover, it has been shown for efficient D:A systems that the internal quantum efficiency (IQE) is generally excitation energy independent, hence EQE \(\propto A\) and the spectral line shape of the EQE follows \(\alpha\)29,30,31,32. Based on this underlying premise, EQE measurements have been frequently employed in the past to determine \({E}_{{\rm{U}}}\). However, a previous lack of sensitivity in the EQE measurements has led to speculative assumptions about the spectral range of trap state absorptions, exponential tails and associated \({E}_{{\rm{U}}}\) in organic semiconductors. By choosing a small fitting range, exponential fits can be forced on to EQE spectra resulting in a rather arbitrary \({E}_{{\rm{U}}}\) dependent on the spectral range of the fitting. More insight can instead be gained from the apparent Urbach energy (\({E}_{{\rm{U}}}^{{\rm{app}}}\)) here defined as:
For a true exponential tail in the form of Eq. 1, \({E}_{{\rm{U}}}^{{\rm{app}}}\) is constant in the sub-gap spectral region and given by \({E}_{{\rm{U}}}\).
Determination of the apparent Urbach energy in organic solar cells
Figure 1 shows the EQE and \({E}_{{\rm{U}}}^{{\rm{app}}}\) for a wide range of organic semiconductor D:A blends. The material systems studied here can be generally grouped according to the D:A energy offset or, in other words, the difference between the CT state energy \({E}_{{\rm{CT}}}\) and the optical bandgap \({E}_{{\rm{opt}}}\). Here, \({E}_{{\rm{opt}}}\,\)is typically equal to the local exciton (LE) energy of the lower-bandgap component (between D and A). Large offset systems, such as PCDTBT:PC70BM, BQR:PC70BM and PBDB-T:PC70BM (Fig. 1a), show three distinct spectral ranges for energies below the gap: (i) mid-gap trap state absorption at energies below 1.2âeV, (ii) CT state absorption in the mid-range (\(\sim\)1.2â1.5âeV) and (iii) LE absorption above 1.5âeV. For mid-gap and CT state absorption, Gaussian-shaped EQE features are observed14. These features are consistent with Marcus charge-transfer between states that are distributed in accordance with a Gaussian density of sub-gap states. Details of the Gaussian fits are provided in the Methods section. Because of the Gaussian line-shape (see Equation 5 in the Methods), the corresponding \({E}_{{\rm{U}}}^{{\rm{app}}}\) in the spectral range of CT absorption is expected to be strongly energy-dependent and given by \({E}_{{\rm{U}}}^{{\rm{app}}}\approx 2{kT}{\lambda }_{{\rm{CT}}}/({E}_{{\rm{CT}}}+{\lambda }_{{\rm{CT}}}-E)\); this is highlighted by the black solid line in Fig. 1a. Here, \({E}_{{\rm{CT}}}\) and \({\lambda }_{{\rm{CT}}}\), as obtained from the Gaussian fits, are to be considered the effective energy and reorganization energy of CT states which includes the effect of static disorder (see Supplementary Note 2)16,18. At higher energies corresponding to the sub-gap LE absorption regime, in turn, \({E}_{{\rm{U}}}^{{\rm{app}}}\) has a narrow parabolic shape with a sharp minimum at roughly 40 â 50âmeV for the large offset blends.
The EQE and \({E}_{{\rm{U}}}^{{\rm{app}}}\) of blends with a smaller offset between \({E}_{{\rm{CT}}}\) and \({E}_{{\rm{opt}}}\) (PTB7-Th:ITIC, PBDB-T:ITIC and PBDB-T:IT-4F) are shown in the Fig. 1b. In this case, the Gaussian CT state line-shape is barely recognizable, and the \(\alpha\) tail in the spectral range of LEs is visible over a wider spectral range. While \({E}_{{\rm{U}}}^{{\rm{app}}}\) retains its parabolic shape in the LE sub-gap absorption regime, because of the wider range, the minimum value of \({E}_{{\rm{U}}}^{{\rm{app}}}\) is significantly reduced to values close to \({kT}\). Finally, the EQE and \({E}_{{\rm{U}}}^{{\rm{app}}}\) of blends in low-offset D:NFA systems (PM6:Y6, PM6:ITIC and PBDB-T:EH-IDTBR) are shown in Fig. 1c. In these blends, the CT state absorption can no longer be discerned from the EQE spectra. Instead, the \(\alpha\) tail in the spectral range of LEs remains dominant down to the energy range of mid-gap state excitation. In this limit, the \({E}_{{\rm{U}}}^{{\rm{app}}}\) in the LE-dominated sub-gap absorption range finally saturates and reaches a broad plateau where \({E}_{{\rm{U}}}^{{\rm{app}}}\approx {kT}\).
The above experimental observation for low offset D:NFA systems suggests the presence of an intramolecular sub-gap absorption regime where \(\alpha \propto {\rm{exp }}\left(E/{kT}\right),\) but that its spectroscopic observation is obstructed by CT absorption in systems with larger offsets. To further clarify this, neat D or A systems (lacking the D:A CT absorption feature in the sub-gap EQE spectrum) were investigated as shown in Fig. 2. In the case of neat PC70BM (Fig. 2a), a narrow spectral range where \({E}_{{\rm{U}}}^{{\rm{app}}}\approx {kT}\) can be identified. We note that this is consistent with photothermal deflection spectroscopy results33 of neat PC60BM (see Supplementary Fig. 3), confirming the underlying assumption that the spectral line shape of the EQE follows α in the sub-gap region. In PC70BM, the spectral range where \({E}_{{\rm{U}}}^{{\rm{app}}}\approx kT\) is limited by deep trap-state absorption at low energies and low signal-to-noise ratio (SNR) due to the poor exciton dissociation in the neat phase (translating into a low Internal Quantum Efficiency [IQE]). However, adding 0.1âmol% of the wide-gap donor m-MTDATA results in enhanced exciton dissociation, which increases the spectral range within which \({E}_{{\rm{U}}}^{{\rm{app}}}\approx kT\) can be observed. Since m-MTDATA:PC70BM is characterized by a low ECT, the LE tail of PC70BM can be clearly distinguished from CT states as shown in Fig. 2a34. By further increasing the donor content, the CT state absorption increases and the parabolic shape of \({E}_{{\rm{U}}}^{{\rm{app}}}\) emerges. This explains the parabolic shapes seen for large offset systems in Fig. 1a, which appear when the \(\alpha\) tail in the spectral range of LE absorption becomes convoluted with CT state absorption. This is further supported by results shown in Fig. 2b for neat NFA devices comprising ITIC and IT-4F active layers, respectively. Here again, exponentially decaying \(\alpha\) tail regions with \({E}_{{\rm{U}}}^{{\rm{app}}}\approx {kT}\) are observed.
To check if \({E}_{{\rm{U}}}^{{\rm{app}}}\approx {kT}\) is also true for donor materials has proven more challenging, since the neat donor absorption is often convoluted with deep trap-state absorption. This convolution effectively increases \({E}_{{\rm{U}}}^{{\rm{app}}}\) above \({kT}\) as shown for PBDB-T in Fig. 2c, where \({E}_{{\rm{U}}}^{{\rm{app}}}\approx\) 30âmeV. The problem can be circumvented by using narrow-gap polymer donors such as PBTQ(OD) or PTTBAI, which show strongly redshifted EQE spectra when blended with PC70BM as previously reported35. As shown in Fig. 2c, \({E}_{{\rm{U}}}^{{\rm{app}}}\) is around 22âmeV in this case.
While deviations of the type \({E}_{{\rm{U}}}^{{\rm{app}}} {\,}> {\,}{kT}+\) 10âmeV can be attributed to the presence of other absorbing species like CT states or deep trap states, small deviations of \({E}_{{\rm{U}}}^{{\rm{app}}}\) around \({kT}\) are likely caused by interference effects, i.e. non-constant \(\widetilde{f}\), as shown previously27,36. Interference effects arise from the spectral dependence of the optical constants of the different layers in the thin-film stack and generally vary with the thickness of the active layer. To demonstrate the presence of interference effects, we therefore fabricated PM6:ITIC devices and PM6:Y6 devices of different active layer thicknesses. Indeed, as illustrated in the Supplementary Fig. 4, \({E}_{{\rm{U}}}^{{\rm{app}}}\) shows a weak thickness dependence leading to deviations around \({kT}\) from \(+\)2.3âmeV to \(-\)5.8âmeV. This is further corroborated by optical transfer-matrix simulations (see Supplementary Fig. 5).
Temperature-dependent EQE measurements
To verify that the exponential \(\alpha\) tails in the spectral range of LEs are characterized by \({E}_{{\rm{U}}}^{{\rm{app}}}\approx {kT}\), we performed \(T\)-dependent EQE measurements on PBDB-T:EH-IDTBR, PM6:Y6 and PM6:ITIC, and the neat materials IT-4F and Y6. Figure 3a shows the normalized EQE and the respective \({E}_{{\rm{U}}}^{{\rm{app}}}\) spectra of three representative systems: PBDB-T:EH-IDTBR, PM6:Y6 and IT-4F. The remaining material systems are provided in the Supplementary Note 8. Depending on \(T,\) two different regimes can be distinguished in the sub-gap absorption tail: At high \(T,\) \({E}_{{\rm{U}}}^{{\rm{app}}}\)spectra show a \(T\)-dependent plateau only influenced by interference effects and the shift in the absorption onset. At low \(T,\) however, the \({E}_{{\rm{U}}}^{{\rm{app}}}\)spectra generally attain a parabolic shape suggesting that the LE tail becomes convoluted with CT and/or mid-gap states. The thermal activation of the exponential \(\alpha\) tail is then finally illustrated in Fig. 3b where \({E}_{{\rm{U}}}^{{\rm{app}}}\) at a constant energy in the plateau region is shown as a function of \({kT}\) for the systems studied. We see that \({E}_{{\rm{U}}}^{{\rm{app}}}\) at higher temperature is linear and equals \({kT}\pm 2.5\) meV, where the offset arises from interference effects. At lower \(T,\) \({E}_{{\rm{U}}}^{{\rm{app}}}\) eventually deviates from linearity, as the spectral shape is increasingly affected by other absorbing species. For example, the low-energy tail of the CT absorption may emerge at low \(T\) owing to its distinctly different \(T\)-dependence compared to LE absorption37. Moreover, the absorption of trap states is expected to play a role as well, while little is yet known about their spectral broadening as a function of \(T\).
For comparison, we also measured \(T\)-dependent EQE spectra of a commercial a-Si:H thin-film solar cell (EQE and \({E}_{{\rm{U}}}^{{\rm{app}}}\) spectra shown in Supplementary Fig. 7). In banded semiconductors such as a-Si:H, the total energetic disorder is \({E}_{{\rm{U}}}\left(T\right)={E}_{{\rm{U}},{\rm{D}}}\left(T\right)+{E}_{{\rm{U}},{\rm{S}}}\), consistent with the presence of an exponential DOS of tail states (defined by the associated width \({E}_{{\rm{U}},{\rm{S}}}\)). In Fig. 3b, two regimes can be distinguished for a-Si:H comprising the low-temperature saturation of \({E}_{{\rm{U}}}\) to \({E}_{{\rm{U}}}\left(0\right)\approx\) 40âmeV and thermal activation at higher temperatures with an offset of roughly 21âmeV from \({kT}\) in agreement with previous reports38,39,40. In contrast, for the organic semiconductor systems in Fig. 3b, an extrapolation to \(T=0\) implies that \({E}_{{\rm{U}},{\rm{S}}}\approx 0\), suggesting the lack of an exponential tail state distribution in these systems. These observations raise the question: what is the role of static disorder in shaping the \({E}_{{\rm{U}}}^{{\rm{app}}}\) spectra?
Model for understanding the sub-gap α
In the context of the Marcus formalism for non-adiabatic charge transfer (high-temperature limit), the presence of \({E}_{{\rm{U}}}^{{\rm{app}}}\approx {kT}\) in the spectral range of LE absorption (described by the absorption coefficient \({\alpha }_{{\rm{LE}}}\)) may be rationalized in terms of non-equal potential energy surfaces for the excited state and the ground state. By assuming a significantly more diffuse or delocalized excited state, compared to the strongly localized ground state, a much smaller reorganization energy is expected for the excited state in comparison to the ground state41,42,43. In this limit, we expect \({\alpha }_{{\rm{LE}}}\left(E\right)\approx {\alpha }_{{\rm{sat}}}{\rm{exp }}\left(\left[E-{E}_{{\rm{opt}}}\right]/{kT}\right)\) for \(E \,<\, {E}_{{\rm{opt}}},\) in the case of negligible static disorder (see Supplementary Note 4). Here, \({\alpha }_{{\rm{sat}}}\) includes a \(1/E\) dependence, however, the sub-gap spectral line-shape is dominated by the Boltzmann factor. For above-gap absorption \(E \,> \, {E}_{{\rm{opt}}}\), we assume \({\alpha }_{{\rm{LE}}}\left(E\right)\approx {\alpha }_{{\rm{sat}}}\). We note that this simplification of \({\alpha }_{{\rm{LE}}}\left(E\right)\) essentially follows a Miller-Abrahams-type charge-transfer formalism, sometimes used to describe exciton migration44, but more commonly used for charge transport45 in organic semiconductors. Accounting for static disorder described by a DOS given by \({g}_{{\rm{DOS}}}\left({E}_{{\rm{opt}}}^{{\prime} }\right)\), the total absorption coefficient is of the form \({\alpha }_{{\rm{LE}}}\left(E\right)=\int {\alpha }_{{\rm{LE}}}\left(E,{E}_{{\rm{opt}}}^{{\prime} }\right){g}_{{\rm{DOS}}}\left({E}_{{\rm{opt}}}^{{\prime} }\right)d{E}_{{\rm{opt}}}^{{\prime} }\). In the Supplementary Note 3, expressions for the associated absorption coefficients are derived for the cases of a Gaussian DOS and an exponential tail DOS. For the case of an exponential tail DOS with the width \(W\), we obtain \({E}_{{\rm{U}}}^{{\rm{app}}}=W\) (see Supplementary Fig. 8). This is indeed not what we observe in experiments. For a Gaussian DOS, in turn, we obtain \(\alpha\) in the form
where \({\sigma}_{s}\) is the standard deviation of the Gaussian DOS and \({E}_{{\rm{opt}}}\) is the associated mean exciton energy. Here, \({E}_{{\rm{opt}}}\) corresponds to the optical gap and is determined by the first excited singlet state of either donor \(({E}_{{\rm{opt}}}={E}_{{\rm{D}},{\rm{LE}}})\) or acceptor \(({E}_{{\rm{opt}}}={E}_{{\rm{A}},{\rm{LE}}})\) in a blend. For \(E\ll {E}_{{\rm{opt}}}\), Eq. (3) reduces to the exponential part and \({E}_{{\rm{U}}}^{{\rm{app}}}\) therefore equals \({kT}\). On the other hand, the error functions govern the spectral shape near the absorption edge as demonstrated in Supplementary Fig. 9.
To validate the model, we applied Eq. (3) to \(T\)-dependent EQE spectra, as shown in Fig. 4a for neat Y6 and IT-4F. Using neat materials avoids the influences of CT states in the low-energy tail although trap states are always present at lower energies. Doing so, we obtain \({\sigma }_{s}\) values of 47.0\(\pm\)0.7âmeV for Y6 and 35.0\(\pm\)2.6âmeV for IT-4F (see Supplementary Fig. 11). More fittings on room temperature EQE spectra of other material systems are shown in the Supplementary Fig. 12). For blend systems, the full experimental sub-gap EQE can be reconstructed assuming \(\alpha \left(E\right)={\alpha }_{{\rm{LE}}}\left(E\right)+{\alpha }_{{\rm{CT}}}\left(E\right)+{\alpha }_{{\rm{t}}}\left(E\right)\) constituting the exponentially decaying \({\alpha }_{{\rm{LE}}}\) according to Eq. (3) and the sum of two Gaussian functions (CT states and deep trap states; see Eq. 5). In Fig. 4b, we demonstrate the model for PBDB-T:PC70BM using the parameters \({E}_{{\rm{opt}}}=\) 1.89âeV and \({\sigma }_{{\rm{s}}}=\) 60âmeV, and other Gaussian fit parameters summarized in Supplementary Table 1. PBDB-T:PC70BM belongs to the group of blends for which \({E}_{{\rm{CT}}}\ll {E}_{{\rm{opt}}}\) with \({E}_{{\rm{CT}}}-{E}_{{\rm{opt}}}\approx\) 0.45âeV. Nevertheless, the simplified model reproduces the strong spectral dependence of the experimental \({E}_{{\rm{U}}}^{{\rm{app}}},\) including the parabolic behavior in the LE absorption dominated spectral range. On the other hand, by keeping all other parameters constant while increasing \({E}_{{\rm{CT}}}\) with respect to \({E}_{{\rm{opt}}}\), the emergence of the \({E}_{{\rm{U}}}^{{\rm{app}}}={kT}\) plateau at around 1.6\(\pm\)0.1âeV can be reproduced.
In the absence of CT states, the above model predicts three sub-gap regimes of the \(\alpha \left(E\right)\) based on the dominant disorder mechanism as illustrated in Fig. 4c. For photon energies close to the gap (\(E \,> \, {E}_{{\rm{opt}}}-{\sigma }_{s}^{2}/{kT}\)), \({E}_{{\rm{U}}}^{{\rm{app}}}\) is dominated by Gaussian static disorder. In this regime, the absorption resembles a Gaussian-like shape, where the steepness close to the gap is determined by \({\sigma }_{{\rm{s}}}\), causing a redshift of the effective energy gap with increasing \({\sigma }_{{\rm{s}}}\). While \({E}_{{\rm{U}}}^{{\rm{app}}}\) does not directly represent \({\sigma }_{s}\) at any energy, a positive correlation between \({E}_{{\rm{U}}}^{{\rm{app}}}\) and \({\sigma }_{{\rm{s}}}\) is observed in this regime (see Supplementary Fig. 10). A closer inspection suggests that \({E}_{{\rm{U}}}^{{\rm{app}}}\) is strongly energy dependent in this regime, scaling with \({\sigma }_{{\rm{s}}}\) as \({E}_{{\rm{U}}}^{{\rm{app}}}\left(E\right)\approx 2{\sigma }_{{\rm{s}}}^{2}/\left({E}_{X}-E\right)\), for \({\sigma }_{{\rm{s}}} < 4{kT}\), where \({E}_{X}\) is a constant independent of \(E\). This explains the previously reported correlation between morphological disorder and \({E}_{{\rm{U}}}^{{\rm{app}}}\) close to the absorption onset22,46,47,48. For \(E \,<\, {E}_{{\rm{opt}}}-{\sigma }_{s}^{2}/{kT}\), however, the exponential term in Eq. (3) eventually starts to dominate the spectral line shape of \({\alpha }_{{\rm{LE}}}\). In this regime, the sub-gap absorption of LEs is dominated by thermal broadening with \({E}_{{\rm{U}}}^{{\rm{app}}}\to {kT}\), thus becoming independent of \({\sigma }_{{\rm{s}}}\). The transition between the static disorder dominated and the thermal broadening dominated regime depends on \({\sigma }_{{\rm{s}}}\) (which is a measure of the Gaussian static disorder) and the temperature. This is simulated for the two cases \({\sigma }_{{\rm{s}}}\) \(=\) 70âmeV and \({\sigma }_{{\rm{s}}}\) \(=\) 100âmeV as shown in Supplementary Fig. 10: For \({\sigma }_{{\rm{s}}}\) \(=\) 70âmeV, \({E}_{{\rm{U}}}^{{\rm{app}}}\) reaches \({kT}\) at \(\alpha \left(E\right)\) values 3 to 4 orders of magnitudes below \({\alpha }_{{\rm{sat}}}\), while this is 6 orders below \({\alpha }_{{\rm{sat}}}\) for \({\sigma }_{{\rm{s}}}=\) 100âmeV at room temperature. Finally, at energies well below the optical gap, \(\alpha \left(E\right)\) eventually becomes dominated by deep trap state absorption, resulting in an artificial increase in \({E}_{{\rm{U}}}^{{\rm{app}}}\) and a concomitant deviation from \({E}_{{\rm{U}}}^{{\rm{app}}}={kT}\). Deep trap state absorption is typically observed 6 orders of magnitude below \({\alpha }_{{\rm{sat}}}\) limiting the spectral range dominated by thermal broadening.
Based on these considerations, we estimate that it is possible to observe \({E}_{{\rm{U}}}^{{\rm{app}}}\approx {kT}\) only when \({\sigma }_{{\rm{s}}}\) \(<\) 100âmeV (see Supplementary Fig. 10). Other conditions that must be met are: (i) the neat phase absorption of one component is spectrally separated from the other neat phase absorption, as well as from the CT states and trap states; (ii) the dynamic range of the EQE (or \(\alpha\)) measurement is sufficiently wide to measure photocurrent at wavelengths well below the absorption onset of the neat material (\(E \,<\, {E}_{{\rm{opt}}}-{\sigma }_{{\rm{s}}}^{2}/{kT}\)); and (iii) optical cavity effects are not significant. Exponential sub-gap EQE spectra previously reported in literature for BHJs often do not fulfill these requirements, explaining reported Urbach energies much larger than \({kT}\). Importantly, the \({E}_{{\rm{U}}}^{{\rm{app}}}\) spectra introduced here do not suffer from the short fitting ranges of previous of \({E}_{{\rm{U}}}\) measurements. In contrast, we have shown that an exponential distribution of tail states cannot explain the sub-gap spectral line-shape associated with singlet absorption, nor CT or trap state absorption in organic semiconductors.
Recombination losses and the radiative V OC limit set by excitons
The presence of sub-gap absorption is known to induce radiative recombination losses additional to those predicted by the Shockley and Queisser (SQ) model. While the short-circuit current density (\({J}_{{\rm{SC}}}\)) remains largely unaffected, the presence of sub-gap absorption mainly translates into enhanced losses in the open-circuit voltage (\({V}_{{\rm{OC}}}\)). In general, the \({V}_{{\rm{OC}}}\) can be expressed as \({V}_{{\rm{OC}}}={V}_{{\rm{OC}}}^{{\rm{SQ}}}-\Delta {V}_{{\rm{OC}}}^{{\rm{RAD}}}-\Delta {V}_{{\rm{OC}}}^{{\rm{NR}}}\), where \({V}_{{\rm{OC}}}^{{\rm{SQ}}}\) represents the upper thermodynamic limit of the \({V}_{{\rm{OC}}}\) based on the SQ model, assuming perfect above-gap absorption and no sub-gap absorption. Furthermore, \(\Delta {V}_{{\rm{OC}}}^{{\rm{RAD}}}={V}_{{\rm{OC}}}^{{\rm{SQ}}}-{V}_{{\rm{OC}}}^{{\rm{RAD}}}\) is the radiative loss induced by sub-gap absorption, while \(\Delta {V}_{{\rm{OC}}}^{{\rm{NR}}}={V}_{{\rm{OC}}}^{{\rm{RAD}}}-{V}_{{\rm{OC}}}\) is the non-radiative loss. Here, \({V}_{{\rm{OC}}}^{{\rm{RAD}}}\) is the radiative limit of the \({V}_{{\rm{OC}}}\) corresponding to the expected \({V}_{{\rm{OC}}}\) in the absence of non-radiative losses. In accordance with detailed balance, we expect \({V}_{{\rm{OC}}}^{{\rm{RAD}}}=\frac{{kT}}{q}{\rm{ln}}\left(\frac{{J}_{{\rm{SC}}}}{{J}_{0}^{{\rm{RAD}}}}+1\right)\), where \({J}_{{\rm{SC}}}=q{\int }_{0}^{\infty }{\rm{EQE}}\left(E\right){\Phi }_{{\rm{sun}}}\left(E\right){dE}\) and \({J}_{0}^{{\rm{RAD}}}=q{\int }_{0}^{\infty }{\rm{EQE}}\left(E\right){\Phi }_{{\rm{BB}}}\left(E\right){dE}\) is the dark saturation current density in the radiative limit; \(q\) is the elementary charge and \({\Phi }_{{\rm{sun}}}\) (\({\Phi }_{{\rm{BB}}}\)) the solar spectrum (black body spectrum)49. The SQ limit (\({V}_{{\rm{OC}}}^{{\rm{RAD}}}={V}_{{\rm{OC}}}^{{\rm{SQ}}}\)) corresponds to the case when \({\rm{EQE}}=1\) for \(E \,> \, {E}_{{\rm{opt}}}\), while \({\rm{EQE}}=0\) for \(E \,<\, {E}_{{\rm{opt}}}\). We note that, because of the different ideality factors, the contribution from deep trap states will be negligible under 1-sun conditions, as shown previously14.
As we have shown above, in low offset systems such as state-of-the-art NFA-based blends (or neat material systems), the CT absorption is mostly overshadowed by the stronger NFA absorption (or is absent). In such systems the sub-gap absorption tail is dominated by excitons (\(\alpha ={\alpha }_{{\rm{LE}}}\)). The expected radiative limit \({V}_{{\rm{OC}}}^{{\rm{RAD}}}\) set by excitons can be calculated by assuming \({\rm{EQE}}={{\rm{EQE}}}_{{\rm{max }}}\times \left({\alpha }_{{\rm{LE}}}/{\alpha }_{{\rm{sat}}}\right)\), with \({\alpha }_{{\rm{LE}}}\) given by Eq. (3). The concomitant radiative loss \(\Delta {V}_{{\rm{OC}}}^{{\rm{RAD}}}\), arising from the imperfect spectral line shape of \({\alpha }_{{\rm{LE}}}\) (leading to deviations from the ideal box-like EQE profile), is shown in Fig. 5a at different \({E}_{{\rm{opt}}}\) and \({\sigma }_{{\rm{s}}}\). As shown, the voltage loss \(\Delta {V}_{{\rm{OC}}}^{{\rm{RAD}}}\) increases with increasing \({\sigma }_{{\rm{s}}}\). For example, when \({\sigma }_{{\rm{s}}}=100\) meV, \(\Delta {V}_{{\rm{OC}}}^{{\rm{RAD}}}\) is larger than 0.2âV. Furthermore, at large \({\sigma }_{{\rm{s}}}\), the voltage loss becomes more prominent at larger \({E}_{{\rm{opt}}}\). We note that these losses are caused by a drastic increase in \({J}_{0}^{{\rm{RAD}}}\), while \({J}_{{\rm{SC}}}\) is essentially unchanged. As \({\sigma }_{{\rm{s}}}\) is reduced, in turn, \(\Delta {V}_{{\rm{OC}}}^{{\rm{RAD}}}\) is correspondingly decreased. For small static disorder (i.e., small \({\sigma }_{{\rm{s}}}\)), an analytical approximation relating \(\Delta {V}_{{\rm{OC}}}^{{\rm{RAD}}}\) and \({\sigma }_{{\rm{s}}}\) can be obtained as
A good agreement between Eq. (4), as indicated by the dashed lines in Fig. 5a, and the full simulation is indeed obtained for small \({\sigma }_{{\rm{s}}}\). Importantly, even in the limit of vanishing static disorder, \({\sigma }_{{\rm{s}}}\to 0\), a voltage loss \(\Delta {V}_{{\rm{OC}}}^{{\rm{RAD}}}\)of roughly 70â80âmeV is expected due to non-vanishing sub-gap absorption induced by thermal broadening (manifested by \(\alpha \propto {\rm{exp }}\left(E/{kT}\right)\) in the gap). Based on the \({E}_{{\rm{opt}}}\) and \({\sigma }_{{\rm{s}}}\) extracted (see Supplementary Fig. 12) for the low-offset D:A systems PM6:Y6 (\({E}_{{\rm{opt}}}=1.44\) eV, \({\sigma }_{{\rm{s}}}=42\) meV), PM6:ITIC (\({E}_{{\rm{opt}}}=1.69\) eV, \({\sigma }_{{\rm{s}}}=37\) meV) and PBDB-T:EH-IDTBR (\({E}_{{\rm{opt}}}=1.77\) eV, \({\sigma }_{{\rm{s}}}=46\) meV), the radiative voltage loss induced by sub-gap absorption can be estimated. Subsequently, we find \(\Delta {V}_{{\rm{OC}}}^{{\rm{RAD}}}\) of 108âmV for PM6:Y6, 107âmV for PM6:ITIC, and 120âmV for PBDB-T:EH-IDTBR. We note that the corresponding would-be \(\Delta {V}_{{\rm{OC}}}^{{\rm{RAD}}}\) (if CT states were absent, see Fig. 4b) for PBDB-T:PC70BM is 150âmV.
Finally, Fig. 5b shows the corresponding effect of \({\sigma }_{{\rm{s}}}\) on the radiative limit of the PCE, assuming \({{\rm{EQE}}}_{{\rm{max }}}=1\) and ideal charge collection. Compared to the SQ limit of the PCE (indicated by black solid line), the presence of the sub-gap absorption results in a PCE peak loss of around 1.5 % for \({\sigma }_{{\rm{s}}}=0\). For \({\sigma }_{{\rm{s}}}\ne 0\), the radiative PCE limit is further lowered, with the PCE peak deceasing with increasing \({\sigma }_{{\rm{s}}}\). Hence, to minimize \(\Delta {V}_{{\rm{OC}}}^{{\rm{RAD}}}\), and thus maximize PCE, it is important to minimize \({\sigma }_{{\rm{s}}}\). However, an unavoidable radiative loss, relative to the SQ limit, will still be present due to thermal broadening. Based on our findings, this loss is inherent to all organic solar cells, and needs to be taken into account in low offset systems where excitons dominate the sub-gap absorption.
Discussion
In summary, we have shown that the exciton sub-gap absorption in organic semiconductors is dominated by Gaussian static disorder near the onset, while thermal broadening dominates at lower photon energies. Furthermore, we find that for a large number of organic semiconductors, the Urbach energy in the sub-gap absorption regime dominated by thermal broadening equals the thermal energy \({kT}\) within the variations caused by optical interference. While exponential absorption tails have been previously observed mainly in non-fullerene systems, we show that this property is universal for organic semiconductors and consistent with a Gaussian density of excitonic states undergoing Boltzmann-like thermally activated optical transitions. The static disorder is shown to arise from the width of the Gaussian DOS and to broaden the absorption onset at energies close to the optical gap. Using our model for the sub-gap excitonic absorption tail, it is possible to discriminate spectral regimes dominated by static disorder (where the Urbach energy is strongly dependent on the energy) and thermal broadening (\({E}_{{\rm{U}}}={kT}\)), as well as to reproduce the temperature dependence of absorption coefficient due to excitons. This modified view of the sub-gap absorption coefficient in disordered organic semiconductors clarifies a longstanding debate concerning the shape of the DOS and the relevance of an Urbach description in these important and intriguing materials. Finally, we demonstrate the implications of exciton static disorder and thermal broadening on the radiative open-circuit voltage losses in organic solar cells based on low offset D:A blends.
Methods
Materials
Poly(3,4-ethylenedioxythiophene) polystyrene sulfonate (PEDOT:PSS) was purchased from Heraeus. Zinc acetate dihydrate and PCDTBT (Poly[N-9â²-heptadecanyl-2,7-carbazole-alt-5,5-(4â²,7â²-di-2-thienyl-2â²,1â²,3â²-benzothiadiazole)]) were purchased from Sigma Aldrich. PC70BM ([6,6]-Phenyl-C71-butyric acid methyl ester) and EH-IDTBR were purchased from Solarmer (Beijing). BQR (benzodithiophene-quaterthiophene-rhodanine) was provided by Prof. David. J Jones (University of Melbourne). m-MTDATA (4,4â²,4â-Tris[(3-methylphenyl) phenylamino]triphenylamine) was purchased from Ossila. IT-4F (3,9-bis(2-methylene-((3-(1,1-dicyanomethylene)-6,7-difluoro)-indanone))-5,5,11,11-tetrakis(4-hexylphenyl)-dithieno[2,3-d:2â²,3â²-dâ]-s-indaceno[1,2-b:5,6-bâ]dithiophene), PM6 (Poly[(2,6-(4,8-bis(5-(2-ethylhexyl-3-fluoro)thiophen-2-yl)-benzo[1,2-b:4,5-bâ]dithiophene))-alt-(5,5-(1â²,3â²-di-2-thienyl-5â²,7â²-bis(2-ethylhexyl)benzo[1â²,2â²-c:4â²,5â²-câ]dithiophene-4,8-dione)]), Y6 (2,2â²-((2Z,2â²Z)-((12,13-bis(2-ethylhexyl)-3,9-diundecyl-12,13-dihydro-[1,2,5]thiadiazolo[3,4-e]thieno[2â,3â:4â²,5â²]thieno[2â²,3â²:4,5]pyrrolo[3,2-g]thieno[2â²,3â²:4,5]thieno[3,2-b]indole-2,10-diyl)bis(methanylylidene))bis(5,6-difluoro-3-oxo-2,3-dihydro-1H-indene-2,1-diylidene))dimalononitrile), ITIC (3,9-bis(2-methylene-(3-(1,1-dicyanomethylene)-indanone))-5,5,11,11-tetrakis(4-hexylphenyl)-dithieno[2,3-d:2â²,3â²-dâ]-s-indaceno[1,2-b:5,6-bâ]dithiophene), PBDB-T (Poly[(2,6-(4,8-bis(5-(2-ethylhexyl)thiophen-2-yl)-benzo[1,2-b:4,5-bâ]dithiophene))-alt-(5,5-(1â²,3â²-di-2-thienyl-5â²,7â²-bis(2-ethylhexyl)benzo[1â²,2â²-c:4â²,5â²-câ]dithiophene-4,8-dione)]) and PTB7-Th (Poly[4,8-bis(5-(2-ethylhexyl)thiophen-2-yl)benzo[1,2-b;4,5-bâ]dithiophene-2,6-diyl-alt-(4-(2-ethylhexyl)-3-fluorothieno[3,4-b]thiophene-)-2-carboxylate-2â6-diyl)]) were purchased from Zhi-yan (Nanjing) Inc.
Device fabrication
Solar cells were fabricated with either a conventional architecture Indium tin oxide (ITO)/PEDOT:PSS/active layer/Ca/Al or inverted architecture ITO/ZnO/active layer/MoO3/Ag. Commercial ITO coated glass substrates from Ossila were cleaned in an aqueous solution of Alconox at 60â°C, followed by an ultrasonic bath in deionize water, acetone and isopropanol. The cleaned substrates were dried with nitrogen, followed by a UV/O3 treatment (Ossila, L2002A2-UK). For the conventional device architecture, 30ânm of PEDOT:PSS was spin-coated at 6000ârpm for 30âs onto precleaned ITO substrates and annealed at 155â°C for 15âmin. As top electrode, 20ânm of calcium (Ca) and 100ânm of Aluminum (Al) were vacuum deposited at 10â6 Tor defining an active area of 4 mm2. For the inverted device architecture with 30ânm of ZnO, a solution of 200âmg of zinc acetate dihydrate in 2-methoxyethanol (2âml) and ethanolamine (56âµl) was prepared and stirred overnight under ambient conditions. The ZnO layer formed upon spin-coating the solution at 4000ârpm followed by thermal annealing at 200â°C for 60âmin. As top electrode, 7ânm of MoO3 and 100ânm of Ag were vacuum deposited at 10â6 Tor defining an active area of 4 mm2.
Devices with conventional structure
BQR:PC70BM devices were prepared using with the conventional architecture ITO/PEDOT:PSS/BQR:PC70BM/Ca/Al. BQR and PC70BM were dissolved in toluene (24âmg/ml with the donor:acceptor ratio of 1:1) and stirred at 60â°C for 3âh. Next, the BQR:PC70BM solution was spin-coated at 1000ârpm on the PEDOT:PSS layer to form a 100ânm thick film.
Devices with inverted structure
m-MTDATA:PC70BM devices: Equimolar solutions of PC70BM and m-MTDATA in dichloromethane (DCM) with a concentration of 19.4âmmol/l were prepared. To obtain a series of solutions with different molar ratios of m-MTDATA:PC70BM (5âmol%, 1âmol%, 0.1âmol% and 0âmol% of m-MTDATA), 50 µl, 10 µl, 1 µl and 0 µl of m-MTDATA in DCM were added to 1âml of PC70BM in DCM. The solutions were spin-coated at a spin rate of 800ârpm to get an active layer thickness of around 90ânm. PM6:Y6 devices: PM6:Y6 was dissolved in chloroform (CF) solution (14âmgâmlâ1 with 0.5 vol.% 1-Chloronaphthalene [CN]) with a donor:acceptor ratio of 1:1.2, and spin-coated (3000ârpm) on ZnO to form a 100ânm thick film. The as cast active layers were thermally annealed at 110â°C for 10âmin. PM6:ITIC devices: PM6:ITIC was dissolved in chlorobenzene (CB) solution (18âmgâmlâ1 with 0.5 vol.% DIO) with a donor:acceptor ratio of 1:1, and spin-coated (1000ârpm) on ZnO to form a 100ânm thick film. The active layers were further treated with thermal annealing at 100â°C for 10âmin. PBDB-T:EH-IDTBR devices: PBDB-T:EH-IDTBR was dissolved in CB solution (14âmgâmlâ1) with a donor:acceptor ratio of 1:1, and spin-coated (800ârpm) on ZnO to form a 100ânm thick film. PBDB-T:ITIC devices: PBDB-T:ITIC was dissolved in CB solution (14âmgâmlâ1 with 0.5 vol.% DIO) with a donor:acceptor ratio of 1:1, and spin-coated (800ârpm) on ZnO to form a 100ânm thick film. The active layers were further treated with thermal annealing at 100â°C for 10âmin. PTB7-Th:ITIC devices: PTB7-Th:ITIC was dissolved in CB solution (14âmgâmlâ1 with 1 vol.% 1,8-diiodooctane [DIO]) with a donor:acceptor ratio of 1:1.4, and spin-coated (1000ârpm) on ZnO to form a 100ânm thick film. PBDB-T:PC70BM devices: PBDB-T:PC70BM was dissolved in CB solution (14âmgâmlâ1 with 3 vol.% DIO) with a donor:acceptor ratio of 1:1.4, and spin-coated (1000ârpm) on ZnO to form a 100ânm thick film. Then the as-cast films were rinsed with 80âμL of methanol at 4000ârpm for 20âs to remove the residual DIO. PBDB-T:IT-4F devices: PBDB-T:IT-4F was dissolved in CB solution (14âmgâmlâ1 with 0.5 vol.% DIO) with a donor:acceptor ratio of 1:1, and spin-coated (800ârpm) on ZnO to form a 100ânm thick film. The active layers were further treated with thermal annealing at 100â°C for 10âmin. Neat ITIC devices: ITIC was dissolved in CF solution (10âmgâmlâ1) and spin-coated on ZnO (2000 rpm) to form a 70ânm thick film. Neat IT-4F devices: IT-4F was dissolved in CF solution (10âmgâmlâ1) and spin-coated on ZnO (2000 rpm) to form a 70ânm thick film. Neat PBDB-T devices: PBDB-T was dissolved in CF solution (10âmgâmlâ1) and spin-coated on ZnO (3000ârpm) to form a 70ânm thick film. Neat Y6 devices: Y6 was dissolved in CF solution (16âmgâmlâ1) and spin-coated on ZnO (3000ârpm) to form a 70ânm thick film.
EQE measurements
A homebuilt setup including a Perkin Elmer UV/VIS/NIR spectrometer (LAMBDA 950) as a source for monochromatic light was used. The light was chopped at 273âHz and directed onto the device under test (DUT). The resulting photocurrent was amplified by a low noise current amplifier (FEMTO DLPCA-200) and measured with the Stanford SR860 lock-in amplifier. To decrease the noise floor of the setup, the DUT was mounted in an electrically shielded and temperature controlled Linkam sample stage. An integration time up to 1000âs was used for detecting wavelengths above 1500ânm. NISTâcalibrated silicon and Ge photodiodes from Newport were used as a calibration reference. The temperature inside the Linkam sample stage was set to â120 to 60â°C by the Linkam T96 temperature controller in combination with an LNP96 liquid nitrogen pump. The commercial amorphous silicon thin film solar cell, used for temperature dependent EQE measurements, was manufactured by TRONY with the part number sc80125s-8.
Gaussian fits
The EQE associated with CT states and mid-gap states were routinely fitted in accordance with the standard Marcus charge-transfer formalism. The associated EQEs are given by \({{\rm{EQE}}}_{{\rm{CT}}}\left(E\right)=g(E,{E}_{{\rm{CT}}},{\lambda }_{{\rm{CT}}},{f}_{{\rm{CT}}})\) and \({{\rm{EQE}}}_{{\rm{t}}}\left(E\right)=g(E,{E}_{{\rm{t}}},{\lambda }_{{\rm{t}}},{f}_{{\rm{t}}})\) for CT and mid-gap state absorption, respectively, where
Here, \({f}_{j}\), \({E}_{j}\) and \({\lambda }_{j}\) are fitting parameters. All fit parameters obtained in this work are summarized in Supplementary Table 1. The EQEs are assumed to be related to their respective absorption coefficients as\(\,{\alpha }_{{\rm{CT}}}\left(E\right)\propto {{\rm{EQE}}}_{{\rm{CT}}}\left(E\right)\) and \({\alpha }_{{\rm{t}}}\left(E\right)\propto {{\rm{EQE}}}_{{\rm{t}}}\left(E\right)\) (neglecting interference effects).
Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
References
Zhao, W. et al. Molecular optimization enables over 13% efficiency in organic solar cells. J. Am. Chem. Soc. 139, 7148â7151 (2017).
Jin, K., Xiao, Z. & Ding, L. 18.69% PCE from organic solar cells. J. Semicond. 42, 060502 (2021).
Urbach, F. The long-wavelength edge of photographic sensitivity and of the electronic Absorption of Solids. Phys. Rev. 92, 1324 (1953).
Cody, G. D., Tiedje, T., Abeles, B., Brooks, B. & Goldstein, Y. Disorder and the optical-absorption edge of hydrogenated amorphous silicon. Phys. Rev. Lett. 47, 1480â1483 (1981).
Bonalde, I. et al. Urbach tail, disorder, and localized modes in ternary semiconductors. Phys. Rev. B 69, 195201 (2004).
Bansal, B., Dixit, V. K., Venkataraman, V. & Bhat, H. L. Alloying induced degradation of the absorption edge of InAsxSb1-x. Appl. Phys. Lett. 90, 101905 (2007).
Pankove, J. I. Absorption edge of impure gallium arsenide. Phys. Rev. 140 (1965).
Kurik, M. V. Urbach rule. Phys. Status Solidi 8, 9â45 (1971).
Grein, C. H. & John, S. Temperature dependence of the Urbach optical absorption edge: a theory of multiple phonon absorption and emission sidebands. Phys. Rev. B 39, 1140â1151 (1989).
Corkish, R. & Green, M. A. Erratum: âBand edge optical absorption in intrinsic silicon: assessment of the indirect transition and disorder modelsâ [J. Appl. Phys. 73, 3988 (1993)]. J. Appl. Phys. 74, 6462â6462 (1993).
Dow, J. D. & Redfield, D. Toward a unified theory of Urbachâs rule and exponential absorption edges. Phys. Rev. B 5, 594â610 (1972).
Shklovskii, B. I & Efros, A. L. Electronic properties of doped semiconductors. (Springer, Berlin Heidelberg, 1984).
Rahimi, K. et al. Light absorption of poly(3-hexylthiophene) single crystals. RSC Adv. 4, 11121â11123 (2014).
Zarrabi, N. et al. Charge-generating mid-gap trap states define the thermodynamic limit of organic photovoltaic devices. Nat. Commun. 11, 5567 (2020).
Vandewal, K., Tvingstedt, K., Gadisa, A., Inganäs, O. & Manca, J. V. Relating the open-circuit voltage to interface molecular properties of donor:acceptor bulk heterojunction solar cells. Phys. Rev. B. 81, 1â8 (2010).
Burke, T. M., Sweetnam, S., Vandewal, K. & Mcgehee, M. D. Beyond Langevin recombination: how equilibrium between free carriers and charge transfer states determines the open-circuit voltage of organic solar cells. Adv. Energy Mater. 1500123, 1â12 (2015).
Kahle, F. J., Rudnick, A., Bässler, H. & Köhler, A. How to interpret absorption and fluorescence spectra of charge transfer states in an organic solar cell. Mater. Horiz. 5, 837â848 (2018).
Tvingstedt, K., Benduhn, J. & Vandewal, K. Temperature dependence of the spectral line-width of charge-transfer state emission in organic solar cells; static vs. dynamic disorder. Mater. Horiz. 7, 1888â1900 (2020).
Felekidis, N., Melianas, A. & Kemerink, M. The role of delocalization and excess energy in the quantum efficiency of organic solar cells and the validity of optical reciprocity relations. J. Phys. Chem. Lett. 11, 3563â3570 (2020).
Qian, D. et al. Design rules for minimizing voltage losses in high-efficiency organic solar cells. Nat. Mater. 17, 703â709 (2018).
Ran, N. A. et al. Harvesting the full potential of photons with organic solar cells. Adv. Mater. 28, 1482â1488 (2016).
Liu, S. et al. High-efficiency organic solar cells with low non-radiative recombination loss and low energetic disorder. Nat. Photonics 14, 300â305 (2020).
Jain, N. et al. Interfacial disorder in efficient polymer solar cells: the impact of donor molecular structure and solvent additives. J. Mater. Chem. A 5, 24749â24757 (2017).
Bässler, H. Charge transport in disordered organic photoconductors a Monte Carlo simulation study. Phys. Status Solidi 175, 15â56 (1993).
Beenken, W. J. D. et al. Sub-bandgap absorption in organic solar cells: experiment and theory. Phys. Chem. Chem. Phys. 15, 16494â16502 (2013).
Blakesley, J. C. & Neher, D. Relationship between energetic disorder and open-circuit voltage in bulk heterojunction organic solar cells. Phys. Rev. B 84, 075210 (2011).
Kaiser, C., Zeiske, S., Meredith, P. & Armin, A. Determining ultralow absorption coefficients of organic semiconductors from the subâbandgap photovoltaic external quantum efficiency. Adv. Opt. Mater. 8, 1901542 (2019).
Armin, A. et al. Limitations of charge transfer state parameterization using photovoltaic external quantum efficiency. Adv. Energy Mater. 10, 2001828 (2020).
Armin, A. et al. Spectral dependence of the internal quantum efficiency of organic solar cells: Effect of charge generation pathways. J. Am. Chem. Soc. 136, 11465â11472 (2014).
Vandewal, K. et al. Efficient charge generation by relaxed charge-transfer states at organic interfaces. Nat. Mater. 13, 63â68 (2014).
Kurpiers, J. et al. Probing the pathways of free charge generation in organic bulk heterojunction solar cells. Nat. Commun. 9, 1â11 (2018).
Lee, J. et al. Charge transfer state versus hot exciton dissociation in polymer-fullerene blended solar cells. J. Am. Chem. Soc. 132, 11878â11880 (2010).
Buchaca-Domingo, E. et al. Direct correlation of charge transfer absorption with molecular donor:acceptor interfacial area via photothermal deflection spectroscopy. J. Am. Chem. Soc. 137, 5256â5259 (2015).
Sandberg, O. J., Zeiske, S., Zarrabi, N., Meredith, P. & Armin, A. Charge carrier transport and generation via trap-mediated optical release in organic semiconductor devices. Phys. Rev. Lett. 124, 128001 (2020).
Gielen, S. et al. Intrinsic detectivity limits of organic nearâinfrared photodetectors. Adv. Mater. 32, 2003818 (2020).
Nikolis, V. C. et al. Strong light-matter coupling for reduced photon energy losses in organic photovoltaics. Nat. Commun. 10, 1â8 (2019).
Göhler, C. et al. Temperature-dependent charge-transfer-state absorption and emission reveal the dominant role of dynamic disorder in organic solar cells. Phys. Rev. Appl. 15, 064009 (2021).
Searle, T. M. & Jackson, W. A. Static versus electron-phonon disorder in amorphous Si: H and its alloys. Philos. Mag. B 60, 237â255 (1989).
Cody, G. D. Urbach edge of crystalline and amorphous silicon: a personal review. J. Non Cryst. Solids 141, 3â15 (1992).
Leblanc, F., Maeda, Y., Ando, M., Wakagi, M. & Minemura, T. Accurate determination of the Urbach energy of a-Si:H thin films by correction for the interference effect. Jpn. J. Appl. Phys. 33, L1755âL1758 (1994).
Tang, J. Electron-transfer reactions involving non-linear spin-boson interactions. Chem. Phys. 188, 143â160 (1994).
Casado-Pascual, J., Morillo, M., Goychuk, I. & Hänggi, P. The role of different reorganization energies within the Zusman theory of electron transfer. J. Chem. Phys. 118, 291â303 (2003).
Keil, T. H. Theory of the Urbach rule. Phys. Rev. 144, 582â587 (1966).
Hood, S., Zarrabi, N., Meredith, P., Kassal, I. & Armin, A. Measuring energetic disorder in organic semiconductors using the photogenerated charge-separation efficiency. J. Phys. Chem. Lett. 10, 3863â3870 (2019).
Miller, A. & Abrahams, E. Impurity conduction at low concentrations. Phys. Rev. 120, 745â755 (1960).
Cui, Y. et al. Single-junction organic photovoltaic cells with approaching 18% efficiency. Adv. Mater. 32, 1â7 (2020).
Song, B., Rolin, C., Zimmerman, J. D. & Forrest, S. R. Effect of mixed layer crystallinity on the performance of mixed heterojunction organic photovoltaic cells. Adv. Mater. 26, 2914â2918 (2014).
Venkateshvaran, D. et al. Approaching disorder-free transport in high-mobility conjugated polymers. Nature 515, 384â388 (2014).
Shockley, W. & Queisser, H. J. Detailed balance limit of efficiency of p-n junction solar cells. J. Appl. Phys. 32, 510â519 (1961).
Acknowledgements
This work was funded through the Welsh Governmentâs Sêr Cymru II Program âSustainable Advanced Materialsâ (Welsh European Funding Office â European Regional Development Fund). C.K. is recipient of a UKRI EPSRC Doctoral Training Program studentship. N.Z. is funded by a studentship through the Sêr Cymru II Program. P.M. is a Sêr Cymru II Research Chair and A.A. is a Rising Star Fellow also funded through the Welsh Governmentâs Sêr Cymru II âSustainable Advanced Materialsâ Program (European Regional Development Fund, Welsh European Funding Office and Swansea University Strategic Initiative). This work was also funded by UKRI through the EPSRC Program Grant EP/T028511/1 Application Targeted Integrated Photovoltaics.
Author information
Authors and Affiliations
Contributions
A.A. and P.M. provided the overall leadership of the project. O.J.S. and A.A. conceptualized the idea. C.K. and A.A. designed the experiments. C.K. performed most measurements, analyzed the data, and performed optical modeling. OJS developed the theoretical model. C.K., O.J.S. and A.A. interpreted the data. N.Z. assisted with \({{\rm{EQE}}}_{{\rm{PV}}}\) measurements. W.L. fabricated the devices. All co-authors contributed in the development of the manuscript which was initially drafted by C.K.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks Martijn Kemerink and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kaiser, C., Sandberg, O.J., Zarrabi, N. et al. A universal Urbach rule for disordered organic semiconductors. Nat Commun 12, 3988 (2021). https://doi.org/10.1038/s41467-021-24202-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-021-24202-9
This article is cited by
-
Physical insights into non-fullerene organic photovoltaics
Nature Reviews Physics (2024)
-
Highly sensitive water pollution monitoring using colloid-processed organic photodetectors
Nature Water (2024)
-
Direct bandgap quantum wells in hexagonal Silicon Germanium
Nature Communications (2024)
-
A case study of comparing two dimerized acceptor molecules built by different branch-connected and terminal-connected approaches
Science China Chemistry (2024)
-
Enhancing Photocatalyst Performance of Magnetic Surfaces Covered by Carbon Clouds for Textile Dye Degradation
Arabian Journal for Science and Engineering (2024)