A microscopic calculation of the solvation dynamics of an ion in liquid water is presented. The calculated solvation time correlation function shows an ultrafast Gaussian decay which carries about 70%-90% of the strength followed by a biexponential decay with time constants equal to 250 fs and 1 ps. These results are in excellent agreement with the computer simulations of Maroncelli and Fleming and also with the experimental findings of Barbara and Jarzeba. In addition, we find that both the rotational librations and the intermolecular translational vibrational modes of water contribute significantly to the initial Gaussian decay.