Numerical drift-diffusion analysis of current flow and photoelectric effect in metal halide perovskite-silicon tandem solar cell considering temperature, structure and doping effects for device optimization.
Juan Alberto Gonzalez Cuevas a, Luis Bernal a, Giovanni Secchia a, Juan De Egea a, Fernando Pio Barrios a
a Facultad de ingeniería, Universidad Nacional de Asunción, Paraguay
Proceedings of MATSUS Spring 2026 Conference (MATSUSSpring26)
A5 From halide perovskites to perovskite-inspired materials – Synthesis, Modelling and Application
Barcelona, Spain, 2026 March 23rd - 27th
Organizers: Gustavo de Miguel, Lorenzo Malavasi and Isabella Poli
Oral, Juan Alberto Gonzalez Cuevas, presentation 618
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

© FUNDACIO DE LA COMUNITAT VALENCIANA SCITO
We use our own and third party cookies for analysing and measuring usage of our website to improve our services. If you continue browsing, we consider accepting its use. You can check our Cookies Policy in which you will also find how to configure your web browser for the use of cookies. More info