Coupling Semiconductor and Heat Transfer

Greetings!

First some information on my model. I have a model of 1D PN diode. I want to measure 3 features:

1. Lattice temperatureincreasement due to forward current.

2. Impact Ionization Generation Rate

3. SRH Recombination Rate

To simulate this I've included Semiconductor and Heat Transfer in Solids modules. I've attached 3 snapshots of all included models and boundries. So here is what I mean by two physics coupling:

Cases 1 - One way couplying semi--->ht

In Semiconductor Material Model, Model Input Temperature is set to constant - 300[K] Heat Source is set to General source: semi/smm1

Time Study Result: Success Simulation completed.

Cases 2 - Two way couplying semi<--->ht

In Semiconductor Material Model, Model Input Temperature is set to - Temperature (ht) Heat Source is set to General source: semi/smm1 (No other change is introduced to the model)

Time Study Result: ERROR - Division by zero. failed to evaluate variable comp1. semi.iig1.alpha p

When I try to debugg this error, I trace the cause to Impact Ionization Generation model, where iig1.Epn and iig1.Epp parallel electric field variables are calculated in the following manner: Epn = 0 if Epn<1E5, max(Epn)=1E6 ; Epp = 0 if Epp<1E5, max(Epp)=1E6

So my question would be, why am I seeing this error when I try a two-way coupling? Is there a initial condition that should be altered in Case 2 ? Alternatively, is there a way to alter the value of iig1.Epn and iig1.Epp parallel electric field variables? Any wisdom on this would be much appreaciated. Thanks!