Dear Mark and Chris,
Thanks a lot for your suggestions which gave me idea to solve the problem.
I have run the minimization again, without the FLEXIBLE water molecules.
After the minimization, structure looks perfectly ok to me. please have a
look at the new snapshoot.
http://i269.photobucket.com/albums/jj58/gromacs/after_em.jpg
It is just my speculation that because I am using the TIP4P water model,
which has dummy atom (-1.04 charge), it might be possible that due to this
water hydrogen's (+0.52 charge) are attracted to the dummy atom of nearby
water molecules. That's why if I closely observed the structures (after
minimization using flexible water model) then, some water molecules
aggregates, and leave some vacuum at other places, that's why structure was
looking weird. It is just my speculation please correct me if I am missing
something.
I have also checked the bond and angle energy values and when I was using
the flexible water model it deviate a lot. For bond energy it increases
8.88556e+04 to 4.10925e+07 and for angle energy 3.60214e+04 to 1.09253e+06
these energy compensated by short range Coulomb which changes -1.65898e+07
to -1.71941e+08 .
Without flexible water molecules there is no such huge deviation from the
original value.Also during the equilibration run I have not observed any
LINCS warnings. I am hoping my problem is solved now :-)
Thanks again to your valuable suggestions.
Best regards,
Alok Jain
----- Original Message -----
From: "Mark Abraham" <[EMAIL PROTECTED]>
To: "Discussion list for GROMACS users" <gmx-users@gromacs.org>
Sent: Friday, February 01, 2008 5:04 AM
Subject: Re: [gmx-users] Weird structureafter minimization (membrane
proteinsimulation
[EMAIL PROTECTED] wrote:
Dear Mark, & Chris,
Thanks for your useful suggestions.
Yes Mark define = -DFLEXIBLE is the only difference in both the cases.
Actually in my previous post I have pasted only small section of
tip4p.itp
file,(sorry if it created unnecessary confusions/doubt). This time I am
pasting full tip4p.itp file. As you expected, it has more #ifdef
statements.
OK
I have not modified this file. this is the original file and comes with
the
gromacs-3.2.1 distribution. I have cross checked it with the clean
distribution.
Is this distribution the only one you have installed? If not, you want to
take care that GMXLIB is pointing to the right spot so that you're
actually including the .itp files you think you are.
I realise that make_hole requires some early version of GROMACS, but once
it has been run, perhaps you should take its output structure and use a
3.3.x version from there.
I have also compared the order of TIP4P atoms in my structure file
(gro/pdb) and the same is pasted below as TIP4P.itp file. Order is same
in
both the
files.
OK, well there goes that idea.
Mark
_______________________________________________
gmx-users mailing list gmx-users@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 [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
_______________________________________________
gmx-users mailing list gmx-users@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 [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
