Analytical and Numerical Study of Photocurrent Transients in Organic Polymer Solar Cells

Carlo de Falco,Riccardo Sacco,Maurizio Verri

DOI: https://doi.org/10.1016/j.cma.2010.01.018

2012-06-28

Abstract:This article is an attempt to provide a self consistent picture, including existence analysis and numerical solution algorithms, of the mathematical problems arising from modeling photocurrent transients in Organic-polymer Solar Cells (OSCs). The mathematical model for OSCs consists of a system of nonlinear diffusion-reaction partial differential equations (PDEs) with electrostatic convection, coupled to a kinetic ordinary differential equation (ODE). We propose a suitable reformulation of the model that allows us to prove the existence of a solution in both stationary and transient conditions and to better highlight the role of exciton dynamics in determining the device turn-on time. For the numerical treatment of the problem, we carry out a temporal semi-discretization using an implicit adaptive method, and the resulting sequence of differential subproblems is linearized using the Newton-Raphson method with inexact evaluation of the Jacobian. Then, we use exponentially fitted finite elements for the spatial discretization, and we carry out a thorough validation of the computational model by extensively investigating the impact of the model parameters on photocurrent transient times.

Numerical Analysis,Mathematical Physics,Analysis of PDEs,Computational Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the mathematical modeling and numerical simulation of the photocurrent transient behavior in organic polymer solar cells (OSCs). Specifically, the authors hope to provide a self - consistent framework, including:

1. **Existence analysis**: Prove the existence of solutions of the proposed mathematical model under steady - state and transient conditions.
2. **Numerical solution algorithm**: Develop effective numerical methods to solve these complex nonlinear partial differential equation systems.

### Mathematical model

The mathematical model of OSCs consists of a set of nonlinear diffusion - reaction partial differential equations (PDEs), which describe the transport of charge carriers (electrons and holes) and are coupled with a kinetic ordinary differential equation (ODE) that describes the time evolution of the exciton concentration. The key equations of the model are as follows:

- **Continuity equation**:
\[
\begin{cases}
\frac{\partial n}{\partial t}-\nabla\cdot J_{n}=G_{n}-R_{nn}\\
\frac{\partial p}{\partial t}-\nabla\cdot J_{p}=G_{p}-R_{pp}
\end{cases}
\]
where \(n\) and \(p\) represent the densities of electrons and holes respectively, \(J_{n}\) and \(J_{p}\) are the corresponding current densities, \(G_{n}\) and \(G_{p}\) are the carrier generation rates, and \(R_{nn}\) and \(R_{pp}\) are the recombination rates.

- **Current density expression**:
\[
\begin{cases}
J_{n}=D_{n}\nabla n-\mu_{n}n\nabla\phi\\
J_{p}=D_{p}\nabla p+\mu_{p}p\nabla\phi
\end{cases}
\]
where \(D_{n}\) and \(D_{p}\) are the diffusion coefficients, \(\mu_{n}\) and \(\mu_{p}\) are the mobilities, and \(\phi\) is the electrostatic potential.

- **Poisson equation**:
\[
-\nabla\cdot(\varepsilon\nabla\phi)=q(p - n)
\]
where \(\varepsilon\) is the dielectric constant and \(q\) is the magnitude of the electron charge.

- **Exciton kinetic equation**:
\[
\frac{\partial X}{\partial t}=g - r
\]
where \(X\) is the volume density of geminate pairs, \(g\) is the generation rate, and \(r\) is the recombination rate.

### Existence analysis

To simplify the problem and prove the existence of solutions, the authors introduced some assumptions (H1 - H4) and proved the existence of solutions under steady - state and transient conditions under these assumptions.

### Numerical method

The authors adopted an implicit adaptive time - discretization method for time discretization and used the Newton - Raphson method to linearize the nonlinear problem. For spatial discretization, the exponential fitting finite element method was used to ensure stable approximation of the internal and boundary layers.

### Research significance

Through these works, the authors aim to fill the gap between experimental research and computational models and provide theoretical support for optimizing the design and performance of organic polymer solar cells. In addition, they also pay special attention to the influence of exciton dynamics on the device turn - on time, which is of great significance for understanding material properties.