drift diffusion simulations