Bond-order based empirical force field models (EFF) (e.g, Tersoff-type BOP, Reactive force field model (ReaxFF) ) are yet another class of interatomic potentials that attempt to account for bond directionality in materials via an angular dependence. The existing set of parameters for these pre-defined models perform well for bulk structures. But in recent years with the advancements in synthesis and applications of nano system materials, the need for having a computational route to assist with the design of novel low-dimensional materials has become increasingly crucial. But these popular and widely available force-field models which have been successful for bulk systems fail to capture the properties of these clusters accurately (e.g. see Fig 1 shows comparison of errors in prediction of energies and forces between some best stateof-the-art potential models vs. DFT energies for a wide size range of cluster configurations sampled from near to far-from equilibrium). Most of the elemental systems have mean absolute errors (MAE) around » 0.5 eV/atom illustrating the need for better models that can describe the interactions in lowdimensional systems. (meV/Å) using force field models available in literature such as EAM/fs , EAM/alloy , MEAM , ADP , SC/EAM )