A Stable, High-Order Time-Stepping Scheme for the Drift-Diffusion Model in Modern Solar Cell Simulation

Abstract

This paper presents a one-dimensional transient drift-diffusion simulator for advanced solar cells, integrating a structure-preserving finite-volume spatial discretization with Scharfetter-Gummel-type fluxes and a high-order, L-stable implicit Runge-Kutta (Radau IIA) temporal integrator. The scheme ensures local charge conservation, handles sharp material interfaces, achieves second-order spatial and fifth-order temporal convergence, and exhibits excellent computational efficiency and numerical stability. Its accuracy is verified against the classical depletion approximation in p-n junction and validated against four independent software implementations for an organic photovoltaic device under a common protocol. The framework's extensibility is demonstrated by incorporating exciton kinetics in organic solar cells, capturing multi-timescale dynamics, and by modeling mobile ions in perovskite solar cells, reproducing characteristic J-V hysteresis without empirical parameters. This work provides a robust, high-order numerical foundation for simulating coupled charge, exciton, and ion transport in next-generation photovoltaic devices.

Abstract: arXiv descriptive metadata (CC0). Thank you to arXiv for use of its open access interoperability. This service was not reviewed or approved by, nor does it necessarily express or reflect the policies or opinions of, arXiv.