Re: [gmx-users] trjconv and ions

2010-03-25 Thread Joe Joe
t; > Hope it helps, > > Tsjerk > > On Mar 25, 2010 7:22 PM, "Joe Joe" wrote: > > I am seeing weird behavior when I trjconv my trajectory containing NaCl, > tip3p, and protein. I use the OPLSAA force field and octahedral PBCs. > > When I run > > trjconv -f

[gmx-users] trjconv and ions

2010-03-25 Thread Joe Joe
I am seeing weird behavior when I trjconv my trajectory containing NaCl, tip3p, and protein. I use the OPLSAA force field and octahedral PBCs. When I run trjconv -f traj.trr -s topol.tpr -pbc nojump -o trajout.xtc The ions in the output form clusters of three similarly charged ions in the shaped

[gmx-users] problem with using pdb2gmx on a single amino acid

2010-03-12 Thread Joe Joe
I think there is a tdb problem when using pdb2gmx on a single amino acid. The tdb databases replaces the CA with different charges depending on whether its a c or n termini. When its a single amino acid you get the c terminal charge and thus the overall charge is incorrect. Am I missing something?

Re: [gmx-users] Regarding defining termini in new rtp entry

2010-03-12 Thread Joe Joe
ory. You can also copy it to your > working directory and change it there, as it will then take precedence > over the other one. Do mind to increment the number at the top of the > file for every amino acid you add. > > Hope it helps, > > Tsjerk > > On Fri, Mar 12, 2010 a

[gmx-users] Regarding defining termini in new rtp entry

2010-03-12 Thread Joe Joe
How do I tell pdb2gmx what the termini are in my molecule. I added a non standard amino acid to the ffoplsaa.rtp database but pdb2gmx says "No N- or C-terminus found: this chain appears to contain no protein". thanks, Ilya -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.o

Re: [gmx-users] topolbuild 1.3, revision to topolbuild that incorporates OPLS-AA support

2010-03-10 Thread Joe Joe
Bruce, This is great. I was able to easily make a molecule that runs with OPLS-AA. Would it be possible to add a feature which creates an *.rtp entry? It would be really useful for making non standard amino acids and incorporating it with pdb2gmx. Thanks, Ilya On Wed, Mar 10, 2010 at 11:07 AM,

Re: [gmx-users] g_dist and vsites

2010-03-10 Thread Joe Joe
-09 20.01, Joe Joe wrote: > >> I also just tried on a different ALA and got the same problem. >> > > can you please submit a bugzilla and upload files to reproduce the problem? > > >> Ilya >> >> >> On Tue, Mar 9, 2010 at 10:58 AM, Joe Joe > <mail

Re: [gmx-users] g_dist and vsites

2010-03-09 Thread Joe Joe
done, bug #402 On Tue, Mar 9, 2010 at 11:44 AM, David van der Spoel wrote: > On 2010-03-09 20.01, Joe Joe wrote: > >> I also just tried on a different ALA and got the same problem. >> > > can you please submit a bugzilla and upload files to reproduce the problem? > >

Re: [gmx-users] g_dist and vsites

2010-03-09 Thread Joe Joe
I also just tried on a different ALA and got the same problem. Ilya On Tue, Mar 9, 2010 at 10:58 AM, Joe Joe wrote: > I am looking the the gmxdump output and everything seems consistent. The > only difference is that the numbering in the index file starts at 1 whereas > in the gm

Re: [gmx-users] g_dist and vsites

2010-03-09 Thread Joe Joe
I am looking the the gmxdump output and everything seems consistent. The only difference is that the numbering in the index file starts at 1 whereas in the gmxdump the arrays are indexed starting at 0. The coordinates look just fine. Thanks, Ilya On Tue, Mar 9, 2010 at 10:30 AM, Joe Joe

Re: [gmx-users] g_dist and vsites

2010-03-09 Thread Joe Joe
On Tue, Mar 9, 2010 at 10:22 AM, David van der Spoel wrote: > On 2010-03-09 19.16, Joe Joe wrote: > >> yep. >> > Have you gmxdump/ed the xtc to check the coordinates are right? > >> >> On Tue, Mar 9, 2010 at 10:15 AM, David van der Spoel >> mailto:sp...@x

[gmx-users] g_dist and vsites

2010-03-09 Thread Joe Joe
Hi I am trying to post process and xtc trajectory using g_dist. I am trying to calculate the CA-CB distance of an Alanine residue but I get NAN in all the distance columns. It works for the other residues I've tried (i.e. SER, VAL). I am using vsites in my simulation and I think it may have some th

Re: [gmx-users] g_dist and vsites

2010-03-09 Thread Joe Joe
yep. On Tue, Mar 9, 2010 at 10:15 AM, David van der Spoel wrote: > On 2010-03-09 19.09, Joe Joe wrote: > >> Hi I am trying to post process and xtc trajectory using g_dist. I am >> trying to calculate the CA-CB distance of an Alanine residue but I get >> NAN in all the di

[gmx-users] g_dist and vsites

2010-03-09 Thread Joe Joe
Hi I am trying to post process and xtc trajectory using g_dist. I am trying to calculate the CA-CB distance of an Alanine residue but I get NAN in all the distance columns. It works for the other residues I've tried (i.e. SER, VAL). I am using vsites in my simulation and I think it may have some th

Re: [gmx-users] Replica Exchange Error

2009-04-12 Thread Joe Joe
Will do thanks. Ilya On Sun, Apr 12, 2009 at 7:50 PM, Mark Abraham wrote: > Joe Joe wrote: > >> It works just fine without RE. >> > > OK well that sounds like after a replica exchange attempt, some data > structure isn't being (re-)constructed correctly. See if yo

Re: [gmx-users] Replica Exchange Error

2009-04-12 Thread Joe Joe
It works just fine without RE. On Sat, Apr 11, 2009 at 8:49 PM, Mark Abraham wrote: > Joe Joe wrote: > >> Hi Justin, >> >> Thanks for helping me with this. >> On Sat, Apr 11, 2009 at 4:46 AM, Justin A. Lemkul > jalem...@vt.edu>> wrote: >> >

Re: [gmx-users] Replica Exchange Error

2009-04-11 Thread Joe Joe
On Sat, Apr 11, 2009 at 8:13 PM, Justin A. Lemkul wrote: > > > Joe Joe wrote: > >> Hi Justin, >> >> Thanks for helping me with this. >> On Sat, Apr 11, 2009 at 4:46 AM, Justin A. Lemkul > jalem...@vt.edu>> wrote: >> >> >> >>

Re: [gmx-users] Replica Exchange Error

2009-04-11 Thread Joe Joe
Hi Justin, Thanks for helping me with this. On Sat, Apr 11, 2009 at 4:46 AM, Justin A. Lemkul wrote: > > > Joe Joe wrote: > >> How do I determine what those values should be? How does it affect the >> simulation? >> >> > The header of your log file shoul

Re: [gmx-users] Replica Exchange Error

2009-04-10 Thread Joe Joe
How do I determine what those values should be? How does it affect the simulation? thanks, Ilya On Fri, Apr 10, 2009 at 8:00 PM, Justin A. Lemkul wrote: > > > Joe Joe wrote: > >> Hi All, >> >> I get the following error when running replica exchange. I have 32

[gmx-users] Replica Exchange Error

2009-04-10 Thread Joe Joe
Hi All, I get the following error when running replica exchange. I have 32 replica's with each replica running on 4 processors (128 total). How do I decrease the domain decomposition cell size, -dss? What should I set it to? Thanks, Ilya DD cell 2 0 0: Neighboring cells do not have atoms: 22217

Re: [gmx-users] Pressure Coupling Problem

2009-04-10 Thread Joe Joe
I finally figured it out. I went through every parameter step by step and it turns our I had epsilon-r set to 80. Not sure why I had that. Wish gromacs would have given me a warning (hint hint). That explains why my P.E. was 10^-5 instead of 10^-6. Thanks everyone for trying!!!. Sometime the most o

Re: [gmx-users] Pressure Coupling Problem

2009-04-09 Thread Joe Joe
Nope not an A/nm problem. As a simple test I take spc.gro from share/top. I reconfigure the box (i.e. editconf -f spc.gro -d 1.0 -c -bt cubic -o water_center. I then solvate with genbox, minimize and run using the mdp file I provided earlier. No matter how much minimization I do the volume of the

Re: [gmx-users] Pressure Coupling Problem

2009-04-09 Thread Joe Joe
On Thu, Apr 9, 2009 at 6:36 AM, Justin A. Lemkul wrote: > > > Joe Joe wrote: > >> So I got my small water box (800 waters) to behave stably with pressure >> coupling after more minimization but I still can't get my large system to >> work with pressure cou

Re: [gmx-users] Release of R.E.D. Server

2009-04-09 Thread Joe Joe
What about use for commercial purposes? On Thu, Apr 9, 2009 at 1:31 AM, FyD wrote: > Dear All, > > I am pleased to announce the release of R.E.D. Server available @ > http://q4md-forcefieldtools.org/REDS/. > > R.E.D. Server provides the software and hardware (i. e. a cluster of > computers) requ

Re: [gmx-users] Pressure Coupling Problem

2009-04-09 Thread Joe Joe
t) and for all terms an amplitude (real) ; and a phase angle (real) E-x = E-xt = E-y = E-yt = E-z = E-zt = On Wed, Apr 8, 2009 at 1:00 PM, Joe Joe wrote: > > > On

Re: [gmx-users] Pressure Coupling Problem

2009-04-08 Thread Joe Joe
On Wed, Apr 8, 2009 at 11:31 AM, Roland Schulz wrote: > > > On Wed, Apr 8, 2009 at 7:53 AM, Joe Joe wrote: > >> HI Chris, >> >> On Tue, Apr 7, 2009 at 9:31 PM, wrote: >> >>> Hi Ilya, >>> >>> First thing that comes to mind is that

Re: [gmx-users] Pressure Coupling Problem

2009-04-08 Thread Joe Joe
I tried .1, and 10 ps tau_p values. I guess I can try smaller values. On Wed, Apr 8, 2009 at 10:44 AM, Justin A. Lemkul wrote: > > > Joe Joe wrote: > >> Hi Chris, >> >> When I create the topology for the 4fs timestep I use pdb2gmx -vsite h. I >> set up the

Re: [gmx-users] Pressure Coupling Problem

2009-04-08 Thread Joe Joe
Yeah I only gave a partial. Tried to remove the QM params. I do use constraints = all-bonds. On Wed, Apr 8, 2009 at 9:18 AM, wrote: > You say "I run all he sims with constraints=all-bonds", but I don't see > that in the mdp options that you provided. I even put your text in a file > and grepped

Re: [gmx-users] Pressure Coupling Problem

2009-04-08 Thread Joe Joe
Hi Chris, When I create the topology for the 4fs timestep I use pdb2gmx -vsite h. I set up the correct constraints. I've tested it and it conserves energy in NVE. I run all he sims with constraints=all-bonds. I am now running a single water box (800 water molecules) with 1s time steps and the volum

Re: [gmx-users] Pressure Coupling Problem

2009-04-08 Thread Joe Joe
HI Chris, On Tue, Apr 7, 2009 at 9:31 PM, wrote: > Hi Ilya, > > First thing that comes to mind is that it is strange to couple a coulombic > switching function with PME. While this could possibly be done correctly, I > doubt that it is in fact done in the way that you expect (i.e. correctly) in

[gmx-users] Pressure Coupling Problem

2009-04-07 Thread Joe Joe
Hi I am having some pressure coupling issues. I have a fairly large protein/water system 400K+ atoms. It minimizes just fine (F < 1000). If I run NVE it conserves energy with appropriate parameter settings. If I run NVT it is stable. When I turn on Pcoupl (i.e. Berendsen or Parinello Rahman), the s

Re: [gmx-users] Energy Conservation with 4fs timestep

2009-04-01 Thread Joe Joe
Conservation with 4fs timestep > > > > Joe Joe wrote: > > > Hi, > > > > > > I get good conservation when running NVE in gromacs with 4 fs when I > use > > > PME-switch for electrostatics but not so good when I use switch. Any > > > thought

[gmx-users] Energy Conservation with 4fs timestep

2009-03-31 Thread Joe Joe
Hi, I get good conservation when running NVE in gromacs with 4 fs when I use PME-switch for electrostatics but not so good when I use switch. Any thoughts why that would be? Params shown below. Thanks, Ilya ; ; File 'mdout.mdp' was generated ; By user: ichorny (502) ; On host:

Re: [gmx-users] Elastic Network Model

2009-03-31 Thread Joe Joe
On Mon, Mar 30, 2009 at 3:23 PM, Joshua Adelman wrote: > I am interested in potentially using gromacs to do some calculations on an > C-alpha based elastic network model of the protein that I'm working on. I'm > curious if other users have done similar types of calculations, and if so, > what prot