Having verified that our simulation protocol gives a good description of liquid water's structure, we turn next to the simulated properties of e - (aq) . Fig. shows the RDF for e - (aq) based on e --oxygen and e --hydrogen distances (panels (a) and (b) respectively). Encouragingly, QM/MM separation is maintained as shown by the red and blue filled regions of Fig. . Comparing to Fig. , we note that there appears to be more overlap in the distributions of QM e --O and MM e --O distances compared to the respective QM and MM distributions of O-O distances relative to a central fixed oxygen atom in the pure water simulations. This arises from diffusion of the excess electron's centroid position away from the center of the QM region (discussed above in Section III), rather than a breakdown of QM/MM separation. In Fig. we see the presence of a clear cavity in which the electron resides, with an e --O turn-on distance of 1.7 Å, and two strong peaks in the e --O RDF at 2.5 Å and 4.3 Å, indicating the first and second solvation shells. The e --H RDF displays a single strong peak at 1.5 Å, consistent with a predominantly O-H bond aligned solvation of the excess electron. Integrating to the first minimum of the e --O RDF at 3.3 Å provides a coordination number of 4.0, which can be compared to the same procedure for the first minimum of the e -H RDF at 2.4 Å which yields a coordination number of 3.9. The structure of e - (aq) is thus essentially the same as found in previous QM/MM simulations at the self-interaction corrected BLYP-D level, and consistent with the overall picture provided by cavity-forming MQC simulations, and recent PIMD simulations with ML potentials, although the latter's prediction of quasi-stable dual cavity structures is not observed here, perhaps due to the neglect of NQE.