For stage 1, we used 8 radial symmetry functions G 1 each with different η values. This decision to use only the radial symmetry functions at this stage was taken to leverage the radially symmetric nature of the Sutton Chen EAM forcefields. For stage 2 along with the above radial symmetry functions 22 angular symmetry functions were also used with a distinct set of parameters. This allowed us to map the increased complexities of the potential energy surface described by DFT for each of the elements in this study.