<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Dear Patrick Gono,<div class=""><br class=""></div><div class="">Thank you so much for your help. What you write makes perfect sense. However, I am still stymied by the example, as it seems to include two exchanges and two correlations.</div><div class="">XC_GGA_X_RPW86 includes both exchange and correlations</div><div class="">and so does VWN.</div><div class="">So what does the code do? It seems, at least in the energy, to ADD GGA and LDA xc (both x and c), </div><div class="">and that is certainly not what is wanted.</div><div class=""><br class=""></div><div class="">There is also the issue of using pseudopotentials that correspond to the xc for the bands.There are limitless combinations of XC possible, but only a few are represented in the cp2k pseudopotentials/\</div><div class=""><br class=""></div><div class="">Sincerely,</div><div class=""><br class=""></div><div class="">Ron</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><blockquote type="cite" class=""><div class="gmail_quote"><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex;"><div dir="ltr" class=""><div class=""><p class="gmail-m_457967070791147240p1">&XC</p></div></div></blockquote></div></blockquote><blockquote type="cite" class=""><div class="gmail_quote"><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex;"><div dir="ltr" class=""><div class=""><p class="gmail-m_457967070791147240p1">&XC_FUNCTIONAL NO_SHORTCUT</p><p class="gmail-m_457967070791147240p1">&LIBXC T</p><p class="gmail-m_457967070791147240p1">FUNCTIONAL XC_GGA_X_RPW86</p><p class="gmail-m_457967070791147240p1">&END LIBXC </p><p class="gmail-m_457967070791147240p1">&VWN T </p><p class="gmail-m_457967070791147240p1">&END VWN</p><p class="gmail-m_457967070791147240p1">&END XC_FUNCTIONAL &VDW_POTENTIAL</p><p class="gmail-m_457967070791147240p1">POTENTIAL_TYPE NON_LOCAL &NON_LOCAL</p><p class="gmail-m_457967070791147240p1">TYPE LMKLL</p><p class="gmail-m_457967070791147240p1">VERBOSE_OUTPUT T KERNEL_FILE_NAME ./vdW_kernel_table.dat CUTOFF 160</p><p class="gmail-m_457967070791147240p1">&END NON_LOCAL &END VDW_POTENTIAL</p><p class="gmail-m_457967070791147240p1">&END XC</p></div></div></blockquote></div></blockquote><div class=""> <div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">---<br class="">Ron Cohen<br class=""><a href="mailto:reco...@gmail.com" class="">reco...@gmail.com</a><br class="">skypename: ronaldcohen</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">twitter: @recohen3<br class=""><div class=""><br class=""></div></div><br class="Apple-interchange-newline"><br class="Apple-interchange-newline"> </div> <div><br class=""><blockquote type="cite" class=""><div class="">On Feb 12, 2019, at 6:12 PM, Patrick Gono <<a href="mailto:patric...@gmail.com" class="">patric...@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div class="">Dear Ronald,</div><div class=""><br class=""></div><div class="">it is true that documentation can be sparse. I suggest following the articles cited in the manual for more information.</div><div class=""><br class=""></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=""><p class="gmail-m_457967070791147240p1">I note some related problems with the documentation. I find examples like for LIBXC:</p><div class=""> <br class="webkit-block-placeholder"></div><p class="gmail-m_457967070791147240p1">FUNCTIONAL XC_MGGA_X_M06_L XC_MGGA_C_M06_L</p><p class="gmail-m_457967070791147240p1">but the docs at </p><p class="gmail-m_457967070791147240p1">say</p><ul class="gmail-m_457967070791147240disc" style="font-family: Times; font-size: inherit;"><li style="margin-left:0em;padding-left:0em;text-indent:0em" class=""><a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/XC/XC_FUNCTIONAL/LIBXC.html#list_FUNCTIONAL" id="gmail-m_457967070791147240FUNCTIONAL" target="_blank" class="">FUNCTIONAL</a></li><li style="margin-left:0em;padding-left:0em;text-indent:0em" class="">This keyword cannot be repeated and it expects precisely one word.<br class=""></li></ul><p class="gmail-m_457967070791147240p1">Yet there are "two words" </p><p class="gmail-m_457967070791147240p1">and including VWN afterwards also seems like more than one functional, as I said above.</p></div></blockquote><div class=""><br class=""></div><div class="">In older versions (prior to 6.1), the FUNCTIONAL keyword accepted a list of words. In the newest version (6.1), this keyword only accepts a single word. Note, however, that the section &LIBXC can be repeated. So, instead of:</div><div class=""><br class=""></div><div class="">&LIBXC<br class=""></div><div class=""> FUNCTIONAL XC_MGGA_X_M06_L XC_MGGA_C_M06_L</div><div class="">&END LIBXC</div><div class=""><br class=""></div><div class="">you can do:</div><div class=""><br class=""></div><div class=""><div class="">&LIBXC<br class=""></div><div class=""> FUNCTIONAL XC_MGGA_X_M06_L</div><div class="">&END LIBXC</div><div class=""><div class="">&LIBXC<br class=""></div><div class=""> FUNCTIONAL XC_MGGA_C_M06_L</div><div class="">&END LIBXC</div><div class=""><br class=""></div><div class="">You can also mix different correlation or exchange functionals, like so:</div><div class=""><br class=""></div><div class="">&XC_FUNCTIONAL<br class=""></div><div class=""> &LIBXC<br class=""> FUNCTIONAL XC_GGA_X_RPW86 <br class=""> &END LIBXC</div><div class=""><div class=""> &LIBXC<br class=""> FUNCTIONAL GGA_C_PBE<br class=""> &END LIBX</div><div class="">&END XC_FUNCTIONAL<br class=""></div><div class="">&vdW_POTENTIAL<br class=""> DISPERSION_FUNCTIONAL NON_LOCAL<br class=""> &NON_LOCAL<br class=""> TYPE RVV10<br class=""> PARAMETERS 9.3 0.0093<br class=""> KERNEL_FILE_NAME ./rVV10_kernel_table.dat<br class=""> CUTOFF 400<br class=""> &END NON_LOCAL<br class="">&END vdW_POTENTIAL<br class=""><br class=""></div><div class=""><br class=""></div><br class=""></div><div class="">where I used RPW86 exchange and PBE correlation. The reason behind that is that I am trying to use the (modified) rVV10 nonlocal van der Waals functional. In the original paper: <a href="https://journals.aps.org/prb/abstract/10.1103/PhysRevB.87.041108" class="">https://journals.aps.org/prb/abstract/10.1103/PhysRevB.87.041108</a> the authors define the total exchange-correlation as:</div><div class="">XC_rVV10 = X_RPW86 + C_PW_LDA + C_non-local_vdW</div><div class=""><br class=""></div><div class="">However, in later work it was found that the non-local vdW correction from the rVV10 works quite well with other combinations of correlation and exchange functionals. Hence, I am now using<br class=""></div></div><div class="">XC_rVV10 = X_RPW86 + C_PBE + C_non-local_vdW</div><div class="">with a specific choice of the b parameter in the vdW correction.</div><div class=""><br class=""></div><div class="">What I am trying to say is that you can mix and match the non-local vdW correction with different exchanges and correlations. The specific choice will always depend on the system you are trying to study. There are many papers testing and comparing various combinations, such as this one <a href="https://pubs.acs.org/doi/10.1021/ct4003527" class="">https://pubs.acs.org/doi/10.1021/ct4003527</a> where the VV10 corrections is added to various exchanges and correlations.</div><div class=""><br class=""></div><div class="">In your input file you use the Lee-Murray-Kong-Lundqvist-Langreth (LMKLL) nonlocal van der Waals density functional (via the keyword TYPE LMKLL). If you look up the relevant publication, you will find what correlation and exchange functional it is supposed to be used with. You may also find benchmark papers that use this non-local correction on top of other-than-originally-intended exchange and correlation functionals. This allows one to choose the best combination for the system at hand (large system vs small; periodic vs isolated; metals vs semiconductors vs small molecules; etc.).</div><div class=""><br class=""></div><div class="">Hope this helps.</div><div class="">Yours sincerely,</div><div class="">Patrick Gono<br class=""></div> </div></div></div></div></div></div></div><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, 12 Feb 2019 at 15:02, Ronald Cohen <<a href="mailto:reco...@gmail.com" class="">reco...@gmail.com</a>> wrote:<br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr" class="">I am quite confused about the (lack of) documentation for the XC functionals and non-local vdW. I find examples such as:<div class=""><br class=""></div><div class=""><p class="gmail-m_457967070791147240p1">&XC</p><p class="gmail-m_457967070791147240p1">&XC_FUNCTIONAL NO_SHORTCUT</p><p class="gmail-m_457967070791147240p1">&LIBXC T</p><p class="gmail-m_457967070791147240p1">FUNCTIONAL XC_GGA_X_RPW86</p><p class="gmail-m_457967070791147240p1">&END LIBXC </p><p class="gmail-m_457967070791147240p1">&VWN T </p><p class="gmail-m_457967070791147240p1">&END VWN</p><p class="gmail-m_457967070791147240p1">&END XC_FUNCTIONAL &VDW_POTENTIAL</p><p class="gmail-m_457967070791147240p1">POTENTIAL_TYPE NON_LOCAL &NON_LOCAL</p><p class="gmail-m_457967070791147240p1">TYPE LMKLL</p><p class="gmail-m_457967070791147240p1">VERBOSE_OUTPUT T KERNEL_FILE_NAME ./vdW_kernel_table.dat CUTOFF 160</p><p class="gmail-m_457967070791147240p1">&END NON_LOCAL &END VDW_POTENTIAL</p><p class="gmail-m_457967070791147240p1">&END XC</p><p class="gmail-m_457967070791147240p1"><br class=""></p><p class="gmail-m_457967070791147240p1">(from <a href="http://archer.ac.uk/training/course-material/2014/08/CP2K/Slides/NSCCS_2014_Sanliang_Ling.pdf" style="font-family:Arial,Helvetica,sans-serif;font-size:13px" target="_blank" class="">http://archer.ac.uk/training/course-material/2014/08/CP2K/Slides/NSCCS_2014_Sanliang_Ling.pdf</a><span style="color:rgb(34,34,34);font-family:Arial,Helvetica,sans-serif;font-size:13px" class=""> )</span></p><p class="gmail-m_457967070791147240p1"><span style="color:rgb(34,34,34);font-family:Arial,Helvetica,sans-serif;font-size:13px" class=""><br class=""></span></p><p class="gmail-m_457967070791147240p1"><span style="color:rgb(34,34,34);font-family:Arial,Helvetica,sans-serif;font-size:13px" class="">and similar even in the regtests.</span></p><p class="gmail-m_457967070791147240p1"><span style="color:rgb(34,34,34);font-family:Arial,Helvetica,sans-serif;font-size:13px" class="">This seems to call both </span>XC_GGA_X_RPW86 and VWN (lda) .</p><p class="gmail-m_457967070791147240p1">What does it mean to call two different xc potentials in cp2k, and why would one do it?</p><p class="gmail-m_457967070791147240p1">Tracing this in the code is not so easy either.</p><p class="gmail-m_457967070791147240p1"><br class=""></p><p class="gmail-m_457967070791147240p1">Ins eom docs it seems to suggest that LDA only should be used for correlation. Is that because GGA is included </p><p class="gmail-m_457967070791147240p1">in the vDW function in cp2k? Also there is the question which pseudopotentials to use.</p><p class="gmail-m_457967070791147240p1"><br class=""></p><p class="gmail-m_457967070791147240p1">It seems to me that XC_GGA_X_RPW86 is correct without the VMN, and this would correapond to what is done for DF2 in</p><p class="gmail-m_457967070791147240p1">quantum espresso it seems.</p><p class="gmail-m_457967070791147240p1"><br class=""></p><p class="gmail-m_457967070791147240p1">I note some related problems with the documentation. I find examples like for LIBXC:</p><div class=""> <br class="webkit-block-placeholder"></div><p class="gmail-m_457967070791147240p1">FUNCTIONAL XC_MGGA_X_M06_L XC_MGGA_C_M06_L</p><p class="gmail-m_457967070791147240p1">but the docs at </p><p class="gmail-m_457967070791147240p1">say</p><ul class="gmail-m_457967070791147240disc" style="font-family: Times; font-size: inherit;"><li style="margin-left:0em;padding-left:0em;text-indent:0em" class=""><a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/XC/XC_FUNCTIONAL/LIBXC.html#list_FUNCTIONAL" id="gmail-m_457967070791147240FUNCTIONAL" target="_blank" class="">FUNCTIONAL</a></li><li style="margin-left:0em;padding-left:0em;text-indent:0em" class="">This keyword cannot be repeated and it expects precisely one word.<br class=""></li></ul><p class="gmail-m_457967070791147240p1">Yet there are "two words" </p><p class="gmail-m_457967070791147240p1">and including VWN afterwards also seems like more than one functional, as I said above.</p><p class="gmail-m_457967070791147240p1"><br class=""></p><p class="gmail-m_457967070791147240p1">Thank you!</p><p class="gmail-m_457967070791147240p1"><br class=""></p><p class="gmail-m_457967070791147240p1">Ronald Cohen</p><p class="gmail-m_457967070791147240p1"><br class=""></p><div class=""><span style="color:rgb(34,34,34);font-family:Arial,Helvetica,sans-serif;font-size:13px" class=""> </span><br class="webkit-block-placeholder"></div> <div class=""></div></div></div><div class=""><br class="webkit-block-placeholder"></div> -- <br class=""> You received this message because you are subscribed to the Google Groups "cp2k" group.<br class=""> To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+uns...@googlegroups.com" target="_blank" class="">cp2k+uns...@googlegroups.com</a>.<br class=""> To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank" class="">cp...@googlegroups.com</a>.<br class=""> Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" class="">https://groups.google.com/group/cp2k</a>.<br class=""> For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" class="">https://groups.google.com/d/optout</a>.<br class=""> </blockquote></div><div class=""><br class="webkit-block-placeholder"></div> -- <br class=""> You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.<br class=""> To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/_1FupbyRLPA/unsubscribe" class="">https://groups.google.com/d/topic/cp2k/_1FupbyRLPA/unsubscribe</a>.<br class=""> To unsubscribe from this group and all its topics, send an email to <a href="mailto:cp2k+uns...@googlegroups.com" class="">cp2k+uns...@googlegroups.com</a>.<br class=""> To post to this group, send email to <a href="mailto:cp...@googlegroups.com" class="">cp...@googlegroups.com</a>.<br class=""> Visit this group at <a href="https://groups.google.com/group/cp2k" class="">https://groups.google.com/group/cp2k</a>.<br class=""> For more options, visit <a href="https://groups.google.com/d/optout" class="">https://groups.google.com/d/optout</a>.<br class=""> </div></blockquote></div><br class=""></div></body></html>