where I and J run over the atom indices in the system, and R IJ is the distance between pair of atoms. V rep is generally considered to be an empirical function that is determined by fitting data from a higher level electronic structure calculations such as DFT. In this work, the adjustment of this potential was carried out using the TANGO code where the repulsive potential is approximated by the following function: