We developed a simple pair-additive Lennard-Jones plus Coulomb potential for molecular simulations of the trivalent cation Al(3+) in water which accounts reasonably well for the behavior of aluminum aqueous solutions. The model predicts an octahedral first hydration shell containing 6 water molecules and a trigonal second shell with 12 molecules on average, in good agreement with the available experimentally determined structure. The peak positions of the cation-oxygen radial distribution function are only slightly compressed compared to the x-ray structure, the hydration enthalpy is 10% too low, and the cation self-diffusion coefficient and the single-particle second rank reorientational time are in excellent agreement with inelastic neutron scattering and NMR spectroscopy data, respectively. The model also captures the essential vibrational features of the hydrated [Al(H(2)O)(6)](3+) complex. It predicts the main O-Al-O bending mode frequency to within approximately 5%, but significantly overestimates the frequency of the totally symmetric Al-O stretching mode. Overall, the accuracy of the proposed model is as good as the best available classical potentials, if not better in some aspects, with a much simpler functional form, which makes it an attractive alternative for computer simulations of Al(3+) in more complex aqueous and biomolecular systems.