Propensities in the solvation of M+-Benzene systems (MÂ =Â Na, K, Rb) investigated by cluster dynamics

► A semiempirical interaction potential is used to model intermolecular interactions. ► The model is based on the decomposition of the molecular polarizability. ► The model uses an improved Lennard Jones function. ► The solvation characteristics depend on the size of the cation.