Gaussian wavefunctions in a semicontinuum model for solvated electrons