The interaction between the particles would be simple LJ potential with cut off of 2.5. In that case I think I should only modify the “pair_lj_cut.cpp”. As you see in the above expression there is an extra coefficient which I would like to introduce. I am looking for “pair_coeff.cpp/h” file in src directory to do this but I can’t find it.
SD