Free charge photogeneration in a single component high photovoltaic efficiency organic semiconductor

Nature Communications volume 13, Article number: 2827 (2022) Cite this article

8811 Accesses

37 Citations

42 Altmetric

Metrics details

Organic photovoltaics (OPVs) promise cheap and flexible solar energy. Whereas light generates free charges in silicon photovoltaics, excitons are normally formed in organic semiconductors due to their low dielectric constants, and require molecular heterojunctions to split into charges. Recent record efficiency OPVs utilise the small molecule, Y6, and its analogues, which – unlike previous organic semiconductors – have low band-gaps and high dielectric constants. We show that, in Y6 films, these factors lead to intrinsic free charge generation without a heterojunction. Intensity-dependent spectroscopy reveals that 60–90% of excitons form free charges at AM1.5 light intensity. Bimolecular recombination, and hole traps constrain single component Y6 photovoltaics to low efficiencies, but recombination is reduced by small quantities of donor. Quantum-chemical calculations reveal strong coupling between exciton and CT states, and an intermolecular polarisation pattern that drives exciton dissociation. Our results challenge how current OPVs operate, and renew the possibility of efficient single-component OPVs.

When light is absorbed by inorganic semiconductors, like silicon or gallium arsenide, a combination of free charges and excitons (bound electrons and holes) are created, and their ratio depends upon the material’s dielectric constant. For organic semiconductors, the paradigm holds that only excitons are intrinsically photogenerated because the low dielectric constant1,2 of these materials (ε ~ 3–4) means the electron and hole Coulombic interaction is not efficiently screened, leading to Frenkel excitons with high binding energies (EB). For the last thirty years3,4, this constraint has guided the development of organic photovoltaics (OPVs) - which promise a step-change in flexible, lightweight, non-toxic solution-processed solar energy production, but are yet to be widely commercialised. Splitting bound excitons into free charges has required sharp molecular heterojunctions between donor and acceptor materials. Rather than purely optimising charge harvesting, OPV devices have been optimised for exciton splitting at interfaces, and demand complex interpenetrating networks of donor and acceptor materials. The molecular heterojunction approach limits device efficiencies5, introduces inherent voltage losses and instabilities from interfaces6, and complicates research progress7. Reports have previously shown some field-assisted exciton dissociation, and ‘extrinsic’ charge formation in neat homopolymers and small molecules8,9,10,11, as well as significant charge-transfer (or polaron-pair) state formation in neat “push-pull” copolymers and small molecules12,13,14. However no organic material (until now) has exhibited substantial, let alone majority, free charge (rather than polaron pair) formation – which undergoes bimolecular (rather than geminate), recombination.

Recently, small molecule non-fullerene fused ring electron acceptors (FREAs)15,16,17,18,19,20,21 have driven a rapid uptick in the power conversion efficiency (PCE) of OPVs. Alongside this advance has arisen unexpected observations, namely: (1) Barrierless free charge generation in PM6:Y6 blends22 (full chemical names and structures shown in Supplementary Fig. 1). (2) Charge generation efficiency in blends of PM6:Y6 increases with increasing incident light intensity23. (3) In Y611,24 and IDIC25, excitons are delocalised, or form an ‘intra-moiety’ intermediate state with likely charge-transfer (CT) like character. These observations have all been explained by invoking CT-states, however, with additional data, we show here that a more profound explanation is required.

FREAs have recently been measured to have very high refractive indices26,27, and hence high complex dielectric constants. Of the FREAs, Y6 and its derivatives are present in both the binary and ternary OPVs that hold the highest PCEs16,28,29,30.

In inorganic semiconductors, the number of bound excitons versus free charges can be described by the Saha-Langmuir relation.

where the free charge fraction, \({x}\), is dependent upon the temperature, T, excitation density, n, effective mass, \({m}^{* }\) and the exciton binding energy,\(\,{E}_{B}\) (which is related to the inverse of the dielectric constant31). From comparison of the highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) gap to the optical gap measured by Karuthedath et al.6, the exciton binding energy of Y6 is EB = 100–250 meV - much lower than, for example, P3HT (EB ~ 700 meV)32. Based on a midpoint, 175 meV binding energy, estimated for Y6 above, and assuming an effective mass between33,34 \({m}^{* }\)= 1.7–0.2 \({m}_{e}\), this simple calculation would predict a high fraction (20–80%) of intrinsically generated free charges.

Here, we show that the high optical frequency dielectric constant of Y6 leads to majority free charge generation upon optical excitation within neat films. To prove that the aforementioned ‘intra-moiety’ TA signature is in fact due to free charges, we demonstrate that recombination follows bimolecular kinetics, rather than monomolecular decay as would occur for CT states. We use a combination of ultrafast transient absorption, photoluminescence up-conversion and time-resolved terahertz spectroscopy, and identify the spectral signatures and decay rates of singlet excitons and free charges (polarons). Intensity-dependent photoluminescence (PL) measurements prove the existence of radiative bimolecular charge recombination, which would not be the case if only CT states were formed. We present simple kinetic models that accurately reproduce our experimental data, wherein a free charge fraction of between 0.7–0.99 exists in equilibrium with bound charges under steady-state illumination, and 60–90% of photoexcited excitons dissociate to charges. Bimolecular, and rapid minority carrier, charge recombination inhibit the efficiency of single-component devices, but this recombination can be slowed by ‘doping’ our Y6 films with only a small amount of ‘p-type’ polymer electron donor. Using quantum-chemical calculations, we show that free charge generation is expected due to strong coupling between equienergetic exciton and CT states, and reveal that asymmetric polarisation of the different lattice positions results in a “donor/acceptor” system derived purely from intermolecular interactions.

Fig. 1a shows the normalised excitation density-dependent external photoluminescence quantum efficiency (PLQE) of two neat films of Y6 (BTP-4F) of differing thicknesses. As the intensity of pulsed excitation is increased, the PLQE increases steadily, peaks, and then decreases for both films. This peak in PL efficiency has not been observed before in similar organic systems, where it is expected that the PL remains steady with increasing fluence before dropping in efficiency as non-radiative bimolecular processes kick-in. Three explanations for the phenomena of a rise in PL with increasing pump intensity are possible: saturation of exciton (or charge) traps, triplet-triplet annihilation (TTA), or bimolecular radiative recombination of intrinsically generated free charges (bulk ionisation).

a Normalised photoluminescence efficiency, as a function of excitation density, of 130 nm thick (blue dots) and 80 nm thick (orange dots) neat films of Y6 on a glass substrate. The sample was excited by 600 ps pulses of 532 nm light. The sample was measured in vacuum, and the relative photoluminescence efficiency was determined from the ratio of the 2nd harmonic scatter peak at 1064 nm to the integrated PL intensity between 850 and 1050 nm. The error bars shown represent the relative error of the change in PL at different excitation densities, rather than the total error (details in the Supplementary Information). b Simulation of intensity-dependent PLQE values showing effect of increasing (orange) or decreasing (blue) radiative bimolecular recombination (from prompt free charge fraction, calculated using the model detailed in the text and SI. c Graphic representation of exciton-free charge equilibrium in Y6.

Exciton trap saturation is considered in detail in the Supplementary Information. We rule out this explanation due to the unphysical rate constants required for it to match the transient absorption and PL data, and due to the lack of any other evidence of exciton trap saturation in the TA spectra/kinetics. TTA is also considered in the Supplementary Information. Due to a number of considerations, such as the large numbers of triplets required for this effect to be solely responsible for the PL rise, which conflicts with experimental evidence, we also rule this out. We therefore conclude that a substantial fraction of photogenerated free charges must be responsible for the PL rise shown in Fig. 1a. As detailed in the below sections, and the Supplementary Information, charge trap saturation, and TTA are possible, and likely present to some extent in our system, but only in conjunction with significant intrinsic charge generation.

Figure 1b illustrates the expected behaviour of excitation density-dependent PLQE for varying proportions of radiative monomolecular to bimolecular processes. The rise in efficiency observed in Fig. 1a is not possible for a system undergoing monomolecular (geminate) radiative recombination, even from intra-moiety CT states, and exciton-exciton annihilation can only produce a bimolecular contribution with reduced radiative efficiency. Charges must be present in the neat film, and these charges are formed ‘intrinsically’, rather than through exciton annihilation processes. Under steady-state illumination, this free charge (polaron) population will exist concurrently with a population of bound Frenkel or CT-like excitons, as illustrated by Fig. 1c.

To quantify the relative proportions of charges versus excitons, we perform ultrafast spectroscopic measurements on films of neat Y6. We identify spectral signatures of the singlet exciton, and free charge/CT-state species and their kinetics to estimate first and second order kinetic rates/lifetimes. Figure 2a shows the area-normalised transient absorption (TA) spectrum of a neat Y6 film pumped with a 150 fs, 700 nm laser pulse at a moderate excitation density of ~1017 excitations/cm3. The spectrum evolves from one initial photoexcited species, a sharp negative peak at 915 nm, to another consisting of both positive and negative features, apparent after ~500 ps. Singular value decomposition, multivariate curve resolution analysis and application of a genetic algorithm (GA)35 consistently resolve two major species present, shown as the blue and orange spectra in Fig. 2b. We confirm that the orange species is that of the polaron – referred to as free charge (FC) for consistency in the rest of the manuscript – by comparing the TA spectra of Y6 blended with two hole accepting species (PTB7-Th and poly-TPD), to reveal the Y6 electron signature. The charge spectrum in neat Y6 closely resembles charge spectra obtained from these measurements (we note that the donor hole signature is also present in the PTB7-Th:Y6 measurement, which must be accounted for when making this comparison). We note that the charge is also likely to have a very similar spectrum to an intraspecies charge-transfer36 state/excimer/weakly-bound exciton delocalised over 2 or more molecules in the Y6 crystal, and therefore designate the spectra here as representing the free charge plus CT state populations. We stress that we cannot explain all of our data unless there is a significant fraction of free charges, as detailed below.

a Transient absorption spectra, normalised to the absolute area integral, of a neat film of Y6 at different time delays. 700 nm, 200 fs pump pulse, excitation density ~ 1016 excitations/cm3. The spectrum evolves from the exciton, blue line, with a negative peak at 915 nm, to the red-shifted charge/CT state spectrum, orange line. b Confirmation of charge and exciton transient absorption spectra by comparison of neat species (extracted using genetic algorithm) with species resolved from blends of Y6 with hole accepting materials (PTB7-Th, polyTPD). The green dashed line shows the hole spectra of PTB7-Th blended with PCBM. c Singlet component from GA of transient absorption (dark blue line), at an excitation density of ~8 × 1016 cm−3, compared with up-converted PL signal (cyan), at an excitation density of ~2 × 1016 cm−3. d Initial kinetics of the exciton (blue) versus charge/CT (orange) TA kinetics in the first 2.5 ps after photoexcitation. The charge species shows a rise concomitant with the fast initial decay of the singlet species. 550 nm, 200 fs pump pulse, excitation density ~ 1018 excitations/cm3.

We confirm that the blue species is that of the singlet exciton in 2 ways: (1) through TA measurements of a ‘solid-solution’ of Y6 diluted in polystyrene (Supplementary Fig. 2), and (2) by comparing the TA spectra of the ‘blue’ species to the transient kinetics of upconverted PL signal from the same film at roughly similar excitation density. As shown in Fig. 2c, the singlet TA matches the transient PL signal. The singlet, and free charge/CT species show different kinetic decays to each other across multiple different excitation density regimes (see below). As an initial indication of bimolecularity, we note that the square of the charge population closely matches the derivative of the charge, even at low fluence.

The singlet and charge species, as shown in Fig. 2d, show a prompt, fluence independent (Supplementary Fig. 4), interconversion within the first 2.5 picoseconds, which also has no clear excitation-energy dependence (observable within our excitation energy resolution – Supplementary Fig. 5). These observations are consistent with observations from Wang et al.11, though our other measurements (terahertz spectroscopy, ultrafast transient PL, and fluence dependent PLQE) necessitate a model beyond pure CT generation. As this effect appears not to be due to excess energy (based on TA measurements pumped with different excess energies), the rapid formation is likely an effect of increased exciton delocalisation at early times37,38. In line with Gillet et al.39, we also observe a small 3rd spectral component in the TA which is consistent with a triplet signature (Supplementary Fig. 6). As further evidence that there are significant free charges present, the triplet signal peak is significantly delayed compared to the charge/CT signal, which is consistent with a non-geminate triplet generation pathway rather than geminate triplet formation from inter-system crossing. Additionally, optical-pump-terahertz-probe spectroscopy (Supplementary Fig. 17) shows a conductivity spectrum consistent with a fraction of free charges present at ~3 ps pump delay time (although we note that this spectrum could also be consistent with a pure CT-state population – insufficient signal-noise ratio in the measurement prevents more precise interpretation).

We fit our transient absorption kinetics and intensity-dependent PLQE data based on a basic kinetic model of free charges and excitons present in a neat film of Y6 upon illumination. Fig. 3a graphically represents this model (with further details in SI and Fig. 3b, c). The key processes illustrated and modelled are: Photogeneration of an initial, delocalised singlet exciton, S1*, which rapidly evolves into localised excitons, S1, and free electrons and holes, FC. Free charges bimolecularly recombine to form singlets and triplets, or non-radiatively via trap/defect states. This basic model is all that is required to accurately recreate our experimental data. We can optionally include further processes, such as triplet-triplet annihilation, or specify the traps as most likely being hole traps40, to improve our fits to the data further. Our data fails to fit to a model that explicitly includes a substantial and long-lived CT state population (see SI for further details).

a Jablonski diagram highlighting key excited-state pathways – rapid initial interconversion (solid red arrows), charge recombination (solid orange arrows), non-radiative (red dashed arrows), and radiative rates (green arrows). b Normalised transient absorption kinetics of excitons (blue circles) and (c) charge states (orange circles) at different excitation densities, fitted with a global fit to the basic kinetic model described in the text and SI (solid blue and orange lines). d Internal (blue crosses) and external (red crosses) PLQE values of Y6 as a function of excitation density, with corresponding simulated values from fits to the transient absorption and intensity dependent PLQE. Internal PLQE is calculated from external PLQE as per Richter et al.41. The solid blue line is the simulated PLQE from the basic model, the dashed blue line shows the PLQE calculated when TTA is included, and the dotted blue line shows the estimated PLQE from an explicit treatment of hole traps (but no TTA). e Steady-state free charge fraction as a function of total excitation density calculated from the rate constants gathered from the basic kinetic model (red line), and from an estimate from the Saha equation (blue line). The Saha equation is calculated based on an excitonic binding energy of 270 meV, corresponding to an exciton effective mass of ~0.5 me.

Figure 3b–d show the results of a global fit of this basic model to a multiple fluence series of transient absorption, and PLQE data (with TTA and trap effects also shown in Fig. 3d). We calculate internal PLQE41 (taking photon recycling into account – see below) as the ratio of radiative decay events to total decay events, integrated over a microsecond, with a simulated 600 ps input photon pulse. The radiative and non-radiative singlet decay rates are obtained from comparison of the PLQE values of Y6 in a solid-state solution of polystyrene.

These models simultaneously fit the transient absorption data, and crucially, fit the large increase, or ‘hump’, in PLQE with only a small number of free parameters (as low as 4). The fitted parameters correspond to physically meaningful rate constants in agreement with literature. We thus view them as instructive for estimating free charge yields and transition rates in Y6. Error analysis of these models is shown in the Supplementary Information. Though there is interdependency between the fitted rates, and hence the individual error in each rate constant is significant (and model dependent), we can use these models to put a lower and upper bound on steady-state free charge yields at given excitation intensities. Setting the conservative criterion that there is a 10% rise in PLQE with excitation intensity, the models show that the exciton dissociation probability (the proportion of excitons that dissociate to charges within their lifetime) lies in the range 0.6–0.9, and the corresponding steady-state free charge fraction lies in the range 0.65–0.99 at 1 Sun excitation intensity. Our total carrier density under AM1.5 lies in the range 1013–1014 cm−3 depending on the model used.

The red line in Fig. 3e illustrates the fraction of steady-state free charges present based on a basic kinetic model with no TTA or hole trap filling. Almost completely overlapping this line is a blue curve given by the Saha equation, with an excitonic binding energy of 270 meV (corresponding to an exciton effective mass of ~0.5 me). The very close match between the two curves shows the equivalence of the kinetic model to the Saha equation.

As illustrated in Fig. 3d, it is important to take into account photon reabsorption in PLQE measurements of Y6 thin films. This is a further consequence of their extremely high refractive index. We use the same process as Richter et al.41 to approximate photon reabsorption, and show that the true PLQE of Y6 is much higher (~6%) than what would be ascertained from a direct reading of the external PLQE (~2%). This has significant importance for solar cell material optimisation efforts, as maximising PLQE is a useful approach for assessing photovoltaic efficiency potential42, and ultimately, being able to harness radiative carrier recombination through photon recycling is key to approaching the Shockley-Quiesser limit. Although we estimate photon recycling is small for these measurements (Supplementary Fig. 16), the high refractive index means that it will become significant if internal PLQE can be increased.

High intrinsic charge yields suggest a new path for OPV design by creating single component devices, or devices with very low donor contents (which may also be beneficial for semi-transparent OPV). However, the results from the transient absorption, PLQE, and terahertz spectroscopy data explain why, even though charges are intrinsically generated, single component, monolayer devices are not efficient. High bimolecular charge recombination, combined with a significant hole trap population, mean that free charges recombine before arriving at electrodes in a device. With our measured charge recombination rates, and assuming (for the sake of illustration), our monomolecular charge recombination term to give the minority carrier lifetime, the hole diffusion length would be less than 15 nm. We fabricated single component Y6 devices. The highest PCE based on a ITO/PEDOT:PSS/Y6/LiF/Al structure was 0.09%. By replacing PEDOT:PSS with PCP-Na, whose HOMO better aligns with that of Y6, the highest PCE obtained was 0.63% (Supplementary Fig. 7). This value is low (though there is no junction in the active layer to separate charges and prevent minority carrier recombination). This loss of photocurrent and photovoltage is due in part to high bimolecular, and space-charge induced recombination, which we show by performing intensity dependent short-circuit current (JSC) measurements. Previous reports have seen power law dependence of JSC on illumination intensity (I) according to:43,44 \({{{{{J}}}}}_{{{{{{\rm{SC}}}}}}}\propto {I}^{\alpha }\). Deviations of the fitting parameter, α, from unity signify bimolecular recombination or space-charge effects due to unbalanced charge mobilities43,44,45,46. Electron and hole mobilities are relatively balanced in Y6 (being 1.8 − 40 × 10−4 and 0.5 − 50 × 10−3 cm2 V−1 s−1 respectively24,47), and at intensities below 0.01 suns, the exponent, α, approaches 1 (Supplementary Fig. 8). Fig. 4a shows that for the region between 0.2 and 2 suns incident photo-intensity, the fitted exponent is α = 0.84, indicating the device suffers from significant bimolecular recombination and/or space-charge effects at 1 sun.

a Intensity dependent JSC of a monolayer Y6 device. The power exponent fitted to the data between 0.2 and 2 suns is 0.84, indicating a high number of bimolecular recombination events. b External quantum efficiencies (EQEs) of devices with active layers consisting of small ratios of donor, PTB7-Th, to acceptor, Y6. For such low donor content, EQEs are unexpectedly high. c Normalised transient absorption kinetics of exciton and charge species in blends of PTB7-Th:Y6 with ratios of 1:50 and 1:100. Pumped with 800 nm, 150 fs pulses at excitation density 1 × 1017 excitations/cm3. The exciton kinetics of both films, dashed lines, are matching, as are the rise times of the charge kinetics (orange and red lines), but the 1:50 blend shows slower charge recombination after 20 ps. The green-blue line shows the hole kinetics in the 1:50 blend, showing a gradual increase in hole population as excitons and charges in Y6 diffuse to the PTB7-Th.

To further emphasise the importance of high charge recombination, rather than exciton splitting, as a limiting factor in device efficiency, we investigated the effect of ‘doping’ our Y6 material with very small amounts of donor polymer, PTB7-Th48. These very dilute blends allow us to temporally separate the processes of charge separation, transport, and recombination in transient absorption measurements. We measured the photovoltaic external quantum efficiencies of the corresponding devices, which both show surprisingly high efficiencies. As shown in Fig. 4b, when blended in a 1:50 ratio of PTB7-Th:Y6, the EQE is nearly twice as high as for the 1:100 parts blend. However, Fig. 4c shows that this is not due to differences in exciton splitting. The transient absorption kinetics show that for both blends the exciton species decay, and the charge species appear, with almost identical kinetics in the first 20 ps, but the charges recombine faster in the 1:100 blend, showing that the polymer donor’s key role is to decrease the rate of recombination. For both bimolecular charge recombination, and hole-trap (or n-dopant) induced recombination, once holes are removed from the Y6 - as they are, rapidly, in normal bulk heterojunction (BHJ) blends (as shown in Supplementary Fig. 10) - then the left-over electrons have nothing with which to recombine, and their recombination kinetics are slowed after electron and hole separation.

Density functional theory calculations were used to gain an understanding of the structural features responsible for free charge generation in Y6. The first requirement is significant electronic coupling between S1 and CT1. This parameter was calculated for all molecular pairs exhibiting π-π stacking interactions in the Y6 crystal structure24 by projecting the quasi-diabatic exciton/CT states calculated for a pair of separated molecules onto the adiabatic states of the molecular pair in the crystal geometry49,50. This approach enables us to simultaneously predict the energies and couplings for exciton/CT states (Fig. 5a, Supplementary Fig. 22, and Supplementary Table 1).

a Exciton-CT energy offsets and electronic couplings (in meV) for selected π-stacked molecular pairs extracted from the Y6 crystal structure (alkyl chains have been removed for clarity, remaining structures can be found in the Supporting Information)24. b Density of states (DOS) for holes (−IE) in a 10 nm thick model thin film based on molecular dynamics equilibration of the Y6 crystal structure, normalized to unit area. Dashed line indicates the DOS onset determined by photoelectron spectroscopy in air6. c Average ionisation energy as a function of π-stacking layer showing differential hole stabilisation for the different lattice positions, which occur in an “…ABBA…” sequence.

The computed couplings and energetics point to a kinetic pathway for exciton break-up. The exciton-CT coupling is significant (10–75 meV), with the CT states either equal or lower in energy to the excitonic states (except one pair in which the stacking interaction is minimal). These two factors, combined with the low reorganisation energy (estimated as ~0.3 eV based on the individual molecules) imply that exciton dissociation is expected, even according to incoherent charge transfer theories. The same electronic coupling parameters/state energies are also consistent with the possibility of coherent CS, as we observe spectroscopically. Indeed, as previously noted, the coupling between excitonic states suggests a high degree of exciton delocalisation24, which is expected to lower the reorganisation barrier for charge formation. Coupling of the CT states with the ground state is also strong (30–85 meV), which is consistent with the rapid recombination observed experimentally. We note that the CT-GS coupling strengths do not correlate with the exciton-CT coupling in a simple manner. This apparent independence of VEx-CT and VCT-GS indicates that the recombination can, at least in principle, be overcome by crystal structure engineering to maximise the VEx-CT/VCT-GS ratio. Evidence of exciton/CT hybridisation may be seen in the broader red tail of the UV-vis absorption spectra.

Our calculations also offer a clue regarding the origin of charge formation in Y6. In pairs involving molecules from different lattice positions, there is a clear energetic preference (~0.15–0.3 eV) for which of the two possible CT states is formed. These observations indicate that the packing geometry results in two distinct polarisation environments, creating a donor-acceptor system through supramolecular, rather than synthetic means. This interpretation also explains the significant recombination observed experimentally: the alternation of 'donor' and 'acceptor' positions on the scale of the unit cell is akin to an extremely intermixed blend system – meaning that both charge formation and recombination are extremely rapid.

The energy levels calculated for free charges also reveal the presence of distinct polarisation environments. Fig. 5b shows the ionisation energy (IE) distributions calculated for a Y6 thin film using the long-range electrostatic embedding procedure available in VOTCA (electron affinity (EA) is shown in Supplementary Fig. 13)51,52,53,54,55,56,57. Figure 5c shows that the different energetics of charges occupying distinct crystal positions results in a bimodal density of states, particularly for holes. These distributions confirm the presence of a driving force for charge formation inherent to the Y6 packing structure, on the order of the exciton binding energy. This situation bears some similarity to single component OPVs in which crystalline domains with different orientations possess distinct charge energetics, resulting in charge formation at interfacial boundaries58. The present work is distinguished from the prior literature in that the different energetics are inherent to the crystal structure, resulting in bulk - rather than interfacial - charge generation.

In summary, our calculations reveal equal energy exciton-CT manifolds (as a result of differential polarisation), strong exciton-CT coupling, and a flat energetic landscape for electron transport. Taken together with previous theoretical work showing appreciable electronic coupling for electron and hole transfer, these observations strongly suggest the existence of a kinetic pathway for exciton-to-free charge conversion in Y6.

The observation of intrinsic charge generation in a single organic material necessitates an expanded understanding of the photovoltaic process in organic materials. The implications for device design are far-reaching, showing that a new approach to reaching maximum efficiency may be to maximise intrinsic charge generation yield, whilst reducing rapid trap/dopant-assisted monomolecular and bimolecular charge recombination, rather than focus on exciton splitting in a donor-acceptor blend. Indeed, we suggest that the highest efficiency OPV blends (such as PM6:Y6) have in fact been optimised for the extraction of intrinsically formed charges.

We propose a mechanism for how recombination is suppressed in current Y6 blends, drawing on studies invoking energy level bending at the donor/acceptor interface as a result of long-range interactions of charges with molecular quadrupole moments6,22,51,59,60,61. In an excitonic picture – where charges only form at the interface - these fields assist CT dissociation and suppress non-geminate recombination. However, these benefits come at the expense of exciton dissociation rates, resulting in an intrinsic limit for the ionisation energy offset between the neat materials6. When free charges are intrinsically generated in Y6, not only is the interfacial exciton dissociation barrier bypassed, but the level bending defines an attractive potential for holes, which are pulled towards the interface before leaving the acceptor material, while electrons are repelled. Under this interpretation, champion efficiency systems, such as PM6/Y6, benefit from both exciton- and charge-funnelling effects: excitons in PM6 are funnelled to Y6 via resonant energy transfer6,27, while level bending pulls holes from Y6 into the polymer domain, and suppresses recombination thereafter. We can test whether quadrupolar fields affect charge recombination in Y6 by comparing charge and exciton kinetics of neat Y6, to a film of Y6 blended with 20% weight ratio PCBM (a ratio present in highly efficient ternary devices28). The energy levels of PCBM mean that both electron and hole transfer from Y6 to PCBM are prohibited. Any reduction in charge recombination rate will therefore be due to interfacial level bending (one would expect the effect of increased crystallinity in Y6 due to the PCBM would increase recombination, rather than decrease). Supplementary Fig. 11 shows that, as predicted, there is a small reduction in charge recombination rate, and concurrent reduction of exciton population, when PCBM is introduced.

Efficient intrinsic charge generation, in the FREA small molecule Y6, has broad implications for OPV device design. While previous works (on materials such as P3HT62 and PCBM63) have generated debate, showing evidence of a proportion of (10–15%) charges generated in neat films, this work differs significantly in that the free charges are not formed ‘extrinsically’ (in the case of P3HT64), or from field-induced exciton separation10. Considering a population of excitons and free charges in neat Y6 necessitates a more nuanced understanding of the simplified four step process – exciton generation, exciton transport, charge separation at a heterojunction, charge transport – that has previously dominated OPV discussion. Exciton diffusion length is no longer a simple function of how fast a lone exciton travels and how long it lives - interconversion into free charges and the reverse process must also be considered. The ease of charge separation, but concurrent rapid charge recombination, necessitate a renewed focus on reducing such recombination. This hints that for current Y6-based devices the key purpose of the bulk heterojunction is more to reduce charge recombination, than to split excitons. Emphasis may shift from not only reducing interfacial and CT state bimolecular recombination, but to also focussing on a reduction in minority-carrier recombination. Simultaneously maximising intrinsic PLQE and charge generation should improve OPV device open circuit voltages – which have yet to achieve parity with inorganic photovoltaics. Specific directions of study include: improving crystal packing to further enhance dielectric constants and hence free charge fraction; tuning energy levels through packing structures and/or molecular orientation rather than different compounds; and crystal structure engineering to minimise CT-ground state coupling.

We have shown, through intensity-dependent PLQE measurements, that free charges exist in neat films of the molecular semiconductor Y6, at illumination intensities below 1 sun. We have corroborated this measurement with optically-pumped time-resolved terahertz, ultrafast photoluminescence up-conversion and transient absorption spectroscopy. Our data can be explained by simple kinetic models, along with a new quantum mechanical picture of exciton-CT state coupling and free charge energetics. We have outlined the implications of this finding for our understanding of the photovoltaic process, and the mechanism of action for current champion efficiency devices, and have suggested further avenues of study for enhancing efficiency in future devices. These findings open up the potential to think beyond the bulk heterojunction, and revisit the possibility of fabricating efficient doped organic p-n junctions.

UV-Vis absorption spectra were obtained using Carry 5000 (Agilent) and a Cary 50 Bio UV-vis spectrometer in the range 190–1100 nm. Photoluminescence spectra were obtained using a Cary Eclipse (Varian). Photoluminescence quantum efficiency measurements were obtained in an integrating sphere using the method of de Mello et al.65. A Q-switched frequency doubled ND:YAg laser with output wavelength 532 nm, pulse-length 600 ps, and rep-rate set to 25 kHz was used to excite the sample for the intensity dependent PLQE measurements. Sample film thicknesses were performed using a Dektak profilometer. Thin films of Y6, and Y6 blends were spin-coated onto quartz spectrosil substrates in a glovebox from 11 mg/ml chloroform solutions, at a spin-speed of 3000 rpm. Solid solutions of Y6:polystyrene were prepared in a 1:50 weight ratio, and spincoated in a glovebox at 800 rpm.

Indium tin oxide-patterned glass substrates were cleaned by sequential sonication in acetone and propan-2-ol and subjected to ozone treatment. A 35 nm film of PEDOT:PSS was spin coated in air from aqueous solution and baked on a hotplate at 110 °C for 15 min. A 30 nm layer of Y6 was spin coated from a chlorobenzene solution in a nitrogen-filled glovebox (0.1 ppm O2, 0.0 ppm H2O). The cathode of 1 nm LiF and 90 nm Al was deposited by vacuum thermal evaporation at pressure ~10−6 mbar through a shadow mask to define 6 mm2 pixels. The devices based on PTB7-Th:Y6 were fabricated with a normal structure as ITO glass/PEDOT:PSS/PTB7-Th:Y6/PNDIT-F3N/Ag. Patterned ITO glass was precleaned in an ultrasonic bath with deionized water, acetone and isopropanol, and treated in an ultraviolet–ozone chamber (Jelight Company, USA) for 20 min. PEDOT:PSS was spin-coated on the pre-cleaned ITO at 5000 rpm, followed by baking at 150 oC for 15 min. Then, the active layers were spin-coated on PEDOT:PSS (PTB7-Th:Y6, X:12 mg mL−1 in chloroform, 2500 rpm). Afterwards, the PNDIT-F3N solution (0.5 mg mL−1 in methanol) was spin-coated on the active layer at 2000 rpm. Finally, Ag electrode (ca. 80 nm) was slowly evaporated onto the surface of the underneath layer under vacuum (ca. 10–5 Pa). The devices were not masked and the active area of devices were 4 mm2.

Photovoltaic device JV characteristics were measured in a nitrogen-filled glove box at room temperature using a solar simulator and source-measure unit. Light intensity was determined using a calibrated silicon diode. We tested 20 devices for each kind and choose the best performing devices. Light intensity was varied by adjusting the distance between source and device under test and also by use of neutral density filters. EQE spectra of PTB7-Th:Y6 devices were obtained using a Solar Cell Spectral Response Measurement System QE-R3011 (Enlitech Co.).

Measurements were performed using a homebuilt experimental setup illuminated by an amplified Ti:sapphire laser, with pulse durations of 100–150 fs, centred around 800 nm and at a repetition rate of 3 kHz. The excitation pulses are generated either using this fundamental, or by using an optical parametric amplifier (TOPAS) with the 800 nm fundamental input and then chopped at \(f/2\) (1.5 kHz). Photoexcitations in the materials were probed via a broadband white light continuum generated by focusing a portion of the fundamental to an undoped YAG (Yttrium Aluminium Garnet) crystal. Pump-probe polarizations were kept under magic angle (54.70) configuration to avoid orientational dynamics. After passing through the photoexcited sample, the probe pulses were spectrally dispersed using a prism spectrometer and are then collected using a CMOS camera (visible components) or an InGaAs photodiode array (IR components). The time resolution is obtained via introducing time delays in the pump path which is achieved using a retroreflector connected to a motorized translational stage. The differential transmission signals at various time delays are calculated from the sequential probe shots corresponding to the pump on versus off. For typical measurements, 8000 shots were averaged at each time point and were repeated at least four times. To avoid degradation, all the samples were measured under a vacuum, or inert nitrogen environment.

Optical pump terahertz probe spectroscopy was performed using a dual lock-in technique, similar to that described by Tiwana et al.66. THz pulses are generated by an amplified Ti:sapphire laser, with pulse durations of 100–150 fs, centred around 800 nm, generated at a repetition rate of 3 kHz chopped to 1.5 kHz, incident onto a 1 mm thick ZnTe crystal. THz pulses are measured by balanced photodetectors using electrooptic sampling in a separate ZnTe crystal. The 800 nm ‘optical’ pump is chopped at 750 Hz, and incident on the sample with a spot-size diameter of 4–5 mm. The THz probe diameter was measured as ~1–2 mm. The entire THz beam path is enclosed in a cavity pumped with dry air, with humidity monitored to ensure % water remains below 2%. The sample is enclosed in a nitrogen environment to avoid laser degradation and is excited from the glass substrate side. THz conductivity was calculated66 as:

where \({\varepsilon }_{0}\) is the vacuum permittivity, c is the speed of light, nA and nB are the refractive indices of the media on both sides of the sample in the THz region – quartz, and air, L is the sample thickness, and \(\frac{\triangle T}{T}\) the measured change in THz transmittance.

Ultrafast photoluminescence dynamics were measured using a photoluminescence up-conversion technique. The setup’s light source is a Ytterbium fiber laser (Tangerine SP, Amplitude Systemes) operating at 44 kHz and generating 150 fs pulses. The laser output was split into pump and gate parts. The pump was frequency-doubled to 515 nm using a BBO crystal and focused to a 50 µm2 spot to excite the sample. The sample fluorescence was collimated and refocused on a nonlinear mixing crystal (1 mm BBO, θ = 32°) by a pair of aluminium-coated off-axis parabolic mirrors. The crystal angle is set for the phase-matching condition (type I) to produce sum frequency signals at 480 nm from mixing the 1030 nm gate and 900 nm photoluminescence. Two achromatic lenses collimated and focussed the upconverted signals onto the spectrometer (Princeton Instruments SP 2150). Signals were detected by an intensified CCD camera (Princeton Instruments, PIMAX3). A combination of 550 nm long‐ and 800 nm short‐pass filters were used to remove residual excitation and 1030 nm gate light, respectively. The time delay between pump and gate beams was controlled via a motorized optical delay line on the excitation beam path. For each delay time, 240,000 shots were accumulated. Samples were measured under a vacuum environment.

Density functional theory calculations were performed using Gaussian 0967 with an ultrafine integration grid, and the accuracy of two-electron integrals increased to 1011. Optimized geometries were confirmed by calculation of the associated vibrational modes, which revealed no imaginary frequencies in all cases.

Electronic coupling calculations were conducted on pairs of molecules extracted from the X-ray crystal structure, with alkyl chains truncated to ethyl. These calculations used the long-range corrected CAM-B3LYP68 exchange–correlation functional in combination with the 6–31+G(d,p) basis set. Linear response TD-DFT calculations for excited states were performed within the Tamm–Dankoff approximation (TDA). Dielectric stabilisation was treated using a polarisable continuum model with \(\epsilon\) = 5.0, based on the experimentally-determined value26,69,70,71,72. To treat CT states, the exchange–correlation functional was range-tuned according to the established non-empirical procedure69,70,71,72, resulting in an optimal range-separation parameter of 0.010 Bohr−1. Electronic couplings and exciton/CT energetics were accessed by performing calculations in which the molecules were separated from one another by 10 Å to define localized/quasi-diabatic states. At this distance, interactions involving orbital overlap are negligible, resulting in the formation of localized CT states, and pure exciton states that can be localized by diabatization73. The energies and couplings between these states were calculated by projection onto the states of the molecular pair in the crystal packing geometry49,50.

Density of states calculations for electrons and holes were performed using the long-range embedding procedure available in VOTCA-CTP51,52,53,54,55,56,57. Molecular charge densities and polarizabilities for the ground state, cation, and anion were calculated using B3LYP/6–311 G(d,p). The charge densities were used to fit distributed atomic multipoles (specifically, atomic charges, dipoles, and quadrupoles) using GDMA (version 2.3)74. Polarizabilities for atoms in alkyl chains were taken directly from the set of Thole polarizabilities73. Atomic polarizabilities for atoms in the Y6 chromophore were obtained by scaling Thole polarizabilities to match the DFT-calculated polarisable volume of the core51. Solid state corrections to the gas-phase IE/EA were calculated using the long-range Ewald summation scheme in VOTCA-CTP. These calculations employed a molecular dynamics equilibrated 8 × 8 × 6 supercell constructed from the Y6 crystal structure24, with the π-π stacking direction aligned with the z-axis, to match the experimentally observed face-on packing preference relative to the substrate. This simulation cell was used to mimic a ~10 nm thin film by including a ~20 nm vacuum buffer in the z direction and applying the shape term for a periodically repeated supercell in the xy plane. The interaction cut-offs for electrostatic and polarization interactions were 8 and 6 nm, respectively.

Molecular dynamics simulations were conducted using the Y6 crystal structure with the missing atoms of the alkyl side chains built in. The system containing 8 × 8 × 6 unit cells was generated using Maestro75 in the Schrödinger Material Science Suite, and then prepared for molecular dynamics simulation with the OPLS3e force field using the System Builder in Desmond76,77,78. The prepared system first underwent 100 ps of Brownian minimization, followed by 1 ns of equilibration using NPT molecular dynamics simulation at 1 fs timestep, 300 K and 1.01325 bar. A further 10 ps molecular dynamics simulation was conducted to collect 200 frames of trajectories at 50 fs time interval for analysis.

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

The raw data and code for data processing can be found at https://github.com/Mikebprice/Free-charge-generation-in-Y6. Questions regarding this data set are encouraged to be sent to the corresponding authors.

Code that was used to support the findings of this study (and which has not been previously made available) can be accessed from https://github.com/PaulAlexanderHume/Free-charge-generation-in-Y6-theory. Questions regarding this data set are encouraged to be sent to the corresponding authors.

Few, S., Frost, J. M. & Nelson, J. Models of charge pair generation in organic solar cells. Phys. Chem. Chem. Phys. 17, 2311–2325 (2015).

Article CAS PubMed Google Scholar

Chen, S., Tsang, S. W., Lai, T. H., Reynolds, J. R. & So, F. Dielectric effect on the photovoltage loss in organic photovoltaic cells. Adv. Mater. 26, 6125–6131 (2014).

Article CAS PubMed Google Scholar

Halls, J. J. M. et al. Efficient photodiodes from interpenetrating polymer networks. Nature 376, 498–500 (1995).

Article ADS CAS Google Scholar

Tang, C. W. Two-layer organic photovoltaic cell. Appl. Phys. Lett. 48, 183 (1986).

Benduhn, J. et al. Intrinsic non-radiative voltage losses in fullerene-based organic solar cells. Nat. Energy 2, 17053 (2017).

Karuthedath, S. et al. Intrinsic efficiency limits in low-bandgap non-fullerene acceptor organic solar cells. Nat. Mater. 20, 378–384 (2020).

Article ADS PubMed Google Scholar

Xiao, F., Li, C., Sun, J. & Zhang, L. Knowledge domain and emerging trends in organic photovoltaic technology: A scientometric review based on citespace analysis. Front. Chem. 5, 1–12 (2017).

Article ADS Google Scholar

Kozlov, O. V. et al. Ultrafast exciton-to-polaron conversion in densely packed small organic semiconducting molecules. Adv. Opt. Mater. 5, 1–7 (2017).

Article ADS Google Scholar

Chandran, H. T. et al. Direct free carrier photogeneration in single layer and stacked organic photovoltaic devices. Adv. Mater. 29 (2017).

Nikolis, V. C. et al. Field Effect versus driving force: charge generation in small-molecule organic solar cells. Adv. Energy Mater. 10 (2020).

Wang, R. et al. Charge separation from an intra-moiety intermediate state in the high-performance PM6:Y6 organic photovoltaic blend. J. Am. Chem. Soc. 142, 12751–12759 (2020).

Article CAS PubMed Google Scholar

Reid, O. G., Pensack, R. D., Song, Y., Scholes, G. D. & Rumbles, G. Charge photogeneration in neat conjugated polymers. Chem. Mater. 26, 561–575 (2014).

Article CAS Google Scholar

Wang, R. et al. Ultrafast hole transfer mediated by polaron pairs in all-polymer photovoltaic blends. Nat. Commun. 10, 398 (2019).

Tautz, R. et al. Structural correlations in the generation of polaron pairs in low-bandgap polymers for photovoltaics. Nat. Commun. 3, 970 (2012).

Hou, J., Inganas, O., Friend, R. H. & Gao, F. Organic solar cells based on non-fullerene acceptors. Nat. Mater. 17, 119–128 (2018).

Article ADS CAS PubMed Google Scholar

Cui, Y. et al. Single-junction organic photovoltaic cells with approaching 18% efficiency. Adv. Mater. 32, 1–7 (2020).

ADS Google Scholar

Yuan, J. et al. Single-junction organic solar cell with over 15% efficiency using fused-ring acceptor with electron-deficient core. Joule 3, 1140–1151 (2019).

Article CAS Google Scholar

Lin, Y. et al. An electron acceptor challenging fullerenes for efficient polymer solar cells. Adv. Mater. 27, 1170–1174 (2015).

Article CAS PubMed Google Scholar

Yan, C. et al. Non-fullerene acceptors for organic solar cells. Nat. Rev. Mater. 3, 18003 (2018).

Wang, J. & Zhan, X. Fused-ring electron acceptors for photovoltaics and beyond. Acc. Chem. Res. 54, 132–143 (2021).

Cheng, P., Li, G., Zhan, X. & Yang, Y. Next-generation organic photovoltaics based on non-fullerene acceptors. Nat. Photonics 12, 131–142 (2018).

Perdigón-Toro, L. et al. Barrierless Free Charge Generation in the High-Performance PM6:Y6 Bulk Heterojunction Non-Fullerene Solar Cell. Adv. Mater. 32 (2020).

Zhu, X. et al. Self-stimulated dissociation in non-fullerene organic bulk-heterojunction solar cells self-stimulated dissociation in non-fullerene organic bulk-heterojunction solar cells. Joule 4, 1–15 (2020).

Article Google Scholar

Zhang, G. et al. Delocalization of exciton and electron wavefunction in non-fullerene acceptor molecules enables efficient organic solar cells. Nat. Commun. 11, 1–10 (2020).

Article ADS Google Scholar

Liu, J. et al. Charge separation boosts exciton diffusion in fused ring electron acceptors. J. Mater. Chem. A 8 (2020).

Kerremans, R. et al. The optical constants of solution-processed semiconductors—new challenges with perovskites and non-fullerene acceptors. Adv. Opt. Mater. 2000319 (2020).

Park, S. Y. et al. Photophysical pathways in efficient bilayer organic solar cells: The importance of interlayer energy transfer. Nano Energy 84, 105924 (2021).

Wen, Z. C., Yin, H. & Hao, X. T. Recent progress of PM6:Y6-based high efficiency organic solar cells. Surf. Interfaces 23, 100921 (2021).

Article CAS Google Scholar

Lin, Y. et al. 17.1% Efficient single-junction organic solar cells enabled by n-type doping of the bulk-heterojunction. Adv. Sci. 1903419, 1–9 (2020).

Google Scholar

Wang, Y. et al. High-efficiency perovskite quantum dot hybrid nonfullerene organic solar cells with near-zero driving force. Adv. Mater. 32, (2020).

Lin, Q. et al. Electro-optics of perovskite solar cells. Nat. Photonics 9, 106–112 (2014).

Article ADS Google Scholar

Li, H. et al. On the study of exciton binding energy with direct charge generation in photovoltaic polymers. Adv. Electron. Mater. 2 (2016).

Mondelli, P. et al. Meta-analysis: The molecular organization of non-fullerene acceptors. Mater. Horiz. 7, 1062–1072 (2020).

Article CAS Google Scholar

Northrup, J. E. Mobility enhancement in polymer organic semiconductors arising from increased interconnectivity at the level of polymer segments. Appl. Phys. Lett. 106 (2015).

Gélinas, S. et al. Ultrafast long-range charge separation in organic semiconductor photovoltaic diodes. Sci. (80-.) 343, 512–517 (2014).

Article ADS Google Scholar

Westenhoff, S. et al. Charge recombination in organic photovoltaic devices with high open-circuit voltages. J. Am. Chem. Soc. 130, 13653–13658 (2008).

Article CAS PubMed Google Scholar

Bakulin, A. A. et al. The role of driving energy and delocalized states for charge separation in organic semiconductors. Science 335, 1340–1344 (2012).

Article ADS CAS PubMed Google Scholar

Gélinas, S. et al. Ultrafast long-range charge separation in organic semiconductor photovoltaic diodes. Science 343, 512–516 (2014).

Article ADS PubMed Google Scholar

Gillett, A. J. et al. The role of charge recombination to spin-triplet excitons in non-fullerene acceptor organic solar cells. Nature 597, 666–671 (2020).

Article ADS Google Scholar

Yao, N. et al. Efficient charge transport enables high efficiency in dilute donor organic solar cells. J. Phys. Chem. Lett. 12, 5039–5044 (2021).

Article CAS PubMed PubMed Central Google Scholar

Richter, J. M. et al. Enhancing photoluminescence yields in lead halide perovskites by photon recycling and light out-coupling. Nat. Commun. 7, 13941 (2016).

Article ADS CAS PubMed PubMed Central Google Scholar

Yablonovitch, E., Miller, O.D. & Kurtz, S.R. A great solar cell also needs to be a great LED: External fluorescence leads to new efficiency record. in AIP Conference Proceedings vol. 1519 (2013).

Kyaw, A. K. K. et al. Intensity dependence of current-voltage characteristics and recombination in high-efficiency solution-processed small-molecule solar cells. ACS Nano 7, 4569–4577 (2013).

Koster, L. J. A., Kemerink, M., Wienk, M. M., Maturová, K. & Janssen, R. A. J. Quantifying bimolecular recombination losses in organic bulk heterojunction solar cells. Adv. Mater. 23, 1670–1674 (2011).

Koster, L. J. A., Mihailetchi, V. D., Xie, H. & Blom, P. W. M. Origin of the light intensity dependence of the short-circuit current of polymer/fullerene solar cells. Appl. Phys. Lett. 87, 203502 (2005).

Mihailetchi, V. D., Wildeman, J. & Blom, P. W. M. Space-charge limited photocurrent. Phys. Rev. Lett. 94, 126602 (2005).

Firdaus, Y. et al. Long-range exciton diffusion in molecular non-fullerene acceptors. Nat. Commun. 11, 5220 (2020).

Liao, S. H., Jhuo, H. J., Cheng, Y. S. & Chen, S. A. Fullerene derivative-doped zinc oxide nanofilm as the cathode of inverted polymer solar cells with low-bandgap polymer (PTB7-Th) for high performance. Adv. Mater. 25, 4766–4771 (2013).

Hume, P. A. & Hodgkiss, J. M. A projective method for the calculation of excited-state electronic coupling: isolating charge transfer/recombination processes in organic photovoltaics. J. Phys. Chem. A 124, 591–600 (2020).

Article CAS PubMed Google Scholar

Hume, P. A., Jiao, W. & Hodgkiss, J. M. Long-range exciton diffusion in a non-fullerene acceptor: approaching the incoherent limit. J. Mater. Chem. C. 9, 1419–1428 (2021).

Article CAS Google Scholar

Poelking, C. & Andrienko, D. Long-range embedding of molecular ions and excitations in a polarizable molecular environment. J. Chem. Theory Comput. 12, 4516–4523 (2016).

Article CAS PubMed Google Scholar

Kirkpatrick, J. An approximate method for calculating transfer integrals based on the ZINDO Hamiltonian. Int. J. Quantum Chem. 108, 51–56 (2008).

Rühle, V., Junghans, C., Lukyanov, A., Kremer, K. & Andrienko, D. Versatile object-oriented toolkit for coarse-graining applications. J. Chem. Theory Comput. 5, 3211–3223 (2009).

Baumeier, B., Kirkpatrick, J. & Andrienko, D. Density-functional based determination of intermolecular charge transfer properties for large-scale morphologies. Phys. Chem. Chem. Phys. 12 (2010).

Lukyanov, A. & Andrienko, D. Extracting nondispersive charge carrier mobilities of organic semiconductors from simulations of small systems. Phys. Rev. B - Condens. Matter Mater. Phys. 82, 193202 (2010).

Rühle, V. et al. Microscopic simulations of charge transport in disordered organic semiconductors. J. Chem. Theor. Comput. 7, 3335–3345 (2011).

Kordt, P. & Andrienko, D. Modeling of spatially correlated energetic disorder in organic semiconductors. J. Chem. Theory Comput. 12, 36–40 (2016).

Dong, Y. et al. Orientation dependent molecular electrostatics drives efficient charge generation in homojunction organic solar cells. Nat. Commun. 11, 4617 (2020).

Poelking, C. & Andrienko, D. Design rules for organic donor-acceptor heterojunctions: Pathway for charge splitting and detrapping. J. Am. Chem. Soc. 137, 6320–6326 (2015).

Article CAS PubMed Google Scholar

Schwarze, M. et al. Impact of molecular quadrupole moments on the energy levels at organic heterojunctions. Nat. Commun. 10, 1–9 (2019).

Article CAS Google Scholar

Poelking, C. et al. Impact of mesoscale order on open-circuit voltage in organic solar cells. Nat. Mater. 14, 434–439 (2015).

Article ADS CAS PubMed Google Scholar

Piris, J. et al. Photogeneration and ultrafast dynamics of excitons and charges in P3HT/PCBM blends. J. Phys. Chem. C. 113, 14500–14506 (2009).

Article CAS Google Scholar

Zou, Y. & Holmes, R. J. The role of exciton ionization processes in bulk heterojunction organic photovoltaic cells. Adv. Energy Mater. 5, 1500019 (2015).

Paquin, F. et al. Charge separation in semicrystalline polymeric semiconductors by photoexcitation: Is the mechanism intrinsic or extrinsic? Phys. Rev. Lett. 106, 1–4 (2011).

Article Google Scholar

de Mello, J. C., Wittmann, H. F. & Friend, R. H. An improved experimental determination of external photoluminescence quantum efficiency. Adv. Mater. 9, 230–232 (1997).

Article Google Scholar

Tiwana, P., Parkinson, P., Johnston, M. B., Snaith, H. J. & Herz, L. M. Ultrafast terahertz conductivity dynamics in mesoporous TiO2: Influence of dye sensitization and surface treatment in solid-state dye-sensitized solar cells. J. Phys. Chem. C. 114, 1365–1371 (2010).

Article CAS Google Scholar

Frisch, M. J. et al. Gaussian 09 Revision D.01 (2010).

Yanai, T., Tew, D. P. & Handy, N. C. A new hybrid exchange-correlation functional using the Coulomb-attenuating method (CAM-B3LYP). Chem. Phys. Lett. 393, 51–57 (2004).

Kronik, L., Stein, T., Refaely-Abramson, S. & Baer, R. Excitation gaps of finite-sized systems from optimally tuned range-separated hybrid functionals. J. Chem. Theory Comput. 8, 1515–1531 (2012).

Article CAS PubMed Google Scholar

Sun, H. et al. Ionization energies, electron affinities, and polarization energies of organic molecular crystals: quantitative estimations from a Polarizable Continuum Model (PCM)-Tuned Range-Separated Density Functional Approach. J. Chem. Theory Comput. 12, 2906–2916 (2016).

Refaely-Abramson, S. et al. Gap renormalization of molecular crystals from density-functional theory. Phys. Rev. B - Condens. Matter Mater. Phys. 88, 081204 (2013).

Chen, X. K. & Brédas, J. L. Voltage losses in organic solar cells: understanding the contributions of intramolecular vibrations to nonradiative recombinations. Adv. Energy Mater. 8, 1702227 (2018).

Hsu, C. P., You, Z. Q. & Chen, H. C. Characterization of the short-range couplings in excitation energy transfer. J. Phys. Chem. C. 112, 1204–1212 (2008).

Article CAS Google Scholar

Stone, A. J. Distributed multipole analysis: Stability for large basis sets. J. Chem. Theory Comput. 1, 1128–1132 (2005).

Van Duijnen, P. T. & Swart, M. Molecular and atomic polarizabilities: Thole’s model revisited. J. Phys. Chem. A 102, 2399–2407 (1998).

Schrödinger Release 2020-2: Maestro. (Schrödinger, LLC, New York, 2020).

Shaw, D. E. et al. Millisecond-scale molecular dynamics simulations on Anton. Proc. Conf. High. Perform. Comput. Netw., Storage Anal., SC ’09 65, 1–11 (2009).

Google Scholar

Schrödinger Release 2020-2: Desmond Molecular Dynamics System. (Schrödinger, LLC, New York, 2020).

Download references

M.B.P., P.A.H., A.I., I.W., R.R.T., K.E.T., N.J.L.K.D., K.C. and J.M.H. acknowledge support from the MacDiarmid Institute of New Zealand. P.A.H., A.I., I.W., R.R.T., K.E.T., N.J.L.K.D., K.C. and J.M.H. acknowledge funding from the Marsden Fund of NZ. M.B.P. and N.J.L.K.D. acknowledge support from the Royal Society of New Zealand. A.G, P.J.C. and G.L. acknowledge support from ARC Centre of Excellence in Exciton Science (funding grant number CE170100026). P.A.H. would like to acknowledge the support of Professor Denis Andrienko, Dr Anastasia Markina and Dr. Carl Pölking for support in utilising VOTCA software. All calculations were performed using the VUW High Performance ComputingCluster ‘‘Rāpoi’’.

These authors contributed equally: Michael B. Price, Paul A. Hume.

School of Chemical and Physical Sciences, Victoria University of Wellington, Wellington, New Zealand

Michael B. Price, Paul A. Hume, Aleksandra Ilina, Isabella Wagner, Ronnie R. Tamming, Karen E. Thorn, Nathaniel J. L. K. Davis, Kai Chen & Justin M. Hodgkiss

MacDiarmid Institute for Advanced Materials and Nanotechnology, Wellington, New Zealand

Michael B. Price, Paul A. Hume, Aleksandra Ilina, Isabella Wagner, Ronnie R. Tamming, Karen E. Thorn, Nathaniel J. L. K. Davis, Kai Chen & Justin M. Hodgkiss

Wellington UniVentures, Victoria University of Wellington, Wellington, New Zealand

Ronnie R. Tamming & Kai Chen

Robinson Research Institute, Faculty of Engineering, Victoria University of Wellington, Wellington, New Zealand

Ronnie R. Tamming & Kai Chen

The Dodd-Walls Centre for Photonic and Quantum Technologies, Dunedin, New Zealand

Ronnie R. Tamming & Kai Chen

Ferrier Research Institute, Victoria University of Wellington, Wellington, New Zealand

Wanting Jiao

ARC Centre of Excellence in Exciton Science, School of Chemistry, University of Sydney, Sydney, NSW, Australia

Alison Goldingay, Patrick J. Conaghan & Girish Lakhwani

School of Materials Science and Engineering, Peking University, Beijing, China

Yifan Wang, Peiyao Xue, Heng Lu & Xiaowei Zhan

College of Materials Science and Engineering, Qingdao University, Qingdao, China

Yifan Wang

You can also search for this author in PubMed Google Scholar

M.B.P., P.A.H.: conceptualisation, investigation, experiments, supervision, writing. M.B.P., P.A.H. contributed equally. A.I., I.W., R.R.T., K.E.T., A.G., P.J.C., N.J.L.K.D., K.C., P.X., H.L.: experiments, investigation. W.J.: molecular dynamics simulations. G.L.: supervision. Y.W.: supervision, investigation, conceptualisation, X.Z.: supervision, conceptualisation, funding acquisition, J.M.H.: conceptualisation, investigation, supervision, writing, funding acquisition.

Correspondence to Michael B. Price, Paul A. Hume or Justin M. Hodgkiss.

The authors declare no competing interests.

Nature Communications thanks Chunfeng Zhang, Yufei Zhong and the other anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

This manuscript was previously reviewed at another journal within the Nature Portfolio. This document only contains reviewer comments and rebuttal letters for versions considered at Nature Communications.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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/.

Reprints and Permissions

Price, M.B., Hume, P.A., Ilina, A. et al. Free charge photogeneration in a single component high photovoltaic efficiency organic semiconductor. Nat Commun 13, 2827 (2022). https://doi.org/10.1038/s41467-022-30127-8

Download citation

Received: 29 November 2021

Accepted: 15 April 2022

Published: 20 May 2022

DOI: https://doi.org/10.1038/s41467-022-30127-8

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Nature Communications (2023)

Nature Energy (2022)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

News

Free charge photogeneration in a single component high photovoltaic efficiency organic semiconductor