Issue 
EPJ Photovolt.
Volume 8, 2017
Topical Issue: Theory and modelling



Article Number  85503  
Number of page(s)  8  
Section  Theory and modelling  
DOI  https://doi.org/10.1051/epjpv/2017004  
Published online  16 June 2017 
https://doi.org/10.1051/epjpv/2017004
Quantum modeling of twolevel photovoltaic systems
^{1} Institut NÉEL, CNRS and Université Grenoble Alpes, 38042 Grenoble, France
^{2} Research Institute for Applied Physics and Astronomy, University of Tabriz, 5166614766 Tabriz, Iran
^{3} School of Electrical, Electronic and Computer Engineering, The University of Western Australia, WA 6009 Crawley, Australia
^{4} Département de Chimie, Université de Montréal, C.P. 6128 Succursale A, Montréal, Québec H3C 3J7, Canada
^{a}
email: tahereh.nematiaram@neel.cnrs.fr
Received: 21 November 2016
Accepted: 3 May 2017
Published online: 16 June 2017
We present a quantum formalism that provides a quantitative picture of the fundamental processes of charge separation that follow an absorption event. We apply the formalism to twolevel photovoltaic cells and our purpose is to pedagogically explain the main aspects of the model. The formalism is developed in the energy domain and provides detailed knowledge about existence or absence of localized states and their effects on electronic structure and photovoltaic yield.
T. Nemati Aram et al., published by EDP Sciences, 2017
This is an Open Access article distributed under the terms of the Creative Commons Attribution License
(http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Nanostructured photovoltaic cells including organic solar cells, dye and quantum dot sensitized solar cells have received great attention as promising photovoltaic technologies because of their prospective applications and interesting operation principle [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. In nanostructured solar cells, following the photon absorption, electronhole pair creation occurs in a confined zone [13, 14, 15, 16] and hence, the Coulomb interaction between the charge carriers and their possible recombination plays an essential role. Fundamentally different from conventional siliconbased photovoltaic devices, in nanostructured solar cells the charge carriers separation and transport rely on a high density of nanoscale interfaces.
While understanding the performance of nanostructured solar cells has been a central effort of the scientific community for many years [17, 18, 19, 20, 21], theoretical approaches still are needed to facilitate the understanding of electronhole interaction and recombination effects on the cell performance. In view of their nanometric size, the semiclassical approaches are expected to fail and fully quantum treatments are needed.
We develop a new quantum formalism, which is based on quantum scattering theory and in particular on the LippmannSchwinger equation [22, 23, 24]. A nonequilibrium wave function, which describes the state of the working cell, is computed and allows determining all parameters of the cell and its performance. As an example, we consider twolevel quantum model and show that the formalism provides a comprehensive framework to understand the effects of electronhole interaction and recombination on the performance of photovoltaic systems.
Fig. 1 The twolevel model in the permanent (left panel) and transitory (right panel) regime of illumination. (Left panel) A molecular photocell with one HOMO and one LUMO orbitals attached to the electrodes. The red line represents the electronhole interaction and recombination inside the molecule and the hopping integrals of electron and hole are denoted by C and J. E is the energy of absorbed photon. (Right panel) Schematic representation of a donoracceptor system. Charge evacuation leads are considered as semiinfinite chains. Here, J_{e} and J_{h} are the coupling energies between two adjacent sites in the electron and hole chains, respectively. Also, C_{e} and C_{h} represent the first coupling energies between the interface states and charge evacuation leads. 
2 Theoretical method
2.1 Twolevel photovoltaic systems
The basic idea of our methodology is described through the example of twolevel photovoltaic systems with the electronhole interaction and nonradiative recombination. The twolevel system characterized by the HOMO (highest occupied molecular orbital) and the LUMO (lowest unoccupied molecular orbital) can be in the permanent or in the transitory regime of illumination.
As an example, for the photovoltaic devices in the permanent regime of illumination, we consider the molecular photocells where the energy conversion process takes place in a single molecular donoracceptor complex attached to electrodes. In the absence of photons, the whole system is in the ground state with filled valence bands and empty conduction bands. Under the flux of photons, absorption by the molecule creates electron and hole in LUMO and HOMO, respectively. Both charge carriers interact via the Coulomb potential and can be recombined in the molecule or can be transferred to their respective channels where they produce permanent photovoltaic current (see Fig. 1left).
The BulkHetero Junction (BHJ) structure is an example of systems in the transitory regime. Indeed in the BHJ structure, following the absorption of a photon an exciton is created in the donor side of the cell at one negative time, and then it diffuses up to the interface (the twolevel system). This exciton arrives at the interface at time t = 0 and our aim is to evaluate the total charge injected in each contact due to the exciton dissociation after a sufficiently large time (see Fig. 1right).
2.2 The electronhole pair Hamiltonian
In both systems of Figure 1, we have to determine the possible quantum states of the photogenerated electronhole pair. The Hilbert space of these electronhole pairs can be mapped onto a square lattice (Fig. 2). For the PV devices in the permanent regime, x(y) represents the position of the electron (hole), in the molecule or in the attached leads. Site x = 0 (y = 0) corresponds to the LUMO (HOMO) orbital in the molecule and x> 0 (y> 0) represents the electron (hole) position in its respective lead. Therefore, site (x = 0, y = 0) is the electronhole pair position on the molecule that is simply the initial state just after the photon absorption, i.e., the excited state.
Furthermore, for the PV devices in the transitory regime, x(y) represents the position of the electron (hole), on the interface or in the attached leads. Site x = 0 (y = 0) corresponds to the LUMO (HOMO) orbital on the interface and x> 0 (y> 0) represents the electron (hole) position in its respective lead. Therefore, site (x = 0, y = 0) is the electronhole pair position on the DA interface. The energy difference between LUMO and HOMO levels on the absorber molecule (permanent regime) or on the interface (transient regime) is equal to Δ. The coupling matrix elements between the molecular or interface states and the possible evacuation leads is denoted by C. The hopping matrix element inside each evacuation lead is considered uniform and denoted by J. The onsite energies of the electron at site (x) and the hole at site (y) are assumed to be ε_{e}(x) and ε_{h}(y), respectively.
Fig. 2 The Hilbert space of the electronhole pair by considering just one evacuation channel for each charge carrier, with one state at each point (x, y) of the lattice. The coordinates x and y of a given state represent the position of electron and hole in their respective leads. ε(x,y) is the onsite energy of each site of the square lattice and the hopping integrals (C and J) are along the bonds of the square lattice. 
Fig. 3 Weight of the scattered wave function ( Ψ_{p}(E) ⟩) on different sites (x,y) of the square lattice normalized to its value on the initial site ( ⟨ x,y  Ψ_{P}(E) ⟩  ^{2}/  ⟨ 0,0  Ψ_{P}(E) ⟩  ^{2}) for various absorbed photon energies (E). The energy continuum (EC) is between 0.8 and 3.2 eV. 
The effective Hamiltonian of the system is of the tightbinding type additionally including the electronhole interaction term (1)where the first term indicates the total onsite energy of each square lattice basis state, which is defined as a summation over the electron onsite energy, the hole onsite energy and the Coulomb interaction energy between them: (2)I(x,y) represents the Coulombtype interaction between the charge carriers and is modeled by: (3)Since I(x,y) is an attractive Coulomb interaction, therefore U and V have negative values.
In equation (3),U represents the strength of the local electronhole interaction, i.e., when they are in the same place (either on the absorber molecule or on the DA interface) and V is the strength of the longrange electronhole interaction. In this formalism, the probability of photogenerated electronhole pair localrecombination inside the absorber molecule or on the DA interface is taken into account by adding an imaginary part − iΓ_{R}/ 2 to the onsite energy of the site (0, 0), where Γ_{R} is the recombination rate.
Finally, the second term in equation (1) represents the coupling energy between two adjacent basis states on the square lattice. As pointed above, the coupling parameters between molecular states or interface states (i.e., site (0, 0)) and their first neighbors are taken different form the other coupling energies. The realistic values of hopping and interaction energies could be obtained from ab initio calculations or from experiments, however in this study they are taken as parameters.
2.3 Nonequilibrium wave function of solar cell
Let us first discuss the formalism for the case in the permanent regime. We consider a photovoltaic cell as a system submitted to an incident flux of photons and assume that the whole system (PV cell and electromagnetic field) is in a stationary state. By applying quantum scattering theory, an open quantum system method, in particular the LippmannSchwinger equation, the photovoltaic system is described by a wave function. The incoming state of the theory Φ_{inc}⟩ represents the photon field with the PV cell in its ground state. By the dipolar interaction between the photovoltaic system and the electromagnetic field this incident state Φ_{inc}⟩ is coupled to a state where one photon is absorbed and one electronhole pair is created.
Based on the LippmannSchwinger equation and applying the Born approximation, the total wave function of the system is (4)The second term in the right hand side of the above equation,  Ψ_{p}(E) ⟩, is called the scattered wave function which represents the charge carriers photogenerated by absorption of a photon with energy E and plays an important role.  Ψ_{p}(E) ⟩ represents the whole system i.e., the cell with one electron hole pair plus the photon field with one less photon. Yet for the rest of this paper we focus on the current in the material system therefore we do not need to specify the part of  Ψ_{p}(E) ⟩ which represents the photon field. With this convention  Ψ_{p}(E) ⟩ is defined by the following equation: (5)where α^{2}(E) is the product of the square of the dipole matrix element of the molecular transition d times the electromagnetic energy density ρ(E) of photons with energy E, i.e., α^{2}(E) = d^{2}ρ(E) / 2ε_{0}. ε_{0} is the vacuum permittivity, ε is an infinitesimal positive energy and  0,0 ⟩ is the excited state after the photon absorption. Furthermore, H is the Hamiltonian of the electronhole pair, which is defined by the tightbinding model on the mentioned square lattice.
To have an idea about the behavior of scattered wave function  Ψ_{p}(E) ⟩, Figure 3 represents the weight of  Ψ_{p}(E) ⟩ on different sites (x,y) of the square lattice normalized to its value on the initial site  ⟨ x,y  Ψ_{P}(E) ⟩  ^{2}/  ⟨ 0,0  Ψ_{P}(E) ⟩  ^{2} by considering two different absorbed photon energies for a cell with energy continuum (EC) lying between 0.8 and 3.2 eV. It has to be noted that the energy continuum is the all allowed energies for the electronhole pair. As can be seen, for the photon energies outside the continuum, the scattered wave function is nearly localized on the site  0,0 ⟩ which is simply the absorber molecule. As the photon energy increases such that lies into the continuum, the scattered wave function extends on the various sites of the square lattice.
2.4 Fluxes
In this formalism, the cell performance can be described through the definition of a series of fluxes and knowledge of scattered wave function  Ψ_{p}(E) ⟩ enables one to compute all the essential fluxes. The main three fluxes are: (1) the flux of absorbed photons Φ_{Ph}(E) which is the number of absorbed photons per unit time; (2) the fluxes of electronhole pairs that recombine in the molecule Φ_{R}(E); and (3) the flux of pairs that escape from the molecule and results in the photovoltaic current Φ_{C}(E). The flux conservation implies that (6)The flux of absorbed photons is related to the local DOS n(E) through the Fremi’s golden rule (7)when the electron and the hole are not interacting n(E) is nothing but the joint density of states between electrons and holes which is a key quantity for lightmatter interaction. If the electron and the hole are interacting n(E) is defined as the local density of states of electronhole pair. This equation presents the direct relation between the photon absorption capability and the local density of states (LDOS) on site  0,0 ⟩ of the square lattice. n(E) can be expressed based on the diagonal matrix elements of the retarded Green’s function (8)As will be discussed in the numerical method section, can be expressed by: (9)where Σ_{0} and Γ_{R} represents the selfenergy of the system and recombination rate, respectively.
It has to be noted that the LDOS indicates all the possible energy states for the electronhole pair and implicitly presents the capability of the photovoltaic cell to absorb the light. Additionally, the current and recombination fluxes can be expressed by Finally, the fluxes that allow describing the cell performance depend on recombination rate Γ_{R} and selfenergy Σ_{0}. The imaginary part of the selfenergy Σ_{0} plays a central role for the injection of current as shown by the above equation. The real part of the selfenergy Σ_{0} plays a central role to determine the position of the resonance in the density of states n(E) and therefore in the frequency of photons that are most efficiently absorbed by the twolevel system. The selfenergy can be determined from Hamiltonian of the square lattice (see numerical method).
Fig. 4 Propagation of different orders of recursion wavevectors. The coupling parameters are J_{e} = J_{h} = 0.3 eV, C_{e} = 0.2 eV and C_{h} = 0.15 eV. 
2.5 Expression of the quantum yield
Let us first consider a photovoltaic cell in the permanent regime of illumination. At a given photon energy E, the yield Y(E) is proportional to the ratio of photogenerated electrons or holes that arrive at the electrodes and the total number of absorbed photons at this given energy (12)Through the assumption that α^{2}(E) = α^{2} in the region where n(E) (i.e., photon absorption) is important, the average yield or in other words, the charge separation yield Y, which is the proportion of the all electronhole pairs, generated by different photons and giving rise to the photovoltaic current can be defined as (13)Here, we draw the reader’s attention to an interesting point: in the limit of small recombination rate, if photon energy E lies into the energy continuum then there is a transformation into an electronhole pair and therefore yield Y(E) is one. Whereas, if E corresponds to a bound electronhole state, no current can be injected in the leads and finally the electronhole pair recombine so that Y(E) is zero. As a consequence, the average yield Y given by equation (12) can be expressed by Y = 1 − P, with P representing the total weight of localized states. This means that the efficiency of the cell depends critically on the existence or absence of the bound states as illustrated bellow. For the PV devices in the transitory regime, we consider the total charges injected in a bind (in units of the electron charge) Q(14)where Â is the current operator on that bind and (15)where  ψ ⟩ is the electronhole pair wave function and U(t) is the time evolution operator.
Through introducing the total electron number Q_{e} or hole number Q_{h} injected in the contacts one obtains (16)Interestingly, through the Fourier transform between the time and energy domain one obtains an expression for the charge separation yield in the transitory regime the same as equation (13).
3 Numerical method
3.1 Recursion method
In this section, we discuss the recursion method, which presents a central role in the computation of the scattered state and fluxes for a given Hamiltonian. An initialstate vector, on whose projection one going to compute a density of states, is selected. One can build a new basis set and a series of recursion coefficients starting from the mentioned initial state. In this new basis, the Hamiltonian becomes tridiagonal and the recursion coefficients are the matrix elements of the Hamiltonian [25, 26, 27, 28].
Fig. 5 LDOS as a function of the energy of the absorbed photon under the influence of shortrange electronhole interaction. 
We consider a system described by a tightbinding Hamiltonian H. To a given normalized state  ψ_{0} ⟩ we can always associate a recursion basis which is constructed by a Schmidt orthogonalization procedure starting from the set of states  ψ_{0} ⟩, Hψ_{0}⟩, H^{2} ψ_{0}⟩, ..., H^{n − 1}ψ_{0}⟩, ...Let us consider Hψ_{0}⟩ and decompose it into two components: first component parallel to ψ_{0}⟩ and the second one orthogonal to ψ_{0}⟩. Therefore, we can write (17)In the next step, we consider Hψ_{1}⟩ which can be decomposed into a component parallel to the space spanned by ψ_{0}⟩, ψ_{1}⟩ and a component orthogonal to this space. We obtain (18)The process can be repeated and leads to the construction of a set of states ψ_{n}⟩ which are orthonormal and satisfy the below equation (19)The a_{n} and b_{n} parameters are called recursion coefficients and ψ_{n}⟩ are the recursion wave vectors. At any step n, the recursion coefficient a_{n} is determined by ⟨ ψ_{n}  H  ψ_{n} ⟩ which is a real number. Then, the recursion coefficient b_{n} and recursion wave vector ψ_{n + 1}⟩ are obtained by the normalization condition for ψ_{n + 1}⟩ and the choice of a real positive b_{n}. An important property of the states ψ_{n}⟩ is that they spread progressively from an initial state. In Figure 4, the propagation of recursion wave vectors on a square lattice is shown.
The significant achievement of the recursion method is that the Hamiltonian H on the basis (ψ_{0}⟩, ψ_{1}⟩, ...) can be written in a tridiagonal form (20)Therefore, can be written in the following form (21)and in a shortened notation (22)where, Σ_{0}(E) is the selfenergy as pointed previously.
Finally, let us recall that the recursion method allows the computation of the scattered wave function and of all fluxes. Technically, one can compute  Ψ_{p}(E) ⟩ based on the product of recursion coefficients.
4 Results and discussions
4.1 Study of occurrence of bound states under the influence of shortrange electronhole interaction
In order to analyze the cell performance we consider the spectral information and in particular, the local DOS on site (0,0). For the numerical calculation, parameters Δ = 2eV, J_{e} = J_{h} = 0.3eV, C_{e} = C_{h} = 0.3eV have been used. With the chosen parameters the electronhole pair energy continuum lies between 0.8 and 3.2 eV.
In Figure 5, the LDOS under the influence of shortrange Coulomb interaction conditions is plotted as a function of the absorbed photon energy. As can be seen for small values of  U , there is a single peak which becomes narrower as  U  increases and all weight of the DOS is in the continuum part (left panel). As  U  increases, the DOS lineshape splits into two parts, a part inside the energy continuum and a sharp peak outside the energy continuum (right panel).
The narrow peak outside the continuum is called excitonic state, which suppresses the charge carrier injection to the energy continuum and expedite the charge carriers recombination. The creation of excitonic state can be understood in full details through the discussion given in reference [21], which shows that there is a critical interaction energy to create localized states. Mathematically an excitonic peak appears at a given energy E outside the energy continuum, if (23)Indeed excitonic states are the poles of the retarded Gree’s function. It is expected that in a boundstate the charge carriers will ultimately recombine since they cannot escape into the leads. The appearance of excitonic states decreases the number of charge carriers arrived at the electrodes and therefore lowers the photocell yield.
For the local interacting electronhole pair, the dependence of the yield Y on the electronhole interaction U for two different recombination rate (Γ_{R}) and different set of coupling parameters (C = C_{e} = C_{h}) is investigated in Figure 6. As can be seen, in all cases, for small values of interaction energy, the yield remains 1 for Γ_{R} = 0. As Uincrease, the yield decreases. This behaviour can be understood based on the information provided in Figure 5. For larger values of  U , the charge carriers will stay on the molecule (or interface) to form a localized state because their energy does not lie in the energy continuum of the contacts. Furthermore, for large values of the coupling parameters (C), charge carriers will transfer more efficiently to the evacuation channels and hence the cell remains efficient over a wider range of the recombination parameter Γ_{R}.
Fig. 6 Yield of a twolevel photovoltaic cell as a function of interaction energy (U) for different recombination rate (Γ_{R}) and coupling parameter (C). 
4.2 Occurrence of bound states under the influence of longrange Coulomb interaction
Now we turn to the case of the longrange electronhole interaction and examine its impact on the spectral properties (see Fig. 7). We find that under the influence of longrange electronhole interaction a series of excitonic peaks appears below the lower band of energy continuum. This is expected, as it is known that the longrange Coulomb attraction creates localized states. In this figure, all the excitonic peaks close to lower band edge cannot be resolved.
We see that, as the interaction strength increases, the weight of localized states increases as well. Therefore, under the influence of longrange interaction less charge carriers exit through the contacts because of localizedstate formation and hence it is expected that the yield decreases.
Fig. 7 LDOS as a function of the energy of the absorbed photon under the influence of longrange electronhole interaction. 
4.3 Occurrence of bound states under the influence of lattice distortion and polaron formation
The charge carrier mobility can be influenced by the electronphonon vibration interaction and therefore in the strong coupling regime, the charge carrier transport can be described by a polaron hopping from molecule to molecule where a polaron is a freely moving charge (electron or hole), surrounded by virtual cloud of phonons [29].
To have a clear idea, suppose in a DA system, the exciton arrives at the interface at time t = 0 coming from the donor side. After the exciton dissociation at the DA interface, the electron either will be recombined with the hole or will move through a set of acceptor sites. To investigate the effects of coupling to the phonon modes, we suppose that electron on each acceptor site can be coupled to one single phonon mode. The model considered here is exactly that described in [29].
Fig. 8 Electronic structure as a function of the energy under the influence of coupling to the phonon modes. (Left panel) Far from the DA interface. (Right panel) On the DA interface. The hopping integral on the linear chain is J = 1 eV [29]. 
Figure 8 represents the typical electronic structure of such a system far from interface (left panel) and on the interface (right panel). Based on the Holstein Hamiltonian model β represents the strength of coupling to the phonon modes and ε is the LUMOLUMO offset [29]. As can be seen from the left panel, the electronic structure is composed of a series of bands called polaronic bands and an energy gap separates the different bands. Based on the right panel, for a given ε the electronic structure may contain the energy states on the allowed polaronic bands as well localized states in the energy gap. The charge carriers lying in a polaronic band evacuate and arrive at the electrodes. On the other hand, the charge carriers localized in the bound state in the gap recombine quickly and cannot lead to photovoltaic current
5 Conclusions
To conclude, we have presented a nonequilibrium quantum formalism, which allows describing the fundamental phenomena related to charge separation or charge recombination after the absorption of a photon. This formalism should be useful in particular for nanostructured photovoltaic cells.
In this paper, we considered the simplest case of a twolevel quantum system This simple study illustrates the fact that localized excitonic or polaronic states can play a central role in the limitation of the quantum yield. This study represents that the formalism can be applied to include in particular the role of electronhole interaction and the role of phonon on exciton propagation and charge separation.
Acknowledgments
The authors would like to thank Simone Fratini and KevinDavis Richler for fruitful discussions.
References
 C.J. Brabec, Sol. Energ. Mat. Sol. Cells 83, 273 (2004) [CrossRef] (In the text)
 C. Lungenschmied, G. Dennler, H. Neugebauer, S.N. Sariciftci, M. Glatthaar, T. Meyer, A. Meyer, Sol. Energ. Mat. Sol. Cells 91, 379 (2007) [CrossRef] (In the text)
 M.T. Lloyd, J.E. Anthony, G.G. Malliaras, Mater. Today 10, 34 (2007) [CrossRef] (In the text)
 M. Niggemann, B. Zimmermann, J. Haschke, M. Glatthaar, A. Gombert, Thin Solid Films 516, 7181 (2008) [CrossRef] (In the text)
 J. Roncali, Acc. Chem. Res. 42, 1719 (2009) [CrossRef] [PubMed] (In the text)
 F.C. Krebs, J. Fyenbo, M. Jørgensen, J. Mater. Chem. 20, 8994 (2010) [CrossRef] (In the text)
 Y. Sun, G.C. Welch, W.L. Leong, C.J. Takacs, G.C. Bazan, A.J. Heeger, Nat. Mater. 11, 44 (2012) [CrossRef] (In the text)
 Y. Li, Acc. Chem. Res. 45, 723 (2012) [CrossRef] [PubMed] (In the text)
 A.J. Nozik, Physica E 14, 115 (2002) [CrossRef] (In the text)
 A.J. Labelleet, Nano Lett. 15, 1101 (2015) [CrossRef] [PubMed] (In the text)
 I. Chung, Nature 485, 486 (2012) [CrossRef] [PubMed] (In the text)
 M. Liang, C. Jun, Chem. Soc. Rev. 42, 3453 (2013) [CrossRef] [PubMed] (In the text)
 D. Cahen, G. Hodes, M. Gratzel, J.F. Guillemoles, I. Riess, J. Phys. Chem. B 104, 2053 (2000) [CrossRef] (In the text)
 J.L. Brédas, J. E. Norton, J. Cornil, V. Coropceanu, Acc. Chem. Res. 42, 1691 (2009) [CrossRef] [PubMed] (In the text)
 S. Ten Cate et al., J. Phys. Chem. Lett. 4, 1766 (2013) [CrossRef] (In the text)
 F. Yang, S.R. Forrest, ACS Nano 2, 1022 (2008) [CrossRef] [PubMed] (In the text)
 S. Ajisaka, B. Žunkovič, Y. Dubi, Sci. Rep. 5, 8312 (2015) [CrossRef] [PubMed] (In the text)
 M. Einax, M. Dierl, A. Nitzan, J. Chem. Phys. 115, 21396 (2011) (In the text)
 M. Einax, A. Nitzan, J. Chem. Phys. 145, 014108 (2016) [CrossRef] (In the text)
 M. Ernzerhof, M.A. Bélanger, D. Mayou, T. Nemati Aram, J. Chem. Phys. 144, 134102 (2016) [CrossRef] (In the text)
 T. Nemati Aram, A. Asgari, D. Mayou, Europhys. Lett. 115, 18003 (2016) [CrossRef] [EDP Sciences] (In the text)
 J.J. Sakurai, Modern Quantum Mechanics (Addison Wesley, 1994) (In the text)
 T. Nemati Aram, P. AnghelVasilescu, A. Asgari, M. Ernzerhof, D. Mayou, J. Chem. Phys. 145, 124116 (2016) [CrossRef] (In the text)
 T. Nemati Aram, M. Ernzerhof, A. Asgari, D. Mayou, J. Chem. Phys. 146, 034103 (2017) [CrossRef] (In the text)
 F. Triozon, J. Vidal, R. Mosseri, D. Mayou, Phys. Rev. B 65, 220202 (2002) [CrossRef] (In the text)
 P. Darancet, V. Olevano, D. Mayou, Phys. Rev. B 81, 155422 (2010) [CrossRef] (In the text)
 D. Mayou, Europhys. Lett. 6, 549 (1988) [CrossRef] (In the text)
 S. Roche, D. Mayou, Phys. Rev. B 60, 322 (1999) [CrossRef] (In the text)
 S. Bera, N. Gheeraert, S. Fratini, S. Ciuchi, S. Florens, Phys. Rev. B 91, 041107 (2015) [CrossRef] (In the text)
Cite this article as: Tahereh Nemati Aram, Asghar Asgari, Matthias Ernzerhof, Pascal Quémerais, Didier Mayou, Quantum modeling of twolevel photovoltaic systems, EPJ Photovoltaics 8, 85503 (2017).
All Figures
Fig. 1 The twolevel model in the permanent (left panel) and transitory (right panel) regime of illumination. (Left panel) A molecular photocell with one HOMO and one LUMO orbitals attached to the electrodes. The red line represents the electronhole interaction and recombination inside the molecule and the hopping integrals of electron and hole are denoted by C and J. E is the energy of absorbed photon. (Right panel) Schematic representation of a donoracceptor system. Charge evacuation leads are considered as semiinfinite chains. Here, J_{e} and J_{h} are the coupling energies between two adjacent sites in the electron and hole chains, respectively. Also, C_{e} and C_{h} represent the first coupling energies between the interface states and charge evacuation leads. 

In the text 
Fig. 2 The Hilbert space of the electronhole pair by considering just one evacuation channel for each charge carrier, with one state at each point (x, y) of the lattice. The coordinates x and y of a given state represent the position of electron and hole in their respective leads. ε(x,y) is the onsite energy of each site of the square lattice and the hopping integrals (C and J) are along the bonds of the square lattice. 

In the text 
Fig. 3 Weight of the scattered wave function ( Ψ_{p}(E) ⟩) on different sites (x,y) of the square lattice normalized to its value on the initial site ( ⟨ x,y  Ψ_{P}(E) ⟩  ^{2}/  ⟨ 0,0  Ψ_{P}(E) ⟩  ^{2}) for various absorbed photon energies (E). The energy continuum (EC) is between 0.8 and 3.2 eV. 

In the text 
Fig. 4 Propagation of different orders of recursion wavevectors. The coupling parameters are J_{e} = J_{h} = 0.3 eV, C_{e} = 0.2 eV and C_{h} = 0.15 eV. 

In the text 
Fig. 5 LDOS as a function of the energy of the absorbed photon under the influence of shortrange electronhole interaction. 

In the text 
Fig. 6 Yield of a twolevel photovoltaic cell as a function of interaction energy (U) for different recombination rate (Γ_{R}) and coupling parameter (C). 

In the text 
Fig. 7 LDOS as a function of the energy of the absorbed photon under the influence of longrange electronhole interaction. 

In the text 
Fig. 8 Electronic structure as a function of the energy under the influence of coupling to the phonon modes. (Left panel) Far from the DA interface. (Right panel) On the DA interface. The hopping integral on the linear chain is J = 1 eV [29]. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.