Ab Initio Molecular Dynamics: Propagating the Density Matrix with Gaussian Orbitals