Hi Francesco,

If you use lammps you can calculate g(r) directly.

Also, if you have the trajectories you can have it using OVITO.

http://www.ovito.org/

Open the trajectorie

go to "Add modification" and then to "coordination analysis": there you have it.

best

A

I used GROMACS to generate the trajectories.

I was able to obtain some radial distribution functions with VMD, however these are the so-called h(r) (basically a counting of the contacts), while I would like to have true g(r).

use gromacs tools gmx_rdf:

http://manual.gromacs.org/programs/gmx-rdf.html

In case you are still wondering how to do it with VMD:

http://www.ks.uiuc.edu/Research/vmd/plugins/gofrgui/