from ase.vibrations import Vibrations
vib = Vibrations(atoms)
vib.run()