Hi, all
I am now using MEAM in lammps. When calculating the pair potential for element, lammps calculates the total energy of the reference structure as
E=-E0* { 1+a(rij) }* exp{-a(rij)}
I wonder whether I can rewrite this equation with the other universal equation proposed by Rose ( Phys Rev.B (1984) ) with two more adjustable parameters?? How to do this modification?
one more question, how can I just change the value of ‘attrac’ and ‘repuls’ parameters?
Is it the only way that I input them in the ‘meam_setup_global.f’, where the defaults of these two parameters are zero