Design of high efficiency organic solar cell with light ... - OSA Publishing

11 downloads 3020 Views 2MB Size Report
*[email protected]. Abstract: We have designed a high efficiency organic solar cell with light trapping structure on transference cylindrical substrate.
Design of high efficiency organic solar cell with light trapping L. Song and A. Uddin* School of Photovoltaic and Renewable Energy Engineering, the University of New South Wales, Sydney, Australia * [email protected]

Abstract: We have designed a high efficiency organic solar cell with light trapping structure on transference cylindrical substrate. An electrical and optical simulation of the light trapping structure has been performed on the basis of finite element and transfer matrix formalism methods. Absorption spectrum, internal quantum efficiency, external quantum efficiency, maximum power output and efficiency of the organic solar cell are simulated and presented in terms of three variables: the height, diameter of the glass substrate and the thickness of the organic active layer. The efficiency of the proposed organic solar cell with light trapping structure is enhanced by a factor of 2 than the similar structure on the flat plain glass substrate. The optimum organic active layer thickness to achieve the highest efficiency is shifted from 65 to 20 nm. Finally, we have investigated the effect of light incident angle on the performance of the proposed cell structure. © 2012 Optical Society of America OCIS codes: (040.5350) Photovoltaic; (220.0220) Optical design and fabrication.

References and links 1.

2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12.

13. 14. 15.

M. Campoy-Quiles, T. Ferenczi, T. Agostinelli, P. G. Etchegoin, Y. Kim, T. D. Anthopoulos, P. N. Stavrinou, D. D. C. Bradley, and J. Nelson, “Morphology evolution via self-organization and lateral and vertical diffusion in polymer:fullerene solar cell blends,” Nat. Mater. 7(2), 158–164 (2008). Konarka, “Single junction solar cell by Konarka with an efficiency of 8.3% on an area of 1 cm2”. G. Dennler, M. C. Scharber, and C. J. Brabec, “Polymer-fullerene bulk-heterojunction solar cells,” Adv. Mater. (Deerfield Beach Fla.) 21(13), 1323–1338 (2009). W. Ma, C. Yang, X. Gong, K. Lee, and A. J. Heeger, “Thermally stable, efficient polymer solar cells with nanoscale control of the interpenetrating network morphology,” Adv. Funct. Mater. 15(10), 1617–1622 (2005). M. S. Ryu, H. J. Cha, and J. Jang, “Effects of thermal annealing of polymer:fullerene photovoltaic solar cells for high efficiency,” Curr. Appl. Phys. 10(2), S206–S209 (2010). A. J. Moulé and K. Meerholz, “Controlling morphology in polymer–fullerene mixtures,” Adv. Mater. (Deerfield Beach Fla.) 20(2), 240–245 (2008). T. Kirchartz, K. Taretto, and U. Rau, “Efficiency limits of organic bulk heterojunction solar cells,” J. Phys. Chem. C 113(41), 17958–17966 (2009). A. C. Mayer, S. R. Scully, B. E. Hardin, M. W. Rowell, and M. D. McGehee, “Polymer-based solar cells,” Mater. Today 10(11), 28–33 (2007). C. Heine and R. H. Morf, “Submicrometer gratings for solar energy applications,” Appl. Opt. 34(14), 2476–2482 (1995). J. E. Cotter, “Optical intensity of light in layers of silicon with rear diffuse reflectors,” J. Appl. Phys. 84(1), 81– 98 (1998). Y. Yi, L. Zeng, C. Hong, J. Liu, N. Feng, X. Duan, L. C. Kimerling, and B. A. Alamariu, “Efficiency enhancement in Si solar cells by textured photonic crystal back reflector,” Appl. Phys. Lett. 89, 111111 (2006). H. Hoppe, M. Niggemann, C. Winder, J. Kraut, R. Hiesgen, A. Hinsch, D. Meissner, and N. S. Sariciftci, “Nanoscale morphology of conjugated polymer/fullerene-based bulk- heterojunction solar cells,” Adv. Funct. Mater. 14(10), 1005–1011 (2004). S. D. Zilio, K. Tvingstedt, O. Inganäs, and M. Tormen, “Fabrication of a light trapping system for organic solar cells,” Microelectron. Eng. 86(4-6), 1150–1154 (2009). H. Huang, Y. Li, M. Wang, W. Nie, W. Zhou, E. D. Peterson, J. Liu, G. Fang, and D. L. Carroll, “Photovoltaic– thermal solar energy collectors based on optical tubes,” Sol. Energy 85(3), 450–454 (2011). Y. Li, W. Nie, J. Liu, A. Partridge, and D. L. Carroll, “The optics of organic photovoltaics: fiber-based devices,” IEEE J. Sel. Top. Quantum Electron. 16(6), 1827–1837 (2010).

#164841 - $15.00 USD (C) 2012 OSA

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A606

16. Y. Li and E. D. Peterson, H. Huang, M. Wang, D. Xue, W. Nie, W. Zhou and D.L. Carroll, “Tube-based geometries for organic photovoltaics,” Appl. Phys. Lett. 96, 243503 (2010). 17. M. R. Lee, R. D. Eckert, K. Forberich, G. Dennler, C. J. Brabec, and R. A. Gaudiana, “Solar power wires based on organic photovoltaic materials,” Science 324(5924), 232–235 (2009). 18. J. Liu, M. A. G. Namboothiry, and D. L. Carroll, “Optical geometries for fiber-based organic photovoltaics,” Appl. Phys. Lett. 90(13), 133515 (2007). 19. L. A. A. Pettersson, L. S. Roman, and O. Inganäs, “Modeling photocurrent action spectra of photovoltaic devices based on organic thin films,” J. Appl. Phys. 86(1), 487–496 (1999).

1. Introduction Organic photovoltaic (OPV) has attracted much attention in recent years due to the rapid increasing of its power conversion efficiency (PCE) and the potential benefits of low cost production methods [1]. The highest reported PCE of an OPV device is ~8.0% [2], which is approaching to the required efficiency for commercialization [3]. At present, OPV devices with the highest performance are based on the bulk heterojunction (BHJ) type in which an electron-donor polymer material and an electron-acceptor fullerene material are blended to form an active photovoltaic layer. The most commonly studied polymer-fullerene system is based on a blend of donor poly (3-hexylthiophene-2,5-diyl) and acceptor [6,6]-phenyl-C61 butyric acid methyl ester(P3HT: PCBM), in which the highest reported PCE values are in the range of 4 to 6% [4–6]. The theoretical PCE limits of BHJ are over 23% [7]. The difference between the theory and experimental PCE can be explained by five reasons: namely the energy level misalignment at the donor/acceptor interfaces, insufficient light trapping, low exciton diffusion lengths, nonradiative recombination and low charge carrier mobilities. The OPV device structure needs a thin active layer due to the low exciton diffusion length (~10 nm) for the high efficiency device [8], but the thin active layer is insufficient to absorb most of the incident light. The light trapping structure is one of the potential ways to improve the PCE of BHJ by increasing the effective optical length within the active layer without the need of altering its physical thickness. The conventional light trapping structures for inorganic solar cells can’t be used on the OPV due to their larger size typically in a micrometer scale, compared to the thickness of OPV active layer in nanometer scale. Some light trapping structures for OPV, have been proposed such as submicron or micron structures, buried Nano-electrodes, multireflection structures and distributed Bragg reflectors [9–12]. However, light trapping structures in direct contact with thin active layer is most likely to introduce defects and contamination into the organic solar cell. Other light trapping structures like micro-lenses structure can avoid direct contact with the active layer, but it is only valid for the normally incident light [13]. Some researchers have also proposed the light trapping structure based on the tube-like geometries [14]. However, the bottom of the structure is closed in a hemispherical cap rather than a proposed back diffuse reflector. Besides, all the parameters associated with the light trapping structure, such as the height of the tube, the thickness of the OPV cell and etc, have not been optimized due to the lack of the detailed theoretical analysis. Finally, the effect of the different incident angles on the efficiency of OPV has not been fully discussed. Other geometries based on horizontal tube-like structure requires very complicated structure and it is still lack of fully theoretical analysis on the effect of all parameters [15–18]. The concept of the proposed structure of OPV device with light trapping on cylindrical glass substrate is shown in Fig. 1 and Fig. 2. The proposed OPV device structure has consisted a columnar glass substrate, indium tin oxide (ITO) layer, Poly (3,4ethylenedioxythiophene) -poly (styrenesulfonate) (PEDOT:PSS) layer, Copper phthalocyanine (CuPc) donor, thiophene-C60 acceptor, bathocuproine (BCP) layer, aluminum (Al) layer as electron collecting electrode (100 nm thick) and lambertian diffuse reflector. The incident light on the columnar glass substrate passes through the glass and is scattered in all directions in the upper hemisphere by the back diffuse reflector (96% lambertian reflector). The scattered lights bounce left and right through the active layer of the device. The

#164841 - $15.00 USD (C) 2012 OSA

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A607

unabsorbed photons in the first pass have additional opportunities to be absorbed in the subsequent passes. In this way, the effective optical length is dramatically increased.

Fig. 1. Schematic diagram of the proposed OPV device structure on cylindrical glass substrate. The device structure contains layers of Al, BCP, CuPc/PCBM active layer, PEDOT:PSS, ITO and cylindrical glass substrate.

Fig. 2. It is the cross sectional view of the proposed OPV structure on cylindrical glass substrate. The random pattern plain at the back of the glass substrate represents the diffuse reflector. When the light is incident on the diffuse reflector, it will be scattered into different directions.

#164841 - $15.00 USD (C) 2012 OSA

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A608

In order to compare the performance of the proposed OPV device structure with light trapping structure, a similar OPV device cell on flat plain glass substrate without any light trapping structure is also simulated and analyzed in the similar way. The OPV device without light trapping is shown in Fig. 3. It contains same structure of the proposed device such as an Al electrode, BCP layer, C60 acceptor, CuPc donor, PEDOT: PSS layer, ITO layer and plain glass substrate.

Fig. 3. Schematic diagram of the conventional OPV device structure on flat plain glass substrate, called control-Planar OPV device. The device structure contains layers of Al, BCP, CuPc/PCBM active layer, PEDOT:PSS, ITO and flat glass substrate with similar layer thicknesses of the proposed OPV device.

2. Theory The proposed OPV device structure described above involves an optical modeling of the light path within the light trapping structure and a comprehensive electrical modeling of the exciton generation and transfer, exciton dissociation and charge transport and collection. After the modeling of optical and electrical properties of the proposed device, we could derive the mathematical expressions for optical absorption curve, external quantum efficiency (EQE), internal quantum efficiency (IQE), short circuit current and power conversion efficiency. The needed assumptions for the model are: (1) the glass substrate and all other layers are considered as homogeneous and isotropic; (2) interfaces are flat and parallel compared to the wavelength of the light [19]; (3) exciton diffusion conforms to the diffusion equation; (4) excitons contribute to the photocurrent are dissociated at the D-A interfaces; (5) the surface of the organic solar cell acts as the exciton sink with large recombination velocity; (6) all the simulations are carried out in the standard test conditions: AM1.5 spectrum, 1000 W/m2 irradiance at 25 C cell temperature; (7) The open circuit voltage of this organic solar cell is 0.5 V and fill factor (FF) is 0.6. 2.1 Optical modeling 2.1.1 Energy dissipation and diffuse reflector It is considering a plane light wave incident from top on a glass substrate with an incident angle θ. Part of the incident light would be reflected by the air/glass interface and the reflection loss could be expressed by 1 sin 2 (θ − θ ) tan 2 (θ1 − θ 2 ) ] Rn = [ 2 1 2 + 2 sin (θ1 + θ 2 ) tan 2 (θ1 + θ 2 )

(1)

where θ1 is the angle of incidence, θ 2 is the angle of refraction, R is the reflection coefficient.

#164841 - $15.00 USD (C) 2012 OSA

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A609

When the incident light is normal to the air/glass interface, the reflection loss is reduced to:  n −n  R= 2 1   n2 + n1 

2

(2)

where, n1 and n2 are the refractive index of the air and glass substrate respectively. The light absorption by the glass substrate is considered to be negligible. Incident light scattered into different directions after reflection from the bottom of glass substrate. The scattering light follow the Lambert's cosine law. The differential radiant flux at any point could be described by the Lambertian distribution [10]:

d Φ L (θ , φ ) = I max cos(θ )sin(θ )dθ dφ

(3)

where θ and φ are the altitude angle defined with respect to the surface normal and azimuth angle and φL is the radiant flux contained within a differential solid angle of sin(θ )d θ d φ ,

I max is the peak radiant flux per steradian emitted out of the surface at normal angle. The total radiant flux from the diffuse reflector for a specific wavelength is equal to the incident radiant flux, which is expressed as: Φ total (λ ) = π I max = Wf (λ )(1 − R ) βπ r 2

(4)

where W is the incident radiant intensity, f (λ ) is the energy percentage of a particular wavelength in AM1.5 spectrum, β is the reflectance efficiency of the diffuse reflector, r is the radius of the glass substrate. Substituting Eq. (1), (2), (3) into Eq. (4), one can get the radiant flux contained within a differential solid angle for a particular wavelength as follow:

d Φ L (θ , φ , λ ) = Wf (λ )(1 − R) β r 2 cos(θ ) sin(θ )dθ dφ

(5)

where λ is the wavelength of the light. At this stage, the s and p polarized lights have the equal intensities, thus one can get the following expressions:

d Φ s _ polarized (θ , φ , λ ) = 0.5Wf (λ )(1 − R) β r 2 cos(θ )sin(θ )dθ dφ

(6)

d Φ p _ polarized (θ , φ , λ ) = 0.5Wf (λ )(1 − R)β r 2 cos(θ )sin(θ )dθ dφ

(7)

2.1.2 The distribution of exciton generation The modeling of the exciton generation distribution in principle is following the method described by Pettersson [19], however, this treatment is modified and expanded for nonnormal incidence case considering the effect of increasing polarization of light after each reflection. When the light with specific solid angle is incident on the multilayer, the propagation of light is generally described by the matrices. Since the tangential component of the electric field is continuous according to the theory of electromagnetic, the interface property could be expressed by: I jk =

#164841 - $15.00 USD (C) 2012 OSA

1 t jk

1   rjk

rjk   1

(8)

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A610

where I jk is the optical transfer matrix of the interface, t jk is the transmission coefficient of the interface, r jk is the reflection coefficient of the interface. The t jk and r jk could be expressed as follows:

t sjk =

2sin θ k cos θ j sin(θ j + θ k )

rjks = − t jkp = rjkp =

sin(θ j − θ k ) sin(θ j + θ k )

(9)

2 sin θ k cos θ j

sin(θ j + θ k ) cos(θ j − θ k ) tan(θ j − θ k ) tan(θ j + θ k )

where θ j and θ k are angle of incidence and angle of refraction respectively; t sjk and rjks are the transmission coefficient and reflection coefficient for s polarized light respectively; t jkp and rjkp are the transmission coefficient and reflection coefficient for p polarized light respectively. According to the transfer matrix method listed in the reference [19]. The total electrical field in the layer j can be expressed as follows:

E j ( x) = E +j ( x) + E −j ( x) = t +j E0+ + t −j E0+

(10)

The generation rate of the excitons is roughly equal to the time average of the energy dissipated in the layer j.

Q j ( x) = Q( x) = α TI 0 [e−α x + ρ ''j 2 ⋅ e

2 1 cε 0α jη j E j ( x) 2

−α (2 d j − x )

+ 2 ρ '' ⋅ e−α d ⋅ cos(

(11)

4πη

λ

(d − x) + δ )]

(12)

The detailed derivations of Eq. (12) are available in the reference [19]. Where α is absorption coefficient, T is internal intensity transmittance, I 0 is the intensity of the incident light, ρ '' is the absolute value of the complex reflection coefficient. 2.2 Electrical modeling 2.2.1 Exciton transport The electrical modeling is carried out in linear coordinates as opposed to cylindrical coordinates, since the active layer thickness, 20 nm, is orders of magnitude smaller than mm diameter cylindrical glass substrate. The exciton transport is dominated by the diffusion equation. At steady state condition if n is the exciton density, the continuity equation is:

∂n ∂ 2 n n Q( x) =D 2 − + =0 ∂t ∂x τ hν

#164841 - $15.00 USD (C) 2012 OSA

(13)

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A611

where D is the diffusion coefficient of the exciton, τ is the mean lifetime of the exciton, h is Planck constant, ν is the frequency of the incident light, and Q(x) is exciton generation rate in Eq. (12). The boundary conditions for Eq. (13) are: (i) the surfaces of the active layer have infinite recombination velocity and act as perfect excitons’ sinks; (ii) all excitons can either recombine or dissociate into free charge carriers at the donor-accepter interface. By applying two boundary conditions into the general solution, one could get the exciton distribution as follows: n( x ) =

θ1α TN  − β x 4πη  (d − x) + δ '' )  (14) Ae + Be β x + e −α x + C1eα x + C2 cos( 2 2  λ D( β − α )  

The photo-generated current density due to the charge transfer process at the donor/acceptor interface is defined as:

I =q

L2 ∂n ( ) τ ∂x x = xDA

(15)

where L is the diffusion length of the material and is the product of D and τ , q is elemental charge. The sum of the photocurrent generated at the donor and acceptor can be written as:

I cell = I D _ s + I D _ p + I A _ s + I A _ p

(16)

where I D _ s and I A _ s are the current generated by the donor and acceptor respectively for s polarized light, I D _ p and I A _ p are the current generated by the donor and acceptor respectively for p polarized light. These currents are expressed by the Eq. (15). 2.2.2 Absorption spectrum, EQE, IQE and efficiency The absorption spectrum, external quantum efficiency (EQE), internal quantum efficiency (IQE) and efficiency are the primary parameters reflecting the fundamental performance of the light trapping system. In this Section, we would derive the expressions for these significant parameters on the basis of the model described in the Section 2.1.2 and Section 2.2.1. The absorption spectrum is generally defined as the ratio of the number of the photons absorbed at the active layer to the number of photons incident on the system. For this light trapping system, the number of photons absorbed could be estimated by:

nabsorbed (λ ) = ∫

π /2

0





0

dA m

dD m

[ ∫ ∑ GA _ k ( x)dx + ∫ ∑ GD _ k ( x) dx]dθ d ε 0 k =1

(17)

0 k =1

where n absorbed (λ ) is the total photons absorbed at the active layer, dA and dD are the thickness of donor material and acceptor material, respectively, G A _ k (x ) and G D _ k (x ) are the exciton generation rate in the donor and acceptor material, respectively and they are calculated by the m

Eq. (12),

∑G

A _k

( x ) stands for the total sum of the exciton generation rate in the acceptor

k =1

because the light will bounce left and right in the light trapping structure until the light emits out from the system. The number of incident photons can be defined as:

nincident (λ ) =

#164841 - $15.00 USD (C) 2012 OSA

Wf (λ ) Aglass hυ

(18)

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A612

Therefore, the absorption spectrum could be expressed as follows:

AS (λ ) =

nabsorbed (λ ) = nincident (λ )

hυ i ∫

π /2

0





0

dA m

dD m

[ ∫ ∑ GA _ k ( x) dx + ∫ 0 k =1

∑G

D_k

( x)dx]dθ d ε

0 k =1

Wf (λ ) Aglass

(19)

The EQE is defined as the number of carriers collected by the electrode per photon incident on the OPV cell. Typically the external quantum efficiency is the product of five steps [Peter, et al., 2003]: (i) optical absorption efficiency of the active layer ηa , (ii) optical enhancement factor due to light trapping Γ , (iii) exciton collection efficiency at the donor-acceptor interface (D-A interface) ηc , (iv) exciton dissociation efficiency at the D-A interface by charge transfer process ηd , and (v) charge collection efficiency by the electrodes ηe . However, the efficiency of charge transfer process typically approaches 100% [Arkhipov, et al., 2003]. The carrier collection efficiency of bilayer OPV device at the short circuit condition is nearly 100% due to the extreme thin film and internal electric field [Peter, et al., 2003]. Therefore, it can be calculated by:

hυ EQE (λ ) = ηaηc Γ =

π / 2 2π m

I k (λ , θ , φ ) dθ dφ q k =1

∫ ∫∑ 0

0

WAglass f (λ )

(20)

where θ and φ are the altitude angle and azimuth angle, A glass is the cross section area of the I k (λ , θ , φ ) stands for the sum of q k =1 the current collected by the organic solar cell, until the light emits out of the light trapping structure. Since we have already derived the EQE and absorption spectrum, the IQE is easily defined as:

glass substrate, W is the incident irradiance (1 kW/m2),

IQE =

m



EQE (λ ) Absorption(λ )

(21)

Now, Short circuit current, Isc, could be easily calculated as the product of External quantum efficiency EQE and AM1.5 spectrum: I sc = q ∫ EQE (λ ) ⋅ W ⋅ Aglass ⋅ f (λ ) d λ

(22)

Since the organic solar cell is in nano-scale, the cross section area of the whole light trapping system is approximately equal to the cross section area of the glass. The efficiency of the whole system is defined as:

η=

I sc ⋅ Voc ⋅ FF = I sc = q ∫ EQE (λ ) ⋅ f (λ )d λ ⋅ Voc ⋅ FF W ⋅ Aglass

(23)

However, the packing density of the proposed cylindrical structure is less than the packing density of the planar structure, thus the effective efficiency compared with the control, in the module level, should be multiplied by a factor of π/4.

η mod =

#164841 - $15.00 USD (C) 2012 OSA

πη 4

(24)

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A613

3. Results and discussions Figure 1 shows the three dimensional (3-D) image of the proposed OPV device structure on cylindrical substrate which we have discussed above. Figure 2 is the cross sectional view of the proposed OPV device structure. An example of perpendicular incident light and the reflected light are shown in Fig. 2. Some reflected light would escape the OPV structure directly, some light would pass the organic solar cell layers several times and some light would be totally trapped by this structure. Figure 3 shows the 3-D view of an OPV structure with similar materials layers on a conventional flat plain glass substrate for comparison with the proposed structure. Figure 4 shows the simulated absorption spectrums of the proposed OPV structure and the conventional OPV on flat plain substrate without light trapping for active layer thicknesses of 20, 40 and 80 nm (d = dD + dA). The equal layer thicknesses of donor dD and acceptor dA are considered. The light absorption in the proposed OPV structure is significantly enhanced due to the cylindrical substrate structure, especially when the organic active layer is thin enough such as 20 nm. The effect of light trapping is greatly reduced for the thicker active layers as thick layer itself can absorb much more light than the thin layers. When the active layer thickness reaches 80 nm, the conventional OPV structure on plain substrate can absorb more light in the wavelength between 300 nm to 450 nm than the proposed structure due to the more reflection loss at the interfaces. The proposed OPV structure could greatly enhance the absorption of light and the performance of device for the ultra-thin organic active layer.

Fig. 4. Simulated absorption spectrum of proposed OPV device structure for 20(a), 40 (b) and 80 (c) nm thick active layers on 5 cm long and 1 mm diameter glass substrate. The absorption spectrum of conventional OPV device on flat plain glass substrate is also compared with the proposed structure. The structure layers thickness was considered as 160 nm ITO/32 nm PEDOT: PSS/(CuPc/PCBM)/12 nm BCP/100 nm Al.

#164841 - $15.00 USD (C) 2012 OSA

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A614

Figure 5 illustrates absorption spectrum of the proposed OPV cell for height 1, 3, 5 and 7 cm of the cylindrical glass substrate with an active layer thickness 20 nm. The absorption spectrum of the OPV structure on flat plain substrate is also compared with 20 nm thick active layer. The light absorption increases with the increase of cylindrical substrate height. As the substrate height increases, the absorption spectrum starts to saturate due to the inevitable reflection losses at the multi-interfaces. However, the increased height would definitely enhance the total cost of the proposed OPV system. Figure 6 shows the absorption spectrum of the proposed OPV cell for various diameters 1, 3, 4 and 6 mm of the glass substrate with an active layer thickness 20 nm. The light absorption increases with the decrease of cylindrical glass substrate for the proposed OPV structure.

Fig. 5. Simulated absorption spectrum of proposed OPV device structure for 1, 3, 5 and 7 cm heights of glass substrate with 1 mm diameter and 20 nm thick active layer. The absorption spectrum of conventional OPV device on flat plain glass substrate is also compared with the proposed structure. The structure layers thickness was considered as 160 nm ITO/32 nm PEDOT: PSS/20 nm (CuPc/PCBM)/12 nm BCP/100 nm Al.

#164841 - $15.00 USD (C) 2012 OSA

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A615

Fig. 6. Simulated absorption spectrum of proposed OPV device structure for 1, 3, 4 and 6 mm diameters of the glass substrate for fixed height 5 cm and 20 nm thick active layer. The absorption spectrum of conventional OPV device on flat plain glass substrate is also compared with the proposed structure. The structure layers thickness was considered as 160 nm ITO/32 nm PEDOT: PSS/20 nm (CuPc/PCBM)/12 nm BCP/100 nm Al.

The EQE is the comprehensive indicator of both the efficiency of the exciton absorption and collection. The IQE is the inherent property of the OPV solar cell and is therefore independent of the light trapping system. With perfect light trapping system, the EQE will approach to the IQE with almost complete absorption. We present the EQE and IQE curve of the proposed system to check its performance. From the Fig. 7, it is clear that the thinner active layer (20 nm) has the higher IQE of the proposed OPV device structure. The generated exciton in thin active layer has the highest rate of dissociation at the D-A interfaces as the distance of D-A interface within the exciton diffusion length. The charge collection probability will increase for thin active layer. The EQE depends on both light absorption ability and IQE. The highest EQE occurs for the active layer thickness 20 nm. The light trapping structure could effectively enhance the light absorption in the proposed OPV structure for thin active layer close to the diffusion length of exciton. The light trapping structure is not very effective for thicker active layer such as 80 nm as active layer itself absorb most of the incident light.

#164841 - $15.00 USD (C) 2012 OSA

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A616

Fig. 7. Simulated EQE spectrum of proposed OPV device structure for active layer thicknesses 20 (a), 40 (b) and 80 nm (c) with glass substrate height 5 cm and diameter 1mm. The EQE spectrum of conventional OPV device on flat plain glass substrate is also compared with the proposed structure. The structure layers thickness was considered as 160 nm ITO/32 nm PEDOT: PSS/CuPc/PCBM/12 nm BCP/100 nm Al.

Figure 8 shows the EQE and IQE curve for different height of the glass substrate of diameter 3 mm cm with active layer thickness 20 nm. Figure 9 shows the EQE and IQE curves for different diameter of the glass substrate with 5 cm height and active layer thickness 20 nm. Similar to the absorption curve, the EQE increases with the higher substrate height and lower diameter. However, when the height reaches 7 cm and the diameter decreases to 1 mm, the EQE begins to saturate due to the limitation of IQE. The highest EQE is far below IQE because of the inevitable reflection loss at the interfaces and the absorption loss by ITO and BCP.

#164841 - $15.00 USD (C) 2012 OSA

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A617

Fig. 8. Simulated EQE and IQE spectrum of proposed OPV device structure for 1, 3, 5, and 7 cm heights of glass substrate with fixed diameter 1mm. The EQE and IQE spectrum of conventional OPV device on flat plain glass substrate is also compared with the proposed structure. The structure layers thickness was considered as 160 nm ITO/32 nm PEDOT: PSS/20 nm (CuPc/PCBM)/12 nm BCP/100 nm Al.

Fig. 9. Simulated EQE and IQE spectrum of proposed OPV device structure for 1, 3, and 6 mm diameters of glass substrate with fixed height 5 cm. The EQE and IQE spectrum of conventional OPV device on flat plain glass substrate is also compared with the proposed structure. The structure layers thickness was considered as 160 nm ITO/32 nm PEDOT: PSS/20 nm (CuPc/PCBM)/12 nm BCP/100 nm Al.

The effective conversion efficiency, considering the effect of the different packing densities, of proposed OPV device versus active layer thickness curves for the different heights of glass substrate is shown in Fig. 10. The light trapping device improves the

#164841 - $15.00 USD (C) 2012 OSA

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A618

efficiency of OPV by a factor of 2 compared to the conventional OPV device structure on flat plain substrate. As expected, longer glass substrate could have the higher conversion efficiency of the proposed OPV structure, although the efficiency would saturate after the height increased to 9 cm. The active layer thickness 20 nm shows the highest efficiency of the OPV devices. The efficiency curves for different diameters of glass substrate as shown in Fig. 11 are similar to Fig. 10. It is very obvious that there is a discontinuity in the slope of the Figs. 10 and 11, due to the active layer thickness is equal to or less than 20 nm. Most of the excitens generated by photons absorption would be collected and dissociated at the D-A interface as the active layer thickness is nearly equal to the sum of the diffusion length of excitons in donor and acceptor materials. The diffusion probability of excitons to the D-A interface are almost equal to 1. When the active layer thickness is more than 20 nm, the collection probability of excitons to the D-A interface is sharply decreased. Figure 12 shows the effect of incident angle of light on the efficiency of proposed OPV device structure. The incident angle of light has a little effect on the performance of device. Most of the incident light on OPV device could be trapped by the cylindrical structure of substrate.

Fig. 10. The simulated effective power conversion efficiency, in the level of module, of the proposed OPV device structure for 1, 3, 5, 7 and 9 cm height of glass substrate with 1 mm diameter. The efficiency of conventional OPV device on flat plain glass substrate is also compared with the proposed structure. The structure layers thickness was considered as 160 nm ITO/32 nm PEDOT: PSS/20 nm (CuPc/PCBM)/12 nm BCP/100 nm Al.

#164841 - $15.00 USD (C) 2012 OSA

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A619

Fig. 11. The simulated effective power conversion efficiency of proposed OPV device structure, in the level of module, for 1, 3, 5, and 7 mm diameters glass substrate with 5 cm height. The conversion efficiency of conventional OPV device on flat plain glass substrate is also compared with the proposed structure. The structure layers thickness was considered as 160 nm ITO/32 nm PEDOT: PSS/20 nm (CuPc/PCBM)/12 nm BCP/100 nm Al.

Fig. 12. The simulated effective conversion efficiency, in the level of module, of proposed OPV device structure for different incident angle of light for 7, 9 and 10 cm height of glass substrate with 1 mm diameter. Positive and negative means the incident angle at left and right side of the normal line. The structure layers thickness was considered as 160 nm ITO/32 nm PEDOT: PSS/20 nm (CuPc/PCBM)/12 nm BCP/100 nm Al.

From the above analysis one can make the best efficiency OPV cells if it is fabricated on 5 cms long and 5 mm diameter cylindrical glass substrate with 20 nm thick active layers. For the same cross sectional areas of intersections light, cylindrical rods will have more surface area and more than double cell efficiency compare to flat surface cells. Cost effects of using long cylindrical rods will be much better than the flat substrate cells because of high efficiency cells.

#164841 - $15.00 USD (C) 2012 OSA

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A620

4. Conclusions We have designed and analyzed a proposed OPV device structure on cylindrical glass substrate to enhance the power conversion efficiency of organic solar cell. An optical and electrical analysis has been carried out for the proposed light trapping structure. Optical interference effects, increased polarization and non-normal incident light are considered for the analysis of the proposed device by expending the Pettersson’s method. Some significant device parameters such as IQE, EQE, efficiency and absorption spectrum are derived. The light absorption curve, EQE and efficiency are significantly enhanced by the proposed light trapping structure. The efficiency of proposed device has the potential to increase by a factor of 2 compared to the OPV device on flat plain glass substrate. The incident angle of light has a little effect on the performance of proposed OPV structure. The proposed OPV device structure can operate efficiently for all incident angle of light. The simulation program is developed in MATLAB by the author.

Acknowledgments The authors would like to thank the School of Photovoltaic and Renewable Energy for its support. We are all grateful for Dr. Patrick Campbell, Dr. Ivan Perez-Wurfl and all the members in PV school for their suggestions and support.

#164841 - $15.00 USD (C) 2012 OSA

Received 26 Mar 2012; revised 27 Jun 2012; accepted 30 Jun 2012; published 10 Jul 2012 10 September 2012 / Vol. 20, No. S5 / OPTICS EXPRESS A621