Thank you Lucas. <div><br></div><div>I've applied  geo-opt before AIMD, temperature still rised a lot. But now I know the temperature rising is because the initial geometry from your advice. I'll try to improve initial geo. And it is also a good advice to add D3 dispersion correction.</div><div><br></div><div>I have another question, when we run AIMD, the outer sphere of water molecules sometimes run too far away from our system, and actually, the box is in a right size. What should I do to limit the motion of water to make it not run too far.</div><div><br></div><div><br></div><div>Best,</div><div><br></div><div>Jessie</div><div><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, August 24, 2022 at 2:48:59 AM UTC+8 Lucas Lodeiro 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">As a first glance, the problem is probably located on your initial geometry... Packmol generates *good* solvation spheres, but they are far away from the minimal energy geometry, and probably the forces on atoms at the first step are very high. This explains why your kinetic energy rised a lot, a huge amount of potential energy is changing to kinetic energy. When I run this type of AIMD, usually relax the structure first with classic MD with OPLS-AA forcefield, to get a better geometrical guess to start the AIMD. Otherwise, you can do a geometrical relaxation (a non-tight relaxation) in ORCA at the same level to get a better geometrical guess. <div>Also, I would add D3(BJ) dispersion correction to your calculation.</div><div><br></div><div>Regards - Lucas</div></div><br><div class="gmail_quote"></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El mar, 23 ago 2022 a las 13:57, Jessie wang (<<a href data-email-masked rel="nofollow">wangch...@gmail.com</a>>) escribió:<br></div></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">Dear all,<div><br></div><div>I'm trying to do NVT AIMD for a system (a cluster surrounded by water molecules). The water molecules was generated by packmol with tolerance 2.0.</div><div>I was set temperature to 298.15K, but after several steps, it rised to hundreds even thousands K.  I was used NOSE-HOOVER thermosat, and tried CVSR also. it does not work.</div><div><br></div><div> Here I attached my input file, and any help would really appreciated. </div><div><br></div><div>Thanks</div><div><br></div><div><br></div><div>Best,</div><div>Jessie</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 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">cp2k+uns...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/7856dc7f-d236-477b-8c20-98fdd1d3037bn%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/7856dc7f-d236-477b-8c20-98fdd1d3037bn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1661392961913000&usg=AOvVaw0R1ren7uoVeY_Nilj-LnNp">https://groups.google.com/d/msgid/cp2k/7856dc7f-d236-477b-8c20-98fdd1d3037bn%40googlegroups.com</a>.<br>
</blockquote></div>
</blockquote></div>

<p></p>

-- <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="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/5dd3975a-13ef-451e-9b73-ff023bb05605n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/5dd3975a-13ef-451e-9b73-ff023bb05605n%40googlegroups.com</a>.<br />