R Signerski - Modeling of photovoltaic properties of bilayer organic systems - страница 1
ВІСНИК ЛЬВІВ. УН-ТУ
Серія фізична. 2007. Вип.40. С.159-167
VISNYKLVIV UNIV. Ser.Physic. 2007. N40. P. 159-167
PACS number(s): 73.50.Pz, 72.80.Le
MODELING OF PHOTOVOLTAIC PROPERTIES OF BILAYER ORGANIC SYSTEMS
Department of Physics of Electronic Phenomena
Gdansk University of Technology G. Narutowicza 11/12, 80-952 Gdansk, Poland e-mail: firstname.lastname@example.org
The work presents a theoretical description of current-voltage dependence and photovoltaic effect of a system formed from two organic layers and electrodes. The model is based on the assumption that in one of organic layers only an electron current can flow, while in the other only a hole current. Taking into account a constant-field approximation simple expressions describing the current-voltage relationship and short-circuit current have been obtained. The influence of space charge in organic layers related to trapped charge carriers is analysed on the basis of numerical calculations. The obtained results enable to estimate the role of electrodes, charge carrier traps and generation-recombination processes in photovoltaic properties of an organic bilayer system.
Key words: organic photovoltaic cells, organic heterojunctions.
Photovoltaic properties of organic heterojunctions have been widely studied due to the prospect of obtaining stable, efficient and low cost solar cells. The power-conversion efficiency in dye-sensitized solar cells has already overcome 10% and in the polymer bulk heterojunction and small molecular weight multilayer solar cells 5% [1-3]. The further progress in organic photovoltaic technology will be determined among others by the progress in the theoretical description of photovoltaic phenomenon. However, due to a wide variety of organic systems [1, 2, 4, 5] it is not possible to elaborate the universal model of this phenomenon.
There are many works, in which authors focused their attention on different processes determining photovoltaic parameters. As an example, in the works [6-8] a short-circuit current generated as a result of exciton dissociation at the interface of two organic layers is calculated, while in the works [9-11] numerical calculations of current-voltage relationships of either bilayer [9, 10] or bulk photovoltaic  devices are presented. Moreover, electric equivalent circuit of organic solar cells is discussed in the works [12-14] and factors limiting the power conversion efficiency of molecular photovoltaic devices are analysed in [15, 16].
This work deals with modeling of current-voltage characteristics of bilayer organic systems. In particular, a concentration of charge carriers at electrodes, generation-
© Signerski R., 2007
recombination processes in junction and trapping of charge carrier in bulk of layers have been taken into account.
Fig. 1. An energy-level diagram of bilayer organic system
Fig. 1 shows an energy-level diagram of the analysed system. The system consists of two organic layers and electrodes. One of organic layers is an electron-transporting layer (ETL) of dn-thickness, while the other is a hole-transporting layer (HTL) of dp-thickness. Bands transporting electrons and holes are schemed as straight lines. The band-bending or potential jumps resulting from dipole layers, which can appear at electrodes or organic-organic interface, are not respected. Free charge carriers are generated only at the ETL-HTL interface either by thermal generation or by exciton dissociation. Thickness of the region, where the charge carrier generation occurs, is small, compared to the thickness of organic layers.
The main assumption of this model supposes that in ETL only an electron current flows, while in HTL only a hole current flows [9, 10].
Let us assume that in each organic layer the electric field is uniform. This assumption is true at low current intensity, low concentration of trapping states and strong built-in electric field. The effect of non-uniform of electric field occurring in thin near-electrode regions is assumed to be negligible.
Now the ETL will be considered as an example. Current density flowing through this layer is described by the expression:
• ( ) d9n (x) + , T dn(x)
j = -e ^nn(x)—--+\ьпкТ——, (1)
where \in - electron mobility, n(x) - concentration of electrons at x, <pn(x) - electric potential at x.
Solving this equation with respect to n(x) and integrating the obtained experssion over <0,dn> as well as appllying the constant-field approximation (e.g.), the following curent-voltage relationship can be derived:
. e nj . \ j J njexp
where Un- voltage across ETL, no - concentration of electrons at the electrode, nj , nj- concentrations of electrons in heterojunction either in thermal equilibrium (njh^
in the presence of current or illumination (nj) .
Furthermore, it will be assumed that the concentration no is independent of voltage and illumination (ohmic contact), and the concentration of electrons in heterojunction is determined by the monomolecular-kinetic equations:
for unilluminated system, (3)
Ghn - Rnnj = 0
Gex + Gtl
for illuminated system.
In the above equations Gth n is a rate of thermal generation, Gex is a rate of photogeneration, Rn is a rate constant of recombination and a is a thickness of generation-recombination region. It is usually assumed that this region is limited to two monomolecular layers located on both sides of a heterojunction interface.
The rate of thermal generation is determined by an energy gap between a valence band of HTL and a conduction band of ETL. The rate constant of photogeneration of charge carriers, Gex, is determined by a flux of exciton reaching the heterojunction interface and by quantum efficiency of exciton dissociation. Usually, Gex is proportional to a light intensity Io, i.e. Gex~Io [7, 8]. The monomolecular recombination assumed in Eqs. (3) and (4), and characterized by the rate constant Rn, is noticeable in the presence of adequately high concentration of trapping-recombination centers in the junction.
If now Eqs. (3) and (4) are taking into account then the current-voltage relationship given by Eq. (2) takes the form:
- +1 - exp
e(Un + Ubin )
Ub,n =— In— = — b-
is the built-in potential for ETL. In particular, Eq. (5) leads to the following formula of: • open-circuit current (when j=0):
short-circuit current (when Un=0):
e V n Ge-.
1 - exp
VJ Since Gex~Io it can be concluded that both current densities, namely j (Eq. (5)) and jsc, (Eq. (8)) are proportional to light intensity Io, while the voltage Uocn (Eq.(7)) exhibits logarithmic dependence on Io. Such relationships are often noticed in organic heterojunctions [1,4,10,18]. Analysis of Eq. (5) indicates that for high positive values of Un the current yields:
j = ea Gn + ) (9) and it is the sum of saturation-current density of the junction:
jsat = eaGthn (10)
and the maximal value of short-circuit-current density:
jScmax = eaGex .
From Eq. (8) it can be noticed that short-circuit current reaches the maximal value when the built-in potential is enough high:
kT 1 Jh^
40 20 0 -20 -40 -60 -80 -100 -120
1 - G =3*1020cm-3s-1
2 - G =1021cm-3s-1
Fig. 2. Photocurrent-voltage relationships for the constant-field approximation
The analysis presented above deals with the current j flowing through ETL. Obviously, analogous relations can be ascribed to HTL. If the flow of the same current j through HTL requires Up then the voltage Ubetween electrodes yields:
U = Un + Up (13)
and j(U) is the current-voltage relationship of the organic bilayer system. As an example, two current-voltage curves calculated according to Eqs. (5) and (13) for different values
of the photogeneration rate constant are presented in Fig. 2. For simplification it was assumed that Un=Up=U/2. In calculation the following values were taken: 7>290K, Jun=^p=10\m2(Vs)'1, d„=d=W5cm, a=2-10-7cm, Gthn=Gthp = 1015cm-V, Rn=Rp=108s-1, n0=p0=1012cm-3. It is seen that for positive voltage the current density is constant and is governed by Eq.(9). Within the photovoltaic range (the second quarter) the fill factor:
FF = V J '™* (14)
can be estimated. It yields 0.53. In the case of organic hererojunctions it is quite high value. In the third quarter the linear dependence of current densisty versus voltage is observed. The analysis of Eq. (5) indicates that the slop of the curves j(Un=U/2) within
this range of voltage yields e ^n"0 and it is determined by mobility of charge carriers,
thickness of organic layers and concentration of charge carriers at electrodes.
Let us analyse now the effect of space charge on the shape of current-voltage curves of the organic bilayer system. If the space-charge concentration is enough high then assumption of uniform field is not true any more. Similarly, as in the previous chapter the attention will be focused on ETL. The space dependence of electric field En(x) in this layer is governed by the Poisson equation:
dEn (x) = d 2cpn (x) = e ( n( x) + nt (x))
in which є0 denotes vacuum permittivity, &rn is dielectric constant and nt(x) is the concentration of trapped charge carriers.
Let us assume that charge carriers are trapped by a discrete energy level of N - concentration and Et - depth. The equation describing processes of trapping and thermal releasing of charge carriers takes the form:
Ynn(x) (N - nt (x)) - v exp ^-kT j nt (x) = 0, (16)
where yn is a rate constant of trapping and v is a thermal collision factor. The set of Eqs. (1), (15), (16), together with the boundary conditions of the following forms: n(x=0)=no, (3) and (4), are solved numerically with respect to n(x) [19,20]. Afterwards the voltage across ETL is calculated:
j d?dx kT, n,
Un = — J — + — to-h . (17)
eц 0 n e nj
The same procedure is applied for HTL and as a result of this the voltage Up is obtained. The potential difference between electrodes of the system is equal to the sum of Up and Un, Eq. (13).
^ 0 О
*|_ -20 -40 -60 -80
-4-3-2-10 1 2 3 4
Fig. 3. Photocurrent-voltage relationships calculated numerically at various densities of charge carrier traps: 1 - N=1015cm-3, 2 - N=1017cm-3, 3 - N=1018cm-3
Current-voltage curves calculated according to the above procedure for the sample under illumination at three different concentrations of traps are presented in Fig.3. Similarly as in the chapter 2.1 it was assumed that Un=Up=U/2. Values of parameters are the same as those taken for Fig.2 and additionally sr=3,5, no=po=1012cm-3, i^1012s-1, E=0.5eV, 7n=10-9cm3s-1.
It is easy to notice that the curve 1 is the same as the curve 2 in Fig.2. It means that charge carrier traps with the concentration not higher than 1015cm-3 (and the depth 0.5eV) do not limit the current flowing through the system. For higher values of the trap concentration the short-circuit current (jsc) and the fill factor (FF) are lower than in the approximation of constant field. The curve 3 departs quite significantly from the curves 1 and 2. Such a shape of current-voltage curves both within and in the vicinity of the photovoltaic range is often recognized as the result of an effect of electric field on the efficiency of charge carrier photogeneration [10,11]. However, in this work such an effect was not respected.
Analysis of current-voltage relationships as well as photovoltaic properties of bilayer organic system presented in this work was based on two main assumptions, namely that exciton dissociation occurs in very thin region near the organic-organic interface and that electrons are transported only by one of organic layers and holes by the other. Simple analytical expressions can be obtained when approximations of constant electric field within organic layers and of monomolecular recombination at organic-organic interface are applied. On the basis of them the role of such parameters as the rates of thermal- and photo-generation of charge carriers, the rate constant of recombination, charge carrier concentrations at electrodes and thickness of organic layers in electric properties of the system can be estimated. In particular, the obtained relations predict the linear dependence of photocurrent and short-circuit current versus light intensity. However, if bimolecular recombination at organic-organic interface is assumed [11, 18]:
• 1 ."
' 3/x '/;'
/ / / ■ ■ ■ ■
Gth + Gex -anjPj = J- (18) ea
with a denoting the rate constant of bimolecular recombination and nj, pj denoting concentrations of electrons and holes in heterojunction, then applying the same procedure as in the chapter 2,1 it can be easy proved that the dependence of short-circuit current versus light intensity (jsc(I0)) is a square root function but within the range of low values of light intensity this dependence can be approximated by a linear function.
If the approximation of constant-electric field is not fulfilled then the numerical calculations have to be performed. Results presented in the chapter 2,2 show the effect of charge carrier traps on the shape of current-voltage curves within the photovoltaic range. The trap concentration enough high results in the same effect as field dependence of photogeneration of charge carriers, namely leads to the decrease of values of short-circuit current jsc and fill factor FF. In the analysis of the influence of charge carrier traps it was assumed only a discrete trapping level of the depth Et, however the same procedure can be applied also for a different trap distribution, e.g. the Gaussian one . Spectral dependences of photocurrent and other quantities characterizing the photovoltaic effect were not considered in this work. These relationships are additionally determined by the relation between the rate of charge carrier photogeneration and the wavelength, i.e. Gex(k), hence they depend on absorption coefficient, length of exciton diffusion, quantum efficiency of exciton dissociation and thickness of layers.
The theoretical description of photovoltaic effect presented in this work was based first of all on parameters, which values can be modified experimentally when designing a photovoltaic system.
This work has been supported by means for a science granted for years 2006-2009 by Polish Ministry of Education and Science under Program No. 3T11B06530.
1. Brabec C.J., Dyakonov V., Parisi J., Sariciftci N. S. (Eds) Organic Photovoltaics. Berlin Heidenberg: Springer-Verlag, 2003.
2. Hoppe H., Sariciftci N. S.Organic solar cells: An overview // J. Mater. Res. 2004. Vol. 83. P. 1924-1945.
3. Xue J., Uchida S., RandB. P., Forrest S. R. Asymetric tandem organic photovoltaic cells with hybrid planer-mixed molecular heterojunctions // Appl. Phys. Lett. 2004. Vol. 85. N 23. P. 5757-5759.
4. Peumans P., Yakimov A., Forrest S.R. Small molecular weight organic thin-film photodetectors and solar cells // J. Appl. Phys. 2003. Vol. 93. N 7. P. 3693-3723.
5. Waldauf Ch., Scharber M. C., Schilinsky P., Huch J. A., Brabec Ch. J. Physics of organic bulk heterojunction devices for photovoltaic applications // J. Appl. Phys. 2006. Vol. 99. Р. 1045031.
6. Wagner J., Fritz T., Bottcher H. Computer Modelling of Organic Thin Film Solar Cells // Phys. stat. sol. (a). 1993. Vol. 136. P. 423-432.
7. Stubinger T., Brutting J. Exciton diffusion and optical interference in organic
donor-acceptor photovoltaic cells // J. Appl. Phys. 2001. Vol. 90. N 7. P. 36323641.
8. Petteron L. A. A., Roman L. S., Inganas O. Modeling photocurrent action spectra of photovoltaic devices based on organic thin films // J. Appl. Phys. 1999. Vol. 86. N 1. P. 487-496.
9. Gregg B. A., Hanna M. C. Comparing organic to inorganic photovoltaic cells: Theory, experiment, and simulation // J. Appl. Phys. 2003. Vol. 93. N 6. P. 36053614.
10. Barker J. A., Ramsdale C. M., Greenham N. C. Modeling the current-voltage characteristics of bilayer polymer photovoltaic devices // Phys. Rev. B. 2003.
Vol. 67. P. 0752051.
11. Koster L. J. A., Smits E. C. P., Mihailetchi V. D., Blom P. W. M. Device model for the operation of polymer/fullerene bulk heterojunction cells // Phys. Rev. B. 2005.
Vol. 72. P. 085201.
12. Rostalski J., Meissner D. Photocurrent spectroscopy for the investigation of charge carrier generation and transport mechanisms in organic p/n junction solar cells // Solar Energy Mat. & Solar Cells. 2000. Vol. 63. P. 37-47.
13. Yoo S., Domercq B., Kippelen B. Intensity-dependent equivalent circuit parameters of organic solar cells // J. Appl. Phys. 2005. Vol. 97. P. 103707061.
14. Mazhari B. An improved solar cell circuit model for organic solar cells // Solar Energy Mat. & Solar Cells. 2006. Vol. 90. P. 1021-1033.
15. Nelson J., Kirkpatrick J., Ravirajan P. Factors limiting the efficieny of molecular
photovoltaic devices // Phys. Rev. B. 2004. Vol. 69. P. 033371.
16. Sun S., Fan Z., Wang Y., Haliburton J. Organic solar cell optimizations // J. Mater. Sci. 2006. Vol. 40. P. 1429-1443.