In our previous work, a perturbed hard-trimersphere equation of state (PHTS EOS) was developed for modeling the phase equilibria of pure ionic liquids (ILs) (M.M. Alavianmehr et al., Ionics 22 (2016) 2447-2459. In this work, we have successfully extended the model to the mixtures of IL + IL and IL + solvent. Two temperaturedependent parameters appearing in the EOS are correlated with two microscopic scaling constants σ, the effective hardsphere diameter, and ε, the non-bonded interaction energy. The overall average absolute deviation (AAD) of the estimated densities from the literature data using the proposed model with and without non-additivity parameter (λ ij ) was found to be 0.44 and 0.79%, respectively. A modified Enskog equation and rough hard-sphere (RHS) theory are combined with our proposed equation of state to calculate the viscosity coefficient of ionic liquids and their mixtures. Finally, from the results obtained, a linear relation between logarithm of surface tension and viscosity property of ionic liquid was developed.