The integrated sorption and diffusion (ISD) model has been developed to quantify radionuclide transport in compacted bentonite. The current ISD model, based on averaged pore aperture and the Gouy-Chapman electric double layer (EDL) theory can quantitatively account for diffusion of monovalent cations and anions under a wide range of conditions (e.g., salinity, bentonite density). To improve the applicability of the current ISD model for multivalent ions and complex species, the excluded volume effect and the dielectric saturation effect were incorporated into the current model, and the modified Poisson-Boltzmann equations were numerically solved. These modified models had little effect on the calculation of effective diffusivity of Sr2+/Cs+/I−. On the other hand, the model, modified considering the effective electric charge of hydrated ions, calculated using the Gibbs free energy of hydration, agreed well with the diffusion data including those of Sr2+.