New tools for high-dimensional, expensive, black-box global optimization functions applied to ReaxFF parameterizations