<div dir="ltr"><div>Thanks for the feedback.<br></div><div><br></div><div>It was possible to go around this problem at an early stage of SIRIUS integration to cp2k but it is not possible to do this anymore. I can come up with a solution for the vdw correction but it will require modifying cp2k accordingly. I would suggest you to open an issue on cp2k github site (as feature to have) so I can track this. It should not take long to include this feature to cp2k.For the time being i have no easy solution to give you unless you are willing to add one line of code for testing and recompiling.</div><div><br></div><div>As a temporary solution add this (it turns it on permanently)<br></div><div><br></div><div>CALL sirius_add_xc_functional(sctx, "XC_FUNC_VDWDF")</div><div><br></div><div>in sirius_interface.f after the line 186 and recompile. Remove the FUNC_VDWXC keyword from the XC section and run it again. <b><br></b></div><div><b><br></b></div><div><b>Warning: </b>it is a temporary solution for <b>testing </b>but it will give me time to actually implement the real fix. <br></div><div><br></div><div>Best</div><div>M.<br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le lun. 5 sept. 2022 à 10:09, Dialogicus <<a href="mailto:dialogicus@gmail.com">dialogicus@gmail.com</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>Thank you for your advice. U<span dir="ltr"><span>nfortunately</span></span>, it doesn't work. <span dir="ltr"><span></span></span> CP2K can't identify FUNC_VDWDF as a keyword or subsection. The out file says</div><div><br></div><div><br></div><div><font size="2">Possible matches for unknown subsection <br><br>FUNC_VDWDF<br><br> enum BEEFVDW in section %__ROOT__%FORCE_EVAL%DFT%KG_METHOD%XC%XC_FUNCTIONAL for keyword _SECTION_PARAMETERS_ score: 103<br> subsection LDA_C_VWN in section %__ROOT__%FORCE_EVAL%DFT%KG_METHOD%XC%XC_FUNCTIONAL score: 103<br> subsection MGGA_C_DLDF in section %__ROOT__%FORCE_EVAL%DFT%KG_METHOD%XC%XC_FUNCTIONAL score: 103<br> subsection GGA_C_CCDF in section %__ROOT__%FORCE_EVAL%DFT%KG_METHOD%XC%XC_FUNCTIONAL score: 103<br> enum BEEFVDW in section %__ROOT__%FORCE_EVAL%DFT%ENERGY_CORRECTION%XC%XC_FUNCTIONAL for keyword _SECTION_PARAMETERS_ score: 103<br><br> *******************************************************************************<br> * ___ *<br> * / \ *<br> * [ABORT] *<br> * \___/ unknown subsection FUNC_VDWDF of section XC_FUNCTIONAL *<br> * | *<br> * O/| *<br> * /| | *<br> * / \ input/input_parsing.F:187 *<br> *******************************************************************************</font></div><div><br></div><div><br></div><div>As far as I understand, this error means invalid syntax for CP2K, but not for SIRIUS. So, is there any way to send information about vdW functional directly to SIRIUS run (e.g. SIRIUS includes this functional automatically in each run)?</div><div><br></div><div>--</div><div>D.<br></div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr"><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>Hi,</div><div><br></div><div>CP2K implementation of vdw is completely independent of libvdwxc which is used in SIRIUS. It means in practice that <br></div><div><ul><li>there is no direct path in place in cp2k to parse the VDW_POTENTIAL section and find the relevant functional supported in libvdwxc (and SIRIUS).It should probably be added if possible.<br></li><li>the VDW_POTENTIAL section is simply ignored when SIRIUS is initialized. It means that the calculations are correct but simply ignore the vdw corrections.<br></li><li>CP2K documentation is misleading</li></ul><div>If you wish to try a vdwxc functional with SIRIUS, then have a look at the functional supported by libvdwxc (<a href="https://libvdwxc.gitlab.io/libvdwxc/background.html#van-der-waals-forces" rel="nofollow" target="_blank">https://libvdwxc.gitlab.io/libvdwxc/background.html#van-der-waals-forces</a>) that fit your needs and then add it to the list of functional XC_FUNC_VDWDF, XC_FUNC_VDWDF2, or XC_FUNC_VDWDFCX to the libxc section of your input file. <br></div><div><br></div><div>something like that might actually work.<br></div></div><div><br></div><div>&XC_FUNCTIONAL<br> &GGA_X_PBE_SOL<br> &END GGA_X_PBE_SOL<br> &GGA_C_PBE_SOL<br> &END GGA_C_PBE_SOL<br></div><div> &FUNC_VDWDF</div><div> &END FUNC_VDWDF</div><div>&END XC_FUNCTIONAL</div><div><br></div><div>let me know if it works or not.<br></div><div><br></div><div>Best</div><div><br></div><div>M.</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Le vendredi 2 septembre 2022 à 11:02:29 UTC+2, Dialogicus a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span>Dear CP2K users,</span>
<p><span>As
CP2K manual says it’s possible to use Van der Waals functionals in
the plane-wave calculations with SIRIUS.</span></p>
<p><span>I tried to run <span>Fm3m gold cell optimization using</span> pseudopotential <span><span>PW method</span></span>. VdW interactions were taken into account in
VDW_POTENTIAL subsection with NON_LOCAL type of dispersion functional. However, it seems there’s no difference between
calculations with and without dispersion correction. I suspect the
vdW functionals don’t work in my case.</span></p><p><span><span>My CP2K and SIRIUS are both linked with libvdwxc
library. </span></span></p>
<p><span>I
don’t know what I’m doing wrong. Is any ideas?</span></p>
<p><span>My
input and output files (with and without vdW correction) are attached
below.<br></span></p>
</blockquote></div></blockquote></div>
<p></p>
-- <br>
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.<br>
To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/abi2t_F_mQ4/unsubscribe" target="_blank">https://groups.google.com/d/topic/cp2k/abi2t_F_mQ4/unsubscribe</a>.<br>
To unsubscribe from this group and all its topics, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/f2c2235c-79f9-48dd-ae85-a90c7c2034c3n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/f2c2235c-79f9-48dd-ae85-a90c7c2034c3n%40googlegroups.com</a>.<br>
</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/CA%2B%3DybGbPc-_gd4yCrF7ACYgF-8pg5iw7tv9T83z6HXaC%3Di8tSQ%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CA%2B%3DybGbPc-_gd4yCrF7ACYgF-8pg5iw7tv9T83z6HXaC%3Di8tSQ%40mail.gmail.com</a>.<br />