A semidiscretization in time of a fourth order nonlinear parabolic system in several space dimensions arising in quantum semiconductor modelling is studied. The system is numerically treated by introducing an additional nonlinear potential. Exploiting the stability of the discretization, convergence is shown in the multi-dimensional case. Under some assumptions on the regularity of the solution, the rate of convergence proves to be optimal.