Move runInterpolation from GaussianTS into parent QMReaction class

Most of this should be general to Gaussian and MOPAC
Make tsScript try to run Nudged Elastic Band
fixup! Move runInterpolation from GaussianTS into parent QMReaction class
Making NEB more general for Mopac and Gaussian.

Reading the start and end geometries directly from the MOL files,
instead of from (optimized) log files.
Running calculations with EMT
Running calculations with MOPAC
Try using PM7
Try to use FIRE optimizer - wasn't as good.
Turn on log file
Try turning on climbing after initial NEB optimization.

Doesn't work very well :-/
Try FIRE with Climbing after initial NEB worked.

THis still takes for ever
Give up on climbing. It's not working (or takes a long time to converge)
Find the highest image after the NEB and save that geometry.

This should be our guess of the TS geometry.