<div dir="ltr"><div>Dear CP2K experts,</div><div><br></div><div>Running a charged calculation with the self-interaction correction by Mauri et al. (as defined in CP2K_INPUT / FORCE_EVAL / DFT / SIC) leads to a large <i>negative</i> Hartree energy contribution -- something that doesn't look physically possible:</div><div><br></div><div>==========================================================<br></div><div>  *** SCF run converged in    37 steps ***<br><br><br>  Electronic density on regular grids:      -1054.9999988859        0.0000011141<br>  Core density on regular grids:             1055.9999999431       -0.0000000569<br>  Total charge density on r-space grids:        1.0000010572<br>  Total charge density g-space grids:           1.0000010572<br><br>  Overlap energy of the core charge distribution:               0.00001597916987<br>  Self energy of the core charge distribution:              -7659.99497268730556<br>  Core Hamiltonian energy:                                   3104.15741917188643<br><b>  Hartree energy:                                          -12743.39184120392383</b><br>  Exchange-correlation energy:                               -820.98172868636971<br><br>  Total energy:                                            -18120.21110742654128<br></div><div>==========================================================</div><div><br></div><div>After this section of the output file, a warning is printed ( <b>WARNING : did not converge in ot_eigensolver</b> ) and the HOMO-LUMO gap for this spin channel is nonsense (64 eV). Finally, a <b>CPASSERT</b> on <b>line 528</b> of the file <b>qs_vxc.F</b> aborts the calculation (see the .out file). I assume this is the consequence (or at least related to) the large negative Hartree energy in the first place.</div><div><br></div><div>We are running this using CP2K version 4.1, but the same happens with version 6.1. The input file, initial structure, and the full output file are attached to this email. The initial structure is sensible and the system behaves well under regular PBE exchange & correlation (no SIC, all other settings kept identical), as well as under the hybrid PBE0 functional.<br></div><div><br></div><div>Are we doing something wrong? I would appreciate any and all input.</div><div><br></div><div>Yours sincerely,</div><div>Patrick Gono<br></div></div>