All of our ab initio molecular dynamics (AIMD) simulations utilized the QUICKSTEP module in the CP2K software package with the unrestricted Kohn-Sham formalism to capture the dynamics of the hydrated electron. All of our NVT simulations used the Nose-Hoover chain thermostat at 300 K, and the AIMD equations of motions were propagated with a 0.5 fs time step. GTH pseudopotentials and the TZV2P basis set were employed to describe the atomic core and valence electrons, respectively. All of our simulations (both hydrated and excess electron calculations) utilized a self-interaction corrected PBE exchange-correlation functional with Grimme's D3 54 dispersion corrections. Specifically, we incorporated a self-interaction correction (SIC) for all orbitals in our AIMD simulations using the average density SIC, as implemented in the CP2K package. Following earlier studies, we set a = 0.2 and b = 0.25 for the scaling parameters in the SIC exchange-correlation functional, which was previously shown to give accurate electronic properties of PFOA/PFOS in water. For the auxiliary plane-wave (PW) basis used in the Gaussian and Plane-wave (GPW) scheme in CP2K, 600 and 60 Ry was used for the PW energy and reference grid cutoffs, respectively.