<div>UPD on this issue: apparently, the Cranc-Nicolson is the less stable one. I did some additional tests with a very small time-step (0.0005 fs) and the CN propagator converged to the ETRS result, the same one that ETRS gives even with dt = 0.006 fs. So, the answer is, one should use the ETRS propagator to be able to use larger time step.</div><div><br /></div><div>Natalia<br /></div><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, October 9, 2023 at 2:06:33 PM UTC+2 Natalia K 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>Dear CP2K team,</div><div><br></div><div>I performed simple tests for H2 with cp2k-9.1 using an input file from <b>cp2k/tests/QS/regtest-rtp-2</b> and modifying the external field (the input file is attached). I then performed 4 tests, using ETRS or CN propagators and two different time steps, dt=0.25 a.u. and 0.125 a.u. With the ETRS propagator, the energies as well as the H atom dynamics are the same for both dt values. However, the results with the Crank-Nicolson propagator (CN) are not only different for two time steps, but also are completely different from when using ETRS. I attach the results here. I tried refining some other parameters, but the two propagators always lead to different EMD.<br></div><div><br></div><div>My question: when performing Ehrenfest dynamics with applied external field, which propagator should I trust? Why they give such different results?</div><div><br></div><div>Best regards,</div><div><br></div><div>Natalia<br></div><div><br></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/45a76931-2f65-4883-81b0-4951c93e9028n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/45a76931-2f65-4883-81b0-4951c93e9028n%40googlegroups.com</a>.<br />