Modelling of hydrogen transport in silicon solar cell structures under equilibrium conditions

This paper presents a model for the introduction and redistribution of hydrogen in silicon solar cells at temperatures between 300 and 700 °C based on a second order backwards difference formula evaluated using a single Newton-Raphson iteration. It includes the transport of hydrogen and interactions with impurities such as ionised dopants. The simulations lead to three primary conclusions: (1) hydrogen transport across an n-type emitter is heavily temperature dependent; (2) under equilibrium conditions, hydrogen is largely driven by its charged species, with the switch from a dominance of negatively charged hydrogen (H−) to positively charged hydrogen (H+) within the emitter region critical to significant transport across the junction; and (3) hydrogen transport across n-type emitters is critically dependent upon the doping profile within the emitter, and, in particular, the peak doping concentration. It is also observed that during thermal processes after an initial high temperature step, hydrogen preferentially migrates to the surface of a phosphorous doped emitter, drawing hydrogen out of the p-type bulk. This may play a role in several effects observed during post-firing anneals in relation to the passivation of recombination active defects and even the elimination of hydrogen-related defects in the bulk of silicon solar cells.