Analytical potential functions have been constructed for Mg m+Ar n , m)1, 2, clusters by combining ab initio calculations with the electrostatic model of ion-induced dipole and induced dipole-induced dipole interactions. The most stable structures and dynamic properties are investigated by molecular dynamics techniques. For Mg + Ar 12 an icosahedral geometry is predicted. For Mg + Ar n clusters with n ) 7-12 two series of minima with small differences in energy were found; in one the magnesium ion stays inside the cluster (solvated) and in the other on the surface of the system. For n > 9, clusters with Mg on the surface are slightly more stable and have a distinguishable high vibrational frequency relative to those with the magnesium ion inside. We find no such trends for Mg 2+ Ar n clusters, in which the Mg cation is always solvated. The absolute minimum of Mg 2+ Ar 4 is a regular tetrahedron, and that of Mg 2+ Ar 6 , a regular octahedron. The stability of all these clusters is investigated by studying caloric curves, root mean square bond length fluctuations, radial distributions, and power spectra.