In this study, we utilised molecular dynamics simulations to investigate the electronic properties of ammonia solutions containing varying concentrations of solvated dielectrons and lithium cations. In order to generate initial configurations we employed force-field molecular dynamics (FFMD) simulations using a liquid-adapted ammonia force field. These configurations then served as the starting point for production runs of AIMD utilising the VASP program package. Within AIMD simulations, the electronic structure was described using the hybrid revPBE0-D3 functional. For each system, multiple trajectories were generated, involving an equilibration phase followed by a production run to provide adequate sampling (for more details see SI).