[gmx-users] Re: Umbrella sampling simulations using make_edi with -restrain and -harmonic

2013-09-12 Thread Poker Chen
I think you should add -tau 0 when using flooding as harmonic restraint (since you do not want the flooding potential to change). Okay, adding -tau 0 does turn it properly into a harmonic potential, according to edsam.xvg. However, the average position of the resulting (gaussian) is near but not

[gmx-users] RE: Umbrella sampling question

2013-07-29 Thread Christopher Neale
Why not do two umbrella sampling simulations: one with initial conformations from your faster pulling and one with initial conformations from your slower pulling. Then you can run them both as regular US simulations until (a) neither US PMF is drifting systematically with increasing simulation t

[gmx-users] Re: Umbrella Sampling _ pulled ion

2013-07-25 Thread Thomas Schlesier
Think your .mdp-file looks reasonable. If you are totaly unsure you could determine the PMF of two water molecules. As a reference one can use the radial distribution function g(r) and calculate the PMF as V_pmf(r) = -kT LN(g(r)) One side note: Since the force constant you are using is relativ

Re: [gmx-users] Re: Umbrella Sampling settings

2013-07-12 Thread Justin Lemkul
Friday, July 12, 2013 9:04 PM Subject: [gmx-users] Re: Umbrella Sampling settings In GROMACS groups are called via the *.ndx file (default: index.ndx) Be aware that 'pull_dim' determines in which diretions (x,y,z) the umbrella potential acts. So use N N Y , if you want that the ion can

Re: [gmx-users] Re: Umbrella Sampling settings

2013-07-12 Thread Shima Arasteh
s" (Y). So what should I do?!   Sincerely, Shima - Original Message - From: Thomas Schlesier To: gmx-users@gromacs.org Cc: Sent: Friday, July 12, 2013 9:04 PM Subject: [gmx-users] Re: Umbrella Sampling settings In GROMACS groups are called via the *.ndx file (default: inde

[gmx-users] Re: Umbrella Sampling settings

2013-07-12 Thread Thomas Schlesier
In GROMACS groups are called via the *.ndx file (default: index.ndx) Be aware that 'pull_dim' determines in which diretions (x,y,z) the umbrella potential acts. So use N N Y , if you want that the ion can move freely (onsidering the pull) in the xy-plane and Y Y Y if you want to also restrit th

[gmx-users] Re: Umbrella Sampling settings

2013-07-12 Thread Thomas Schlesier
But you need each of these lines for both cases (SMD and US). Probably one could skip two lines and use the default values, but it's better to set them manually. See below for comments (comments are under the related entry): Thanks for your reply. But when I don't understand why these extra

[gmx-users] Re: Umbrella sampling- force vs time plots

2013-07-10 Thread Thomas Schlesier
Generally: Using a higher force constant and / or pulling velocity drives the system faster out of equilibrium, which results in higher rupture forces. Varying the force constant has two effects. The softer the potential is, the larger are the fluctuations in the coordinates but the lower are th

[gmx-users] Re: umbrella sampling for two polymer interaction

2013-06-03 Thread Thomas Schlesier
Sorry, my writing was not really excat. If you use a reference group, the force / potential acting on the pulled group is always relative to the reference group. If you use 'pull_geometry = distance' the origin potential is always in the distance 'pull_init1' along the vector from the reference

Re: [gmx-users] Re: umbrella sampling for two polymer interaction

2013-06-01 Thread gromacs query
Dear Thomas, >> In AMBER I remember I did once methane-methane interaction, just distance based umbrella sampling. But there I did not >> provide any direction. So should it not be N N N in Gromacs if I want to allow them freely move in xyz. Sorry am wrong that should be Y Y Y in Gromcacs! for su

Re: [gmx-users] Re: umbrella sampling for two polymer interaction

2013-05-31 Thread gromacs query
Dear Thomas, Thanks a lot again for great reply. Please clarify this also, >> If one uses only 'Y N N' B would only move along the x-axis due to the pull, but could move freely in the yz-plane >>You never want to use the pull-code and 'pull_dim = N N N' >>This would mean that there is no force ac

[gmx-users] Re: umbrella sampling for two polymer interaction

2013-05-31 Thread Thomas Schlesier
For comments to your questions see below. More general: (somewhat longer than i wanted. Hope you find some answers here) Imagine two interacting particles A and B which are alinged to the x-axis. We take A as the reference group, B as pulled group and put the origin of the umbrella potential

Re: [gmx-users] Re: umbrella sampling for two polymer interaction

2013-05-31 Thread gromacs query
Dear Thomas, Thanks a lot for your time and nice explanation. I was not able to get specially the pull_start flag but now its quite clear. I feel sorry, that should be pull_dim = N N N in my case. Also I will be much thankful if you please can help me to make understand following: 1) >> If you d

[gmx-users] Re: umbrella sampling for two polymer interaction

2013-05-31 Thread Thomas Schlesier
Look also into the manual. But the tutorial is a nice place to start. For further comments see below: Dear Lloyd, I have read that but my system is different regards, On Thu, May 30, 2013 at 8:28 PM, lloyd riggs wrote: >Dear Jiom, > >Look at justines tutorial, there's example pull .mdp. >

[gmx-users] Re: umbrella sampling

2013-05-13 Thread Thomas Schlesier
Hi, it seems that you've only coupled your protein to the thermostat, but not the solvent, hence the error message. Generally one would couple both domains of the protein to one thermostat and the solvent (inluding ions) to another thermostat. Side note: If you want to use the WHAM method for

[gmx-users] Re: Umbrella sampling with large pulling distance (larger than half of the box size) (Thomas Schlesier)

2013-03-20 Thread DeChang Li
>-- > >Message: 1 >Date: Wed, 20 Mar 2013 12:55:45 +0100 >From: Thomas Schlesier >Subject: [gmx-users] Umbrella sampling with large pulling distance >(larger than half of the box size) >To: >Message-ID: <5149a3c1.2000...@

[gmx-users] Re: Umbrella sampling and free energy estimation

2013-01-15 Thread James Starlight
Restating of my question about ussage of essential eigenvectors as the pulling coordinates I'd like to introduce some example. The case study is the open to close transition in some protein. For instance by means of essential dynamics sampling technique (in targeted-contraction mode) I've simulated

Re: [gmx-users] Re: Umbrella sampling question

2012-11-16 Thread Gmx QA
Thanks Erik! /PK 2012/11/16 Erik Marklund > Hi, > > Blindly defining the center of mass for a group of atoms is not possible > in a periodic system such as a typical simulation box. You need some clue > as to which periodic copy of every atom that is to be chosen. By providing > pull_pbcatom0 y

Re: [gmx-users] Re: Umbrella sampling question

2012-11-16 Thread Erik Marklund
Hi, Blindly defining the center of mass for a group of atoms is not possible in a periodic system such as a typical simulation box. You need some clue as to which periodic copy of every atom that is to be chosen. By providing pull_pbcatom0 you tell gromacs to, for every atom in grp0, use the pe

[gmx-users] Re: Umbrella sampling question

2012-11-15 Thread Gmx QA
Hi Chris Seems my confusion was that I assumed that the distances in the profile.xvg-file should correspond to something I could measure with g_dist. Turns out it does not. Thank you for helping me sorting out this, I got it now :-) About pull_pbcatom0 though. My box is >> 2*1.08 nm in all direct

Re: [gmx-users] Re: Umbrella sampling question

2012-11-15 Thread Gmx QA
Hi Erik, Thank you for your answer. I see your point now. Went and had a look in gmx_wham.c to see how things are calculated, and that makes sense. I was looking for an easy way of relating different parts of the resulting PMF to my original starting frames, as a means to understand exactly what

Re: [gmx-users] Re: Umbrella sampling question

2012-11-14 Thread Erik Marklund
Hi, See below. 14 nov 2012 kl. 15.06 skrev Gmx QA: > Hi Chris, > > and thank you for your reply. I should have included my g_dist command in > my first mail. Here goes: > > I first run trjconv to extract the individual frames from my pulling > trajectory > $ trjconv -f pull.xtc -s pull.tpr -o

[gmx-users] Re: Umbrella sampling question

2012-11-14 Thread Christopher Neale
What you reported is not what you did. It appears that grompp, gtraj, and g_dist report the same value. Please also report the value that you get from your pullx.xvg file that you get from mdrun, which I suspect will also be the same. The difference that you report is actually between the first

[gmx-users] Re: Umbrella sampling question

2012-11-14 Thread Gmx QA
Hi Chris, and thank you for your reply. I should have included my g_dist command in my first mail. Here goes: I first run trjconv to extract the individual frames from my pulling trajectory $ trjconv -f pull.xtc -s pull.tpr -o conf.gro -sep -pbc mol -ur compact Then, g_dist like so: $ g_dist -s

RE: [gmx-users] Re: Umbrella Sampling Pull code Problem - Follow-Up

2012-09-14 Thread Steinbrecher, Thomas (IPC)
Dear Gromacs users, a while ago I have sent out a message regarding different distance measurements in the pull code vs. those made by g_dist. Thanks to everybody who responded and observed the same problem. I am still somewhat puzzled. Several users suggested that this is connected to the 'sig

[gmx-users] Re: umbrella sampling (PMF) position discrepancy

2012-09-10 Thread Sheeba
Hi, In the output pmf that you get from wham, do you have negative distances? You mentioned your reaction coordinate ranges from -2 to 2 and your pullx files have negative distances but from what I have observed, wham takes only the magnitude of the distances and I presume your wham output will ran

[gmx-users] Re: Umbrella Sampling Pull code Problem

2012-08-31 Thread Sheeba
As Bogdan suggested this could be an issue with the sign of the distance. When you run grompp you get >Pull group natoms pbc atom distance at start reference at t=0 > 0 994 497 > 1 2 46347 1.224 1.950 which suggests your startin

Re: [gmx-users] Re: umbrella sampling with distances larger than half box size

2012-06-25 Thread lloyd riggs
the positive number should be in the calculations used inside WHAM. Grüess Stephan Original-Nachricht > Datum: Mon, 25 Jun 2012 02:32:40 -0700 (PDT) > Von: anaome > An: gmx-users@gromacs.org > Betreff: [gmx-users] Re: umbrella sampling with distances larg

[gmx-users] Re: umbrella sampling with distances larger than half box size

2012-06-25 Thread anaome
Thanks Stephan for your suggestions, A finally found a solution for distance larger than half the box size: 1) Regenerate the tpr files with pbc=no 2) Perform WHAM analysis with the umbrella forces (-if) -- View this message in context: http://gromacs.5086.n6.nabble.com/umbrella-sampling-with-d

[gmx-users] Re: Umbrella Sampling - ligand

2012-03-01 Thread Justin A. Lemkul
Steven Neumann wrote: Dear Gmx Users, Dear Justin, I run Umbrella Sampling for my ligand pulling using Justin's Tutorial. Please, see attached my histogram and PMF profile. I know it requires more sampling around distance from 0.8 to 2.5 nm. 1) Will more sampling affect the final deltaG wh

Re: [gmx-users] Re: Umbrella Sampling - Justin tutorial

2011-11-17 Thread Steven Neumann
Thank you Justin, I will read it for sure and come back to my simulations! Steven On Thu, Nov 17, 2011 at 1:25 PM, Justin A. Lemkul wrote: > > > Steven Neumann wrote: > >> Hi Justin, >> I am sorry for so many questions but I do not understand something. >> First we run the simulation of pullin

[gmx-users] Re: Umbrella Sampling - Justin tutorial

2011-11-17 Thread Justin A. Lemkul
Steven Neumann wrote: Hi Justin, I am sorry for so many questions but I do not understand something. First we run the simulation of pulling Chain A from ChainB with constant force (pull_k1=1000) and constant velocity of pulling (pull_rate1=0.01) We extract windows as we discussed and then r

[gmx-users] Re: Umbrella Sampling - Justin tutorial

2011-11-17 Thread Steven Neumann
Hi Justin, I am sorry for so many questions but I do not understand something. First we run the simulation of pulling Chain A from ChainB with constant force (pull_k1=1000) and constant velocity of pulling (pull_rate1=0.01) We extract windows as we discussed and then run simulations with those con

Re: [gmx-users] Re: umbrella sampling with "pull=constraint"

2011-11-07 Thread Justin A. Lemkul
From: "Justin A. Lemkul" mailto:jalem...@vt.edu>> Subject: Re: [gmx-users] Re: umbrella sampling with "pull=constraint" To: Discussion list for GROMACS users mailto:gmx-users@gromacs.org>> Message-ID: <4eb7f727.3000...@vt.edu <mailto:4eb7f727.3000...

[gmx-users] Re: umbrella sampling with "pull=constraint"

2011-11-07 Thread Vijayaraj
sampling windows? this observation is from last 10ns of 20ns simulation. > > Message: 3 > Date: Mon, 07 Nov 2011 10:20:07 -0500 > From: "Justin A. Lemkul" > Subject: Re: [gmx-users] Re: umbrella sampling with "pull=constraint" > To: Discussion list for GROMA

Re: [gmx-users] Re: umbrella sampling with "pull=constraint"

2011-11-07 Thread Justin A. Lemkul
Vijayaraj wrote: and also I tried to plot the PMF curve from this data, I got few warning messages like "no data point in bin 20, You may not get a reasonable profile. Check your histograms!". the histogram shows poor sampling from 0.6 nm to 0.7 nm, and the windows in this position completely

[gmx-users] Re: umbrella sampling with "pull=constraint"

2011-11-07 Thread Vijayaraj
and also I tried to plot the PMF curve from this data, I got few warning messages like "no data point in bin 20, You may not get a reasonable profile. Check your histograms!". the histogram shows poor sampling from 0.6 nm to 0.7 nm, and the windows in this position completely overlaps at 0.5-0.6 nm

[gmx-users] Re: umbrella sampling with "pull=constraint"

2011-11-07 Thread Vijayaraj
Hello, Thanks for the suggestions. I did umbrella simulation for 20ns and the pull force is as below. 19960.-6.78192 19962.33.3579 19964.-3.1808 19966.-15.0304 19968.-55.4436 19970.-38.9422 19972.-9.41927 19974.-5.95981 19976.

[gmx-users] Re: umbrella sampling

2011-10-28 Thread Vijayaraj
Hi Justin, Here is my complete code for umbrella sampling title = Umbrella pulling simulation define = -DPOSRES_2 integrator = md dt = 0.002 tinit = 0 nsteps = 250 ; 5 ns nstcomm = 10 nstxout = 5 ; every 100 ps nstvout = 5 nstfout

Re: [gmx-users] Re: Umbrella sampling with temperature and pressure coupling method problem (Justin A. Lemkul)

2010-10-19 Thread Justin A. Lemkul
DeChang Li wrote: Message: 6 Date: Tue, 19 Oct 2010 09:30:47 -0400 From: "Justin A. Lemkul" mailto:jalem...@vt.edu>> Subject: Re: [gmx-users] Umbrella sampling with temperature and pressurecoupling method problem To: Discussion list for GROMACS users ma

[gmx-users] Re: Umbrella sampling with temperature and pressure coupling method problem (Justin A. Lemkul)

2010-10-19 Thread DeChang Li
> Message: 6 > Date: Tue, 19 Oct 2010 09:30:47 -0400 > From: "Justin A. Lemkul" > Subject: Re: [gmx-users] Umbrella sampling with temperature and >pressurecoupling method problem > To: Discussion list for GROMACS users > Message-ID: <4cbd9d87.6080...@vt.edu> > Content-Type: text/p

[gmx-users] Re: Umbrella Sampling

2008-08-04 Thread Alexandr Malafeev
2008/8/4 Vitaly Chaban <[EMAIL PROTECTED]> > Hi Alexandr, > > Have you thoroughly examined this thread as for using umbrella sampling: > http://www.gromacs.org/pipermail/gmx-users/2005-December/018615.html ? > > > -- > Vitaly V. Chaban > School of Chemistry > National University of Kharkiv > Svobo

[gmx-users] RE: Umbrella Sampling

2008-08-04 Thread Vitaly Chaban
Hi Alexandr, Have you thoroughly examined this thread as for using umbrella sampling: http://www.gromacs.org/pipermail/gmx-users/2005-December/018615.html ? -- Vitaly V. Chaban School of Chemistry National University of Kharkiv Svoboda sq., 4, Kharkiv 61077, Ukraine email: [EMAIL PROTECTED] sky

[gmx-users] Re: Umbrella Sampling

2008-01-30 Thread Ilya Chorny
I am attempting to pull a ligand through a channel protein. I have setup the umbrella sampling code and it seems to be working just fine except I get large jumps in the deviation of each pulled group from its restrained position at 1 ps intervals. I am pulling three groups all relative to the backb

[gmx-users] Re: umbrella sampling pull.pdo output

2006-09-20 Thread Manohar Murthi
> > is this perchance a bug or am i doing something > wrong? > sorry, i rtfm. __ Do You Yahoo!? Tired of spam? Yahoo! Mail has the best spam protection around http://mail.yahoo.com ___ gmx-users maili