<div dir="ltr"><div>Dear CP2K users,</div><div><br></div><div>I am interested in using the LC-wPBE functional in CP2K with a different value of omega from the original functional.</div><div>However, my test results in CP2K turned out to be quite different to those obtained using another program when I changed the value of omega.</div><div><br></div><div>Using the default value of omega for LC-wPBE (0.4), I obtained SCF energies of water as:<br>CP2K:               -76.41700469511423<br>Other program: -76.41700042778284</div><div><span style="display: inline !important; float: none; background-color: rgb(255, 255, 255); color: rgb(34, 34, 34); font-family: "Arial","Helvetica",sans-serif; font-size: 13px; font-style: normal; font-variant: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: left; text-decoration: none; text-indent: 0px; text-transform: none; -webkit-text-stroke-width: 0px; white-space: normal; word-spacing: 0px;">(Basis set: cc-pVTZ, taken from the Basis Set Exchange for CP2K)</span><b></b><i></i><u></u><strike></strike><br></div><div><b></b><i></i><u></u><strike></strike><br></div><div>That's OK. But when I changed the value of omega to 0.3 for both, my water energies differ significantly:<br>CP2K:               -75.96879714659988<br>Other program: -76.40681085213058</div><div><br></div><div>Then, I tried to reproduce some of the results from this paper (pubs.acs.org/doi/10.1021/acs.jctc.5b00820) using CP2K (LC-wPBE, omega=0.309).</div><div><br></div><div>HOMO/LUMO energies of uracil (eV):<br>From the paper: -9.62 and 0.51<br>CP2K:                -8.79 and 1.34<br>Other program:  -9.51 and 0.67<br>(I used the cc-pVTZ basis set in CP2K and the other program again.)</div><div><br></div><div>Hence, I think I could be doing something wrong with the omega parameters in CP2K?</div><div><br></div><div>More details:<br>I compiled CP2K 6.1.0 using the toolchain; that uses Libxc version 4.0.3. <br><span style="display: inline !important; float: none; background-color: rgb(255, 255, 255); color: rgb(34, 34, 34); font-family: "Arial","Helvetica",sans-serif; font-size: 13px; font-style: normal; font-variant: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: left; text-decoration: none; text-indent: 0px; text-transform: none; -webkit-text-stroke-width: 0px; white-space: normal; word-spacing: 0px;">All regtests passed without any errors after compilation.</span><br><br></div><div>This is the relevant section of my CP2K input file under DFT:</div><div><br></div><div>    &XC<br>      &XC_FUNCTIONAL<br>        &LIBXC<br>          FUNCTIONAL HYB_GGA_XC_LC_WPBE<br>        &END LIBXC<br>      &END XC_FUNCTIONAL</div><div>      &HF<br>        &INTERACTION_POTENTIAL<br>          POTENTIAL_TYPE MIX_CL<br>          OMEGA 0.4000<br>          SCALE_COULOMB 0.0<br>          SCALE_LONGRANGE 1.0<br>        &END<br>        FRACTION 1.0<br>      &END HF<br>    &END XC</div><div><br></div><div>My input files are attached as well.</div><div><br></div><div>Thanks for any help or suggestions!</div><div><br></div><div>Best,</div><div>Zi Cheng</div></div>