<div dir="ltr">After computing the total electronic energy for a neon atom using wB97X-V and comparing it to QChem reference results, it seems that the DFT XC energy in CP2K is not being computed correctly. Does anyone have any idea what might be wrong here? Perhaps the parametrization?<div><br><div>CP2K Results:</div><div><div class="im" style="color: rgb(80, 0, 80);"><div>  Overlap energy of the core charge distribution:               0.00000000000000</div><div>  Self energy of the core charge distribution:               -148.47094303888323</div></div><div>  Core Hamiltonian energy:                                      0.61509680993492</div><div>  Hartree energy:                                              14.01686743627672</div><div>  Exchange-correlation energy:                                 -4.11859163513813</div><div>  Hartree-Fock Exchange energy:                                -3.38266780588875</div><div>  Dispersion energy:                                            0.02876776737550</div><div><br></div><div>  GAPW| Exc from hard and soft atomic rho1:                    -6.47114106845743</div><div>  GAPW| local Eh = 1 center integrals:                         17.39257016010844</div><div><br></div><div>  Total energy:                                              -130.39004137467197</div><div><br></div><div>  outer SCF iter =    2 RMS gradient =   0.85E-08 energy =       -130.3900413747</div><div>  outer SCF loop converged in   2 iterations or  435 steps</div><div><br></div><div><br></div><div> ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):             -130.390041374671966</div></div><div><br></div><div>QChem results (aug-cc-pVQZ/250,590):</div><div><br></div><div><div> Nonlocal correlation =     0.0441808008</div><div>    7    -128.9329200166      6.16E-09 Convergence criterion met</div><div> -----------------------------<wbr>----------</div><div> One-Electron    Energy =  -182.3744736608</div><div> Total Coulomb   Energy =    65.8986384358</div><div> Alpha Exchange  Energy =    -1.6882204462</div><div> Beta  Exchange  Energy =    -1.6882204462</div><div> DFT   Exchange  Energy =    -8.<a href="tel:6509462294" value="+16509462294" target="_blank" style="color: rgb(17, 85, 204);">6509462294</a></div><div> DFT Correlation Energy =    -0.4296976698</div><div> Nuclear Repu.   Energy =     0.0000000000</div><div> Nuclear Attr.   Energy =  -311.0040993334</div><div> Kinetic         Energy =   128.6296256727</div></div><div><br><div><br></div><div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME /global/u1/a/amelcher/soft/cp2k/cp2k/tests/QS/EMSL_BASIS_SETS</div><div>    POTENTIAL_FILE_NAME /global/u1/a/amelcher/soft/cp2k/cp2k/tests/QS/POTENTIAL</div><div>    &MGRID</div><div>      CUTOFF 1000</div><div>      NGRIDS 5</div><div>      REL_CUTOFF 50</div><div>    &END MGRID</div><div>    &QS</div><div>      METHOD GAPW</div><div>      EPS_PGF_ORB 1.0E-12</div><div>      EPS_FILTER_MATRIX 0.0e0</div><div>    &END QS</div><div>    &POISSON</div><div>      PERIODIC NONE</div><div>      PSOLVER MT</div><div>    &END</div><div>    &SCF</div><div>      EPS_SCF 1.0E-8</div><div>      SCF_GUESS ATOMIC</div><div>      MAX_SCF 400</div><div>      &OUTER_SCF</div><div>        EPS_SCF 1.0E-8</div><div>        MAX_SCF 50</div><div>      &END OUTER_SCF</div><div>      &OT ON</div><div>      &END</div><div>    &END SCF</div><div>    &XC</div><div>      &XC_FUNCTIONAL</div><div>        &BECKE97</div><div>            PARAMETRIZATION  wB97X-V</div><div>            SCALE_X 1.0</div><div>            SCALE_C 1.0</div><div>        &END</div><div>      &END XC_FUNCTIONAL</div><div>      &HF</div><div>        &SCREENING</div><div>           EPS_SCHWARZ 1.0E-14</div><div>        &END</div><div>        &INTERACTION_POTENTIAL</div><div>           POTENTIAL_TYPE MIX_CL</div><div>           SCALE_COULOMB   0.167</div><div>           SCALE_LONGRANGE 0.833</div><div>           OMEGA 0.30</div><div>        &END</div><div>        &MEMORY</div><div>          MAX_MEMORY 256</div><div>        &END</div><div>      &END</div><div>      &vdW_POTENTIAL</div><div>         DISPERSION_FUNCTIONAL NON_LOCAL</div><div>         &NON_LOCAL</div><div>           TYPE RVV10</div><div>           PARAMETERS 6.3 0.0093</div><div>           VERBOSE_OUTPUT</div><div>           KERNEL_FILE_NAME /global/u1/a/amelcher/soft/cp2k/cp2k/tests/QS/rVV10_kernel_table.dat</div><div>           CUTOFF  150</div><div>         &END NON_LOCAL</div><div>      &END vdW_POTENTIAL</div><div>   &END XC</div><div>  &END DFT</div><div>  &SUBSYS</div><div>    &CELL</div><div>      ABC 10.0 10.0 10.0</div><div>      PERIODIC NONE</div><div>    &END CELL</div><div>    &COORD</div><div>Ne       0.0000 0.0000 0.0000</div><div>    &END COORD</div><div>    &KIND Ne</div><div>      BASIS_SET aug-cc-pVQZ</div><div>      LEBEDEV_GRID 590</div><div>      RADIAL_GRID 250</div><div>      POTENTIAL ALL</div><div>    &END KIND</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div>&GLOBAL</div><div>  PROJECT wb97xv25</div><div>  PRINT_LEVEL LOW</div><div>  RUN_TYPE ENERGY</div><div>  &TIMINGS</div><div>    THRESHOLD 0.000000001</div><div>  &END</div><div>&END GLOBAL</div></div><div><br></div></div></div></div>