<p class="MsoNormal" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: tahoma, sans-serif; border-collapse: collapse; font-size: 13px; "><font size="2">Forwarding this email from my group colleague:</font></p>

<p class="MsoNormal" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: tahoma, sans-serif; border-collapse: collapse; font-size: 13px; "><font size="2"><br></font></p><p class="MsoNormal" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: tahoma, sans-serif; border-collapse: collapse; font-size: 13px; ">

<font size="2">Dear Gromacs users,</font></p><p class="MsoNormal" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: tahoma, sans-serif; border-collapse: collapse; font-size: 13px; ">

<font size="2"> </font></p><p class="MsoNormal" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: tahoma, sans-serif; border-collapse: collapse; font-size: 13px; "><font size="2">I am trying to simulate a cellulose fiber in an ionic liquid solution in the NPT ensemble.<span>  </span>During the simulation, the entire system is coupled to a thermostat.<span>  </span>Yet, I observe an inhomogeneous temperature distribution throughout my system (hot-solvent/cold-solute) when I use Parrinello-Rahman pressure coupling but NOT when I employ Berendsen pressure coupling.<span>  </span>I have tested velocity-rescaling and the Nose-Hoover scheme to keep the temperature constant and in both cases Parrinello-Rahman pressure coupling seems to cause the solute’s temperature to become significantly lower than the solvent’s (to decompose temperatures, I am using “mdrun -rerun” with a run input that defines tc_grps separately).</font></p>

<p class="MsoNormal" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: tahoma, sans-serif; border-collapse: collapse; font-size: 13px; "><font size="2"> </font></p><p class="MsoNormal" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: tahoma, sans-serif; border-collapse: collapse; font-size: 13px; ">

<font size="2">I was wondering whether there were any known algorithmic reasons for this unphysical temperature gradient when using Parrinello-Rahman pressure coupling.</font></p><p class="MsoNormal" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: tahoma, sans-serif; border-collapse: collapse; font-size: 13px; ">

<font size="2">Thank you.</font></p><p class="MsoNormal" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: tahoma, sans-serif; border-collapse: collapse; font-size: 13px; "><font size="2">Barmak</font></p>

<p class="MsoNormal" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-family: tahoma, sans-serif; border-collapse: collapse; font-size: 13px; "><font size="2"><br></font></p><p class="MsoNormal" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">

<font size="2"><font class="Apple-style-span" face="tahoma, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">Comment from me: The effect is large. The ionic liquid is 5 degrees higher and the cellulose is 50 degrees lower (after 50ps, after that it stays constant). With Berendsen pressure both parts fluctuate around the same target temperature (as one would expect). Any reason why one doesn&#39;t get the correct temperature with rerun? Or is their a better way to get the temperature for different groups(for a simulation with just one tc-group)? Any reason why Parrinello-Rahman pressure coupling would have this effect on the temperature?</span></font></font></p>

<p class="MsoNormal" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font size="2"><font class="Apple-style-span" face="tahoma, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br>

</span></font></font></p><p class="MsoNormal" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font size="2"><font class="Apple-style-span" face="tahoma, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">Roland</span></font></font></p>

<br>-- <br>ORNL/UT Center for Molecular Biophysics <a href="http://cmb.ornl.gov">cmb.ornl.gov</a><br>865-241-1537, ORNL PO BOX 2008 MS6309<br>