<html><body><div style="color:#000; background-color:#fff; font-family:bookman old style, new york, times, serif;font-size:14pt"><div id="yiv1015436642"><div><div style="color:#000;background-color:#fff;font-family:bookman old style, new york, times, serif;font-size:14pt;"><div id="yiv1015436642yui_3_2_0_18_132807401390048"><span id="yiv1015436642yui_3_2_0_18_1328074013900157">Dear Prof.</span><span></span></div><div id="yiv1015436642yui_3_2_0_18_1328074013900164"><br><span id="yiv1015436642yui_3_2_0_18_1328074013900157"></span></div><div id="yiv1015436642yui_3_2_0_18_1328074013900165"><span id="yiv1015436642yui_3_2_0_18_1328074013900157">Thank you very much from your reply.</span></div><div id="yiv1015436642yui_3_2_0_18_1328074013900168"><span id="yiv1015436642yui_3_2_0_18_1328074013900157">Is your mean that "density=(N/V)*g(r)", where N is the number of total atoms in box not in the shell and V is the volume of box not in the shell?<br><br>Best
Regards<br>Sara<br></span></div><div id="yiv1015436642yui_3_2_0_18_1328074013900179"><span id="yiv1015436642yui_3_2_0_18_1328074013900157"><br></span></div><div><br id="yiv1015436642yui_3_2_0_18_132807401390051"></div> <div class="yiv1015436642yui_3_2_0_18_132807401390054" id="yiv1015436642yui_3_2_0_18_132807401390056" style="font-family:bookman old style, new york, times, serif;font-size:14pt;"> <div class="yiv1015436642yui_3_2_0_18_132807401390097" style="font-family:times new roman, new york, times, serif;font-size:12pt;"> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1"> <b><span style="font-weight:bold;">From:</span></b> Dallas Warren <Dallas.Warren@monash.edu><br> <b><span style="font-weight:bold;">To:</span></b> Discussion list for GROMACS users <gmx-users@gromacs.org> <br> <b><span style="font-weight:bold;">Sent:</span></b> Wednesday, February 1, 2012 2:12 AM<br> <b><span style="font-weight:bold;">Subject:</span></b> RE:
[gmx-users] g_rdf<br> </font> </div> <br>
<div id="yiv1015436642">
<style><!--
#yiv1015436642
filtered {font-family:Calibri;panose-1:2 15 5 2 2 2 4 3 2 4;}
#yiv1015436642 filtered {font-family:Tahoma;panose-1:2 11 6 4 3 5 4 4 2 4;}
#yiv1015436642 filtered {font-family:"Bookman Old Style";panose-1:2 5 6 4 5 5 5 2 2 4;}
#yiv1015436642
p.yiv1015436642MsoNormal, #yiv1015436642 li.yiv1015436642MsoNormal, #yiv1015436642 div.yiv1015436642MsoNormal
{margin:0cm;margin-bottom:.0001pt;font-size:12.0pt;font-family:"serif";}
#yiv1015436642 a:link, #yiv1015436642 span.yiv1015436642MsoHyperlink
{color:blue;text-decoration:underline;}
#yiv1015436642 a:visited, #yiv1015436642 span.yiv1015436642MsoHyperlinkFollowed
{color:purple;text-decoration:underline;}
#yiv1015436642 span.yiv1015436642EmailStyle17
{font-family:"sans-serif";color:#1F497D;}
#yiv1015436642 .yiv1015436642MsoChpDefault
{font-size:10.0pt;}
#yiv1015436642 filtered {margin:72.0pt 72.0pt 72.0pt 72.0pt;}
#yiv1015436642 div.yiv1015436642WordSection1
{}
--></style>
<div>
<div class="yiv1015436642WordSection1">
<div class="yiv1015436642MsoNormal"><span style="font-size:11.0pt;color:#1F497D;">An RDF is normalised to the density for the entire box, so you should simply be using that.</span></div>
<div class="yiv1015436642MsoNormal"><span style="font-size:11.0pt;color:#1F497D;"> </span></div>
<div>
<div class="yiv1015436642MsoNormal"><span style="font-size:10.0pt;color:#1F497D;">Catch ya,<br>
<br>
Dr. Dallas Warren</span></div>
<div class="yiv1015436642MsoNormal"><span style="font-size:10.0pt;color:#1F497D;">Medicinal Chemistry and Drug Action</span></div>
<div class="yiv1015436642MsoNormal"><span style="font-size:10.0pt;color:#1F497D;">Monash Institute of Pharmaceutical Sciences, Monash University<br>
381 Royal Parade, Parkville VIC 3010<br>
dallas.warren@monash.edu</span></div>
<div class="yiv1015436642MsoNormal"><span style="font-size:10.0pt;color:#1F497D;">+61 3 9903 9304<br>
---------------------------------<br>
When the only tool you own is a hammer, every problem begins to resemble a nail.</span><span style="font-size:11.0pt;color:#1F497D;">
</span></div>
</div>
<div class="yiv1015436642MsoNormal"><span style="font-size:11.0pt;color:#1F497D;"> </span></div>
<div style="border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt;">
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm;">
<div class="yiv1015436642MsoNormal"><b><span style="font-size:10.0pt;">From:</span></b><span style="font-size:10.0pt;"> gmx-users-bounces@gromacs.org [mailto:gmx-users-bounces@gromacs.org]
<b>On Behalf Of </b>mohammad agha<br>
<b>Sent:</b> Tuesday, 31 January 2012 7:26 PM<br>
<b>To:</b> gmx-users@gromacs.org<br>
<b>Subject:</b> [gmx-users] g_rdf</span></div>
</div>
</div>
<div class="yiv1015436642MsoNormal"> </div>
<div>
<div>
<div class="yiv1015436642MsoNormal" style="background:white;"><span style="font-size:14.0pt;color:black;">Dear Prof.</span></div>
</div>
<div>
<div class="yiv1015436642MsoNormal" style="background:white;"><span style="font-size:14.0pt;color:black;"> </span></div>
</div>
<div>
<div class="yiv1015436642MsoNormal" style="background:white;"><span style="font-size:14.0pt;color:black;">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:</span></div>
</div>
<div>
<div class="yiv1015436642MsoNormal" style="background:white;"><span style="font-size:14.0pt;color:black;">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.</span></div>
</div>
<div>
<div class="yiv1015436642MsoNormal" style="background:white;"><span style="font-size:14.0pt;color:black;">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.</span></div>
</div>
<div>
<div class="yiv1015436642MsoNormal" style="background:white;"><span style="font-size:14.0pt;color:black;">Please help me.</span></div>
</div>
<div>
<div class="yiv1015436642MsoNormal" style="background:white;"><span style="font-size:14.0pt;color:black;"> </span></div>
</div>
<div>
<div class="yiv1015436642MsoNormal" style="background:white;"><span style="font-size:14.0pt;color:black;">Best Regards</span></div>
</div>
<div>
<div class="yiv1015436642MsoNormal" style="background:white;"><span style="font-size:14.0pt;color:black;">Sara</span></div>
</div>
</div>
</div>
</div>
</div>
</div><br>-- <br>gmx-users mailing list <a rel="nofollow" ymailto="mailto:gmx-users@gromacs.org" target="_blank" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>http://lists.gromacs.org/mailman/listinfo/gmx-users<br>Please search the archive at <a rel="nofollow" target="_blank" href="http://www.gromacs.org/Support/Mailing_Lists/Search">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>Please don't post (un)subscribe requests to the list. Use the <br>www interface or send it to <a rel="nofollow" ymailto="mailto:gmx-users-request@gromacs.org" target="_blank" href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>Can't post? Read http://www.gromacs.org/Support/Mailing_Lists<br><br> </div> </div> </div></div></div></div></body></html>