Alternative scheme for vibrations, better handling of scipy optimizations