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