Variational calculations employing explicitly correlated Gaussian functions have been performed for the ground states of 9Be and 9Be+ including the nuclear motion i.e., without assuming the Born-Oppenheimer BO ...