The remainder of the paper is as follows. In section II we benchmark DFT for electron-water interactions. In section III we describe the computational details of our condensed-phase model of e - (aq) . Results are presented in section IV, where we first apply FlexiBLE to a QM/MM treatment of neat water at the BH&HLYP-D3/6-31+G* level in section IV A and then to a QM/MM description of e - (aq) in section IV B. Dynamics of e - (aq) are explored in section IV C, and the NBO analysis is discussed in section IV D. Finally, conclusions are drawn in section V.