Hi Rod,
Thanks for the quick response! We suspect that the breakdown in PSC's is caused by holes tunneling from the ETL into the perovskite's valence band at sufficient -V's (thanks to the mobile ions screening the field at the perovskite interfaces!). I'm interested in modeling this behavior of tunneling at the ETL/PVSK interface that causes current to flow. This type of modeling would be huge to accomplish.
I've lately been using the oghmanano software to look at band diagrams and calculate electric fields at the interface for now since I can't seem to get a PSC to "break down" in reverse bias (for now...

).
There are obviously other electrochemical processes occurring in the perovskite under reverse bias, so i dont expect the software to be able to accurately capture this behavior. There are also likely thermal breakdowns occurring in our real cells and it would be incredible to use this software to look at the Joule heating across a pn diode under reverse bias (is this already possible?).
Basically, modeling the tunneling of holes from the ETL into the perovskite would be amazing if i can see a PSC breakdown and pass current in reverse bias. My second goal is to study the breakdown of normal PN diodes via Zener tunneling at sufficient doping densities or show that the electric field strengths are strong enough to cause Joule heating/dielectric breakdown. I'm struggling to get current to flow in revB with the basic pn diodes for now, and I'd like to know your thoughts on whether modeling this breakdown is possible! I can provide a DOI on a paper for the breakdown of PSC's if that would help guide your thought process!
Thanks