Both shape and ordering optimization tasks use surrogate models to evaluate the NP energies at each optimization step. These models rely on a graph partitioning scheme and a topological energy approximation to compute (local) atomic contributions to the energy of single-metal and bimetallic NPs. In addition to providing sufficiently accurate NP total energies, the separation into local contributions allows us to implement a optimalexchange schemes to accelerate established global optimization methods.