In this research, the solubility of two single gases, hydrogen sulfide and carbon dioxide, in the sixteen ionic liquids ([emim][eFAP], [emim][EtSO4], [emim][OTf], [emim][Tf2N], [bmim][BF4], [bmim][OTf], [bmim][PF6], [bmim][Tf2N], [hmim][PF6], [hmim][Tf2N], [omim][PF6], [omim][Tf2N], [HOemim][BF4], [HOemim][OTf], [HOemim][PF6] and [HOemim][Tf2N]) are modeled using the Peng-Robinson and the PC-SAFT equations of state. For the ionic liquids, the pure component parameters of the PC-SAFT EoS have been obtained by the predictions of the model to both of experimental H2S and CO2 solubility data in ionic liquids as well as experimental liquid density data. Also, the critical properties and acentric factor of the ionic liquids were estimated using the modified Lydersen-Joback-Reid method. In the next step, the temperature-independent interaction parameters between the ionic liquids and the gases for both equations of state were calculated by fitting to the experimental VLE data and the obtained results were analyzed. In most of the studied binary systems, the percent average absolute deviation was lower than 10% that shows good agreement between the experimental and predicted values.