My MD system consists of several LARGE particles and a number of SMALL particles.
(say, 10 Large + 100000 Small particles)
I use LJ potentials for Large-Large and Large-Small pairs, but I like to impose no interactions between Small-Small pairs.
I could do the simulation by using “pair_style lj/cut” and setting epsilon=0 for Small-Small pairs.
But it is very inefficient because LAMMPS should calculate all Small-Small pairs for nothing.
If you know any way to set no interaction between certain pairs, could you please let me know?
Thank you very much.