Re: [gmx-users] Genion error

2013-10-24 Thread Justin Lemkul
On 10/24/13 10:21 AM, felipe vasquez wrote: Hi, I chose group 0 (System), but I also tried others like 1 (Protein) or 2 (Protein+H) with the same result. You don't want to embed ions into your protein or haphazardly into the system. You'll start deleting random segments of molecules, or th

Re: [gmx-users] Genion error

2013-10-24 Thread felipe vasquez
Hi, I chose group 0 (System), but I also tried others like 1 (Protein) or 2 (Protein+H) with the same result. Regards, Andrés F. *Andrés Felipe Vásquez J., MSc.* Grupo de Fisiología Molecular Subdirección de Investigación Científica y Tecnológica Dirección de Investigación en Salud Pública In

Re: [gmx-users] Genion error

2013-10-24 Thread Justin Lemkul
On 10/24/13 9:34 AM, felipe vasquez wrote: Hi, I have been trying to add ions to my system so as to make it neutral. However, I always obtain the same result: Solvent Group size (950) is not multiple of 9. What group did you choose at the genion prompt? The command I entered in this step

[gmx-users] Genion error

2013-10-24 Thread felipe vasquez
Hi, I have been trying to add ions to my system so as to make it neutral. However, I always obtain the same result: Solvent Group size (950) is not multiple of 9. The command I entered in this step was: genion -s protein_w.tpr -o protein_solv.gro -conc 0.15 -neutral -pname NA+ -nname CL- How c

Re: [gmx-users] genion doesn't recognize SOL in top file

2013-10-17 Thread Justin Lemkul
On 10/17/13 1:01 PM, sunyeping wrote: Dear Gromacs users, I am trying to add ions to my system using:genion -s ions.tpr -o solv_ions.gro -p topol.top -pname NA -nname CL -np 8 but it returns a error message saying: Fatal error: No line with moleculetype 'SOL' found the [ molecules ] section

[gmx-users] genion doesn't recognize SOL in top file

2013-10-17 Thread sunyeping
Dear Gromacs users, I am trying to add ions to my system using:genion -s ions.tpr -o solv_ions.gro -p topol.top -pname NA -nname CL -np 8 but it returns a error message saying: Fatal error: No line with moleculetype 'SOL' found the [ molecules ] section of file 'topol.top' I checked the topol.

Re: [gmx-users] Genion command not working

2013-08-27 Thread Tsjerk Wassenaar
Hi Deepak, You have to set the minimal distance between ions lower. Check the help of genion. Cheers, Tsjerk On Aug 28, 2013 8:04 AM, "Deepak Ojha" wrote: Dear Gmxers, I want to simulate 5M nabr aqueous solution.To add the ions I use genion command and add 22 na and br ions respectively to t

[gmx-users] Genion command not working

2013-08-27 Thread Deepak Ojha
Dear Gmxers, I want to simulate 5M nabr aqueous solution.To add the ions I use genion command and add 22 na and br ions respectively to the system of pure 108 water molecules.But on running the command genion I get error message Processing topology Back Off! I just backed up topol.top to ./#top

Re: [gmx-users] genion

2013-05-03 Thread Justin Lemkul
On 5/3/13 11:27 AM, Ewaru wrote: Hi, Just wondering. My system shows this: "System has non-zero total charge: 11.98" and thus, I use the command as below to add 12 Cl to the protein: "genion -s ions.tpr -o m4_solv_ions.gro -p topol.top -pname NA -nname CL -nn 12" However, when I checked

[gmx-users] genion

2013-05-03 Thread Ewaru
Hi, Just wondering. My system shows this: "System has non-zero total charge: 11.98" and thus, I use the command as below to add 12 Cl to the protein: "genion -s ions.tpr -o m4_solv_ions.gro -p topol.top -pname NA -nname CL -nn 12" However, when I checked the genion.log file, it shows that th

Re: [gmx-users] Genion cannot read from input

2013-04-27 Thread Justin Lemkul
On 4/27/13 12:28 PM, Justin Lemkul wrote: On 4/27/13 11:21 AM, Nathan Bullock wrote: I am running genion neutralize the charge of my protein. I got an error that wasn't very descriptive Fatal error: Cannot read from input The part of the script that called this program read as follows /usr/l

Re: [gmx-users] Genion cannot read from input

2013-04-27 Thread Justin Lemkul
On 4/27/13 11:21 AM, Nathan Bullock wrote: I am running genion neutralize the charge of my protein. I got an error that wasn't very descriptive Fatal error: Cannot read from input The part of the script that called this program read as follows /usr/local/gromacs-4.5.5/bin/grompp -f em.mdp -c fi

[gmx-users] Genion cannot read from input

2013-04-27 Thread Nathan Bullock
I am running genion neutralize the charge of my protein. I got an error that wasn't very descriptive Fatal error: Cannot read from input The part of the script that called this program read as follows /usr/local/gromacs-4.5.5/bin/grompp -f em.mdp -c filled_box.pdb -p protein_mod_filled5.top -o ions

RE: [gmx-users] genion doesn't work in 4.6.1?

2013-04-05 Thread Abhishek Acharya
cience, into this forum. With Regards Abhishek Acharya > > > From: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] on > behalf of Abhishek Acharya [aacha...@iitk.ac.in] > Sent: Friday, April 05, 2013 9:35 PM > To: Discussion list for GROMACS use

RE: [gmx-users] genion doesn't work in 4.6.1?

2013-04-05 Thread Luis Felipe Pineda de Castro
ndian "science" Best, Felipe From: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] on behalf of Abhishek Acharya [aacha...@iitk.ac.in] Sent: Friday, April 05, 2013 9:35 PM To: Discussion list for GROMACS users Subject: Re: [gm

Re: [gmx-users] genion doesn't work in 4.6.1?

2013-04-05 Thread Abhishek Acharya
> On 04/05/2013 09:09 PM, Warren Gallin wrote: >> Don't you have to give tyne name of the positive and negative ions that >> will be added by genion? Otherwise, how would it know? >> >> Warren Gallin > > No, it doesn't work either: > > genion_mpi -s ions.tpr -o solv_ions.gro -p complex.top -neutra

Re: [gmx-users] genion doesn't work in 4.6.1?

2013-04-05 Thread Justin Lemkul
On Fri, Apr 5, 2013 at 3:09 PM, Warren Gallin wrote: > Don't you have to give tyne name of the positive and negative ions that > will be added by genion? Otherwise, how would it know? > > Default names are used (see genion -h), and the default names are correct for all force fields as of version

Re: [gmx-users] genion doesn't work in 4.6.1?

2013-04-05 Thread Albert
On 04/05/2013 09:09 PM, Warren Gallin wrote: Don't you have to give tyne name of the positive and negative ions that will be added by genion? Otherwise, how would it know? Warren Gallin No, it doesn't work either: genion_mpi -s ions.tpr -o solv_ions.gro -p complex.top -neutral -norandom -p

Re: [gmx-users] genion doesn't work in 4.6.1?

2013-04-05 Thread Justin Lemkul
On Fri, Apr 5, 2013 at 3:05 PM, Albert wrote: > Hello: > > I found that there would be some problem of genion in 4.6.1: > > grompp_mpi -f em.mdp -c solv.gro -p complex.top -o ions.tpr > > NOTE 1 [file complex.top, line 44008]: > System has non-zero total charge: -11.00 > Total charge sho

Re: [gmx-users] genion doesn't work in 4.6.1?

2013-04-05 Thread Justin Lemkul
On Fri, Apr 5, 2013 at 3:05 PM, Albert wrote: > Hello: > > I found that there would be some problem of genion in 4.6.1: > > grompp_mpi -f em.mdp -c solv.gro -p complex.top -o ions.tpr > > NOTE 1 [file complex.top, line 44008]: > System has non-zero total charge: -11.00 > Total charge sho

Re: [gmx-users] genion doesn't work in 4.6.1?

2013-04-05 Thread Warren Gallin
Don't you have to give tyne name of the positive and negative ions that will be added by genion? Otherwise, how would it know? Warren Gallin On 2013-04-05, at 1:05 PM, Albert wrote: > Hello: > > I found that there would be some problem of genion in 4.6.1: > > grompp_mpi -f em.mdp -c solv.gr

[gmx-users] genion doesn't work in 4.6.1?

2013-04-05 Thread Albert
Hello: I found that there would be some problem of genion in 4.6.1: grompp_mpi -f em.mdp -c solv.gro -p complex.top -o ions.tpr NOTE 1 [file complex.top, line 44008]: System has non-zero total charge: -11.00 Total charge should normally be an integer. See http://www.gromacs.org/Docum

[gmx-users] genion error { problem only with PROTEIN-LIGAND complex} - No line with moleculetype 'SOL' found the [ molecules ] section of file ‘topol.top’

2013-02-08 Thread Ashima
Dear All, I need to add 6 CL ions to make the system neutral but when I type the command line genion -s ions.tpr -o solv_ions.gro -p topol.top -pname NA -nname CL -nn 6 I am able to add the ion easily to a protein or protein-peptide complex by using the above command line where I replace

[gmx-users] genion error - No line with moleculetype 'SOL' found the [ molecules ] section of file ‘topol.top’ [only with PROTEIN-LIGAND complex]

2013-02-08 Thread Ashima
Dear All, I need to add 6 CL ions to make the system neutral but when I type the command line genion -s ions.tpr -o solv_ions.gro -p topol.top -pname NA -nname CL -nn 6 I am able to add the ion easily to a protein or protein-peptide complex by using the above command line where I replace

[gmx-users] genion error - No line with moleculetype 'SOL' found the [ molecules ] section of file ‘topol.top’

2013-02-08 Thread Ashima
Dear All, I need to add 6 CL ions to make the system neutral but when I type the command line genion -s ions.tpr -o solv_ions.gro -p topol.top -pname NA -nname CL -nn 6 and choose group 15 (SOL), the error is Fatal error: No line with moleculetype 'SOL' found the [ molecules ] section of fi

Re: [gmx-users] genion problem with divalent ions

2012-09-09 Thread Leandro Bortot
Hi there, In addition to setting -pq to 2 as Justin said, you can set -conc to a very low value such as 0.1 for neutralizing the system. One more thing... genion apparently uses the box volume when determining how many ions are necessary for getting a given ionic strength. However,

Re: [gmx-users] genion problem with divalent ions

2012-09-09 Thread Justin Lemkul
On 9/9/12 7:53 AM, Matthias Ernst wrote: Hi, I just encountered (another) strange thing with genion. I wanted to neutralized a system by adding Mg-Ions (charge in ions.itp: 2, forcefield amber99sb), so I expected half of the system charge (-52) to be added as ions. BUT it added the same amount

[gmx-users] genion problem with divalent ions

2012-09-09 Thread Matthias Ernst
Hi, I just encountered (another) strange thing with genion. I wanted to neutralized a system by adding Mg-Ions (charge in ions.itp: 2, forcefield amber99sb), so I expected half of the system charge (-52) to be added as ions. BUT it added the same amount, i.e. 52 Mg ions instead of 26 which wou

Re: [gmx-users] genion

2012-03-31 Thread Justin A. Lemkul
e Março de 2012 14:25 *Assunto:* Re: [gmx-users] genion Eduardo Oliveira wrote: > Hi all, > > I' ve stumbled upon a problem when i used genion to prepare my files. Here's the comand line: > > genion_d -s sistema_stpr.tpr -o sistema_ions.gro -nname CL- -nn 3 -nq

Re: [gmx-users] genion

2012-03-31 Thread Justin A. Lemkul
Eduardo Oliveira wrote: Hi all, I' ve stumbled upon a problem when i used genion to prepare my files. Here's the comand line: genion_d -s sistema_stpr.tpr -o sistema_ions.gro -nname CL- -nn 3 -nq -1 -p sistema.toI I've got this in return: Fatal error: Your solvent group size (73282) is

[gmx-users] genion

2012-03-31 Thread Eduardo Oliveira
Hi all, I' ve stumbled upon a problem when i used genion to prepare my files. Here's the comand line: genion_d -s sistema_stpr.tpr -o sistema_ions.gro -nname CL- -nn 3 -nq -1 -p sistema.toI I've got this in return: Fatal error: Your solvent group size (73282) is not a multiple of 19. How to

Re: [gmx-users] genion -pname no NA molecule (sodium) type

2012-03-16 Thread Mark Abraham
On 17/03/2012 4:57 PM, Wholly Peach wrote: Will you please tell me where can I find the ions.itp file? http://www.gromacs.org/Documentation/Terminology/Environment_Variables Mark -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please sea

Re: [gmx-users] genion -pname no NA molecule (sodium) type

2012-03-16 Thread Wholly Peach
Will you please tell me where can I find the ions.itp file?   I am looking forward to getting your reply.   Wholly From: Mark Abraham To: Discussion list for GROMACS users Sent: Saturday, 17 March 2012 3:55 PM Subject: Re: [gmx-users] genion -pname no NA

Re: [gmx-users] genion -pname no NA molecule (sodium) type

2012-03-16 Thread Mark Abraham
On 17/03/2012 4:48 PM, Wholly Peach wrote: Dear All, I use the following to neutralize the charge "genion -s ions.tpr -o 1AKI_solv_ions.gro -p topol.top -pname NA -nname CL -nn 8", Then I do the following "grompp -f minim.mdp -c 1AKI_solv_ions.gro -p topol.top -o em.tpr"". However the

[gmx-users] genion -pname no NA molecule (sodium) type

2012-03-16 Thread Wholly Peach
 Dear All, I use the following to neutralize the charge "genion -s ions.tpr -o 1AKI_solv_ions.gro -p topol.top -pname NA -nname CL -nn 8", Then I do the following "grompp -f minim.mdp -c 1AKI_solv_ions.gro -p topol.top -o em.tpr"". However the feedback says there is no Na (sodium) moleculetype. I

[gmx-users] genion writes a log file?

2010-11-18 Thread Justin A. Lemkul
Hi, When using versions 4.5.2 and 4.5.3, I have noticed that genion now outputs a file called "genion.log" (or "genion_node0.log" with 4.5.2) that appears to contain the header information normally found in an md.log file produced by mdrun (see below). Is this normal? There really isn't any

[gmx-users] genion

2010-07-21 Thread Neha Gandhi
Hi List, I was just wondering what is the relationship between the molarity and the number of molecules. I have a system with 20 atoms of ligand, 761 water molecules, the box volume of 2.86 nanometer cube. I added 0.5 M NaCl salt concentration using genion. I need to add. Genion put 7 Na+ and 7 Cl

Re: [gmx-users] genion doubt

2010-05-15 Thread Justin A. Lemkul
Bharath.K. Chakravarthi wrote: HI I'm not sure about this kind of error but as per my knowledge the problem might be with the parameter cpp in mdp file. In my mdp file it is as follows The problem is unrelated to cpp, since the error complained about the use of "warnings" in the .mdp

Re: [gmx-users] genion doubt

2010-05-15 Thread Justin A. Lemkul
manjula kasinathan wrote: hi snip Thank u. as per ur suggestion i added sodium ions to my protein. then i carried out mdrun for energy minimization nsteps 3000 but i runs only for 1108 steps. then i went to position restrain steps the command i gave grompp -f pr.mdp -c em

Re: [gmx-users] genion doubt

2010-05-15 Thread Bharath.K. Chakravarthi
HI I'm not sure about this kind of error but as per my knowledge the problem might be with the parameter cpp in mdp file. In my mdp file it is as follows title = ecoli_nusG cpp = /usr/bin/cpp define = -DPOSRES constraints = all-bonds integrator = md dt = 0.001 ; ps ! nsteps = 10 ; total 10

[gmx-users] genion doubt

2010-05-15 Thread manjula kasinathan
hi snip Thank u. as per ur suggestion i added sodium ions to my protein. then i carried out mdrun for energy minimization nsteps 3000 but i runs only for 1108 steps. then i went to position restrain steps the command i gave grompp -f pr.mdp -c em.gro -p top.top -o pr.tpr i got a

Re: [gmx-users] genion: can you place the ions at one location in the box?

2010-03-29 Thread Justin A. Lemkul
Justin A. Lemkul wrote: Smiruthi Ramasubramanian wrote: Justin, I tried to use an index file with atom numbers of the water molecules in the required volume. But when I try to ionize the box, I get this error: Fatal error: The solvent group SOL is not continuous: index[717]=21689, ind

Re: [gmx-users] genion: can you place the ions at one location in the box?

2010-03-29 Thread Justin A. Lemkul
Smiruthi Ramasubramanian wrote: Justin, I tried to use an index file with atom numbers of the water molecules in the required volume. But when I try to ionize the box, I get this error: Fatal error: The solvent group SOL is not continuous: index[717]=21689, index[718]=21696 This makes s

Re: [gmx-users] genion: can you place the ions at one location in the box?

2010-03-29 Thread Smiruthi Ramasubramanian
Justin, I tried to use an index file with atom numbers of the water molecules in the required volume. But when I try to ionize the box, I get this error: Fatal error: The solvent group SOL is not continuous: index[717]=21689, index[718]=21696 This makes sense, but is there a way get around this

Re: [gmx-users] genion: can you place the ions at one location in the box?

2010-03-25 Thread Justin A. Lemkul
Smiruthi Ramasubramanian wrote: I am simulating an ion channel and require help placing a concentration of ions at one side of the box (below the membrane). "Sidedness" is irrelevant in a periodic system; the water layers are continuous with PBC. SETUP: membrane + embedded protein in

[gmx-users] genion: can you place the ions at one location in the box?

2010-03-25 Thread Smiruthi Ramasubramanian
I am simulating an ion channel and require help placing a concentration of ions at one side of the box (below the membrane). SETUP: membrane + embedded protein in the center of the cube. I tried to generate a smaller cube (of the required volume below the cube), solvate, add ions and then *cat* w

Re: [gmx-users] genion

2009-11-15 Thread Justin A. Lemkul
leila karami wrote: Hi total charge of my system (pr-dna) is (-10). which is better in genion command: 1) 10 Na or 2) 20 NA + 10 Cl http://www.gromacs.org/Documentation/Gromacs_Utilities/genion Additional ions model increased ionic strength. Whether or not that is "better" is up to yo

[gmx-users] genion

2009-11-15 Thread leila karami
dear mark my object is neutralization of system. please guide more -- gmx-users mailing listgmx-users@gromacs.org 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

Re: [gmx-users] genion

2009-11-15 Thread Mark Abraham
leila karami wrote: Hi total charge of my system (pr-dna) is (-10). which is better in genion command: 1) 10 Na or 2) 20 NA + 10 Cl Only you can answer that, because only you know your objective. Mark -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/list

[gmx-users] genion

2009-11-15 Thread leila karami
Hi total charge of my system (pr-dna) is (-10). which is better in genion command: 1) 10 Na or 2) 20 NA + 10 Cl -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please

Re: [gmx-users] genion question

2009-08-25 Thread Justin A. Lemkul
Paymon Pirzadeh wrote: Well, if I have understood your words properly, apparently there is a problem in the first residue in amino terminal which I treated the NH2 as a zwitterionic piece(+1), and gave the side chain -1 charge. I thought it should add up to zero, but apparently it does not:

Re: [gmx-users] genion question

2009-08-25 Thread Paymon Pirzadeh
Well, if I have understood your words properly, apparently there is a problem in the first residue in amino terminal which I treated the NH2 as a zwitterionic piece(+1), and gave the side chain -1 charge. I thought it should add up to zero, but apparently it does not: 1 opls_287 1ASP

Re: [gmx-users] genion question

2009-08-25 Thread Justin A. Lemkul
Paymon Pirzadeh wrote: Well, when I look into my .top file, almost no where I see closeness to an integer (as close as you said 0.9998). It gets up to -1.05 or -0.91 sometimes. But I can not see a clear trend on where things go wrong. My protein has only 37 amino acids with 451 atoms (including

Re: [gmx-users] genion question

2009-08-25 Thread Paymon Pirzadeh
Well, when I look into my .top file, almost no where I see closeness to an integer (as close as you said 0.9998). It gets up to -1.05 or -0.91 sometimes. But I can not see a clear trend on where things go wrong. My protein has only 37 amino acids with 451 atoms (including hydrogens). I am not sure

Re: [gmx-users] genion question

2009-08-25 Thread Justin A. Lemkul
Paymon Pirzadeh wrote: Thanks for the tips! But I still have problem with my added ion. Justin had told me that ion names are force-field specific! Since I am using OPLSaa, I checked the itp file for Na ion name, and it was Na+. I used this name with the No, it's not. It's NA+. From ions.it

Re: [gmx-users] genion question

2009-08-25 Thread Paymon Pirzadeh
Thanks for the tips! But I still have problem with my added ion. Justin had told me that ion names are force-field specific! Since I am using OPLSaa, I checked the itp file for Na ion name, and it was Na+. I used this name with the -pname switch, but again during the grompp I ran into the following

Re: [gmx-users] genion question

2009-08-25 Thread Mark Abraham
Paymon Pirzadeh wrote: Can I change the charge of Na ion added from +1 to +0.99 to cancel the negative charge exactly? Does that hurt science or simulation? It's probably irrelevant. The representation of "decimal" numbers on computers can be inexact, such that things like 0.02 added to itself

Re: [gmx-users] genion question

2009-08-25 Thread Paymon Pirzadeh
Can I change the charge of Na ion added from +1 to +0.99 to cancel the negative charge exactly? Does that hurt science or simulation? Payman On Tue, 2009-08-25 at 17:22 -0400, Justin A. Lemkul wrote: > > Paymon Pirzadeh wrote: > > I realized sth about the previous error message I posted. In th

Re: [gmx-users] genion command with PDB files for amber03 force field

2009-03-09 Thread Tsjerk Wassenaar
Hi, To set things straight a bit. >> "If you are using ion-related GROMACS tools, such as genion, you will need >> to enter the AMBER ion definition to the ions.itp file in the "top" >> directory of the GROMACS distribution." > > This isn't too tough, as all the parameters for the ions are includ

Re: [gmx-users] genion command with PDB files for amber03 force field

2009-03-08 Thread Joshua Ballanco
Andrew, Just ran into this myself. If you read the notes on http://chemistry.csulb.edu/ffamber/#usage closely, you'll find this little tid-bit: "If you are using ion-related GROMACS tools, such as genion, you will need to enter the AMBER ion definition to the ions.itp file in the "top" di

Re: [gmx-users] genion command with PDB files for amber03 force field

2009-03-05 Thread Justin A. Lemkul
drugdesign wrote: Dear gromacs users, I am using genion command to add ions to .tpr file which was made for protein by use of amber03 force field. In intial pdb file N-terminal amino acid is NHIE, but genion changes it to NHI, which isnot recognized by next pdb2gmx. PDB format dictates 3 ch

[gmx-users] genion command with PDB files for amber03 force field

2009-03-05 Thread drugdesign
Dear gromacs users, I am using genion command to add ions to .tpr file which was made for protein by use of amber03 force field. In intial pdb file N-terminal amino acid is NHIE, but genion changes it to NHI, which isnot recognized by next pdb2gmx. This is not a big problem , but it looks is a l

[gmx-users] genion and a few other problems

2008-12-27 Thread Manik Mayur
Hi all, I have a few probIems/confusions and will be extremely thankful for any help. I added some "Cl- ions" by replacing water using "genion (and -nn)". It generated a "topol.top" file with Cl listed in [ molecules ]. Now when I try to run grompp, it gives an error saying "no such molecule type

[gmx-users] genion problem with large systems

2008-09-09 Thread Justin A. Lemkul
Hi all, I'm experiencing a strange problem with genion that I have never seen before. I am preparing some membrane systems, many of which involve a large number of charged lipids (several hundred). When using lipids with PC headgroups (no net charge), the following works fine: genion -s i

Re: [gmx-users] Genion failing in parallel runs.

2008-08-26 Thread Florian Haberl
Hi, On Tuesday, 26. August 2008, Bhanu wrote: > Hi, > I'm using Gromacs 3.3.3, installed with double precision. When I run genion > with this command: > > genion_d -s *.tpr -o *.gro -g genion.log -pot pot.pdb -p *.top -neutral > > the following error message was generated: > --

Re: [gmx-users] Genion failing in parallel runs.

2008-08-26 Thread Justin A. Lemkul
Bhanu wrote: Hi, I'm using Gromacs 3.3.3, installed with double precision. When I run genion with this command: genion_d -s *.tpr -o *.gro -g genion.log -pot pot.pdb -p *.top -neutral the following error message was generated: --

[gmx-users] Genion failing in parallel runs.

2008-08-26 Thread Bhanu
Hi, I'm using Gromacs 3.3.3, installed with double precision. When I run genion with this command: genion_d -s *.tpr -o *.gro -g genion.log -pot pot.pdb -p *.top -neutral the following error message was generated: Program g

[gmx-users] genion strips +, -, digits. Bug or feature?

2008-07-23 Thread Hans Martin Senn
Dear all I have added the new AMBER ion definitions from Joung and Cheatham to my ffamber* force field files. Following the nomenclature these ions have in AMBER, I use the same name for the atom name, the residue name, and the molecule type, which includes the charge. E.g., I have In ion

[gmx-users] genion strips +, -, digits. Bug or feature?

2008-07-19 Thread Hans Martin Senn
I sent the same message already yesterday but it doesn't seem to have come through. Sorry if it did, and it was just me who for whatever reason didn't get it. Dear all I have added the new AMBER ion definitions from Joung and Cheatham to my ffamber* force field files. Following the nomenc

Re: [gmx-users] genion

2008-06-18 Thread Tsjerk Wassenaar
Hi sh-karbalaee, Well, running an MD simulation has nothing to do with genion per se. But it's better not to have a net chare in your system. Therefore you'd best add an ion (at the least) to your system with unit charge. Whether you do that using genion or do it by hand, is entirely up to you. Ch

[gmx-users] genion

2008-06-18 Thread shahrbanoo karbalaee
Dear justin thanks for reply about g_dist. about genion for peptides.Do I must use command genion for MD peptides with one charge(13 aa)?Is it necessary for md? best sh-karbalaee ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org

Re: [gmx-users] genion problem: is replacing 84 water molecules instead of 200

2008-04-03 Thread David van der Spoel
Subhrangshu Supakar wrote: Hi All !! I have defined my own system which contains 200 water molecules, which I can see in the .gro file. Then I created the .tpr file for 200 steps to steep. In the confout file I can again see the 200 water molecules. But When I am trying to replace the 200 wate

[gmx-users] genion problem: is replacing 84 water molecules instead of 200

2008-04-03 Thread Subhrangshu Supakar
Hi All !! I have defined my own system which contains 200 water molecules, which I can see in the .gro file. Then I created the .tpr file for 200 steps to steep. In the confout file I can again see the 200 water molecules. But When I am trying to replace the 200 water molecules with 200 Na+ GROMAC

Re: [gmx-users] Genion: specify minimum distance between counter-ion and solute

2007-11-06 Thread Mark Abraham
Narayanan Veeraraghavan wrote: Hi, I was wondering if there was a way to specify a minimum distance between the counter-ion and solute, using Genion. There is a rmin option, buts thats for the dist between the ions. But, I am interested in the distance between ions and solute. Would appreciate an

[gmx-users] Genion: specify minimum distance between counter-ion and solute

2007-11-06 Thread Narayanan Veeraraghavan
Hi, I was wondering if there was a way to specify a minimum distance between the counter-ion and solute, using Genion. There is a rmin option, buts thats for the dist between the ions. But, I am interested in the distance between ions and solute. Would appreciate any help. Thanks, Narayan

RE: [gmx-users] genion - hangs up

2007-07-26 Thread Rigden, LucianeVMello
Building Tel:(+44) 151 795 5140 FAX:(+44) 151 794 5130 University of Liverpool Crown St.,Liverpool L69 7ZB, U.K. -Original Message- From: [EMAIL PROTECTED] on behalf of Tsjerk Wassenaar Sent: Thu 7/19/2007 7:53 AM To: Discussion list for GROMACS users Subject: Re: [gmx-users] genion - hangs

Re: [gmx-users] genion - hangs up

2007-07-18 Thread Tsjerk Wassenaar
Hi Luciane, David van der Spoel has replied saying: there is a problem compiling gromacs using gcc 4.x that may cause this. see a recent bugzilla entry about it. Did you compile using the 4.x gcc? In that case, recompile with an older version (3.x). That solution was suggested in the same di

[gmx-users] genion - hangs up

2007-07-18 Thread Rigden, LucianeVMello
Dear all, I am having problem with genion. I have already search the gmx archive, and got the following message: While using genion with -np 3, i get 2 lines of Replacing sol. xxx with Na and it hangs up. With -np 1 it works fine and replaces 1 SOL with Na. I thought it might be my system, so

Re: [gmx-users] genion - why shall we neutralize the charges?

2007-06-11 Thread Yang Ye
list. It is not easy to get to know what you are asking. Regards, Yang Ye - Original Message From: Elias santos <[EMAIL PROTECTED]> To: gmx-users Sent: Tuesday, June 12, 2007 1:34:39 AM Subject: [gmx-users] genion Hi!! I want know because we must neutralize the total load

Re: [gmx-users] genion

2007-06-11 Thread David van der Spoel
Elias santos wrote: Hi!! I want know because we must neutralize the total load of the system (protein + water) into of the box with genion. You should neutralize the system and you also may want to add extra ions to get up to physiological salt concentration. --

[gmx-users] genion

2007-06-11 Thread Elias santos
Hi!! I want know because we must neutralize the total load of the system (protein + water) into of the box with genion. ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at htt

RE: [gmx-users] genion

2007-06-10 Thread Dallas B. Warren
riginal Message- From: [EMAIL PROTECTED] on behalf of Elias santos Sent: Sun 6/10/2007 4:21 AM To: gmx-users@gromacs.org Subject: [gmx-users] genion Hi!! I want know because we must neutralize the total load of the system (protein + water) in the interior of the box wit

[gmx-users] genion

2007-06-09 Thread Elias santos
Hi!! I want know because we must neutralize the total load of the system (protein + water) in the interior of the box with genion. Elias ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please search t

Re: [gmx-users] genion halted

2007-03-29 Thread andrea spitaleri
Hi yes I do use genion and g_rdf properly with 4.0.2 Regards and Sampo Karkola wrote: > Hi, > > so you got genion working with 4.0.2? > > Sampo > > andrea spitaleri wrote: >> gcc-4.0.2 seems ok for gmx-3.3.1 >> I have also submitted in bugzilla a similar problem with g_rdf >> >> >> Regards >>

Re: [gmx-users] genion halted

2007-03-29 Thread Sampo Karkola
Hi, so you got genion working with 4.0.2? Sampo andrea spitaleri wrote: gcc-4.0.2 seems ok for gmx-3.3.1 I have also submitted in bugzilla a similar problem with g_rdf Regards andrea Sampo Karkola wrote: Hi, the same one, gcc-4.1. Have you tried with other versions? Sampo andrea spital

Re: [gmx-users] genion halted

2007-03-29 Thread andrea spitaleri
gcc-4.0.2 seems ok for gmx-3.3.1 I have also submitted in bugzilla a similar problem with g_rdf Regards andrea Sampo Karkola wrote: > Hi, > > the same one, gcc-4.1. Have you tried with other versions? > > Sampo > > andrea spitaleri wrote: >> Hi >> which gcc have you used for compiling gromac

Re: [gmx-users] genion halted

2007-03-29 Thread Sampo Karkola
Hi, the same one, gcc-4.1. Have you tried with other versions? Sampo andrea spitaleri wrote: Hi which gcc have you used for compiling gromacs? I have the same problem in one pc in which I have build gromacs using gcc-4.1 Regards andrea Sampo Karkola wrote: Dear list, the genion program g

Re: [gmx-users] genion halted

2007-03-29 Thread andrea spitaleri
Hi which gcc have you used for compiling gromacs? I have the same problem in one pc in which I have build gromacs using gcc-4.1 Regards andrea Sampo Karkola wrote: > Dear list, > > the genion program gets halted when I'm trying to replace the solvent > molecules with ions. The command > > gen

[gmx-users] genion halted

2007-03-29 Thread Sampo Karkola
Dear list, the genion program gets halted when I'm trying to replace the solvent molecules with ions. The command genion -s protein_in_waterbox.tpr -o protein_out.gro -p topology.top -pname NA+ -nname CL- -random -neutral -conc 0.1 asks for the group of solvent molecules and when I tell him

Re: [[gmx-users] genion error!]

2007-03-22 Thread Pär Bjelkmar
Maybe you forgot to delete the waters that were removed in the ionization process? /Pär > Dear users: > I want to prepare protein simulation on 0.15mol/l salt > concentration,I use > genion command to add ions,After genion,I change cln.top file last > line to NA+ > 61,CL- 41. >When I

RE: [gmx-users] genion error!

2007-02-28 Thread Triguero, Luciano O
Memorial Drive, Room 146 P.O Box 249118 Coral Gables, FL 33124-0431 Cellular: 305-904-2419 Office: 305-284-3938 -Original Message- From: [EMAIL PROTECTED] on behalf of Qiang Zhong Sent: Wed 2/28/2007 3:23 AM To: gmx-users@gromacs.org Subject: [gmx-users] genion error! Dear users: I want

Re: [gmx-users] genion error!

2007-02-28 Thread Yang Ye
On 2/28/2007 7:52 PM, Mark Abraham wrote: Qiang Zhong wrote: Dear users: I want to prepare protein simulation on 0.15mol/l salt concentration,I use genion command to add ions,After genion,I change cln.top file last line to NA+ 61,CL- 41. When I execute grompp ,there are warning : War

Re: [gmx-users] genion error!

2007-02-28 Thread Mark Abraham
Qiang Zhong wrote: Dear users: I want to prepare protein simulation on 0.15mol/l salt concentration,I use genion command to add ions,After genion,I change cln.top file last line to NA+ 61,CL- 41. When I execute grompp ,there are warning : Warning: atom names in cln.top and cln_b4em.gro d

Re: [gmx-users] genion error!

2007-02-28 Thread Erik Marklund
Something seems wrong with your topology. It seems grompp gets confused about CL in one file being Na in the other. Maybe you added the last lines to cln.top in the wrong order. It doesn't necessarily cause the LINCS error later on, but suggests that your topology may need some more work.

[gmx-users] genion error!

2007-02-28 Thread Qiang Zhong
Dear users: I want to prepare protein simulation on 0.15mol/l salt concentration,I use genion command to add ions,After genion,I change cln.top file last line to NA+ 61,CL- 41. When I execute grompp ,there are warning : Warning: atom names in cln.top and cln_b4em.gro don't match (CL - Na)

Re: [gmx-users] genion question

2006-12-11 Thread Mark Abraham
> Hi Guys, > > > > I have a question. I am using genion to add two ZN2+ ions to a protein > that > already has two. I thought I would just use genion with –np 2. However > it > seemed to recognize that there were already two in the file and added > none. > So instead I used the parameter –np 4.

[gmx-users] genion question

2006-12-11 Thread Samuel C Flores
Hi Guys, I have a question. I am using genion to add two ZN2+ ions to a protein that already has two. I thought I would just use genion with –np 2. However it seemed to recognize that there were already two in the file and added none. So instead I used the parameter –np 4. This seemed to a

Re: [gmx-users] GENION: Ions inside the protein

2006-10-10 Thread Tsjerk Wassenaar
Hi, The error message tries to make clear that genion can only replace water molecules in an index group without gaps. It has to be continuous. So the point that was made earlier on the list was that the water molecules should be deleted in the .gro file. Alternatively, they can also be moved to

Re: [gmx-users] GENION: Ions inside the protein

2006-10-10 Thread Viswanadham Sridhara
you might have deleted just one atom of water molecule. Ex: If you are using SPC water , you have to delete all 3 atoms in index file.     Best way is to make changes in .gro file.. Renumber the resulting .gro file using genconf, make a new index file and start running.  On 10/10/06, Una Bjarnadott

[gmx-users] GENION: Ions inside the protein

2006-10-10 Thread Una Bjarnadottir
Dear users, In some of my simulations after neutralizing the systems with genion 2 of the ions were inside the active site cleft of the protein. In the discussion on the list it was commented to either delete the solvent within the protein which I did I deleted the water molecule (11798) and

  1   2   >