[gmx-users] About the example case A toy peptide in water and other protein case (CG martini force field )

2008-10-24 Thread wang
Dear Sir , I used your CG parameters and gromacs V3.3 to simulate protein aggregation. And my command is : grompp_mpi_d -np 2 -f aggre.mdp -c aggre.gro -p aggre.top -o aggre.trp The error message is : double-checking input for internal consistency... Velocities were taken from a Maxwell

[gmx-users] About the example case A toy peptide in water and other protein case (CG martini force field )

2008-10-24 Thread wang
Dear Sir , I used your CG parameters and gromacs V3.3 to simulate protein aggregation. And my command is : grompp_mpi_d -np 2 -f aggre.mdp -c aggre.gro -p aggre.top -o aggre.trp The error message is : double-checking input for internal consistency... Velocities were taken from a Maxwell

[gmx-users] About the example case A toy peptide in water and other protein case (CG martini force field )

2008-10-24 Thread wang
Dear Sir , I used your CG parameters and gromacs V3.3 to simulate protein aggregation. And my command is : grompp_mpi_d -np 2 -f aggre.mdp -c aggre.gro -p aggre.top -o aggre.trp The error message is : double-checking input for internal consistency... Velocities were taken from a Maxwell

[gmx-users] cannot generate .top file

2008-10-24 Thread Kwee Hong
Hi. I seem having some problem in running ./pdb2gmx. Here is the output file. Can you help me in giving me some idea to deal with this? bcsb09s-imac52:bin bcsb09$ pdb2gmx -f 1aml.pdb -p 1aml.top -o 1aml.gro :-) G R O M A C S (-: Getting the Right Output

RE: [gmx-users] stochastic dynamics , langevin

2008-10-24 Thread Berk Hess
From: [EMAIL PROTECTED] To: gmx-users@gromacs.org Subject: RE: [gmx-users] stochastic dynamics , langevin Date: Thu, 23 Oct 2008 15:32:02 -0700 Along these same lines... (a)is there an thermostat setting in GROMACS that comes close to the DPD thermostat as described in (T. Soddemann,

[gmx-users] table_bonded.xvg

2008-10-24 Thread A.Rzepiela
Dear Users In gromacs 4.0, in file force.c in function make_bonded_tables line (1007) sprintf(tabfn + strlen(basefn) - strlen(ftp2ext(efXVG)) - 1,_%s%d.%s, tabext, i, ftp2ext(efXVG)); produces table_bonded_bonded.xvg instead of table_bonded.xvg Greetings Andrzej Rzepiela

Re: [gmx-users] cannot generate .top file

2008-10-24 Thread Mark Abraham
Kwee Hong wrote: Hi. I seem having some problem in running ./pdb2gmx. Here is the output file. Can you help me in giving me some idea to deal with this? bcsb09s-imac52:bin bcsb09$ pdb2gmx -f 1aml.pdb -p 1aml.top -o 1aml.gro snip ---

[gmx-users] Warning with L-BFGS minimization

2008-10-24 Thread Jochen Hub
Hi, with l-bfgs minimzation, grompp always reports WARNING 1 [file em.mdp, line unknown]: For efficient BFGS minimization, use switch/shift/pme instead of cut-off. even if the mdp gives ; Method for doing electrostatics coulombtype = PME rcoulomb_switch = 0 rcoulomb

Re: [gmx-users] Warning with L-BFGS minimization

2008-10-24 Thread Mark Abraham
Jochen Hub wrote: Hi, with l-bfgs minimzation, grompp always reports WARNING 1 [file em.mdp, line unknown]: For efficient BFGS minimization, use switch/shift/pme instead of cut-off. even if the mdp gives ; Method for doing electrostatics coulombtype = PME rcoulomb_switch

[gmx-users] Gromacs 4 SRC RPM?

2008-10-24 Thread Mike Hanby
Is there any eta on when the gromacs-4.0-1.src.rpm is going to show up on: ftp://ftp.gromacs.org/pub/gromacs/rpm Thanks, Mike ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search

Re: [gmx-users] Warning with L-BFGS minimization

2008-10-24 Thread Jochen Hub
Mark Abraham wrote: Jochen Hub wrote: Hi, with l-bfgs minimzation, grompp always reports WARNING 1 [file em.mdp, line unknown]: For efficient BFGS minimization, use switch/shift/pme instead of cut-off. even if the mdp gives ; Method for doing electrostatics coulombtype =

Re: [gmx-users] Warning with L-BFGS minimization

2008-10-24 Thread David van der Spoel
Jochen Hub wrote: Mark Abraham wrote: Jochen Hub wrote: Hi, with l-bfgs minimzation, grompp always reports WARNING 1 [file em.mdp, line unknown]: For efficient BFGS minimization, use switch/shift/pme instead of cut-off. even if the mdp gives ; Method for doing electrostatics coulombtype

RE: [gmx-users] Warning with L-BFGS minimization

2008-10-24 Thread Berk Hess
Hi, The problem is you vdwtype setting, not coulombtype. Berk Date: Fri, 24 Oct 2008 16:08:14 +0200 From: [EMAIL PROTECTED] To: gmx-users@gromacs.org Subject: Re: [gmx-users] Warning with L-BFGS minimization Mark Abraham wrote: Jochen Hub wrote: Hi, with l-bfgs minimzation,

Re: [gmx-users] Gromacs 4 SRC RPM?

2008-10-24 Thread Jussi Lehtola
On Fri, 2008-10-24 at 09:06 -0500, Mike Hanby wrote: Is there any eta on when the gromacs-4.0-1.src.rpm is going to show up on: ftp://ftp.gromacs.org/pub/gromacs/rpm You may be able to use the one in Fedora, e.g.

RE: [gmx-users] Gromacs 4 SRC RPM?

2008-10-24 Thread Mike Hanby
Thanks for the link, I'll take a look at that one. I have to build my own RPM to compile it with Intel compilers and OpenMPI. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jussi Lehtola Sent: Friday, October 24, 2008 9:44 AM To: Discussion list for

Re: [gmx-users] Warning with L-BFGS minimization

2008-10-24 Thread Tsjerk Wassenaar
Hi Jochen, It's probably the rvdw cutoff (from readir.c): if (ir-eI == eiLBFGS (ir-coulombtype==eelCUT || ir-vdwtype==evdwCUT) ir-rvdw != 0) { warning(For efficient BFGS minimization, use switch/shift/pme instead of cut-off.); } Bummer! Tsjerk On 10/24/08, Jochen Hub [EMAIL

RE: [gmx-users] table_bonded.xvg

2008-10-24 Thread Berk Hess
Hi, According to me it does not. I never change the default names. But with the default names the table file name option is table.xvg and for a tabulated bond table number 0, it correctly opens table_b0.xvg Berk From: [EMAIL PROTECTED] To: gmx-users@gromacs.org Date: Fri, 24 Oct 2008

[gmx-users] atoms not part of any of the T-coupling groups

2008-10-24 Thread garry
I am a novice at Gromacs (although I used Gromos 20+ years ago!)I have a dimer with two zincs per monomer, and I get the following complaint when running the grompp with pr.dmp, which I assume is referring to the 4 zincs. There are: 6789 OTHER residues There are: 262

Re: [gmx-users] atoms not part of any of the T-coupling groups

2008-10-24 Thread Justin A. Lemkul
garry wrote: I am a novice at Gromacs (although I used Gromos 20+ years ago!)I have a dimer with two zincs per monomer, and I get the following complaint when running the grompp with pr.dmp, which I assume is referring to the 4 zincs. There are: 6789 OTHER residues There are:

[gmx-users] How to assign/make ionion bonds ?

2008-10-24 Thread Chih-Ying Lin
Hi How to assign / make ionion bonds? Thank you Lin ___ 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

[gmx-users] Does gromacs simulate rotation correctly?

2008-10-24 Thread Arthur Roberts
Hi, all, Just curious. Does gromacs simulate rotation correctly? From my observation, it would appear that every object in the simulation both small and large has a similar rotational correlation time, which is not reality. I would expect water for example to have an extremely short

Re: [gmx-users] How to assign/make ionion bonds ?

2008-10-24 Thread Justin A. Lemkul
Chih-Ying Lin wrote: Hi How to assign / make ionion bonds? Assign them yourself in the [ bonds ] section of the topology. As for what parameters to assign to such a species, that's up to you; they probably don't exist in the standard force fields within Gromacs... -Justin Thank you

Re: [gmx-users] Does gromacs simulate rotation correctly?

2008-10-24 Thread Xavier Periole
On Fri, 24 Oct 2008 10:42:26 -0700 Arthur Roberts [EMAIL PROTECTED] wrote: Hi, all, Just curious. Does gromacs simulate rotation correctly? From my observation, it would appear that every object in the simulation both small and large has a similar rotational correlation time, which is

[gmx-users] The structure property

2008-10-24 Thread Chih-Ying Lin
Hi From the following tutorial, the protein pdb file is downloaded from the Protein Data Bank. Before running the MD simulation, we have to make sure the structure property. I have some questions about this tutorials. http://nmr.chem.uu.nl/~tsjerk/course/md-tutorial/ 1. Make sure that there are

Re: [gmx-users] The structure property

2008-10-24 Thread Justin A. Lemkul
Chih-Ying Lin wrote: Hi From the following tutorial, the protein pdb file is downloaded from the Protein Data Bank. Before running the MD simulation, we have to make sure the structure property. I have some questions about this tutorials. http://nmr.chem.uu.nl/~tsjerk/course/md-tutorial/

Re: [gmx-users] atoms not part of any of the T-coupling groups

2008-10-24 Thread Justin A. Lemkul
Wouldn't it be better to couple the Protein and Zn2+ together (i.e., index group Protein_Zn2+)? Coupling ions separately often leads to weird behavior/explosions when they're free in the solvent. I've never dealt with proteins that have ions bound, so I'm curious to see what everyone else

Re: [gmx-users] The structure property

2008-10-24 Thread Tsjerk Wassenaar
Hi Lin, Justin is right on all points. The first is indeed a check you should always perform, but in this case (1LW9), there are no missing residues/atoms. The point is that one should be aware that there could. The same thing goes with the hint for the histidines. I could have not written you

RE: [gmx-users] f(x) g(x) h(x)in the user defined potential functions

2008-10-24 Thread He, Yang
Hi Mark, Thank you for your reply. From your point of view, I wonder whether you mean I need to write a table based on different interaction types and then change the source code for recognition when running my case in the gromacs .I just have no idea about that .Can you give me further

[gmx-users] gromacs 3.3.1/4.0 and fftw 3.1.2

2008-10-24 Thread chris . neale
Hello, can anybody confirm if fftw-3.1.2 is compatible with gromacs 3.3.1 and gromacs 4.0? A colleague has recently indicated that I must use fftw 2.1.5 for gromacs mpi. However, I have been using fftw-3.1.2 for quite a while now. It compiles fine and my searching today turned up some

[gmx-users] Is there a way to measure the shortest distance between two groups of atoms?

2008-10-24 Thread Arthur Roberts
Hi, all, I was wondering if there is a way to measure the shortest distance between two groups of atoms? g_dist seems to only measure the average distance between the two groups. Your input would be greatly appreciated. Sincerely, Art Roberts

[gmx-users] Is there a way to measure the shortest distance between two groups of atoms?

2008-10-24 Thread chris . neale
g_mindist (with or without -pi depending on what you want) Hi, all, I was wondering if there is a way to measure the shortest distance between two groups of atoms? g_dist seems to only measure the average distance between the two groups. Your input would be greatly appreciated.

Re: [gmx-users] gromacs 3.3.1/4.0 and fftw 3.1.2

2008-10-24 Thread Nicolas
[EMAIL PROTECTED] a écrit : Hello, can anybody confirm if fftw-3.1.2 is compatible with gromacs 3.3.1 and gromacs 4.0? A colleague has recently indicated that I must use fftw 2.1.5 for gromacs mpi. However, I have been using fftw-3.1.2 for quite a while now. It compiles fine and my

[gmx-users] The code : from .pdb to .gro

2008-10-24 Thread Chih-Ying Lin
Hi Could anyone please to direct me the code / the command (pdb2gmx) , which can transfer .pdb file to .gro file? Thank you Lin ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the

Re: [gmx-users] gromacs 3.3.1/4.0 and fftw 3.1.2

2008-10-24 Thread Mark Abraham
[EMAIL PROTECTED] wrote: Hello, can anybody confirm if fftw-3.1.2 is compatible with gromacs 3.3.1 and gromacs 4.0? A colleague has recently indicated that I must use fftw 2.1.5 for gromacs mpi. However, I have been using fftw-3.1.2 for quite a while now. It compiles fine and my searching

Re: [gmx-users] f(x) g(x) h(x)in the user defined potential functions

2008-10-24 Thread Mark Abraham
He, Yang wrote: Hi Mark, Thank you for your reply. From your point of view, I wonder whether you mean I need to write a table based on different interaction types and then change the source code for recognition when running my case in the gromacs. Yes that's what I said. I just have no

Re: [gmx-users] The code : from .pdb to .gro

2008-10-24 Thread Mark Abraham
Chih-Ying Lin wrote: Hi Could anyone please to direct me the code / the command (pdb2gmx) , which can transfer .pdb file to .gro file? It's unlikely you actually need a .gro file for your purpose. pdb2gmx is a tool for making topology files. editconf is the easiest tool for making simple

[gmx-users] gromacs 3.3.1/4.0 and fftw 3.1.2

2008-10-24 Thread chris . neale
Hi Mark, did you post something on this topic? There is a link indicating that you did but it is broken for some reason. The links above and below work fine. If you wouldn't mind, could you please repost? Thanks, Chris. ___ gmx-users mailing