Molecular dynamics (MD) simulations based on ab initio evaluated two- and three-body potentials were performed for Ti3+, Cr3+, and Co3+ ions in water. The ions' hydration structure was evaluated in terms of radial distribution functions, coordination numbers, and angular distributions. The first solvation shell shows an exact coordination number of 6 for, all ions and average M3+-O distances of 2.08, 2.05, and 2.03 Angstrom for Ti3+, Cr3+, and Co3+, respectively. The structural parameters obtained after inclusion of three-body correction terms are in good agreement with experimental values. The energies of hydration obtained from three-body corrected molecular dynamics simulations are also in good agreement with experiment, except for Till where the classical simulations are failing to reproduce the Jahn-Teller effect.