A polarizable force field for the sulfate anion SO$_4^{2–}$ has been developed and extended from nonpolarizable force fields in order to reproduce its structural and thermodynamics properties in aqueous solution. Two force fields with different atomic partial charges on S and O have been tested and used with molecular dynamics with explicit polarization. The results obtained with our developed force field are in good agreement with the experimental hydration properties of the sulfate anion. In addition to molecular dynamics simulations of the sulfate anion in aqueous solution, potentials of mean force of sulfate electrolytes have been calculated via umbrella-sampling molecular dynamics simulations, i.e., MgSO$_4$, EuSO$_4$$^+$, and UO$_2$SO$_4$. These potentials allow for calculating pair association constants directly comparable to the experimental ones. In the case of monoatomic cations such as Mg$^{2+}$ and Eu$^{3+}$, the association constants calculated are in very good agreement with the experimental values, i.e., p$K_{calc}$ = 2.21 (vs 2.21 experimentally) and 3.86 (vs 3.56–3.78 experimentally) for MgSO$_4$ and EuSO$_4$$^+$, respectively. In the case of purely molecular electrolyte (UO$_2$SO$_4$), the association constant calculated (p$K_{calc}$ = 1.58–2.07) is in agreement with the range of values available in the literature (p$K_{exp}$ = 1.17–3.14).