Dear 

Marcella<div>Thank you very much for your reply, this problem has troubled me for a long time. When I set EPS_SCF to 1.0E-6 during NVT BOMD simulation, the conserved quantity always shows a straight line drift, while the time spent is basically stable. Is it true that the conserved quantity will increase stably in NVT simulation? By the way, can SGCPMD improve the sampling efficiency of this system?</div><div><br></div><div>Kind regards</div><div>X. J  Xu</div><div><img alt="4.png" data-iml="599293.400000006" width="892px" height="683px" src="cid:94f094cc-7e49-47f6-be8c-c7fdce6c4827"><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2022年10月14日星期五 UTC+8 15:26:16<Marcella Iannuzzi> 写道：<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear X.J.Xu<div><br></div><div>Such total energy drifts indicate that the results of the simulation are not reliable, the errors in the integration of the equation of motion are to big. </div><div>The conserved quantity must be conserved.</div><div>The increase in time per MD step is most probably related to an increase in the number of SCF iterations, if you are running BOMD. </div><div>If this is the case, i.e., if it becomes more and more difficult to optimise the wave function, it might indicate that the computed forces are rather wrong and lead to the sampling of weird geometries.</div><div>By the way, have you considered that graphene is a semi-metal? This makes the convergence of the wave function even more problematic. </div><div>Before worrying on how to speed up the sampling, it is recommendable to find settings that produce reliable results, for instance further reducing EPS_SCF, adding virtual MOS, increasing the cutoff, etc.</div><div><br></div><div>Kind regards</div><div> Marcella</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, October 14, 2022 at 4:22:59 AM UTC+2 <a href data-email-masked rel="nofollow">pxp11...@gmail.com</a> 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">Hi cp2k groups,<br><div>I'm a new cp2k user and recently i simulated the contact Angle of water on graphene using AIMD. I found that different values of EPS_SCF lead to different changes in cons qty, thus leading to different increases in used time. As shown in the following three figures, when cons qty is almost unchanged, the time taken by each MD step is the same, but when cons qty rises, usedtime rises sharply, may I ask what is the problem? How can I improve it? The following is my input file.<br></div><div>EPS_SCF = 1.0*10^-5, cutoff = 320 Ry</div><div><img alt="1E-5.png" width="534px" height="409px" src="https://groups.google.com/group/cp2k/attach/1027e835dd231/1E-5.png?part=0.3&view=1"><br></div><div>5.0*10^-6, cutoff =  380Ry</div><div><img alt="5E-6.png" width="534px" height="409px" src="https://groups.google.com/group/cp2k/attach/1027e835dd231/5E-6.png?part=0.2&view=1"></div><div>2*10-6,  cutoff = 380 Ry</div><div><img alt="2E-6.png" width="534px" height="409px" src="https://groups.google.com/group/cp2k/attach/1027e835dd231/2E-6.png?part=0.4&view=1"><br></div><div>I used NVT 330 K and nose-hoover chain, and here is my input file.<br></div><div>My questions are: 1. Does cons Qty remain constant over time once NVT reaches equilibrium? 2. Why does the cons Qty in 2*10-6 fluctuate wildly after it has remained unchanged for a long time? 3. How to accelerate AIMD. <br></div><div><br></div><div>Best regards, </div><div>X.J Xu</div><div><br></div><div><br></div><div><br></div></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/2121f47f-ed89-4fa8-8fc9-48e85189dd9an%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/2121f47f-ed89-4fa8-8fc9-48e85189dd9an%40googlegroups.com</a>.<br />