<div dir="ltr">Dear Frederick, <div>Thanks a lot for your clear answers. Indeed, they help a lot!</div><div><br></div><div>Regarding question 1 above, now it makes sense. However, my confusion originated after looking at the PURE SUBROUTINE wigner_slater_functional in src/atoms_util.F:</div><div><br></div><div><b><i> PURE SUBROUTINE wigner_slater_functional(rho, vxc)<br> REAL(KIND=dp), DIMENSION(:), INTENT(IN) :: rho<br> REAL(KIND=dp), DIMENSION(:), INTENT(OUT) :: vxc<br><br> INTEGER :: i<br> REAL(KIND=dp) :: ec, ex, rs, vc, vx<br><br> vxc = 0._dp<br> DO i = 1, SIZE(rho)<br> IF (rho(i) > 1.e-20_dp) THEN<br> ! 3/4 * (3/pi)^{1/3} == 0.7385588<br> ex = -0.7385588_dp*rho(i)**0.333333333_dp<br> vx = 1.333333333_dp*ex<br> rs = (3._dp/fourpi/rho(i))**0.333333333_dp<br> ec = -0.88_dp/(rs + 7.8_dp)<br> vc = ec*(1._dp + rs/(3._dp*(rs + 7.8_dp)))<br> vxc(i) = vx + vc<br> END IF<br> END DO<br><br> END SUBROUTINE wigner_slater_functional</i><br></b></div><div><i><br></i></div><div>You can see this subroutine computes both the correlation energy (ec, vc) as well as the exchange energy (ex,vx), which contradicted the flag<b> LDA_C_WIGNER</b> to compute <b>only the correlation bit.</b> Now your explanation helps to unpuzzle my concern. Since there was no date recorded in the description of each sections in the online manual, I could not even guess when these flags were introduced.</div><div><br></div><div>I shall continue investigating further. Thanks a lot for your help.<br></div><div>Best wishes</div><div><br></div><div>Ivan</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, 26 Apr 2021 at 13:25, 'Frederick Stein' via cp2k <<a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>> wrote:<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>Dear Ivan,</div><div><br></div><div>Short answer:<br></div><div>All section names you are referring to are actually functionals from LibXC. That is why you do not find the respective code in CP2K.</div><div>A pseudo potential optimized for a functional of a given class usually provides reasonable results for other functionals of the same class. For metals, you should check that separately.<br></div><div><br></div><div>Long answer:<br></div><div>To 1+2) Do you mean by accident the master branch instead of 8.1 ? Your given functional names are functional names from LibXC. In the recent trunk, we have added separate sections for the functionals from LibXC. That is why you will not find the code in CP2K but in LibXC and why these sections do not work with older versions of CP2K. Some more information about the the naming scheme of LibXC (also consult the manual of LibXC):</div><div>- Names starting with "HYB_" refer to hybrid functionals for which you have to set up a separate Hartree-Fock section on your own.</div><div>- Then, you will find the class of the underlying DFT functional ("LDA", "GGA", "MGGA").</div><div>- Then, the type of functional ("K" for kinetic energy, "X" for exchange-only, "C" for correlation-only, "XC" for exchange-correlation functional).</div><div>- Finally, the functional name (have a look at the manual of LibXC for you LibXC distribution in-use).<br></div><div><br></div><div>
There is one functional which sticks to the above convention but is not provided by LibXC but by CP2K (LDA_HOLE_T_C_LR).</div><div>Some functionals from LibXC are natively supported by CP2K like LYP (in LibXC: GGA_C_LYP) or PZ81 (in LibXC: LDA_C_PZ).<br></div><div><br></div><div>To 3+4)</div><div>Basically yes, but you should always check the suitability of the pseudopotentials, especially for metals.</div><div><br></div><div>To 5)</div><div>Probably yes. In case of doubt, check the references in the LibXC manual (I do not know these functionals)</div><div><br></div><div>To 6)</div><div>As stated above, "_C_" tells you that the corresponding functional is just a correlation functional and you have to add a suitable exchange functional (or Hartree-Fock). For LDA-based functionals, it is LDA_X (from LibXC) or XALPHA (from CP2K).</div><div><br></div><div>I hope it helps you. If you have more questions, just ask.</div><div><br></div><div>Best,</div><div>Frederick<br></div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr"><a href="mailto:iva...@gmail.com" target="_blank">iva...@gmail.com</a> schrieb am Montag, 26. April 2021 um 13:10:26 UTC+2:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br><div><div>Dear CP2K community,</div><div><br></div><div>I attach the test.png file where I have tested different XC functionals of my interest in a table format. I have tried CP2K versions from 6.1 to 8.1. Questions:</div><div><br></div><div>1)<span style="white-space:pre-wrap"> </span>Functional with directives LDA_C_HL, LDA_C_PZ, LDA_C_WIGNER, GGA_C_AM05, GGA_X_AM05 and GGA_X_RPBE do not seem to be recognised for version 7.1 and 6.1, but they work well for 8.1. Is this correct? If not, what I am doing wrong? The only two XC subroutines I have found in the src folder for version 8.1 related to the above functionals is wigner_slater_functional (for WIGNER defined in src/atom_utils.F, created in 2008) and (for PZ, in src/xc/xc_perdew_zunger.F, created in 2004). So I wonder why LDA_C_WIGNER nor LDA_C_PZ are recognised for versions 6.1 and 7.1.</div><div>2)<span style="white-space:pre-wrap"> </span>This is related to 1), somehow. Apart from WIGNER, PZ, PBE, REVPBE and PBESOL, I cannot find the subroutines that correspond to the other functionals using <b>grep</b>. Could you please instruct me how/where to find the relevant subroutines?</div><div>3)<span style="white-space:pre-wrap"> </span>Shall one use PADE potentials for HL, PZ and Wigner? </div><div>4)<span style="white-space:pre-wrap"> </span>Shall one use PBE potentials for AM05, RPEB, PBESOL?</div><div>5)<span style="white-space:pre-wrap"> </span>Regarding AM05: if I want to use this functional, shall I specify both GGA_C_AM05 and GGA_X_AM05 together as follows?</div><div><br></div><div>&XC</div><div> &XC_FUNCTIONAL</div><div> &GGA_X_AM05<br></div><div> &END GGA_X_AM05</div><div><div> &GGA_C_AM05<br></div><div> &END GGA_C_AM05</div></div><div> &END XC_FUNCTIONAL</div><div>&END XC</div><div><br></div><div>6)<span style="white-space:pre-wrap"> </span>For LDA_C_HL, “_C_” refers that the HL only proposed a correlation term, and used the Slater functional for the exchange part (from the paper Hedin and Lundquist wrote). Does LDA_C_HL (and the rest LDA_C_ of the table) include the exchange part already? Or one also need to include directive LDA_X (Slater exchange) in the input file? I mean, shall LDA_X block below be defined or not?</div><div><div><br></div><div>&XC</div><div> &XC_FUNCTIONAL</div><div><b style="background-color:red"> &LDA_X </b></div><div><b style="background-color:red"> &END LDA_X</b></div><div><div> &LDA_C_HL<br></div><div> &END LDA_C_HL</div></div><div> &END XC_FUNCTIONAL</div><div>&END XC</div><div><br></div></div><div>Thanks a lot in advance for your help with these questions<br></div></div><div><br></div><div><br></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:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/d54e315d-468a-4e25-92fd-d4fd8950501cn%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/d54e315d-468a-4e25-92fd-d4fd8950501cn%40googlegroups.com</a>.<br>
</blockquote></div>