Re: [gmx-users] Energy minimization has stopped....

2013-11-05 Thread Justin Lemkul
On 11/5/13 7:31 AM, Kalyanashis Jana wrote: Dear Justin, Can you please tell me, how I can solve this problem?? If I will change the coordinate of atom 15461, will it help me? But you know, I did this step changing the position of drug molecule and I got same error. You should first do wh

Re: [gmx-users] Energy minimization has stopped....

2013-11-05 Thread Kalyanashis Jana
Dear Justin, Can you please tell me, how I can solve this problem?? If I will change the coordinate of atom 15461, will it help me? But you know, I did this step changing the position of drug molecule and I got same error. On Tue, Nov 5, 2013 at 5:44 PM, Justin Lemkul wrote: > > > On 11/5/13

Re: [gmx-users] Energy minimization has stopped....

2013-11-05 Thread Justin Lemkul
On 11/5/13 6:28 AM, Kalyanashis Jana wrote: Hi, Whenever I am trying to do position retrained MD run, It has been stopped at middle of the MD run. I have given the following error. Can you please suggest me something to resolve this error? Energy minimization has stopped, but the forces have

Re: [gmx-users] Energy minimization has stopped....

2013-11-05 Thread jkrieger
What does your curve look like? What parameters are you using in the mdp? How big is your system and what kind of molecules are in there? Providing this kind of information would help people work out what the problem is. Then again it may be ok that the minimisation has converged without reaching

[gmx-users] Energy minimization has stopped....

2013-11-05 Thread Kalyanashis Jana
Hi, Whenever I am trying to do position retrained MD run, It has been stopped at middle of the MD run. I have given the following error. Can you please suggest me something to resolve this error? Energy minimization has stopped, but the forces havenot converged to the requested precision Fmax < 1

Re: [gmx-users] energy minimization problem

2013-11-03 Thread Justin Lemkul
On 11/3/13 6:20 AM, kiana moghaddam wrote: Dear Justin Further to your previous email, I want to calculate free energy for DNA-ligand interaction. According to your answer, for more sensitive calculations like free energy simulations and normal modes, emtolshould be lower than 1. As I understa

Re: [gmx-users] energy minimization problem

2013-11-02 Thread Justin Lemkul
On 11/2/13 3:38 PM, kiana moghaddam wrote: Hi GMX Users I 'm running MD Simulation for DNA-ligand interaction. I ran energy minimization with emtol=10, 100 and 250 in 1 steps (steepest descent followed by the conjugate gradient algorithm). The lowest value for energy obtained for emtol=

[gmx-users] energy minimization problem

2013-11-02 Thread kiana moghaddam
Hi GMX Users I 'm running MD Simulation for DNA-ligand interaction. I ran energy minimization with emtol=10, 100 and 250 in 1 steps (steepest descent followed by the conjugate gradient algorithm). The lowest value for energy obtained for emtol=10.  I have some questions about energy minimiz

Re: [gmx-users] energy minimization

2013-09-23 Thread Mark Abraham
On Sep 23, 2013 9:08 AM, "marzieh dehghan" wrote: > > Hi every body > in order to protein- ligand docking, energy minimization was done by > GROMACS. I did the following steps for insulin pdb file: > > 1- pdb2gmx -ignh -f 3inc.pdb -o test.pdb -p topol.top -water spce > 2- grompp -f em.mdp -c test.

Re: [gmx-users] energy minimization

2013-09-23 Thread Justin Lemkul
On 9/23/13 3:08 AM, marzieh dehghan wrote: Hi every body in order to protein- ligand docking, energy minimization was done by GROMACS. I did the following steps for insulin pdb file: 1- pdb2gmx -ignh -f 3inc.pdb -o test.pdb -p topol.top -water spce 2- grompp -f em.mdp -c test.pdb -p topol.top

[gmx-users] energy minimization

2013-09-23 Thread marzieh dehghan
Hi every body in order to protein- ligand docking, energy minimization was done by GROMACS. I did the following steps for insulin pdb file: 1- pdb2gmx -ignh -f 3inc.pdb -o test.pdb -p topol.top -water spce 2- grompp -f em.mdp -c test.pdb -p topol.top -o em.tpr 3-mdrun -v -deffnm em 4- editconf -f

Re: [gmx-users] Energy minimization failure

2013-03-19 Thread Justin Lemkul
On Tue, Mar 19, 2013 at 9:16 AM, Jeff Woodford < jwoodf...@missouriwestern.edu> wrote: > Greetings, > I am relatively new to MD, and I am attempting to simulate a metal-organic > framework but I don't seem to be even able to get past the energy > minimization phase. I could appreciate any insight

[gmx-users] Energy minimization failure

2013-03-19 Thread Jeff Woodford
Greetings, I am relatively new to MD, and I am attempting to simulate a metal-organic framework but I don't seem to be even able to get past the energy minimization phase. I could appreciate any insight into what I might be doing wrong. Here is a brief summary of what I have done: - constructe

Re: [gmx-users] Energy minimization

2013-02-23 Thread Kieu Thu Nguyen
Thank Justin so much ! On Sat, Feb 23, 2013 at 1:50 AM, Justin Lemkul wrote: > > > On 2/21/13 9:10 PM, Kieu Thu Nguyen wrote: > >> Dear all, >> >> I can not distinguish when the Fmax < 10 kJmol-1nm-1, < 100, < 1000 ? is >> appropriate >> >> Can anyone tell me which value is appropriate for li

Re: [gmx-users] Energy minimization

2013-02-22 Thread Justin Lemkul
On 2/21/13 9:10 PM, Kieu Thu Nguyen wrote: Dear all, I can not distinguish when the Fmax < 10 kJmol-1nm-1, < 100, < 1000 ? is appropriate Can anyone tell me which value is appropriate for lipid membrane, for protein in water, for protein in lipid membrane ? There is no absolute rule. Th

RE: [gmx-users] Energy Minimization/ NpT settling problem

2012-12-20 Thread pcl
inity? From: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] on behalf of emmanuelle [e.a.y.m...@student.tudelft.nl] Sent: Thursday, December 20, 2012 2:31 PM To: gmx-users@gromacs.org Subject: [gmx-users] Energy Minimization/ NpT settling problem Hi,

[gmx-users] Energy Minimization/ NpT settling problem

2012-12-20 Thread emmanuelle
Hi, I am trying to simulate a system consisting of 1790 molecules of monoethanolamine (MEA), 203 molecules of CO2 and 18210 molecules of water. After having generated my box and set the density to 1000g/l using editconf, I am facing a problem with the energy minimization step. It does not converge

Re: [gmx-users] Energy minimization

2012-12-18 Thread Peter C. Lai
sion list for GROMACS users > Sent: Sunday, December 16, 2012 11:18 AM > Subject: Re: [gmx-users] Energy minimization > > Probably remove the overlapping lipid then. Once you run MD it will repack. > > On 2012-12-15 09:19:49PM -0800, Shima Arasteh wrote: > > Thanks

Re: [gmx-users] Energy minimization

2012-12-17 Thread Shima Arasteh
d the overlapping atoms and the reason of clashe among a vast number of atoms? Sincerely, Shima - Original Message - From: Peter C. Lai To: Shima Arasteh Cc: Discussion list for GROMACS users Sent: Sunday, December 16, 2012 11:18 AM Subject: Re: [gmx-users] Energy minimization Probably remov

Re: [gmx-users] Energy minimization

2012-12-16 Thread Shima Arasteh
It's printed in terminal, exactly before showing the values of infinite potential and Max Force and . Sincerely, Shima From: Dariush Mohammadyani To: Shima Arasteh Sent: Sunday, December 16, 2012 5:01 PM Subject: Re: [gmx-users] Energy minimiz

Fw: [gmx-users] Energy minimization

2012-12-16 Thread Shima Arasteh
  Sincerely, Shima - Forwarded Message - From: Shima Arasteh To: Dariush Mohammadyani Sent: Sunday, December 16, 2012 5:03 PM Subject: Re: [gmx-users] Energy minimization It's printed in terminal, exactly before showing the values of infinite potential and Max Forc

Re: [gmx-users] Energy minimization

2012-12-15 Thread Peter C. Lai
Cc: > Sent: Sunday, December 16, 2012 8:44 AM > Subject: Re: [gmx-users] Energy minimization > > It depends on what the atom is overlapping with and some conjecture as to > what might be causing the overlap: > > You can always manually move it, either by editing the .gro fil

Re: [gmx-users] Energy minimization

2012-12-15 Thread Shima Arasteh
- From: Peter C. Lai To: Shima Arasteh ; Discussion list for GROMACS users Cc: Sent: Sunday, December 16, 2012 8:44 AM Subject: Re: [gmx-users] Energy minimization It depends on what the atom is overlapping with and some conjecture as to what might be causing the overlap: You can always manually

Re: [gmx-users] Energy minimization

2012-12-15 Thread Peter C. Lai
t; users > Sent: Saturday, September 29, 2012 3:01 PM > Subject: Re: [gmx-users] Energy minimization > > > > On 9/29/12 3:19 AM, Shima Arasteh wrote: > > > > Dear all, > > > > My system contains lipids, protein and water. > > I want to energy mi

Re: [gmx-users] Energy minimization

2012-12-15 Thread Shima Arasteh
Subject: Re: [gmx-users] Energy minimization On 9/29/12 3:19 AM, Shima Arasteh wrote: > > Dear all, > > My system contains lipids, protein and water. > I want to energy minimize it, so ran grompp: > > > # grompp -f em.mdp -c system_solv_ions.gro -p topol.top -o em.tpr &g

[gmx-users] Energy minimization with walls

2012-11-23 Thread harshaljain950
Hi all, I am simulating a box filled only with water and having walls at Z=0 and Z=Z_box. I am confused that whether energy minimization process removes the bad contacts between wall atoms and solvent or not ? any kind of help would be really appreciated -- View this message in context: http

Re: [gmx-users] Energy minimization

2012-09-29 Thread Justin Lemkul
On 9/29/12 3:19 AM, Shima Arasteh wrote: Dear all, My system contains lipids, protein and water. I want to energy minimize it, so ran grompp: # grompp -f em.mdp -c system_solv_ions.gro -p topol.top -o em.tpr and then: # mdrun -v -deffnm em The output is: Steepest Descents: Tolerance

[gmx-users] Energy minimization

2012-09-29 Thread Shima Arasteh
Dear all, My system contains lipids, protein and water. I want to energy minimize it, so ran grompp: # grompp -f em.mdp -c system_solv_ions.gro -p topol.top -o em.tpr and then: # mdrun -v -deffnm em The output is: Steepest Descents:    Tolerance (Fmax)   =  1.0e+03    Number of steps   

Re: [gmx-users] energy minimization output

2012-09-26 Thread Erik Marklund
It seems that your input structure has too severe clashes for energy minimization, that there are problems with your topology, or that your mdp options are not adequate. One possible solution, if your topology is correct, is to first perform an energy minimization using soft-core potentials to g

Re: [gmx-users] energy minimization output

2012-09-26 Thread Archana Sonawani
Dear Erik, I also encountered similar problem. So you mean it cannot be further energy minimized? How can we proceed with NPT if the system is not in equilibrium? On Wed, Sep 26, 2012 at 1:11 PM, Erik Marklund wrote: > > 26 sep 2012 kl. 08.16 skrev Shima Arasteh: > > > Dear all, > > > > My syst

Re: [gmx-users] energy minimization output

2012-09-26 Thread Erik Marklund
26 sep 2012 kl. 08.16 skrev Shima Arasteh: > Dear all, > > My system contains lipids, protein and water. > I want to energy minimize it, so ran grompp: > > > # grompp -f em.mdp -c system_solv_ions.gro -p topol.top -o em.tpr > > and then: > # mdrun -v -deffnm em > > > The output is: > Steepe

[gmx-users] energy minimization output

2012-09-25 Thread Shima Arasteh
Dear all, My system contains lipids, protein and water. I want to energy minimize it, so ran grompp: # grompp -f em.mdp -c system_solv_ions.gro -p topol.top -o em.tpr and then: # mdrun -v -deffnm em The output is: Steepest Descents:    Tolerance (Fmax)   =  1.0e+03    Number of steps    =

Re: [gmx-users] Energy Minimization of Alanin

2012-09-19 Thread Lara Bunte
:12 Mittwoch, 19.September 2012 Betreff: Re: [gmx-users] Energy Minimization of Alanin for EM you can probably ignore this, but note that these are the wrong cutoffs for CHARMM27. (rlist=1.2 rlistlong=1.4 rcoulomb=1.2 rvdw=1.2 rvdw_switch=0.8 and vdwtype=switch). On 2012-09-19 11:04:47AM +0100

Re: [gmx-users] Energy Minimization of Alanin

2012-09-19 Thread Peter C. Lai
for EM you can probably ignore this, but note that these are the wrong cutoffs for CHARMM27. (rlist=1.2 rlistlong=1.4 rcoulomb=1.2 rvdw=1.2 rvdw_switch=0.8 and vdwtype=switch). On 2012-09-19 11:04:47AM +0100, Lara Bunte wrote: > Hello > > I want to do md simulations with the amino acid alanin fo

Re: [gmx-users] Energy Minimization of Alanin

2012-09-19 Thread Justin Lemkul
On 9/19/12 6:04 AM, Lara Bunte wrote: Hello I want to do md simulations with the amino acid alanin for practice. I choose alanin because it is parametrized in my forcefield. I use charmm27. I created with pdb2gmx my topology and I use water model tip3p. I choosed a dodecahedron box with dis

[gmx-users] Energy Minimization of Alanin

2012-09-19 Thread Lara Bunte
Hello I want to do md simulations with the amino acid alanin for practice. I choose alanin because it is parametrized in my forcefield. I use charmm27. I created with pdb2gmx my topology and I use water model tip3p. I choosed a dodecahedron box with distance of 0.5 between the solute and this b

Re: [gmx-users] Energy Minimization - not getting correct lipid area

2012-06-08 Thread Justin A. Lemkul
On 6/8/12 3:15 PM, Hicks, Erica wrote: Hi, I am working through the KALP-15 tutorial and having difficulties getting the correct area per lipid (~71 A^2). I first scaled the lipid by a factor of 4 (perl inflategro.pl system.gro 4 DPPC 14 system_inflated.gro 5 area.dat) then did energy minimiza

[gmx-users] Energy Minimization - not getting correct lipid area

2012-06-08 Thread Hicks, Erica
Hi, I am working through the KALP-15 tutorial and having difficulties getting the correct area per lipid (~71 A^2). I first scaled the lipid by a factor of 4 (perl inflategro.pl system.gro 4 DPPC 14 system_inflated.gro 5 area.dat) then did energy minimization (mdrun -v -deffnm em) and then scal

RE: [gmx-users] Energy minimization of crystal structure with water

2012-03-18 Thread NG HUI WEN
Hi Mark, Thanks - I completely forgot about that! Cheers, Huiwen From: gmx-users-boun...@gromacs.org [mailto:gmx-users-boun...@gromacs.org] On Behalf Of Mark Abraham Sent: Sunday, March 18, 2012 2:03 PM To: Discussion list for GROMACS users Subject: Re: [gmx-users] Energy minimization of

Re: [gmx-users] Energy minimization of crystal structure with water

2012-03-17 Thread Mark Abraham
On 18/03/2012 1:01 PM, NG HUI WEN wrote: Hi all, I have a seemingly simple task that turned quite tricky here. I am trying to perform energy minimization on a crystal structure after adding hydrogen to the protein using pdb2gmx (OPLS). The problem I am facing now is that after energy minimi

[gmx-users] Energy minimization of crystal structure with water

2012-03-17 Thread NG HUI WEN
Hi all, I have a seemingly simple task that turned quite tricky here. I am trying to perform energy minimization on a crystal structure after adding hydrogen to the protein using pdb2gmx (OPLS). The problem I am facing now is that after energy minimization, the protein and water (origina

Re: [gmx-users] Energy minimization Error

2012-02-10 Thread Mark Abraham
On 11/02/2012 4:15 AM, francesca vitalini wrote: Ok now I have tryied to restart it all and the problem seems to be that the system has some overlapping. In fact, no matter what I freeze, water, CA or nothing, I get this error message VERSION 3.3.1 This program is free software; you

Re: [gmx-users] Energy minimization Error

2012-02-10 Thread Mark Abraham
On 11/02/2012 2:48 AM, francesca vitalini wrote: In fact in the reverse transformation I'm feeding the CG structure information. Yes, but you need a source of atomistic intra-residue knowledge about protein structure in order to convert a CG structure to something that is vaguely plausible at

Re: [gmx-users] Energy minimization Error

2012-02-10 Thread francesca vitalini
Ok now I have tryied to restart it all and the problem seems to be that the system has some overlapping. In fact, no matter what I freeze, water, CA or nothing, I get this error message VERSION 3.3.1 This program is free software; you can redistribute it and/or modify it under

Re: [gmx-users] Energy minimization Error

2012-02-10 Thread francesca vitalini
In fact in the reverse transformation I'm feeding the CG structure information. Once I look through VMD to the FG structure I notice that the backbones are not in planes so definitely I need to run some minimization there. In the tutorial they were using simulating annealing, but I don't think I ne

Re: [gmx-users] Energy minimization Error

2012-02-10 Thread Mark Abraham
On 11/02/2012 12:41 AM, francesca vitalini wrote: In order to overcome the problem I tried to fix everything except the backbone (solvent, sidechain and CA, as I want the structure to be maintained). However, if I do then I have problems with the energy minimization as the force on the 15300 is

Re: [gmx-users] Energy minimization Error

2012-02-10 Thread francesca vitalini
In order to overcome the problem I tried to fix everything except the backbone (solvent, sidechain and CA, as I want the structure to be maintained). However, if I do then I have problems with the energy minimization as the force on the 15300 is infinite. Getting Loaded... Reading file EM1-1.tpr,

Re: [gmx-users] Energy minimization Error

2012-02-10 Thread Justin A. Lemkul
francesca vitalini wrote: I achieve Steepest Descents converged to machine precision in 205 steps, but did not reach the requested Fmax < 10. Potential Energy = -1.49131940478719e+07 Maximum force = 1.09664530279637e+06 on atom 1520(parto of the protein) Norm of force = 2.2836

Re: [gmx-users] Energy minimization Error

2012-02-10 Thread francesca vitalini
I achieve Steepest Descents converged to machine precision in 205 steps, but did not reach the requested Fmax < 10. Potential Energy = -1.49131940478719e+07 Maximum force = 1.09664530279637e+06 on atom 1520(parto of the protein) Norm of force = 2.28369808518165e+06 and the system l

Re: [gmx-users] Energy minimization Error

2012-02-09 Thread Justin A. Lemkul
francesca vitalini wrote: Tryed the double precision without any success... Again the same type of error message. Might need to try other minimization algorithm. Any suggestions on that? What is the Fmax you achieve, and what are your intents for this process? Minimizing very low is only

Re: [gmx-users] Energy minimization Error

2012-02-09 Thread francesca vitalini
Tryed the double precision without any success... Again the same type of error message. Might need to try other minimization algorithm. Any suggestions on that? 2012/2/9 Justin A. Lemkul > > > francesca vitalini wrote: > >> Hi all! >> I'm trying to run some energy minimization with gromacs 3.3.1

Re: [gmx-users] Energy minimization Error

2012-02-09 Thread Justin A. Lemkul
francesca vitalini wrote: Hi all! I'm trying to run some energy minimization with gromacs 3.3.1 on a system that I have transformed from cg to aa. After the reverse transformation I solvate the system and I try to minimize just the proteins (my system is a dimer) keeping fixed the water. The

[gmx-users] Energy minimization Error

2012-02-09 Thread francesca vitalini
Hi all! I'm trying to run some energy minimization with gromacs 3.3.1 on a system that I have transformed from cg to aa. After the reverse transformation I solvate the system and I try to minimize just the proteins (my system is a dimer) keeping fixed the water. Then the idea would be to minimize t

Re: [gmx-users] Energy Minimization

2012-01-10 Thread Mark Abraham
On 10/01/2012 10:46 PM, Steven Neumann wrote: Thank you. Imagine I would like to put ligands in further simulations. Should I then copy coordinates of my ligands in smaller box (not to overlap my tube) and then copy both coordinates of ions and ligands plus protein? Any way you want to constr

Re: [gmx-users] Energy Minimization

2012-01-10 Thread Steven Neumann
Thank you. Imagine I would like to put ligands in further simulations. Should I then copy coordinates of my ligands in smaller box (not to overlap my tube) and then copy both coordinates of ions and ligands plus protein? On Tue, Jan 10, 2012 at 11:31 AM, Mark Abraham wrote: > On 10/01/2012 10:13

Re: [gmx-users] Energy Minimization

2012-01-10 Thread Mark Abraham
On 10/01/2012 10:13 PM, Steven Neumann wrote: Thank you. Should I also copy and paste coordinates of my ions or just my protein? The randomly-placed ions will be immaterial for EM. Mark -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Ple

Re: [gmx-users] Energy Minimization

2012-01-10 Thread Steven Neumann
Thank you. Should I also copy and paste coordinates of my ions or just my protein? Steven On Tue, Jan 10, 2012 at 11:03 AM, Mark Abraham wrote: > On 10/01/2012 9:54 PM, Steven Neumann wrote: > >> Dear Gmx Users, >> I am setting up my simulations of carbon tube with protein. I solvated my >> syst

Re: [gmx-users] Energy Minimization

2012-01-10 Thread Mark Abraham
On 10/01/2012 9:54 PM, Steven Neumann wrote: Dear Gmx Users, I am setting up my simulations of carbon tube with protein. I solvated my system, added ions and I would like to run EM of my system. My carbons of the tube in MD will be restrained. In this case should I run EM of my protein in wate

[gmx-users] Energy Minimization

2012-01-10 Thread Steven Neumann
Dear Gmx Users, I am setting up my simulations of carbon tube with protein. I solvated my system, added ions and I would like to run EM of my system. My carbons of the tube in MD will be restrained. In this case should I run EM of my protein in water (and with ions) separately and the copy coordin

Re: [gmx-users] energy minimization error

2011-05-21 Thread David van der Spoel
On 2011-05-21 09.09, sreelakshmi ramesh wrote: Dear all, I am using tabulated potential do simulate nacl in water .when i was doing my energy minimization it gives me the following error.Any help please.. Do you mean to say your table is correct? Check the derivative numerically

[gmx-users] energy minimization error

2011-05-21 Thread sreelakshmi ramesh
Dear all, I am using tabulated potential do simulate nacl in water .when i was doing my energy minimization it gives me the following error.Any help please.. WARNING: For the 1498 non-zero entries for table 2 in table_NA_CL.xvg the forces deviate on average 39% from minus the numeri

Re: [gmx-users] Energy minimization error

2011-04-15 Thread Kavyashree M
Ok thanks On Fri, Apr 15, 2011 at 9:56 PM, Justin A. Lemkul wrote: > > > Kavyashree M wrote: > >> Dear gromacs users, >> >> While doing energy minimization for a protein (from pdb), with oplsaa >> force field >> and tip4p water model, there was an error and em stopped - >> >> Fatal error: >> A c

Re: [gmx-users] Energy minimization error

2011-04-15 Thread Justin A. Lemkul
Kavyashree M wrote: Dear gromacs users, While doing energy minimization for a protein (from pdb), with oplsaa force field and tip4p water model, there was an error and em stopped - Fatal error: A charge group moved too far between two domain decomposition steps This usually means that your

[gmx-users] Energy minimization error

2011-04-15 Thread Kavyashree M
Dear gromacs users, While doing energy minimization for a protein (from pdb), with oplsaa force field and tip4p water model, there was an error and em stopped - Fatal error: A charge group moved too far between two domain decomposition steps This usually means that your system is not well equilib

Re: [gmx-users] energy minimization of a charged system in vacuum

2011-02-02 Thread devicerandom
On 30/01/11 16:13, David van der Spoel wrote: On 2011-01-30 17.08, ms wrote: On 30/01/11 15:41, David van der Spoel wrote: Since I have exactly the same needs (charged system in vacuum) I jump in... In http://www.gromacs.org/Documentation/Errors it says: Note for PME users: It is possible to

Re: [gmx-users] energy minimization of a charged system in vacuum

2011-02-02 Thread devicerandom
On 30/01/11 16:13, David van der Spoel wrote: On 2011-01-30 17.08, ms wrote: Since I have exactly the same needs (charged system in vacuum) I jump in... In http://www.gromacs.org/Documentation/Errors it says: Note for PME users: It is possible to use a uniform neutralizing background charge in

Re: [gmx-users] energy minimization of a charged system in vacuum

2011-01-30 Thread David van der Spoel
On 2011-01-30 17.08, ms wrote: On 30/01/11 15:41, David van der Spoel wrote: My question is what effect does running a simulation with a charged system have? I recall reading that something related to PME calculations assumes the system is neutral, but it did not specify whether it was referrin

Re: [gmx-users] energy minimization of a charged system in vacuum

2011-01-30 Thread ms
On 30/01/11 15:41, David van der Spoel wrote: My question is what effect does running a simulation with a charged system have? I recall reading that something related to PME calculations assumes the system is neutral, but it did not specify whether it was referring to MD or EM. From the mailing

Re: [gmx-users] energy minimization of a charged system in vacuum

2011-01-30 Thread David van der Spoel
On 2011-01-30 16.02, ms wrote: On 29/01/11 05:08, Matt Chan wrote: Perfect. This is great reading. Thanks for the pointers Mark. Matt On 01/28/2011 06:50 PM, Mark Abraham wrote: On 29/01/2011 9:57 AM, Matthew Chan wrote: Hi, I'm a first time GROMACS user. I've got 2 questions, which I'll

Re: [gmx-users] energy minimization of a charged system in vacuum

2011-01-30 Thread ms
On 29/01/11 05:08, Matt Chan wrote: Perfect. This is great reading. Thanks for the pointers Mark. Matt On 01/28/2011 06:50 PM, Mark Abraham wrote: On 29/01/2011 9:57 AM, Matthew Chan wrote: Hi, I'm a first time GROMACS user. I've got 2 questions, which I'll ask in separate emails. The firs

Re: [gmx-users] energy minimization of a charged system in vacuum

2011-01-28 Thread Matt Chan
Perfect. This is great reading. Thanks for the pointers Mark. Matt On 01/28/2011 06:50 PM, Mark Abraham wrote: On 29/01/2011 9:57 AM, Matthew Chan wrote: Hi, I'm a first time GROMACS user. I've got 2 questions, which I'll ask in separate emails. The first is about running EM on a charged p

Re: [gmx-users] energy minimization of a charged system in vacuum

2011-01-28 Thread Mark Abraham
On 29/01/2011 9:57 AM, Matthew Chan wrote: Hi, I'm a first time GROMACS user. I've got 2 questions, which I'll ask in separate emails. The first is about running EM on a charged protein in vacuum. I'm presently walking through some of the tutorials and trying to simplify them for my purpose

[gmx-users] energy minimization of a charged system in vacuum

2011-01-28 Thread Matthew Chan
Hi, I'm a first time GROMACS user. I've got 2 questions, which I'll ask in separate emails. The first is about running EM on a charged protein in vacuum. I'm presently walking through some of the tutorials and trying to simplify them for my purposes. One is the energy minimization of the 1AKI lys

Re: [gmx-users] eNERGY mINIMIZATION

2011-01-17 Thread Mark Abraham
On 18/01/2011 3:12 PM, bharat gupta wrote: Sir, That's what I want to know how u found out that the P.E. is ok and then what about the PE graph ?? A value of large magnitude and negative sign is probably OK. Exactly how large varies widely with system composition, force field, etc. So nobody

Re: [gmx-users] eNERGY mINIMIZATION

2011-01-17 Thread bharat gupta
Sir, That's what I want to know how u found out that the P.E. is ok and then what about the PE graph ?? Does it has to be some what parallel to x-axis in case of every minimized structure as in my case its not parallel to x-axis , rather its getting down towards the x-axis.. ... generally for a m

Re: [gmx-users] eNERGY mINIMIZATION

2011-01-17 Thread Mark Abraham
On 18/01/2011 2:50 PM, bharat gupta wrote: Hi, I did the energy minimization of my protein (230 amino acids) after adding ions .. and after checking the potential energy of the system I was getting thee graph which is different from the one given in the lysozyme tutorial i.e. the the potentia

[gmx-users] eNERGY mINIMIZATION

2011-01-17 Thread bharat gupta
Hi, I did the energy minimization of my protein (230 amino acids) after adding ions .. and after checking the potential energy of the system I was getting thee graph which is different from the one given in the lysozyme tutorial i.e. the the potential line in graph is not parallel to the x axix..

Re: [gmx-users] Energy minimization

2011-01-17 Thread Justin A. Lemkul
bharat gupta wrote: The potential energy graph is not coming parallel to x-axis is getting down towars the x-axis .. normally it should be paralle to x-axis.. Depending on how quickly the potential energy converges, it may or may not. The most important outcome of EM is that you achieve

Re: [gmx-users] Energy minimization

2011-01-17 Thread Amit Choubey
Check your log file to see if the minimization criterion was satisfied ie the max force is less than emtol. On Mon, Jan 17, 2011 at 12:29 AM, bharat gupta wrote: > > The potential energy graph is not coming parallel to x-axis is getting down > towars the x-axis .. normally it should be paralle to

Re: [gmx-users] Energy minimization

2011-01-17 Thread bharat gupta
The potential energy graph is not coming parallel to x-axis is getting down towars the x-axis .. normally it should be paralle to x-axis.. -- Bharat Ph.D. Candidate Room No. : 7202A, 2nd Floor Biomolecular Engineering Laboratory Division of Chemical Engineering and Polymer Science Pusan National

Re: [gmx-users] Energy minimization

2011-01-17 Thread Amit Choubey
OK what seemed wrong to you? On Mon, Jan 17, 2011 at 12:24 AM, bharat gupta wrote: > I am doin the first minimization step.. here are the details of minim.mdp > file > > ; minim.mdp - used as input into grompp to generate em.tpr > ; Parameters describing what to do, when to stop and what to save

Re: [gmx-users] Energy minimization

2011-01-17 Thread bharat gupta
I am doin the first minimization step.. here are the details of minim.mdp file ; minim.mdp - used as input into grompp to generate em.tpr ; Parameters describing what to do, when to stop and what to save integrator = steep ; Algorithm (steep = steepest descent minimization) emtol = 1000.0 ; St

Re: [gmx-users] Energy minimization

2011-01-17 Thread Amit Choubey
Can you provide more details ? Which graph are you talking about ? Amit On Mon, Jan 17, 2011 at 12:05 AM, bharat gupta wrote: > Hi all, > > I am minimizing my protein (230 amino acids) with the minim.mdp file given > in the lysozyme tutorial but I am not getting a stable graph .. I checked > th

[gmx-users] Energy minimization

2011-01-17 Thread bharat gupta
Hi all, I am minimizing my protein (230 amino acids) with the minim.mdp file given in the lysozyme tutorial but I am not getting a stable graph .. I checked the manual also to look for parameters but I am not able to find such information there apart from the theory about the minimzation algos an

[gmx-users] Energy minimization - single/double precision

2011-01-13 Thread Kavyashree M
Dear Gromacs users, I wanted to know whether it is correct doing minimization in double precision and rest all - pdb2gmx, editconf, genbox, grompp and mdrun itself in single precision. I had got convergence during energy minimization using double precision but not single precision (input

Re: [gmx-users] Energy minimization query

2010-12-28 Thread TJ Mustard
Nikhil, I had a similar problem, and was getting the same error. In my case the order of my atoms was different in my pdb/gro and my top/itp file. This made my system do a weird inversion, due to incorrect placement of atoms. And my energies were terrible.  

[gmx-users] Energy minimization query

2010-12-28 Thread nikhil damle
Hi, I am trying to simulate a dimeric structure and I have already simulated its monomeric form. I use exactly same parameters which i use for monomer during dimer minimisation in solvent. System Details: # Atoms = 4488 (# Res = 582) octahedron box with d = 0.9 Solvent (spce water model) den

RE: [gmx-users] Energy Minimization

2010-10-18 Thread Dallas Warren
er 2010 7:26 AM To: Discussion list for GROMACS users Subject: [gmx-users] Energy Minimization Hi all, I am using CHARMM forcefield with tables for Protein SOL interactions for alanine dipeptide. Tables supplied: table.xvg, table_Protein_SOL.xvg, tablep.xvg (Standard tables for 6-12 in

[gmx-users] Energy Minimization

2010-10-18 Thread Sai Pooja
Hi all, I am using CHARMM forcefield with tables for Protein SOL interactions for alanine dipeptide. Tables supplied: table.xvg, table_Protein_SOL.xvg, tablep.xvg (Standard tables for 6-12 interactions used) Combination rule changed from '2' to '1' in forcefield.itp file. In the energy minimizati

Re: [gmx-users] Energy minimization by keeping heavy atom fixed

2010-08-26 Thread Mark Abraham
- Original Message - From: priyabrata panigrahi Date: Thursday, August 26, 2010 17:09 Subject: [gmx-users] Energy minimization by keeping heavy atom fixed To: gmx-users@gromacs.org > I have a protein structure determined by x-ray crystallography method where > hydrogen atom

[gmx-users] Energy minimization by keeping heavy atom fixed

2010-08-26 Thread priyabrata panigrahi
I have a protein structure determined by x-ray crystallography method where hydrogen atoms were missing. I have added hydrogen atom. Now i need to do energy minimization of this protein by keeping the heavy atom fixed at their crystallographic position. How can I manage which atoms to be fixed and

[gmx-users] Energy minimization: problems with ramachanrad score

2010-06-16 Thread chris . neale
ay from the native fold although stereochemically the model look good? -- next part -- An HTML attachment was scrubbed... URL: http://lists.gromacs.org/pipermail/gmx-users/attachments/20100617/aafe7f09/attachment.html Previous message: [gmx-users] Electric field, p

Re: [gmx-users] Energy minimization: problems with ramachanrad score

2010-06-16 Thread Lucio Ricardo Montero Valenzuela
The implicit water models were built in a time when computing time was limiting because of the computing power if this epoch. But now is feasible to do simulations with explicit water. Maybe one should argue that at the first minimization of a protein waterbox the waters are not in the natural posi

Re: [gmx-users] Energy minimization: problems with ramachanrad score

2010-06-16 Thread Thomas Evangelidis
> > Without having actually read their paper... there's could easily be a kind > of apples-and-oranges comparison going on. The local optimum resulting from > a single EM using MM forcefields optimized for non-vacuum conditions on > near-native folds in vacuo are being compared with (say) crystal s

Re: [gmx-users] Energy minimization: problems with ramachanrad score

2010-06-16 Thread Mark Abraham
- Original Message - From: Thomas Evangelidis Date: Thursday, June 17, 2010 10:48 Subject: [gmx-users] Energy minimization: problems with ramachanrad score To: gmx-users@gromacs.org > I 'm trying to relax my homology models from steric clashes, and while > searching for the

[gmx-users] Energy minimization: problems with ramachanrad score

2010-06-16 Thread Thomas Evangelidis
I 'm trying to relax my homology models from steric clashes, and while searching for the appropriate minimization scheme, I came across this old thread: http://lists.gromacs.org/pipermail/gmx-users/2007-April/027043.html The authors in the cited paper have created near-native structures as a test

[gmx-users] energy minimization + freeze atoms

2010-05-03 Thread Nilesh Dhumal
Hello, I am trying to do energy minimization with normal mode analysis for my ststem. I used freezegrps and freezedim optitions to freeze the atoms. I found that after minimization the atoms are moved. Can you tell me how to freeze atoms for energy minimzation. I need this geometry for normal

[gmx-users] energy minimization + freeze atoms

2010-05-03 Thread Nilesh Dhumal
Hello, I am trying to do energy minimization with normal mode analysis for my ststem. I used freezegrps and freezedim optitions to freeze the atoms. I found that after minimization the atoms are moved. Can you tell me how to freeze atoms for energy minimzation. I need this geometry for normal

Re: [gmx-users] energy minimization - infinite force on atom

2010-04-07 Thread Mark Abraham
On 8/04/2010 4:23 PM, Evelyne Deplazes wrote: Hi gmx-users I am somehow famliar with gromacs but only used it with coarse grained force fields (Martini) so far, not with all atom ones. I am now trying to run a simulation on a membrane protein system (protein embedded in a POPC bilayer solvated w

[gmx-users] energy minimization - infinite force on atom

2010-04-07 Thread Evelyne Deplazes
Hi gmx-users I am somehow famliar with gromacs but only used it with coarse grained force fields (Martini) so far, not with all atom ones. I am now trying to run a simulation on a membrane protein system (protein embedded in a POPC bilayer solvated with water) using the gromos force field (Gromacs

  1   2   3   >