<div>Dear Guillaume,</div><div><br /></div><div>thank you for your instructions. After applying the patch, and making sure the test files are correct, I ran calculations comparing the latest version 2023.2 (with the patch) with the version 9.1, that I was previously using, for the system Ag6-N2. The results are different. Namely, the kinetic energy is significantly different (and consequently the atomic motion). This calculation is similar to the one from https://dx.doi.org/10.1021/acs.jpcc.0c02979?ref=pdf, where on Fig.6 (middle left) they show the N-N bond change with time. Their result is obtained with the Gaussian code. With all the parameters similar to theirs (except the XC functional, they use hybrid and I use GGA, but I don't think it is important), CP2K-v9.1 gives similar result, while v2023.2 does not. I attach my results and the input file I used in both calculations below.</div><div><br /></div><div>Best regards,</div><div><br /></div><div>Natalia<br /></div><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, October 12, 2023 at 10:54:18 AM UTC+2 Guillaume Le Breton 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;">Dear CP2K community, <br><br>On 24 February 2023 I made some changes in the Real-Time calculation part of the code which turned out to break one specific behavior: Ehrenfest Dynamics with a time-dependent field. This has been fixed the 9 October 2023. <br><br>Thus, if you are using Ehrenfest Dynamics with a time-dependent field, you are NOT impacted if:<br><br>- You are using a version before 24.02.2023. This includes V2023.1 and the previous ones. <br>   - The master branch from 09.10.2023. <br><br>You ARE impacted if:<br>- You are using a version between 24.02.2023 and 09.10.2023. This includes V2023.2. <br><br>If you are impacted, I recommend you rerun your simulations (involving Ehrenfest + Field) with the newest version of the code as the issue is quite important. If you want to keep using version 2023.2, here is a patch to solve this problem. <br><br>1) Download the patch and untar it.<br>2) Apply the patch using the patch command in the shell: <br><br>$> patch cp2k/src/qs_force.F qs_force.patch<br>$> patch cp2k/tests/QS/regtest-rtp-2/TEST_FILES TEST_FILES.patch<br><br>To check that the patch has been applied correctly, open the cp2k/tests/QS/regtest-rtp-2/TEST_FILES file. Line 11 should be:<br><br>H2-emd-efield.inp                                      2      4e-11                -0.894818188949<br><br>3) recompile the code and run regtests. <br><br>I sincerely apologize for this problem. <br><br>Guillaume Le Breton</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/10dea554-5950-4416-992d-49eec123ce12n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/10dea554-5950-4416-992d-49eec123ce12n%40googlegroups.com</a>.<br />