I grab force and dist using compute pair/local e.g.
compute fd all pair/local force dist
I now want to multiply the force per pair by the corresponding distance per pair.
[I am after the derivative of the LJ potential with respect to sigma to get its corresponding conjugate variable λ).
dU_lj/ds = -(dU_lj/dr_ij) (rij/s)
I can dump them and do it outside lammps, but it would be nice to do it within as the data set is huge for large simulations and I need the average over a run.
Thanks,
Rene