We propose a novel optical design of organic solar cell with a hybrid plasmonic system, which comprises a plasmonic cavity coupled with a dielectric core-metal shell nanosphere. From a rigorous solution of Maxwell’s equations, called volume integral equation method, optical absorption of the active polymer material has a four-fold increase. The significant enhancement mainly attributes to the coupling of symmetric surface wave modes supported by the cavity resonator. The dispersion relation of the plasmonic cavity is characterized by solving an 1D eigenvalue problem of the air/metal/polymer/metal/air structure with finite thicknesses of metal layers. We demonstrate that the optical enhancement strongly depends on the decay length of surface plasmon waves penetrated into the active material. Furthermore, the coherent interplay between the cavity and the dielectric core-metal shell nanosphere is undoubtedly confirmed by our theoretical model. The work offers detailed physical explanations to the hybrid plasmonic cavity device structure for enhancing the optical absorption of organic photovoltaics.
© 2011 OSA
Critically different from the thin-film polycrystalline or amorphous silicon solar cells (SCs) with active layer thickness of a few microns , the active polymer layer of thin-film organic solar cells (OSCs) only has few hundreds nanometers or even thinner thickness due to an extreme short exciton diffusion length . Such thin active layer with low refractive index (n ≈ 1.8) induces not only the weak optical absorption of OSCs but also fundamental (half-wavelength) limits of the optical design. On one hand, the strong Fabry-Pérot photonic mode or waveguide mode [3, 4] cannot be expected in the ultra-thin active layer. On the other hand, the physical mechanism of near-field concentration (not far-field scattering) should be taken into account in the design. Having unique features of tunable subwavelength resonance and engineered near-field concentration, plasmon is an enabling technique for light manipulation and management [5–7]. Plasmonic effects allow us to unprecedentedly improve the optical absorption of thin-film OSCs  and promote emerging SC technology meeting clean energy demands. Plasmonic nanostructures can offer three principles to enhance the optical absorption of OSCs. The first one is surface plasmon resonance by metal gratings fabricated on the top or bottom of the active layer [9–13]. The second one is local plasmon resonance by metal particles incorporated into or near the active layer [14–18]. The third one is plasmon coupling and hybridization, such as surface plasmon resonance coupled with local plasmon resonance or plasmon resonance coupled with photonic resonance .
In the paper, we propose a novel optical design of OSC with a hybrid plasmonic system, which comprises a plasmonic cavity coupled with a dielectric core-metal shell nanosphere. It has been investigated that optical absorption of the active polymer material has a four-fold increase. With the help of rigorous electrodynamic approaches, we unveil the fundamental physics of the significant enhancement, which mainly attributes to the coupling of symmetric surface wave modes supported by the cavity resonator with a long penetration length. Moreover, coherent interaction between the cavity and the dielectric core-metal shell nanosphere is definitely confirmed by our theoretical model. A distribution of polarization charges on the surface of the cavity indicates a bonding and antibonding coupling modes in the hybrid plasmonic system. The work introduces a new hybrid plasmonic cavity device structure to enhance the optical absorption of organic photovoltaics with detailed physical explanations.
2. Theoretical Model
As a rigorous solution to Maxwell’s equations, a volume integral equation (VIE) method is developed to characterize the optical absorption of OSCs. Considering non-magnetic optical materials with an arbitrary inhomogeneity profile, the VIE can be written as20]. Due to the hypersingularity of the Green’s tensor and spurious discontinuity of tangential E-field induced by the scalar (piecewise constant) basis functions, the DDA method cannot accurately characterize the subwavelength plasmonic physics  and breaks down in the multilayer device structure with high-contrast metallic materials. Here, we develop an alternate algorithm to bypass the difficulties. In our model, the polarization currents are expanded using the roof-top vector basis functions  and thus the continuity of normal current is naturally satisfied at the material interfaces. Furthermore, the hypersingular Green’s tensor is smoothened by using the finite-difference approximation.
From the VIE [Eq. (1)], the scattered electric field generated by the volumetric polarization current J can be written as
Using the rooftop basis functions to expand the unknown currents, we haveEq. (8) can be defined in the same way.
As a result, the discretized form for the operator in Eq. (5) can be written asEq. (6) can be discretized as
The computations of the discrete convolutions can be performed efficiently by means of cyclic convolutions and fast Fourier transform (FFT) , which is similar to the DDA method. As a traditional iterative solver of the resulting VIE matrix equation, the conjugate-gradient method converges very slowly and will produce the non-physical random errors in the calculation of optical absorption. To tackle the problem, we employ the fast and smoothly converging biconjugate gradient stabilized (BI-CGSTAB) method . The FFT is adopted to accelerate the matrix-vector multiplications encountered in the BI-CGSTAB solver with computational complexity of O(N logN) and memory of O(N).
Through the rigorous VIE solution to Maxwell’s equations, we can access some important physical quantities to reveal the physical mechanism of plasmonic effects in OSCs and optimize device performances. The absorption spectrum of OSCs is calculated byEqs. (14) and (15), the SEF and TEF can be easily accessible. In addition, both the two layered OSC structure and embedded concentrators are meshed by the volumetric cuboid cells with a uniform grid size of 1 nm.
3. Results and Discussions
Figure 1 shows the schematic pattern of a heterojunction OSC. A hybrid plasmonic system, which comprises a plasmonic cavity coupled with a dielectric core-metal shell nanosphere, is employed for improving the optical absorption of the active polymer material. A transparent spacer is inserted to avoid local shunt and extract carriers. The incident light is propagated from the spacer to the active layer at the vertical incident angle with an E-field polarized along the x direction. Figure 2(a) shows the real and imaginary parts of the refractive index of the active material measured from ellipsometry . The complex refractive index of metal silver can be expressed by the Brendel-Bormann model . The SiO2 as a dielectric layer is adopted for the core-shell sphere and its refractive index can be found in the literature .
First, various nanosphere concentrators (excluding the plasmonic cavity) are systematically and comparatively observed. These nanospheres include a dielectric sphere, a metal sphere, a metal core-dielectric shell (MC-DS) sphere, and a dielectric core-metal shell (DC-MS) sphere . The scattering cross section (SCS) of the nanospheres can be obtained from the generalized reflection coefficients of the spherically layered media [29, 30]Fig. 1, respectively, m is the order of the modified spherical Bessel (Hankel) functions, and and are the generalized reflection coefficients of the TM and TE spherical waves in the layer 3 reflected by the layer 2. For small spherical particles, the leading term (m = 1) of determines the value of the SCS. The generalized reflection coefficient can be written as a recursive equation Figure 2(b) and 2(c) show the SEFs and the SCS of the nanospheres, respectively. As seen in Figs. 2(b) and 2(c), the peaks of the SEFs agree with those of the SCS well. The small dielectric nanosphere, although has no loss and large refractive index (n=4), is not a good concentrator for OSCs. The dielectric nanosphere with positive refractive index cannot produce a strong dipole resonance compared to the metal nanosphere. Moreover, in contrast to the DC-MS sphere that has a metal layer adjacent to different materials (SiO2 and polymer), the resonance of the MC-DS sphere is blue-shifted because only one material SiO2 with lower refractive index is adjacent to the metal layer. The near field of the MC-DS sphere confines to the shell layer and cannot sufficiently scatter to the active layer. As a result, the optical enhancement by the MC-DS sphere is very weak. Figure 2(d) shows a tunable plasmon resonance by engineering the geometry of the DC-MS sphere. The resonance is red-shifted and becomes damped as the core radius increases.
Second, we investigate the plasmonic silver cavity structure enclosing the active polymer material. Here the nanosphere is excluded. Figure 3(a) shows the SEF by the plasmonic cavity. A broadband enhancement has a good overlap with the absorption spectrum of the active material as depicted in Fig. 1(a). The dispersion relation of the cavity resonator plays a key role in unveiling the fundamental physics of the significant optical enhancement. Different from previous works assuming infinitely thick metal claddings , we take into account the finite-thickness effect of the metal claddings. The governing equation for an 1D eigenvalue problem of the Air/Ag/polymer/Ag/Air planar layered media is given by Eq. (23)] can be easily solved by the finite-difference method  with perfectly matched layer absorbing boundary conditions [32, 33]. We find two eigenvalues of the TM wave closest to the momentum of surface plasmon polariton (SPP). The SPP propagates at the interface between semi-infinite polymer and Ag half-spaces and its dispersion relation is of the form 31]. These modes also are associated with the guided-wave poles of the generalized reflection coefficient for the planar layered media . Figure 3(b) shows the dispersion relations. At long wavelengths, the dispersion relation of the SPP is located between those of symmetric and asymmetric modes. Then we calculate the decay lengths of the surface plasmon waves penetrated into the active material as plotted in Fig. 3(c). The peaks of the decay length of the symmetric mode strongly coincide with those of the SEF as shown in Fig. 3(a). The incident light excites symmetric surface plasmon waves at the bottom corners of the cavity and the waves propagate in the active layer bounded between the metal claddings. A short decay length away from the metal claddings makes E-fields concentrate at the surfaces of the metals. Contrarily, a long decay length induces a concentrated E-field in the center of the active material. Figure 4 demonstrates the near-field distributions in the active layer at the wavelengths denoted with the arrows of Fig. 3(a). At 800 nm, the slowly decaying near field away from the metal claddings leads to the most significant enhancement.
Third, we study the hybrid plasmonic system, which comprises a plasmonic cavity coupled with a nanosphere. The optical enhancement shows little improvements when the cavity is coupled with the dielectric or MC-DS sphere as illustrated in Fig. 3(a). However, a significant enhancement can be achieved if the cavity is coupled with metal or DC-MS sphere as shown in Fig. 3(d). By evaluating the TEF with the aid of Eq. (15), the optical absorption of the active polymer material has a four-fold increase when the cavity is coupled with the DC-MS sphere. We calculate the algebraic summation of the SEF by the uncoupled single cavity and that by the uncoupled single DC-MS sphere. As shown in Fig. 5, the summation is observably smaller than the SEF by the cavity coupled with the DC-MS sphere at long wavelengths ranging from 620 nm to 800 nm. The coherent interplay between the cavity and the DC-MS sphere is undoubtedly confirmed by the result. Then, we analyze the metal dissipation ratio defined as the metal loss of the coupled concentrator over that of the uncoupled one. The inset of Fig. 5 shows that the metal dissipation ratio of the DC-MS sphere is substantially larger than that of the cavity. The evanescent SPPs from the inner surface of metallic cavity sufficiently penetrate the DC-MS sphere inducing stronger local plasmon resonance as well as larger metallic loss. Thus the DC-MS sphere becomes more effective concentrator when it is coupled with the cavity. Finally, we discuss the coupling modes in the hybrid plasmonic system. Figure 6 depicts a polarization charge distribution on the surface of the cavity by using Eq. (16), indicating a bonding and antibonding coupling modes in the hybrid system. The bonding modes at 620 nm and 750 nm have a denser charge distribution if the distance between the cavity and the shell surfaces becomes closer. However, the antibonding mode at 650 nm reverses the polarity of the polarization charge when the face-to-face distance approaches the minimum. The near fields of the bonding modes show more concentrated field at the gap between the DC-MS sphere and the cavity, which can be seen in the inset of Fig. 5. For the bonding mode, the polarity of its polarization charge is marked in Fig. 1. Due to the in-phase plasmon oscillation and hybridization, the bonding mode is superradiant or strongly radiative, and provides a great help for the optical enhancement.
The hybrid plasmonic system, which comprises the plasmonic cavity coupled with the DC-MS nanosphere, can increase the optical absorption of the OSC by four fold. The significant enhancement mainly results from the coupling of symmetric surface wave modes supported by the cavity resonator and strongly depends on the decay length of surface plasmon waves penetrated into the active layer. Furthermore, the coherent interplay between the cavity and the DC-MS nanosphere is strongly demonstrated by our theoretical model. The bonding coupling mode in the hybrid plasmonic system enhances the optical absorption further. The work provides detailed physical explanations for the hybrid plasmonic cavity device structure to enhance the optical absorption of organic photovoltaics.
The authors acknowledge the support of the Research Grant Council of the Hong Kong Special Administrative Region, China (grants No. 712108, No. 711508, and No. 711609). This project is also supported in part by a Hong Kong UGC Special Equipment Grant ( SEG HKU09).
References and links
1. K. L. Chopra, P. D. Paulson, and V. Dutta, “Thin-film solar cells: an overview,” Prog. Photovoltaics 12, 69–92 (2004). [CrossRef]
2. H. Hoppe and N. S. Sariciftci, “Organic solar cells: an overview,” J. Mater. Res. 19, 1924–1945 (2004). [CrossRef]
3. R. A. Pala, J. White, E. Barnard, J. Liu, and M. L. Brongersma, “Design of plasmonic thin-film solar cells with broadband absorption enhancements,” Adv. Mater. 21, 3504–3509 (2009). [CrossRef]
6. J. A. Schuller, E. S. Barnard, W. S. Cai, Y. C. Jun, J. S. White, and M. L. Brongersma, “Plasmonics for extreme light concentration and manipulation,” Nat. Mater. 9, 193–204 (2010). [CrossRef] [PubMed]
7. S. A. Maier, Plasmonics: Fundamentals and Applications (Springer, 2007).
9. A. J. Morfa, K. L. Rowlen, T. H. Reilly, M. J. Romero, and J. van de Lagemaat, “Plasmon-enhanced solar energy conversion in organic bulk heterojunction photovoltaics,” Appl. Phys. Lett. 92, 013504 (2008). [CrossRef]
10. W. L. Bai, Q. Q. Gan, G. F. Song, L. H. Chen, Z. Kafafi, and F. Bartoli, “Broadband short-range surface plasmon structures for absorption enhancement in organic photovoltaics,” Opt. Express 18, A620–A630 (2010). [CrossRef] [PubMed]
11. M. G. Kang, T. Xu, H. J. Park, X. G. Luo, and L. J. Guo, “Efficiency enhancement of organic solar cells using transparent plasmonic Ag nanowire electrodes,” Adv. Mater. 22, 4378–4383 (2010). [CrossRef] [PubMed]
12. C. J. Min, J. Li, G. Veronis, J. Y. Lee, S. H. Fan, and P. Peumans, “Enhancement of optical absorption in thin-film organic solar cells through the excitation of plasmonic modes in metallic gratings,” Appl. Phys. Lett. 96, 133302 (2010). [CrossRef]
14. D. Duche, P. Torchio, L. Escoubas, F. Monestier, J. J. Simon, F. Flory, and G. Mathian, “Improving light absorption in organic solar cells by plasmonic contribution,” Sol. Energy Mater. Sol. Cells 93, 1377–1382 (2009). [CrossRef]
15. J. Jung, T. G. Pedersen, T. Sondergaard, K. Pedersen, A. N. Larsen, and B. B. Nielsen, “Electrostatic plasmon resonances of metal nanospheres in layered geometries,” Phys. Rev. B 81, 125413 (2010). [CrossRef]
16. S. J. Tsai, M. Ballarotto, D. B. Romero, W. N. Herman, H. C. Kan, and R. J. Phaneuf, “Effect of gold nanopillar arrays on the absorption spectrum of a bulk heterojunction organic solar cell,” Opt. Express 18, A528–A535 (2010). [CrossRef] [PubMed]
17. I. Diukman, L. Tzabari, N. Berkovitch, N. Tessler, and M. Orenstein, “Controlling absorption enhancement in organic photovoltaic cells by patterning Au nano disks within the active layer,” Opt. Express 19, A64–A71 (2011). [CrossRef] [PubMed]
18. J. Jung, T. Sondergaard, T. G. Pedersen, K. Pedersen, A. N. Larsen, and B. B. Nielsen, “Dyadic Green’s functions of thin films: Applications within plasmonic solar cells,” Phys. Rev. B 83, 085419 (2011). [CrossRef]
19. V. E. Ferry, L. A. Sweatlock, D. Pacifici, and H. A. Atwater, “Plasmonic nanostructure design for efficient light coupling into solar cells,” Nano Lett. 8, 4391–4397 (2008). [CrossRef]
20. B. T. Draine and P. J. Flatau, “Discrete-dipole approximation for scattering calculations,” J. Opt. Soc. Am. A 11, 1491–1499 (1994). [CrossRef]
21. A. M. Kern and O. J. F. Martin, “Surface integral formulation for 3D simulations of plasmonic and high permittivity nanostructures,” J. Opt. Soc. Am. A 26, 732–740 (2009). [CrossRef]
22. A. W. Glisson and D. R. Wilton, “Simple and efficient numerical methods for problems of electromagnetic radiation and scattering from surfaces,” IEEE Trans. Antennas Propag. 28, 593–603 (1980). [CrossRef]
23. M. F. Catedra, E. Gago, and L. Nuno, “A numerical scheme to obtain the RCS of three-dimensional bodies of resonant size using the conjugate gradient method and the fast Fourier transform,” IEEE Trans. Antennas Propag. 37, 528–537 (1989). [CrossRef]
24. H. A. Vandervorst, “Bi-CGSTAB: A fast and smoothly converging variant of Bi-CG for the solution of nonsymmetric linear systems,” SIAM J. Sci. Stat. Comput. 13, 631–644 (1992). [CrossRef]
25. W. C. H. Choy and H. H. Fong, “Comprehensive investigation of absolute optical properties of organic materials,” J. Phys. D: Appl. Phys. 41, 155109 (2008). [CrossRef]
26. A. D. Rakic, A. B. Djurisic, J. M. Elazar, and M. L. Majewski, “Optical properties of metallic films for vertical-cavity optoelectronic devices,” Appl. Opt. 37, 5271–5283 (1998). [CrossRef]
27. E. D. Palik, Handbook of Optical Constants of Solids (Academic Press, 1998).
29. W. C. Chew, Waves and Fields in Inhomogenous Media (Wiley-IEEE Press, 1999). [CrossRef]
30. L. Tsang, J. A. Kong, and K. H. Ding, Scattering of Electromagnetic Waves: Theories and Applications (Wiley, 2000). [CrossRef]
31. B. Prade, J. Y. Vinet, and A. Mysyrowicz, “Guided optical waves in planar heterostructures with negative dielectric-constant,” Phys. Rev. B 44, 13556–13572 (1991). [CrossRef]
32. J. P. Berenger, “A perfectly matched layer for the absorption of electromagnetic-waves,” J. Comput. Phys. 114, 185–200 (1994). [CrossRef]
33. W. C. Chew and W. H. Weedon, “A 3-D perfectly matched medium from modified Maxwell’s equations with stretched coordinates,” Microw. Opt. Technol. Lett. 7, 599–604 (1994). [CrossRef]