Recent application of SR for physical models can be found in civil engineering and material science . Although successfully proven, the proposed SR was based on a heuristic search that could terminate the optimization in local minima solutions, producing potentially less suitable models than possible. Additionally, as the structure of a model reflects the actual mechanistic interactions within the system studied, these approximate solutions cannot be used reliably to infer any mechanistic information about the system, i.e. to use it to identify chemical reactions mechanisms with certainty. Acknowledging this disadvantage, SR is formulated as a mixed-integer nonlinear programming (MINLP) in Ref. and solved to global optimum.