We built a spherical nanodroplet model of the condensed-phase e - (aq) following a similar protocol to the companion paper. A cavity for the excess electron was carved out by placing a chloride ion at the origin of a 25-Å radius droplet of 2017 water molecules created using the SolvateCap command of tleap in the Amber 2018 package. The water forcefield was described by the SPC/Fw model, and chloride Lennard-Jones parameters were taken from Ref. 51. A half-harmonic cap potential with a force constant of 10 kcal/mol/ Å2 prevented evaporation of water at the droplet/vacuum interface, and the chloride ion was restrained to the origin using a harmonic potential with a force constant of 1000 kcal/mol/ Å2 . Following minimization, heating, and equi-libration, a 100-ps NVT production run was performed to generate ten independent snapshots separated by 10 ps, from which the chloride ion was then removed to generate initial conditions for the e - (aq) simulations. To model e - (aq) , 32 water molecules nearest the origin of the droplet were treated at a QM level, with the remaining 1985 water molecules treated at an MM level. Since the switch from MM to QM introduced large forces on the QM water molecules, we quenched the largest components of their forces by performing 5 steps of minimization using the L-BFGS algorithm, implemented in DL-FIND. During dynamics, QM/MM separation was maintained with FlexiBLE implemented in a development version of Terachem, and interfaced with OpenMM 7.3. We used a FlexiBLE exponent parameter of α = 15 Å-1 and a convergence parameter of γ = 10 -3 k B T . The QM region was described at the BH&HLYP-D3/6-31++G* level of theory, where D3 represents dispersion corrections with Becke-Johnson damping. The QM and MM subsystems were coupled with the standard approach of elecrostatically embedding the QM particles in the field of the MM point charges. Mechanical embedding forces were also included via Lennard-Jones pair potentials between the QM and MM particles according to the SPC/Fw forcefield. AIMD was propagated with the velocity Verlet algorithm 59 using the Niklasson time-reversible Born-Oppenheimer integrator with a 0.5-fs timestep. Temperature was maintained at 298 K with the Bussi-Parrinello thermostat. We generated ten trajectories of 7 ps in length. The first 2 ps of data of each trajectory was discarded as an equilibration period.