Development and validation of an integrated computational approach for the study of ionic species in solution by means of effective two-body potentials. The case of Zn2+, Ni2+, and Co2+ in aqueous solutions