No one can ever deny the significance of calculations of solubilities of industrial solid compounds in supercritical CO2 in separation processes. In this work, the Peng-Robinson (PR) and the Esmaeilzadeh-Roshanfekr (ER) equations of state (EoS) along with several mixing rules including the Wong-Sandler (WS), the covolume dependent (CVD) and the van der Waals one (VDW1) and two (VDW2) fluid mixing rules are applied to evaluate the solubilities of 52 mostly used solid compounds in supercritical carbon dioxide. Besides, the Van-Laar excess Gibbs energy (G(ex)) model is applied in phase behavior calculations by the WS mixing rule. The optimal values of the proposed thermodynamic model parameters are evaluated using the DE (differential evolution) optimization strategy. The absolute average deviations of the model results from 1776 experimental data points and the optimal values of the adjustable parameters of the model are reported to investigate the capabilities of combinations of each equation of state with different mixing rules in calculations of the solubilities. The results indicate that the combination of the ER EoS with the WS mixing rule leads to more accurate results (MD = 9.0%) compared with other ones. (C) 2010 Elsevier B.V. All rights reserved.