[gmx-users] hydrogen bond calculation problem
Dear gromacs users, I have NaCl-water system of finite concentration .I want to calculate water-water hydrogen bond in the first solvation shell around ion..So I used g_hbond -shell option.For single ion in water,its easy to calculate..But say,for water-water hydrogen bonds around two Na+ ions,it is saying to Select one atom for shell. Although I have more than 1 Na+ ions and I want to see water-water hydrogen bond around all Na+ ions simultaneously.Is it possible to take solvation shell of each ion at a same time? Thanks in advance. Regards, *Sushma* -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
[gmx-users] Re: Reaction field zero and ions
Hello If I did the MD simulation using PME and neutralized with ions, and I want to rerun this time with reaction field zero, is there any problem if I keep the ions? This is for LIE calculation. I am using AMBER99SB. Thanks 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
Re: [gmx-users] installation error under openSuse 12.2
On 11/11/13 11:27 AM, kolnkempff wrote: Dear gmx-users, I am trying to install gromacs-4.6.3 on an older Dell that is running openSuse 12.2 Using "DGMX_BUILD_OWN_FFTW=ON" failed for me so to get through cmake I had to compile fftw from scratch and I followed the recommendation of going directly to their website. Now I am at the make stage and get the following message: [ 67%] Built target gmxfftw make[2]: *** No rule to make target `//home/koln/bin/gromacs-4.6.3/build/src/contrib/fftw/gmxfftw-prefix/lib/libfftw3f.a', needed by `src/gmxlib/libgmx.so.8'. Stop. make[1]: *** [src/gmxlib/CMakeFiles/gmx.dir/all] Error 2 make: *** [all] Error 2 When I check for libfftw, I get: koln@linux-5bim:~/bin/gromacs-4.6.3/build2> rpm -qf /usr/lib/libfftw3f.* libfftw3-3-32bit-3.3.2-1.1.2.x86_64 Any suggestions for how to get past this obstacle would be greatly appreciated! Are you building in a clean installation directory? It looks to me that cmake is still trying to build the internal FFTW code, rather than looking for your manually installed FFTW. Providing your full cmake command would also be very useful (hint: always provide it). -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
[gmx-users] segmentation fault on gromacs 4.5.5 after mdrun
Hello evryone, I doing a simulation of a ligand-protein interaction with gromacs 4.5.5. Everything looks fine after I equilibrate the protein-ligand complex. I'm running these commands: grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr mdrun -deffnm nvt Nevertheless, I got this error: "Reading file nvt.tpr, VERSION 4.5.5 (double precision) Segmentation fault" What should I do? Carlos Javier Alméciga Díaz, QF., PhD. Profesor Asistente Pontificia Universidad Javeriana Facultad de Ciencias Instituto de Errores Innatos del Metabolismo Tel: 57-1-3208320 Ext 4140-4099 Fax: 57-1-3208320 Ext 4099 Bogotá. D.C. - COLOMBIA cjalmec...@javeriana.edu.co http://www.javeriana.edu.co/ieim AVISO LEGAL: El presente correo electronico no representa la opinion o el consentimiento oficial de la PONTIFICIA UNIVERSIDAD JAVERIANA. Este mensaje es confidencial y puede contener informacion privilegiada la cual no puede ser usada ni divulgada a personas distintas de su destinatario. Esta prohibida la retencion, grabacion, utilizacion, aprovechamiento o divulgacion con cualquier proposito. Si por error recibe este mensaje, por favor destruya su contenido y avise a su remitente. En este aviso legal se omiten intencionalmente las tildes. Este mensaje ha sido revisado por un sistema antivirus, por lo que su contenido esta libre de virus. This e-mail has been scanned by an antivirus system, so its contents is virus free. -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
Re: [gmx-users] segmentation fault on gromacs 4.5.5 after mdrun
On 11/11/13 11:24 AM, Carlos Javier Almeciga Diaz wrote: Hello evryone, I doing a simulation of a ligand-protein interaction with gromacs 4.5.5. Everything looks fine after I equilibrate the protein-ligand complex. I'm running these commands: grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr mdrun -deffnm nvt Nevertheless, I got this error: "Reading file nvt.tpr, VERSION 4.5.5 (double precision) Segmentation fault" What should I do? Instantaneous failure typically indicates that the forces are nonsensically high and the constraint algorithm immediately fails. Likely the previous energy minimization did not adequately complete. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
[gmx-users] installation error under openSuse 12.2
Dear gmx-users, I am trying to install gromacs-4.6.3 on an older Dell that is running openSuse 12.2 Using "DGMX_BUILD_OWN_FFTW=ON" failed for me so to get through cmake I had to compile fftw from scratch and I followed the recommendation of going directly to their website. Now I am at the make stage and get the following message: [ 67%] Built target gmxfftw make[2]: *** No rule to make target `//home/koln/bin/gromacs-4.6.3/build/src/contrib/fftw/gmxfftw-prefix/lib/libfftw3f.a', needed by `src/gmxlib/libgmx.so.8'. Stop. make[1]: *** [src/gmxlib/CMakeFiles/gmx.dir/all] Error 2 make: *** [all] Error 2 When I check for libfftw, I get: koln@linux-5bim:~/bin/gromacs-4.6.3/build2> rpm -qf /usr/lib/libfftw3f.* libfftw3-3-32bit-3.3.2-1.1.2.x86_64 Any suggestions for how to get past this obstacle would be greatly appreciated! Koln -- View this message in context: http://gromacs.5086.x6.nabble.com/installation-error-under-openSuse-12-2-tp5012430.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
Re: [gmx-users] help about logfile
On 11/11/13 5:04 AM, guozhicheng222 wrote: Hi I am confusing with the output file (.log) about the sample frequency (frames) in my simulation. The average information in .log file is 'Statistics over 31 steps using 3001 frames' where nstxout =4000 and nstlog =4000. While, 'Statistics over 31 steps using 20001 frames', appeared in .log file, where nstxout =15 and nstlog =15. In my opinion, the former sample frequency should be 75 frames rather than 3001 frames and the latter is right. I want to know how can I control the sample frequency, arbitrarily. Check to make sure you set nstlog correctly in the first run. Look for the value in the .log file itself. It is highly unlikely that something so fundamental and simple is being executed incorrectly in one run, but correctly in another. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
Re: [gmx-users] Re: g_analyze
On 11/11/13 4:06 AM, bharat gupta wrote: In addition to my previous question, I have another question about g_analyze. When I used the hbond.xvg file to get the average and plotted the average.xvg file I found that the average value is round 4 to 5 according to the graph. But g_analyze in its final calculation gives 7.150 as the average values... Here's the link for the graph and result of average value calculated by g_analyze :- std. dev.relative deviation of standard - cumulants from those of set average deviation sqrt(n-1) a Gaussian distribition cum. 3 cum. 4 SS1 * 7.150740e+00 * 8.803173e-01 1.760635e-02 0.0620.163 SS2 1.490604e+00 1.164761e+00 2.329523e-02 0.4950.153 https://www.dropbox.com/s/1vqixenyerha7qq/115-water.png Here's the link hbond.xvg file and its averaged file https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg Neither of these files produce output that corresponds to the PNG image above. Both files have values in 6-9 H-bond range and thus agree with the g_analyze output, which I can reproduce. I suspect you're somehow getting your files mixed up. -Justin On Mon, Nov 11, 2013 at 3:30 PM, bharat gupta wrote: 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: 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 that may or may not be an informed decision - with an RDF it is clear where the hydration layers are. -Justin 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. 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 trjorder and g_hbond. If you want to measure H-bonds, there's only one tool for that. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
Re: [gmx-users] Re: g_analyze
On 11/11/13 1:30 AM, bharat gupta wrote: 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 ?? g_select -Justin On Mon, Nov 11, 2013 at 10:40 AM, Justin Lemkul wrote: 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 that may or may not be an informed decision - with an RDF it is clear where the hydration layers are. -Justin 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. 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 trjorder and g_hbond. If you want to measure H-bonds, there's only one tool for that. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
[gmx-users] help about logfile
Hi I am confusing with the output file (.log) about the sample frequency (frames) in my simulation. The average information in .log file is 'Statistics over 31 steps using 3001 frames' where nstxout =4000 and nstlog =4000. While, 'Statistics over 31 steps using 20001 frames', appeared in .log file, where nstxout =15 and nstlog =15. In my opinion, the former sample frequency should be 75 frames rather than 3001 frames and the latter is right. I want to know how can I control the sample frequency, arbitrarily.-- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
Re: [gmx-users] Re: g_analyze
In addition to my previous question, I have another question about g_analyze. When I used the hbond.xvg file to get the average and plotted the average.xvg file I found that the average value is round 4 to 5 according to the graph. But g_analyze in its final calculation gives 7.150 as the average values... Here's the link for the graph and result of average value calculated by g_analyze :- std. dev.relative deviation of standard - cumulants from those of set average deviation sqrt(n-1) a Gaussian distribition cum. 3 cum. 4 SS1 * 7.150740e+00 * 8.803173e-01 1.760635e-02 0.0620.163 SS2 1.490604e+00 1.164761e+00 2.329523e-02 0.4950.153 https://www.dropbox.com/s/1vqixenyerha7qq/115-water.png Here's the link hbond.xvg file and its averaged file https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg https://www.dropbox.com/s/4n0m47o3mrjn3o8/hbond_115-water.xvg On Mon, Nov 11, 2013 at 3:30 PM, bharat gupta wrote: > 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: > >> >> >> 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 that may or may not be an informed decision - with an RDF it is >> clear where the hydration layers are. >> >> -Justin >> >> >> >>> 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. > > 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 trjorder and g_hbond. If you want to measure H-bonds, there's only one tool for that. -Justin -- == Justin A. Lemkul, Ph.D. Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 601 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 == -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists >> -- >> == >> >> Justin A. Lemkul, Ph.D. >> Postdoctoral Fellow >> >> Department of Pharmaceutical Sciences >> School of Pharmacy >> Health Sciences Facility II, Room 601 >> University of Maryland, Baltimore >> 20 Penn St. >> Baltimore, MD 21201 >> >> jalem...@outerbanks.umaryland.edu | (410) 706-7441 >> >> == >> -- >> 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 (un)subscribe requests to the list. Use the www >> interface or send it to gmx-users-requ...@gromacs.org. >> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists >> > > > > -- 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 (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
