Publication date: 15th December 2025
A drift-diffusion model is developed to calculate the spectral response, current-voltage characteristics and conversion efficiency of a thin film metal halide perovskite-silicon tandem solar cell. The model is based on solving carrier continuity partial differential equations to obtain the spatial distribution and time evolution of electron and hole densities, coupled with Poisson’s equation in order to get the electric field and potential influence on carrier drift considering the doped charges in layers of the device.
The photogeneration of carriers is dependent upon the absorption coefficient and normal incident reflectivity. Drift and diffusion currents correspond to the flow of carriers due to the electric field and concentration gradient, which depend on carrier mobilities, recombination mechanisms both at the surface and within the device [1-3], and absorption-reflection coefficients derived from the dielectric function of the perovskite and silicon materials [4-6], with temperature effects arising from the variation of energy levels due to lattice expansion and contraction [7].
The structure of the modeled two-terminal solar cell consists of spiro-OMeTAD/CH3NH4PbI3/SnO2 layers in tandem over a p++/n/ n++ silicon cell. Ohmic contacts are assumed at the front and back of the device with corresponding boundary conditions for carrier densities and potential. Results have been compared to similar perovskite-silicon tandem solar cells with high conversion efficiency [8-11].
The equations were numerically solved employing both finite difference and finite element discretization methods using MATLAB and COMSOL Multiphysics software.
For performance enhancement, the solar cell quantum efficiency, fill factor, conversion efficiency and open circuit voltage have been evaluated for different values of structural parameters such as layer depth and doping, as well as operating variables such as temperature and photon energy, in order to decrease charge accumulation and carrier recombination and approach the Shockley-Queisser efficiency limit.
The authors acknowledge the financial support from CONACYT Paraguay Consejo Nacional de Ciencia y Tecnología grant number PINV01-623
