We carried out AIMD simulations within the Born-Oppenheimer approximation using the FHI-vibes package interfaced to the FHI-aims code for the computation of energies and forces. FHI-aims is an all-electron code based on numerical atom-centered orbitals (NAOs). The simulations were performed at 300 K with a time step of 0.5 fs, using the NVT ensemble with a Langevin thermostat. The electronic structure calculations were performed with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional in combination with the Tkatchenko-Scheffler (TS) correction to account for long-range dispersion interactions. NAOs of tier1 quality in combination with "light" numerical settings were used. The AIMD simulations were equilibrated for approximately 5.0 ps followed by a production run of at least 25.0 ps.