In solar cell applications I need to model PN junction, In lumerical's FDTD we can create new material in material database, there we can make any material with desired np density (electron hole density), but after that my simulation diverges due to material gain I suppose.