Re: [gmx-users] Fatal error when using Hbond restraints

2007-03-17 Thread David van der Spoel

Robert Johnson wrote:

Hello everyone,
I'm performing a simulation that includes a 2D infinite SiO2
substrate. I'm not interested in the high frequency motion of the
substrate bonds, so I'm applying position restraints to all the SiO2
atoms. I have other molecules in the system such as water and DNA that
require hydrogen bond restraints. However, I'm getting the following
error when trying to run the system across several nodes:
Fatal error:
Shake block crossing node boundaries
constraint between atoms (1680,1716)

The atoms listed are in the substrate. Since I'm already restraining
the substrate atoms, I don't really need to constrain the Hbonds in
the SiO2. Is there any way to apply Hbond restraints only to a group
of atoms?

Alternatively, how does Gromacs decide which bonds contain hydrogen? I
could just alter the topology file for the substrate to trick Gromacs
into thinking that the hydrogens in the SiO2 are really a different
type of atom.
Thanks,
Bob Johnson
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
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

Bob,

try using
constraints = all-bonds,
or make sure that your charge groups overlap with the constraint groups.

--
David.

David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,  75124 Uppsala, Sweden
phone:  46 18 471 4205  fax: 46 18 511 755
[EMAIL PROTECTED]   [EMAIL PROTECTED]   http://folding.bmc.uu.se

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
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


Re: [gmx-users] Fatal error when using Hbond restraints

2007-03-17 Thread David van der Spoel

Robert Johnson wrote:

Hello everyone,
I'm performing a simulation that includes a 2D infinite SiO2
substrate. I'm not interested in the high frequency motion of the
substrate bonds, so I'm applying position restraints to all the SiO2
atoms. I have other molecules in the system such as water and DNA that
require hydrogen bond restraints. However, I'm getting the following
error when trying to run the system across several nodes:
Fatal error:
Shake block crossing node boundaries
constraint between atoms (1680,1716)

The atoms listed are in the substrate. Since I'm already restraining
the substrate atoms, I don't really need to constrain the Hbonds in
the SiO2. Is there any way to apply Hbond restraints only to a group
of atoms?

Alternatively, how does Gromacs decide which bonds contain hydrogen? I
could just alter the topology file for the substrate to trick Gromacs
into thinking that the hydrogens in the SiO2 are really a different
type of atom.
Thanks,
Bob Johnson
I'm at this very moment looking at 
http://bugzilla.gromacs.org/show_bug.cgi?id=110
which reports the same error. You might try to use pbc=full and 
contraint_algorithm = lincs to circumvent the problem.


The problem you have does not arise from the position restraints though, 
and, by the way, you probably don't need those restraints at all. I 
don't know what you mean with hbond restraints though...


--
David.

David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,  75124 Uppsala, Sweden
phone:  46 18 471 4205  fax: 46 18 511 755
[EMAIL PROTECTED]   [EMAIL PROTECTED]   http://folding.bmc.uu.se

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
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] Fatal error when using Hbond restraints

2007-03-17 Thread Robert Johnson

Hello everyone,
I'm performing a simulation that includes a 2D infinite SiO2
substrate. I'm not interested in the high frequency motion of the
substrate bonds, so I'm applying position restraints to all the SiO2
atoms. I have other molecules in the system such as water and DNA that
require hydrogen bond restraints. However, I'm getting the following
error when trying to run the system across several nodes:
Fatal error:
Shake block crossing node boundaries
constraint between atoms (1680,1716)

The atoms listed are in the substrate. Since I'm already restraining
the substrate atoms, I don't really need to constrain the Hbonds in
the SiO2. Is there any way to apply Hbond restraints only to a group
of atoms?

Alternatively, how does Gromacs decide which bonds contain hydrogen? I
could just alter the topology file for the substrate to trick Gromacs
into thinking that the hydrogens in the SiO2 are really a different
type of atom.
Thanks,
Bob Johnson
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
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