So should we lose all hope in computing the atomic forces in the way I had imagined?<div> </div><div>Best regards,</div><div>Riccardo Piombo </div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Il giorno martedì 10 ottobre 2023 alle 12:02:47 UTC+2 Jürg Hutter ha scritto: </div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">No, there are random contributions to the propagation that you will not be able to retrace. JH ________________________________________ From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of riccardo piombo <<a href data-email-masked rel="nofollow">riccard...@gmail.com</a>> Sent: Tuesday, October 10, 2023 10:37 AM To: cp2k Subject: Re: [CP2K:19335] Extracting Atomic Forces from XYZ Trajectories in CP2K Hi Prof. Hutter, my simulations are also based on the Langevin ensamble. do we have to give up on calculating atomic forces from the xyz trajectory? Kind regards, Riccardo Piombo Il giorno martedì 10 ottobre 2023 alle 10:17:17 UTC+2 Simone Ritarossi ha scritto: Hi users, I have the same problem and I used Langevin (second-generation Car-Parrinello). In this case? Simone Ritarossi Il Mar 10 Ott 2023, 08:58 Jürg Hutter <<a href data-email-masked rel="nofollow">hut...@chem.uzh.ch</a>> ha scritto: Hi all integrators in CP2K are based on the Velocity Verlet algorithm. The different routines can be found in the src/motion directory. If your trajectory was generated for a NVE ensemble without any constraints, you can easily extract the forces from inverting the Velocity Verlet algorithm. If you were using any type of thermostat or any geometrical constraints it will be very difficult or impossible. At least I don't see a straightforward way. Depending on you final goal, recalculating the forces at a subset of your trajectory will be your only option. regards JH ________________________________________ From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of riccardo piombo <<a href data-email-masked rel="nofollow">riccard...@gmail.com</a>> Sent: Monday, October 9, 2023 5:54 PM To: cp2k Subject: [CP2K:19326] Extracting Atomic Forces from XYZ Trajectories in CP2K Dear CP2K Community, We have encountered a challenge while working with MD trajectories in XYZ format. Specifically, we have trajectories containing atomic positions and velocities for each frame, but unfortunately, we do not have the corresponding atomic forces. In theory, it is possible to calculate atomic forces based on atomic positions and velocities. However, we are unsure about the specific algorithm and where it can be found within the CP2K codebase. We are reaching out to the community to inquire if anyone can provide guidance on: 1. 1. Locating the algorithm or module responsible for calculating atomic forces within CP2K. 2. 2. Assessing the feasibility of extracting atomic forces from our XYZ trajectories in a reasonably short amount of time. Any assistance or insights on this matter would be greatly appreciated. We understand that this may be a non-trivial task, but we are hopeful that the CP2K community's expertise can help us navigate this challenge. Thank you in advance for your support and expertise. Best regards Riccardo Piombo -- 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 data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/311dfc7b-f7ec-4d45-a7c4-32e9203b8ee1n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://groups.google.com/d/msgid/cp2k/311dfc7b-f7ec-4d45-a7c4-32e9203b8ee1n%2540googlegroups.com&source=gmail&ust=1697025479275000&usg=AOvVaw2xnqOrytnPPnXpGJ7EUQCf">https://groups.google.com/d/msgid/cp2k/311dfc7b-f7ec-4d45-a7c4-32e9203b8ee1n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/311dfc7b-f7ec-4d45-a7c4-32e9203b8ee1n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://groups.google.com/d/msgid/cp2k/311dfc7b-f7ec-4d45-a7c4-32e9203b8ee1n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1697025479275000&usg=AOvVaw3Z64TAbajxZA4zgKWvCobu">https://groups.google.com/d/msgid/cp2k/311dfc7b-f7ec-4d45-a7c4-32e9203b8ee1n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>. -- 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 data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ZR0P278MB0759B06368A3D5E374E749D39FCDA%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://groups.google.com/d/msgid/cp2k/ZR0P278MB0759B06368A3D5E374E749D39FCDA%2540ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM&source=gmail&ust=1697025479275000&usg=AOvVaw1iGJJLnCsgnQoDE4TiuKHO">https://groups.google.com/d/msgid/cp2k/ZR0P278MB0759B06368A3D5E374E749D39FCDA%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM</a>. -- 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 data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/b9230076-a40b-4e74-8b91-4d90d9ae180dn%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://groups.google.com/d/msgid/cp2k/b9230076-a40b-4e74-8b91-4d90d9ae180dn%2540googlegroups.com&source=gmail&ust=1697025479275000&usg=AOvVaw2wvJGySJRlCnlxCICaSPD0">https://groups.google.com/d/msgid/cp2k/b9230076-a40b-4e74-8b91-4d90d9ae180dn%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/b9230076-a40b-4e74-8b91-4d90d9ae180dn%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://groups.google.com/d/msgid/cp2k/b9230076-a40b-4e74-8b91-4d90d9ae180dn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1697025479275000&usg=AOvVaw0G8_27XaC3l-Ud_G63Dd1j">https://groups.google.com/d/msgid/cp2k/b9230076-a40b-4e74-8b91-4d90d9ae180dn%40googlegroups.com?utm_medium=email&utm_source=footer</a>>. </blockquote></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/05ca489e-163b-4bc7-98f6-2f4206b5e7d0n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/05ca489e-163b-4bc7-98f6-2f4206b5e7d0n%40googlegroups.com</a>.