Thanks JJ and Mark. I shall revise my algorithm/technique. This was just a lazy way of implementing a more robust technique and I shall use the more rigorous form. <br>However, if one wants to use the current technique with a smaller temperature range, the dynamics of the system with say Protein at 500K and water at 300K would be unphysical or can that be corrected by using a smaller timestep and thermostat type/parameters?<br>
<br>Pooja<br><br><div class="gmail_quote">On Sun, Oct 31, 2010 at 10:15 PM, <span dir="ltr"><<a href="mailto:thompsjj@purdue.edu">thompsjj@purdue.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
I concur with Mark - you're pumping heat into the solvent which is physically unsound, particularly at the solute-solvent boundary in your situation. This would be analogous to making a block of metal at 900K instantly "appear" in water at 300K: what would happen?<br>
Technically the Nose-Hoover thermostat is going to be better than Berendsen at reducing the noncanonical bias in the replica exchange, but of course using a smaller range of temperatures is going to reduce the pressure on the thermostat to keep up with velocity scaling. This is a different problem entirely, though. To be honest, I see the difference in your management of protein and solvent temperatures as the biggest problem here.<div>
<div></div><div class="h5"><br>
<br>
<br>
Quoting Mark Abraham <<a href="mailto:mark.abraham@anu.edu.au" target="_blank">mark.abraham@anu.edu.au</a>>:<br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
<br>
----- Original Message -----<br>
From: Sai Pooja <<a href="mailto:saipooja@gmail.com" target="_blank">saipooja@gmail.com</a>><br>
Date: Monday, November 1, 2010 10:24<br>
Subject: Re: [gmx-users] LINCS vs SHAKE<br>
To: Discussion list for GROMACS users <<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Some information:<br>
Except for this note, I do not get any other warning when I generate the .tpr files. The cut-offs are in line with the forcefield recommendations and<br>
The largest charge group contains 12 atoms.<br>
Since atoms only see each other when the centers of geometry of the charge<br>
groups they belong to are within the cut-off distance, too large charge<br>
groups can lead to serious cut-off artifacts.<br>
For efficiency and accuracy, charge group should consist of a few atoms.<br>
For all-atom force fields use: CH3, CH2, CH, NH2, NH, OH, CO2, CO, etc<br>
</blockquote>
<br>
Probably not critical, but IIRC the latest GROMACS imposes single-atom charge groups for CHARMM, in line with CHARMM developer advice.<br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
@JJ<br>
<br>
I am currently using Nose-hoover and it does happen close to a replica-exchange. Since I am attempting an algorithm which only scales the temperature of the protein and not the solvent across replicas, your pointer on the thermostat seems very pertinent.<br>
</blockquote>
<br>
Yep. That's a prime source for problems - IIRC post-exchange velocities are rescaled globally, not group-wise. Even if the latter were implemented, you'd still have to reconcile the fact that your simulations constantly have heat flowing from the solute heat bath to the solvent heat bath via your system. Is that sound?<br>
<br>
Also, I'm not sure which temperature GROMACS is taking as the replica temperature, whether that is dependent on the order of the temperature-coupling groups in the .mdp file, and whether it shouldn't be regarding your situation as questionable/illegal and complaining at the grompp stage.<br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Do you think a weaker coupling algo like Berendsen might help resolve this problem?<br>
<br>
I can also try using a smaller temperature range. I am currently using 300-900K (with solvent set at 300 for all 6 replicas). Perhaps if I use a smaller temperature range like 300-550, this could be prevented?<br>
</blockquote>
<br>
If that's the cause, then the answer is "probably", but you should address the underlying problem/question as above.<br>
<br>
Mark<br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
Pooja<br>
<br>
On Sun, Oct 31, 2010 at 7:13 PM, <<a href="mailto:thompsjj@purdue.edu" target="_blank">thompsjj@purdue.edu</a>> wrote:<br>
Have you looked into changing the thermostat yet? Can you determine if the malformed residues are a result of a recent replica exchange? This continues to sound like an overheat that the thermostat cannot handle. ><br>
</blockquote>
><br>
><br>
><br>
><br>
><br>
><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Quoting Mark Abraham <<a href="mailto:mark.abraham@anu.edu.au" target="_blank">mark.abraham@anu.edu.au</a>>:<br>
</blockquote>
><br>
><br>
><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
----- Original Message -----<br>
From: Sai Pooja <<a href="mailto:saipooja@gmail.com" target="_blank">saipooja@gmail.com</a>><br>
Date: Monday, November 1, 2010 8:12<br>
Subject: Re: [gmx-users] LINCS vs SHAKE<br>
To: Discussion list for GROMACS users <<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>
</blockquote>
><br>
> Yes, I understand that. But then this is a protein in water simulation with CHarmm forcefield and I cannot find a reason as to why does this happen(esp after variable number of steps)?<br>
><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
MD integrates equations of motion. If the starting configuration, or integration procedure is not well-formed, then atoms will experience large forces and motions, which lead to more, which eventually lead to a smoking pile of debris. See <a href="http://www.gromacs.org/Documentation/Terminology/Blowing_Up" target="_blank">http://www.gromacs.org/Documentation/Terminology/Blowing_Up</a> for discussion and suggestions<br>
</blockquote>
><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Mark<br>
</blockquote>
><br>
> I have checked the .pdb files generated for the steps where Lincs fails and it seems that H atoms for atleast 1 sidechain are in incorrect positions(usually clashing). There doesn't seem to be a good way of analyzing this error. If you have any suggestions I would be eager to try them.<br>
> ><br>
> Pooja<br>
><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
On Sat, Oct 30, 2010 at 3:40 PM, Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>> wrote:<br>
</blockquote>
> ><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
><br>
</blockquote>
> Sai Pooja wrote:<br>
> ><br>
> For Replica exchange, is there any advantage of in using SHAKE over LINCS(other than the stepsize)?<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
I am running an REM simulation and the simulation stops after running for variable number of steps (100000, 10000000 etc.) because some bond moves more than 30 degrees and LINCS gives a warning.<br>
</blockquote>
> ><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
><br>
</blockquote>
> It has been said (check the archive) that LINCS is more stable. Just because you're getting LINCS warnings does not mean the constraint algorithm is to blame, it's just the first algorithm that fails when your model physics implodes.<br>
> ><br>
> -Justin<br>
> ><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
><br>
</blockquote>
> Pooja<br>
> ><br>
> --<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Quaerendo Invenietis-Seek and you shall discover.<br>
</blockquote>
> ><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
><br>
</blockquote>
> --<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
========================================<br>
</blockquote>
> ><br>
> Justin A. Lemkul<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
</blockquote>
> > <a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
><br>
</blockquote>
> ========================================<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
</blockquote>
> > <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
> Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</blockquote>
> ><br>
> ><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
--<br>
Quaerendo Invenietis-Seek and you shall discover.<br>
</blockquote>
> > --<br>
> gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
> <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Please search the archive at<br>
</blockquote>
> <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</blockquote>
><br>
><br>
><br>
><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
</blockquote>
> <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use thewww interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>. ><br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</blockquote>
><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
><br>
--<br>
Quaerendo Invenietis-Seek and you shall discover.<br>
</blockquote>
> --<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at<br>
<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</blockquote>
<br>
</blockquote>
<br>
<br>
<br>
-- <br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use thewww interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Quaerendo Invenietis-Seek and you shall discover.<br>