Alternatively, one thing that has come to my mind as a solution inspired by a post by Shern (Force between two groups as a function of coordinates - #2 by srtee) and by the idea you proposed yourself would be is to run a simulation saving a trajectory file. Then you could do a rerun (see rerun command — LAMMPS documentation) using your potential but setting the cutoff to be whatever custom distance you are interested on. And then you could use the pair/local command to output the information you want, since you said it would do the job if you could set a custom cutoff. The pairwise energies would be messed up though since you are using an unrightful cutoff for the force-fields you have, but if you indeed only want to output distances, I think this should work also.
1 Like