The sequential one-parameter parabolic interpolation (SOPPI) method has been used as one of the first techniques for the development of ReaxFF force fields . Here, the total error is calculated for three different values of a selected parameter to derive the optimal value using a parabolic fit. However, with an increasing set of parameters, the required number of one-parameter optimisation steps to converge increases dramatically, as only a small portion of the search space is accessible. Also, it requires a reasonable initial guess and a well-considered choice of the order for the parameters to be optimised in each step to obtain an accurate solution .