Dear lammps user,
I am interested in a simulation which contains 2 carbon nanotubes located one after the other. The ends of the nanotubes (facing each other) are fixed by using set_force = 0 command. I would like to prevent the nanotubes for interacting with each other. i.e. if I apply a 5 N load to the first nanotube, the sum of all the carbon nanotube atoms + the reaction at the fixed boundary should be 5 N for the 1st nanotube and 0 for the second one.
Does anyone know how I can achieve this in LAMMPS?