changed md-procedure to multiple p and VASP reference energy