Drift-diffusion model

Hieoretical works have focused on morphology predictions in order to address the property charaaeristics of photovoltaic devices consisting of rod-coil block copolymers. For this purpose, a continuum formalism termed the drift diffusion model has been adopted. This characterizes the transport of excitons, electrons, and holes under photoadsorption to render predictions regarding the device characteristics. ... 

Figure 14.13 Response of the short-circuit photocurrent of P3HT F8TBT solar cells to a square pulse illumination of varying intensity. Turn-on of the source is at 0 is and turn-off at 100 is. (a) The experimental photocurrent transients (b) the photocurrent transients simulated using a time-dependent drift diffusion model...

The capability of drift-diffusion models can be increased by introducing localized states into the band gap. Recently, several studies have introduced single trap levels [83] as well as distributions of localized states in order to describe the results of transient and steady state experiments on polymer fullerene solar cells [54, 84-89]. Most of these models use a Shockley-Read-Hall type occupation statistics for the localized states, which we will discuss in more detail in Sect. 2.3 and the Appendix 2 before discussing some of the implications of this model in the case studies in Sect. 4.1. 

To implement the drift-diffusion models for organic photovoltaic devices, we need a way to determine the parameters describing transport, generation and recombination and their functional form. It is not yet possible to predict these parameters from the chemical structure of the materials they should therefore be determined experimentally for the materials under investigation. Here, we give a brief overview of the experimental methods that are used to determine input parameters for optical or electrical modelling and to validate the model. 

Fig. 20 Determination of the coefficient of bimolecular recombination by performing TDCF experiments with variable delay between the excitation pulse and application of the collection bias, (a) Scheme of the experiment, (b) Experimental TDCF photocurrent transients (open squares) measured on a 200 nm thick layer of slow-dried POFTTiPCBM (1 1) during application of different collection biases Fcoii- The collection bias was applied 150 ns after the laser pulse (t = 0 in this graph). Solid lines show fits to the data using a numerical drift diffusion model with constant electron and hole mobilities. A noteworthy observation is that charges can be fully extracted from these layers within a few hundreds of nanoseconds for a sufficiently high collection bias [171]. (c-f) Q-pre, 2coii> and 2,o, plotted as a function of the delay time for as-prepared and thermally annealed chloroform-cast P3F1T PCBM, and with the pre-bias Fpre set either to 0.55 V (near open circuit) or to 0 V (short-circuit conditions) [172]. Solid lines show fits with an iterative model that considers bimolecular recombination of free charges in competition with their extraction...

The high aspect ratio of nanorods can facilitate charge transport, while the handgap can he tuned by vaiying the nanorod radius. This enables the absorption spectmm of the devices to be tailored to overlap with the solar emission spectmm, whereas traditionally polymer absorption has been limited to only a small fraction of the incident solar irradiation. At present, the nanorods in polymer solar cells are typically incorporated into a homopolymer matrix. An alternative to this approach is to incorporate the nanorods into either a polymer blend or diblock copolymer system. The photovoltaic properties of nanorod polymer composites could potentially be improved due to the percolation of nanorods, and the presence of continual electrical pathways, from the DA interfaces to the electrodes. To test this hypothesis, we use the distribution of nanorods from the self-assembled stmcture in Figure 1(b) as the input into a drift-diffusion model of polymer photovoltaics. 

We now couple the GH-BD model of nanorod polymer composites with the above photovoltaic model. In particular, we feed the morphology from Figure 1(b) into the drift-diffusion model and elucidate the photovoltaic properties of this self-assembled stmcture. Figure 4 shows the exciton, electron, and hole concentrations in the device. 

The device model describes transport in the organic device by the time-dependent continuity equation, with a drift-diffusion form for the current density, coupled to Poisson s equation. To be specific, consider single-carrier structures with holes as the dominant carrier type. In this case,... 

In the fluid model the momentum balance is replaced by the drift-diffusion approximation, where the particle flux F consists of a diffusion term (caused by density gradients) and a drift term (caused by the electric field ) ... 

In order to be able to explain the observed results plasma modeling was applied. A one-dimensional fluid model was used, which solves the particle balances for both the charged and neutral species, using the drift-diffusion approximation for the particle fluxes, the Poisson equation for the electric field, and the energy balance for the electrons [191] (see also Section 1.4.1). 

In reality, the slip velocity may not be neglected (except perhaps in a microgravity environment). A drift flux model has therefore been introduced (Zuber and Findlay, 1965) which is an improvement of the homogeneous model. In the drift flux model for one-dimensional two-phase flow, equations of continuity, momentum, and energy are written for the mixture (in three equations). In addition, another continuity equation for one phase is also written, usually for the gas phase. To allow a slip velocity to take place between the two phases, a drift velocity, uGJ, or a diffusion velocity, uGM (gas velocity relative to the velocity of center of mass), is defined as... 

This is the traditional diffusion model given in Eq. (6.305) with the diffusion coefficient D proportional to 1/er. For this case, the so-called spurious drift term vanishes because the effects of a and a cancel each other out in the stationary state. The stationary distribution is proportional to the Boltzmann distribution exp(-I 7/c7 ) and independent of D. 2. a = constant. Then, Eq. (6.308) becomes... 

To determine the dispersed phase velocities as occurring in the phasic continuity equations in both formulations, the momentum equation of the dispersed phases are usually approximated by algebraic equations. Depending on the concept used to relate the phase k velocity to the mixture velocity the extended mixture model formulations are referred to as the algebraic slip-, diffusion- or drift flux models. 

Drift-diffusion for ions The drift-diffusion approximation (Eq. 25) is made for ions, e.g., ion inertia is neglected. This appears to be a good approximation for Kn < 0.2. Eq. (37) provides the field that drives ions in the plasma. Also, a constant mobility is assumed for ions. Actually, when collisions are infrequent (high Kn number), it is better to use a variable mobility model, for which the ion mobility is inversely proportional to the magnitude of the directed velocity of the ion [6]. This model results by assuming a constant mean-free-path of ion-neutral collisions. The constant mobility model results by assuming a constant mean-free-time for ion-neutral collisions. 

Lastly, the walls of the drifts are subjected to a convection condition with air circulating at a temperature T2 of 20°C in the lower drift and Tj of 60°C in the upper. From the geometry of each drift and air flow circulating in each well, the exchange coefficients hi and h2 are respectively about 11 W/mVK for the upper drift and 8 WW/K for the lower. Using these boundary conditions and a diffusion model taking of account the conduction of heat in the rock mass, the temperature field has been calculated by CAST3M (see fig.4). 

This observation is fundamental, reveling the fact that the characteristic energies of the Fokker-Planck equation are reactive energies, and in the final, nonequilibrium energies. This aspect is directly correlated with the nonequilibrium character specific for the Fokker-Planck equation while modeling open systems (driven by drift diffusion and factors, stochastic noise, etc.). Moreover, if the analytical solution of the eigen-values for the Schrodinger equation with the potential ) is considered, the consecrated expression is obtained ... 

---