optimized hermiteintegrals, overloaded for bomd
energy calc runs with optimized R, gradient with old R