“…The design of numerical schemes for drift-diffusion models is also an mature but still very active field of research (see for instance [9,11,25,28,29,34,35,36,40]). In order to ensure the quality of the numerical simulation and the stability of the numerical method, efforts have been made towards the design of structure preserving schemes [6,12,24,31,34]. Their aim is to preserve physical features of the original model such as the decay of free-energy or non-negativity of solutions.…”