[gmx-users] density profile of bilayers

[Zoltan Varga]

Hello,

priyanka srivastava <priyankaps4 at yahoo.com> wrote: Dear all,

I wish to calculate the density profile of a lipid bilayer i.e the electron density and this is how I am using the g_density command:

g_density -f final.trr -s final.tpr -ei electrons.dat -o density.xvg -ed -symm

where electrons.dat contains # of electrons  for each and every atomtype defined in the system e.g.
C2A = 6
NTM = 7
OA   = 8 etc....

However it reports electron density for 50bins (setup by -sl) which ranges from 0 to 50. However I want it to vary from, -30 to +30. 

Maybe I'm wrong, but it reports the electron density along the z axis of the box with a grid what contains 50 bins. So to get an electron density which ranges from -30 to +30 (in Angstroms and with the bilayer center at 0, I guess) you should only shift the z values with half of your box size. 

Another point is, I want the density profile separately for P, N, acyl, water etc. Should I define a separate electrons.dat for this purpose? 
No, you should make an index file for that which contains P, N, etc. groups and run g_density -n new_index.ndx ...

Kindly suggest me as to how should I proceed,
eagerly waiting for your reply,
Priyanka S...
        

---------------------------------
Looking for a deal? Find great prices on flights and hotels with Yahoo! FareChase._______________________________________________
gmx-users mailing list    gmx-users at gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-request at gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Hope that helped,

Zoltan Varga

       
---------------------------------
Need a vacation? Get great deals to amazing places on Yahoo! Travel. 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20071017/d72a724f/attachment.html>