In the present study, a variational path integral molecular dynamics method developed by the Miura [Chem. Phys. Lett. 482, (2009) p. 165] is applied to molecular vibrations on adiabatic potential energy surfaces. This method is a molecular dynamics algorithm for a variational path integral method, which can be used to generate the exact ground state numerically. Calculations on a Morse oscillator and a water molecule have been performed. The ground-state energies and the groundstate distributions were accurately calculated.