<div dir="ltr"><div><div><div><br>Hi Berk,<br><br></div>About the dynamics, you are right. As we modified the force, the dynamics become incorrect. Of cause, we have some special way to correct it, but it is too complicate so we haven&#39;t do it yet. In fact, this is one of disadvantage of our method. However, we usually use our method to calculate the thermodynamics.<br><br></div>In fact, our method have already published for many years (Yang and Gao. JCP 2009: <span>214109</span><span></span>). This method is similar to the REST method, but replica exchange is not needed. This method is much more efficient than REST. We perform only one trajectory, but this trajectory includes the thermodynamics of all the replicas of REST. As we only realized this method in AMBER (Each time when the new edition of AMBER was released, we have to modified the codes once again), just a few people had used it. Now I want to transplant this method to Gromacs. I think it will help us to make this method more popular. <br><br></div><div>Thank you and best regards,<br></div><div>Isaac<br></div><div><br></div><div><div class="gmail_extra"><br><div class="gmail_quote">2016-09-07 14:36 GMT+02:00 David van der Spoel <span dir="ltr">&lt;<a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a>&gt;</span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On 07/09/16 14:15, Yi Isaac Yang wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Berk,<br>
<br>
I am sorry to send so many same mails. Because this is the first time to<br>
use the developer list. I didn&#39;t receive my mail after I sent it, so I<br>
send it again and again.<br>
<br>
Thank you for reply! Our method is similar as REST (Replica exchange<br>
with solute tempering) which just to enhance the solute part of the<br>
system but keep the rest of the system as usual. However, our method<br>
can&#39;t be realized by modified the force field file like in REST<br>
(Terakawa et al. JCC 2010: 1228). We must read the intra-energy of<br>
solute (marked as U_slu-slu) and the interaction energy between solute<br>
and solvent (marked as U_slu-sol) and to modified them. In order to put<br>
the modified potential energy into the system, we need to know the<br>
component forces caused by the component energy F_slu-slu=-(\partial<br>
U_slu-slu)/(\partial r) and F_slu-sol=-(\partial U_slu-sol)/(\partial r).<br>
<br>
Of cause, I know we cannot divide the whole potential energy into<br>
different parts (solute, solvent or interaction energy). Even in AMBER,<br>
we just draw the short range interaction of the different parts.<br>
However, even if we can only read a part of component energy, at least<br>
we can still enhance the sampling of this part.<br>
</blockquote></span>
I think the problem is ill defined and if you just do the short-range parts you are missing all the reaction field due to the solvent at long range. By selecting just part of the forces and doing something funny with it your dynamics may also be incorrect (that is if you do not strictly have F_i = dU/dr_i your dynamics becomes incorrect which you probably need to compensate for with a strong T-coupling).<br>
<br>
With REST there is a similar problem that only one of the replica&#39;s is physical, the rest are not. However, even the ensemble at the &quot;physical replica&quot; in that case may be skewed by the ensembles at different temperatures, especially for short  exchange times I would think.<br>
<br>
Does your method work with replica&#39;s as well? Can you ascertain that at least one replica is completely physical? If not I would advise you to reconsider your project. Sorry for the harsh words, but if the physics is not consistent the method is not useful.<br>
<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">
<br>
Thank you and best regrads,<br>
Isaac<br>
<br></span>
2016-09-07 12:17 GMT+02:00 Berk Hess &lt;<a href="mailto:hess@kth.se" target="_blank">hess@kth.se</a> &lt;mailto:<a href="mailto:hess@kth.se" target="_blank">hess@kth.se</a>&gt;&gt;:<div><div class="h5"><br>
<br>
    Hi,<br>
<br>
    Please don&#39;t send the same mail multiple times.<br>
<br>
    Your request is not clear enough. Do you want forces on different<br>
    parts of the system or between different parts of the system? The<br>
    latter is difficult to obtain.<br>
    Note that with PME, which is what you should be using in most cases,<br>
    it is difficult, if not impossible, to decompose forces and energies<br>
    between group (pairs).<br>
<br>
    Cheers,<br>
<br>
    Berk<br>
<br>
<br>
    On 2016-09-07 11:46, Yi Isaac Yang wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
    Dear developers,<br>
<br>
    I am developing an enhanced sampling method, and I have already<br>
    realized it in AMBER. Now I want to add the code of this method<br>
    into Gromacs. This method need to read the component energy and<br>
    force of different part of the system (like protein, solution,<br>
    etc.) during the simulation. However, I am not very familiar withe<br>
    the codes of Gromacs, so I don&#39;t know how to read it. I know when<br>
    I perform MD simulation using Gromacs I can set the groups in<br>
    &quot;energygrps&quot;, then I can read the energy between those groups in<br>
    .edr file. So I think it shouldn&#39;t be too difficult to read the<br>
    component energy and force in the codes. Can you tell me how to<br>
    read the component energy and force from the codes?<br>
<br>
    Thank you very much and best regards,<br>
    Isaac<br>
<br>
<br>
</blockquote>
<br>
<br>
    --<br>
    Gromacs Developers mailing list<br>
<br>
    * Please search the archive at<br>
    <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List" rel="noreferrer" target="_blank">http://www.gromacs.org/Support<wbr>/Mailing_Lists/GMX-developers_<wbr>List</a><br>
    &lt;<a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List" rel="noreferrer" target="_blank">http://www.gromacs.org/Suppor<wbr>t/Mailing_Lists/GMX-developers<wbr>_List</a>&gt;<br>
    before posting!<br>
<br>
    * Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" rel="noreferrer" target="_blank">http://www.gromacs.org/Support<wbr>/Mailing_Lists</a><br>
    &lt;<a href="http://www.gromacs.org/Support/Mailing_Lists" rel="noreferrer" target="_blank">http://www.gromacs.org/Suppor<wbr>t/Mailing_Lists</a>&gt;<br>
<br>
    * For (un)subscribe requests visit<br>
    <a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers" rel="noreferrer" target="_blank">https://maillist.sys.kth.se/ma<wbr>ilman/listinfo/gromacs.org_gmx<wbr>-developers</a><br>
    &lt;<a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers" rel="noreferrer" target="_blank">https://maillist.sys.kth.se/m<wbr>ailman/listinfo/gromacs.org_gm<wbr>x-developers</a>&gt;<br>
    or send a mail to <a href="mailto:gmx-developers-request@gromacs.org" target="_blank">gmx-developers-request@gromacs<wbr>.org</a><br>
    &lt;mailto:<a href="mailto:gmx-developers-request@gromacs.org" target="_blank">gmx-developers-request<wbr>@gromacs.org</a>&gt;.<br>
<br>
<br>
<br>
<br>
--<br>
Yesterday Young<br>
College of Chemistry and Molecular Engineering<br>
Peking University<br>
<br>
<br>
</div></div></blockquote><span class="HOEnZb"><font color="#888888">
<br>
<br>
-- <br>
David van der Spoel, Ph.D., Professor of Biology<br>
Dept. of Cell &amp; Molec. Biol., Uppsala University.<br>
Box 596, 75124 Uppsala, Sweden. Phone:  <a href="tel:%2B46184714205" value="+46184714205" target="_blank">+46184714205</a>.<br>
<a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a>    <a href="http://folding.bmc.uu.se" rel="noreferrer" target="_blank">http://folding.bmc.uu.se</a></font></span><div class="HOEnZb"><div class="h5"><br>
-- <br>
Gromacs Developers mailing list<br>
<br>
* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List" rel="noreferrer" target="_blank">http://www.gromacs.org/Support<wbr>/Mailing_Lists/GMX-developers_<wbr>List</a> before posting!<br>
<br>
* Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" rel="noreferrer" target="_blank">http://www.gromacs.org/Support<wbr>/Mailing_Lists</a><br>
<br>
* For (un)subscribe requests visit<br>
<a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers" rel="noreferrer" target="_blank">https://maillist.sys.kth.se/ma<wbr>ilman/listinfo/gromacs.org_gmx<wbr>-developers</a> or send a mail to <a href="mailto:gmx-developers-request@gromacs.org" target="_blank">gmx-developers-request@gromacs<wbr>.org</a>.<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">Yesterday Young<br>College of Chemistry and Molecular Engineering<div>Peking University</div></div>
</div></div></div>