<div dir="ltr">Thanks, the keyword EVAL_FORCES does work.</div> <div class="gmail_quote"><div dir="ltr" class="gmail_attr">Krack Matthias <<a href="mailto:matthias.krack@psi.ch">matthias.krack@psi.ch</a>> 于2023年9月7日周四 21:50写道： </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg-3451602874250218323"> <div lang="en-CH" style="overflow-wrap: break-word;"> <div class="m_-3451602874250218323WordSection1"> Hi Add the keyword EVAL_FORCES to the &REFTRAJ section or use the latest CP2K development version after <a href="https://github.com/cp2k/cp2k/pull/2974" target="_blank">PR2974</a> has been merged successfully. Alternatively, you can apply the same patch to your local CP2K version. HTH Matthias <div id="m_-3451602874250218323mail-editor-reference-message-container"> <div> <div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(181,196,223);padding:3pt 0cm 0cm"> From: <a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a> <<a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a>> on behalf of Tianhua Wang <<a href="mailto:njuwth@gmail.com" target="_blank">njuwth@gmail.com</a>> Date: Thursday, 7 September 2023 at 09:12 To: cp2k <<a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a>> Subject: [CP2K:19225] Re: The force output of RETRAJ module under a r2scan functional is zero </div> In fact, according to my tests, this bug exists in almost all CP2K versions later than 6.1, i.e., 7.1, 8.1, 8.2, 9.1, 2022.1, and 2023.2 (I did not test 2023.1), for both SCAN and r2SCAN functionals. Instead, CP2K 6.1 can output the forces normally. It seems to be a systemic problem. <div> <div> 在2023年9月7日星期四 UTC+8 11:20:17<qin li> 写道： </div> <blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin-left:4.8pt;margin-right:0cm"> <div> The force output of the RETRAJ module works problematiclly in version cp2k2022.1. </div> <div> <div> qin li <<a href="mailto:flytomoo...@gmail.com" target="_blank">flytomoo...@gmail.com</a>> 于2023年9月7日周四 10:20写道： </div> </div> <div> <blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin-left:4.8pt;margin-right:0cm"> <div> Hi, <div> </div> <div> I'm trying to train a deep potential model of the MgCl2 solution based on r2scan functional, where the positions and forces are needed. I employed the r2scan functional to rerun the trajectory resulting from the BLYPD3 functional, but the force output of the RETRAJ module comes out to be all zero. However, the RETRAJ module works fine for the BLYPD3 functional. Does anyone have an idea of how to resolve this? </div> <div> </div> <div> The cp2k version is cp2k2022.1, inputs of the r2scan and force output are as follows and the attachment: </div> <div> &XC &XC_FUNCTIONAL &MGGA_X_R2SCAN &END MGGA_X_R2SCAN &MGGA_C_R2SCAN &END MGGA_C_R2SCAN &END XC_FUNCTIONAL &XC_GRID XC_DERIV SPLINE3 XC_SMOOTH_RHO NONE &END XC_GRID &END XC </div> <div> </div> <div> &REFTRAJ EVAL_ENERGY_FORCES TRAJ_FILE_NAME ${REFTRAJ} &END REFTRAJ </div> <div> </div> <div> &PRINT &FORCES &EACH MD 1 &END EACH &END FORCES </div> <div> &END PRINT </div> <div> </div> </div> </blockquote> </div> </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" target="_blank">cp2k+unsubscribe@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/760142b3-db53-4c9a-8ca2-882cccc136ban%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank"> https://groups.google.com/d/msgid/cp2k/760142b3-db53-4c9a-8ca2-882cccc136ban%40googlegroups.com</a>. </div> </div> </div> </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" target="_blank">cp2k+unsubscribe@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ZRAP278MB0827DF9B11ACCAF76AD4DE17F4EEA%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/ZRAP278MB0827DF9B11ACCAF76AD4DE17F4EEA%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM</a>. </div></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/CALSfb8o4W20%2BvdUhMvyyiuo16FSz_JNc-TNjMw4P9woiWWk_Dg%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CALSfb8o4W20%2BvdUhMvyyiuo16FSz_JNc-TNjMw4P9woiWWk_Dg%40mail.gmail.com</a>.