Berk Hess wrote:
Hi,
Are you 100% sure you have no other charge groups in your system
that consist of only virtual sites, but at least with 2 virtual sites?
Ah, you're right. One of the CG which caused the crash is the following,
with 2 H's in one CG:
1 amber99_39 1 NPRO N 1 -0.202 16.026 ;
qtot -0.202
2 amber99_17 1 NPRO H1 2 (!) 0.312
0 ; qtot 0.11
3 amber99_17 1 NPRO H2 2 (!) 0.312 0
; qtot 0.422
4 amber99_11 1 NPRO CD 3 -0.012 14.026 ;
qtot 0.41
Note however, that the CG which also caused crashes was different, with
2 vsites and 1 heavy atom in the same group, plus a separate CG with 3
hydrogens. According to the lincs warning and 1-4 distance error, the
hydrogens seemed to be involved. Typically the warning/error mentioned
the hydrogen to have moved/rotated too much:
1 MNH3 1 NASN MN1 1 0 8.517 ;
qtot 0
2 MNH3 1 NASN MN2 1 0 8.517 ;
qtot 0
3 amber99_39 1 NASN N 1 0.1801 0 ;
qtot 0.1801
4 amber99_17 1 NASN H1 2 0.1921 0 ;
qtot 0.3722
5 amber99_17 1 NASN H2 2 0.1921 0 ;
qtot 0.5643
6 amber99_17 1 NASN H3 2 0.1921 0 ;
qtot 0.7564
7 amber99_11 1 NASN CA 3 0.0368 13.018 ;
qtot 0.7932
I just fixed a bug with this situation for 4.0.6 and 4.1.
You can try the git master or 4.0 release branch.
I will get the lates gmx with your bugfix and write back if the crashes
still appear.
I don't understand why only your terminal NH3 becomes a single charge
group,
unless the hydrogens were generated with the .hdb file,
while they were present for all other NH3 groups.
I don't understand that either. The result is the same no matter if I
use -ignh, or if the hydrogens are present in the input gro file or not.
For examle, using
pdb2gmx -f asn[-h].gro -ff amber99sb [-ignh]
with asn-h.gro or asn.gro are just the ASN at the N-terminus (with or
without hydrogens), e.g. asn-h.gro:
Gallium Rubidium Oxygen Manganese Argon Carbon Silicon
8
97NASN N 1423 0.062 -1.317 1.235
97NASN CA 1427 0.109 -1.209 1.154
97NASN CB 1429 0.233 -1.249 1.073
97NASN CG 1432 0.213 -1.368 0.981
97NASN OD1 1433 0.263 -1.478 1.005
97NASN ND2 1434 0.143 -1.345 0.867
97NASN C 1437 0.003 -1.148 1.068
97NASN O 1438 0.033 -1.090 0.964
0.00000 0.00000 0.00000
The written topology is:
[ atoms ]
; nr type resnr residue atom cgnr charge mass
typeB chargeB massB
1 amber99_39 1 NASN N 1 0.1801 14.01 ;
qtot 0.1801
2 amber99_17 1 NASN H1 2 0.1921 1.008 ;
qtot 0.3722
3 amber99_17 1 NASN H2 2 0.1921 1.008 ;
qtot 0.5643
4 amber99_17 1 NASN H3 2 0.1921 1.008 ;
qtot 0.7564
5 amber99_11 1 NASN CA 3 0.0368 12.01 ;
qtot 0.7932
6 amber99_28 1 NASN HA 4 0.1231 1.008 ;
qtot 0.9163
...
Weird. The same thing happens by the way with a proline N-terminus, no
matter if hydrogens are present before pdb2gmx or not.
So thanks a lot for the help,
Jochen
Berk
> Date: Wed, 30 Sep 2009 16:45:39 +0200
> From: [email protected]
> To: [email protected]
> Subject: Re: [gmx-users] vsites, ffamber, charge groups, N-terminus,
and crashes
>
> Hey Berk,
>
> I just found that the simulation with a separate CG on each atom also
> crashed, but much later (after 11ns), that why I just noticed that. I
> have repeated the simulations now to make really sure that the crashes
> are reproducible.
>
> Jochen
>
>
> Jochen Hub wrote:
> >
> > I have only tested on 8 nodes, and 4 x 2 x 1 DD. Should I try
someting
> > else?
> >
> > Jochen
> >
> >
> > Berk Hess wrote:
> >> Hi,
> >>
> >> One question:
> >> Are these crashes single processor, with domain decomposition, or
both?
> >>
> >> Berk
> >>
> >> > Date: Wed, 30 Sep 2009 15:35:42 +0200
> >> > From: [email protected]
> >> > To: [email protected]
> >> > Subject: [gmx-users] vsites, ffamber, charge groups, N-terminus,
> >> and crashes
> >> >
> >> > Hi,
> >> >
> >> > I have been trying to use vsites with AMBER99SB, but the simulation
> >> > frequently cashes after a few hundred ps with the suspects (lincs
> >> > warnings, 1-4 distance error). I could pinpoint the problem to
> >> erroneous
> >> > charge group assignments generated by pdb2gmx at the NH3 group
at the
> >> > N-terminus. The following charge group assignments caused the
error:
> >> >
> >> > ; nr type resnr residue atom cgnr charge mass
> >> > 1 MNH3 1 NASN MN1 1 0 8.517 ;
> >> > qtot 0
> >> > 2 MNH3 1 NASN MN2 1 0 8.517 ;
> >> > qtot 0
> >> > 3 amber99_39 1 NASN N 1 0.1801 0 ;
> >> > qtot 0.1801
> >> > 4 amber99_17 1 NASN H1 2 0.1921 0 ;
> >> > qtot 0.3722
> >> > 5 amber99_17 1 NASN H2 2 0.1921 0 ;
> >> > qtot 0.5643
> >> > 6 amber99_17 1 NASN H3 2 0.1921 0 ;
> >> > qtot 0.7564
> >> > 7 amber99_11 1 NASN CA 3 0.0368 13.018 ;
> >> > qtot 0.7932
> >> > ...
> >> >
> >> > The NH3 groups in the lysines were fine and did not cause any
> >> error. The
> >> > only difference compared to the N-terminus is that each of the
three H
> >> > in the NH3 has its own charge group, but the vsites and the N are
> >> still
> >> > in the same CG:
> >> >
> >> > 30 amber99_11 2 LYP CE 29 -0.0143 14.026 ;
> >> > qtot 0.9937
> >> > 31 amber99_28 2 LYP HE1 30 0.1135 0 ;
> >> > qtot 1.107
> >> > 32 amber99_28 2 LYP HE2 31 0.1135 0 ;
> >> > qtot 1.221
> >> > 33 MNH3 2 LYP MNZ1 32 0 8.517 ;
> >> > qtot 1.221
> >> > 34 MNH3 2 LYP MNZ2 32 0 8.517 ;
> >> > qtot 1.221
> >> > 35 amber99_39 2 LYP NZ 32 -0.3854 0 ;
> >> > qtot 0.8353
> >> > 36 amber99_17 2 LYP HZ1 33 0.34 0 ;
> >> > qtot 1.175
> >> > 37 amber99_17 2 LYP HZ2 34 0.34 0 ;
> >> > qtot 1.515
> >> > 38 amber99_17 2 LYP HZ3 35 0.34 0 ;
> >> > qtot 1.855
> >> >
> >> > Alternatively, I have also tried to give every single atom
(including
> >> > the vsites MN??) a separate charge group, which also worked fine.
> >> >
> >> > 1)
> >> > So, is there any reason why the dummies should be in the same
CG as
> >> the
> >> > nitrogen? Or should it be fine just to assign a separate CG to
> >> every atom?
> >> >
> >> > 2)
> >> > Is a bug report appreciated on that issue?
> >> >
> >> > Best,
> >> >
> >> > Jochen
> >> >
> >> >
> >> > --
> >> > ---------------------------------------------------
> >> > Dr. Jochen Hub
> >> > Molecular Biophysics group
> >> > Dept. of Cell & Molecular Biology
> >> > Uppsala University. Box 596, 75124 Uppsala, Sweden.
> >> > Phone: +46-18-4714451 Fax: +46-18-511755
> >> > ---------------------------------------------------
> >> >
> >> > _______________________________________________
> >> > gmx-users mailing list [email protected]
> >> > http://lists.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
> >>
> >>
------------------------------------------------------------------------
> >> See all the ways you can stay connected to friends and family
> >> <http://www.microsoft.com/windows/windowslive/default.aspx>
> >>
------------------------------------------------------------------------
> >>
> >> _______________________________________________
> >> gmx-users mailing list [email protected]
> >> http://lists.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
> >
> >
>
>
> --
> ---------------------------------------------------
> Dr. Jochen Hub
> Molecular Biophysics group
> Dept. of Cell & Molecular Biology
> Uppsala University. Box 596, 75124 Uppsala, Sweden.
> Phone: +46-18-4714451 Fax: +46-18-511755
> ---------------------------------------------------
>
> _______________________________________________
> gmx-users mailing list [email protected]
> http://lists.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
------------------------------------------------------------------------
Express yourself instantly with MSN Messenger! MSN Messenger
<http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/>
------------------------------------------------------------------------
_______________________________________________
gmx-users mailing list [email protected]
http://lists.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
--
---------------------------------------------------
Dr. Jochen Hub
Molecular Biophysics group
Dept. of Cell & Molecular Biology
Uppsala University. Box 596, 75124 Uppsala, Sweden.
Phone: +46-18-4714451 Fax: +46-18-511755
---------------------------------------------------
_______________________________________________
gmx-users mailing list [email protected]
http://lists.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
