Re: [gmx-users] anisotropic pressure coupling for temperature annealing

2007-10-10 Thread Alan Dodd
Have you simulated the system without the annealing, to see if it does the same 
thing?  Bilayers are very sensitive to minor changes in setup, using someone 
else's equilibrated bilayer coordinates doesn't mean it's equilibrated for your 
system.


- Original Message 
From: Q733 [EMAIL PROTECTED]
To: gmx-users@gromacs.org
Sent: Wednesday, October 10, 2007 6:40:06 AM
Subject: [gmx-users] anisotropic pressure coupling for temperature annealing


Dear gmx-users, I used anisotropic p-coupling for lipid bilayers and annealed 
the system from 330K-300K at the ratio of 2 degree/ns.And the x dimension and y 
dimension increased with time in the whole procedure, but actually the 
area/lipid should decrease when tempreture decreases. I think there may due to 
the P couple parameters I used. Now I write them here:
 
 tau_p   =  2.0  2.0 2.0 0   0   0
compressibility =  4.5e-6   4.5e-6  4.5e-6  0   0   0
ref_p   =  11   1   1   1   1
 
Is there anything wrong with my parameter? or should I use
 
 
ref_p   =  11   1   0   0   0   
 
I saw someone use this ref_p  before.
 
Any suggestion will be appreciated, thanks in advance.


  

Catch up on fall's hot new shows on Yahoo! TV. Watch previews, get listings, 
and more!
http://tv.yahoo.com/collections/3658 ___
gmx-users mailing listgmx-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

Re: [gmx-users] How to change the box size during simuation

2007-10-10 Thread Alan Dodd
You can apply external forces to subsets of atoms, in such a way as to cause 
your box to shrink by itself.


- Original Message 
From: WU Yanbin [EMAIL PROTECTED]
To: gmx-users@gromacs.org
Sent: Tuesday, October 9, 2007 11:40:21 PM
Subject: [gmx-users] How to change the box size during simuation

Hey,
  Now I'm simulating infinite molecules. I want to change the box size a lit 
bit every time step so as to induce some stress in this system.
  Is there any way in gromacs to do this, apart from modifying the code? Thanks 
in advance. 
   Yours Sincerely,
WU Yanbin


  

Fussy? Opinionated? Impossible to please? Perfect.  Join Yahoo!'s user panel 
and lay it on us. http://surveylink.yahoo.com/gmrs/yahoo_panel_invite.asp?a=7 
___
gmx-users mailing listgmx-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

RE: [gmx-users] How to change the box size during simuation

2007-10-10 Thread Berk Hess





From: WU Yanbin [EMAIL PROTECTED]
Reply-To: Discussion list for GROMACS users gmx-users@gromacs.org
To: gmx-users@gromacs.org
Subject: [gmx-users] How to change the box size during simuation
Date: Tue, 9 Oct 2007 17:40:21 -0500

Hey,
  Now I'm simulating infinite molecules. I want to change the box size a 
lit

bit every time step so as to induce some stress in this system.
  Is there any way in gromacs to do this, apart from modifying the code?


You can use pressure coupling to apply a stress and induce a strain
or you can use the deform option to apply a strain and induce stress.

Berk.

_
Talk with your online friends with Messenger 
http://www.join.msn.com/messenger/overview


___
gmx-users mailing listgmx-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] What is this warning on running grompp ? WARNING 1 [file md-new.mdp, line unknown]

2007-10-10 Thread maria goranovic
Hi,
When running grompp I get the following error often

###
checking input for internal consistency...
WARNING 1 [file md-new.mdp, line unknown]:
###

It does not look like a critical error, but I want to make sure. Did not
find any specific discussions on the mailing list immediately.

My md-new.mdp file is attached.

Thank you for the help.

Also, how can one send a reply to a specific message in the list ?

Thanks very much !

-Maria
-- 
Maria G.
Technical University of Denmark
Copenhagen


md-new.mdp
Description: application/mdp
___
gmx-users mailing listgmx-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] xtc error

2007-10-10 Thread luciano.pgomes
Hi All,

I had this problem running MD on a single Debian machine:

Program mdrun_d, VERSION 3.3.1
Source code file: stat.c, line: 257

Fatal error:
XTC error


There is no chance of being size exceeded of file or quota, or disk error.
I compiled gromacs with gcc-3.4...

Does anyone know what this error means? Gcc version may be the cause?

Thanks,

Luciano P Gomes


___
gmx-users mailing listgmx-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


Re: [gmx-users] xtc error

2007-10-10 Thread David van der Spoel

luciano.pgomes wrote:

Hi All,

I had this problem running MD on a single Debian machine:

Program mdrun_d, VERSION 3.3.1
Source code file: stat.c, line: 257

Fatal error:
XTC error


it can also mean that your system is exploding.



There is no chance of being size exceeded of file or quota, or disk error.
I compiled gromacs with gcc-3.4...

Does anyone know what this error means? Gcc version may be the cause?

Thanks,

Luciano P Gomes


___
gmx-users mailing listgmx-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



--
David van der Spoel, Ph.D.
Molec. Biophys. group, Dept. of Cell  Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
[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 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


Re: [gmx-users] Criteria and Energy Conservation

2007-10-10 Thread David van der Spoel

[EMAIL PROTECTED] wrote:

 Dear GROMACS Users,

 I'm trying to use GROMACS to get some MD snapshots for single-point quantum
mechanical calculations, trying to observe HOMO-LUMO gap of a protein immersed
in water with counterions, at a given salt concentration (4600 atoms in total).
I'd like to know if they are reliable enough to send to QM calculations. I've
got two questions: Which are the criteria I should be looking at? And one of
them I've analyzed was energy conservation (I'm doing 300K NVT MD, with 10A
cutoff for vdw and coulomb, with PME electrostatics, and double precision), and
the manual says a Ekin RMSD/Tot Energy RMSD ratio below 0.05 is fine, but I've
got 0.3-0.5 in all simulations I've done. What should I do?

This is because of T coupling. If you turn off T coupling it will be 
better (but you will not have constant T, and the energy will dirft in 
the long run). Your starting structures are fine for QM, but QM 
minimization will take a LONG time.




 Thanks in advance,

 Gustavo Troiano Feliciano



This message was sent using IFUSP Webmail - USP/Sao Paulo/Brazil.
___
gmx-users mailing listgmx-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



--
David van der Spoel, Ph.D.
Molec. Biophys. group, Dept. of Cell  Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
[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 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] simulating hydrocarbons in water

2007-10-10 Thread Adam Fraser
Hello,

I'm very new to Gromacs, and I am interested in simulating the interactions
of 2 hexadecane (C16H34) molecules in water (SPC/E specificall).  I've spent
some time experimenting with tutorials in Gromacs, but I've found little
help in non-protein simulations like this one.

I've built topologies for hexadecane compatible with NAMD, so I figure I
should be able to do the same in Gromacs but I'm fuzzy on how to go about
doing so (what files to edit).  I was also hoping to get some input on which
forcefield would be best for this sort of system.

I'd greatly appreciate it if someone could give me some pointers on how to
get started here.

Thanks very much,
Adam
___
gmx-users mailing listgmx-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

Re: [gmx-users] simulating hydrocarbons in water

2007-10-10 Thread David van der Spoel

Adam Fraser wrote:

Hello,

I'm very new to Gromacs, and I am interested in simulating the 
interactions of 2 hexadecane (C16H34) molecules in water (SPC/E 
specificall).  I've spent some time experimenting with tutorials in 
Gromacs, but I've found little help in non-protein simulations like this 
one.


I've built topologies for hexadecane compatible with NAMD, so I figure I 
should be able to do the same in Gromacs but I'm fuzzy on how to go 
about doing so (what files to edit).  I was also hoping to get some 
input on which forcefield would be best for this sort of system.


I'd greatly appreciate it if someone could give me some pointers on how 
to get started here.


Thanks very much,
Adam


if you have a pdb file you can run x2top (3.3.2 only).





___
gmx-users mailing listgmx-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



--
David van der Spoel, Ph.D.
Molec. Biophys. group, Dept. of Cell  Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
[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 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


Re: [gmx-users] simulating hydrocarbons in water

2007-10-10 Thread Adam Fraser
Yes I tried this and it doesn't work for me.  I get:

Fatal error:
Library file ffG53a6.n2t not found in current dir nor in default
directories.
(You can set the directories to search with the GMXLIB path variable)

Also, doesn't x2top take in a gro file?  (not a pdb)  I wrote a script that
converts pdb to gro, and I still get the above error.

-Adam

On 10/10/07, David van der Spoel [EMAIL PROTECTED] wrote:

 Adam Fraser wrote:
  Hello,
 
  I'm very new to Gromacs, and I am interested in simulating the
  interactions of 2 hexadecane (C16H34) molecules in water (SPC/E
  specificall).  I've spent some time experimenting with tutorials in
  Gromacs, but I've found little help in non-protein simulations like this
  one.
 
  I've built topologies for hexadecane compatible with NAMD, so I figure I
  should be able to do the same in Gromacs but I'm fuzzy on how to go
  about doing so (what files to edit).  I was also hoping to get some
  input on which forcefield would be best for this sort of system.
 
  I'd greatly appreciate it if someone could give me some pointers on how
  to get started here.
 
  Thanks very much,
  Adam
 
 if you have a pdb file you can run x2top (3.3.2 only).

 
  
 
  ___
  gmx-users mailing listgmx-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


 --
 David van der Spoel, Ph.D.
 Molec. Biophys. group, Dept. of Cell  Molec. Biol., Uppsala University.
 Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
 [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 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 listgmx-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

Re: [gmx-users] simulating hydrocarbons in water

2007-10-10 Thread TJ Piggot

Just as a note editconf will covert a pdb to a gro file

Tom

--On Wednesday, October 10, 2007 09:47:03 -0400 Adam Fraser 
[EMAIL PROTECTED] wrote:



Yes I tried this and it doesn't work for me.  I get:


Fatal error:
Library file ffG53a6.n2t not found in current dir nor in default
directories.
(You can set the directories to search with the GMXLIB path variable)

Also, doesn't x2top take in a gro file?  (not a pdb)  I wrote a script
that converts pdb to gro, and I still get the above error.

-Adam


On 10/10/07, David van der Spoel [EMAIL PROTECTED] wrote:

Adam Fraser wrote:

Hello,

I'm very new to Gromacs, and I am interested in simulating the
interactions of 2 hexadecane (C16H34) molecules in water (SPC/E
specificall).  I've spent some time experimenting with tutorials in
Gromacs, but I've found little help in non-protein simulations like this
one.

I've built topologies for hexadecane compatible with NAMD, so I figure I
should be able to do the same in Gromacs but I'm fuzzy on how to go
about doing so (what files to edit).  I was also hoping to get some
input on which forcefield would be best for this sort of system.

I'd greatly appreciate it if someone could give me some pointers on how
to get started here.

Thanks very much,
Adam


if you have a pdb file you can run x2top (3.3.2 only).





___
gmx-users mailing listgmx-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



--
David van der Spoel, Ph.D.
Molec. Biophys. group, Dept. of Cell  Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
[EMAIL PROTECTED][EMAIL PROTECTED]   http://folding.bmc.uu.se
___
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






--
TJ Piggot
[EMAIL PROTECTED]
University of Bristol, UK.

___
gmx-users mailing listgmx-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


Re: [gmx-users] Criteria and Energy Conservation

2007-10-10 Thread gtroiano

 Thank you very much again, David! Then it's fine for me, since what I want to
do is a single point QM in termalized snapshots of the dynamics and not QM
minimization.

 Gustavo

 [EMAIL PROTECTED] wrote:
   Dear GROMACS Users,
 
   I'm trying to use GROMACS to get some MD snapshots for single-point
 quantum
  mechanical calculations, trying to observe HOMO-LUMO gap of a protein
 immersed
  in water with counterions, at a given salt concentration (4600 atoms in
 total).
  I'd like to know if they are reliable enough to send to QM calculations.
 I've
  got two questions: Which are the criteria I should be looking at? And one
 of
  them I've analyzed was energy conservation (I'm doing 300K NVT MD, with 10A
  cutoff for vdw and coulomb, with PME electrostatics, and double precision),
 and
  the manual says a Ekin RMSD/Tot Energy RMSD ratio below 0.05 is fine, but
 I've
  got 0.3-0.5 in all simulations I've done. What should I do?
 
 This is because of T coupling. If you turn off T coupling it will be
 better (but you will not have constant T, and the energy will dirft in
 the long run). Your starting structures are fine for QM, but QM
 minimization will take a LONG time.


   Thanks in advance,
 
   Gustavo Troiano Feliciano
 
 
  
  This message was sent using IFUSP Webmail - USP/Sao Paulo/Brazil.
  ___
  gmx-users mailing listgmx-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


 --
 David van der Spoel, Ph.D.
 Molec. Biophys. group, Dept. of Cell  Molec. Biol., Uppsala University.
 Box 596, 75124 Uppsala, Sweden. Phone:+46184714205. Fax: +4618511755.
 [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 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






This message was sent using IFUSP Webmail - USP/Sao Paulo/Brazil.
___
gmx-users mailing listgmx-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


Re: [gmx-users] simulating hydrocarbons in water

2007-10-10 Thread David van der Spoel

Adam Fraser wrote:

Yes I tried this and it doesn't work for me.  I get:

Fatal error:
Library file ffG53a6.n2t not found in current dir nor in default 
directories.

(You can set the directories to search with the GMXLIB path variable)



did you use 3.3.2 ? The default and only supported FF there is OPLS. For 
GROMOS variants you can try the prodrg webserver.


Also, doesn't x2top take in a gro file?  (not a pdb)  I wrote a script 
that converts pdb to gro, and I still get the above error.


-Adam

On 10/10/07, * David van der Spoel* [EMAIL PROTECTED] 
mailto:[EMAIL PROTECTED] wrote:


Adam Fraser wrote:
  Hello,
 
  I'm very new to Gromacs, and I am interested in simulating the
  interactions of 2 hexadecane (C16H34) molecules in water (SPC/E
  specificall).  I've spent some time experimenting with tutorials in
  Gromacs, but I've found little help in non-protein simulations
like this
  one.
 
  I've built topologies for hexadecane compatible with NAMD, so I
figure I
  should be able to do the same in Gromacs but I'm fuzzy on how to go
  about doing so (what files to edit).  I was also hoping to get some
  input on which forcefield would be best for this sort of system.
 
  I'd greatly appreciate it if someone could give me some pointers
on how
  to get started here.
 
  Thanks very much,
  Adam
 
if you have a pdb file you can run x2top (3.3.2 only).

 
 


 
  ___
  gmx-users mailing listgmx-users@gromacs.org
mailto: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]
mailto:[EMAIL PROTECTED].
  Can't post? Read http://www.gromacs.org/mailing_lists/users.php
http://www.gromacs.org/mailing_lists/users.php


--
David van der Spoel, Ph.D.
Molec. Biophys. group, Dept. of Cell  Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
[EMAIL PROTECTED]
mailto:[EMAIL PROTECTED][EMAIL PROTECTED]
mailto:[EMAIL PROTECTED]   http://folding.bmc.uu.se
___
gmx-users mailing list gmx-users@gromacs.org
mailto: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]
mailto:[EMAIL PROTECTED].
Can't post? Read http://www.gromacs.org/mailing_lists/users.php





___
gmx-users mailing listgmx-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



--
David van der Spoel, Ph.D.
Molec. Biophys. group, Dept. of Cell  Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
[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 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] xtc error

2007-10-10 Thread luciano.pgomes
Hi,

 Hi All,

 I had this problem running MD on a single Debian machine:

 Program mdrun_d, VERSION 3.3.1
 Source code file: stat.c, line: 257

 Fatal error:
 XTC error


 There is no chance of being size exceeded of file or quota, or disk error.
 I compiled gromacs with gcc-3.4...

 Does anyone know what this error means? Gcc version may be the cause?

 Thanks,

 Luciano P Gomes



 it can also mean that your system is exploding.

No.
The xtc file doesn't opened. It crashed VMD.

I forgot saying that I was running simulated annealing protocol with a DNA 
molecule and crystallographic water molecules (HOH) freezed, only moving ions 
and water molecules placed by genbox (tip3p). The crash occured on 2 ps of 
simulation time in a temperature of 50K. I had excluded non-bonded 
interactions: DNA DNA HOH DNA HOH HOH.

I am running now a simple MD at 300K and this error haven't occured yet.
I am trying to understand...


Thanks,
Luciano

___
gmx-users mailing listgmx-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


Re: [gmx-users] Criteria and Energy Conservation

2007-10-10 Thread David van der Spoel

[EMAIL PROTECTED] wrote:

 Thank you very much again, David! Then it's fine for me, since what I want to
do is a single point QM in termalized snapshots of the dynamics and not QM
minimization.
But then you will get systematically incorrect energies due to weird 
(force field dependent) bond lengths etc.


 Gustavo


[EMAIL PROTECTED] wrote:

 Dear GROMACS Users,

 I'm trying to use GROMACS to get some MD snapshots for single-point

quantum

mechanical calculations, trying to observe HOMO-LUMO gap of a protein

immersed

in water with counterions, at a given salt concentration (4600 atoms in

total).

I'd like to know if they are reliable enough to send to QM calculations.

I've

got two questions: Which are the criteria I should be looking at? And one

of

them I've analyzed was energy conservation (I'm doing 300K NVT MD, with 10A
cutoff for vdw and coulomb, with PME electrostatics, and double precision),

and

the manual says a Ekin RMSD/Tot Energy RMSD ratio below 0.05 is fine, but

I've

got 0.3-0.5 in all simulations I've done. What should I do?


This is because of T coupling. If you turn off T coupling it will be
better (but you will not have constant T, and the energy will dirft in
the long run). Your starting structures are fine for QM, but QM
minimization will take a LONG time.



 Thanks in advance,

 Gustavo Troiano Feliciano



This message was sent using IFUSP Webmail - USP/Sao Paulo/Brazil.
___
gmx-users mailing listgmx-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


--
David van der Spoel, Ph.D.
Molec. Biophys. group, Dept. of Cell  Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
[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 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







This message was sent using IFUSP Webmail - USP/Sao Paulo/Brazil.
___
gmx-users mailing listgmx-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



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


Re: [gmx-users] xtc error

2007-10-10 Thread David van der Spoel

luciano.pgomes wrote:

Hi,


Hi All,

I had this problem running MD on a single Debian machine:

Program mdrun_d, VERSION 3.3.1
Source code file: stat.c, line: 257

Fatal error:
XTC error


There is no chance of being size exceeded of file or quota, or disk error.
I compiled gromacs with gcc-3.4...

Does anyone know what this error means? Gcc version may be the cause?

Thanks,

Luciano P Gomes





it can also mean that your system is exploding.


No.
The xtc file doesn't opened. It crashed VMD.

I forgot saying that I was running simulated annealing protocol with a DNA 
molecule and crystallographic water molecules (HOH) freezed, only moving ions 
and water molecules placed by genbox (tip3p). The crash occured on 2 ps of 
simulation time in a temperature of 50K. I had excluded non-bonded 
interactions: DNA DNA HOH DNA HOH HOH.

I am running now a simple MD at 300K and this error haven't occured yet.
I am trying to understand...


did you use pressure coupling? you can not have Pcoupl with freezing and 
constraints.



Thanks,
Luciano

___
gmx-users mailing listgmx-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



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


Re: [gmx-users] limit on energygrp_table

2007-10-10 Thread Adrien Leygue
Thanks for your suggestions,

I have gone for the solution where I change the type of group numbers
from uchars to unsigned short.

It seems however that this solution is not very elegant:

I fail to understand the logic of limiting the number of groups in the
energygrp_table field to 254. Indeed the energygrp_table field is made
of pairs of groups. So if the group numbers are stored as uchar the
number of possible pairwise interactions between those groups is much
greater than 254 (of the order of 2^16 and not 2^8).

I am actually in that case: I have not too many groups (~20) so I can
store all the numbers in a uchar but when I generate all possible
pairwise interactions I have to input something like 400 group names
under energygrp_table and then grompp complains.

 -Is there a way to increase this limit (by changing a constant and
 then recompile) ? or is it more complicated (from the value 254 one
 could guess that there is some one one-byte optimization)
 

 This is indeed a on-byte optimization.
 The group numbers are stored as uchars (include/types/atoms.h).
 You could change this, then you also have to change 2 or more other files
 for the io and communication of these numbers.


I have modified include/types/atoms.h
Then it seemed logical to also change the following lines in src/kernel/readir.c
#define MAXPTR 65534
#define NOGID  65535
to account for the new limits.

I can compile the modified code and run grompp without a warning on my
system. The only strange thing is that grompp complains if I do not
provide tablep*.xvg files even though I want to use the same tabulated
potentials for my inter and intra (defined through some pairs)
molecular non-bonded potentials.

Nevertheless when I try to run some energy minimization it directly
complains that my first 1-4 interaction is beyond table range of 0
nm,... and will ignore all 1-4 (i assume it means pairs).
This message seems really odd,... especially the 0nm part.

I have probably to propagate the change of type to more files. I
welcome all suggestions.

Adrien.

-- 
Adrien Leygue

Department of Materials Science and Engineering
School of Chemical Engineering
National Technical University of Athens
9 Heroon Polytechniou Street
Zografou Campus, Athens 157 80,
GREECE
Tel/Fax: +30210 772-3112
___
gmx-users mailing listgmx-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] Problems prior to running lipid-protein simulation

2007-10-10 Thread adak Nasiripour
Dear All;
when I try to run pdb2gmx command for making .top file of POPC+Protein, an 
error appears: Residue 'POPC' not found in residue topology database.
I have copied all necessary .itp files in share folder and also in the path 
where it is run.
  what could be the problem and how to solve it? (I have passed these steps 
long time ago but now I am not able to reproduce the procedure!)
thanks;


Adak Nasiripour.

   
-
 Check out  the hottest 2008 models today at Yahoo! Autos.___
gmx-users mailing listgmx-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

Re: [gmx-users] simulating hydrocarbons in water

2007-10-10 Thread Adam Fraser
Much appreciated, works beautifully!

On 10/10/07, TJ Piggot [EMAIL PROTECTED] wrote:

 Just as a note editconf will covert a pdb to a gro file

 Tom

 --On Wednesday, October 10, 2007 09:47:03 -0400 Adam Fraser
 [EMAIL PROTECTED] wrote:

  Yes I tried this and it doesn't work for me.  I get:
 
 
  Fatal error:
  Library file ffG53a6.n2t not found in current dir nor in default
  directories.
  (You can set the directories to search with the GMXLIB path variable)
 
  Also, doesn't x2top take in a gro file?  (not a pdb)  I wrote a script
  that converts pdb to gro, and I still get the above error.
 
  -Adam
 
 
  On 10/10/07, David van der Spoel [EMAIL PROTECTED] wrote:
 
  Adam Fraser wrote:
  Hello,
 
  I'm very new to Gromacs, and I am interested in simulating the
  interactions of 2 hexadecane (C16H34) molecules in water (SPC/E
  specificall).  I've spent some time experimenting with tutorials in
  Gromacs, but I've found little help in non-protein simulations like
 this
  one.
 
  I've built topologies for hexadecane compatible with NAMD, so I figure
 I
  should be able to do the same in Gromacs but I'm fuzzy on how to go
  about doing so (what files to edit).  I was also hoping to get some
  input on which forcefield would be best for this sort of system.
 
  I'd greatly appreciate it if someone could give me some pointers on how
  to get started here.
 
  Thanks very much,
  Adam
 
  if you have a pdb file you can run x2top (3.3.2 only).
 
 
 
 
 
  ___
  gmx-users mailing listgmx-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
 
 
  --
  David van der Spoel, Ph.D.
  Molec. Biophys. group, Dept. of Cell  Molec. Biol., Uppsala University.
  Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
  [EMAIL PROTECTED][EMAIL PROTECTED]   http://folding.bmc.uu.se
  ___
  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
 
 



 --
 TJ Piggot
 [EMAIL PROTECTED]
 University of Bristol, UK.

 ___
 gmx-users mailing listgmx-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 listgmx-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

Re: [gmx-users] simulating hydrocarbons in water

2007-10-10 Thread Adam Fraser
Ah, okay thanks very much!  I thought-sure I had the most recent version,
but apparently not.  Still there's one quirk...

x2top -f HexA.gro -o HexA.top -r HexA.rtp -alldih
...This works to generate a topology file, but still, the program is
seg-faulting before writing out the residue type file (*.rtp).

I'm not yet sure if I can get by without the .rtp since hexadecane (C16)
won't likely be a recognized residue type.

Any ideas why this is seg-faulting or how I can get around the fact that
restype (C16) won't be recognized?

Thanks so much,
-Adam


On 10/10/07, David van der Spoel [EMAIL PROTECTED] wrote:

 Adam Fraser wrote:
  Yes I tried this and it doesn't work for me.  I get:
 
  Fatal error:
  Library file ffG53a6.n2t not found in current dir nor in default
  directories.
  (You can set the directories to search with the GMXLIB path variable)
 

 did you use 3.3.2 ? The default and only supported FF there is OPLS. For
 GROMOS variants you can try the prodrg webserver.

  Also, doesn't x2top take in a gro file?  (not a pdb)  I wrote a script
  that converts pdb to gro, and I still get the above error.
 
  -Adam
 
  On 10/10/07, * David van der Spoel* [EMAIL PROTECTED]
  mailto:[EMAIL PROTECTED] wrote:
 
  Adam Fraser wrote:
Hello,
   
I'm very new to Gromacs, and I am interested in simulating the
interactions of 2 hexadecane (C16H34) molecules in water (SPC/E
specificall).  I've spent some time experimenting with tutorials
 in
Gromacs, but I've found little help in non-protein simulations
  like this
one.
   
I've built topologies for hexadecane compatible with NAMD, so I
  figure I
should be able to do the same in Gromacs but I'm fuzzy on how to
 go
about doing so (what files to edit).  I was also hoping to get
 some
input on which forcefield would be best for this sort of system.
   
I'd greatly appreciate it if someone could give me some pointers
  on how
to get started here.
   
Thanks very much,
Adam
   
  if you have a pdb file you can run x2top (3.3.2 only).
 
   
   
 
 
 
   
___
gmx-users mailing listgmx-users@gromacs.org
  mailto: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]
  mailto:[EMAIL PROTECTED].
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
  http://www.gromacs.org/mailing_lists/users.php
 
 
  --
  David van der Spoel, Ph.D.
  Molec. Biophys. group, Dept. of Cell  Molec. Biol., Uppsala
 University.
  Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax:
 +4618511755.
  [EMAIL PROTECTED]
  mailto:[EMAIL PROTECTED][EMAIL PROTECTED]
  mailto:[EMAIL PROTECTED]   http://folding.bmc.uu.se
  ___
  gmx-users mailing list gmx-users@gromacs.org
  mailto: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]
  mailto:[EMAIL PROTECTED].
  Can't post? Read http://www.gromacs.org/mailing_lists/users.php
 
 
 
  
 
  ___
  gmx-users mailing listgmx-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


 --
 David van der Spoel, Ph.D.
 Molec. Biophys. group, Dept. of Cell  Molec. Biol., Uppsala University.
 Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
 [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 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 listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 

[gmx-users] Re:gmx-users Digest, Vol 42, Issue 36

2007-10-10 Thread luciano.pgomes
HI David,

I used -DPOSRE, no constraints but temperature coupling (berendsen), no 
pressure coupling.

thanks,
Luciano

 luciano.pgomes wrote:
  Hi,
 
  Hi All,
 
  I had this problem running MD on a single Debian machine:
 
  Program mdrun_d, VERSION 3.3.1
  Source code file: stat.c, line: 257
 
  Fatal error:
  XTC error
 
 
  There is no chance of being size exceeded of file or quota, or disk error.
  I compiled gromacs with gcc-3.4...
 
  Does anyone know what this error means? Gcc version may be the cause?
 
  Thanks,
 
  Luciano P Gomes
 
 
 
  it can also mean that your system is exploding.
 
  No.
  The xtc file doesn't opened. It crashed VMD.
 
  I forgot saying that I was running simulated annealing protocol with a DNA 
  molecule and crystallographic water molecules (HOH) freezed, only moving 
  ions and water molecules placed by genbox (tip3p). The crash occured on 2 
  ps of simulation time in a temperature of 50K. I had excluded non-bonded 
  interactions: DNA DNA HOH DNA HOH HOH.
 
  I am running now a simple MD at 300K and this error haven't occured yet.
  I am trying to understand...

 did you use pressure coupling? you can not have Pcoupl with freezing and
 constraints.
 
 
  Thanks,
  Luciano
 
  ___
  gmx-users mailing listgmx-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


 --
 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 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] Re: Problems prior to running lipid-protein simulation

2007-10-10 Thread Bruce Milne

Dear Adak,

[EMAIL PROTECTED] wrote:

Dear All;
when I try to run pdb2gmx command for making .top file of POPC+Protein, an 
error appears: Residue 'POPC' not found in residue topology database.
I have copied all necessary .itp files in share folder and also in the path 
where it is run.
  what could be the problem and how to solve it? (I have passed these steps 
long time ago but now I am not able to reproduce the procedure!)
thanks;


Adak Nasiripour.
  
I think the process you might be wanting to perform is to run pdb2gmx 
for your protein only (with suitably adjusted coordinates so that it 
matches your bilayer coordinates) and then add #include lines to the 
resulting topology to add in the POPC part *after* the force field 
#include statement:


#include ffgmx.itp
#include lipids.itp
#include popc.itp

- assuming that these are your .itp file names. Then adjust the [ 
molecules ] section to list the number and type of your components and 
you should have a working .top file.


You don't have to run pdb2gmx for the whole system at any point (it's 
not going to work anyway unless you add the POPC parameters to all the 
force field files and then add the POPC residue to the  .rtp file so 
that  pdb2gmx can recognize it).


Cheers,

Bruce

--
Hamewith - the road that's never dreary

Dr. Bruce F. Milne
CEQOFFUP
Faculdade de Farmácia
Universidade do Porto
Rua Aníbal Cunha - 164
4050-047
Porto
Portugal

email: [EMAIL PROTECTED]

___
gmx-users mailing listgmx-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


Re: [gmx-users] Criteria and Energy Conservation

2007-10-10 Thread gtroiano

 Well, yeah. But I want to look at the HOMO-LUMO gap of the system in the living
environment. Since the system is too big to do MD with the QM code, I tried to
use the best force field description I could, and the smallest box as well, to
get structures that could represent the system at 300K. And the forces from QM
are not so bad so far, within the order I'd expect from MD with the QM code.
Anyway, thanks very much for the hints so far, David, and if you have some
suggestion or anything at all to add, feel free to do so.

 [EMAIL PROTECTED] wrote:
   Thank you very much again, David! Then it's fine for me, since what I want
 to
  do is a single point QM in termalized snapshots of the dynamics and not
 QM
  minimization.
 But then you will get systematically incorrect energies due to weird
 (force field dependent) bond lengths etc.
 
   Gustavo
 
  [EMAIL PROTECTED] wrote:
   Dear GROMACS Users,
 
   I'm trying to use GROMACS to get some MD snapshots for single-point
  quantum
  mechanical calculations, trying to observe HOMO-LUMO gap of a protein
  immersed
  in water with counterions, at a given salt concentration (4600 atoms in
  total).
  I'd like to know if they are reliable enough to send to QM calculations.
  I've
  got two questions: Which are the criteria I should be looking at? And one
  of
  them I've analyzed was energy conservation (I'm doing 300K NVT MD, with
 10A
  cutoff for vdw and coulomb, with PME electrostatics, and double
 precision),
  and
  the manual says a Ekin RMSD/Tot Energy RMSD ratio below 0.05 is fine, but
  I've
  got 0.3-0.5 in all simulations I've done. What should I do?
 
  This is because of T coupling. If you turn off T coupling it will be
  better (but you will not have constant T, and the energy will dirft in
  the long run). Your starting structures are fine for QM, but QM
  minimization will take a LONG time.
 
 
   Thanks in advance,
 
   Gustavo Troiano Feliciano
 
 
  
  This message was sent using IFUSP Webmail - USP/Sao Paulo/Brazil.
  ___
  gmx-users mailing listgmx-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
 
  --
  David van der Spoel, Ph.D.
  Molec. Biophys. group, Dept. of Cell  Molec. Biol., Uppsala University.
  Box 596, 75124 Uppsala, Sweden. Phone: +46184714205. Fax: +4618511755.
  [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 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
 
 
 
 
 
  
  This message was sent using IFUSP Webmail - USP/Sao Paulo/Brazil.
  ___
  gmx-users mailing listgmx-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


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






This message was sent using IFUSP Webmail - USP/Sao Paulo/Brazil.
___
gmx-users mailing listgmx-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 

RE: [gmx-users] Problems prior to running lipid-protein simulation

2007-10-10 Thread Dallas B. Warren
http://wiki.gromacs.org/index.php/Errors#Residue_.27XXX.27_not_found_in_
residue_topology_database

You only run pdb2gmx on molecules that you want to generate a topology
file for.  For the lipids, you will already have the .itp files for it,
or should have.  So you don't run pdb2gmx on it.

Catch ya,

Dr. Dallas Warren
Lecturer
Department of Pharmaceutical Biology and Pharmacology
Victorian College of Pharmacy, Monash University
381 Royal Parade, Parkville VIC 3010
[EMAIL PROTECTED]
+61 3 9903 9524
-
When the only tool you own is a hammer, every problem begins to resemble
a nail. 
___
gmx-users mailing listgmx-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


Re: [gmx-users] Dummies usage doubts

2007-10-10 Thread David van der Spoel

Jones de Andrade wrote:

Hi all.

Well, I'm having a bit of trouble here because the work has decided to 
go in the direction of something I've never used before with gromacs: 
dummy atoms (or now virtual sites).


I assume this is not for a protein, in which case pdb2gmx does it for you.


I've got at least three questions in order to use them:

1 - does the virtual sites need to be included in the .gro files in 
order to make a simulation run?



Yes.


2 - do I need to include them in the exclusion list?

It depends, but usually yes.



3 - do I need to include them in the atoms section?

Yes.




Thanks a lot in advance for any help anyone can provide in these 
subjects. ;)

Check the mannual too. Chapter 5. An dtopologies for TIP4P, TIP5P etc.




Sincerally yours,

Jones




___
gmx-users mailing listgmx-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



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