<div dir="ltr"><div>Dear Frederick,</div><div><br></div><div>Thanks very much for the explanation and the suggestions.</div><div><br></div><div>I've tried both of your suggested methods.</div><div><br></div><div>The first method, reparametrising the HSE functional in Libxc, crashed CP2K.<br>I did a bit of testing, and apparently setting FUNCTIONAL HYB_GGA_XC_HSE03 seems to be the cause of it.</div><div><br></div><div>Fortunately, the second method worked out well.<br>I was able to reproduce my test calculations from the other program, as well as those in the paper, so that solves my problem.</div><div><br></div><div>Thanks again!<br>Zi Cheng<br></div><div><br></div><br>On Tuesday, January 7, 2020 at 5:57:17 PM UTC+8, Frederick Stein wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Zi,<div><br></div><div>You have only changed the omega value in the HF section. But you also need a different omega value for the Libxc functional. These omega values are independent from each other.</div><div><br></div><div>I had a look at the source code. There is no way to change any parameters of the LC-wPBE functional in Libxc. I think you can reparametrize the HSE functional of Libxc</div><div><br></div><div>      &XC_FUNCTIONAL<br>        &LIBXC<br>          FUNCTIONAL HYB_GGA_XC_HSE03</div><div># According to <a onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fgitlab.com%2Flibxc%2Flibxc%2Fblob%2F4.0.3%2Fsrc%2Fhyb_gga_xc_hse.c\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHy4c_0mzIvlblOcQx8kYC9MZYnCA';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fgitlab.com%2Flibxc%2Flibxc%2Fblob%2F4.0.3%2Fsrc%2Fhyb_gga_xc_hse.c\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHy4c_0mzIvlblOcQx8kYC9MZYnCA';return true;" href="https://gitlab.com/libxc/libxc/blob/4.0.3/src/hyb_gga_xc_hse.c" target="_blank" rel="nofollow">https://gitlab.com/libxc/<wbr>libxc/blob/4.0.3/src/hyb_gga_<wbr>xc_hse.c</a> the parameters should be (beta, omega_HF, omega_PBE)</div><div>          PARAMETERS 0.0 0.309 0.309<br>        &END LIBXC<br>      &END XC_FUNCTIONAL<br><br></div><div>Or you might try it with the CP2K implementation of WPBE:</div><div><br>      &XC_FUNCTIONAL</div><div>        &XWPBE<br>          OMEGA 0.309</div><div>          SCALE_X   1.0</div><div>          SCALE_X0 0.0<br>        &END XWPBE</div><div>        &PBE</div><div>          SCALE_X   0.0</div><div>          SCALE_C   1.0<br>        &END PBE</div><div>      &END XC_FUNCTIONAL<br></div><div><br></div><div>I hope I could help you.</div><div><br></div><div>Best,</div><div><br></div><div>Frederick</div></div></blockquote></div>