<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'><div dir='ltr'>
Thanks again for your quick answer!<br>No, I did not prior equilibration, however for the analysis I was considering only the last 0.5 ns. <br>I will try to extend the simulations and see what happens. Is it possible to do it using tpbconv as in a typical simulation? Should I create different files for pullf*.xvg and pullx*.xvg and then join them for the analysis?<br>Best wishes,<br>Rebeca.<br><br><div>&gt; Date: Tue, 12 Jul 2011 17:59:44 -0400<br>&gt; From: jalemkul@vt.edu<br>&gt; To: gmx-users@gromacs.org<br>&gt; Subject: Re: [gmx-users] number of windows in PMF<br>&gt; <br>&gt; <br>&gt; <br>&gt; Justin A. Lemkul wrote:<br>&gt; &gt; <br>&gt; &gt; <br>&gt; &gt; Rebeca García Fandiño wrote:<br>&gt; &gt;&gt; Thanks a lot for your answer.<br>&gt; &gt;&gt; My system have about 29000 atoms, and the simulation time was 1ns.<br>&gt; &gt;&gt; I have generated the error bars, and indeed you were right, they are <br>&gt; &gt;&gt; too big. The histograms looked good, so I thought they should be <br>&gt; &gt;&gt; well-converged...:S<br>&gt; &gt;&gt; What should I do, extending the simulation or generating more windows?<br>&gt; &gt; <br>&gt; &gt; The simulations probably need to be longer.  I've gotten good <br>&gt; &gt; convergence in much larger systems using far fewer windows but somewhat <br>&gt; &gt; longer simulations.<br>&gt; &gt; <br>&gt; <br>&gt; Another thing to consider (aside from my offhand guesses of simulation length) <br>&gt; is that if you're generating velocities at the beginning of each simulation in <br>&gt; each window in the absence of prior equilibration, you're not *really* getting 1 <br>&gt; ns of viable data.  Some should be discarded as equilibration, if you're not <br>&gt; doing that already.  If you've done prior equilibration in each window, then <br>&gt; disregard this thought.<br>&gt; <br>&gt; -Justin<br>&gt; <br>&gt; -- <br>&gt; ========================================<br>&gt; <br>&gt; Justin A. Lemkul<br>&gt; Ph.D. Candidate<br>&gt; ICTAS Doctoral Scholar<br>&gt; MILES-IGERT Trainee<br>&gt; Department of Biochemistry<br>&gt; Virginia Tech<br>&gt; Blacksburg, VA<br>&gt; jalemkul[at]vt.edu | (540) 231-9080<br>&gt; http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin<br>&gt; <br>&gt; ========================================<br>&gt; -- <br>&gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting!<br>&gt; Please don't post (un)subscribe requests to the list. Use the <br>&gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; Can't post? Read http://www.gromacs.org/Support/Mailing_Lists<br></div>                                               </div></body>
</html>