I am using Athena to create a p-n junction diode device, and model the p and n doping such that we have p+ and n+ near the electrodes, then pwell and nwell below them. For a dark current simulation I can only bias the structure to 0.0001V and I see an electric field of ~3e6 V/cm in the aluminum electrode near the silicon p+ and n+ layers.
The breakdown should be much higher and I do not know why there is such a high E-field in the aluminum electrodes, when it should be zero.