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

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.

💡 Research Summary

This paper presents a comprehensive mathematical and computational study of photocurrent transients in organic polymer solar cells (OSCs). The authors formulate a coupled system consisting of nonlinear diffusion‑reaction partial differential equations (PDEs) for charge carriers (electrons and holes) with electrostatic convection, and an ordinary differential equation (ODE) describing exciton dynamics. The exciton ODE provides the generation term for the charge‑carrier PDEs and includes exciton diffusion, separation, and recombination processes, thereby creating a feedback loop between exciton population and charge transport.

To establish a rigorous foundation, the model is reformulated in a variational framework. By selecting appropriate Sobolev spaces for potentials and carrier densities, the authors prove existence (and under certain monotonicity conditions, uniqueness) of solutions for both steady‑state and transient regimes. The proof exploits the monotone operator theory and maximum principle, handling the nonlinear source terms that couple the PDE and ODE.

For numerical simulation, an adaptive implicit time‑discretization is employed. The time step is automatically refined during rapid current rise, ensuring stability without sacrificing efficiency. At each time level, the nonlinear system is linearized via a Newton‑Raphson scheme; however, the Jacobian is evaluated inexactly (e.g., using lagged approximations or finite‑difference estimates) to reduce computational cost while preserving quadratic convergence in practice. Spatial discretization uses exponentially fitted finite elements, which accurately capture steep gradients in electric potential and carrier concentrations that arise near contacts and depletion regions. This fitting mitigates spurious oscillations typical of standard Galerkin methods in convection‑dominated regimes.

Extensive simulations validate the model against experimental transient data. A parametric sensitivity analysis investigates the impact of exciton diffusion coefficient, separation rate, recombination rate, carrier mobilities, and boundary conditions on two key performance metrics: turn‑on time (the time required for the photocurrent to reach a specified fraction of its steady‑state value) and peak current magnitude. Results show that the exciton separation rate exerts the strongest influence on turn‑on time; higher separation accelerates the rise of photocurrent dramatically. Carrier mobility primarily affects the slope of the transient and the steady‑state current level, while recombination rates modulate the peak amplitude. The analysis provides clear guidance for material and device engineering: enhancing exciton dissociation pathways and optimizing carrier transport properties are crucial for fast response and high efficiency.

Finally, the authors discuss the broader applicability of their framework to other organic electronic devices such as OLEDs and photodetectors, and outline future work involving multilayer architectures, non‑ideal contact models, and experimental validation of the predicted parameter trends.