4. HREMD
5. Fourbody interactions - eliminate pp-contact etc. storage (currently 4body
switched off by default, -DFOURBODY flag switches on).
+6. Fix LBFGS - the linesearch (search) subroutine does not handle difficult
+ cases when changing variables can result in nans.