In such a situation, molecular dynamics (MD) simulations can serve as a very useful tool with an arbitrary spatial and temporal resolution, provided the underlying interactions are described with sufficient rigour and accuracy. For the present systems containing solvated electrons, a quantum description of the electronic structure is imperative. The first application of molecular dynamics to the EMT in alkali metal-liquid ammonia systems accomplished 30 years ago had to resort -due to limited methodology and computational resources available at the time -to semi-empirical methods for the description of the interactions of the excess electrons among themselves and with the rest of the system. This severe approximation strongly limited both the range of accessible electronic structure scenarios and the fidelity of the predictions. In this study, we overcome the above limitations by employing ab initio molecular dynamics (AIMD) with an explicit quantum mechanical description of all the excess and valence electrons in the system based on the density functional theory (DFT) and the one-body Green's function (GW ) approximation. The use of a state-of-the-art hybrid density functional for MD ensures that our trajectories are faithful, while the GW method provides accurate electronic structure information via a post-DFT description of Coulombic screening. This simulation methodology is suficiently robust to describe the dynamics and electronic structure of highly concentrated metallic solutions.