Note that MgCl2 is considered as one single element in this case since it cannot be separated in this model. For higher likelihood of convergence and finding of the global minimum, the basin-hopping algorithm scipy.optimize.basinhopping 37 is applied with "temperature" parameter of -1000 and the minimizer keyword arguments maxiter of 100 and ftol of 10 -6 .