<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">CP2K_INPUT</a> /
               <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/MOTION.html">MOTION</a> /
               <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/MOTION/GEO_OPT.html">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">tang...@gmail.com schrieb am Freitag, 26. März 2021 um 15:16:33 UTC+1:<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>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>