Re: [gmx-users] Re: g_analyze

2013-11-10 Thread bharat gupta
thank you informing about g_rdf... Is it possible to dump the structure with those average water molecules interacting with the residues. I generated the hbond.log file which gives the details but I need to generate a figure for this ?? On Mon, Nov 11, 2013 at 10:40 AM, Justin Lemkul wrote: >

Re: [gmx-users] Re: g_analyze

2013-11-10 Thread Justin Lemkul
On 11/10/13 8:38 PM, bharat gupta wrote: But trjorder can be used to calculate the hydration layer or shell around residues ... Right ?? Yes, but I also tend to think that integrating an RDF is also a more straightforward way of doing that. With trjorder, you set some arbitrary cutoff tha

Re: [gmx-users] Re: g_analyze

2013-11-10 Thread bharat gupta
But trjorder can be used to calculate the hydration layer or shell around residues ... Right ?? On Mon, Nov 11, 2013 at 10:35 AM, Justin Lemkul wrote: > > > On 11/10/13 8:30 PM, bharat gupta wrote: > >> Thanks for your reply. I was missing the scientific notation part. Now >> everything is fine

Re: [gmx-users] Re: g_analyze

2013-11-10 Thread Justin Lemkul
On 11/10/13 8:30 PM, bharat gupta wrote: Thanks for your reply. I was missing the scientific notation part. Now everything is fine. Regarding trjorder, it doesn't measure h-bonds but gives the water nearest to protein. I wouldn't try to draw any sort of comparison between the output of trjo

Re: [gmx-users] Re: g_analyze

2013-11-10 Thread bharat gupta
Thanks for your reply. I was missing the scientific notation part. Now everything is fine. Regarding trjorder, it doesn't measure h-bonds but gives the water nearest to protein. On Mon, Nov 11, 2013 at 10:12 AM, Justin Lemkul wrote: > > > On 11/10/13 7:18 PM, bharat gupta wrote: > >> I check

Re: [gmx-users] Re: g_analyze

2013-11-10 Thread Justin Lemkul
On 11/10/13 7:18 PM, bharat gupta wrote: I checked the file hbnum.xvg file and it contains three columns - time, hbonds, hbonds that donot follow the angle criteria. In that case SS1 is The third column is not actually H-bonds, then ;) the average of actual hbonds (2nd column ) and SS2 is t

[gmx-users] Re: Thankful

2013-11-10 Thread Williams Ernesto Miranda Delgado
Justin, thank you very much for your kind help about LIE and PME Williams -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! * Please don't post

Re: [gmx-users] Re: g_analyze

2013-11-10 Thread bharat gupta
I checked the file hbnum.xvg file and it contains three columns - time, hbonds, hbonds that donot follow the angle criteria. In that case SS1 is the average of actual hbonds (2nd column ) and SS2 is the average of 3rd column. Am I right here or not ?? I tried to calculate the h-bond for residues 1

[gmx-users] Bilayer COM removal issue: Large VCM

2013-11-10 Thread rajat desikan
Hi All, I am experiencing a few problems in membrane simulations wrt COM removal. I downloaded a 400 ns pre-equilibrated Slipid-DMPC membrane with all the accompanying files. I then carried out the following steps: 1) energy minimization 2) NVT Eq - 100 ps 3) NPT Eq - 250 ps (Berendsen temp, Pres c

Re: [gmx-users] Umbrella Sampling tutorial

2013-11-10 Thread Justin Lemkul
On 11/10/13 9:14 AM, shahab shariati wrote: Dear Justin Very thanks for your reply. What you described earlier should not be attempted with "distance" geometry. It won't work very well. The use of restraints is almost NEVER necessary, especially in the case where the reference group > is mu

[gmx-users] Umbrella Sampling tutorial

2013-11-10 Thread shahab shariati
Dear Justin Very thanks for your reply. > What you described earlier should not be attempted with "distance" > geometry. It won't work very well. The use of restraints is almost > NEVER necessary, especially in the case where the reference group > is much more massive than the pulling group. I w

Re: [gmx-users] problem in running mdrun command

2013-11-10 Thread Mark Abraham
Hi, There's nothing GROMACS-specific here - something about your MPI installation, configuration or use is pretty wrong, but we can't help work out what. Mark On Sun, Nov 10, 2013 at 12:31 PM, S.Chandra Shekar < chandrashe...@iisertvm.ac.in> wrote: > Dear all > > I encounter a problem while ru

Re: [gmx-users] Umbrella Sampling tutorial

2013-11-10 Thread Justin Lemkul
On 11/10/13 5:31 AM, shahab shariati wrote: Dear Justin Thanks for your reply. You are right. I should not extrapolate too literally from your tutorial to my system. But, I have a general question: There is 2 groups in COM pulling method (reference group + pull group). If I want to use pu

Re: [gmx-users] Re: g_analyze

2013-11-10 Thread Justin Lemkul
On 11/10/13 12:20 AM, bharat gupta wrote: Hi, I used the command g_hbond to find h-bond between residues 115-118 and water. Then I used g_analyze to find out the average and it gives the value for the hbonds like this :- std. dev.relative deviation of

[gmx-users] problem in running mdrun command

2013-11-10 Thread S.Chandra Shekar
Dear all I encounter a problem while running command mdrun_mpi -v -deffnm em in gromacs. I am new to the gromacs. i just ran test calculation, *simulation of lyzozyme in water*. i am able to generate gro, tpr files. But in the final step i got following error. Thanks in advance. [localhost.l

[gmx-users] Re: choosing force field

2013-11-10 Thread pratibha
Thank you Justin for your kind help. The simple reason for considering only gromos parameter sets is that the parameters for the metal ions (in my protein) are not defined in other force fields. On Sat, Nov 9, 2013 at 7:18 PM, Justin Lemkul [via GROMACS] < ml-node+s5086n5012376...@n6.nabble.com>

[gmx-users] Umbrella Sampling tutorial

2013-11-10 Thread shahab shariati
Dear Justin Thanks for your reply. You are right. I should not extrapolate too literally from your tutorial to my system. But, I have a general question: There is 2 groups in COM pulling method (reference group + pull group). If I want to use pull_geometry = distance, so, I should fix referen

Re: mdrun on 8-core AMD + GTX TITAN (was: Re: [gmx-users] Re: Gromacs-4.6 on two Titans GPUs)

2013-11-10 Thread Mark Abraham
On Sun, Nov 10, 2013 at 5:28 AM, Dwey Kauffman wrote: > Hi Szilard, > > Thank you very much for your suggestions. > > >Actually, I was jumping to conclusions too early, as you mentioned AMD > >"cluster", I assumed you must have 12-16-core Opteron CPUs. If you > >have an 8-core (desktop?) AMD CPU

Re: [gmx-users] About Compiler Compatibility for Gromacs 4.6.2 Compilation

2013-11-10 Thread Mark Abraham
On Nov 10, 2013 10:04 AM, "Mark Abraham" wrote: > Yes (unless you are using AMD cpus), per the installation instructions, > although you will probably do slightly better with GCC 4.7, and should not > do a new install of 4.6.2 after 4.6.3 is released. In particular, 4.6.2 has > an affinity-relate