[gmx-users] REMD: A small bug in repl_ex.c and a related question

2014-06-26 Thread Suman Chakrabarty
Hello! 1. It seems I have encountered a minor bug (?) in repl_ex.c for version 4.6.x (for NPT simulations only): Line 880 (in version 4.6.5): fprintf(fplog, dpV = %10.3e d = %10.3e\nb, dpV, delta + dpV); should be: fprintf(fplog, dpV = %10.3e d = %10.3e\n, dpV,

Re: [gmx-users] REMD - step 12005: Water molecule starting at atom 7030 can not be settled

2014-05-14 Thread Mark Abraham
OK, but with Berendsen you're also not observing physics. Your call. Mark On Wed, May 14, 2014 at 3:50 AM, Steven Neumann s.neuman...@gmail.comwrote: Thanks a lot! Berendsen barostat works well with V-rescale using REMD, no crashes observed. On Tue, May 13, 2014 at 2:08 PM, Mark Abraham

Re: [gmx-users] REMD - step 12005: Water molecule starting at atom 7030 can not be settled

2014-05-13 Thread Mark Abraham
On May 13, 2014 7:21 AM, Steven Neumann s.neuman...@gmail.com wrote: Thanks a lot! So which combination of thermostat/barostat would provide a proper ensemble without crashes? Various - but only P-R for the barostat in GROMACS, so you must be conservative with REMD exchange period. That will

[gmx-users] REMD - step 12005: Water molecule starting at atom 7030 can not be settled

2014-05-12 Thread Steven Neumann
Dear Gmx Users, I am running REMD of a small protein in water. I have 36 temperatures. After minimization I equilibrate each with nvt and npt ensemble. Then I start REMD with -multidir. The mdp: ; Run parameters integrator = md; leap-frog integrator nsteps = 5000 dt =

Re: [gmx-users] REMD - step 12005: Water molecule starting at atom 7030 can not be settled

2014-05-12 Thread Steven Neumann
What would you recommend then? Equilibrate with Berendsen barostat and v-rescale thermostat? On Mon, May 12, 2014 at 2:01 PM, Justin Lemkul jalem...@vt.edu wrote: On 5/12/14, 3:58 AM, Steven Neumann wrote: Dear Gmx Users, I am running REMD of a small protein in water. I have 36

Re: [gmx-users] REMD - step 12005: Water molecule starting at atom 7030 can not be settled

2014-05-12 Thread Mark Abraham
I would generally encourage the use of NVT REMD. However, note that P-R does exhibit unstable oscillatory behaviour when not close to equilibrium. REMD is supposed to stay in equilibrium by design, but this gets double crossed if the exchange period is too short with respect to the autocorrelation

Re: [gmx-users] REMD - step 12005: Water molecule starting at atom 7030 can not be settled

2014-05-12 Thread Steven Neumann
Thanks a lot! So which combination of thermostat/barostat would provide a proper ensemble without crashes? I am looking in the literature but did not find anything. On Tue, May 13, 2014 at 6:34 AM, Mark Abraham mark.j.abra...@gmail.comwrote: I would generally encourage the use of NVT REMD.

Re: [gmx-users] REMD, demux after 100 ns per replica

2014-05-02 Thread Christopher Neale
, Chris. From: gromacs.org_gmx-users-boun...@maillist.sys.kth.se gromacs.org_gmx-users-boun...@maillist.sys.kth.se on behalf of aperez a.pervi...@gmail.com Sent: 22 April 2014 12:46 To: gmx-us...@gromacs.org Subject: Re: [gmx-users] REMD, demux after 100 ns

Re: [gmx-users] REMD, demux after 100 ns per replica

2014-05-02 Thread Christopher Neale
. From: gromacs.org_gmx-users-boun...@maillist.sys.kth.se gromacs.org_gmx-users-boun...@maillist.sys.kth.se on behalf of Christopher Neale chris.ne...@alum.utoronto.ca Sent: 02 May 2014 14:28 To: gmx-us...@gromacs.org Subject: Re: [gmx-users] REMD, demux after 100 ns per replica Dear

Re: [gmx-users] REMD, demux after 100 ns per replica

2014-04-22 Thread aperez
Hi, I did an extra modification to the *.pl with the aim of being more general with the time step. I have just added a line that asks the time step. Note: I don't know perl, so I apologize for the rustic way how I modified the file. Thanks, Andrea proof_demux.pl

[gmx-users] REMD, demux after 100 ns per replica

2014-04-16 Thread aperez
Hi all, Today I have experienced a problem with the demux script applied in a replica exchange simulation. I have 16 replicas, and for each replica I have about 120 ns of MD. When I executed the script demux.pl on the log file, I noticed that the first column of the replica_*.xvg files had a

Re: [gmx-users] REMD, demux after 100 ns per replica

2014-04-16 Thread Mark Abraham
Hi, Yes, I have seen this. Unfortunately that script dates from the days when 10ns was a long run, and like most things in GROMACS, hasn't seen much love since one of the developers was actively using it. If you can test your fixed version and want to contribute it back, you can follow the

Re: [gmx-users] REMD, demux after 100 ns per replica

2014-04-16 Thread Christopher Neale
expand upon the bug report. Chris. From: gromacs.org_gmx-users-boun...@maillist.sys.kth.se gromacs.org_gmx-users-boun...@maillist.sys.kth.se on behalf of aperez a.pervi...@gmail.com Sent: 16 April 2014 11:48 To: gmx-us...@gromacs.org Subject: [gmx-users

Re: [gmx-users] REMD Melting Curve

2014-03-05 Thread Justin Lemkul
On 3/4/14, 11:26 PM, atanu das wrote: Thanks Justin. That's reassuring. However, I was wondering if any smoothing function is used in getting the melted fractions at intermediate temperatures with 1-K resolution. I could not find it in the description of the g_kinetics program. I need to know

Re: [gmx-users] REMD Melting Curve

2014-03-04 Thread atanu das
Dear Sir,                I successfully ran demix.pl and generated the xvg files - replica_ndx.xvg and replica_temp.xvg. I am attaching the file that I got as melt.xvg below .. # This file was created Wed Feb 26 10:10:27 2014 # by the following command: # g_kinetics_d -f replica_temp.xvg -d

Re: [gmx-users] REMD Melting Curve

2014-03-04 Thread Justin Lemkul
On 3/4/14, 12:22 PM, atanu das wrote: Dear Sir, I successfully ran demix.pl and generated the xvg files - replica_ndx.xvg and replica_temp.xvg. I am attaching the file that I got as melt.xvg below .. # This file was created Wed Feb 26 10:10:27 2014 # by the following

Re: [gmx-users] REMD Melting Curve

2014-03-03 Thread David van der Spoel
On 2014-03-04 00:20, atanu_das wrote: Dear Gromacs users, I have performed a REMD simulation with 12 replicas exponentially spaced between 260-420K. However, when I applied g_kinetics program to generate the melting curve, I got 161 values starting from 260K to 420K having the

Re: [gmx-users] REMD slow's down drastically

2014-02-25 Thread Mark Abraham
of Mark Abraham mark.j.abra...@gmail.com Sent: 24 February 2014 21:26 To: Discussion list for GROMACS users Subject: Re: [gmx-users] REMD slow's down drastically On Feb 24, 2014 11:01 PM, Christopher Neale chris.ne...@alum.utoronto.ca wrote: Presuming that you have indeed set up

Re: [gmx-users] REMD slow's down drastically

2014-02-24 Thread Dr. Vitaly Chaban
sorry, and what performance did you expect from REMD? Dr. Vitaly V. Chaban On Mon, Feb 24, 2014 at 8:32 AM, Singam Karthick sikar...@yahoo.in wrote: Dear members, I am trying to run REMD simulation for poly Alanine (12 residue) system. I used remd generator to get the range of temperature

Re: [gmx-users] REMD slow's down drastically

2014-02-24 Thread Singam Karthick
Dear Francis we are running in the  Xeon E5-2670 8C 2.60GHz (2 CPUs , 8 cores, 16 threads) for each temperature. and the exchange attempt frequency is 500 steps. The other system with 126 replicas run 30 ns per day ( system size of  ~38000 atoms). could you please help us in solving this

Re: [gmx-users] REMD slow's down drastically

2014-02-24 Thread Mark Abraham
Adding replicas cannot of itself slow things down, though it will increase the cost linearly. Don't try to run them all on the same amount of hardware as a smaller calculation! You are shooting yourself in the foot if you do not have at least one processor per replica (= MPI rank). Mark On Mon,

Re: [gmx-users] REMD slow's down drastically

2014-02-24 Thread Francis Jing
So what's the difference between the two systems? I think you can look through the end of the log file, and find which part slowed down. Maybe your job did not use the right partition of cores (too much comm. time)? I don't know. Francis On 25 Feb, 2014 1:19 am, Singam Karthick sikar...@yahoo.in

Re: [gmx-users] REMD slow's down drastically

2014-02-24 Thread Christopher Neale
. From: gromacs.org_gmx-users-boun...@maillist.sys.kth.se gromacs.org_gmx-users-boun...@maillist.sys.kth.se on behalf of Singam Karthick sikar...@yahoo.in Sent: 24 February 2014 02:32 To: gromacs.org_gmx-users@maillist.sys.kth.se Subject: [gmx-users] REMD slow's down drastically Dear members, I

Re: [gmx-users] REMD slow's down drastically

2014-02-24 Thread Christopher Neale
. From: gromacs.org_gmx-users-boun...@maillist.sys.kth.se gromacs.org_gmx-users-boun...@maillist.sys.kth.se on behalf of Mark Abraham mark.j.abra...@gmail.com Sent: 24 February 2014 21:26 To: Discussion list for GROMACS users Subject: Re: [gmx-users] REMD slow's

[gmx-users] REMD slow's down drastically

2014-02-23 Thread Singam Karthick
Dear members, I am trying to run REMD simulation for poly Alanine (12 residue) system. I used remd generator to get the range of temperature with the exchange probability of 0.3. I was getting the 125 replicas. I tried to simulate 125 replicas its drastically slow down the simulation time (for

Re: [gmx-users] REMD slow's down drastically

2014-02-23 Thread Francis Jing
Because every replica has to be calculated separately, it seems the speed is not that slow (125 * 0.07 ~ 9 ns). Also, to evaluate it, you should post how many CPUs you used, and other information like exchange attempt frequency... Francis On Mon, Feb 24, 2014 at 3:32 PM, Singam Karthick

[gmx-users] REMD

2013-12-03 Thread Shine A
Sir, I want to do an REMD simulation having 16 replicas. But I have only 8 processors. Is it possible to do 16 replicas in 8 processors? How I can do this? -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before

[gmx-users] REMD on GPU

2013-11-28 Thread James Starlight
Dear Gromacs users! I'd like to perform replica exchange simulation For this I've made bash script which create n folders like replica-298 replica-312 replica-323 replica-334 ... replica-N with all files needed for simulation considted of specified mdp file with different ref_t value No I'd

Re: [gmx-users] REMD on GPU

2013-11-28 Thread Mark Abraham
On Thu, Nov 28, 2013 at 3:01 PM, James Starlight jmsstarli...@gmail.comwrote: Dear Gromacs users! I'd like to perform replica exchange simulation For this I've made bash script which create n folders like replica-298 replica-312 replica-323 replica-334 ... replica-N with all files

<    1   2