output dynamic Fermi level in bulk or QD conduction band
skeleton for general relaxation function
added general relaxation RHS function that can pick between RTA and kinetic models; still need to implement RELAX_RTA function
fixed bug in kinetic model where gamma2 was in int instead of double