Molecular rotations and vibrations have been extensively studied by chemists for decades, both experimentally using spectroscopic methods and theoretically with the help of quantum chemistry. However, the theoretical investigation of molecular rotations and vibrations in strong magnetic fields requires computationally more demanding tools. As such, proper calculations of rotational and vibrational spectra were not feasible up until very recently. In this work, we present rotational and vibrational spectra for two small linear molecules, H2 and LiH, in strong magnetic fields. By treating the nuclei as classical particles, trajectories for rotations and vibrations are simulated from ab initio molecular dynamics. Born-Oppenheimer potential energy surfaces are calculated at the Hartree-Fock and MP2 levels of theory using London atomic orbitals to ensure gauge origin invariance. For the calculation of nuclear trajectories, a highly efficient Tajima propagator is introduced, incorporating the Berry curvature tensor accounting for the screening of nuclear charges.