Re: [gmx-users] umbrella sampling tutorial

On 7/6/18 10:42 AM, hosein geraili wrote: Dear all, I have a question regarding the umbrella sampling tutorial. I am using the bash script to create the summary_distances.dat, but the problem is at first I would get the error "Unknown command-line option -o", then I used -oh instead of -o

Re: [gmx-users] Umbrella sampling tutorial WHAM issue

On 11/10/16 10:46 AM, Lindsay, Richard J. wrote: Hi, I am attempting to follow the gromacs umbrella sampling tutorial but have run into an issue at the data analysis step. When I run the gmx wham command I get only one histogram, even though I collected 24 umbrella trajectories. You're

Re: [gmx-users] Umbrella sampling tutorial

Any idea on this topic? Thanks > On 07 Oct 2016, at 14:13, gozde ergin wrote: > > Thanks Justin. > I want infinite dilute bulk that’s why I do not want any surfactant molecule > inside the bulk. > I may put a ‘’wall’’ just below the surface molecules in order to inhibit

Re: [gmx-users] Umbrella sampling tutorial

Thanks Justin. I want infinite dilute bulk that’s why I do not want any surfactant molecule inside the bulk. I may put a ‘’wall’’ just below the surface molecules in order to inhibit their travel to bulk but I do not know how to remove this effect from PMF. > On 07 Oct 2016, at 13:40, Justin

Re: [gmx-users] Umbrella sampling tutorial

On 10/7/16 4:41 AM, gozde ergin wrote: Dear all, I would like to ask a question about Justin’s tutorial. By using the POSRES_B, a restraint is applied on chain B in order to immobile it. Is the effect of this restraint removed during the WHAM? Not to my knowledge, no. The restraints in

Re: [gmx-users] Umbrella sampling Tutorial by Dr. Lemkul

On 6/7/16 3:13 AM, amitbe...@chemeng.iisc.ernet.in wrote: Hello, I was running the simulation for the umbrella sampling tutorial by Dr. Justin Lemkul in GROMACS 5.1.2 and I followed the steps exactly as they were mentioned. But in my summary_distances.dat there is one column. the COM column is

Re: [gmx-users] umbrella sampling tutorial

March 2015 9:45 PM To: gmx-us...@gromacs.org Subject: Re: [gmx-users] umbrella sampling tutorial On 3/17/15 9:06 PM, Ming Tang wrote: Hi Justin, Thanks a lot. I can move forward now. Recently, I am trying to pull a triple helix, and want to fix the center of mass of the three terminal atoms

Re: [gmx-users] umbrella sampling tutorial

On 4/7/15 7:11 AM, Ming Tang wrote: Hi, Justin I am working with a triple helix using umbrella sampling method. When generating configurations, I pulled 1ns (md_pull.mdp). Then I got 15 configurations using 0.03 spacing. They are conf0.gro conf5.gro conf10.gro conf15.gro conf17.gro conf18.gro

Re: [gmx-users] umbrella sampling tutorial

On 3/17/15 9:06 PM, Ming Tang wrote: Hi Justin, Thanks a lot. I can move forward now. Recently, I am trying to pull a triple helix, and want to fix the center of mass of the three terminal atoms. One paper said that it can be done by means of a strong harmonic restrain with a super large

Re: [gmx-users] umbrella sampling tutorial

[mailto:gromacs.org_gmx-users-boun...@maillist.sys.kth.se] On Behalf Of Justin Lemkul Sent: Monday, 16 March 2015 9:27 PM To: gmx-us...@gromacs.org Subject: Re: [gmx-users] umbrella sampling tutorial On 3/15/15 9:57 PM, Ming Tang wrote: Dear Justin, I have tried many cases, and only your code runs smoothly

Re: [gmx-users] umbrella sampling tutorial

-Original Message- From: gromacs.org_gmx-users-boun...@maillist.sys.kth.se [mailto:gromacs.org_gmx-users-boun...@maillist.sys.kth.se] On Behalf Of Justin Lemkul Sent: Monday, 16 March 2015 9:27 PM To: gmx-us...@gromacs.org Subject: Re: [gmx-users] umbrella sampling tutorial On 3/15/15 9:57 PM

Re: [gmx-users] umbrella sampling tutorial

...@gromacs.org Subject: Re: [gmx-users] umbrella sampling tutorial On 3/17/15 8:08 PM, Ming Tang wrote: Hi Justin, Here is the error message. tm@tm-HP-Z420-Workstation:~/umbrella$ perl distances.pl Bareword found where operator expected at distances.pl line 8, near gmx distance -n

Re: [gmx-users] umbrella sampling tutorial

...@maillist.sys.kth.se] On Behalf Of Justin Lemkul Sent: Monday, 16 March 2015 9:27 PM To: gmx-us...@gromacs.org Subject: Re: [gmx-users] umbrella sampling tutorial On 3/15/15 9:57 PM, Ming Tang wrote: Dear Justin, I have tried many cases, and only your code runs smoothly. If I set pull-ngroups = 1

Re: [gmx-users] umbrella sampling tutorial

...@maillist.sys.kth.se [mailto:gromacs.org_gmx-users-boun...@maillist.sys.kth.se] On Behalf Of Ming Tang Sent: Wednesday, 18 March 2015 8:14 AM To: gmx-us...@gromacs.org Subject: Re: [gmx-users] umbrella sampling tutorial Yes, it can not run. But I can extract the COM of different conf.groone by one manually

Re: [gmx-users] umbrella sampling tutorial

March 2015 8:14 AM To: gmx-us...@gromacs.org Subject: Re: [gmx-users] umbrella sampling tutorial Yes, it can not run. But I can extract the COM of different conf.groone by one manually. Due to the time difference, I will update the error message when I get to the office. Thanks. Sent from my Huawei

Re: [gmx-users] umbrella sampling tutorial

On 3/15/15 9:57 PM, Ming Tang wrote: Dear Justin, I have tried many cases, and only your code runs smoothly. If I set pull-ngroups = 1, it will warn that pull-group2-name= Chain_A has not been defined. The code you gave me last time runs well. Though the COM between Chain_A and Chain_B

Re: [gmx-users] umbrella sampling tutorial

...@maillist.sys.kth.se [mailto:gromacs.org_gmx-users-boun...@maillist.sys.kth.se] On Behalf Of Justin Lemkul Sent: Friday, 13 March 2015 10:57 AM To: gmx-us...@gromacs.org Subject: Re: [gmx-users] umbrella sampling tutorial On 3/12/15 8:49 PM, Ming Tang wrote: Dear Justin, I switched Chain_A and Chain_B

Re: [gmx-users] umbrella sampling tutorial

[mailto:gromacs.org_gmx-users-boun...@maillist.sys.kth.se] On Behalf Of Ming Tang Sent: Thursday, 12 March 2015 11:34 AM To: gmx-us...@gromacs.org Subject: Re: [gmx-users] umbrella sampling tutorial Hi Justin, Your code works! Thank you very much. Regards, Ming -Original Message- From

Re: [gmx-users] umbrella sampling tutorial

On 3/11/15 4:17 AM, Ming Tang wrote: Dear all, I was learning the Umbrella Sampling tutorial step by step on Gromacs 5.0.4. Everything was fine until I got to the step five (run the continuous pulling simulation). Then I turned to the 5.0.4 manual, and tried to modify the pull code,

Re: [gmx-users] umbrella sampling tutorial

Message- From: gromacs.org_gmx-users-boun...@maillist.sys.kth.se [mailto:gromacs.org_gmx-users-boun...@maillist.sys.kth.se] On Behalf Of Justin Lemkul Sent: Thursday, 12 March 2015 9:50 PM To: gmx-us...@gromacs.org Subject: Re: [gmx-users] umbrella sampling tutorial On 3/12/15 2:41 AM, Ming

Re: [gmx-users] umbrella sampling tutorial

On 3/12/15 8:49 PM, Ming Tang wrote: Dear Justin, I switched Chain_A and Chain_B in the .mdp file, but got the same result. In the tutorial: pull-ngroups: (1) The number of pull groups, not including the absolute reference group, when used. And we need to assign the Chain_B as the reference

Re: [gmx-users] umbrella sampling tutorial

On 3/12/15 2:41 AM, Ming Tang wrote: Dear Justin, The results turned out that the COM distance between them is the same over time, seems that both Chain_A and Chain_B had been pulled by the same force in the same direction, which is different from the tutorial. Switch Chain_A and

Re: [gmx-users] umbrella sampling tutorial

On 3/11/15 8:58 PM, Ming Tang wrote: Dear Justin, Thanks for your help. I changed the group names to Chain_A and Chain_B in my index file, and modified the COM pulling code. ; COM pulling pull= umbrella pull_geometry = distance ; simple distance increase pull_dim= N N

Re: [gmx-users] umbrella sampling tutorial

On 3/11/15 4:17 AM, Ming Tang wrote: Dear all, I was learning the Umbrella Sampling tutorial step by step on Gromacs 5.0.4. Everything was fine until I got to the step five (run the continuous pulling simulation). Then I turned to the 5.0.4 manual, and tried to modify the pull code, but

Re: [gmx-users] umbrella sampling tutorial

...@gromacs.org Subject: Re: [gmx-users] umbrella sampling tutorial On 3/11/15 8:58 PM, Ming Tang wrote: Dear Justin, Thanks for your help. I changed the group names to Chain_A and Chain_B in my index file, and modified the COM pulling code. ; COM pulling pull= umbrella pull_geometry

Re: [gmx-users] Umbrella Sampling Tutorial Question

On 10/31/14 8:48 PM, Agnivo Gosai wrote: Dear Users Firstly thanks to Dr. Lemkul for his immensely informative and helpful tutorials. I am following the umbrella sampling tutorial for one of my projects. During the placement of the molecule inside the box , the tutorial tells us to visualize

Re: [gmx-users] Umbrella Sampling Tutorial

On 1/11/14, 11:50 AM, Andres Ortega wrote: Hi, Dear Gromacs Users I was running de Umbrella Sampling Tutorial, and when i get the frames from the trajectory, there are some frames like this: 489 5.3776196 490 5.3817739 491 5.3752654 492 5.3619911 493 5.3890083 494

Re: [gmx-users] Umbrella Sampling tutorial

On 12/1/13 8:04 AM, shahab shariati wrote: Dear Justin Very thanks for your reply I will slow pulling by changing pull_rate1 from 0.01 to 0.001 and pull_k1 from 1000 to 100. I have a basic question. Before minimization, where should I put drug molecule? into water molecules or into lipid

Re: [gmx-users] Umbrella Sampling tutorial

On 12/1/13 10:38 AM, shahab shariati wrote: Dear Justin When I placed drug within the lipid bilayer, after Generating Configurations step, I obtained 501 gro files. When I used distance.pl script for these gro files, my summary_distances.dat file is as follows: 00.2756352 10.2535359

Re: [gmx-users] Umbrella Sampling tutorial

The output you got indicates that the pulled group came to the location where it is expected by gromacs to be. The input looks correct, but gromacs' pull code is with surprise inside. I would try to pull in opposite direction, i.e. pull_vec1 = -0 0 1. Turning off barostat during pulling also

Re: [gmx-users] Umbrella Sampling tutorial

On 11/30/13 9:22 AM, shahab shariati wrote: Dear Justin My system contains lipid bilayer + drug + water molecules. I want to calculate Potential of mean force as a function of the distance between the centers of mass of drug and the lipid bilayer. Based on your suggestion I used