Method of improving exciton dissociation at organic donor-acceptor heterojunctions

Information

  • Patent Grant
  • 11302881
  • Patent Number
    11,302,881
  • Date Filed
    Tuesday, May 3, 2011
    13 years ago
  • Date Issued
    Tuesday, April 12, 2022
    2 years ago
Abstract
The present disclosure generally relates to organic photosensitive optoelectronic devices and polaron pair recombination dynamics to impact efficiency and open circuit voltages of organic solar cells. The present disclosure also relates, in part, to methods of making organic photosensitive optoelectronic devices comprising the same.
Description
FIELD OF THE DISCLOSURE

The present disclosure generally relates to organic photosensitive optoelectronic devices and polaron pair recombination dynamics to impact efficiency and open circuit voltages of organic solar cells. The present disclosure also relates, in part, to methods of making organic photosensitive optoelectronic devices comprising the same.


BACKGROUND

Optoelectronic devices rely on the optical and electronic properties of materials to either produce or detect electromagnetic radiation electronically or to generate electricity from ambient electromagnetic radiation.


Photosensitive optoelectronic devices convert electromagnetic radiation into electricity. Solar cells, also called photovoltaic (PV) devices, are a type of photosensitive optoelectronic device that is specifically used to generate electrical power. PV devices, which may generate electrical energy from light sources other than sunlight, can be used to drive power consuming loads to provide, for example, lighting, heating, or to power electronic circuitry or devices such as calculators, radios, computers or remote monitoring or communications equipment. These power generation applications also often involve the charging of batteries or other energy storage devices so that operation may continue when direct illumination from the sun or other light sources is not available, or to balance the power output of the PV device with a specific application's requirements. As used herein the term “resistive load” refers to any power consuming or storing circuit, device, equipment or system.


Another type of photosensitive optoelectronic device is a photoconductor cell. In this function, signal detection circuitry monitors the resistance of the device to detect changes due to the absorption of light.


Another type of photosensitive optoelectronic device is a photodetector. In operation a photodetector is used in conjunction with a current detecting circuit which measures the current generated when the photodetector is exposed to electromagnetic radiation and may have an applied bias voltage. A detecting circuit as described herein is capable of providing a bias voltage to a photodetector and measuring the electronic response of the photodetector to electromagnetic radiation.


These three classes of photosensitive optoelectronic devices may be characterized according to whether a rectifying junction as defined below is present and also according to whether the device is operated with an external applied voltage, also known as a bias or bias voltage. A photoconductor cell does not have a rectifying junction and is normally operated with a bias. A PV device has at least one rectifying junction and is operated with no bias. A photodetector has at least one rectifying junction and is usually but not always operated with a bias. As a general rule, a photovoltaic cell provides power to a circuit, device or equipment, but does not provide a signal or current to control detection circuitry, or the output of information from the detection circuitry. In contrast, a photodetector or photoconductor provides a signal or current to control detection circuitry, or the output of information from the detection circuitry but does not provide power to the circuitry, device or equipment.


Traditionally, photosensitive optoelectronic devices have been constructed of a number of inorganic semiconductors, e.g., crystalline, polycrystalline and amorphous silicon, gallium arsenide, cadmium telluride and others. Herein the term “semiconductor” denotes materials which can conduct electricity when charge carriers are induced by thermal or electromagnetic excitation. The term “photoconductive” generally relates to the process in which electromagnetic radiant energy is absorbed and thereby converted to excitation energy of electric charge carriers so that the carriers can conduct, i.e., transport, electric charge in a material. The terms “photoconductor” and “photoconductive material” are used herein to refer to semiconductor materials which are chosen for their property of absorbing electromagnetic radiation to generate electric charge carriers.


PV devices may be characterized by the efficiency with which they can convert incident solar power to useful electric power. Devices utilizing crystalline or amorphous silicon dominate commercial applications, and some have achieved efficiencies of 23% or greater. However, efficient crystalline-based devices, especially of large surface area, are difficult and expensive to produce due to the problems inherent in producing large crystals without significant efficiency-degrading defects. On the other hand, high efficiency amorphous silicon devices still suffer from problems with stability. Present commercially available amorphous silicon cells have stabilized efficiencies between 4 and 8%.


PV devices may be optimized for maximum electrical power generation under standard illumination conditions (i.e., Standard Test Conditions which are 1000 W/m2, AM1.5 spectral illumination), for the maximum product of photocurrent times photovoltage. The power conversion efficiency of such a cell under standard illumination conditions depends on the following three parameters: (1) the current under zero bias, i.e., the short-circuit current ISC, in Amperes (2) the photovoltage under open circuit conditions, i.e., the open circuit voltage VOC, in Volts and (3) the fill factor, ff.


PV devices produce a photo-generated current when they are connected across a load and are irradiated by light. When irradiated under infinite load, a PV device generates its maximum possible voltage, V open-circuit, or VOC. When irradiated with its electrical contacts shorted, a PV device generates its maximum possible current, I short-circuit, or ISC. When actually used to generate power, a PV device is connected to a finite resistive load and the power output is given by the product of the current and voltage, I×V. The maximum total power generated by a PV device is inherently incapable of exceeding the product, ISC×VOC. When the load value is optimized for maximum power extraction, the current and voltage have the values, Imax and Vmax, respectively.


A figure of merit for PV devices is the fill factor, ff, defined as:

ff={ImaxVmax}/{ISCVOC}  (1)


where ff is always less than 1, as ISC and VOC are never obtained simultaneously in actual use. Nonetheless, as ff approaches 1, the device has less series or internal resistance and thus delivers a greater percentage of the product of ISC and VOC to the load under optimal conditions. Where Pinc is the power incident on a device, the power efficiency of the device, ηP, may be calculated by:

ηP=ff*(ISC*VOC)Pinc


To produce internally generated electric fields that occupy a substantial volume of the semiconductor, the usual method is to juxtapose two layers of material with appropriately selected conductive properties, especially with respect to their distribution of molecular quantum energy states. The interface of these two materials is called a photovoltaic junction. In traditional semiconductor theory, materials for forming PV junctions have been denoted as generally being of either n or p type. Here n-type denotes that the majority carrier type is the electron. This could be viewed as the material having many electrons in relatively free energy states. The p-type denotes that the majority carrier type is the hole. Such material has many holes in relatively free energy states. The type of the background, i.e., not photo-generated, majority carrier concentration depends primarily on unintentional doping by defects or impurities. The type and concentration of impurities determine the value of the Fermi energy, or level, within the gap between the conduction band minimum and valance band maximum energies. The Fermi energy characterizes the statistical occupation of molecular quantum energy states denoted by the value of energy for which the probability of occupation is equal to ½. A Fermi energy near the conduction band minimum energy indicates that electrons are the predominant carrier. A Fermi energy near the valence band maximum energy indicates that holes are the predominant carrier. Accordingly, the Fermi energy is a primary characterizing property of traditional semiconductors and the prototypical PV junction has traditionally been the p-n interface.


The term “rectifying” denotes, inter alia, that an interface has an asymmetric conduction characteristic, i.e., the interface supports electronic charge transport preferably in one direction. Rectification is associated normally with a built-in electric field which occurs at the junction between appropriately selected materials.


The current-voltage characteristics of organic heterojunctions are often modeled using the generalized Shockley equation derived for inorganic diodes. However, since the Shockley equation does not rigorously apply to organic semiconductor donor-acceptor (D-A) heterojunctions (HJs), the extracted parameters lack a clear physical meaning.


Without wishing to be bound by theory, a mathematical derivation for the current density-voltage (J-V) characteristics directed to D-A HJ solar cells is provided. This derivation predicts the general dependence of dark current, open-circuit voltage (VOC), and short-circuit current (JSC) on temperature and light intensity as well as the maximum VOC for a given D-A material pair. Also provided herein is the investigation of kPPr at CuPc/C60 and SubPc/C60 HJs into PP dynamics in archetype planar organic solar cells comprising, for example, copper phthalocyanine and boron subphthalocyanine chloride donor-acceptor pairs.


Due to the disorder at the D-A interface, trap-limited recombination leads to the introduction of two temperature-dependent ideality factors. Such factors describe the dark current of copper phthalocyanine/C60 and boron subphthalocyanine/C60 cells at low temperature, and fits to the generalized Shockley equation break down. As described in detail below, one embodiment of the present disclosure identifies the polaron pair recombination rate as a key factor that determines the J-V characteristics in the dark and under illumination and provides direct measurements of this process. These results provide a general physical framework for interpreting the J-V characteristics and understanding the efficiency of both small molecule and polymer organic, planar and bulk HJ solar cells.


SUMMARY

There is disclosed an organic photosensitive optoelectronic device comprising:


two electrodes comprising in superposed relation;


at least one donor material, and


at least one acceptor material,


the at least one donor material and at least one acceptor material forming at least one photoactive region between the two electrodes, wherein the at least one donor material and at least one acceptor material are ordered in their bulk and comprise a disordered region at their interface.


There is also disclosed a method of making an organic photosensitive optoelectronic device, the method comprising depositing onto a substrate,


at least one donor material,


at least one acceptor material, and


at least one electrode,


wherein the at least one donor material and at least one acceptor material are deposited to form at least one photoactive region between the two electrodes,


wherein the at least one donor material and at least one acceptor material are ordered in their bulk and comprise a disordered region at their interface.


The foregoing and other features of the present disclosure will be more readily apparent from the following detailed description of exemplary embodiments, taken in conjunction with the attached drawings. It will be noted that for convenience all illustrations of devices show the height dimension exaggerated in relation to the width.





BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1. (a) represents an energy level diagram showing the anode and cathode work functions, WFa and WFc, and their associated injection barriers ϕa and ϕc, respectively. The interfacial gap, ΔEHL, is the energy difference between the highest occupied molecular orbital energy of the donor and the lowest unoccupied molecular orbital energy of the acceptor. Current is unipolar in the donor (JP) and acceptor (Jn) layers, and is determined from generation/recombination in the heterojunction (HJ) region, roughly defined by the spatial extent, a0, of the polaron pair distribution at the interface. (b) Represents the processes occurring within the HJ region. Excitons diffuse, with current density, JX, to the HJ and undergo charge transfer to form polaron pairs. These may recombine, at rate kPPr, or dissociate with rate, kPPd, as determined by the Onsager-Braun model (C. L. Braun, J. Chem. Phys. 80, 4157 (1984)). The current density, J, contributes to the interfacial free electron (nI) and hole (pI) densities, which bimolecularly recombine to form polaron pairs at rate krec.



FIG. 2. (a) Shows the calculated dark current density (J)-voltage characteristics over the temperature range from 120K to 300 K in 20 K increments using Eq. (13) and the parameters in Table 2. The slope in reverse bias is due to the field dependent dissociation of thermally generated polaron pairs. In forward bias, three regions are apparent. At Va<0.3 V, trap-limited recombination involving free acceptor electrons and trapped donor holes dominates, and the slope is proportional to the donor ideality factor, nD. At higher bias, the inverse process dominates (i.e. free donor holes recombine with trapped acceptor electrons) and the slope is proportional to the acceptor ideality factor, nA. Series resistance (Rs) limits the current at Va>0.8 V. (b) Represents the diode parameters n and n ln(Js) corresponding to the dark currents in (a). Both ideality factors increase with decreasing temperature.



FIG. 3. (a) Shows the intensity dependence of the open-circuit voltage (Voc), short-circuit current (Jsc) and fill factor (FF) derived from the same model and parameters as in FIG. 2. (b). Temperature dependence of Voc, FF, and polaron pair dissociation efficiency at short circuit, ηPPd (Jsc).



FIG. 4. Shows the dark current density vs. forward voltage for (a) CuPc/C60 and (b) SubPc/C60 devices recorded for T=296, 275, 247, 218, 193, 171, 155, 145, 128, and 114 K. Bold (red) lines indicate fits to Eq. (18) in the description. Thin (black) lines connect the data points and serve as a guide to the eye. Both data sets are refit using the generalized Shockley Equation in (c) and (d), where the difference between data and theory is most pronounced at low voltage and temperature. Several intermediate temperatures have been omitted in (c) and (d) for clarity.



FIG. 5. (a) Represents temperature dependence of the ideality factors derived from the J-V data of FIG. 4. The ideality factors increase with decreasing temperature in all cases. Values of nD and JsD) cannot be extracted for CuPc/C60 at T>220 K because the second exponential term in Eq. (18) is insignificant. (b) Shows a plot of n ln(Js) vs 1/kbT, where kbT is the thermal energy, for the CuPc/C60 and SubPc/C60 devices. Note the similarity in nA and nA 1n(JsA) between the two cells, which share a common C60 acceptor, whereas nD and nD 1n(JsD) clearly differ due to the different donors.



FIG. 6. (a) Represents the dependence of the open circuit voltage, Voc, on the intensity of simulated AM1.5 G solar illumination at room temperature. Solid lines are calculated from Eq. (18) using values of n and Js from FIG. 5 along with the short-circuit current densities measured at each intensity. (b) Shows the temperature dependence of Voc under laser illumination at λ=496 nm, and at an intensity of 30 mW/cm2. Solid lines are calculated using the temperature dependent n and Js of FIG. 5 and the short-circuit current density measured at room temperature. Deviation between the calculation and the data at low temperature indicates a reduced polaron pair dissociation efficiency, ηPPd(T), shown on the right hand axis, and normalized to its value at room temperature. Dashed lines connecting the data points are a guide to the eye.



FIG. 7. (a) Represents the dependence of open circuit voltage, fill factor and short circuit current (Voc, FF, Jsc, respectively), and normalized power conversion efficiency, ηP, on the polaron pair recombination rate, kPPr. An exciton current of qJX=10 mA/cm2 is assumed. (b) Shows the Fourth quadrant J-V characteristics at several values of kPPr resulting in the trends shown in (a).



FIG. 8. (a) Shows the Current-Voltage characteristics calculated for Vbi ranging from 0.1 to 0.7 V in 0.2 V increments using the same model as in FIG. 7b with kPPr=(20 ns)−1. (b) Is a schematic illustrating the physical basis for the ‘S-kink’ observed in (a). The top panel shows the dependence of the internal field, FI, and polaron pair dissociation efficiency, ηPPd, on the applied bias, Va. Dissociation is aided by FI<0. The lower panel shows the corresponding photocurrents, which sum with the dark current to give the total current (dashed lines). When Vbi is small, the internal field reaches zero at small Va, resulting in a rapid change in photocurrent before the dark current becomes appreciable, leading to the ‘S-kink’ as shown.



FIG. 9. (a) Shows the emission spectra of two solar cells as an embodiment of the present disclosure. The red solid line shows SubPc/C60 electroluminescence obtained at a constant drive current of 450 mA/cm2. Black closed and blue open circles denote SubPc/C60 and CuPc/C60 photoluminescence, respectively, excited at λex=442 nm with a modulated intensity of 2 mW/cm2 rms. (b) Shows the absolute SubPc/C60 photoluminescence intensity measured for the following structure: BCP (10 nm)/C60 (X nm)/SubPc (11 nm) deposited on quartz, where X=10, 20, 30, and 40 nm. Inset: Open circles and squares show the intensity for SubPc fluorescence and exciplex/fullerene emission integrated over the wavelength domains λem={575, 670} nm and λem={730, 800} nm, respectively. Lines indicate the predictions of a combined optical interference and diffusion model, where the exciplex emission is taken to be proportional to the exciton diffusive flux reaching the SubPc/C60 interface. The fluorescence intensity of SubPc decreases in spite of its constant thickness due to an interference minimum that develops at its location as the C60 thickness increases.



FIG. 10. (a) Represents the photoluminescence intensity as a function of electrical bias for CuPc/C60 and SubPc/C60 HJ devices obtained for a steady-state optical bias intensity of 10 mW/cm2. The λem=730 nm trends follow the intensity of exciplex emission from the heterojunction of each cell, whereas the λem=615 nm trend follows the SubPc fluorescence intensity. Solid lines are guides for the eye. (b) Shows SubPc/C60 exciplex emission for optical biases of 0.8 and 30 mW/cm2, where emission from the 0.8 W/cm2 bias has lower intensity because the modulated signal excitation beam was decreased from 2 to 0.2 mW/cm2 to prevent it from dominating the optical bias. The current-voltage characteristics for each device, obtained concurrently with the luminescence data, are shown by solid lines, with the scale on the right hand ordinate. Red lines through the intensity data denote linear fits to estimate the emission peak positions, which scale with open-circuit voltage.



FIG. 11. (a) Is a Cole-Cole representation of intensity modulated photocurrent spectra (IMPS) for the SubPc/C60 cell. Near short-circuit, the complex quantum efficiency, {tilde over (Φ)} accumulates a small phase lag with increasing modulation frequency due to the charge carrier transit time from the interface to the contacts. At open-circuit, {tilde over (Φ)} exhibits a phase advance due to superposition of an oppositely directed bimolecular recombination current at the heterojunction. (b) Similar IMPS spectra obtained for CuPc/C60. The contraction |{tilde over (Φ)}|→0 and lack of phase advance at open-circuit indicates that free carrier generation is suppressed due to geminate recombination.



FIG. 12. (a) Shows the simulated IMPS spectra assuming constant kPPr (100 ns)−1 (black solid lines) and 1 ns−1 (blue dashed lines) using Eq. (23) in the text. (b) Shows the simulated IMPS spectra assuming the field-dependent form for kPPr as given by Eq. (24) in text. In this case, the zero-field rates are kPPr0=(100 ns)−1 (black solid lines) and 1 ns−1 (blue dashed lines). The inset shows the predicted polaron pair density in each case (solid lines), that agree with the SubPc/C60 (black circles) and CuPc/C60 (blue squares) exciplex data of FIG. 10a, up to the open circuit voltage.



FIG. 13. Represents a proposed field-dependent model of the recombination rate, kPPr. The interfacial electric field modifies the polaron pair potential, and hence the activation energy required for polaron pair recombination. In reverse bias, FI<0, reducing the PP potential to favor dissociation and increasing the barrier for recombination. As FI increases, the dissociation rate decreases and the recombination rate increases. The polaron pair Coulomb potential is asymmetric due to the offset between donor and acceptor lowest unoccupied molecular orbitals, whose energy scheme is shown at the top of the figure.





The current density vs. voltage (J-V) characteristics of organic semiconductor heterojunctions are often similar to those of inorganic p-n junctions. As a consequence, theoretical treatments based on the generalized Shockley Equation originally derived for inorganic devices, have been extended to model the operation of organic solar cells. See, e.g., W. Shockley, Bell Syst. Tech. J. 28, 435 (1949); W. Shockley, Electrons and Holes in Semiconductors (D. Van Nostrand, Inc., New York, 1950); V. Dyakonov, Appl. Phys. A 79, 21 (2004); W. J. Potscavage, S. Yoo, and B. Kippelen, Appl. Phys. Lett. 93, 193308 (2008); B. P. Rand, D. P. Burk, and S. R. Forrest, Phys. Rev. B 75, 115327 (2007); and N. Li, B. E. Lassiter, R. R. Lunt, G. Wei, and S. R. Forrest, Appl. Phys. Lett. 94, 023307 (2009).


This phenomenological approach often yields a reasonably accurate description, although it fails in some cases and obscures the inherently different physics of organic semiconductors whose understanding would lead to a richer, more rigorous picture of these structures.


The Shockley Equation is derived for inorganic p-n junctions with well-developed energy band structure where thermal and optical excitation results in delocalized free charge carriers. In contrast, organic semiconductors are generally characterized by hopping transport and tightly bound, localized exciton states that require significant energy to dissociate into free charge carriers. Thus, in one aspect, the need exists to develop an ideal diode equation specifically for the case of organic heterojunctions (HJs).


By explicitly treating polaron pair generation, recombination and dissociation at the HJ, provided herein is a current-voltage characteristic similar in form to the Shockley equation as described in S. M. Sze, Physics of Semiconductor Devices (John Wiley & Sons, New York, 2005), but differing in several key aspects. First, it predicts the temperature dependence of the dark current, the intensity and temperature dependence of the open circuit voltage (Voc) and short-circuit current (Jsc), and the maximum Voc attainable for a given small molecule or polymer HJ material pair. Additionally, and without wishing to be bound by theory, the diode ideality factor arose from recombination via disorder-induced traps at the heterointerface, and analytically treats the case in which these traps have an exponential energetic distribution extending into the energy gap between the highest occupied (HOMO) and lowest unoccupied molecular orbitals (LUMO). This model, applies to, for example, data obtained for two archetype, planar HJ organic photovoltaic cells composed of copper phthalocyanine (CuPc) and boron subphthalocyanine chloride (SubPc) donors, and a fullerene (such as, for example, C60) acceptor.


The treatment described herein assumes that current is governed solely by generation and recombination at the heterojunction, and that both processes proceed through the polaron pair (PP) intermediate state. The average polaron pair separation, a0, thus defines the ‘volume’ of the heterojunction region, as shown in FIG. 1a. Current outside of this region is unipolar, with pure hole and electron currents flowing in the donor and acceptor bulk, respectively.


The schematic in FIG. 1a also defines the energetics, where the interfacial gap, ΔEHL, is the difference between donor HOMO and acceptor LUMO along with any shift due to formation of an interface dipole. The hole and electron injection barriers at the anode and cathode are ϕa and ϕc, respectively, again including any interface dipoles, and the built-in potential of the device is given by the corresponding difference in contact work functions: Vbi=WFa−WFc.



FIG. 1b shows the processes that occur within the HJ volume. The recombination of polaron pairs is described via:













J
X


a
0


-


k
PPr



(

ζ
-

ζ
eq


)


-


k
PPd


ζ

+


k
rec



n
I



p
I



=
0

,




(
1
)








and for free carriers:













k
PPd


ζ

-


k
rec



n
I



p
I


+

J

qa
0



=
0

,




(
2
)








where steady-state conditions are assumed. Here, ζ is the PP density, JX is the exciton current density diffusing to the interface, J is the charge current density flowing through the device, q is the electron charge, and nI and pI are the interfacial free electron and hole densities, respectively. Definitions of select variables used in this section are summarized in Table 1.


Polaron pairs recombine to the ground state at rate kPPr, which is also linked to the thermal equilibrium PP population, ζeq, determined by detailed balance. See, e.g., T. Kirchartz, B. E. Pieters, K. Taretto, and U. Rau, J. Appl. Phys. 104, 094513 (2008). Polaron pairs dissociate at rate kPPd, which is a function of temperature and the electric field at the interface according to the Onsager-Braun model (C. L. Braun, J. Chem. Phys. 80, 4157 (1984)), see Calculation of the Interfacial Electric Field below.


Finally, free carriers bimolecularly recombine to form PPs with rate constant, krec, approximated by its bulk Langevin value, qμtot/ε, as described in M. Pope and C. E. Swenberg, Electronic Processes in Organic Crystals and Polymers (Oxford University Press, New York, 1999) and C. Groves and N. C. Greenham, Phys. Rev. B 78, 155205 (2008).


Here, μtot is the sum of the electron and hole mobilities in the acceptor and donor layers, respectively, and ε is the average permittivity.


Solving Eq. (1) for the PP density and substituting the result into Eq. (2) gives:










J
=



qa
0




k
rec



(


k
PPr



k
PPd

+

k
PPr



)




(



n
I



p
I


-



k
PPd


k

PPd
,
eq





n

I
,
eq




)


-


qJ
X



(


k
PPd



k
PPd

+

k
PPr



)




,




(
3
)








using ζeq=krecnI,eqPI,eq/kPPd,eq from Eq. (2). The subscript eq indicates the thermal equilibrium value in the absence of bias or illumination. Similar to the Shockley Equation, quasi-equilibrium is assumed. Hence, the carrier densities at the interface (nI, pI) and contacts (nc, pc) are related via:










n
I

=


n
C


exp






δ
A



q


(


V
a

-

V
bi


)





k
b


T









and





(

4

a

)








p
I

=


p
C


exp






δ
D



q


(


V
a

-

V
bi


)





k
b


T






,




(

4

b

)








where δDA=1 are the fractions of the potential dropped across the donor (D) and acceptor (A) layers, respectively. Here, Va is the applied bias, kb is Boltzmann's constant, and T is the temperature. These relations are strictly valid only when J=0, but are a good approximation at low current when J is much smaller than either of its drift or diffusion components.


Use of Eqs. (4) in Eq. (3) yields:










J
=



qa
0



k
rec



n
C




p
C



(

1
-

η
PPd


)




exp


(



-

qV
bi


/

k
b



T

)




{


exp


(



qV
a

/

k
b



T

)


-


k
PPd


k

PPd
,
eq




}


-

q






η
PPd



J
X




,




(
5
)








where ηPPd=kPPd/(kPPd+kPPr) is the PP dissociation probability. Assuming detailed balance of the charge density adjacent to an injecting contact:

nc=ƒ(Fc,T)NLUMO exp(−ϕc/kbT),  (6)

where NLUMO is the density of states (DOS) at the acceptor LUMO, and Fc is the electric field at the cathode contact. The analogous relation involving the injection barrier, ϕa, (see FIG. 1a) exists for holes at the anode, with NHOMO as the DOS at the donor HOMO. The term, f(Fc,T) is dominated by Schottky barrier lowering; since it is near unity except for the case of high field and/or low temperature, it's neglected here. Thus, the ideal diode equation for an organic HJ in the absence of traps is:










J
=




qa
0



k
rec



N
HOMO




N
LUMO



(

1
-

η
PPd


)




exp


(


-
Δ








E
HL

/

k
b



T

)




{


exp


(



qV
a

/

k
b



T

)


-


k
PPd


k

PPd
,
eq




}


-

q






η
PPd



J
X



=



J

s





0




{


exp


(



qV
a

/

k
b



T

)


-


k
PPd


k

PPd
,
eq




}


-

q






η
PPd



J
X





,




(
7
)








where ΔEHLac+qVbi from FIG. 1a. This equation carries an implicit dependence on temperature and the interfacial electric field through ηPPd (FI,T), and consequently ηPPd. The solution for these terms is provided in Calculation of the Interfacial Electric Field below.


In the absence of illumination (JX=0), Eq. (7) represents the ideal organic HJ current density-voltage relationship. Compared to the Shockley Equation, the prefactor of the bracketed term is analogous to the dark saturation current (called Js0 here), whose temperature dependence is predominantly exponential in ΔEHL. For Va<0, the interfacial field is high and aids dissociation. Thus, kPPd grows larger than kPPd,eq and the reverse saturation current increases with increasing reverse bias, as typically observed in organic HJs. Under forward bias, kPPd is similar to or less than kPPd,eq and the current density increases exponentially with an ideality factor n=1. In this case, Eq. (7) reduces to the familiar:

J=Js0{exp(qVa/kbT)−1}−qηPPdJX  (8)

frequently used to model organic HJ solar cells. As expected, the photocurrent (last term on the right) is directly proportional to the PP dissociation efficiency, which diminishes with increasing forward bias.


Solving Eq. (7) for the open-circuit voltage gives:










qV
oc

=


Δ






E
HL


-


k
b


T





ln






(


k
PPr


k
PPd


)





k
rec



N
HOMO



N
LUMO




J
X

/

a
0






.







(
9
)








Equation (9) predicts that VOC should increase with decreasing temperature, and that it should increase as the logarithm of JX (proportional to intensity) with slope kbT/q. When polaron pairs are strongly coupled to the ground state, kPPr is large and VOC decreases. At VOC, the interfacial field, FI, is small or even reverses direction (i.e. aiding recombination). In this case, kPPd takes on its zero-field value, kPPd≈[3krec/4πa03]exp(−EB/kbT), where EB is the polaron pair binding energy.


Recognizing 4Baa3/3 as the volume occupied by a polaron pair, this relationship may be inserted into Eq. (9) to give:











qV
oc

=


(


Δ






E
HL


-

E
B


)

-


k
b


T





ln






k
PPr



N
HOMO



N
LUMO




ζ
max




J
X

/

a
0









,




(
10
)








where ζmax is the maximum PP density that can be sustained at the interface (i.e., all states are occupied).


According to Eq. (10), the maximum possible open circuit voltage, Vocmax, is equal to the interfacial energy gap, less the polaron pair binding energy. This situation, arrived at empirically by Rand, et al. Phys. Rev. B 75, 115327 (2007), is reached in the limit T→0 or when JX increases to JX=a0kPPrNHOMONLUMOmax. Since ζmax˜NHOMO˜NLUMO are all on the order of the molecular density, this final condition becomes JX=a0kPPrNHOMO, which states that the exciton current is delivering the maximum density of charge that recombination at the interface can sustain. It is not possible to increase JX beyond this limit since, with all states filled, the net radiative transition rate approaches zero and the device becomes transparent. Accordingly, Eq. (10) may be viewed as an organic analog to the band-gap limitation of Vocmax in inorganic solar cells established by Shockley and Queisser, J. Appl. Phys. 32, 510 (1961).









TABLE 1







Definition of variables









Variable
Definition
Units










Trap-free









JX
exciton current density reaching the HJ
cm−2 s−1


ζ
polaron pair density at the HJ
cm−3


a0
polaron pair spatial extent
cm


krec
free carrier bimolecular recombination coefficient
cm3 s−1


kPPd
polaron pair dissociation rate
s−1


kPPr
polaron pair recombination rate
s−1


0PPd
polaron pair dissociation efficiency



δA, δD
fractions of the potential dropped across acceptor




and donor layers



nI, pI
free electron and hole densities at the HJ interface
cm−3


nC, pC
free electron and hole densities adjacent to
cm−3



respective contacts



NHOMO,
donor HOMO and acceptor LUMO edge densities of
cm−3


NLUMO
states








Exponential Trap Distribution









nlt, plt
trapped free electron and hole densities at the HJ
cm−3



interface



Tt,D,
characteristic temperatures for hole and electron trap
K


Tt,A
distributions in the donor and acceptor



HD, HA
band edge trap densities in the donor and acceptor
cm−3


krec,n,
recombination coefficients for recombination of free
cm3 s−1


krec,p
electrons with trapped holes and vice versa



nD, nA
ideality factors due to trap limited recombination in




the donor and acceptor









Calculation of the Interfacial Electric Field


The interfacial electric field is required to calculate the polaron pair dissociation efficiency, ηPPd, in Eqs (7) and (15). Assuming that free carriers determine the potential distributions in their respective layers (e.g., holes (p) in the donor, and electrons (n) in the acceptor), then in quasi-equilibrium, pμF≈Dp′(x), and the Poisson equation gives FF′−(D/μ)F″=0. Here, D is the diffusion coefficient, μ is the mobility of the charge under consideration, and F is the electric field. Assuming the Einstein relation between D and μ, and defining β=q/kbT, these equations may be solved as in D. Gupta, et al., Appl. Phys. Lett. 92, 093301 (2008) or they may be integrated once and solved as a Riccati equation (C. M. Bender and S. A. Orszag, Advanced Mathematical Methods for Scientists and Engineers (McGraw-Hill, New York, 1978)). Defining the constant CI=FI2−2qpI/(βε), the following results:










F


(
x
)


=


[


1

F
1


-


β





x

2


]


-
1






(
A1a
)







p


(
x
)


=



βɛ

2

q




[


1

F
1


-


β





x

2


]



-
2






(
A1b
)







when C1=0 and:















F


(
x
)


=



C
1






1
+


2


(


F
I

-


C
1



)





(


F
I

+


C
1



)



exp


(


-


C
1




β





x

)



-

(


F
I

-


C
1



)












(

A





2

a

)







p


(
x
)


=



2

β





ɛ

q



C
1



exp


(



C
1



β





x

)






(


F
I

-


C
1



)



(


F
I

+


C
1



)




[


(


F
I

+


C
1



)

-


(


F
I

-


C
1



)



exp


(



C
1



β





x

)




]

2







(
A2b
)








when C1≠0, for both positive and negative values of C1. As described herein, FI and pI are the electric field and hole density at the donor-acceptor interface located at position x=0. Since both the hole density at the interface and at the contact (x=xc) are known from Eqs. (4) and (6), Eq. (A2b) can be solved implicitly for FI.


For the case of an exponential trap distribution, assuming pt>>p, trapped carriers dominate the field distribution, and the Poisson equation is F′≈qpt/ε. Since quasi-equilibrium maintains for free carriers, the analog of Eq. (11) for holes is used to obtain FF′−(DID/μ)F″=0. Thus, Eqs. (A1) and (A2) remain valid upon substituting β→β/ID, where Eq. (A2b) now refers to the trapped hole density. Using Eq. (11) to relate the trapped (pc,t) and free (pc) hole densities at the anode, with the latter given by Eq. (6), pIt=pc,t exp[−qδD(Va−Vbi)/lDkbT] and Eq. (A2b) is used again to solve for FI.


After determining FI, the polaron pair dissociation efficiency is calculated using the Onsager-Braun model, summarized here. Given an initial PP separation, a0, the dissociation rate kPPd is:











k
PPd

=


3

4

π






a
0
3





k
rec



exp


(



-

E
B


/

k
b



T

)






J
1



[

2




-
2


b



]





-
2


b





,




(
A3
)








where b=q3FI/(8πεkb2T2) and JI is the Bessel function of order one. To account for disorder, a range of initial separations was taken according to the normalized distribution function ƒ(x)=4/(√{square root over (πa03)})x2 exp(−x2/a02) and integrate it over the entire current density expression, since kPPd appears independently in addition to ηPPd.


The validity of Eq. (A3) for negative fields (i.e., when FI reverses direction under sufficient forward bias) is unclear, and its numerical evaluation becomes unstable for FI<104 V/cm. Analogy to the Onsager theory (M. Pope and C. E. Swenberg, Electronic Processes in Organic Crystals and Polymers (Oxford University Press, New York, 1999)), suggests that kPPd retains its zero field value at FI<0, however, this discontinuity seems unphysical, and Monte-Carlo simulations show that the dissociation efficiency at a D-A heterojunction does indeed decrease for fields directed toward the HJ. Thus, when FI<0, the zero field kPPd is taken and the PP binding energy is augmented with the additional barrier −qfIrc cos θ, where rc=q2/(4πεkbT) is the Onsager radius. Averaging over the forward half-space of angles −π/2<θ<π/2 between the PP separation vector and the HJ normal accounts for the distribution of initial PP orientations.


Most organic solar cells are characterized by significant molecular disorder in both the donor and acceptor layers, leading to a broad density of states with an approximately Gaussian energetic distribution at the HOMO and LUMO levels. Evidence also suggests that ground state interactions between the donor and acceptor further broaden this distribution near the HJ interface. The low energy tail of the DOS, where most charge carriers reside, can be modeled as an exponential distribution of traps. This changes the kinetics of recombination at the interface, and introduces an ideality factor n>1 in the diode equation, as shown below.


Assuming an exponential trap distribution with characteristic trap temperature, Tt,A, in the acceptor, the trapped (nt) and free (n) electron densities are related approximately via:











n
t




H
A



exp


(



E
Fn

-

E
LUMO




k
b



T

t
,
A




)







H
A



(

n

N
LUMO


)



1
/

l
A




,




(
11
)








where HA is the density of trap states at the acceptor LUMO, EFn is its electron quasi-Fermi energy, ELUMO is the LUMO energy of the acceptor, and IA=Tt,A/T. A similar relationship holds for the trapped hole density, pt, in the donor. Assuming that the trapped carrier density significantly exceeds the free carrier density, Eq. (3) becomes:










J
=




qa
0



(


k
PPr



k
PPd

+

k
PPr



)




{



k

rec
,
n




(



n
I



p
It


-



k
PPd


k

PPd
,
eq





n

I
,
eq




p

It
,
eq




)


+


k

rec
,
p




(



p
I



n
It


-



k
PPd


k

PPd
,
eq





p

I
,
eq




n

It
,
eq




)



}


-


qJ
X



(


k
PPd



k
PPd

+

k
PPr



)




,




(
12
)








where recombination now occurs primarily at trap states, and krec,n and krec,p are the rate constants describing recombination at the HJ between a free electron in the acceptor (nI) with a trapped hole in the donor (pIt), and vice versa.


Using Eqs. (4), (11) and (12) gives:










J
=




qa
0



(

1
-

η
PPd


)




{



k

rec
,
n




N
LUMO



H
D




exp


(



-

α
D


/

k
b



T

)




[


exp


(



qV
a

/

n
D




k
b


T

)


-






k
PPd


k

PPd
,
eq




]



+








k

rec
,
p




N
HOMO



H
A




exp


(



-

α
A


/

k
b







T

)


[






exp
(







qV
a

/

n
A








k
b






T

)

-






k
PPd


k

PPd
,
eq




]



}


-

q






η
PPd



J
X




,








where




(

13

a

)













α
D

=



Δ






E
HL



n
D


+




l
D

-
1


l
D




(



δ
D



ϕ
c


-


δ
A



ϕ
a



)




,








and





(

13

b

)












α
A

=



Δ






E
HL



n
A


+




l
A

-
1


l
A





(



δ
A



ϕ
a


-


δ
D



ϕ
c



)

.








(

13

c

)








The ideality factors, nA and nD are given by:











n
A

=


l
A




δ
D



(


l
A

-
1

)


+
1



,




and




(

14

a

)







n
D

=



l
D




δ
A



(


l
D

-
1

)


+
1


.





(

14

b

)








More compactly, Eq. 13 becomes the ideal diode equation in the presence of traps:









J
=



J
sD



[


exp


(



qV
a

/

n
D




k
b


T

)


-


k
PPd


k

PPd
,
eq




]


+


J
sA



[


exp


(



qV
a

/

n
A




k
b


T

)


-


k
PPd


k

PPd
,
eq




]


-

q






η
PPd




J
X

.







(

15

a

)








which simplifies to:

J=JsD[exp(qVa/nDkbT)−1]+JsA[exp(qVa/nAkbT)−1]−PPdJX,  (15b)

when kPPd≤kPPd,eq under forward bias (c.f. Eq. 8). Here JsD and JSA are the dark saturation currents given by the prefactors in Eq. (13).


Thus, in general, there are two sources of dark current stemming from the recombination of free electrons with trapped holes and vice versa. Each pathway produces its own ideality factor (nD and nA, respectively) and dark saturation current, both of which depend on the balance of the voltage drop across the D and A layers as well as their characteristic trap temperatures.


In the case of a perfectly symmetric device, with identical transport properties and injection barriers for electrons and holes, Eq. (13) simplifies to:










J
=




qa
0



(

1
-

η
PPd


)



2


k

rec
,
sym




N
sym



H
sym



exp


(


-
Δ








E
HL

/

n
sym




k
b


T

)




{


exp


(



qV
a

/

n
sym




k
b


T

)


-


k
PPd


k

PPd
,
eq




}


-

q






η
PPd



J
X




,








or




(

16

a

)












J
=



J
sym



{


exp


(



qV
a

/

n
sym




k
b


T

)


-
1

}


-

q






η
PPd



J
X




,





(

16

b

)








where nsym=2Isym/(Isym+1) and the subscript sym indicates parameters that are equivalent for holes in the donor and electrons in the acceptor. Solving for Voc leads to:











qV
oc

=


Δ






E
HL


-


n
sym



k
b


T





ln





(


k
PPr


k
PPd


)





k

rec
,
sym




N
sym
2




J
X

/

a
0









,




(
17
)








resulting in a slope of nsymkbT/q when plotted versus the logarithm of intensity. The maximum open circuit voltage, Vocmax, still reduces to ΔEHL−EB, since in the limit of high light intensity, the original assumption of nt, pt>>n, p no longer holds. In this regime all traps are filled and recombination is no longer trap-limited, hence, Eq. (17) reverts to Eq. (9).



FIG. 2a shows a set of dark J-V characteristics calculated using Eq. (13) over the temperature range from 120 K to 296 K for a device using the parameters listed in Table 2, which are typical of organic photovoltaic cells. Included is a series resistance of Rs=1Ωcm2 (in which case, Va is replaced by Va−JRs) and assumed that most of the potential is dropped across the donor layer, resulting in different ideality factors according to Eq. (14). The ideality factors and their associated dark saturation currents are given in FIG. 2b. Both ideality factors increase with decreasing temperature, and n ln(Js) is nonlinear when plotted vs. I/kbT. This contrasts with previous analyses based on the Shockley equation, which predict a linear relationship for this quantity with a slope equal to −ΔEHL/2.









TABLE 2







Model Parameter Values








Parameter
Value





donor thickness = acceptor thickness
40 nm


ΔEHL
1.2 eV


Vbi
0.5 V


Tt,A = Tt,D
1000 K


HA = HD
1018 cm−3


NHOMO = NLUMO
1021 cm−3


δA
0.1


a0
1.5 nm


kPPr
1 μs−1


krec,n = krec,p = qμ/ε
ε/ε0 = 3, μ = 10−3 cm2/Vs


Rs
1 Ω cm2









In comparison to the dark current characteristics of many reported planar HJ organic solar cells, the following observation were made. In the generalized Shockley Equation, increasing dark current with reverse bias is phenomenologically accounted for by a shunt resistance, whereas here it is due to the field dependent dissociation of thermally generated PPs [i.e. kPPd=kPPd(FI)]. The two slopes routinely observed in semi-log plots under forward bias (the second is only evident at low temperature for some devices) are not the result of a shunt resistance and a single diode ideality factor, but rather of two ideality factors that reflect recombination with trapped carriers at each side of the HJ.


In forward bias, the ideality factor nD is evident in the slope of the J-V characteristics for 0<Va<0.3 V, whereas nA dominates at 0.4<Va<0.7 V, beyond which series resistance from contacts and the layer bulk limits the current. The asymmetry of the voltage dropped across the donor and acceptor layers (δA=0.1 is assumed here) is the primary cause of the difference in ideality factors, and is expected by analogy to the asymmetries in organic light emitting devices.


At low bias, the current is predominantly mediated by PPs formed from the recombination of free holes in the donor with trapped electrons in the acceptor. At higher bias (0.4<Va<0.7), the reverse process dominates, and the slope is determined by nA. As temperature decreases, carriers freeze into the trap states that constitute the tail of the DOS, and recombination becomes further trap-limited, resulting in the increase of both ideality factors with temperature (see FIG. 2b, Eq. 14).


In FIG. 3a, a calculation for the intensity dependence of Voc, Jsc, and the fill factor, FF, from Eq. (13), using the same parameters from the dark J-V characteristics of FIG. 2 is provided. The exciton current density, qJX, is directly proportional to light intensity and is approximately equal to Jsc. The FF increases and then saturates at high light intensity. The open-circuit voltage increases as the log of intensity, with the slopes nDkbT/q and nAkbT/q at low and high intensity, respectively, as expected from the double exponential form of Eq. (13).



FIG. 3b shows the temperature dependence of Voc, FF, and the polaron pair dissociation efficiency at short circuit, ηPPd. The open-circuit voltage increases by ˜0.25 V at low temperature, as compared to its value at T=300K, whereas FF steadily declines for T≤200 K. The polaron pair dissociation efficiency decreases substantially with decreasing T from its room temperature value of ηPPd˜0.6. The theoretical trends in FIG. 3 have been observed in both small molecular and polymer organic solar cells.


The trap-free and trap-limited ideal diode equations derived here (Eqs. 7 and 15, respectively) are similar in form to the Shockley Equation but the interpretation of the fitting parameters (ideality factors, saturation current densities) is different due to the different nature of excited states in organic and inorganic semiconductors. In inorganic semiconductors, excitation produces free carriers directly, whereas in organic semiconductors, the result of photon absorption is a tightly bound exciton that has a very low probability of dissociating unless charge transfer is initiated at a D-A HJ to form a more loosely bound polaron pair. Thus, in contrast to inorganic p-n junctions where current is due to drift-diffusion and/or recombination within the depletion region, the current in organic HJs depends on polaron pair recombination and dissociation that occur over a small volume at the HJ interface.


In the absence of shunt paths or other junction defects, generation and recombination via interfacial polaron pairs is the sole source of current flow. This assumption is reasonable at low bias since PPs provide the lowest energy recombination pathway in D-A HJs with Type II (i.e., staggered) HJ energy level offsets, and hence should be the most heavily populated state at small quasi-Fermi level separations. As described in more detail below, this is not the case at high bias (Va>2 V), where both SubPc and C60 fluorescence are detected, indicating injection beyond the HJ (i.e. of holes into the acceptor bulk, and electrons into the donor), with subsequent exciton formation and recombination in each of the layers. This assumption thus does not apply to organic light emitting devices, which generally operate at high bias (Va>2 V) and are designed to minimize energy level offsets to maximize injection into, and exciton formation in the emissive layer bulk.


As in the Shockley Equation, quasi-equilibrium conditions are assumed, implying that J is negligible compared to either of its drift or diffusion components. This assumption is rigorous at low currents, but breaks down for J>0.1 A/cm, where the voltage drop due to contact and bulk series resistance becomes significant. Beyond these assumptions, the trap-free current-voltage relationship and its consequences for VOC established by Eqs. (9) and (10) are of general validity.


Disorder-induced and other trap states complicate the trap-free picture, as the recombination kinetics depend on the particular trap distribution chosen. Regardless of this choice, the double exponential form of the J-V relation given by Eq. (15) is a general result, as it stems from the two possible recombination pathways of free electrons with trapped holes, and vice versa.


Shockley-Hall-Read (SHR) theory provides an alternate way to describe trap-limited recombination. Although this approach has previously been applied to organic solar cells, it assumes detailed balance of traps active for both carrier types in a single material. In contrast, the organic HJ involves two different sets of traps, active for holes and electrons in the donor and acceptor layers, respectively. The Langevin bimolecular approach of Eq. (12) is well established for organic semiconductors, and its application here is more appropriate than the SHR approach.


The fits of FIG. 4 demonstrate that Eq. (18) is an accurate description of the J-V characteristics, as compared to the generalized Shockley Equation for CuPc/C60 and SubPc/C60 solar cells. In FIG. 5, the similarity in nA and nA ln(JSA) between CuPc/C60 and SubPc/C60 suggests that the current component involving these terms is due to recombination with trapped electrons in the C60 acceptor common to both HJs. Analogously, the differing nD and nD ln(JsD) suggests recombination with trapped holes occupying the different density of states characteristic of CuPc and SubPc.


The values of n and Js from FIG. 5 give rise to the intensity and temperature dependencies of Voc in FIG. 6. Only one slope is evident in FIG. 6a for each device. A second slope is expected, as predicted in FIG. 3a to become apparent at low intensity or low temperature. In FIG. 6b, the deviation in calculated Voc from the data at low temperature is due to the drop in polaron pair dissociation efficiency, as expected from the Onsager-Braun model and replicated in FIG. 3b.


While the trends predicted in FIG. 2b are qualitatively similar to those in FIG. 5, their primary difference is the magnitude and degree of change predicted for the ideality factors. This is partially attributed to the choice of an exponential trap distribution, which, although useful in enabling analytical results, may not accurately represent the trap distributions in the materials provided here. Furthermore, δD and δA do not cancel as in the trap-free case and hence changes in the potential distribution across the D and A layers with bias and/or temperature will also influence the ideality factors [see Eq. (14)].


A further concern is the implicit use of the conventional Einstein relation for the charge carrier mobility, which may not be appropriate for disordered semiconductors. The Einstein relation depends on the particular density of states or trap distribution function and in general increases (i.e., D/μ>kbT/q, where D is the diffusivity) with carrier density and with decreasing temperature. As the ideality factors depend directly on this change, its incorporation into the model results in an increase in the magnitude of nA and nD and leads to a stronger temperature dependence, as observed in the data. Despite this shortcoming, as shown in FIG. 3, the model qualitatively predicts the magnitudes, and intensity and temperature dependences of Voc, Jsc, and FF that are commonly observed in organic HJ solar cells using the typical material parameters listed in Table 2.


Previously, Voc was observed to saturate at T˜175 K for a number of different HJs (V. Dyakonov, Appl. Phys. A 79, 21 (2004) and B. P. Rand, D. P. Burk, and S. R. Forrest, Phys. Rev. B 75, 115327 (2007)), whereas other reports (H. Kumar, P. Kumar, N. Chaudhary, R. Bhardwaj, S. Chand, S. C. Jain, and V. Kumar, J. Phys. D 42 (2009)), and the data in FIG. 6b show a monotonic increase with decreasing temperature. The reason for this discrepancy may lie in the difficulty of accurately measuring device temperature under illumination on a thermally insulating glass substrate. A plateau in Voc at intermediate temperature is not predicted in FIG. 3b using the parameters in Table 2. However, in some asymmetric HJs with different D and A layer thicknesses, trap temperatures, and injection barriers (i.e. ϕa and ϕc), a plateau is found to emerge between T=100 K and 150 K, as observed.


As shown in FIG. 7, the solar power conversion efficiency depends critically on the polaron pair recombination rate, kPPr. The decline in Jsc is due to the decrease in PP dissociation efficiency with increasing kPPr, and the drop in Voc follows Eq. (9). Physically, kPPr describes either the direct recombination of the D+−A charge transfer state, or its indirect recombination, e.g., via a lower-lying triplet state. Estimates for this parameter generally range between 1 ns−1 and 1 μs−1, as supported by polaron transient absorption kinetics, electronic modeling of D-A excited states, and Monte-Carlo simulations of geminate recombination.


Finally, the ‘S-kink’ behavior found many organic HJ cells in FIG. 8a is explained schematically in FIG. 8b. The top panel illustrates the shift in zero-crossing of the internal field (F1<0 favors dissociation) toward positive applied bias with increasing Vbi. It leads to a similar shift in the PP dissociation efficiency, ηPPd, and hence to a concomitant shift in the photocurrent (c.f. Eq. 7).


The lower panel shows that the total current is the sum of the dark and photocurrents. The ‘S-kink’ appears when the slope of the photocurrent, as determined by ηPPd, is large in the fourth quadrant at low bias, and when the dark current is not yet significant. At room temperature, this typically occurs for a combination of large kPPr (e.g. ≥107 s−1) and small Vbi (e.g. ≤0.3 V). In contrast, at low temperature this feature emerges for a much broader range of Vbi and kPPr, since under these conditions, kPPd is small, shifting the θPPd curve (see top panel) toward low applied bias. An ‘S-kink’ is thus expected to be a general feature in the J-V characteristic of organic HJ solar cells under illumination at low temperature.


Understanding and minimizing the PP recombination rate is clearly important for achieving high efficiency organic solar cells. The kPPr employed in the Onsager-Braun model may in fact be the result of multiple recombination pathways, potentially involving intermediate states, and depending on such factors as mutual orientation and orbital overlap between the donor and acceptor molecules at the heterointerface. For example, Perez et. al. linked the degree of π-π overlap and intermolecular interaction between the donor and acceptor to the dark saturation current and Voc in a broad selection of molecules. Both of these factors depend directly on the magnitude of kPPr. Additionally, rapid intersystem crossing between singlet and triplet PP configurations due to the small exchange splitting can make recombination via a triplet exciton favorable if such a state exists below the PP energy.


Based in part on this understanding, the present disclosure further describes the nature of kPPr at CuPc/C60 and SubPc/C60 HJs. For example, in one embodiment, the present disclosure relates to PP dynamics in two archetype planar HJ organic solar cells consisting of copper phthalocyanine (CuPc)/C60 and boron subpthalocyanine chloride (SubPc)/C60 donor-acceptor (D-A) pairs.


In one embodiment, the present disclosure provides for an organic photosensitive optoelectronic device comprising:


two electrodes comprising in superposed relation;


at least one donor material, and


at least one acceptor material,


the at least one donor material and at least one acceptor material forming at least one photoactive region between the two electrodes,


wherein the at least one donor material and at least one acceptor material are ordered in their bulk and comprise a disordered region at their interface.


In one embodiment, the disordered region comprises a non-crystalline region, a steric bulk region, or combinations thereof.


The disordered region may comprise disorder in the donor material, disorder in the acceptor material, or combinations thereof.


In one embodiment, the disordered region comprises less than ten monolayers of the donor material and/or the acceptor material.


In one embodiment, the photo-active region forms at least one of a mixed heterojunction, bulk heterojunction, nanocrystalline-bulk heterojunction, and hybrid planar-mixed heterojunction.


In one embodiment, the order in the bulk comprises at least one of single crystal, nanocrystalline or polycrystalline material.


Non-limiting examples of the donor that may be used herein include subphthalocyanine (SubPc), copper pthalocyanine (CuPc), chloroaluminium phthalocyanine (CIAIPc), tin phthalocyanine (SnPc), pentacene, tetracene, diindenoperylene (DIP), and squaraine (SQ).


Non-limiting examples of the acceptor that may be used herein include C60, C70, 3,4,9,10-perylenetetracarboxylicbis-benzimidazole (PTCBI), PCBM, PC70BM and hexadecafluorophthalocyanine (F16CuPc).


In one embodiment, the at least one electrode comprises transparent conducting oxides or transparent conducting polymers.


Non-limiting examples of conducting oxides that may be used herein include indium tin oxide (ITO), tin oxide (TO), gallium indium tin oxide (GITO), zinc oxide (ZO), and zinc indium tin oxide (ZITO), and the transparent conductive polymers comprise polyanaline (PANI).


In one embodiment, the at least one electrode comprises a metal substitute, a non-metallic material or a metallic material chosen from Ag, Au, Ti, Sn, and Al.


The organic photosensitive optoelectronic device as described herein may further include at least one exciton blocking layer.


Non-limiting examples of the exciton blocking layer that may be used herein include bathocuproine (BCP), bathophenanthroline (BPhen), 3,4,9,10-perylenetetracarboxylicbis-benzimidazole (PTCBI), 1,3,5-tris(N-phenylbenzimidazol-2-yl)benzene (TPBi), tris(acetylacetonato) ruthenium(III) (Ru(acaca)3), and aluminum(III) phenolate (Alq2 OPH).


In one embodiment, the organic photosensitive optoelectronic device is an organic photodetector. In another embodiment, the organic photosensitive optoelectronic device is an organic solar cell.


In one embodiment, the present disclosure also provide for a method of making an organic photosensitive optoelectronic device, the method comprising depositing onto a substrate,


at least one donor material,


at least one acceptor material, and


at least one electrode,


wherein the at least one donor material and at least one acceptor material are deposited to form at least one photoactive region between the two electrodes,


wherein the at least one donor material and at least one acceptor material are ordered in their bulk and comprise a disordered region at their interface.


In one embodiment, the disordered region is formed by forming a non-crystalline region, a steric bulk region, or combinations thereof. In another embodiment, the disordered region is formed by disordering the donor material, disordering in the acceptor material, or disordering both the donor and acceptor materials.


In one embodiment, the disordered, non-crystalline region comprises less than ten monolayers each of the donor material and the acceptor material.


The disordered, non-crystalline region described herein may be formed by changing the deposition conditions used to deposit the donor and acceptor materials. In another embodiment, the disordered, non-crystalline region is formed by modifying the molecules used as the donor and acceptor materials.


In one embodiment, the depositing occurs in a vacuum.


Weak emission from exciplex states at the D-A heterointerface, and intensity modulated photocurrent spectroscopy (IMPS) to monitor changes in the PP density as functions of bias and illumination intensity are also described.


kPPr depends on electric field for both material pairs, and that it is significantly larger in CuPc/C60 than in SubPc/C60. To understand these differences, an expression for kPPr was found, which qualitatively explains the exciplex luminescence and IMPS trends observed. Additionally, further observations support the rate-based description of polaron pair dynamics used to derive the organic HJ ideal diode equation described above.


Intensity modulated photocurrent spectroscopy (IMPS) is a useful tool to explore the balance between photocurrent and recombination current that has previously been applied to electrochemical and dye-sensitized solar cells. It involves superimposing a small optical modulation onto the steady-state excitation incident on a cell, and measuring both the magnitude and phase of the resulting photocurrent. This yields a complex external quantum efficiency, {tilde over (Φ)}, that depends on the optical modulation frequency in addition to, for example, the steady-state light intensity and cell electrical bias.


The time-dependent polaron pair and interfacial free carrier densities are described via:












d





ζ

dt

=



j
X


a
0


-


k
PPr


ζ

-


k
PPd


ζ

+


k
rec



n
1



p
1




,




and




(
19
)








dp
1

dt

=



k
PPd


ζ

-


k
rec



n
1



p
1


-


J

qa
0


.






(
20
)








Here, JX is the exciton current density diffusing to the HJ which undergoes quantitative charge transfer and contributes to ζ the polaron pair density. The electronic charge is q, and nI and pI are the free electron and hole densities at the HJ interface, respectively. The heterojunction volume is defined by the average polaron pair separation, a0, and krec is the bimolecular Langevin recombination coefficient. Polaron pairs dissociate at rate kPPd given by the Onsager-Braun model, and recombine to the ground state at rate kPPr, resulting in a dissociation efficiency ηPPd=kPPd/(kPPd+kPPr). The intrinsic thermal equilibrium PP density, ζeq<<ζ under IMPS illumination, and has therefore been omitted from Eq. (19). Additionally, J is redefined in Eq. (20) as positive in the direction of photocurrent, consistent with IMPS convention.


Due to the harmonic perturbation in IMPS illumination at frequency, ω, the photogenerated exciton diffusive flux reaching the interface becomes JX→JX(1+δeiωτ), which assumes that the exciton lifetime is small compared to ω−1. Similarly, the steady-state and time varying components of Eqs. (19) and (20) are separated as, for example, J→J+{tilde over (j)}eiωτ and pI→pI+{tilde over (p)}Ieiωτ, where the tilde denotes a complex, small-signal quantity. Retaining only first order harmonic terms Eq. (20) gives:










i





ω



p
~

1


=



k
PPd



ζ
~


-


k
rec



n
1




p
~

1


-



j
~


qa
0


.






(
21
)








The layer bulk between the interface and the contact may be approximated by its equivalent, parallel resistance-capacitance (RC) circuit, where the impedance is {tilde over (Z)}=R/(1+iωτ) and τ=RC. Hence, relating the areal charge density at the interface, {tilde over (Q)}I=qa0{tilde over (p)}I, to the small-signal current density, {tilde over (j)} via the usual capacitive relationship:

{tilde over (Q)}1=(C{tilde over (Z)}){tilde over (j)}.  (22)


Equations (19), (21) and (22) are solved for {tilde over (j)} and normalized by qδJX to give the complex IMPS response:











Φ
~

=



n
PPd



(

1
+

i





ωτ


)



1
+


k
rec



n
1



τ


(

1
-

n
PPd


)



+

2

i





ωτ




,




(
23
)








where the modulation frequency is assumed to be low compared to the polaron pair dynamics (i.e. ω<<kPPd+kPPr). Equation (23) predicts a transition from Im({tilde over (Φ)})<0 to Im({tilde over (Φ)})>0 when krecnIτ(1−θPPd)>1. Hence, an increase in the interfacial free carrier density can cause a crossover from phase lag to phase advance in the IMPS spectrum.


The trap-free model detailed above was used to determine the steady-state quantities (e.g. nI and θPPd) in Eq. (23), and assume a symmetric planar HJ for simplicity. Observations discussed below suggest that kPPr depends on the electric field at the interface, FI. As a first approximation, a field-dependent shift (qFIr) in the PP potential is added to existing energy barriers that determine the rate of recombination, viz.:










k
PPr

=


k

PPr





0





exp


(



qF
1


r



k
b


T


)


.






(
24
)







Here, kPPr0 is the zero-field recombination rate, r is the PP separation, and FI is the electric field at the interface, as calculated in Calculation of the Interfacial Electric Field above. As before, kPPr and kPPd over a Gaussian distribution of separations is evaluated, taking the average value a0=1.5 nm. Equation (24) describes a physical situation analogous to field-induced energy barrier lowering such as the Schottky effect, or Frenkel-Poole emission in conventional semiconductors.


EXAMPLES
Example 1

Conventional CuPc/C60 and SubPc/C60 cells were fabricated on glass substrates coated with a transparent indium tin oxide (sheet resistance of ˜15 ohms/square) anode pre-patterned into 1 mm wide stripes. All organic materials were purified by thermal gradient sublimation prior to use. Following a thorough solvent degrease of the substrate, layers were deposited sequentially by thermal evaporation in a chamber with base pressure ˜10−7 Torr. The devices consist of 20 nm thick CuPc or 11 nm thick SubPc, followed by 40 nm thick C60, 10 nm thick bathocuprine, and a 100 nm thick Al cathode, deposited through a shadow mask as 1 mm stripes positioned orthogonally to the patterned anode to form 1 mm2 device areas.


Cell efficiency was characterized at room temperature in air under simulated AM1.5G illumination, and found to be comparable to previously reported devices using these materials combinations K. L. Mutolo, E. I. Mayo, B. P. Rand, S. R. Forrest, and M. E. Thompson, J. Am. Chem. Soc. 128, 8108 (2006); and M. D. Perez, C. Borek, S. R. Forrest, and M. E. Thompson, J. Am. Chem. Soc. 131, 9281 (2009). Temperature dependent measurements were conducted in an evacuated, closed-cycle He cryostat using illumination from the 2=496 nm line of an Ar+ laser with an intensity of 30 mW/cm2. The temperature of each device was measured with a Ge thermistor soldered to the substrate surface. Current-voltage characteristics were obtained using an Agilent 4156 semiconductor parameter analyzer at a voltage sweep rate of 0.2 V/s to minimize unwanted capacitive effects. The current-voltage characteristics were fit using a nonlinear least squares trust-region algorithm with statistical weighting.


In FIGS. 4a and b, the J-V data taken over a range of temperatures are fit for both CuPc/C60 and SubPc/C60 cells to the ideal diode equation of Eq. (15) modified to include the effect of series resistance, Rs:

J=JsD{exp[q(Va−JRs)/nDkbT]−1}+JsA{exp[q(Va−JRs)/nAkbT]−1}−qηPPdJX.  (18)


The results of Eq. (18) matched the data over the entire temperature range, whereas those calculated using the generalized Shockley Equation and plotted in FIGS. 4c and d, show increasingly poor agreement at low temperature and voltage.



FIG. 5 presents the ideality factors and dark saturation currents extracted from the fits in FIGS. 4a and 4b. Both sets of ideality factors increase with decreasing temperature, similar to the prediction of FIG. 2b. nA was assigned to the overlapping ideality factors of CuPc/C60 and SubPc/C60 as above. In FIG. 5b, the quantity n ln(Js) showed a non-linear Arrhenius dependence in all cases, where the trend was also similar to that predicted in FIG. 2b.



FIG. 6a shows the intensity dependence of Voc under simulated AM1.5G solar illumination. Lines are calculated using Eq. (18) along with the ideality factors and dark saturation currents determined from data in FIG. 5. The photocurrent, qηPPdJX, is approximated by the short-circuit current density, Jsc, measured at each intensity. The slopes are dominated by nAkbT/q for both CuPc/C60 and SubPc/C60 devices; that of nDkbT/q only becomes evident at lower intensity.


The temperature dependence of Voc under λ=496 nm laser illumination is shown in FIG. 6b. The solid lines are also calculated from Eq. (18) using the temperature dependent n and Js from FIG. 5, with the photocurrent term given by Jsc measured at T=300 K. Since the diode parameters characterize the dark current temperature dependence, inference provides that the deviation at low temperature is due to a decrease in the polaron pair dissociation efficiency, ηPPd. The ratio, ηPPd(T)/ηPPd (300K), calculated from this difference, is also shown in FIG. 6b (dashed lines), and drops by a factor of ˜40 between room temperature and T=114 K.


In FIG. 7a, the normalized power efficiency, ηp, is calculated along with Voc, Jsc, and FF as a function of the polaron pair recombination rate using the model as in FIGS. 2 and 3. Current-voltage characteristics corresponding to several values of kPPr are shown in FIG. 7b. All of the cell efficiency parameters decrease with increasing recombination rate, and the power conversion efficiency drops most steeply for kPPr>1 μS−1. In contrast, FIG. 8a shows the J-V dependence on built-in potential, Vbi, for the same model as in FIG. 7b, with kPPr=(20 ns)−1. Small Vbi result in an ‘S’ shaped kink as observed in some small molecule and polymer organic solar cells.


Example 2

Conventional CuPc/C60 and SubPc/C60 cells were fabricated and characterized using methods similar to those as described in Example 1. In brief, the cells consisted of a 20 nm thick CuPc or 11 nm thick SubPc layer vacuum-deposited at a chamber base pressure of 10−7 Torr on an indium tin oxide covered glass substrate, followed by a 40 nm thick layer of C60, 10 nm thick bathocuprine (BCP), and a 100 nm thick Al cathode.


Room temperature IMPS was conducted in a sealed, electrically and optically accessible chamber continuously purged by dry N2 to minimize device degradation. Low temperature data were obtained in a closed-cycle, evacuated He cryostat, and temperature was measured using a Ge thermistor soldered to the substrate surface. An Ar+ laser (λex=496) nm with an intensity of 30 mW/cm2 was used for excitation; the trends observed were independent of excitation wavelength over the range from λex=476 nm to 514 nm. A small fraction of the incident beam (<5%) was acousto-optically modulated at frequencies ranging from 1 Hz to 100 kHz, and the photocurrent was measured with a lock-in amplifier. Impedance spectra were used to measure the time constant, τ, for each device. Although τ varies with bias and illumination intensity, it was found that τ˜5 μs is typical for both cells near Voc under the 30 mW/cm2 IMPS illumination. This value was used for all simulations.


Luminescence was characterized from the same N2 purged sample cell using the λex=442 nm line of a HeCd laser for excitation. This wavelength coincides with a minimum in SubPc absorption to prevent its fluorescence from obscuring the weak exciplex emission. The excitation beam was split into a weak signal beam and a variable intensity bias beam. The transverse electric polarized signal beam was acousto-optically modulated at 2 mW/cm2 rms intensity, and was incident on the sample at an angle of 60°. Steady-state bias illumination was incident at 40°, and its intensity was varied in the range from 0.8 to 30 mW/cm2. Emission was collected normal to the sample with an ƒ/2 lens and focused into a multimode fiber where it was directed to a monochromator equipped with a photomultiplier. Luminescence at this frequency was obtained as a function of sample electrical bias and optical bias intensity using a lock-in amplifier.


Luminescence data were obtained by staggering the data points in random fashion to avoid systematic error due to sample degradation over the course of the measurement. All luminescence trends were independent of modulation frequency in the range of 100 Hz to 100 kHz; data reported were obtained at 10 kHz. Electroluminescence intensity was measured using a calibrated Si photodiode, and spectra were obtained with a fiber-coupled spectrometer under DC current drive from a semiconductor parameter analyzer.


Emission spectra for both CuPc/C60 and SubPc/C60 devices are shown in FIG. 9a. The peak at λem=615 nm in both the photoluminescence (PL) and electroluminescence (EL) spectra of SubPc/C60 corresponds to SubPc fluorescence. The electroluminescence peak at λem=715 nm is primarily due to fullerene fluorescence, as is the long wavelength shoulder at λem=800 nm. Overlapping emission from SubPc structural defects or aggregates may also contribute to the peak at λem=715 nm.


In contrast, the PL spectra of both SubPc/C60 and CuPc/C60 showed a significant redshift and broadening from the fullerene reference. This is ascribed to an interfacial Pc/C60 exciplex in each case, noting that it is likely to be superimposed on a bulk fullerene fluorescence background. Similar exciplex emission red-shifted from the fullerene peak has been reported in porphyrin-fullerene systems, and has recently been established for several phthalocyanine-fullerene dyads as well.



FIG. 9b shows PL spectra obtained from a quartz/SubPc (11 nm)/C60 (X nm)/BCP (10 nm) stack with excitation incident on SubPc through the substrate, where X=10, 20, 30, and 40 nm. The inset shows the intensity of SubPc fluorescence and exciplex/fullerene emission integrated over the wavelength intervals λem={575, 670 nm} and λem={730, 800 nm}, respectively, as a function of C60 thickness. Lines indicate predictions based on an optical interference/diffusion model, where the exciton generation profile in each layer is calculated using the transfer matrix approach along with the optical constants of each layer. Exciplex emission intensity is assumed to be proportional to the diffusive exciton flux toward the heterointerface from both the SubPc and C60 layers, whereas the fluorescence intensity from each layer is proportional its integrated exciton density. Literature reported diffusion lengths of 8 nm and 40 nm are used for SubPc and C60, respectively.



FIG. 10a shows the emission intensity at λen=615 nm and λem=730 nm for the CuPc/C60 and SubPc/C60 HJs, respectively, as functions of electrical bias under a steady state illumination intensity of Ibias=10 mW/cm2. Whereas SubPc fluorescence increases only slightly with increasing bias, the exciplex/fullerene emission at λem=730 nm increases linearly, peaks at Va˜0.5 V and then sharply decreases. This contrasts with the behavior of CuPc/C60, whose exciplex/fullerene emission decreases under reverse bias for Va≥−0.4 V, and continues before leveling off at Va>0.5 V.



FIG. 10b shows the bias dependence of SubPc/C60 exciplex/fullerene emission for Ibias=30 and 0.8 mW/cm2. For the low intensity case, the probe modulation was decreased from 2 to 0.2 mW/cm2 rms to maintain it as a small fraction of Ibias, resulting in a proportionately lower luminescence signal. The corresponding current-voltage characteristics, obtained at the same time as the luminescence, are also plotted in FIG. 10b (solid lines). The emission peaks determined by linear fitting of the two slopes evident in each trend, appear to correspond to the open-circuit voltage as identified by the minimum in the log(J) vs. V data.


Intensity modulated photocurrent spectra are shown for SubPc/C60 in the Cole-Cole plot of FIG. 11a, with frequency was the implicit variable. At zero bias (short-circuit), {tilde over (Φ)} is purely real at low frequency, consistent with a steady-state external quantum efficiency of ˜18% at λex=496 nm. Negative phase accumulates with increasing frequency due to the charge carrier transit time delay from the interface to the contacts. This is typical of the IMPS response, similar to that observed in dye sensitized solar cells, although the phase lag here is smaller, and occurs at higher frequency due to the faster charge transport in organic thin films.


This behavior is maintained for increasing forward bias, with |{tilde over (Φ)}| decreasing as expected from the current-voltage characteristics. Near open-circuit conditions (Va˜Voc=0.81 V), the response changes to a phase advance. Although the photocurrent appears to lead the optical signal, this is in fact a result of the superposition of a positive photocurrent (toward the contacts), and a negative recombination current (toward the interface). According to the crossover condition, krecnIτ(1−ηPPd)>1 derived above, this indicates substantial bimolecular recombination (i.e. krecnI becomes large) at open-circuit in the SubPc/C60 cell. FIG. 11b shows that the CuPc/C60 device behaves similarly to SubPc/C60 at short-circuit. At open-circuit (Va=0.44 V), however, the CuPc/C60 cell shows no phase advance, instead maintaining its short-circuit form as |{tilde over (Φ)}|→0.



FIG. 12a shows IMPS spectra calculated from Eq. (23) assuming an initial polaron pair separation a0=1.5 nm and two different constant values for the PP recombination rate, kPPr. When kPPr=(100 ns)−1 (black lines), crossover from phase lag to phase advance is observed upon transition from Jsc to Voc, similar to that shown for SubPc/C60 in FIG. 11a. If kPPr is increased to 1 ns−1 (blue dashed lines), the modulus decreases at all values of bias, becoming small as Va→Voc, where a small phase advance is still predicted.



FIG. 12b shows the change in IMPS signal that results from a field-dependent kPPr as given by Eq. (24), where the zero-field recombination rates of kPPr0=(100 ns)−1 (solid black lines) and kPPr0=1 ns−1 (dashed blue lines) have been maintained. Similar to FIG. 12a, black curves show a crossover in the IMPS response to a phase advance near open-circuit. Whereas the modulus of the blue dashed curves is more substantial at short circuit, it remains small in the vicinity of Voc, similar to FIG. 12a.


Solid lines in the inset of FIG. 12b show the corresponding calculated PP densities, which are proportional to (kPPd+kPPr)−1. The exciplex luminescence data from FIG. 10a are reproduced and overlayed for comparison, showing agreement for V≤Voc. Similar calculations for the case of constant kPPr in FIG. 12a (not shown) fit poorly to the data over the entire range of bias.


According to Morteani et. al., excitons from the donor and acceptor bulk diffuse to the heterojunction where they undergo charge transfer and thermalize within the distance of a few molecular layers to form Coulombically bound geminate polaron pairs. These may dissociate in the electric field at the interface, or relax to neighboring donor and acceptor molecules, at which point they may form an exciplex or recombine to the ground state. Since polaron pairs are precursors to exciplex formation, and exciplex decay is expected to be largely field independent, the exciplex emission intensity provides a means for monitoring changes in the PP population at the HJ, which depends sensitively on the local electric field.


In this context, previous reports have observed exciplex emission for various Pc/C60 molecular dyads that is red-shifted and broadened from the fullerene fluorescence peak at λem=715 nm. The intensity trend of this red-shifted peak with increasing C60 thickness is directly proportional to the number of excitons that diffuse to the SubPc/C60 interface, as shown in FIG. 1b. In contrast, fluorescence from the C60 bulk is predicted to grow monotonically with C60 thickness. Hence, SubPc/C60 and CuPc/C60 interfacial exciplexes constitute a significant component of the broad emission bands at λem˜740-750 nm in their respective PL spectra of FIG. 1a. The exciplex is far more prominent in the PL spectrum of SubPc/C60 than in its EL spectrum. This may be the case because the EL spectrum is measured at high current density (J=450 mA/cm2), where hole injection beyond the heterojunction, leading to recombination in the C60 bulk dominates over that at the interface.


The weak bias dependence of SubPc fluorescence at λ=615 nm in FIG. 10a is consistent with the strongly bound bulk exciton, in contrast to the exciplex emission at λ=730 nm (c.f. FIG. 9a), which varies significantly. As Va→Voc the magnitude of the internal field in the device decreases, reducing dissociation, and thereby leading to increased emission. The exciplex bias dependence is more pronounced than that of the SubPc bulk exciton due to the high field sensitivity of the PP precursor state.


According to FIG. 10, the SubPc/C60 exciplex emission, and hence the PP density at this HJ, peaks at Voc. In contrast, the CuPc/C60 exciplex trend indicates that the geminate PP density decreases long before reaching open-circuit. Since the PP density follows ζ˜(kPPd+kPPr)−1, this results from an increase in kPPr, since kPPd decreases in moving from reverse bias to open-circuit due to the decrease in internal field. Therefore, SubPc/C60 PPs are dominated by the field dependence of dissociation (kPPd), whereas those of CuPc/C60 are dominated by the field dependence of recombination (kPPr).


The IMPS spectra of FIGS. 11a and 11b are also consistent with this interpretation. Both CuPc/C60 and SubPc/C60 show external quantum efficiencies of approximately 20% at short circuit, indicating that PP dissociation yields there are high (kPPd>kPPr). Near open-circuit, the SubPc/C60 IMPS spectrum shows a phase advance, indicating a significant increase in the interfacial free carrier concentration according to the crossover condition derived above. The non-vanishing IMPS modulus also indicates that ηPPd remains substantial (c.f. Eq. (23)), which in turn implies that kPPd is comparable to kPPr.


In contrast, the phase lag and steady decrease in modulus of the CuPc/C60 spectrum from short-circuit to open-circuit indicates that ηPPd becomes very small over this same interval. Hence, the IMPS spectra show that the ratio, kPPd/kPPr, decreases more strongly for CuPc/C60 between short- and open-circuit conditions than it does for SubPc/C60, consistent with the interpretation of the exciplex luminescence trends.


This difference is evident in the simulated IMPS spectra of FIG. 12a, calculated according to Eq. (23) for the two cases of kPPr=(100 ns)−1 and kPPr=1 ns−1. The smaller recombination rate allows for a high dissociation yield at short-circuit, and leads to an IMPS crossover similar to that observed for SubPc/C60. The fast recombination rate kPPr=1 ns−1 leads to a vanishing IMPS modulus at Voc, as seen for CuPc/C60, but also causes a significant decrease in the short-circuit quantum yield, in contrast to that observed for this device. Thus, a constant kPPr cannot reconcile the large short-circuit quantum yield and vanishing open-circuit IMPS modulus observed for CuPc/C60.


Therefore, this evidence suggests that the recombination rate changes with bias, and thus depends on electric field. Such field-dependence is often observed for donor-acceptor dyads in solution, where recombination of the charge-transfer excited state depends strongly on solvent polarity, and hence on the local electric field that is induced.


Equation (24) suggests a simple form for the field dependence of kPPr where the change in polaron pair potential with field modifies the activation energy barrier required for the transition, as shown schematically in FIG. 13. In reverse bias, FI<0 lowers the potential of the PP, decreasing the recombination rate, whereas in forward bias, FI>0 increases it. The potential shift typically ranges from −0.1 to 0.1 eV, and hence can significantly change the recombination rate at room temperature, especially approaching open-circuit where FI changes rapidly with applied voltage.


Marcus theory provides a more rigorous basis for the field-dependence of the recombination rate, where the driving force for the transition, ΔG, is modified by the change in PP potential. In the normal region, the driving force increases with FI and the field dependence of kPPr is similar to that given by Eq. (24). In the inverted region, however, increasing FI slows the recombination rate and vice versa. Since the recombination pathways and their associated free energies are not well established for the two HJs studied here, the Marcus formulation involves too many unknowns, and hence Eq. (24) was used as a reasonable approximation.


Simulating the IMPS spectra and PP density in FIG. 4b with the field-dependent kPPr as given by Eq. (24) results in trends that provide a qualitative match to the data. Using kPPr0=(100 ns)−1 reproduced the observed IMPS spectrum and, in the inset, predicted a PP density that follows the SubPc/C60 exciplex emission trend for V<Voc. Alternatively, when kPPr0=1 ns−1, the calculated PP density starts to decline under reverse bias, in reasonable agreement with the CuPc/C60 exciplex trend, where the emission floor observed for Va>0.5 V probably reflects the base fullerene fluorescence component, and not a leveling off of the exciplex emission. The field-dependent kPPr also expands the modulus of the simulated IMPS spectra at Jsc while maintaining its vanishing modulus at open circuit, in better agreement with the CuPc/C60 data of FIG. 11b.


The model for the PP density breaks down beyond Voc as the predicted density decreases more rapidly than the data in the inset of FIG. 12b. In this region, kPPr>kPPd; hence deviations from the approximation in Eq. (24) result in larger errors in the predicted PP density than when kPPr<kPPd at V<Voc. Additionally, other physical mechanisms may play a role. For example, carriers may undergo drift during thermalization, leading to a distribution of initial PP separations that is field dependent. Alternatively, the rate of exciton dissociation at the interface can also be field-dependent, where FI<0 (i.e., opposite to the direction aiding dissociation) decreases the free energy that drives the initial charge transfer step.


The variation observed in kPPr has consequences for solar cell performance. As discussed above, the decrease in open-circuit voltage from its theoretical Vocmax, depends on ln(kPPr/kPPd). High recombination rate is consistent with the IMPS and exciplex observations for CuPc/C60. Indeed, the approximately 100-fold difference in recombination rate that distinguishes the trends for SubPc/C60 and CuPc/C60 translates to a ˜0.25 V decrease in Voc at room temperature (see Eq. (17), above) in these devices. Hence, the fast recombination rate largely explains why CuPc/C60 cells operate at approximately 0.3 V below their Vocmax, whereas SubPc/C60 cells reach their maximum under standard operating conditions.


There are several possible routes to reduce kPPr and thus improve cell efficiency. For example, the electronic coupling upon which kPPr depends can be decreased by adding steric bulk to D and A molecules at the interface or by disrupting their ordered stacking as suggested previously. Only interfacial molecules within a distance ˜a0 of the HJ participate in the recombination process. Hence, the bulk of the D and A layers need not be altered from compositions or stacking habits that promote charge transport and exciton diffusion


Investigation of geminate polaron pair recombination in archetype CuPc/C60 and SubPc/C60 organic heterojunction solar cells using intensity modulated photocurrent spectroscopy, and measuring the emission from interfacial exciplex states has been provided. Also, provided is the nature of polaron pair recombination in both material systems, and validation of the rate-based model of PP dynamics used to derive the ideal diode equation for organic HJs. The PP density in SubPc/C60 peaks near open circuit, and is dominated by dissociation kinetics between V=0 and Voc, whereas that of CuPc/C60 decreases due to increasing recombination over this interval. An expression for the recombination rate that explains this behavior and qualitatively describes both the IMPS and exciplex luminescence trends observed for each cell are discussed above.


As provided herein, polaron pair recombination dynamics significantly impact the efficiency and open circuit voltages of organic solar cells. While this behavior for two small-molecule based cells has been provided, the techniques employed here are general, and may be applied to non-planar, bulk heterojunction devices based on, for example, small molecules or polymers.


An analytical model that describes the current-voltage characteristics of organic heterojunctions is provided herein. The model was based on polaron pair generation and recombination processes at the donor-acceptor interface, and leads to an ideal diode equation that is generally applicable to both polymer and small molecular weight, organic planar and bulk heterojunctions. The dark current is shown to be directly proportional to the polaron pair recombination rate, and the open-circuit voltage is ultimately limited by the difference between the donor-acceptor interface energy gap and the polaron pair binding energy, as observed previously.


It was proposed that recombination at traps within the disorder-induced density of states tail in each material at the D-A heterointerface results in two, temperature-dependent ideality factors, and a commonly observed double exponential J-V characteristic. Using an exponential trap distribution, the trends observed for the dark current, Voc, and Jsc as a function of temperature and intensity have been reproduced herein.


Thus, the aforementioned disclosure provides, in part, the following interpretation of the current vs. voltage characteristics: (1) The reverse bias slope of the dark current reflects the increasing dissociation probability of thermally generated PPs at the HJ; (2) in forward bias, there are two exponential regions, each with its own ideality factor stemming from recombination with trapped/immobile carriers at each side of the HJ; (3) the ideality factors are interpreted as the extent to which trap-limited recombination (one free and one immobile carrier) dominates over free carrier recombination (both carriers mobile); (4) the ratio of the two exponential current components under forward bias reflects the relative contribution of free donor holes recombining with trapped acceptor electrons vs. the reverse process to the total current density; (5) large ideality factors (>2) are due to a broad trap DOS distribution and/or large asymmetries in the potential dropped across the donor and acceptor layers; and (6) both ideality factors increase with decreasing temperature as more carriers freeze into the DOS tail and trap-limited recombination increases.


Specific examples of the invention are illustrated and/or described herein. However, it will be appreciated that modifications and variations of the invention are covered by the above teachings and within the purview of the appended claims without departing from the spirit and scope of the invention.

Claims
  • 1. An organic photosensitive optoelectronic device comprising: two electrodes in superposed relation;at least one donor material, andat least one acceptor material,said at least one donor material and at least one acceptor material forming at least one photoactive region between the two electrodes,wherein said at least one donor material and at least one acceptor material each comprises at least one of a single crystalline, nanocrystalline, and polycrystalline material in the bulk of the respective donor and acceptor materials,wherein the at least one donor material and the at least one acceptor material form a steric bulk region at the interface of the donor and acceptor materials, andwherein the steric bulk region comprises both molecules of the donor material and the acceptor material that have increased steric bulk compared to molecules of the donor material and the acceptor material outside of the steric bulk region.
  • 2. The device of claim 1, wherein said steric bulk region comprises less than ten monolayers of said donor material and/or said acceptor material.
  • 3. The device of claim 1, wherein the photo-active region forms at least one of a mixed heterojunction, bulk heterojunction, nanocrystalline-bulk heterojunction, and hybrid planar-mixed heterojunction.
  • 4. The device of claim 1, wherein the donor material is chosen from subphthalocyanine (SubPc), copper pthalocyanine (CuPc), chloroaluminium phthalocyanine (CIAIPc), tin phthalocyanine (SnPc), pentacene, tetracene, diindenoperylene (DIP), and squaraine (SQ).
  • 5. The device of claim 1, wherein the acceptor is chosen from C60, C70, 3,4,9,10 -perylenetetracarboxylicbis-benzimidazole (PTCBI), PCBM, PC70BM and hexadecafluorophthalocyanine (F16CuPc).
  • 6. The device of claim 1, wherein at least one electrode comprises transparent conducting oxides or transparent conducting polymers.
  • 7. The device according to claim 6, wherein the conducting oxides are chosen from indium tin oxide (ITO), tin oxide (TO), gallium indium tin oxide (GITO), zinc oxide (ZO), and zinc indium tin oxide (ZITO), and the transparent conductive polymers comprise polyanaline (PANI).
  • 8. The device according to claim 1, wherein at least one electrode comprises a metal substitute, a non-metallic material or a metallic material chosen from Ag, Au, Ti, Sn, and Al.
  • 9. The device of claim 1, further comprising at least one exciton blocking layer.
  • 10. The device of claim 9, wherein said exciton blocking layer is chosen from bathocuproine (BCP), bathophenanthroline (BPhen), 3,4,9,10-perylenetetracarboxylicbis-benzimidazole (PTCBI), 1,3,5-tris(N-phenylbenzimidazol-2-yl)benzene (TPBi), tris(acetylacetonato) ruthenium(III) (Ru(acaca)3), and aluminum(111)phenolate (Alq2 OPH).
  • 11. The device of claim 1, wherein the device is an organic photodetector.
  • 12. The device of claim 1, wherein the device is an organic solar cell.
  • 13. A method of making an organic photosensitive optoelectronic device, said method comprising depositing onto a substrate at least one donor material, at least one acceptor material, and two electrodes, wherein said at least one donor material and at least one acceptor material are deposited to form at least one photoactive region between the two electrodes,wherein said at least one donor material and at least one acceptor material each comprises at least one of a single crystalline, nanocrystalline, and polycrystalline material in the bulk of the respective donor and acceptor materials,wherein the at least one donor material and the at least one acceptor material form a steric bulk region at the interface of the donor and acceptor materials, andwherein the steric bulk region comprises both molecules of the donor material and the acceptor material that have increased steric bulk compared to molecules of the donor material and the acceptor material outside of the steric bulk region.
  • 14. The method of claim 13, wherein said steric bulk region comprises less than ten monolayers each of said donor material and/or said acceptor material.
  • 15. The method of claim 13, wherein said depositing occurs in a vacuum.
  • 16. The method of claim 13, wherein the photo-active region forms at least one of a mixed heterojunction, bulk heterojunction, nanocrystalline-bulk heterojunction, and hybrid planar-mixed heterojunction.
  • 17. The method of claim 13, wherein the donor layer is chosen from subphthalocyanine (SubPc), copper pthalocyanine (CuPc), chloroaluminium phthalocyanine (CIAIPc), tin phthalocyanine (SnPc), pentacene, tetracene, diindenoperylene (DIP), and squaraine (SQ).
  • 18. The method of claim 13, wherein the acceptor is chosen from C6o, C7o, 3,4,9,10-perylenetetracarboxylicbis-benzimidazole (PTCBI), PCBM, PC70BM and hexadecafluorophthalocyanine (F16CuPc).
  • 19. The method of claim 13, wherein the at least one electrode comprises a metal substitute, a non-metallic material or a metallic material chosen from Ag, Au, and Al.
  • 20. The method of claim 13, further comprising depositing at least one exciton blocking layer between the substrate and the electrode, wherein said exciton blocking layer is chosen from bathocuproine (BCP), bathophenanthroline (BPhen), 3,4,9,10-perylenetetracarboxylicbis-benzimidazole (PTCBI), 1,3,5-tris(N-phenylbenzimidazol-2-yl) benzene (TPBi), tris(acetylacetonato) ruthenium(III) (Ru(acaca)3), and aluminum(lll)phenolate (Alq2 OPH).
  • 21. The method of claim 13, wherein at least one electrode comprises transparent conducting oxides or transparent conducting polymers.
  • 22. The method of claim 21, wherein the substrate comprises a conducting oxide chosen from indium tin oxide (ITO), tin oxide (TO), gallium indium tin oxide (GITO), zinc oxide (ZO), and zinc indium tin oxide (ZITO), and the transparent conductive polymers comprises polyanaline (PANT).
CROSS-REFERENCE TO RELATED APPLICATION

This application claims the benefit of U.S. Provisional Application No. 61/331,798, filed May 5, 2010, which is incorporated herein by reference in its entirety. The subject matter of the present disclosure was made by, on behalf of, and/or in connection with one or more of the following parties to a joint university-corporation research agreement: University of Michigan and Global Photonic Energy Corporation. The agreement was in effect on and before the date the subject matter of the present disclosure was prepared, and was made as a result of activities undertaken within the scope of the agreement.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH

This invention was made with U.S. Government support under Contract No. DE-SC0000957 awarded by the Department of Energy and FA9550-07-1-0364 awarded by the Air Force Office of Scientific Research. The government has certain rights to this invention.

US Referenced Citations (6)
Number Name Date Kind
20050224113 Xue et al. Oct 2005 A1
20050227390 Shtein et al. Oct 2005 A1
20060060239 Peumans et al. Mar 2006 A1
20070215868 Forrest Sep 2007 A1
20080116536 Forrest et al. May 2008 A1
20110132440 Sivarajan Jun 2011 A1
Foreign Referenced Citations (2)
Number Date Country
WO 2009029548 Mar 2009 WO
WO 2010036963 Apr 2010 WO
Non-Patent Literature Citations (3)
Entry
Yang et al., “Nanoscale Morphology of High-Performance Polymer Solar Cells”, Nano Letters, 2005.
International Search Report from PCT/US2011/035017.
Perez, M. Dolores et al., “Molecular and Morphological Influences on the Open Circuit Voltages of Organic Photovoltaic Devices,” J. Am. Chem. Soc. 2009, 131, pp. 9281-9286.
Related Publications (1)
Number Date Country
20110297234 A1 Dec 2011 US
Provisional Applications (1)
Number Date Country
61331798 May 2010 US