Hi,<div><br></div><div>Yes, use the configuration from geometry optimization to start MD. Though in your original post, the thermostat coupling constant is much greater than the length of the trajectory you ran, whether massive or global - it's not doing anything here. The velocities are randomized and scaled to the target temperature by default as well so you should expect differences in the initial temperature fluctuations. After running for some time and reaching equilibrium, the average temperature from both thermostat settings should be equal. The massive thermostat controls every degree of freedom so it should get to the target temperature in fewer MD steps, making it ideal for equilibration. You should exercise caution in using MASSIVE for a given thermostat and coupling constant when you are interested in dynamic properties during data collection (aka production) runs. See, J. Chem. Theory Comput. 2013, 9, 2887−2899.</div><div><br></div><div>-T<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, March 26, 2021 at 10:15:18 PM UTC-4 tang...@gmail.com 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 dir="ltr">Hello Julian,<div> Thanks for your suggestion. I will do geometry optimization immediately. However, I'm a little confused about the last sentence in your reply, whether I should use the structure gotten by geometry optimization to run massive and global thermostat firstly and than run others processing? I hope i have expressed my question clearly. </div><div> Tangyu.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Julian Heske <<a href data-email-masked rel="nofollow">j...@gmail.com</a>> 于2021年3月26日周五 下午11:03写道：<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"></blockquote></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>Hello,</div><div><br></div><div>in order to run a MD calculation at a specific tempertare, you have to equilibrate the system first to reach your target temperature. The time you need to equilibrate your system is dependent on your initial configuration. Since you are starting with a structure which has a very high potential energy (-157,5 vs. ~ -168 a. u. of the others, see out.txt), this energy remains as kinetic energy in the system and the thermostat will need a long time to remove the energy from the system. You can first do a geometry optimization ( <span><a href="https://manual.cp2k.org/trunk/CP2K_INPUT.html" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://manual.cp2k.org/trunk/CP2K_INPUT.html&source=gmail&ust=1616907924017000&usg=AFQjCNHJ9fDPFOplXQujOfqaNKNsSjaJWQ">CP2K_INPUT</a> /
<a href="https://manual.cp2k.org/trunk/CP2K_INPUT/MOTION.html" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://manual.cp2k.org/trunk/CP2K_INPUT/MOTION.html&source=gmail&ust=1616907924017000&usg=AFQjCNGpnLIs0A25m5mOb9pGLYm8AYpIoA">MOTION</a> /
<a href="https://manual.cp2k.org/trunk/CP2K_INPUT/MOTION/GEO_OPT.html" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://manual.cp2k.org/trunk/CP2K_INPUT/MOTION/GEO_OPT.html&source=gmail&ust=1616907924017000&usg=AFQjCNHpS_LNhwyZnhFoY05vEJ8nXXTeoA">GEO_OPT</a> )</span>, which leads to an optimized structure. Then start your MD run with your obtained structure as initial configuration. However, you should then still take a look at the temperature and not use the equilibration steps for your analysis.<br></div><div><br></div><div>Kind regards,</div><div>Julian<br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr"><a href data-email-masked rel="nofollow">tang...@gmail.com</a> schrieb am Freitag, 26. März 2021 um 15:16:33 UTC+1:<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>Dears cp2k's users,</div><div> I'm doing a vibration analysis exercise.When I did massive and global thermostat, the temperature I get is different from what I set. Could you tell me what causes this error? I would appreciate you generous help. The input and output file are listed as follow.</div></blockquote></div>
<p></p></blockquote></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
-- <br>
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.<br>
To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/MnB_Ejspnpw/unsubscribe" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/topic/cp2k/MnB_Ejspnpw/unsubscribe&source=gmail&ust=1616907924017000&usg=AFQjCNFbaYas0q90tEdLCmSRidatp1oe7w">https://groups.google.com/d/topic/cp2k/MnB_Ejspnpw/unsubscribe</a>.<br>
To unsubscribe from this group and all its topics, 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/ac99e428-7df6-4e03-8b6e-cbf7c16009dbn%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/ac99e428-7df6-4e03-8b6e-cbf7c16009dbn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1616907924017000&usg=AFQjCNGpXFoEC12f0zrsJMau31_c8VHzTQ">https://groups.google.com/d/msgid/cp2k/ac99e428-7df6-4e03-8b6e-cbf7c16009dbn%40googlegroups.com</a>.<br>
</blockquote></div>
</blockquote></div>