A critical bottleneck for improving the performance of organic solar cells (OSC) is minimising non-radiative losses in the interfacial charge-transfer (CT) state via the formation of hybrid energetic states. This requires small energetic offsets often detrimental for high external quantum efficiency (EQE). Here, we obtain OSC with both non-radiative voltage losses (0.24 V) and photocurrent losses (EQE > 80%) simultaneously minimised. The interfacial CT states separate into free carriers with ≈40-ps time constant. We combine device and spectroscopic data to model the thermodynamics of charge separation and extraction, revealing that the relatively high performance of the devices arises from an optimal adjustment of the CT state energy, which determines how the available overall driving force is efficiently used to maximize both exciton splitting and charge separation. The model proposed is universal for donor:acceptor (D:A) with low driving forces and predicts which D:A will benefit from a morphology optimization for highly efficient OSC.
Understanding the mechanism of non-radiative losses in organic photovoltaics is crucial to improve the performance further. Here, the authors use combined device and spectroscopic data to reveal universal model to maximise exciton splitting and charge separation by adjusting the energy of charge transfer state.
According to the Shockley and Queisser description, an ideal solar cell exhibits only unavoidable radiative losses1. Operated as a light-emitting diode, such a device should deliver electroluminescence (EL) external quantum efficiency (EQEEL) of unity2. For inorganic (GaAs) and hybrid (perovskites) solar cells, EQEEL close to 100% is indeed observed, while in organic photovoltaic (OPV) devices, EQEEL is orders of magnitude lower, indicating the dominance of non-radiative pathways for charge recombination3–6. These losses are mainly caused by recombination of the charge transfer (CT) states generated at the interface between the electron donor (D) and acceptor (A) materials7–11. Because the energy of the CT state (ECT) is lower than the energy (ELE) of the localised excitations (LE), the rate of non-radiative recombination is expected to increase following the energy gap law12–14. The recent introduction of non-fullerene acceptors (NFAs), with tunable energy levels and optical absorption, greatly improved the power conversion efficiencies (PCE) in OPV, which have now reached 17% in single-junction devices15–18. However, the non-radiative open circuit voltage losses () are still substantial and reduction of
Here, we demonstrate a pathway to minimise non-radiative voltage losses in an OPV device without penalizing EQE. We introduce a bulk heterojunction (BHJ) solar cell based on polymer:NFA26 blends that combine relatively high power conversion efficiency (PCE > 12%) with high EQEEL (10−4) and low non-radiative loss
The medium bandgap D–A copolymer WF3 was selected as electron donor material, while the electron acceptor was either the fullerene derivative PC70BM or the small molecule NFA O-IDTBR (Fig. 1a). The optical bandgap of neat WF3 is 1.85 eV and this value is not affected when blended with PC70BM. The optical bandgap of the WF3:O-IDTBR blend is 1.6 eV (Fig. 1b), which lies between the optical bandgaps of neat O-IDTBR annealed and as-deposited films (1.58 and 1.65 eV, respectively, see Supplementary Fig. 2). This shows that the absorption onset of WF3:O-IDTBR blends can be assigned to small aggregates of O-IDTBR in contact with WF3; the absorption spectra show no evidence of an interfacial CT state (Supplementary Fig. 3).


Steady-state characterisation of organic solar cells.
a Chemical structures of the donor and acceptor materials used in this work. b Absorption spectra of WF3:O-IDTBR and WF3:PC70BM blends. c Current density–voltage characteristics under 1 sun illumination for WF3:O-IDTBR and WF3:PC70BM devices. d External quantum efficiency measurements of the same devices as shown in (c). e Sensitive EQE measurements (EQEPV) for the same devices as shown in (c). f Electroluminescence spectra of the pristine materials and the WF3:O-IDTBR blend. The inset depicts a photograph of the WF3:O-IDTBR electroluminescence (the original photograph can be found in the SI). g EQEEL of WF3:O-IDTBR and WF3:PC70BM devices as a function of the driving current.
Figure 1c shows the current-voltage characteristics measured at 1 sun condition for the OSCs (in an inverted configuration, for details see Supplementary Information 9). WF3:PC70BM and WF3:O-IDTBR based devices delivered a PCE of 9% and 12%, respectively, the latter showing a high open-circuit voltage VOC = 1.06 V and high EQE over 80% in the 1.6–1.8 eV region (Fig. 1d). Compared to the optical bandgap Eg, calculated by the intersection between the absorption and PL spectra, the total VOC loss in the WF3:O-IDTBR device is only
Non-radiative voltage losses are strongly reduced when ELE and ECT are very similar22,29. However, different from Ref. 22. where the blends with the lowest
To get insight into the dynamics of the CT and separation processes, we use femtosecond TA spectroscopy. The samples are excited by a 10-fs pulse, mainly in resonance with WF3 excitons (shaded grey in Fig. 2C), and probed with a broadband continuum covering the 1.3–2.2 eV photon energy range. We choose broadband excitation because the resulting 10-fs pulses yield a much better time resolution compared to standard TA measurements with a time resolution of the order of 100–200 fs, allowing us to observe photophysical processes from their very onset. Figure 2C shows TA spectra at 20 ps pump-probe delay, following excitation with a 10 fs pulse (shaded grey), of films of neat WF3 (blue dots) and annealed O-IDTBR (green dots), along with the fits obtained from the global analysis (solid lines, see Supplementary Section 2 for modelling details) and the normalised ground state absorption, shown in dotted lines. The TA of WF3 includes a photobleaching (PB) signal peaking at 1.87 eV, corresponding to the bandgap absorption and its vibronic progression, a broad photoinduced absorption (PA) band in the near-infrared (NIR), and no stimulated emission (SE). These signatures are due to the formation of a dark interchain exciton (also known as a polaron pair) within 1 ps (see Supplementary Section 2), which being an excited state retains PB but suppresses SE, so that the weak steady-state emission from WF3 films can be explained by the Boltzmann equilibrium between bright singlet excitons and the less energetic dark interchain excitons (see Supplementary Fig. 5). The TA spectrum of annealed neat O-IDTBR films, on the other hand, contains a PA band in the NIR, a PB corresponding to the vibronic progression from the absorption spectrum (at 1.68 and 1.9 eV) along with its mirror image SE peak at 1.4 eV, thus corresponding to singlet excitons. The (00) vibronic transitions of the SE and PB bands are nearly completely superposed. Spectral decomposition (Supplementary Table 1) yields a Stokes shift of only about 80 meV, showing that there is very little geometric relaxation along with low energy torsional modes.


Ultrafast spectroscopy of the OPV blend.
A Schematic representation of states contributing to the TA spectra. Open shapes represent organic semiconductors in the ground state. Coloured shapes are organic semiconductors either in the electronically excited state (green: emissive O-IDTBR exciton; blue: WF3 polaron pair), or in a state created by a charge transfer event (gold: interfacial charge transfer state or charge-separated state), or in the electronic ground state but exposed to field lines of a charge pair, exhibiting a Stark shift of the fundamental electronic transition (light blue for WF3, light green for O-IDTBR). B Contributions to the TA spectra of these states. Identical bands occurring in several states are connected by dashed vertical lines. C TA spectra (dots) and fit functions (solid lines): the spectra were recorded 20 ps after excitation with 10 fs broadband pulses at parallel polarization with the probe pulses for films of pristine WF3 (blue) and annealed O-IDTBR (green). Dotted lines correspond to the normalised linear absorption and the shaded grey area shows the excitation spectrum. D TA of the WF3:O-IDTBR blend at 50 fs delays following photoexcitation (black dots), along with the fit (black line) and its components, following the same colour code as in panel (B). E Same as (D) for 200 ps after photoexcitation. F Dynamics of the concentration of each species included in the spectral modelling of the WF3:O-IDTBR blends (see Supplementary Information). Black dots correspond to the ratio EA/(CT + CS) that enables tracking the first jump of a charge away from the interface. The black solid line is an exponential fit with a 40.5 ps rise time.
To discuss the TA of the blend, we refer to Fig. 2A, where we schematically represent an area around the D/A interface and the possible photoinduced states, and Fig. 2B, showing the corresponding simplified cartoons of the TA signatures. Besides the LEs from WF3 and O-IDTBR (shown in blue and green, respectively), charged states are possible both at the interface (CT, shown in gold) or when at least one charge has diffused away from the interface (charge-separated states, CS, also shown in gold). Both CT and CS states are represented by the same colour because they present very similar intrinsic transient signatures, including PB of both charged organic semiconductors as well as a PA band in the NIR, whose spectral sensitivity to the environment (bulk or interface) we cannot probe due to spectral limitations. However, they can be distinguished by consideration of the Stark effect that they introduce upon neighbouring molecules by the electric field between positive and negative charges: as shown in Fig. 2A, the strongest Stark effect is experienced if there is a molecule between the positive and negative charge, which is the case only for a CS state but not for a CT state (for a simple electrostatic calculation, see Supplementary Section 5). These molecules will then produce a TA signal shaped as the derivative of their ground state absorption, which is often referred to as “electroabsorption” (EA)30.
Figure 2D shows the TA spectrum of the WF3:O-IDTBR blend 50 fs after excitation with 10 fs broadband pulses predominantly in resonance with the WF3 exciton (black dots). Pump and probe beams have parallel polarization; we also performed the study at magic angle polarization, qualitatively showing similar results. For a quantitative evaluation, we preferred parallel polarization due to the higher signal quality.
Derivative-shaped EA features of both WF3 and O-IDTBR are already clearly present, resulting in a spectrum that is very different from the one that would result from the superposition of the TA spectra from the pure materials in Fig. 2C. Hence, CT after WF3 excitation is extremely rapid and not diffusion-limited, outcompeting the formation of polaron pairs seen in the pristine WF3 films. This can be associated with high excess energy available for CT generation, allowing the formation of non-equilibrium, band-like states31. The observation is in agreement with a recent demonstration of ultrafast electron transfer across the donor–acceptor interface after exciting donors in NFA-blends32. Since the formation of CT states from donor excitons outperforms all possible deactivation channels, it follows that the photoexcitation dynamics of our systems only depend on the HOMO–HOMO offsets; the LUMO–LUMO offsets do not contribute quantitatively as long as it is above a certain threshold value allowing the observed ultrafast transfer33.
The spectral decomposition, as shown in Fig. 2D, yields a small initial population of acceptor excitons (dark green spectrum). The relative strength of this contribution approximately agrees with the expected amount of resonantly created acceptor excitons (about 30%), obtained by multiplication of the pump pulse spectral density with the pure donor and acceptor absorption spectra (Supplementary Information, Part 10). Therefore, we can associate the initial population of O-IDTBR excitons with resonant creation, showing that ultrafast energy transfer from donor to acceptor excitons does not play a significant role in our samples.
It has been shown that the ratio EA/(CT + CS) between the strength of the EA signal and the time-dependent densities of CT and CS states can be used to trace the separation of CT states into CS states34. To find EA/(CT + CS), we decompose the TA spectra of the WF3:O-IDTBR blend in Fig. 2D, E into contributions from O-IDTBR excitons (green), CT/CS states (gold), and EA contributions from the WF3 (light blue) and O-IDTBR (light green), using all spectral shapes of these states available from ground-state absorption, photoluminescence (PL) and TA spectroscopy (for details see Supplementary Section 2). Repeating this analysis for all pump-probe delays, we obtain the time-dependent spectral weight of these contributions as shown in Fig. 2D, E with the same colour scheme from Fig. 2A, B. Figure 2F shows the normalised kinetics of the different states, and the EA/(CT + CS) ratio (shown in black) presents an exponential build-up with a time constant of 40.5 ps. A simple electrostatic simulation of the electric field caused by two charges across a planar D–A interface (see Supplementary Section 5) shows that this monoexponential rise can be assigned to the first jump away from the interface because further increases in the distance between the charges have a negligible effect on the EA signal (see Supplementary Fig. 20). Note that the first jump already creates charge-separated states for which the back-electron transfer to the ground state is strongly suppressed.
It is worth mentioning that charge separation kinetics in BHJ were previously studied employing pump-push photocurrent techniques25,35. However, it was recently shown that the push pulse does not cause enhanced CT dissociation in optimized blends where CT separation is only weakly activated36. The measured 40.5 ps charge separation time suffices to achieve 90% of charge separation yield, given the exciton lifetime of 530 ps of O-IDTBR LE states. This explains the high EQE of over 80% found in WF3:O-IDTBR, a system with very low HOMO–HOMO offset.
An interesting point is the decay of acceptor excitons on the same time scale as the build-up of the CS state (compare green and black curves in Fig. 2F). This is a direct consequence of the LE ⇔ CT equilibrium: the femtosecond pump pulse predominantly generates excitons in the donor phase. Due to the high driving force for exciton splitting, they convert instantaneously into CT states. However, acceptor LE states are accessible from CT states, and therefore an equilibration must occur. If the equilibration time is much shorter than the charge separation time, then both equilibrated species (LE and CT) are depleted on the same time scale by the charge separation process. An alternative explanation, namely diffusion of residual acceptor excitons followed by dissociation, can be excluded because it would not lead to monoexponential decay kinetics.
We note that triplet states have not been included in our model because the quantitative agreement between fit and experiment was achieved without considering them. A recent study has shown that triplet states do form upon CS recombination, however, they are relatively stable against deactivation (the main channel being triplet-triplet annihilation, which occurs predominantly at high excitation fluences)37. We have performed a comparative study using narrowband excitation exclusively in resonance with the acceptor or the donor exciton (Supplementary material, Part 6). These studies are qualitatively in agreement with the results of Fig. 2 albeit with a much longer instrumental response function, due to the narrowband pump pulses. We also observed a possible diffusion limitation of hole transfer after resonantly exciting the acceptor exciton; this has been found recently and ascribed to a lower wavefunction extension in the small molecule NFA system as compared to the donor polymer32. For this reason, we did not use the data from the narrowband experiment for kinetic modelling.
We complement the above experiments with time-dependent density functional theory (TDDFT) calculations to gain insight into the electronic structure of the WF3:O-IDTBR interface, modelled as a dimer formed by one repeat unit of WF3 and one O-IDTBR molecule (for details see Supplementary Section 4). Supplementary Fig. 15 shows the frontier molecular orbital wavefunctions for isolated WF3 and O-IDTBR monomers and for the blend of the two components (panels a,b, and c, respectively). Since the HOMO energies for the isolated WF3 and O-IDTBR are very similar, we find a high degree of hybridization of these orbitals in the HOMO and HOMO-1 of WF3:O-IDTBR. This hybridization transfers significant oscillator strength from the LE state to the CT state so that the CT state acquires 1/5 of the oscillator strength of the LE, a relatively high value. In consequence, the energy of the CT state is also influenced; the calculation yields an energy offset
In order to understand the combination of high EQE and low VOC losses in WF3:O-IDTBR, we performed a multi-objective optimization, aiming to reproduce as closely as possible the three main results of this work (EQE > 80%, EL of blend indistinguishable from that of the pure acceptor, charge separation time 40 ps) by the simplest possible numerical model, taking into account both equilibria LE < = > CT and CT < = > CS and including radiative and non-radiative loss channels from LE and CT states. The model is given in Supplementary Fig. 12.
In our calculation, we consider five free parameters: the exciton breaking rate


Photophysical model deployed for multi-objective optimization.
Black arrows: known values from the literature; blue: values optimized in this work; orange: driving forces for exciton splitting (
The result is shown in Fig. 4. The best match of the three objectives is found at a position marked by a yellow star, suggesting that the driving force for exciton splitting is close to zero. However, the orange trace in Fig. 4a, b indicates regions in which the match of the three objectives is nearly as good as in the global optimum. Therefore, we can only give an upper limit to


Multi-objective modelling.
a IQE and (b) PCE as a function of the driving force for exciton splitting
This result is counter-intuitive at first sight, given the stronger Coulomb attraction of a CT state as compared to a CS state and suggests that charge separation in low driving force systems should be endothermic because these systems lack excess energy needed for coherent long-range phenomena. Recently, temperature-dependent pump-push-probe transient data were presented in low driving force NFA blends, showing very slow charge separation kinetics in the hundreds of picoseconds time domain, supporting the picture of endothermic charge separation39. However, there is an increasing body of evidence that many systems display exothermic charge separation. For example, in similar IDTBR-based blends36, it was suggested that an intermixed donor–acceptor interface can increase the CT energy thus providing a driving force towards the bulk leading to energy-dependent charge dynamics on the picosecond timescale40. Very recently, it has been shown that interfacial quadrupole moments can provide an additional driving force for CT dissociation33,41. Entropy contributions, on the other hand, can be controlled by optimizing the effective interfacial area42. Finally, by charge extraction experiments, charge separation without activation barrier has been demonstrated, another indication that charge separation is an energetically downhill process41.
Although the multi-objective optimization does not yield a precise value for the driving force of exciton splitting, it does explain why WF3:O-IDTBR has high IQE and low voltage losses simultaneously. The area enclosed by the yellow ellipsis in Fig. 4b, which is our best estimate of the combination of driving force and activation energies for WF3:O-IDTBR, is very close to the maximum achievable PCE for each given value of
Our model is generally valid for donor–acceptor blends of the low driving force so that an equilibrium LE-CT is formed. As such, it has a significant impact on material search and device optimization in OPV. Consider the case of Y6, also boasting high IQE and low voltage losses in some of its blends41. This NFA is designed to be particularly rigid, thus combining good charge transport properties with long exciton lifetimes. It can be shown that long exciton lifetimes shift the normalised IQE (orange curve in Fig. 4c) to the left: the material can thus keep unity exciton splitting efficiencies for even lower driving forces38. This means that the maximum of the resulting normalised PCE (green curve in Fig. 4c) will be much broader for Y6 than for O-IDTBR in a given blend, so that for Y6, adjusting the CT energy to an optimum value is less critical for maximising PCE as it is for O-IDTBR. By predicting whether a dedicated morphology optimization will lead to a PCE gain in a given material combination, the application of our model can be used to avoid a tedious trial and error process. In this way, it broadens the range of possible candidate materials by lowering the PCE threshold for consideration.
Both
In conclusion, we have demonstrated a new blend WF3:O-IDTBR showing very low non-radiative losses and at the same time an EQE of more than 80% in the spectral region of the low driving force. Using femtosecond TA spectroscopy, we have shown that the CT state dissociates with a 40.5 ps time constant, which is slower than in previously measured systems with a high driving force, but sufficient to explain an EQE of 80%, despite the low CT decay rate. Numerical simulations using input from quantum chemistry successfully reproduce all experimental observations. We have shown that the reason for the high EQE, despite a low overall driving force, lies in a near-optimal distribution of this energy between the exciton dissociation and charge separation steps. From general considerations of the early CT processes, we propose that for each blend with given LE and CT lifetimes, there exists an optimal value for the driving forces for exciton dissociation and charge separation, that can be adjusted independently from the HOMO level alignment by exploiting interfacial energy destabilization and CT hybridization.
WF3 and O-IDTBR were synthesized as reported elsewhere8.
Fabrication of photovoltaic devices: pre-structured indium tin oxide (ITO) substrates were cleaned with acetone and isopropyl alcohol in an ultrasonic bath for 10 min each. After drying, the substrates were spin-coated with 40 nm of zinc oxide (ZnO) and different active layer based on 20 g L−1 solution in chlorobenzene (CB). To complete the fabrication of the devices 10 nm of MoOx and 100 nm of Ag were thermally evaporated through a mask (with a 10.4 mm² active area opening) under a vacuum of ~1 × 10−6 mbar.
Current density–Voltage (J–V) measurements: the J–V characteristics were measured using a source measurement unit from BoTest. Illumination was provided by a solar simulator (Oriel Sol 1 A, from Newport) with AM1.5G spectrum at 100 mW/cm2. UV–Vis absorption was performed on a Lambda 950, from Perkin Elmer. EQEs were measured using an integrated system from Enlitech, Taiwan. All the devices were tested in ambient air.
Sensitive EQE measurements (EQEPV): EQEPV was measured using a modified Vertex 70 FTIR spectrometer from Brucker optics, equipped with a QTH lamp, quartz beam splitter, and external detector. A low noise current amplifier (Femto, DLPCA-200) is used to amplify the photocurrent produced upon illumination of the photovoltaic device with light modulated by the FTIR. The output voltage of the current amplifier is fed back to the external detector port of the FTIR, in order to be able to use the FTIR’s software to collect the photocurrent spectrum.
EL: EL measurements were performed by using a chopper and applying a constant current (100 mA/cm2) supplied by an external current/voltage source through the devices which have an active area of 0.104 cm2. The emitted light then collected by a monochromator and detected by a liquid-nitrogen-cooled InGaAs detector. The spectrum was recorded by a standard lock-in technique. The system was wavelength calibrated.
TA: the TA setup is fed by a 100-fs, 1-kHz repetition rate Ti:sapphire laser system (Libra, Coherent) with a central wavelength of 800 nm. The ≈10 fs pump pulses span the spectral region from 1.7 to 2.4 eV and are generated by a non-collinear optical parametric amplifier (OPA), compressed to the near-transform limited duration at the sample position with a pair of chirped mirrors43. The probe pulses consist of a white-light continuum (WLC) generated by tightly focusing the output of an OPA at 1225 nm on a 2 mm thick sapphire plate, generating a WLC extending from 1.3 to 2.2 eV, limited on the NIR side by the short pass filter used to reject the residual 1225 nm OPA light44. Measurements are performed in transmission using an SP2150 Acton, Princeton Instruments spectrometer equipped with a CCD detector. The differential transmission (ΔT/T) spectrum is acquired as a function of probe photon energy and pump-probe delay, and pump and probe beams are collinearly polarized.
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
The online version contains supplementary material available at 10.1038/s41467-021-22032-3.
N.G. acknowledges the Imperial College Research Fellowship scheme. C.J.B. gratefully acknowledge funding from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation), Project no. 182849149-SFB 953. C.J.B. gratefully acknowledges financial support through the “Aufbruch Bayern” initiative of the state of Bavaria (EnCN and SFF) and the Bavarian Initiative “Solar Technologies go Hybrid” (SolTech) and funding from DFG project DFG INST 90/917. C.C.L. thanks the European Union for the financial support. G.C. acknowledges the support from the PRIN 2017 Project 201795SBA3—HARVEST.
N.G. and C.J.B. conceived and developed the ideas. N.G. designed the experiments and performed device fabrication, electrical characterisation and data analysis. N.G. and M.S. performed EL, PL and FTPS measurements. F.V.A.C. performed TA measurements under the supervision of T.N., L.L. and G.C. S.F. performed the molecular quantum-chemical calculation under the supervision of A.G. L.L. developed the modelling of the TA, EQE, EL spectra. A.C. performed transient PL measurements. C.C. synthesized WF3. S.R. performed EL measurements under the supervision of D.N. A.W. synthesized the IDTBR acceptor under the supervision of I.M. N.G., F.V.A.C., and L.L. wrote the manuscript. M.S., D.N., C.J.B., D.B., V.G. and D.N. contributed to the revision of the manuscript. The projects were supervised by G.C., C.J.B. and L.L.
The data that support the findings of this study are available on request from the corresponding author (N.G.).
The authors declare no competing interests.
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
29.
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
40.
41.
42.
43.
44.