where V ≡ y 1 +y 2 2 here, and the parameter V 1 2 = 1 control the steepness of the weight function in the short-range repulsive region. The next step would search the local maximums of the uncertainty function, where the initial points are also generated by the Sobol quasirandom method. These points would then be selected as new points via two criteria. The first criterion is related to the similarity of points in configuration space, which removes the point if there is another point close to it (d(x i , x j ) < d cut ) and has a larger uncertainty.