Dear Lenard, dear Thomas, <div><br></div><div>I am still of the idea that defining these as NVT or NVE simulations is misleading, </div><div>since the generated configurations do not belong to any of these ensembles. </div><div>Moreover, adding some kinetic energy stepwise and running 10/20 ps to equilibrate end explore, at least partially, the corresponding ensemble, is for sure much better than increasing the temperature from 0 to 800K in 2 ps, as was the original intention. </div><div>Best</div><div>Marcella</div><div><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, April 1, 2021 at 12:29:54 PM UTC+2 tkuehne wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div style="word-wrap:break-word;line-break:after-white-space">Dear Lenard, <div><br></div><div>alternatively you may try using TEMEPRATURE_ANNEALING with values slightly larger than 1. </div><div>However, this corresponds to an exponential not linear increase in temperature. </div><div><br></div><div>Cheers, </div><div>Thomas </div><div><br><div><blockquote type="cite"></blockquote></div></div></div><div style="word-wrap:break-word;line-break:after-white-space"><div><div><blockquote type="cite"><div>Am 29.03.2021 um 18:45 schrieb Lenard Carroll <<a href data-email-masked rel="nofollow">len...@gmail.com</a>>:</div><br></blockquote></div></div></div><div style="word-wrap:break-word;line-break:after-white-space"><div><div><blockquote type="cite"><div><div dir="ltr">But adding the kinetic energy all at once would result in the system acting non-physically. I know this as I've done classical MD simulations for years. I like the approach of multiple NVT simulations, gradually increasing the temperature. I think I will also be using the previous temperature's final trajectory/velocities for the new temperature NVT simulation.<br><br>Thanks</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 29, 2021 at 4:20 PM Marcella Iannuzzi <<a href data-email-masked rel="nofollow">ma...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Lenard, <div><br></div><div>Heating the system by constantly rescaling the velocity generates a dynamics that does not correspond to any meaningful thermodynamic ensemble.</div><div>For sure, to set up a thermostat while forcing the temperature to change is not a good idea.</div><div>If you want to use the ANNEALING in cp2k just set NVE, obviously the energy is not going to be conserved.</div><div>Activating the annealing, at each step the actual temperature is multiplied by the scaling factor and the velocities are rescaled accordingly.</div><div><br></div><div>Probably a better way to bring the system from a lower to a higher temperature would be to run a set of NVT simulations at increasing temperature, say 150K, 300K, ...800K</div><div>Each simulation should be several ps long, such that the system has the time to sort of equilibrate at each temperature before ramping up.</div><div><br></div><div>best</div><div>Marcella</div><br><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, March 29, 2021 at 4:19:17 PM UTC+2 <a href data-email-masked rel="nofollow">jos...@gmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="auto">There is no way to do 800K over 2 ps.</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 29, 2021, 15:03 ASSIDUO Network <<a rel="nofollow">len...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Good day, I trust you are doing well.<div><br></div><div>I'm in the process of getting my hands dirty with AIMD using CP2K. I am familiar with MD simulations from AMBER, but not quite yet CP2K.<br><br>I would like set up a gentle heating in my AIMD simulation, gently heating the box up from 0 K to 800 K. I have read on this group that annealing should be used for this, but my issues here are the following:<br>- The heating phase has the volume and number of molecules constant, but allows the energy, the temperature and the pressure to change. None of the ensembles seem to fit this. So, I decided to go with NVT since that is what I have used with the heating simulation with AMBER (many tutorials available on it), but there you could specify the initial and reference temperature.</div><div>- How do I properly use annealing to heat the temperature up from 0 K to 800 K over 2 ps? Annealing just scales the velocities of the particles, so I'm not sure how to best choose this scaling factor.</div><div>- How do I print the temperature throughout the AIMD run?</div><div><br></div><div>I have included my input file. I would appreciate if you could go through the &MOTION section and tell me if it all makes sense.<br><br>Kind regards,<br>Lenard</div><div><br></div>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="noreferrer nofollow">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/e2da3b1a-39d8-45aa-912f-531eb540f12en%40googlegroups.com?utm_medium=email&utm_source=footer" rel="noreferrer nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/e2da3b1a-39d8-45aa-912f-531eb540f12en%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1617365901597000&usg=AFQjCNHKzp9h1rPQxFZGlI4x3aL9cS1FaQ">https://groups.google.com/d/msgid/cp2k/e2da3b1a-39d8-45aa-912f-531eb540f12en%40googlegroups.com</a>.<br>
</blockquote></div>
</blockquote></div><div><br></div>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/5f2e3349-9440-4bde-b87f-cdde5a90b9d3n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/5f2e3349-9440-4bde-b87f-cdde5a90b9d3n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1617365901597000&usg=AFQjCNGaG3L3QAa5amF4vjlbeRPAz4jeWw">https://groups.google.com/d/msgid/cp2k/5f2e3349-9440-4bde-b87f-cdde5a90b9d3n%40googlegroups.com</a>.<br>
</blockquote></div><div><br></div>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>.<br></div></blockquote></div></div></div><div style="word-wrap:break-word;line-break:after-white-space"><div><div><blockquote type="cite"><div>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAKur15-3ngs%3D1_ZxjM9FVNvmSj_Q9oBVwG0OYF%2BV%2BX53oru9BQ%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/CAKur15-3ngs%253D1_ZxjM9FVNvmSj_Q9oBVwG0OYF%252BV%252BX53oru9BQ%2540mail.gmail.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1617365901597000&usg=AFQjCNFQhZuVT4zvJe-AYhddwihZsBhWbg">https://groups.google.com/d/msgid/cp2k/CAKur15-3ngs%3D1_ZxjM9FVNvmSj_Q9oBVwG0OYF%2BV%2BX53oru9BQ%40mail.gmail.com</a>.<br>
</div></blockquote></div><br></div></div></blockquote></div>