Dear All

 I would like to compute the radial density function rho(r) (in g/cm3) of
different parts of micelle formed with DPC molecules (such as headgroup,
alkyl tail and water) relative to the center of mass of the aggregate. I
have constructed an index files containing all the atoms of the surfactant
headgroup and alkyl chain and used g_rdf (4.5.3):


 /work/cont003/abel01/gmx4.5.3/bin/g_rdf_mpi -f *_Center_All.xtc -s
run_1.tpr -o DPC-Self-AMB_100-155ns_HEADGROUP_RDF.xvg -b 100000 -e 155000
-norm -com -n index_Radial_Profile_DPC.ndx



 For water, as expected, I obtained a RDF profile where rho (r) values tend
to the water density (~1g/cm3) near the simulation  edge  box.  For the
headgroup and the alkyl tail of the DPC, it is not the case since I obtained
two peaks with rho (r) values higher than the density a phosphocholine group
and a dodecyl alkyl chain (in my case 26.4 and 16.4 ).


 What is wrong with command used above. Did i forget something ?


 Thank you in advance for your help.

SA
-- 
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to