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