<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>

<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body text="#000000" bgcolor="#ffffff">
Stefane,<br>
<br>
If you are asking for a comparison of simulation time required, then
that is not something that can be answered in the abstract. There are
even cases in which REMD will take more total simulation time to
converge the 300K ensemble than will N simulations all at 300K. The
efficiency of REMD will be highly dependent on the size of the
important energy barriers. REMD becomes more efficient than MD as the
barriers get higher. To take that to an extreme, if you wanted to know
the distribution of x-coordinate values of one specific water molecule
in a box of pure water, REMD with 100 replicas would be 100 times less
efficient than simple MD -- there are no energy barriers and so the 99
replicas at temperatures above 300K are a total waste of time. At the
other extreme, if your energy barrier is 30 Kcal/mol at 300K, then you
will essentially never cross it at 300K and REMD, if it works, will be
infinitely more "efficient" because the replica at 300K will
essentially never sample the other side of that barrier.<br>
<br>
There have been a few comparisons, but I don;t know of any general rule
that you can apply here.<br>
<br>
On 10/13/10 7:42 PM, XAvier Periole wrote:<br>
<pre>&gt;<i>
</i>&gt;<i> Well the exchanges do not cost time ... so basically each replica runs
</i>&gt;<i> at the speed it would run it was a regular simulation on the number of
</i>&gt;<i> cpu you give each replica.
</i>&gt;<i>
</i>&gt;<i> I do not think there any significant difference in NPT/NVT simulations.
</i>&gt;<i> Note NVT is often used for REMD
</i>It us, but that gives high pressure at high T, which is typically not 
what you want.

If you use NPT REMD than the run time between exchanges is in general 
equal to that of the coldest replica, since that will be at the highest 
density (unless you go below the density maximum of water :)).


&gt;<i>
</i>&gt;<i> XAvier.
</i>&gt;<i>
</i>&gt;<i> On Oct 13, 2010, at 11:36 AM, ABEL Stephane 175950 wrote:
</i>&gt;<i>
</i>&gt;&gt;<i> Dear All,
</i>&gt;&gt;<i>
</i>&gt;&gt;<i> For a futur project, I would like to perform the REMD calculations
</i>&gt;&gt;<i> with GROMACS4.5.X. To have an estimation of CPU time required, I have
</i>&gt;&gt;<i> a naive question: What is the speed of REMD compared to a classic NPT
</i>&gt;&gt;<i> MD ? I am aware that the response depends a lot of factor (number of
</i>&gt;&gt;<i> replica, etc.), but since i have no idea and experience with this type
</i>&gt;&gt;<i> of MD, I ask the question ;).
</i>&gt;&gt;<i>
</i>&gt;&gt;<i> The system will have approximatively 25000 atoms. Previously for this
</i>&gt;&gt;<i> system, I got 11 ns/day on 32 CPU for a NPT MD with GROMOS53A6 and
</i>&gt;&gt;<i> GMX4.0.5.
</i>&gt;&gt;<i>
</i>&gt;&gt;<i> Thanks you in advance for your response.
</i>&gt;&gt;<i>
</i>&gt;&gt;<i> Stefane
</i>&gt;&gt;<i>
</i>&gt;&gt;<i> --
</i>&gt;&gt;<i> gmx-users mailing list <a
 href="http://lists.gromacs.org/mailman/listinfo/gmx-users">gmx-users at gromacs.org</a>
</i>&gt;&gt;<i> <a
 href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a>
</i>&gt;&gt;<i> Please search the archive at
</i>&gt;&gt;<i> <a
 href="http://www.gromacs.org/Support/Mailing_Lists/Search">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!
</i>&gt;&gt;<i> Please don't post (un)subscribe requests to the list. Use the
</i>&gt;&gt;<i> www interface or send it to <a
 href="http://lists.gromacs.org/mailman/listinfo/gmx-users">gmx-users-request at gromacs.org.</a>
</i>&gt;&gt;<i> Can't post? Read <a
 href="http://www.gromacs.org/Support/Mailing_Lists">http://www.gromacs.org/Support/Mailing_Lists</a>
</i>&gt;<i>
</i>

-- 
David.
</pre>
</body>
</html>