Nilesh Dhumal wrote:
Hello,

I have a system with a glucose molecule.

I want to calculate the dipole moment of a particular OH in glucose
molecule.I made an index file which have a group containing atoms.

[ O8 ]
  8   18

10 is oxygen no. and 20 is hydroen no.



I'm assuming you mean 8 is O and 18 is H?


But, if I try to use g_dipoles to get the dipole moment contribution from using following command:

g_dipoles -f 6.trr -s 6.tpr  -n index.ndx -corr mol -nonormalize  -c

I get the following error.

Fatal error:
index[1]=8 does not correspond to the first atom of a molecule

I read on gmx-users mailing list.

http://www.mail-archive.com/gmx-users@gromacs.org/msg40091.html

Calculating a dipole moment of a possibly-charged species requires that
there be a reference point, which is conventionally the center of mass of
the molecule. This means all the atoms of the molecule have to be known,
and g_dipoles assumes the index group consists of whole molecules.

So to do the partition you want, you will need to provide the same
molecules, but with zero charges. That will mean making a copy of your
.top and hacking those charges to zero to generate a new such .tpr.

Actually, tpbconv has the ability to set the charges of a group in a .tpr
to zero.

So I tried

tpbconv -f 6.trr -s 6.tpr -n glu-emi-cl-128-no.ndx -e 6.edr  -o zero.tpr


This command doesn't accomplish anything. If you're trying to zero out the charges, you need to use the -zeroq option. But then, if you have two atoms with zero charge, then the dipole will of course be zero!

I got the error

segmentataion fault

How can I solve this problem to calcualte the dipole moment of a bond?


I doubt you can do it with g_dipoles. You've been told already that the group analyzed has to have a net charge of zero. An OH group will not satisfy that need, nor does assigning each of the atoms zero charge. You're better off using g_dist (or maybe g_bond) to find the distance between the O and H atoms over time and plug that into an equation using the partial charges on the atoms to calculate the dipole yourself. The charges are fixed, so the only thing that changes (in the absence of constraints) is the distance between them.

-Justin

I am using gromacs version 4.0.7.

Thanks

Nilesh







--
========================================

Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================
--
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