[gmx-users] g_rdf

[Dallas Warren]

An RDF is normalised to the density for the entire box, so you should simply be using that.

Catch ya,

Dr. Dallas Warren
Medicinal Chemistry and Drug Action
Monash Institute of Pharmaceutical Sciences, Monash University
381 Royal Parade, Parkville VIC 3010
dallas.warren at monash.edu
+61 3 9903 9304
---------------------------------
When the only tool you own is a hammer, every problem begins to resemble a nail.

From: gmx-users-bounces at gromacs.org [mailto:gmx-users-bounces at gromacs.org] On Behalf Of mohammad agha
Sent: Tuesday, 31 January 2012 7:26 PM
To: gmx-users at gromacs.org
Subject: [gmx-users] g_rdf

Dear Prof.

I am confused about generation a radial density graph (density vs distance from center of mass), I know that I should use g_rdf. Then I count the number of atoms in the shells around the COM of special group by " trjorder -com -nshell -r " , next I use from this formula for compute the density:
density = (dN/4*pi*dr*r^2)*g(r), where dN, dr, r and g(r) are the number of atoms in each shell, wide of shell, radius of COM of special group, and averaged RDF in each shell, respectively.
While the plot of g(r) is consistent with published previous results but density has not consistence with them, I think that I am wrong in compute the density.
Please help me.

Best Regards
Sara
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20120131/bbb76efd/attachment.html>