for maximum performance, you may want to try out the
latest alpha test version of VMD, since the gofr code
is now multi-threaded. for instructions to get those
builds see: http://www.ks.uiuc.edu/Research/vmd/alpha/
You can do this, but you can't specify r > pair_cutoff. LAMMPS
doesn't know anything about neighbors further away. You could
reset the pair cutoff to be long right before using these commands
and then use a longer r.