<div dir="ltr">Hi Fernan,<br><blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">I would like to know how CP2k models an excess electron placed in a neutral system with N electrons by writing CHARGE -1.0 in the section CP2K_INPUT / FORCE_EVAL / DFT. Is this excess electron placed in the LUMO N-orbital?<br><br></div></blockquote><div>this isn't really how things work. The KS equations are solved, then filled up obeying the Aufbau principle. If the excess charge doesn't disturb the existing charge density, then it would be in an orbital very similar to the previous LUMO.</div><div> </div><blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">In addition, does cp2k correct the eigenvalues to the vacuum level? Or do they need to be corrected after the simulation is completed? Is there any reference?<br></div></blockquote><div><br></div><div>this depends on you boundary conditions (POISSON solver). In a periodic calculation there is no vacuum, or anything outside the system, and the energy levels are arbitrary. In a surface or molecular system (without dipole) the WAVELET solver will give the levels relative to zero.</div><div><br></div><div>Matt </div><blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><br>Best regards,<br> - Fernan Saiz, PhD<br>Department of chemistry<br>Imperial College London<br></div></blockquote></div>