Quantum molecular dynamics: propagating wavepackets and density operators using the multiconfiguration time-dependent Hartree method