<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dear Mengxu Li<div> </div><div>After reading your inputs, I think there are some parameters you should pay attention to.</div><div>1. REGION MASSIVE # too aggressive. global may be better for equilibrium.</div><div>2. NOSE LENGTH 5 # 3 or 4 may be suitable</div><div>3. REL_CUTOFF 30 # you should test. I suggest 50 or 60.</div><div>4. BASIS_SET DZVP-MOLOPT-SR-GTH # this is good for metals, which may not be suitable for h2o. Please test convergence such as DZVP ,TZVP,TZV2P, etc.</div><div>Lastly, I suggest that you can consider the ORCA program for nonperiodic systems.</div><div> </div><div> Best wishes,</div><div> QJ</div><div> <div class="gmail_quote"><div dir="ltr" class="gmail_attr">'Mengxu Li (Mengxu)' via cp2k <<a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a>> 于2022年11月27日周日 14:54写道： </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"> Dear all, I' m simulating the AIMD of the H+(H2O)3 cluster. I have seen this work, Schran et al., J. Chem. Theory Comput., 16, 88 (2020). The article explicitly mentioned used parameter settings when simulating AIMD in the article: the cutoff of 500Ry, TZV2P basis set, the time step of 0.25fs, temperature of 300K, Nosé chain thermostat, etc. We have used similar parameter settings and simulated 25ps, but the temperature fluctuation range is large and there is no trend of equilibrium . Could you give me some advice? I would appreciate your assistance in my problem. PS: the input and output files are attached.<img alt="temperature.png" width="540px" height="174px" src="cid:60f90d24-246b-458b-8b5e-e78771dc0725" title="" style=""> Yours sincerely, LMX -- You received this message because you are subscribed to the Google Groups "cp2k" group. To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/a338640d-10d7-4740-84f9-a68fdaedbf07n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/a338640d-10d7-4740-84f9-a68fdaedbf07n%40googlegroups.com</a>. </blockquote></div></div></div></div></div></div> -- You received this message because you are subscribed to the Google Groups "cp2k" group. To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAHhCLPHjUtepsG_EeXcZqRnCYWTWPxL3mpg4y35Gfp8QjmTB%2Bw%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAHhCLPHjUtepsG_EeXcZqRnCYWTWPxL3mpg4y35Gfp8QjmTB%2Bw%40mail.gmail.com</a>.