Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread Justin A. Lemkul
Justin A. Lemkul wrote: ERIKSSON, EMMA wrote: I was using those "strange" values of tau_p and compressibility to keep the z box length fixed in order to avoid problems associated with scaling the positions of the molecules in the box when we constrain the distance between DPPC and the smal

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread Justin A. Lemkul
Gromacs Users' List Ämne: Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box ERIKSSON, EMMA wrote: Pcoupl = Parrinello-Rahman pcoupltype = semiisotropic tau_p= 1.0 1.0e-14 compressibility =

RE: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread ERIKSSON, EMMA
jalem...@vt.edu] Skickat: den 29 april 2010 17:20 Till: Gromacs Users' List Ämne: Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box ERIKSSON, EMMA wrote: > Pcoupl = Parrinello-Rahman > pcoupltype

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread Justin A. Lemkul
ERIKSSON, EMMA wrote: Pcoupl = Parrinello-Rahman pcoupltype = semiisotropic tau_p= 1.0 1.0e-14 compressibility = 4.5e-5 4.5e-15 I would bet almost anything that this is the cause of your problem. How did you come up with these bi

RE: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread ERIKSSON, EMMA
alway, Ireland Från: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] för Justin A. Lemkul [jalem...@vt.edu] Skickat: den 29 april 2010 15:44 Till: Gromacs Users' List Ämne: Re: [gmx-users] Constrained simulations crash when bilayer moves in

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread XAvier Periole
...@gromacs.org [gmx-users- boun...@gromacs.org] för XAvier Periole [x.peri...@rug.nl] *Skickat:* den 27 april 2010 16:32 *Till:* Discussion list for GROMACS users *Ämne:* Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box you should remove the water and li

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread Justin A. Lemkul
] för Justin A. Lemkul [jalem...@vt.edu] Skickat: den 29 april 2010 15:28 Till: Gromacs Users' List Ämne: Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box ERIKSSON, EMMA wrote: The last section of my mdp file:

RE: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread ERIKSSON, EMMA
- Galway Galway, Ireland Från: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] för Justin A. Lemkul [jalem...@vt.edu] Skickat: den 29 april 2010 15:28 Till: Gromacs Users' List Ämne: Re: [gmx-users] Constrained simulations crash when bilayer

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread Justin A. Lemkul
tional University of Ireland - Galway Galway, Ireland Från: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] för Justin A. Lemkul [jalem...@vt.edu] Skickat: den 29 april 2010 14:33 Till: Gromacs Users' List Ämne: Re: [gmx-users] Constrain

RE: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread ERIKSSON, EMMA
eland Från: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] för XAvier Periole [x.peri...@rug.nl] Skickat: den 29 april 2010 14:28 Till: jalem...@vt.edu; Discussion list for GROMACS users Ämne: Re: [gmx-users] Constrained simulations crash when bilayer moves in the z directi

RE: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread ERIKSSON, EMMA
y National University of Ireland - Galway Galway, Ireland Från: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] för Justin A. Lemkul [jalem...@vt.edu] Skickat: den 29 april 2010 14:33 Till: Gromacs Users' List Ämne: Re: [gmx-users] Constrai

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread Justin A. Lemkul
. -Justin Emma Från: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] för Justin A. Lemkul [jalem...@vt.edu] Skickat: den 29 april 2010 14:09 Till: Discussion list for GROMACS users Ämne: Re: [gmx-users] Constrained simulations crash when

RE: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread ERIKSSON, EMMA
pbc. Emma Från: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] för Justin A. Lemkul [jalem...@vt.edu] Skickat: den 29 april 2010 14:09 Till: Discussion list for GROMACS users Ämne: Re: [gmx-users] Constrained simulations crash when bilayer

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread XAvier Periole
romacs.org] för XAvier Periole [x.peri...@rug.nl] *Skickat:* den 27 april 2010 16:32 *Till:* Discussion list for GROMACS users *Ämne:* Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box you should remove the water and lipid bilayer COM separately. I am

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread Justin A. Lemkul
ce. Emma *Från:* gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] för XAvier Periole [x.peri...@rug.nl] *Skickat:* den 27 april 2010 16:32 *Till:* Discussion list for GROMACS users *Ämne:* Re: [gmx-users] Constrained simulations crash when bilayer moves

RE: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread ERIKSSON, EMMA
] för Justin A. Lemkul [jalem...@vt.edu] Skickat: den 27 april 2010 15:37 Till: Discussion list for GROMACS users Ämne: Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box ERIKSSON, EMMA wrote: > Hi again, > Thanks Xavier for your reply. I

RE: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-29 Thread ERIKSSON, EMMA
XAvier Periole [x.peri...@rug.nl] Skickat: den 27 april 2010 16:32 Till: Discussion list for GROMACS users Ämne: Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box you should remove the water and lipid bilayer COM separately. I am not sure what you should

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-27 Thread XAvier Periole
idn't understand you explanation. Emma Från: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] för XAvier Periole [x.peri...@rug.nl] Skickat: den 27 april 2010 12:41 Till: Discussion list for GROMACS users Ämne: Re: [gmx-users] Constrained simulations crash when bilayer moves in t

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-27 Thread Justin A. Lemkul
n't understand you explanation. Emma *Från:* gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] för XAvier Periole [x.peri...@rug.nl] *Skickat:* den 27 april 2010 12:41 *Till:* Discussion list for G

RE: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-27 Thread ERIKSSON, EMMA
Sorry that I didn't understand you explanation. Emma Från: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] för XAvier Periole [x.peri...@rug.nl] Skickat: den 27 april 2010 12:41 Till: Discussion list for GROMACS users Ämne: Re: [gmx-users] Cons

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-27 Thread XAvier Periole
The crash seems to be expected as by crossing the pbc the distance will change significantly and in way the algorithm can not handle. Note that the overall translational motion of your system should always be removed. The removal of the COM motion of your bilayer should be sufficient to

[gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-27 Thread ERIKSSON, EMMA
Hi all! I have a problem regarding lipid bilayer simulations in Gromacs 4. During some of my simulations the whole system is moving in the z direction in the box, meaning that after some time the lipids are going out in the bottom of the box and coming in in the top of the box, since I'm using

[gmx-users] CONSTRAINED SIMULATIONS

2010-01-24 Thread chris . neale
Dear Giuseppe: This is a pretty good example of how it can be useful to completely describe what you are doing at first posting. What you mention could easily be a problem for a couple of reasons: 1. The manual indicates that you may need more than one shake iteration in this case (or at

Re: [gmx-users] CONSTRAINED SIMULATIONS

2010-01-24 Thread Giuseppe Pellicane
Dear Chris: I don't think it is something related to pbc because my solutes are well inside the simulation box (as I told you before I visualize them via a movie while the simulation proceeds). Nor it is something related to the barostat because I saw that warning and in order to try to see it wo

[gmx-users] CONSTRAINED SIMULATIONS

2010-01-24 Thread chris . neale
Dear Giuseppe: OK, those are the forces. They seem pretty huge and massively fluctuating, although I use umbrella sampling myself so this might be quite normal for constraint sampling (something for you to look into). The only thing that comes to mind is that your running up against a pbc

Re: [gmx-users] CONSTRAINED SIMULATIONS

2010-01-24 Thread Giuseppe Pellicane
Got it! Ok, first lines: # CONSTRAINT3.0 # Component selection: 1 1 1 # nSkip 1 # Ref. Group 'protein_b' # Nr. of pull groups 1 # Group 1 'protein_a' Pos. 3.099513 2.982803 3.377482 # 0.00-787.501479 0.002000-271.159600 0.004000-1147.35

[gmx-users] CONSTRAINED SIMULATIONS

2010-01-24 Thread chris . neale
Dear Giuseppe: *** What I am looking for is raw -pd pull.pdo data from the first run. *** It looks like you did not define a name for -pd: mpiexec -n $NSLOTS mdrun_mpi -np 1 -v -s start_res.tpr -o output_res -c output_conf -pi pull.ppa -pn mdgrp1.ndx -po pullout -pdo pullout1 * so based on

Re: [gmx-users] CONSTRAINED SIMULATIONS

2010-01-24 Thread Giuseppe Pellicane
Chris: you're right man, let me be more concise. Here there is the execution line I used: mpiexec -n $NSLOTS mdrun_mpi -np 1 -v -s start_res.tpr -o output_res -c output_conf -pi pull.ppa -pn mdgrp1.ndx -po pullout -pdo pullout1 * As for the pull.ppa file, here you go: ; GENERAL verbose

[gmx-users] CONSTRAINED SIMULATIONS

2010-01-24 Thread chris . neale
Dear Giuseppe: I don't think your method of showing the change is very good since it introduces unnecessary variables (e.g. did you use the correct files for the second run). What I am looking for is raw -pd pull.pdo data from the first run. I found no such file in the body of the email

Re: [gmx-users] CONSTRAINED SIMULATIONS

2010-01-24 Thread Giuseppe Pellicane
Hi Chris, the pdo file was attached at the bottom of my previous e-mail. As for the output, you may read the information about the COMs in the md.log file. Here there is an extended portion of it at the beginning of the simulation: **

[gmx-users] CONSTRAINED SIMULATIONS

2010-01-24 Thread chris . neale
Please provide actual gromacs output and tell us where it is from. I know it's sad, but not all of us can recall what the gromacs 3.3.3 .pdo file format looked like. So please include a sufficiently large portion of the file to help us recall. If, on the other hand, these values that you pl

[gmx-users] CONSTRAINED SIMULATIONS

2010-01-24 Thread Giuseppe Pellicane
Hello, I'm trying to do constrained simulations of a system consisting of two solvated proteins separated at a given distance. Inside the box there are also around 14.000 water molecules. At the bottom there is the .pdo file (I'm using gromacs 3.3.3). However, it seems as it does not work. At the