<div>Dear all and dear author,</div><div><br /></div><div>Did the Input suggetions with "
FORCE_EVAL / PROPERTIES / ATOMIC ENERGY .TRUE." work for you? If so, would you mind sharing your input as an example?<br /></div><div><br /></div><div>I am also trying to run an AIMD simulation and would like to print the atomic forces and energies too. When I use the block:<br /></div><div>&PROPERTIES<br /> &ATOMIC<br /> ENERGY TRUE<br /> &END ATOMIC<br />&END PROPERTIES</div><div><br /></div><div>the calculation the SCF converges fine and the atomic forces are printed for the first timestep but then the job crashes with the error:</div><div> *******************************************************************************<br /> * ___ *<br /> * / \ *<br /> * [ABORT] *<br /> * \___/ CPASSERT failed *<br /> * | *<br /> * O/| *<br /> * /| | *<br /> * / \ force_env_methods.F:430 *<br /> *******************************************************************************<br /><br /><br /> ===== Routine Calling Stack =====<br /><br /> 2 qs_mol_dyn_low<br /> 1 CP2K<br /></div><div><br /></div><div>This behavior does not seem to be dependent on the basisset (I tested it for DZVP and TZVP MOLOPT) or on the CP2K version (I tested it with the 2023.1, 2023.2 and 2024.1 versions).</div><div><br /></div><div>I hope this is the right place for this question, otherwise I will open a new topic.</div><div>Thanks in advance!</div><div>Jakob<br /></div><div><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Nathalie Smith schrieb am Dienstag, 27. Februar 2024 um 13:34:39 UTC+1:<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>Hi,</div><div><br></div>Thank you so much for your answer and the explanation!<br><br><div>Best wishes</div><div>Nathalie Smith</div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Jürg Hutter schrieb am Mittwoch, 21. Februar 2024 um 17:17:01 UTC+1:<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
<br>
<br>atomic energies are not uniquely defined in QM calculations. They are not
<br>real observables and are usually not considered.
<br>However, there is a simple atomic energy decomposition available in CP2K
<br>that works for QM systems using Mulliken populations and some other
<br>CP2K specific assumptions for the energy decomposition. You can
<br>get it with
<br>
<br>FORCE_EVAL / PROPERTIES / ATOMIC ENERGY .TRUE.
<br>
<br>No warenty that the numbers are of any real use.
<br>
<br>regards
<br>JH
<br>
<br>________________________________________
<br>From: <a rel="nofollow">cp...@googlegroups.com</a> <<a rel="nofollow">cp...@googlegroups.com</a>> on behalf of Nathalie Smith <<a rel="nofollow">nathali...@gmail.com</a>>
<br>Sent: Wednesday, February 21, 2024 4:46 PM
<br>To: cp2k
<br>Subject: [CP2K:19943] Re: Energies of the atoms of my system
<br>
<br>Hi, thank you so much for your reply!
<br>
<br>I tested to set the input section like this:
<br>
<br>&MOTION
<br> &MD
<br> &PRINT
<br> &ENERGY
<br> FILENAME energy.dat
<br> &EACH
<br> MD 1
<br> &END EACH
<br> &END ENERGY
<br> &END PRINT
<br> [...]
<br> &END MD
<br>&END MOTION
<br>
<br>The output I got was a regular .ener-file with the system energies, not the energies of each atom. If anyone has an idea how I can get the energies of each atom in my simulation - or if someone knows whether this option is simply not implemented in CP2K - I'd love you to let me know.
<br>
<br>Best wishes
<br>Nathalie Smith
<br>Muhammad Saleh schrieb am Mittwoch, 21. Februar 2024 um 15:16:57 UTC+1:
<br>Hi..
<br>
<br>You may have to check section CP2K_INPUT / MOTION / MD / PRINT / ENERGY, on CP2K documentation, set FILENAME and EACH. Hope it will help
<br>
<br>Best
<br>MuS
<br>
<br>On Tuesday, February 20, 2024 at 2:26:59 PM UTC+1 Nathalie Smith wrote:
<br>Hello everyone,
<br>
<br>I am trying to run a MD simulation and want to print the energies of each atom, similarly to the force output of MOTION/PRINT/FORCES. Could someone kindly tell me where to find such an option in CP2K? Thanks in advance!
<br>
<br>Best wishes,
<br>Nathalie Smith
<br>
<br>--
<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 rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a 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/f0e4c7fe-db7d-493c-91e6-25012439177an%40googlegroups.com" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=de&q=https://groups.google.com/d/msgid/cp2k/f0e4c7fe-db7d-493c-91e6-25012439177an%2540googlegroups.com&source=gmail&ust=1725961870161000&usg=AOvVaw1lSyQbZrxE6YbRoOB32ZlH">https://groups.google.com/d/msgid/cp2k/f0e4c7fe-db7d-493c-91e6-25012439177an%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/f0e4c7fe-db7d-493c-91e6-25012439177an%40googlegroups.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=de&q=https://groups.google.com/d/msgid/cp2k/f0e4c7fe-db7d-493c-91e6-25012439177an%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1725961870161000&usg=AOvVaw2OQX8XtTXYsFG3Jo8unAPS">https://groups.google.com/d/msgid/cp2k/f0e4c7fe-db7d-493c-91e6-25012439177an%40googlegroups.com?utm_medium=email&utm_source=footer</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/f01483cb-150b-4fa2-90df-2445fd4af40bn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/f01483cb-150b-4fa2-90df-2445fd4af40bn%40googlegroups.com</a>.<br />