  We establish the well-posedness of the transient van Roosbroeck system in
three space dimensions under realistic assumptions on the data: non-smooth
domains, discontinuous coefficient functions and mixed boundary conditions.
Moreover, within this analysis, recombination terms may be concentrated on
surfaces and interfaces and may not only depend on charge-carrier densities,
but also on the electric field and currents. In particular, this includes
Avalanche recombination. The proofs are based on recent abstract results on
maximal parabolic and optimal elliptic regularity of divergence-form operators.
