The ability to separate fission products by electrodeposition from molten salts depends, in part, on differences between the interactions of the different fission product cations with the ions present in the molten salt "solvent". These differences may be expressed as ratios of activity coefficients, which depend on the identity of the solvent and other factors. Here. we demonstrate the ability to calculate these activity coefficient ratios using molecular dynamics simulations with sufficient precision to guide the choice of suitable solvent systems in practical applications. We use polarizable ion interaction potentials which have previously been shown to give excellent agreement with structural, transport, and spectroscopic information of the molten salts, and the activity coefficients calculated in this work agree well with experimental data. The activity coefficients are shown to vary systematically with cation size for a a set of trivalent cations.