×

Derivation and solution of effective medium equations for bulk heterojunction organic solar cells. (English) Zbl 1380.82057

Summary: A drift-diffusion model for charge transport in an organic bulk heterojunction solar cell, formed by conjoined acceptor and donor materials sandwiched between two electrodes, is formulated. The model accounts for (i) bulk photogeneration of excitons, (ii) exciton drift and recombination, (iii) exciton dissociation (into polarons) on the acceptor-donor interface, (iv) polaron recombination, (v) polaron dissociation into a free electron (in the acceptor) and a hole (in the donor), (vi) electron/hole transport and (vii) electron-hole recombination on the acceptor-donor interface. A finite element method is employed to solve the model in a cell with a highly convoluted acceptor/donor interface. The solutions show that, with physically realistic parameters, and in the power generating regime, the solution varies little on the scale of the micro-structure. This motivates us to homogenise over the micro-structure; a process that yields a far simpler one-dimensional effective medium model on the cell scale. The comparison between the solution of the full model and the effective medium (homogenised) model is very favourable for applied voltages less than the built-in voltage (the power generating regime) but breaks down as the applied voltages increases above it. Furthermore, it is noted that the homogenisation technique provides a systematic way to relate effective medium modelling of bulk heterojunctions to a more fundamental approach that explicitly models the full micro-structure and that it allows the parameters in the effective medium model to be derived in terms of the geometry of the micro-structure. Finally, the effective medium model is used to investigate the effects of modifying the micro-structure geometry, of a device with an interdigitated acceptor/donor interface, on its current-voltage curve.

MSC:

82D37 Statistical mechanics of semiconductors
65N22 Numerical solution of discretized equations for boundary value problems involving PDEs
35B27 Homogenization in context of PDEs; PDEs in media with periodic structure

Software:

WIAS-TeSCA
Full Text: DOI

References:

[1] [1]AllaireG. (1992) Homogenization and two-scale convergence. SIAM J. Math. Anal.23, 1482-1518.10.1137/0523084 · Zbl 0770.35005 · doi:10.1137/0523084
[2] [2]AllsopN., NürnbergR., Lux-SteinerM. Ch. & Schedel-NiedrigTh. (2009) Three-dimensional simulations of a thin film heterojunction solar cell with a point contact/defect passivation structure at the heterointerface. Appl. Phys. Lett.95, 122108-1-122108-3.
[3] [3]BarkerJ. A., RamsdaleC. M. & GreenhamN. C. (2003) Modelling the current-voltage characteristic of bilayer polymer devices. Phys. Rev. B67, 075205.10.1103/PhysRevB.67.075205 · doi:10.1103/PhysRevB.67.075205
[4] [4]BarrettJ. W. & ElliottC. M. (1987) Fitted and unfitted finite-element methods for elliptic equations with smooth interfaces. IMA J. Numer. Anal.7283-300.10.1093/imanum/7.3.283 · Zbl 0629.65118 · doi:10.1093/imanum/7.3.283
[5] [5]BraunC. L. (1984) Electric-field assisted dissociation of charge-transfer states as a mechanism for photocarrier production. Chem. Phys.80, 4157-4161.
[6] [6]BrunaM. & ChapmanS. J. (2015) Diffusion in spatially varying porous media. SIAM J.Appl. Maths.75, 1648-1674.10.1137/141001834 · Zbl 1320.35038 · doi:10.1137/141001834
[7] [7]BrinkmanD., FellnerK., MarkowichP. A. & WolframM.-T. (2013) A drift-diffusion-reaction model for excitonic photovoltaic bilayers: Asymptotic analysis and a 2-D HDG finite-element scheme. Math. Models Methods Appl. Sci.23, 839-872.10.1142/S0218202512500625 · Zbl 1264.35118 · doi:10.1142/S0218202512500625
[8] [8]BuxtonG. A. & ClarkeN. (2006) Predicting structure and property relations in polymeric photovoltaic devices. Phys. Rev. B74, 085207.10.1103/PhysRevB.74.085207 · doi:10.1103/PhysRevB.74.085207
[9] [9]BuxtonG. A. & ClarkeN. (2007) Computer simulation of polymer solar cells. Model. Simul. Mater. Sci. Eng.15, 13-26.10.1088/0965-0393/15/2/002 · doi:10.1088/0965-0393/15/2/002
[10] [10]ChenJ.-D., CuiC., LiY.-Q., ZhouL., OuQ.-D., LiC., LiY. & TangJ.-X. (2015) Single-junction polymer solar cells exceeding 10 · doi:10.1002/adma.201404535
[11] [11]ColeJ. D. (1995) Limit process expansions and homogenization. SIAM J. Appl. Math.55, 410-424.10.1137/S0036139993269230 · Zbl 0819.35012 · doi:10.1137/S0036139993269230
[12] [12]ClarkeT. M. & DurrantJ. R. (2010) Charge photogeneration in organic solar cells. Chem. Rev.110, 6736-6767.10.1021/cr900271s · doi:10.1021/cr900271s
[13] [13]CloverI. (2016) Heliatek raises bar for OPV efficiency to 13.2. pv magazine.
[14] [14]CioranescuD. & DonatoP. (1999) An Introduction to Homogenization, Oxford Lecture Series in Mathematics and its Applications, Oxford, Oxford University Press. · Zbl 0939.35001
[15] [15]CredgingtonD., KimY., LabramJ., AnthopoulosT. D. & DurrantJ. (2011) Analysis of recombination in polymer C60 solar cells. J. Phys. Chem. Lett.2, 2759.10.1021/jz201338d · doi:10.1021/jz201338d
[16] [16]CredgingtonD., JamiesonF. C., WalkerB., NguyenT.-Q. & DurrantJ. R. (2012) Quantification of geminate and non-geminate recombination losses within a solution-processed small-molecule bulk heterojunction solar cell. Adv. Mater.24, 2135-2141.10.1002/adma.201104738 · doi:10.1002/adma.201104738
[17] [17]CroneB. K., DavidsP. S., CampbellI. H. & SmithD. L. (2000) Device model investigation of bilayer organic light emitting diodes. J. Appl. Phys.87, 1974.10.1063/1.372123 · doi:10.1063/1.372123
[18] [18]DavidsP. S., CampbellI. H. & SmithD. L. (1997) Device model for single carrier organic diodes. J. Appl. Phys.82, 6319.10.1063/1.366522 · doi:10.1063/1.366522
[19] [19]de FalcoC., SaccoR. & VerriM. (2010) Analytical and numerical study of photocurrent transients in organic polymer solar cells. Comput. Methods Appl. Mech. Eng.199, 1722-1732.10.1016/j.cma.2010.01.018 · Zbl 1231.78038 · doi:10.1016/j.cma.2010.01.018
[20] [20]DeibelC. & DyakonovV. (2010) Polymer-fullerene bulk heterojunction solar cells. Rep. Prog. Phys.73, 096401.10.1088/0034-4885/73/9/096401 · doi:10.1088/0034-4885/73/9/096401
[21] [21]FosterJ. M., KirkpatrickJ. & RichardsonG. (2013) Asymptotic and numerical prediction of current-voltagencurves for an organic bilayer solar cell under varying illumination and comparison to the Shockley equivalent circuit. J. Appl. Phys.114, 104501.10.1063/1.4820567 · doi:10.1063/1.4820567
[22] [22]GajewskiH., KaiserH. Chr., LangmachH., NürnbergR. & RichterR. H. (2003) Mathematical modelling and numerical simulation of semiconductor detectors. In: W.Jäger & H. J.Krebs (editors), Mathematics? Key Technology for the Future, Springer, Berlin, Heidelberg, pp. 355-364.10.1007/978-3-642-55753-8 · Zbl 1016.00014 · doi:10.1007/978-3-642-55753-8
[23] [23]GajewskiH.et al.TeSCA Two- and Three-Dimensional Semi-Conductor Analysis Package, Weierstrass Institute for Applied Analysis and Stochastics, Berlin.
[24] [24]GünesS., NeugebauerH. & SariciftciN. S. (2007) Conjugated polymer-based organic solar cells. Chem. Rev.103, 1324.
[25] [25]GreggK. A. & HannaM. C. (2003) Comparing organic to inorganic photovoltaic cells: Theory, experiment, and simulation. J. Appl. Phys.93, 3605-3614.10.1063/1.1544413 · doi:10.1063/1.1544413
[26] [26]GrovesC., BlakesleyJ. C. & GreenhamN. C. (2010) Effect of charge trapping on geminate recombination and polymer solar cell performance. Nano Lett.10, 1063-1069.10.1021/nl100080r · doi:10.1021/nl100080r
[27] [27]GrovesC., KimberR. G. E. & WalkerA. B. (2010) Simulation of loss mechanisms in organic solar cells. J. Chem. Phys.133, 144110.10.1063/1.3483603 · doi:10.1063/1.3483603
[28] [28]HoppeH. & SariciftciN. S. (2004) Organic solar cells: An overview. J. Mater. Res.19, 1924-1945.10.1557/JMR.2004.0252 · doi:10.1557/JMR.2004.0252
[29] [29]de JonghP. E. & VanmaekelberghD. (1996) Trap-limited transport in assemblies of nanometer-size TiO_2 particles. Phys. Rev. Lett.77, 3427-3430.10.1103/PhysRevLett.77.3427 · doi:10.1103/PhysRevLett.77.3427
[30] [30]KellerJ. B. (1980) Darcy’s law for flow in porous media and the two-space method. In: Lecture Notes in Pure and Applied Mathematics vol. 54, Dekker, New York. · Zbl 0439.76017
[31] [31]KellerJ. B. (1977) Effective behavior of heterogeneous media. In: U.Landman (editor), Statistical Mechanics and Statistical Methods in Theory and Application, Plenum, New York, pp. 631-644.10.1007/978-1-4613-4166-6 · doi:10.1007/978-1-4613-4166-6
[32] [32]KimberR. G. E., WrightE. N., O’KaneS. E. J., WalkerA. B. & BlakesleyJ. C. (2012) Mesoscopic kinetic Monte Carlo modeling of organic photovoltaic device characteristics. Phys. Rev. B86, 235206.10.1103/PhysRevB.86.235206 · doi:10.1103/PhysRevB.86.235206
[33] [33]KirchartzT., PietersB. E., KirkpatrickJ., RauU. & NelsonJ. (2011) Recombination via tail states in polythiophene: Fullerene solar cells. Phys. Rev. B83, 115209.10.1103/PhysRevB.83.115209 · doi:10.1103/PhysRevB.83.115209
[34] [34]KirkpatrickJ., MarconV., KremerK., Nelson, J. & AndrienkoD. (2007) Charge mobility in discotic mesophases: A multiscale quantum and classical study. Phys. Rev. Lett.98, 227402.10.1103/PhysRevLett.98.227402 · doi:10.1103/PhysRevLett.98.227402
[35] [35]KodaliH. K. & GanapathysubramanianB. (2012) Computer simulation of heterogeneous polymer photovoltaic devices. Model. Simul. Mater. Sci. Eng.20, 035015.10.1088/0965-0393/20/3/035015 · doi:10.1088/0965-0393/20/3/035015
[36] [36]KosterL. J. A., SmitsE. C. P., MihailetchiV. D. & BlomP. W. M. (2005) Device model for the operation of polymer/fullerene bulk heterojunction solar cells. Phys. Rev. B.72, 085205.10.1103/PhysRevB.72.085205 · doi:10.1103/PhysRevB.72.085205
[37] [37]KotlarskiJ. D., BlomP. W., KosterL. J., Lenes, M. & SloofL. H. (2008) Combined optical and electrical modeling of polymer: Fullerene bulk heterojunction solar cells. J. Appl. Phys.103, 084502.10.1063/1.2905243 · doi:10.1063/1.2905243
[38] [38]MartinC. M., BurlakovV. M. & AssenderH. E. (2006) Modellng charge transport in composite solar cells. Sol. Energy Mater. Sol. Cells90, 900-915.10.1016/j.solmat.2005.05.009 · doi:10.1016/j.solmat.2005.05.009
[39] [39]MartinC. M., BurlakovV. M., AssenderH. E. & BarkhouseD. A. R. (2007) A numerical model for explaining the role of the interface morphology in composite solar cells. J. Appl. Phys.102, 104506.10.1063/1.2817603 · doi:10.1063/1.2817603
[40] [40]MarkovD. E., AmsterdamE., BlomP. W. M., SievalA. B. & HummelenJ. C. (2005) Accurate measurement of the exciton diffusion length in a conjugated polymer using a heterostructure with a side-chain cross-linked fullerene layer. J. Phys. Chem. A109, 5266-5274.10.1021/jp0509663 · doi:10.1021/jp0509663
[41] [41]McNeillC. R., WestenhoffS., GrovesC., FriendR. H. & GreenhamN. C. (2007) Influence of nanoscale phase separation on the charge generation dynamics and photovoltaic performance of conjugated polymer blends: Balancing charge generation and separation. J. Phys. Chem. C111, 19153-19160.10.1021/jp075904m · doi:10.1021/jp075904m
[42] [42]NelsonJ. (2003) Diffusion-limited recombination in polymer-fullerene blends and its influence on photocurrent collection. Phys. Rev. B67, 155209.10.1103/PhysRevB.67.155209 · doi:10.1103/PhysRevB.67.155209
[43] [43]NelsonJ. (2003) The Physics of Solar Cells, London, Imperial College Press.
[44] [44]OffermansT., MeskersS. C. J. & JanssenR. A. J. (2005) Monte-Carlo simulations of geminate electron-hole pair dissociation in a molecular heterojunction: A two-step dissociation mechanism. Chem. Phys.308, 125-133.10.1016/j.chemphys.2004.08.015 · doi:10.1016/j.chemphys.2004.08.015
[45] [45]PautmeierL., RichertR. & BässlerH. (1990) Poole-Frenkel behaviour of charge transport in organic solids with off-diagonal disorder studied by Monte Carlo simulation. Synth. Met.37, 271.10.1016/0379-6779(90)90158-H · doi:10.1016/0379-6779(90)90158-H
[46] [46]PeumansP., UchidaS. & ForrestS. R. (2003) Efficient bulk heterojunction photovoltaic cells using small-molecular-weight organic thin films. Nature425, 158-162.10.1038/nature01949 · doi:10.1038/nature01949
[47] [47]PotscavageW. J., YooS. & KippelenB. (2008) Origin of the open-circuit voltage in a multilayer heterojunction organic solar cells. Appl. Phys. Lett.93, 193308.10.1063/1.3027061 · doi:10.1063/1.3027061
[48] [48]RichardsonG., DenuaultG. & PleaseC. P. (2012) Multiscale modelling and analysis of lithium-ion battery charge and discharge. J. Eng. Math.72, 41-72.10.1007/s10665-011-9461-9 · Zbl 1254.78026 · doi:10.1007/s10665-011-9461-9
[49] [49]RichardsonG., PleaseC. P., FosterJ. & KirkpatrickJ. A. (2012) Asymptotic solution of a model for bilayer organic diodes and solar cells. SIAM J. Appl. Math.72, 1792-1817.10.1137/110825807 · Zbl 1276.82064 · doi:10.1137/110825807
[50] [50]RichardsonG. & ChapmanS. J. (2011) Derivation of the bidomain equations for a beating heart with a general microstructure. SIAM J. Appl. Math.71, 657-675.10.1137/090777165 · Zbl 1243.35015 · doi:10.1137/090777165
[51] [51]ScottJ. C. & MalliarasG. G. (1999) Charge injection and recombination at the metal-organic interface. Chem. Phys. Lett.299, 115-119.10.1016/S0009-2614(98)01277-9 · doi:10.1016/S0009-2614(98)01277-9
[52] [52]SeunhyupY., PotscavageW. J., DomercquaB., LicT. D., JonesS. C., SzozskiewiczR., LevibD., RiedocE., MarderS. R. & KillenB. (2007) Analysis of improved photovoltaic properties of pentacene/C60 organic solarcells: Effects of exciton blocking layer thickness and thermal annealing. Solid-State Electron.51, 1367.10.1016/j.sse.2007.07.038 · doi:10.1016/j.sse.2007.07.038
[53] [53]RaoA., WilsonM. W. B., HodgkissJ. M., Albert-SeifriedS., BässlerH. & FriendR. H. (2010) Exciton fission and charge generation via triplet excitons in pentacene/C_60 bilayers. J. Am. Chem. Soc.132, 12698-12703.10.1021/ja1042462 · doi:10.1021/ja1042462
[54] [54]ScharfetterD. L. & GummelH. K. (1969) Large-signal analysis of a silicon Read diode oscillator. IEEE Trans. Electron. Dev.16, 64-77.10.1109/T-ED.1969.16566 · doi:10.1109/T-ED.1969.16566
[55] [55]SzeS. M. & KwokK. Ng (2006) Physics of Semiconductor Devices, 3rd ed., Wiley-Interscience, New York.
[56] [56]TansaeC., BlomP. W. M., de LeeuwD. M. & MeijerE. J. (2004) Charge carrier density dependence of mobility in poly-p-phenylene vinylene. Phys. Status Solidi B201, 1236.10.1002/pssa.200404340 · doi:10.1002/pssa.200404340
[57] [57]VerlaakS., BeljonneD., CheynsD., RolinC., LinaresM., CastetF., CornilJ. & HeremansP. (2009) Electronic structure and geminate pair energetics at organic-organic interfaces: The case of pentacene/C_60 heterojunctions. Adv. Funct. Mater.19, 3809-3814.10.1002/adfm.200901233 · doi:10.1002/adfm.200901233
[58] [58]WilliamsJ. & WalkerA. B. (2008) Two-dimensional simulations of bulk heterojunction solar cell characteristics. Nanotechnology19, 424011.10.1088/0957-4484/19/42/424011 · doi:10.1088/0957-4484/19/42/424011
[59] [59]WagenpfahlA., RauhD., BinderM., DeibelC. & DyakonovV. (2010) S-shaped current-voltage characteristics of organic solar devices. Phys. Rev. B82, 115306.10.1103/PhysRevB.82.115306 · doi:10.1103/PhysRevB.82.115306
[60] [60]YangF., ShteinM. & ForrestS. R. (2005) Controlled growth of a molecular bulk heterojunction photovoltaic cell. Nat. Mat.4, 37-41.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.