[gmx-users] How can I determine how Gromacs was compiled?
Hi, I'm running Gromacs 4.5.5. Is there a way to determine how Gromacs was compiled on my system? I would like to know whether the person who compiled Gromacs used gcc or icc (both compilers are available and were presumably available when the compilation occurred) and, if possible, what version of FFTW they used. Is there a file or a command that has this information? Thank you so much! Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] Does Gromacs take advantage of ECC memory?
Hi, I'm looking into a workstation which I will use to run some Gromacs simulations. I'm a chemist, not a computer person, so I'm not very well-versed in computer hardware. I have the option of getting non-ECC or ECC memory, with the ECC being at extra cost. In asking around about whether I need ECC, I've been advised to "research your simulation apps and see if they can take advantage of ECC." I don't see ECC mentioned in the manual. Does Gromacs take advantage of ECC? The simulations I will be running will not be terribly large: on the order of ~10,000 atoms for a few tens of ns. My system will likely have 8 cores. Thanks so much for your time. Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] Pausing and resuming mdrun
Hi, I'm running an MD simulation using mdrun (specifically version 4.5.5 -- an old version). I would like to pause this simulation -- perhaps for days or weeks -- so that I can run a different one on this node. Of course, one way is to just kill the simulation altogether: kill $ProcessID where $ProcessID is the process ID of the simulation. Then, when I want to resume the simulation, I can just pass the last checkpoint file to mdrun. Checkpoint files have been written every 15 minutes (i.e., the default setting), so with this method I will lose at most 15 minutes of computation time. But, is there any way to literally _pause_ the simulation and resume it a few days or weeks later? A Unix/Linux question and answer site ( http://unix.stackexchange.com/questions/2107/how-to-suspend-and-resume-processes ) says that one can pause/resume a process with this method: kill -SIGSTOP $ProcessID kill -SIGCONT $ProcessID or: kill -SIGSTP $ProcessID kill -SIGCONT $ProcessID Another site ( http://www.cyberciti.biz/faq/unix-kill-command-examples/ ) says to just use: kill -STOP $ProcessID kill -CONT $ProcessID My question is, do you think that these Linux methods will work with mdrun? I did a test with mdrun on another node, and it seems to work, but I'm just wondering if there are any dangers in using these methods. (I am accessing the Linux machine remotely, by SSH. Sometimes my SSH connection gives out, so when starting a simulation I always use "nohup mdrun -s topol.tpr" so that the mdrun process is not terminated when my SSH connection flakes out. I'm not sure if this will affect the viability of the "kill -STOP/CONT" method...) Thanks so much, Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
Re: [gmx-users] Pausing and resuming mdrun
Szilard, Thanks for your time; it is really helpful. Just in case people are reading this thread in the future, I want to point out a mistake I made in my original message. I said that using the signal "-SIGSTP" was an alternative to -SIGSTOP. It is actually "-SIGTSTP", not "-SIGSTP". -SIGTSTP is short for "_t_erminal stop," as this reference page explains: http://www.ucs.cam.ac.uk/docs/course-notes/unix-courses/Building/files/signals.pdf Thanks again, Andrew -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] What is Mnbf/s?
Hi, I feel like this has probably been asked and answered before, but I can't seem to find it. Even if I Google for Mnbf/s, I can't seem to find a definition. In the .log output, what is Mnbf/s? Thanks so much for your time, Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] Formatting of spacing in .top and .itp files
Hi, I am running Gromacs 4.5.5. In that rather old version, there is a binary called g_x2top -- I think it is now called something else -- to generate "a primitive topology from a coordinate file." When I use g_x2top with my (custom) force field, g_x2top generates a nicely formatted topology file (either .top or .itp). Here is a snippet from the nicely formatted topology (assuming you display the text in a fixed width font): [ dihedrals ] ; aiajakal functc0c1 c2c3c4c5 8197 113 129 1 11297 113 114 1 11297 113 129 1 82989983 1 829899 115 1 (Only atom indices and a function type are printed because I used the option -noparam in g_x2top.) But, what if I want to generate my own .top or .itp text file, using for example Fortran, Perl, or Python? Are .top and .itp fixed-spaced file formats? In other words, do I have to nicely align the columns as g_x2top does? Or can I just use any number of spaces to delimit the columns, like in the following example? [ dihedrals ] 81 97 113 129 1 112 97 113 114 1 112 97 113 129 1 82 98 99 83 1 82 98 99 115 1 I have tried this ugly formatting, and it seems to give the same results as the nice formatting. But can anyone confirm that .top and .itp are NOT fixed-spaced file formats (unlike .pdb and .gro for example)? Thanks for your time, Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] Multiple parameter entries in itp with the same name
Hi, If I have a force field parameter file (for example, ffnonbonded.itp) with multiple entries with the same name (for example, multiple [ atomtypes ] entries with the name "CA"), which one is used by grompp: the first one or the last one? My guess is that the last one listed is used, but I guess I can do a gmxdump of the .tpr file to find out for sure. Thanks so much for your time. Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] Questions about 1-4 interactions in alkanes
Hi, Section 4.2 of the manual has a discussion of dihedrals. The section contains, in part, the following text: "With the periodic GROMOS potential a special 1-4 LJ-interaction must be included; with the Ryckaert-Bellemans potential _for alkanes_ the 1-4 interactions must be excluded from the non-bonded list. Note: Ryckaert-Bellemans potentials are also used in e.g. the OPLS force field in combination with 1-4 interactions. You should therefore not modify topologies generated by pdb2gmx in this case." I would like to ask in particular about the statement that "with the Ryckaert-Bellemans potential _for alkanes_ the 1-4 interactions must be excluded from the non-bonded list." Here are two questions about that statement: (1) Does that statement apply ONLY to the GROMOS force field, or does it apply to the OPLS-AA force field as well? (2) What is the meaning of "alkanes" in that statement? - (2a) Does it refer ONLY to saturated hydrocarbons, which contain only hydrogen and carbon and contain only single covalent bonds? (Examples of such molecules are butane, isobutane, pentane, isopentane, and hexane.) - (2b) Or, does it refer to ANY molecule which contains an alkyl group? (An example of such a molecule is the amino acid isoleucine. Isoleucine contains an alkyl group -- isoleucine's side chain is an isopentyl group. Does this mean that I need to EXCLUDE 1-4 interactions among hydrogens and carbons of the isopentyl group from the non-bonded list?) (3) In the OPLS-AA force field, nexcl is usually set to 3, which _initially_ excludes 1-4 interactions. However, in OPLS-AA, 1-4 interactions are explicitly specified in the [ pairs ] section; they are taken at half-strength by setting fudgeQQ = 0.5 and fudgeLJ = 0.5 in the topology. (See, for example, this helpful explanation: https://mailman-1.sys.kth.se/pipermail/gromacs.org_gmx-users/2012-December/0 77005.html .) If the statement applies to OPLS-AA in addition to GROMOS, and if the statement applies to the alkyl atoms in, e.g. isoleucine, then do I need to EXCLUDE 1-4 interactions among isoleucine's hydrogens and carbons from the [ pairs ] section? Thank you for your time! Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] What does PHE(OL) mean?
Hi, This might be more of a biochemistry question rather than a Gromacs question, so I apologize in advance... In the "top" directory (e.g., on my installation, it's /usr/local/gromacs/share/gromacs/top/) of the Gromacs installation, several force fields are included. I'm looking at the file "ffbonded.itp" in the "oplsaa.ff" directory. I see, for example, this entry under [ angles ]: CA CA CA 1 120.000527.184 ; PHE(OL) PHE, I'm sure, refers to phenylalanine, where CA is an aromatic carbon. But to what does (OL) refer? I see other (OL) entries in the topology: for example, TYR(OL), TRP(OL), and ARG(OL)? Does (OL) somehow refer to the protonation state of the amino acid? Thanks so much for your time! Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] How can I tell whether or not Gromacs was compiled with MPI?
Hi, I am running Gromacs version 4.5.5. However, I did not compile it myself; the system administrator did many months ago. I have been running "standard" MD simulations (i.e., no replica exchange) over the 24 cores that are available on the node. To run "standard" MD, I just call mdrun directly -- without referencing MPI or mpirun (at least not explicitly) -- using a command like this: mdrun -v -s topol.tpr -nt 24 >From looking at the old release notes, it looks like this "thread-based parallelization" was introduced sometime in version 4.5 (maybe from the beginning of version 4.5?). As the page http://www.gromacs.org/About_Gromacs/Release_Notes/Versions_4.5.x says, "Running on a multi-core node now uses thread-based parallelization to automatically spawn the optimum number of threads in the default build. MPI is now only required for parallelization over the network." Now, however, I want to run replica exchange simulations. I have no experience with such simulations, but the page http://www.gromacs.org/Documentation/How-tos/REMD says: "REMD with any version of GROMACS requires mdrun compiled with MPI (i.e. not threading), and that the number of processors be a multiple of the number of replicas." But, my question is, how can I tell whether Gromacs was compiled with MPI? Probably it was not compiled with MPI, since I am using threading, but how can I be sure? Thanks so much for your time! Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
Re: [gmx-users] About g_potential
Hi Guangwei, The simulation box is binned, or "sliced," along z in planes parallel to the xy plane. A greater number of slices means the box is sliced "finer" -- the bins or "slices" are smaller. The charges are binned (like making a histogram) into those slices, and the potential is calculated by integrating the charge density. So using more slices means more resolution, but using too many slices can lead to irregularities if you haven't saved configurations very often in your trajectory. Andrew DeYoung -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] Mailing list archive
Hi, Does anyone know if the gmx-users Archive has been moved or is temporarily down? When I do to: http://lists.gromacs.org/pipermail/gmx-users/ it redirects me to: https://mailman-1.sys.kth.se/pipermail/gmx-users/ but that page says, "The requested URL /pipermail/gmx-users/ was not found on this server." Thanks! Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
Re: [gmx-users] Mailing list archive
Hi again, I just realized that Mark sent out an email earlier today announcing a power outage for most of the Gromacs infrastructure, starting at midnight. I think it is almost midnight in Stockholm right now, so probably the outage has started and that's why I'm not able to access http://lists.gromacs.org/pipermail/gmx-users/ . Sorry for not noticing this earlier. I will try accessing the archive in a couple of days and will be in touch if I still have trouble. Thanks! Andrew DeYoung Carnegie Mellon University -Original Message----- From: Andrew DeYoung Sent: Tuesday, November 26, 2013 4:08 PM To: gmx-us...@gromacs.org Subject: Mailing list archive Hi, Does anyone know if the gmx-users Archive has been moved or is temporarily down? When I do to: http://lists.gromacs.org/pipermail/gmx-users/ it redirects me to: https://mailman-1.sys.kth.se/pipermail/gmx-users/ but that page says, "The requested URL /pipermail/gmx-users/ was not found on this server." Thanks! Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
Re: [gmx-users] Mailing list archive
Hi Justin, Thank you so much for your help with the mailing list archive URLs! This really helps me, because I was trying to access an old mailing list conversation. With your tip, now I can; thanks! Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] What does it mean that particle decomposition does not support checkpointing?
Hi, I am a Gromacs user (I'm running version 4.5.5). I am not _at all_ an experienced programmer, but I was looking through this redmine page about "mdrun features to deprecate for 5.0": http://redmine.gromacs.org/issues/1292 Near the top of that page (for issue #1292), it mentions particle decomposition as a feature to deprecate for 5.0: "1) Particle decomposition Doesn't support checkpointing. Doesn't scale because it does not balance load. Removing it will simplify a lot of conditionally executed code paths. The main caveat is whether there are important algorithms whose implementations still require it (see below)." What does it mean that particle decomposition doesn't support checkpointing? I use particle decomposition for my Gromacs simulations (in version 4.5.5), and I also save checkpoints from mdrun every 15 minutes using the -cpo flag. The resulting .cpt files seem to work fine -- at least "on the surface"; I have not checked this in detail -- for doing a restart whenever my system crashes. As a test, I took the final checkpoint file of a run. Then, using the -t flag in mdrun, I started a new simulation using that final checkpoint file. After that new simulation finished, I used trjconv -dump 0 to dump the first frame into a configuration file (in my case, a .g96 file). I compared this to the coordinates in the checkpoint file (as read by gmxdump -cp checkpoint.cpt), and the coordinates in the dumped configuration and in the checkpoint were the same up to the 0.1 nm decimal place (since .g96 files and gmxdump output have different precision). So it seems that, at least in some ways, checkpointing is working even when particle decomposition is used. If so, then what does it mean in the above redmine post that "particle decomposition doesn't support checkpointing"? If not, why does my very quick, superficial check described in the previous paragraph appear to be successful? Thanks so very much for your time! Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
Re: [gmx-users] electrical potential using Gromacs
Hi, I was able to "fix" the drift in the electric potential by using particle decomposition parallelization instead of domain decomposition when running the simulations using mdrun (although I still do not know why this "fixes" the problem). Domain decomposition is the default parallelization in Gromacs, but particle decomposition can be used with the -pd switch of mdrun. However, I think that particle decomposition is slated to be removed in future versions of Gromacs. (I also did not use the -correct switch of g_potential, but that alone did not fix the problem.) Happy New Year. Andrew -Original Message- From: Fereshte Taherian [mailto:tab...@csi.tu-darmstadt.de] Sent: Friday, December 27, 2013 8:04 PM To: adeyo...@andrew.cmu.edu Subject: electrical potential using Gromacs Dear Andrew, I am Fereshte Taherian, and doing my Phd in Technische Universität Darmstadt, Germany. Sorry to bother you. Currently I am using Gromacs to calculate the electrical potential of a ionic liquid confined between two charged surfaces (ewald_geometry = 3dc). Using g_potential I tried to calculate the potential as following: g_potential -correct Looking at the electrical potential, I found a drift in the potential as you mentioned in the Gromacs users forum (http://gromacs.5086.x6.nabble.com/Pseudo-2-D-Ewald-Summation-and-Slab-Geome try-td4880471.html). I really appreciate, if you could let me know how did you fix this problem. Thank you very much. I wish you a happy new year. Best regards, Fereshte -- Fereshte Taherin Tabasi PhD student of Computational Physical Chemistry Center of Smart Interfaces Technische Universität Darmstadt Petersenstraße 32 64287 Darmstadt Tel:(+49)6151-1675040 (Office) Fax:(+49)6151-162048(university) Email: tab...@csi.tu-darmstadt.de http://www.csi.tu-darmstadt.de -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
Re: [gmx-users] electrical potential using Gromacs
Hi Fereshte, I used PME. Andrew -Original Message- From: Fereshte Taherian [mailto:tab...@csi.tu-darmstadt.de] Sent: Saturday, January 04, 2014 5:29 PM To: Andrew DeYoung Subject: RE: electrical potential using Gromacs Hi Andrew, Many Thanks for your reply. I am checking my system based on your suggestion. just one additional question, could you please let me know how did you calculate the coulombic interaction in your system, with PME or the other methods? Thanks again, Fereshte Quoting "Andrew DeYoung" : > Hi, > > I was able to "fix" the drift in the electric potential by using particle > decomposition parallelization instead of domain decomposition when running > the simulations using mdrun (although I still do not know why this "fixes" > the problem). Domain decomposition is the default parallelization in > Gromacs, but particle decomposition can be used with the -pd switch of > mdrun. However, I think that particle decomposition is slated to be removed > in future versions of Gromacs. > > (I also did not use the -correct switch of g_potential, but that alone did > not fix the problem.) > > Happy New Year. > > Andrew > > -Original Message- > From: Fereshte Taherian [mailto:tab...@csi.tu-darmstadt.de] > Sent: Friday, December 27, 2013 8:04 PM > To: adeyo...@andrew.cmu.edu > Subject: electrical potential using Gromacs > > > > Dear Andrew, > > I am Fereshte Taherian, and doing my Phd in Technische Universität > Darmstadt, Germany. Sorry to bother you. Currently I am using Gromacs > to calculate the electrical potential of a ionic liquid confined > between two charged surfaces (ewald_geometry = 3dc). > > Using g_potential I tried to calculate the potential as following: > > g_potential -correct > > Looking at the electrical potential, I found a drift in the potential > as you mentioned in the Gromacs users forum > (http://gromacs.5086.x6.nabble.com/Pseudo-2-D-Ewald-Summation-and-Slab-Geome > try-td4880471.html). > > I really appreciate, if you could let me know how did you fix this problem. > > Thank you very much. > > I wish you a happy new year. > > Best regards, > Fereshte > > > > > > > -- > Fereshte Taherin Tabasi > PhD student of Computational Physical Chemistry > Center of Smart Interfaces > Technische Universität Darmstadt > Petersenstraße 32 > 64287 Darmstadt > Tel:(+49)6151-1675040 (Office) > Fax:(+49)6151-162048(university) > > Email: tab...@csi.tu-darmstadt.de > http://www.csi.tu-darmstadt.de > > > -- Fereshte Taherin Tabasi PhD student of Computational Physical Chemistry Center of Smart Interfaces Technische Universität Darmstadt Petersenstraße 32 64287 Darmstadt Tel:(+49)6151-1675040 (Office) Fax:(+49)6151-162048(university) Email: tab...@csi.tu-darmstadt.de http://www.csi.tu-darmstadt.de -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] Supported features in Gromacs on GPU
Hi, It seems that Gromacs support for GPU began with OpenMM (which works with Gromacs 4.5). But starting with Gromacs 4.6, GPU support was native. I am interested in which features are currently implemented for GPU. (I know that, given the staggering amount of code and years and years of work that have gone into Gromacs, it is a huge task to implement a lot of features on GPU!) The page below lists some features supported in the OpenMM GPU-enabled Gromacs (i.e., 4.5): http://www.gromacs.org/Documentation/Installation_Instructions_4.5/GROMACS-O penMM#Supported_features But this page refers specifically to the older implementation (i.e., OpenMM). My question is, have any of the features unsupported in the OpenMM been implemented in the native GPU implementation (i.e., 4.6)? The following page describes the native GPU implementation in general, but I don't see specific references to features supported: http://www.gromacs.org/GPU_acceleration I am particularly interested in these features: (1) ability to compile in double precision (2) leapfrog integration (i.e., integrator = md; in OpenMM, only velocity Verlet is available) (3) Nose-Hoover thermostat (in OpenMM, only Andersen is available) (4) Ewald sum with dipole correction (i.e., ewald-geometry = 3dc) (5) free energy calculations Do you know if there have been any changes in support for any of (1)-(5) above between versions 4.5 and 4.6? But, again, I know that the developers are working hard and that it must be an extreme amount of work to implement things on the GPU. I'm not complaining! Thanks so much, Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
Re: [gmx-users] Supported features in Gromacs on GPU
Hi Mark, Thank you so much for your kind response! If you have time, may I please ask a follow-up question? In your response below, you said that three of the features that I asked about (leapfrog integration, Nose-Hoover thermostat, and Ewald sum with dipole correction) are, in fact, now supported on the GPU. Were you referring to version 4.6.x or to 5.0? Thanks so much! Andrew DeYoung Carnegie Mellon University > Hi, > > It seems that Gromacs support for GPU began with OpenMM (which works with > Gromacs 4.5). But starting with Gromacs 4.6, GPU support was native. > > I am interested in which features are currently implemented for GPU. (I > know that, given the staggering amount of code and years and years of work > that have gone into Gromacs, it is a huge task to implement a lot of > features on GPU!) The page below lists some features supported in the > OpenMM GPU-enabled Gromacs (i.e., 4.5): > > > http://www.gromacs.org/Documentation/Installation_Instructions_4.5/GROMACS-O penMM#Supported_features > > But this page refers specifically to the older implementation (i.e., > OpenMM). My question is, have any of the features unsupported in the > OpenMM > been implemented in the native GPU implementation (i.e., 4.6)? The > following page describes the native GPU implementation in general, but I > don't see specific references to features supported: > > http://www.gromacs.org/GPU_acceleration > > I am particularly interested in these features: > > (1) ability to compile in double precision > No - see archive for post from Szilard Pall on this topic, but in short the combination of needing a high-end card for actual good performance and significant developer work that wouldn't be used by the core team is a bit of deal breaker. > (2) leapfrog integration (i.e., integrator = md; in OpenMM, only velocity > Verlet is available) > (3) Nose-Hoover thermostat (in OpenMM, only Andersen is available) > (4) Ewald sum with dipole correction (i.e., ewald-geometry = 3dc) > Yes - the needs of an OpenMM implementation made life awkward for most algorithms, but those limitations are history because we are no longer working with OpenMM. > (5) free energy calculations > Berk has been working hard for the upcoming 5.0 release, which should fill in most or all of the blanks here - patch is in Gerrit awaiting us organizing some test cases. I don't actually do such calculations, but if you want to see if your needs are or might be met, please try out the code from https://gerrit.gromacs.org/#/c/2803/ right away! Do you know if there have been any changes in support for any of (1)-(5) > above between versions 4.5 and 4.6? > The only thing the 4.5 and 4.6 implementations have in common is running on GPUs and being called GROMACS. But, again, I know that the developers are working hard and that it must be > an extreme amount of work to implement things on the GPU. I'm not > complaining! > Sure. But if you're in the position to see if a test calculation does something that seems correct, this is exactly the feedback we need during the lead up to a release! Cheers, Mark -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
Re: [gmx-users] Supported features in Gromacs on GPU
Hi Mark, Thanks so much for your time in answering my questions. It helped me a lot. Andrew DeYoung Carnegie Mellon University These are all standard in 4.6. The GPU implementation there only does non-bonded pair interactions, which affects none of the above. The limitation in 4.5 was that the GPU did lots more parts of the computation, and not all kinds of algorithms are amenable to that. Mark -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] How can I concatenate multiple g96 configurations into a single trajectory?
Hi, I have been creating several .g96 configuration files "by hand" (actually, using a scripting language). I like .g96 as a configuration file format because it provides 9 decimal places of precision for the x, y, and z coordinates of the position and velocity of atoms in my system -- whereas .gro and .pdb configuration file formats offer much lower precision. My system is crystalline, and I found many months ago that the lower precision offered by .gro and .pdb is not quite enough to achieve the exact bond lengths that I want in my crystal. Many months ago, Professor van der Spoel suggested on this list that I switch to .g96 for my configurations. So for many months I have been successfully using .g96 configurations as input to grompp to make the .tpr files necessary for simulation. Today, however, I am trying to do something somewhat unusual, and I run into an unusual problem. Suppose I have made three .g96 configurations "by hand": conf1.g96, conf2.g96, and conf3.g96. I want to compute the single-point energy of all three configurations. (Why would I want to do this? I have written code in a scripting language to generate frames in which a particular angle is "scanned" from one extreme to the other; I want to use Gromacs to simply compute the single-point energy at each frame in the "scan," without taking any MD moves.) One way to do this might be to run three different zero-step MD simulations, by passing the three configurations separately to grompp, obtaining three separate .tpr files, and passing those three .tpr files to mdrun. But this documentation page http://www.gromacs.org/Documentation/How-tos/Single-Point_Energy warns that zero-step MD simulations have "complications related to catering to possible restarts" and is not recommended. Doing it that way is also a little inconvenient because eventually I will have dozens or hundreds of separate "handmade" .g96 configurations for which to calculate the single-point energy, and running hundreds of zero-step simulations would probably require some sort of bash script. So, my thought is to concatenate the three (or more) .g96 configurations into a single trajectory, and then feed that trajectory to the -rerun switch of mdrun to compute the energy of each frame. My question is, how can I concatenate three (or more) .g96 configurations into a single trajectory of any type? I have tried: trjcat -f conf1.g96 conf2.g96 conf3.g96 -o output.trr But I get this error message: --- Program trjcat, VERSION 4.5.5 Source code file: trxio.c, line: 693 Fatal error: Reading trajectories in .g96 format is broken. Please use a different file format. For more information and tips for troubleshooting, please check the GROMACS website at http://www.gromacs.org/Documentation/Errors --- Reading trajectories in .g96 format may be broken, but reading CONFIGURATIONS in .g96 is not broken, as far as I know (I routinely pass g96 configurations to the -c switch of grompp). Is there any other way I can concatenate three (or more) .g96 configurations into a single trajectory of any type (.trr, .xtc, etc) -- WITHOUT going through .gro or .pdb (e.g., by converting from .g96 using editconf), which have lower precision? Thanks so much for your time! Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
Re: [gmx-users] How can I concatenate multiple g96 configurations into a single trajectory?
Hi Tsjerk, Thanks so much for your time. Unfortunately, when I tried to convert a .g96 configuration to a .trr trajectory, with this command: trjconv -f conf1.g96 -o conf1.trr I get this error message: --- Program trjconv_d, VERSION 4.5.5 Source code file: trxio.c, line: 693 Fatal error: Reading trajectories in .g96 format is broken. Please use a different file format. For more information and tips for troubleshooting, please check the GROMACS website at http://www.gromacs.org/Documentation/Errors --- I am running Gromacs 4.5.5. Do you know if there is another way to convert a .g96 configuration to a .trr trajectory, without using trjconv? Or do you know if a different version of Gromacs would not have this problem? Thanks for your time. Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] Saving velocities of one atom at every timestep
Hi, I have a rather large system (~10,000 atoms) -- at least, I consider it rather large since I'm using a cluster, not a supercomputer. I want to save the velocities of a single atom of interest at every timestep, so that I can use g_velacc to calculate the velocity autocorrelation function. I only care about the single atom of interest; I'm not interested in the coordinates, velocities, and forces of the other many thousands of atoms in my system. In other words, I have enough disk space to save the velocities of the atom of interest for every timestep for, say, several millions of steps. But I don't have enough space to save the velocities of all 10,000 atoms for every timestep for several millions of steps. I'm running Gromacs version 4.5.5. Setting nstxtcout = 1 allows one to save the coordinates of the groups specified in xtc-grps to the xtc trajectory. Is there any built-in way to save velocities -- rather than coordinates -- for only a specified group (a single atom, in this situation)? Or would it be possible to use the xtc to save the coordinates of the atom of interest at every timestep, and then somehow convert these coordinates to velocities so that I can use g_velacc? I'm using the leapfrog integrator (integrator = md). Thanks so much for your time! Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] Specifying the number of OpenMP threads that g_hbond uses
Hi, >From the web page http://www.gromacs.org/About_Gromacs/Release_Notes/Versions_4.6.x , it looks like parallelization for g_hbond (via OpenMP) was added in version 4.6.x. What is a typical command to run g_hbond in parallel? Since I mostly have experience with Gromacs 4.5.5 (which for mdrun has auto MPI threading as the default parallelization, which means I only have to specify the number of threads with -nt), I am not sure how to launch OpenMP for g_hbond in 4.6.x. When I just run g_hbond in 4.6.1 with no extra flags, the output says, "Frame loop parallelized with OpenMP using 8 threads." Is there any way to change the number of OpenMP threads that g_hbond uses? Thanks so much! Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] Embedded view of the users list
Hi, For the past several years and until a few weeks ago, the GMX-Users List webpage here ( http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List ) had an embedded "live"/"interactive" view of the list. Messages were organized in threads such that when I clicked on a thread, I could see all posts in the thread on one page. It also had a nice way to search old messages and sort them by date, etc. Could someone please confirm that the embedded "live"/"interactive" view of the list has been taken down from http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List ? It is fine if it has been taken down, since I can still browse through the archives here ( https://mailman-1.sys.kth.se/pipermail/gromacs.org_gmx-users/ ) and I can search the archives by using google with site:https://mailman-1.sys.kth.se as Mark suggested. But I just wanted to make sure that it is not just me and my browser that no longer sees the embedded "live"/"interactive" view of the list. Thank you very much! Andrew DeYoung Carnegie Mellon University -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.