Molecular dynamics (MD) simulations have been used as a powerful computational tool by researchers in the fields of physics, chemistry, biology, and materials engineering, to model dynamical behavior of several material systems, including gases, liquids, solids, surfaces, and clusters. The accuracy of all these simulations and subsequent calculations is determined by the quality of the potential energy surface (PES) description, in the form of a force field model used in the simulation. The force field model describes the inter particle interactions in the system in question, and it is the most computationally expensive aspect of the MD simulation. High fidelity techniques, for instance, density functional theory (DFT) and ab-initio molecular dynamics (AIMD) are computationally expensive and are limited to much shorter time (∼ ps) and lengthscales (∼ nm). Hence the developement of an accurate classical description of the PES would be beneficial. Empirical potentials such as the Embedded atom model (EAM) which are derived from concepts based on the density functional theory, and its variants like Suttona) Also at Center for Nanoscale Materials, Argonne National Laboratory, Argonne, IL 60439 Chen (SC) and Gupta potentials, , are computationally efficient and have been successful in predicting much of the bulk properties including the elastic coefficients, thermal expansion coefficients , melting points , free energies in both solid and liquid phases , description of defects and many more for metals and metal alloys.