Re: [gmx-users] C-terminus residue name in Gromos43a1

2014-02-26 Thread Francesca Vitalini
Dear Justin 2014-02-25 19:07 GMT+01:00 Justin Lemkul : > > > On 2/25/14, 12:56 PM, Francesca Vitalini wrote: > >> Dear Justin, >> >> Thanks for your answer. >> However I noticed that I had made a mistake and there is no definition of >> NAC in GROMOS53a6, it was the .rtp file of OPLS-AA I was lo

Re: [gmx-users] User-defined Potentials in Gromacs

2014-02-26 Thread Dr. Vitaly Chaban
>> I have looked up "tabulated potential" in the manual and also the suggested >> document, neither of them actually helps the situation. Why do they not? Dr. Vitaly V. Chaban -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_L

[gmx-users] please guide me through the confusing gromacs results!!!

2014-02-26 Thread delara aghaie
Dear Gromacs users, we want to simulateHSA protein using8 processors. Usually with our available system 10ns simulation on 8 processors lasts 2-3 days. This time we submitted 10 ns simulation, after almost 6 days it has finished but when we look at md.log file, only 179411 steps has been complete

Re: [gmx-users] C-terminus residue name in Gromos43a1

2014-02-26 Thread Justin Lemkul
On 2/26/14, 4:08 AM, Francesca Vitalini wrote: Dear Justin 2014-02-25 19:07 GMT+01:00 Justin Lemkul : On 2/25/14, 12:56 PM, Francesca Vitalini wrote: Dear Justin, Thanks for your answer. However I noticed that I had made a mistake and there is no definition of NAC in GROMOS53a6, it was

Re: [gmx-users] difference in potenital energy on desktop and server run

2014-02-26 Thread Justin Lemkul
On 2/26/14, 12:51 AM, gupta.rakesh082 wrote: Dear allI performed MD simulation of lipid bilayer on my desktop it is giving correct results and negative potential energy (~ -10^5) but when I ran the same simulation on a server using multi thread (-nt 8), It is giving positive potential energy (~

Re: [gmx-users] please guide me through the confusing gromacs results!!!

2014-02-26 Thread Justin Lemkul
On 2/26/14, 4:31 AM, delara aghaie wrote: Dear Gromacs users, we want to simulateHSA protein using8 processors. Usually with our available system 10ns simulation on 8 processors lasts 2-3 days. This time we submitted 10 ns simulation, after almost 6 days it has finished but when we look at md

Re: [gmx-users] please guide me through the confusing gromacs results!!!

2014-02-26 Thread Dr. Vitaly Chaban
if the system sends a termination request to MDRUN, it stops the main cycle and writes down the summary of the run. Dr. Vitaly V. Chaban On Wed, Feb 26, 2014 at 10:31 AM, delara aghaie wrote: > Dear Gromacs users, > we want to simulateHSA protein using8 processors. Usually with our available

Re: [gmx-users] Adding TPO and SEP to G53a6 Forcefield

2014-02-26 Thread lalithkumar
Dear Dr. Justin, Thank you for pointing out the details of error. I have missed out the version of GROMACS I have been using and so I would like to let you know that my version is 4.5.7. In this version ffnonbonded.itp file of the gromos53a6.ff forcefield doesnt have pairtypes fr LJ-14 interaction

Re: [gmx-users] Adding TPO and SEP to G53a6 Forcefield

2014-02-26 Thread Justin Lemkul
On 2/26/14, 7:05 AM, lalithkumar wrote: Dear Dr. Justin, Thank you for pointing out the details of error. I have missed out the version of GROMACS I have been using and so I would like to let you know that my version is 4.5.7. In this version ffnonbonded.itp file of the gromos53a6.ff forcefiel

Re: [gmx-users] difference in potenital energy on desktop and server run

2014-02-26 Thread gupta.rakesh082
I am using gromacs 4.6.5. Surprisingly system is stable at such high potential energy. I created bilayer system using packmol then solvated it with water using genbox. I did energy minimization using following command grompp -f em.mdp -c cerblayer.gro -p cer.itp -o min.tpr mdrun -c -v -deffnm min

Re: [gmx-users] difference in potenital energy on desktop and server run

2014-02-26 Thread Justin Lemkul
On 2/26/14, 10:50 AM, gupta.rakesh082 wrote: I am using gromacs 4.6.5. Surprisingly system is stable at such high potential energy. I created bilayer system using packmol then solvated it with water using genbox. I did energy minimization using following command grompp -f em.mdp -c cerblayer.

[gmx-users] how to increase vacuum thickness in a periodic system

2014-02-26 Thread decaiyu
Dear All, I ran a NVT simulation with a layer of oil+vaccuum. Because of periodic boundary condition, there are some molecules on the top of vaccuum layer. Now I want to increase the vaccum thickness so that I can add some water in it. Is there a simple way to do that? The molecules on the top of

Re: [gmx-users] how to increase vacuum thickness in a periodic system

2014-02-26 Thread Dr. Vitaly Chaban
editconf -box X Y Z Dr. Vitaly V. Chaban On Wed, Feb 26, 2014 at 7:58 PM, decaiyu wrote: > Dear All, > > I ran a NVT simulation with a layer of oil+vaccuum. > Because of periodic boundary condition, there are some molecules on the top > of vaccuum layer. > Now I want to increase the vaccum thi

[gmx-users] Martini Lo-phase seems gel-like

2014-02-26 Thread David Ackerman
Hello Everyone, I have been using the Martini CG forcefield to simulate phase separation in a bilayer containing DPPC/DUPC/Chol at a ratio of 1:1:.5. Though this does phase separate at 295 K, when I use mdp files based on http://md.chem.rug.nl/cgmartini/images/parameters/exampleMDP/martini_v2.x_ex

[gmx-users] Martini Lo-phase has gel-like properties

2014-02-26 Thread David Ackerman
Hello Everyone, I have been using the Martini CG forcefield to simulate phase separation in a bilayer containing DPPC/DUPC/Chol at a ratio of 1:1:.5. Though this does phase separate at 295 K, when I use mdp files based on http://md.chem.rug.nl/cgmartini/images/parameters/exampleMDP/martini_v2.x_ex

Re: [gmx-users] C-terminus residue name in Gromos43a1

2014-02-26 Thread Francesca Vitalini
Residue 2 is just alanine,the most standard amino acid ever. Any help? Francesca. Il 26/feb/2014 11:58 "Justin Lemkul" ha scritto: > > > On 2/26/14, 4:08 AM, Francesca Vitalini wrote: > >> Dear Justin >> >> >> 2014-02-25 19:07 GMT+01:00 Justin Lemkul : >> >> >>> >>> On 2/25/14, 12:56 PM, Franc

Re: [gmx-users] difference in potenital energy on desktop and server run

2014-02-26 Thread gupta.rakesh082
EM RUNSteepest Descents converged to machine precision in 2820 steps,but did not reach the requested Fmax < 10.Potential Energy = -2.9268122e+05Maximum force = 1.9479546e+02 on atom 123Norm of force = 6.3931437e+00grompp -f em.mdp -c cerblayer.gro -p cer.itp -o min.tprmdrun -c -v -deffnm

Re: [gmx-users] difference in potenital energy on desktop and server run

2014-02-26 Thread gupta.rakesh082
EM RUN Steepest Descents converged to machine precision in 2820 steps, but did not reach the requested Fmax < 10. Potential Energy = -2.9268122e+05 Maximum force = 1.9479546e+02 on atom 123 Norm of force = 6.3931437e+00 grompp -f em.mdp -c cerblayer.gro -p cer.itp -o min.tpr mdrun -c -