# Rigorous Connectionist Models to Predict Carbon Dioxide Solubility in Various Ionic Liquids

Hocine Ouaer, Amir Hossein Hosseini, Menad Nait Amar, Mohamed El Amine

Ben Seghier, Mohammed Abdelfetah Ghriga, Narjes Nabipour, Pål Østebø

Andersen, Amir Mosavi, Shahaboddin Shamshirband

► To cite this version:

Hocine Ouaer, Amir Hossein Hosseini, Menad Nait Amar, Mohamed El Amine Ben Seghier, Mohammed Abdelfetah Ghriga, et al.. Rigorous Connectionist Models to Predict Carbon Dioxide Solubility in Various Ionic Liquids. Applied Sciences, 2020, 10 (1), pp.304. 10.3390/app10010304. hal-02482008

HAL Id: hal-02482008

https://hal.science/hal-02482008

Submitted on 29 Oct 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Rigorous Connectionist Models to Predict Carbon Dioxide Solubility in Various Ionic Liquids

Hocine Ouaber ¹, Amir Hossein Hosseini ², Menad Nait Amar ³,*,
Mohamed El Amine Ben Seghier ⁴,⁵,*, Mohammed Abdelfetah Ghriga ¹,⁶,*
Narjes Nabipour ⁷,*, Pål Østebø Andersen ⁸,⁹,*, Amir Mosavi ¹⁰,¹¹,¹²,¹³ and
Shahaboddin Shamshirband ¹⁴,¹⁵,*,

¹ Laboratoire Génie Physique des Hydrocarbures, Faculté des Hydrocarbures et de la Chimie, Université M'Hamed Bougara de Boumerdes, Avenue de l'Indépendance, Boumerdes 35000, Algeria; ouaerhocine@univ-boumerdes.dz (H.O.); ghriga.ma@univ-pau.fr (M.A.G.)

2 Petroleum Department, Semnan University, Semnan 3513119111, Iran; s.amirhosein69@yahoo.com

³ Département Etudes Thermodynamiques, Division Laboratoires, Sonatrach, Boumerdes 35000, Algeria; m.naitamar@univ-boumerdes.dz

⁴ Division of Computational Mathematics and Engineering, Institute for Computational Science, Ton Duc Thang University, Ho Chi Minh City 729000, Vietnam; benseghier@tdtu.edu.vn

5 Faculty of Civil Engineering, Ton Duc Thang University, Ho Chi Minh City 729000, Vietnam

6 Institut des Sciences Analytiques et de Physico-chimie pour l'Environnement et les Matériaux, IPREM, UMR 5254, CNRS Université de Pau et des Pays de l'Adour/E2S, 2 avenue P. Angot, Technopôle Hélioparc, 64000 Pau, France

Institute of Research and Development, Duy Tan University, Da Nang 550000, Vietnam

8 Department of Energy Resources, University of Stavanger, 4036 Stavanger, Norway; pal.andersen@uis.no

9 The National IOR Centre of Norway, University of Stavanger, 4036 Stavanger, Norway

10 Institute of Structural Mechanics, Bauhaus University Weimar, D-99423 Weimar, Germany; amir.mosavi@uni-weimar.de

11 Institute of Automation, Kalman Kando Faculty of Electrical Engineering, Obuda University, 1034 Budapest, Hungary; amir.mosavi@kvk.uni-obuda.hu

12 Faculty of Health, Queensland University of Technology, 130 Victoria Park Road, Brisbane City, QLD 4059, Australia; a.mosavi@qut.edu.au

13 School of the Built Environment, Oxford Brookes University, Oxford OX3 0BP, UK; a.mosavi@brookes.ac.uk

14 Department for Management of Science and Technology Development, Ton Duc Thang University, Ho Chi Minh City, Vietnam

15 Faculty of Information Technology, Ton Duc Thang University, Ho Chi Minh City, Vietnam

* Correspondence: narjesnabipour@duytan.edu.vn (N.N.); Shahaboddin.shamshirband@tdtu.edu.vn (S.S.)

Received: 16 November 2019; Accepted: 23 December 2019; Published: 31 December 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Estimating the solubility of carbon dioxide in ionic liquids, using reliable models, is of paramount importance from both environmental and economic points of view. In this regard, the current research aims at evaluating the performance of two data-driven techniques, namely multilayer perceptron (MLP) and gene expression programming (GEP), for predicting the solubility of carbon dioxide (CO₂) in ionic liquids (ILs) as the function of pressure, temperature, and four thermodynamical parameters of the ionic liquid. To develop the above techniques, 744 experimental data points derived from the literature including 13 ILs were used (80% of the points for training and 20% for validation). Two backpropagation-based methods, namely Levenberg-Marquardt (LM) and Bayesian Regularization (BR), were applied to optimize the MLP algorithm. Various statistical and graphical assessments were applied to check the credibility of the developed techniques. The results were then compared with those calculated using Peng-Robinson (PR) or Soave-Redlich-Kwong (SRK) equations of state (EoS). The highest coefficient of determination (R² = 0.9965) and the lowest root mean square error (RMSE = 0.0116) were recorded for the MLP-LMA model on the full dataset (with a negligible difference to the MLP-BR model). The comparison of results from this model with
the vastly applied thermodynamic equation of state models revealed slightly better performance, but the EoS approaches also performed well with $R^2$ from 0.984 up to 0.996. Lastly, the newly established correlation based on the GEP model exhibited very satisfactory results with overall values of $R^2$ = 0.9896 and RMSE = 0.0201.

**Keywords:** CO₂ solubility; ionic liquids; carbon dioxide; multilayer perceptron; gene expression programming; prediction; equation of state; machine learning

# 1. Introduction

The natural gas produced from the subterranean gas fields and subsequently transported through pipelines should meet certain specifications such as environmental and safety standards as well as those of sale gas sectors. The products destined for sale should be free of undesirable contaminants, e.g., carbon dioxide (CO₂) and hydrogen sulfide (H₂S), which are both toxic and unfriendly from an environmental point of view. For instance, CO₂ is considered the main contributor to global warming and climate change. Several treatments are employed to remove acidic gases from natural gas. The most famous are the alkanolamine-based treatment operations [1–3]. This technique was firstly introduced for carbon dioxide removal in 1991 [4,5]. The used alkanolamines are organic compounds, such as monoethanolamine (MEA), diethanolamine (DEA), and triethanolamine (TEA). By far, MEA was the most preferred alkanolamine compared to DEA and TEA because of its reactivity, low molecular weight, and lower required circulation to maintain a given amine to acid-gas mole ratio.

The gas-liquid absorption in amine-based solvents is an efficient process in gas sweetening. Nevertheless, some imperfections have been observed, such as the creation of corrosive byproducts due to the amine degradation, water transfer to the gas stream during the desorption stage, and loss of the feedstock (amine), making the treatment operations expensive [6–10]. As an alternative, a new class of non-aqueous and environmentally friendly innovating fluids, known as ionic liquids (ILs), has emerged. ILs have many industrial applications such as catalysis for clean technology [11] and the removal of contaminants from refinery feedstock [12]. Furthermore, Ion Engineering Company is intended to use the know-how of ionic liquids for industrial-scale sweetening of natural gas and flue gas CO₂ separation [13,14], as stated by Hasib-ur-Rahman et al. [15].

Ionic liquids are molten salts, which are liquid (non-volatile) at room temperature. They are comprised exclusively of positively and negatively charged ions. Due to their bulky and asymmetrical cation structure, ILs have a low affinity to constitute crystals [16]. Manipulation of the cation and/or anion allows designing ILs adaptable to any particular application requirements [17]. Moreover, ILs are a perfect medium for acid gas solubilization over wide ranges of temperature and pressure. Thus, great attention was paid to evaluate the performance of ILs as a gas-cleaning agent in gas refinery plants [10]. The only and the most discussed disadvantage related to the use of ILs is their high viscosity, but it can be bypassed as the viscosity can be regulated over a reasonable range of about <50 cP to >10,000 cP by selecting an adequate mixture of cation and anion [15].

In the past few years, many experimental studies were conducted to estimate the solubility of acid gases in ILs [1, 6, 8, 10, 18, 19], especially the carbon dioxide solubility [20–24]. The obtained results confirmed the ILs to be very efficient in carbon dioxide removal. Unfortunately, the experimental studies require many laboratory tests, which are expensive, difficult, tedious, and time-consuming. As an alternative strategy, the solubility of acid gases in ILs has been modeled using thermodynamical laws and the equation of state (EoS). The thermodynamic laws used for modeling the solubility of acid gases in ILs can be divided into four groups that include cubic equations, quantum mechanics-based methods, activity coefficient methods, and statistical mechanics-based molecular approaches [25]. The most used models are the Peng–Robinson equation of state (EoS), the generic Van der Waals EoS, the generalized Redlich–Kwong cubic EoS, the law and extended law of Henry, and the equation of
Krichevsky-Kasarnovsky [8,26-32]. It was noticed that the models describe well the systems at low and moderate pressures [33], however, the equations of state suffer from many weaknesses. The equation of state can be reliable only for an individual system and not for more interestingly, multiple systems. They require various adjustable parameters, which should be optimized based on real data within a particular and limited range of thermodynamic conditions. Consequently, developing more general and powerful models to predict the solubility of acid gases, especially carbon dioxide, in ILs is of paramount importance.

Recently, many soft computing methods have been applied to model gas solubility and phase equilibrium. One of these methods is the artificial neural network (ANN), which represents an important embranchment computational intelligence method that can be used without any pre-assumption of the input-output relationship [34,35]. Multilayer perceptron (MLP), radial basis function network (RBF), multi-layer feed-forward network, and gene expression programming (GEP) are the general categories of ANN. Fuzzy logic (FL) is also one of the computational intelligence methods, which can model complicated nonlinear relations [36]. Likewise, the support vector machine (SVM) proposed by Vapnik was shown to be a very performant smart model [37]. Several researchers have used these smart models in the petroleum industry and to predict the solubility of acid gases in ILs. Bagban et al. [36] have predicted CO₂ solubility in the presence of various ILs using MLP and Adaptive Neuro-Fuzzy Inference System (ANFIS). Amedi et al. [38] have used MLP, RBF, and ANFIS to predict H₂S solubility in the presence of various ILs. Otherwise, in 2017, Rostami et al. [39] have applied the GEP method to model CO₂ solubility in crude oil during carbon dioxide enhanced oil recovery.

The current work aims at developing highly robust and easy-to-use machine learning models that can be applied for forecasting the solubility of CO₂ in 13 different ionic liquids at different temperature and pressure conditions. Two rigorous connectionist techniques, namely multilayer perceptron (MLP), and gene expression programming (GEP) are applied on a set of experimental data that was gathered from different literature sources [40–44]. The MLP method was optimized using either Levenberg-Marquardt (LMA) or Bayesian Regularization (BR) techniques. The results obtained using three methods (MLP-LMA, MLP-BR, GEP) are then compared with results calculated using thermodynamic models based on Peng-Robinson (PR) and Soave-Redlich-Kwong (SRK) equation of states. Statistical indicators, including the determination coefficient (R²) and Root Mean Square Error (RMSE), are used to evaluate the accuracy of the methods, in addition to graphical assessments using cross plots and bar plots. In the end, outliers detection is performed to test and analyze the validity of the best-developed model and quantify the doubtful experimental points from the database. It is worth noting that the two backpropagation-based learning algorithms (LMA and BR) that were employed in the training process of MLP, alongside the explicit correlations established to predict the CO₂ solubility in ILs, make the current work different from previously published works in the literature.

This paper is constructed as follows; Section 2 depicts the data used in the study and the input and output parameters in the models. Section 3 describes in detail the rigorous connectionist models and optimization techniques. An overview of the PR and SRK equations of state is presented in Section 4. In Section 5, the results are presented and discussed, and in Section 6, the conclusions of the study are summarized.

## 2. Data Collection and Preparation

In order to develop reliable models, it was crucial to use a large number of experimental data points representing a variety of conditions. In the current work, 744 experimental data points for carbon dioxide solubility in 13 different ILs (36 to 80 points for each IL) were collected from the literature [40–44]. Table 1 presents the references from where the data were gathered and the observed range of experimental measurements. The dissolved mole fraction (termed x) of CO₂ in each IL was the parameter we wished to match and predict based on the corresponding input of temperature *T*, pressure *P*, and the given IL. To represent the IL, we specified the thermodynamic properties of the IL, including molecular weight (*M*<sub>w</sub>), critical temperature (*T*<sub>c</sub>), critical pressure (*P*<sub>c</sub>), and acentric factor
$(w)$, as summarized in Table 2. These 6 input parameters were hence used for predicting the dissolved mole fraction $x$ of CO$_2$. All the 744 experimental data points were randomly divided into a training data set (80% of the database) and a test data set (20% of the database).

<table><caption>**Table 1.** Ranges of observed temperature, pressure, and CO$_2$ solubility measurements of the ionic liquids used in this study.</caption><thead><tr><th>No.</th><th>Ionic Liquid</th><th>Temperature Range (K)</th><th>Pressure Range (MPa)</th><th>CO<sub>2</sub> Solubility Range (Mole Fraction)</th><th>No. of Data Points</th><th>References</th></tr></thead><tbody><tr><td>1</td><td>[C3mpy][Tf2N]</td><td>303.15–373.15</td><td>0.52–47.1</td><td>0.186–0.787</td><td>56</td><td>[40]</td></tr><tr><td>2</td><td>[bmmim][tf2N]</td><td>298.15–343.15</td><td>0.01–1.9</td><td>0.002–0.382</td><td>36</td><td>[41]</td></tr><tr><td>3</td><td>[P(5)mpyrr][Tf2N]</td><td>298.15</td><td>0.01–1.9</td><td>0.002–0.406</td><td>36</td><td>[41]</td></tr><tr><td>4</td><td>[HMIM][Tf2N]</td><td>303.15–373.15</td><td>0.42–45.28</td><td>0.165–0.824</td><td>64</td><td>[42]</td></tr><tr><td>5</td><td>[HMIM][TfO]</td><td>303.15–373.15</td><td>1.42–100.12</td><td>0.267–0.816</td><td>64</td><td>[42]</td></tr><tr><td>6</td><td>[HMIM][BF4]</td><td>303.15–373.15</td><td>1.2–41.69</td><td>0.212–0.622</td><td>48</td><td>[42]</td></tr><tr><td>7</td><td>[HMIM][MeSO4]</td><td>303.15–373.15</td><td>0.87–50.14</td><td>0.158–0.602</td><td>48</td><td>[42]</td></tr><tr><td>8</td><td>[HMIM][PF6]</td><td>303.15–373.15</td><td>0.3–55.63</td><td>0.216–0.691</td><td>48</td><td>[42]</td></tr><tr><td>9</td><td>[C2mim][SCN]</td><td>303.15–373.15</td><td>1.3–95.34</td><td>0.169–0.474</td><td>72</td><td>[43]</td></tr><tr><td>10</td><td>[C2mim][N(CN)2]</td><td>303.15–373.15</td><td>0.88–96.2</td><td>0.171–0.585</td><td>80</td><td>[43]</td></tr><tr><td>11</td><td>[C2mim][C(CN)3]</td><td>303.15–373.15</td><td>0.59–88.29</td><td>0.17–0.503</td><td>80</td><td>[43]</td></tr><tr><td>12</td><td>[BMP][Tf2N]</td><td>303.15–373.15</td><td>0.68–62.77</td><td>0.2276–0.8029</td><td>72</td><td>[44]</td></tr><tr><td>13</td><td>[BMP][MeSO4]</td><td>303.15–373.15</td><td>3.07–97.3</td><td>0.2871–0.6049</td><td>40</td><td>[44]</td></tr></tbody></table>

<table><caption>**Table 2.** Thermodynamic properties of ionic liquids (ILs) used in this study.</caption><thead><tr><th>No.</th><th>Ionic Liquid</th><th>M<sub>w</sub> (kg/kmole)</th><th>T<sub>c</sub> (K)</th><th>P<sub>c</sub> (MPa)</th><th>Acentric Factor w (-)</th><th>References</th></tr></thead><tbody><tr><td>1</td><td>[C3mpy][Tf2N]</td><td>408.38</td><td>1196.86</td><td>2.6749</td><td>0.2753</td><td>[40]</td></tr><tr><td>2</td><td>[bmmim][tf2N]</td><td>433.4</td><td>1255.8</td><td>2.031</td><td>0.3193</td><td>[41]</td></tr><tr><td>3</td><td>[P(5)mpyrr][Tf2N]</td><td>436.4</td><td>1221.9</td><td>1.828</td><td>0.2603</td><td>[41]</td></tr><tr><td>4</td><td>[HMIM][Tf2N]</td><td>447.92</td><td>1292.78</td><td>2.3888</td><td>0.3893</td><td>[42]</td></tr><tr><td>5</td><td>[HMIM][TfO]</td><td>316.34</td><td>1055.6</td><td>2.4954</td><td>0.489</td><td>[42]</td></tr><tr><td>6</td><td>[HMIM][BF4]</td><td>278.37</td><td>1110.84</td><td>2.9611</td><td>0.4899</td><td>[42]</td></tr><tr><td>7</td><td>[HMIM][MeSO4]</td><td>254.08</td><td>716.61</td><td>1.7941</td><td>0.6589</td><td>[42]</td></tr><tr><td>8</td><td>[HMIM][PF6]</td><td>312.24</td><td>759.16</td><td>1.5499</td><td>0.9385</td><td>[42]</td></tr><tr><td>9</td><td>[C2mim][SCN]</td><td>169.25</td><td>1013.63</td><td>2.226</td><td>0.3931</td><td>[43]</td></tr><tr><td>10</td><td>[C2mim][N(CN)2]</td><td>177.21</td><td>998.96</td><td>2.9108</td><td>0.7661</td><td>[43]</td></tr><tr><td>11</td><td>[C2mim][C(CN)3]</td><td>201.23</td><td>1149.26</td><td>2.4591</td><td>0.8509</td><td>[43]</td></tr><tr><td>12</td><td>[BMP][Tf2N]</td><td>422.41</td><td>1209.16</td><td>2.48</td><td>0.32</td><td>[44]</td></tr><tr><td>13</td><td>[BMP][MeSO4]</td><td>253.36</td><td>1023.74</td><td>3.09</td><td>0.42</td><td>[44]</td></tr></tbody></table>

# 3. Modeling Techniques

## 3.1. Multilayer Perceptron (MLP)

Artificial neural network (ANN) is one of the well-known machine learning tools, which exhibits reliable ability for recognizing and identifying relationships between input and output parameters in complex systems [45]. The conception and the mathematical formulation of an ANN model were inspired by the human brain and its strategy for processing information. Multilayer perceptron (MLP) is one of the most applied types of ANN for modeling purposes.

An MLP model includes three kinds of layers, an input layer: From which the inputs enter into the model. The *i*th input parameter of the model is noted $z_i$.

One or more hidden layers: Their role is to transform the inputs into a higher feature in order to capture the non-linearity of the system. This mechanism is done by means of the activation functions. Tansig and logsig are among the frequently applied transfer functions and they are defined as follows:

$$ \text{Logsig: } g(z) = \frac{1}{e^z + 1} \qquad (1) $$

$$ \text{Tansig: } g(z) = \frac{e^z - e^{-z}}{e^z + e^{-z}} \qquad (2) $$

An output layer: The outputs of the model are returned from this layer. Pureline is generally considered the transfer function for this layer.

$$ \text{Pureline: } f(z) = z \qquad (3) $$
If the total number of hidden layers is considered to be *M*, the input layer is defined as layer 0. Suppose that the number of neurons in the hidden layer *m* is $N_m$, *m* = 1, 2, ..., *M*. Considered $w_{ij}^m$ the weight of the linkage between the *j*th neuron of the *m* – 1 th hidden layer and *i*th neuron of the *m*th hidden layer, and $b_i^m$ the bias term of *i*th neuron of the *m*th hidden layer. Hence, the output ($\bar{y}_i^m$) of *i*th neuron of the *m*th hidden layer can be determined using the following equation:

$$ \bar{y}_i^m = f\left(\sum_{j=1}^{N_{m-1}} w_{ij}^m \bar{y}_i^{m-1} + b_i^m\right) \quad (4) $$

$$ \bar{y}_i^0 = z_i, \quad i = 1, \dots, N_z, \quad N_z = N_0 \quad (5) $$

Each of the above-mentioned layers includes neurons. The number of neurons in the input and output layers corresponds to the number of input and output parameters of the system, respectively. The number of hidden layers and their neurons depends on the complexity of the system. Generally, one hidden layer is sufficient to model systems having moderate complexity, while more than one hidden layer is recommended for highly complex systems. The trial and error method is usually applied for the proper determination of the number of hidden layers and their numbers of neurons.

The training phase of an MLP model consists of determining appropriate values for the weights and bias terms that result in the minimization of the difference between the experimental data (measured dissolved CO₂ mole fraction) and the predictions of the model (calculated CO₂ mole fraction as a function of the input parameters). Backpropagation learning techniques were applied for the training of the MLP model, namely the Bayesian Regularization (BR) and Levenberg-Marquardt (LM) algorithms. More details about these two algorithms can be found in prior published works [46–48].

## 3.2. Gene Expression Programming (GEP)

In spite of computer-based artificial intelligent models, based genetic calculations, such as genetic programming (GP) can construct formulas, which have high precision. It has been demonstrated that GP can be more accurate than traditional empirical models. Furthermore, models developed by genetic calculations and simulators and software systems are easy to use [49].

Using correction and development of the GP method [50] led to a more relevant version, which was GEP [51]. Indeed, the latter removed the found weaknesses in GP, i.e., the big number of possible programs that can be built by the algorithm and the high sensitivity to small changes in fitness of the created solution. In general, the GEP [51] had 2 significant parts: The chromosome and the symbolic expression trees (ETs). The chromosome encoded the potential solutions, and after that, these solutions were converted to the particular applicant solution named ET [52–54]. Reliable ETs continuously were reproduced by unlimited genetic conversions because of the conversion method of the structural organization of the chromosome to the ET [52]. In fact, the unique genes’ placements in the GEP method produced impressive development of the solutions leading to encryption of any downside [55]. Preceding researches recognized that the speed of convergence in the GEP scheme was much more than the GP scheme [53,54].

The chromosome included unchanging factors and variables as terminals and predetermined functions with equal length in one or more genes [52]. By the user, the constants were produced by the GEP in a specific range, however, the functions and variables were the input information. Each gene consisted of a tail, which comprised terminals and also a head that contained constants and variables in the role of terminals and functions [52]. The given value of the gene’s head is taken into account as an input of the GEP method, and the value of gene’s tail (t) is given by the following equation:

$$ t = h(r - 1) + 1 \quad (6) $$

with *t* denotes the length of the tail, *h* is the length of the head, and *r* is the number of arguments of the function with the most arguments.
The essential stages of the GEP method are illustrated in Figure 1. These steps can be elucidated briefly as follows: Firstly, a generation of chromosomes of an initial population was created randomly. In a second step, the chromosomes were expressed with ET, and their fitness was computed. Then, the individuals were selected to perform reproduction based on their fitness values. This reproduction created offspring with fresh signs. The same procedure was followed by individuals of the resulted generation, i.e., genomes expression, selection, and reproduction with adjustment. The stages expressed above were repeated until reaching the stopping criteria. More explanation of the GEP technique can be found in the literature [51,55].

<figure><img src="image_4.png" /><figcaption>Figure 1. Flowchart of gene expression programming (GEP) algorithm procedure.</figcaption></figure>

In the current work, the mathematical expression used to predict the CO₂ solubility in ILs was generated in the following form: xCO₂ = f(T, P, Tc, IL, Pc, IL and wIL).

### 3.3. Techniques Evaluation

Statistical indicators alongside the graphical techniques were used to assess the performance and robustness of the proposed models. The statistical indicators include root mean square error (RMSE) and the coefficient of determination (R²), which is defined mathematically as follows:

$$RMSE = \sqrt{\frac{1}{n} \sum_{i=1}^{n} (x_{i \text{exp}} - x_{i \text{pred}})^2} \quad (7)$$

$$R^2 = 1 - \frac{\sum_{i=1}^{n} (x_{i \text{exp}} - x_{i \text{pred}})^2}{\sum_{i=1}^{n} (x_{i \text{exp}} - \bar{x}_{i \text{pred}})^2} \quad (8)$$
The credibility of the developed models was tested using the leverage statistical approach. In this approach, the standardized residuals, which refer to the difference between the predicted results and the experimental points, were presented versus the statistical Hat matrix leverage values; this graph is called the Williams plot. The Hat can be calculated using the Hat matrix of the following form [56,57]:

$$H = X(X^T X)^{-1}X^T \qquad (9)$$

where X is an (k × d) matrix, k and d represent the dimension and the data points number, respectively, and X^t denotes the X transpose matrix. The limit leverage value (H*) depicted on the Williams plot with a purple line is a constant calculated as 3(d + 1)/k. The selection of data points is in the range of ± 3 of the standard deviation from the mean, where the cut-off value of 3 ensures covering 99% of normally distributed data. The validity of the developed model and the credibility of the predictions were conditioned by the values of data points, which must be situated in the range of 0 ≤ H ≤ H* and -3 ≤ R ≤ 3.

## 4. Equations of State for Modeling the Solubility of Acid Gases

Many equations of state have been used by researchers to model the solubility of acid gases in ILs. Two of the most widely used EoS are the Soave–Redlich–Kwong (SRK) and the Peng–Robinson (PR) EoS, which are defined with Equations (10) and (15), respectively [35]. The results of these two EoS are used to contrast those obtained by using the proposed models.

Noting that the calculation of solubility using EoS is related to the calculation of the mole fraction, the equations below express how to determine the mole fraction.

Soave–Redlich–Kwong (SRK) EoS:

$$P = \frac{RT}{v-b} - \frac{a(T)}{v(v+b)} \qquad (10)$$

where T, P, v, and R indicate the temperature, pressure, molar volume, and gas constant respectively, and a and b represent the EoS variables.

For computing the variables, a and b in the case of mixtures, the classic van der Waals one-fluid rules of mixing is used [58]:

$$a = \sum_{i=1}^{N} \sum_{j=1}^{N} x_i x_j \sqrt{a_i a_j (1 - K_{ij}(T))}; \quad b = \sum_{i=1}^{N} x_i b_i \qquad (11)$$

with x_i and x_j denote the mole fractions of components i and j, and N indicates the number of the mixture's components. The parameter K_{ij}(T) is a parameter of binary interaction, that enlarges the molecular interactions between molecules i and j. For pure materials, the parameters can be computed as follows:

$$a_i(T) = 0.42747 \frac{R^2 T_{c,i}^2}{P_{c,i}} \alpha_i(T_R) \qquad (12)$$

$$b_i = \frac{0.08664RT_{c,i}}{P_{c,i}} \qquad (13)$$

$$\alpha_i(T_R) = [1 + (0.480 + 1.574w_i - 0.176w_i^2)(1 - \sqrt{T_{Ri}})]^2 \qquad (14)$$

with T_{Ri} = T/T_{ci}, denotes the reduced temperature. We note that the above description is valid for a system of arbitrarily many components. In our system, it was assumed that we had two components where CO₂ was one and the IL the other.
Peng–Robinson (PR) EoS:

$$P = \frac{RT}{v-b} - \frac{a(T)}{v(v+b)+b(v-b)} \qquad (15)$$

with:

$$a_i(T) = 0.457235 \frac{R^2 T_{c,i}^2}{P_{c,i}} * [1 + (0.37464 + 1.54226w_i - 0.26992w_i^2)(1 - \sqrt{T_{Ri}})]^2 \qquad (16)$$

$$b_i = \frac{0.077796RT_{c,i}}{P_{c,i}} \qquad (17)$$

For mixtures, the parameters a and b are determined using the same formulas as in SRK EoS. The parameter of binary interaction ($K_{ij}$) is the only adjustable parameter for both EoS, SRK, and PR. This parameter is determined using a genetic algorithm (GA) with the next objective function [58]:

$$f = \sum_{i=1}^{N} \left( \frac{x_{\text{solute}}^{\text{exp}} - x_{\text{solute}}^{\text{cal}}}{x_{\text{solute}}^{\text{exp}}} \right)^2 \qquad (18)$$

with $x_{\text{solute}}^{\text{exp}}$ and $x_{\text{solute}}^{\text{cal}}$ indicate the experimental and EoS-calculated mole fractions of the solute respectively.

## 5. Computational Procedure

For the training phase, the mean square error (MSE) was used as the assessment criterion, which is defined mathematically as follows:

$$MSE = \frac{\sum_{i=1}^{n} (x_{i \text{exp}} - x_{i \text{pre}})^2}{n} \qquad (19)$$

where x stands for mole fraction of CO₂, exp and pre indicate the experimental and the predicted values, respectively, and n represents the number of samples. The model tuning parameters giving the lowest MSE on the training set were considered the choice for the trained model.

For the modeling task using MLP, the data points were normalized between -1 and 1. To select appropriate topologies for the MLP approach, trial and error were used. The obtained models were designated MLP-LMA, and MLP-BR, respectively, and both included 3 hidden layers with 11, 11, and 9 neurons, respectively. The suitable activation functions in all the hidden layers and for the output layer were Tansig and Pureline, respectively.

## 6. Results and Discussion

To obtain an accurate and trustworthy expression, the calculation procedure reported in Section 3.2 has been followed. Thus, taking into account the series of five independent parameters ($T$, $P$, $T_c$, $P_c$, $w$), the GEP-based model makes the mathematical based correlation for the dependent parameter $xCO_2$ as follows:

$$x = 0.0001011 \times A_1 + A_2 \tanh(P) - A_3 \sqrt{P} + A_4 \times \ln(P) + 10^{-6} \times M_w^2 \times A_5 + A_6 \times P_c^3 + A_7 \times T_c + A_8 \qquad (20)$$

where $P$ (MPa) and $P_c$ (MPa) are pressure and critical pressure, respectively, $M_w$ (Kg/Kmole) is the molecular weight, $T_c$ (K) indicates the critical temperature and $A_1, A_2, A_3, \dots$ and $A_8$ are expressed as below:

$$A_1 = P_c \times w + P \times (P_c + w + P) \qquad (21)$$

$$A_2 = 3.362 - 0.1832 \times \sqrt{T} \qquad (22)$$
$$A_3 = 0.4802 + \frac{0.08072 \times T}{M_w} + 0.03085 \times \sqrt{T} + 5.14 \times 10^{-6} \times M_w \times \sqrt{P} \times \ln(T) \quad (23)$$

$$A_4 = \frac{0.3091}{\tanh(P_c)} - \frac{0.4344}{P_c^3} - 0.002638 \times \exp(P_c) \quad (24)$$

$$A_5 = 4.047 - 5.481 \times 10^{-3} \times M_w \quad (25)$$

$$A_6 = 2.33 \times 10^{-5} \times T_c \times \tanh(M_w) - 0.0004489 \times P_c^3 \quad (26)$$

$$A_7 = 4.047 \times 10^{-6} - 0.000305 \times w \quad (27)$$

$$A_8 = -5.213 \times 10^{-7} \times P^3 - 0.146 \times P_c + 0.5632 \times w + 0.04124 \times \exp\left(\frac{T}{M_w}\right) - 0.2381 \times \ln(\tanh(P)) - 0.339 \quad (27)$$

## 6.1. Performance Evaluation

Table 3 summarizes the statistical parameters (RMSE and R²) of the three proposed models in predicting CO₂ solubility for training, test, and overall data. The graphical representation of the overall RMSE and R² through bar plots, as shown in Figure 2. This graphical comparison allows an alternative visualization for the performance of the different models. A concise view of the figure shows that the values of RMSE and R² were closest to 0 and 1, respectively, pointing out that the models were more credible and qualified for the forecast.

<table><caption>**Table 3.** Statistical parameters of the proposed models in predicting CO₂ solubility.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Training Data</th><th colspan="2">Test Data</th><th colspan="2">Overall</th></tr><tr><th>RMSE</th><th>R<sup>2</sup></th><th>RMSE</th><th>R<sup>2</sup></th><th>RMSE</th><th>R<sup>2</sup></th></tr></thead><tbody><tr><td>MLP-LMA</td><td>0.0107</td><td>0.9971</td><td>0.0150</td><td>0.9941</td><td>0.0116</td><td>0.9965</td></tr><tr><td>MLP-BR</td><td>0.0117</td><td>0.9966</td><td>0.0138</td><td>0.9951</td><td>0.0121</td><td>0.9963</td></tr><tr><td>GEP</td><td>0.0192</td><td>0.9907</td><td>0.0238</td><td>0.9854</td><td>0.0201</td><td>0.9896</td></tr></tbody></table>

Form a deep comparative standing point of view, Table 3 and Figure 2 show that the reliability of the developed models took the following order: MLP-LMA > MLP-BR > GEP. The MLP model optimized using LMA had the highest R² coefficients and the lowest RMSE values for both the training and overall data, followed by the MLP-BR model. On the other hand, very satisfactory performances and high predictions potential of the GEP model were asserted by the values reported in Table 3 for training, test, and overall data.

Cross plots in Figure 3 were used to compare the models' predictions against the experimental data further. In these plots, the predicted values using the three methods were represented versus the previously established experimental data for the solubility of CO₂ in ILs. The criterion of precision for the different models in typical cross plots can be indicated by the compressed gathering of data points cloud nearby the 45° line, otherwise, all points falling on the unit slope would correspond to a perfect model. The accumulation of data cloud nearby the 45° line was attributed to the elevated value of R² close to the unit and to the low value of RMSE near zero.

It can be observed in Figure 3, that the MLP-LMA, MLP-BR, and GEP models exhibited distributions of predictions near the unit slope line for both training and test datasets, confirming their high accuracy in estimating the solubility of CO₂ in the considered ionic liquids.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Performance comparison on the overall dataset for the proposed models: (a) Root mean square error (RMSE), and (b) R².</figcaption></figure>
## 6.2. Comparison with PR and SRK EoS Results

The results obtained by the MLP-LMA model were compared with values calculated using thermodynamic models based on PR and SRK EoS.

Tables 4 and 5 provide the CO₂ solubility and thermodynamic properties respectively. Furthermore, Table 6 summarizes a comparison, where the RMSEs and R² of the best-proposed model, i.e., MLP-LMA, and PR and SRK EoS were mentioned. To perform this comparison, two samples of ILs that had not been included in the training phase of the intelligent models were considered. These samples included [bmim][Tf2N] and [bmim][PF6] ILs. Further description is given in Tables 4 and 5. All three models (MLP-LMA, and the PR and SRK EoS) gave accurate descriptions with R² above 0.98. For the two ILs, the MLP-LMA model had a higher coefficient of determination and lower RMSE compared to the PR and SRK EoS, in the following order of reliability: MLP-LMA > SRK > PR in the case of [bmim][PF6], and MLP-LMA > PR > SRK in the case of [bmim][Tf2N].

<table><caption>**Table 4.** Ranges of temperature, pressure, and CO₂ solubility in [bmim][Tf2N] and [bmim][PF6] ILs.</caption><thead><tr><th>No.</th><th>Ionic Liquid</th><th>T (K)</th><th>P (MPa)</th><th>x CO<sub>2</sub></th><th>References</th></tr></thead><tbody><tr><td>1</td><td>[bmim][PF6]</td><td>298.15</td><td>1.33–5.37</td><td>0.2054–0.5604</td><td>[59]</td></tr><tr><td>2</td><td>[bmim][Tf2N]</td><td>298.15</td><td>0.42–5.91</td><td>0.0973–0.7117</td><td>[60]</td></tr></tbody></table>

<table><caption>**Table 5.** Thermodynamic properties of [bmim][TF2N] and [bmim][PF6] ILs.</caption><thead><tr><th>No.</th><th>Ionic Liquid</th><th>M<sub>w</sub><br>(g/mole)</th><th>T<sub>c</sub> (K)</th><th>P<sub>c</sub> (MPa)</th><th>ω</th><th>References</th></tr></thead><tbody><tr><td>1</td><td>[bmim][PF6]</td><td>284.18</td><td>708.9</td><td>1.73</td><td>0.7553</td><td>[59]</td></tr><tr><td>2</td><td>[bmim][Tf2N]</td><td>419.36</td><td>1265.0</td><td>2.76</td><td>0.2656</td><td>[60]</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 3.** *Cont.*</figcaption></figure>
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 3.** Cross plots of the proposed predictive models.</figcaption></figure>

Figures 4 and 5 exhibited the comparison performed between PR and SRK EoS and the MLP-LMA model, and at the same time, the evolution of CO₂ solubility with pressure in the two ILs at a constant temperature. As it is clear in the figures, the solubility of CO₂ increased with the increase in pressure for the two investigated ILs. Henry's law can explain this phenomenon, which supposes that the gas
solubility throughout a liquid is directly a function of the gas pressure above the solution's surface [28]. The increase in pressure forces gas molecules into the solution, thus the number of gas molecules dissolved raises. It can also be seen from the figures that the solubility values predicted by the MLP-LMA model are in good agreement with the experimental data, whereas those established using the EoS differs slightly.

<table><caption>**Table 6.** Comparison of R² and RMSE of EoS and the best proposed model.</caption><thead><tr><th rowspan="2">ILs</th><th colspan="2">PR EoS</th><th colspan="2">SRK EoS</th><th colspan="2">MLP-LMA Model</th></tr><tr><th>R²</th><th>RMSE</th><th>R²</th><th>RMSE</th><th>R²</th><th>RMSE</th></tr></thead><tbody><tr><td>[bmim][PF6]</td><td>0.9842</td><td>0.0158</td><td>0.9931</td><td>0.0104</td><td>0.9994</td><td>0.0030</td></tr><tr><td>[bmim][Tf2N]</td><td>0.9960</td><td>0.0125</td><td>0.9951</td><td>0.0138</td><td>0.9980</td><td>0.0088</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Comparison of Peng–Robinson (PR) and Soave–Redlich–Kwong equations of state (SRK EoS), and MLP-LMA model for prediction of CO₂ solubility in [bmim][PF6] at T = 298.15 K.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 5.** Comparison of PR and SRK EoS, and MLP-LMA model for prediction of CO₂ solubility in [bmim][Tf2N] at T = 298.15 K.</figcaption></figure>
### 6.3. Outlier Detection

A Williams plot of the MLP-LMA model is depicted in Figure 6. It can be seen from the figure that the majority of data points are situated in the range of -3 ≤ R ≤ 3 and 0 ≤ H ≤ 0.03, which validated and verified the implemented MLP-LMA model statistically, where 13 points are considered suspect. In summary, the leverage approach states that the developed model is reliable and can be used with high accuracy to predict CO₂ solubility in ILs.

<figure><img src="image_11.png" /><figcaption>Figure 6. Williams plot of the MLP-LMA model.</figcaption></figure>

## 7. Conclusions

In the current paper, based on different experimental data gathered from the literature, MLP-LMA, GEP, PR, and SRK EoS methods were successfully used to predict carbon dioxide solubility in 13 various ionic liquids under the (298.15–373.15 K) temperature and (0.01–100.12 MPa) pressure range. Multilayer perceptron optimized with two back-propagation algorithms, viz. LMA and BR, and a reliable white-box technique, namely GEP, were the applied data-driven methods. The various statistical indicators obtained in this study showed that MLP optimized with LMA was the most advantageous paradigm, which can be used with high precision to predict the CO₂ solubility in ILs for natural gas sweetening and purification process, which may be of great interest for experts and gas engineers. Furthermore, the established correlation using the GEP-based model was reliable and gave good results. Finally, the comparison made between the MLP-LMA model and PR and SRK EoS demonstrated the high performance of the proposed model against the cubic EoS.

**Author Contributions:** Data curation, H.O., N.N., and A.M.; formal analysis, H.O., P.Ø.A., and A.M.; funding acquisition, A.M.; investigation, A.H.H., M.E.A.B.S., P.Ø.A., and S.S.; methodology, A.H.H. and M.E.A.B.S.; project administration, M.N.A., M.A.G., and N.N.; resources, M.N.A., M.A.G., and N.N.; software, M.N.A., M.A.G., and N.N.; supervision, S.S.; writing—original draft, H.O.; writing—review and editing, P.Ø.A. and A.M. All authors have read and agreed to the published version of the manuscript.

**Funding:** We acknowledge the financial support of this work by the Hungarian State and the European Union under the EFOP-3.6.1-16-2016-00010 project.

**Acknowledgments:** We acknowledge the support of the German Research Foundation (DFG) and the Bauhaus-Universität Weimar within the Open-Access Publishing Programme. Furthermore, we acknowledge the Research Council of Norway and the industry partners, ConocoPhillips Skandinavia AS, Aker BP ASA, Vår Energi AS, Equinor ASA, Neptune Energy Norge AS, Lundin Norway AS, Halliburton AS, Schlumberger Norge AS, and Wintershall DEA, of The National IOR Centre of Norway for support.

**Conflicts of Interest:** The authors declare no conflict of interest.
## Abbreviations and Symbols

<table><tr><td>ANN</td><td>Artificial neural network</td></tr><tr><td>BR</td><td>Bayesian Regularization</td></tr><tr><td>CO<sub>2</sub></td><td>Carbon dioxide</td></tr><tr><td>d</td><td>Data points number</td></tr><tr><td>DEA</td><td>Diethanolamine</td></tr><tr><td>EoS</td><td>Equations of state</td></tr><tr><td>ETs</td><td>Expression trees</td></tr><tr><td>exp</td><td>Experimental</td></tr><tr><td>FL</td><td>Fuzzy logic</td></tr><tr><td>GEP</td><td>Gene expression programming</td></tr><tr><td>GP</td><td>Genetic programming</td></tr><tr><td>H</td><td>Hat matrix</td></tr><tr><td>H<sub>2</sub>S</td><td>Hydrogen sulfide</td></tr><tr><td>H*</td><td>Limit Leverage value</td></tr><tr><td>ILs</td><td>Ionic liquids</td></tr><tr><td>LM</td><td>Levenberg-Marquardt</td></tr><tr><td>MEA</td><td>Monoethanolamine</td></tr><tr><td>MLP</td><td>Multilayer perceptron</td></tr><tr><td>MSE</td><td>Mean square error</td></tr><tr><td>M<sub>w</sub></td><td>Molecular weight (kg/kmole)</td></tr><tr><td>k</td><td>Dimension</td></tr><tr><td>n</td><td>Number of samples</td></tr><tr><td>P<sub>c</sub></td><td>Critical pressure (MPa)</td></tr><tr><td>pred</td><td>Predicted</td></tr><tr><td>PR</td><td>Peng-Robinson</td></tr><tr><td>RMSE</td><td>Root mean square error</td></tr><tr><td>R<sup>2</sup></td><td>Coefficient of determination</td></tr><tr><td>SRK</td><td>Soave-Redlich-Kwong</td></tr><tr><td>SVM</td><td>Support vector machine</td></tr><tr><td>t</td><td>Gene's tail</td></tr><tr><td>T<sub>c</sub></td><td>Critical temperature (K)</td></tr><tr><td>TEA</td><td>Triethanolamine</td></tr><tr><td>w</td><td>Acentric factor</td></tr><tr><td>x</td><td>Mole fraction of CO<sub>2</sub></td></tr><tr><td>X</td><td>(k × d) matrix</td></tr><tr><td>X<sup>t</sup></td><td>X transpose matrix</td></tr></table>

## References

1. Taib, M.M.; Murugesan, T. Solubilities of CO₂ in aqueous solutions of ionic liquids (ILs) and monoethanolamine (MEA) at pressures from 100 to 1600 kPa. *Chem. Eng. J.* **2012**, *181*, 56–62. [CrossRef]

2. Kidnay, A.J.; Parrish, W.R.; McCartney, D.G. *Fundamentals of Natural Gas Processing*; CRC Press: Boca Raton, FL, USA, 2011.

3. Astaria, G.; Savage, D.W.; Bisio, A. *Gas Treating with Chemical Solvents*; John Wiley: New York, NY, USA, 1983.

4. Romeo, L.M.; Bolea, I.; Escosa, J.M. Integration of power plant and amine scrubbing to reduce CO₂ capture costs. *Appl. Therm. Eng.* **2008**, *28*, 1039–1046. [CrossRef]

5. Mulukutla, T.; Obuskovic, G.; Sirkar, K.K. Novel scrubbing system for post-combustion CO₂ capture and recovery: Experimental studies. *J. Membr. Sci.* **2014**, *471*, 16–26. [CrossRef]

6. Jalili, A.H.; Rahmati-Rostami, M.; Ghotbi, C.; Hosseini-Jenab, M.; Ahmadi, A.N. Solubility of H₂S in ionic liquids [bmim][PF₆], [bmim][BF₄], and [bmim][Tf₂N]. *J. Chem. Eng. Data* **2009**, *54*, 1844–1849. [CrossRef]

7. Shiflett, M.B.; Yokozeki, A. Separation of CO₂ and H₂S using room-temperature ionic liquid [bmim][PF₆]. *Fluid Phase Equilib.* **2010**, *294*, 105–113. [CrossRef]
8. Sakhaeinia, H.; Taghikhani, V.; Jalili, A.H.; Mehdizadeh, A.; Safekordi, A.A. Solubility of H₂S in 1-(2-hydroxyethyl)-3-methylimidazolium ionic liquids with different anions. *Fluid Phase Equilib.* **2010**, *298*, 303–309. [CrossRef]

9. Sánchez, L.G.; Meindersma, G.; De Haan, A. Solvent properties of functionalized ionic liquids for CO₂ absorption. *Chem. Eng. Res. Des.* **2007**, *85*, 31–39. [CrossRef]

10. Safavi, M.; Ghotbi, C.; Taghikhani, V.; Jalili, A.H.; Mehdizadeh, A. Study of the solubility of CO₂, H₂S and their mixture in the ionic liquid 1-octyl-3-methylimidazolium hexafluorophosphate: Experimental and modelling. *J. Chem. Thermodyn.* **2013**, *65*, 220–232. [CrossRef]

11. Han, D.; Tian, M.; Park, D.-W.; Choi, D.K. Application of ionic liquids as mobile phase additives and surface-bonded stationary phase in liquid chromatography. *Korean J. Chem. Eng.* **2009**, *26*, 1353–1358. [CrossRef]

12. Martínez-Paloua, R.; Luqueb, R. Applications of Ionic liquids in contaminants removal from Refinery Feedstocks: An industrial perspective. *Energy Environ. Sci.* **2014**, *7*, 2414–2447. [CrossRef]

13. Bemani, A.; Baghban, A.; Shamshirband, S.; Mosavi, A.; Csiba, P.; Varkonyi-Koczy, A.R. Applying ANN, ANFIS, and LSSVM Models for Estimation of Acid Solvent Solubility in Supercritical CO₂. *arXiv* **2019**, arXiv:1912.05612. [CrossRef]

14. Shamshirband, S.; Hadipoor, M.; Baghban, A.; Mosavi, A.; Bukor, J.; Várkonyi-Kóczy, A.R. Developing an ANFIS-PSO Model to Predict Mercury Emissions in Combustion Flue Gases. *Mathematics* **2019**, *7*, 965.

15. Hasib-ur-Rahman, M.; Siaj, M.; Larachi, F. Ionic liquids for CO₂ capture—Development and progress. *Chem. Eng. Process.* **2010**, *49*, 313–322. [CrossRef]

16. Bara, J.E. Potential for hydrogen sulfide removal using ionic liquid solvents. In *Green Solvents II*; Springer: Berlin/Heidelberg, Germany, 2012; pp. 155–167.

17. Davis, H.; James, J. Task-specific ionic liquids. *Chem. Lett.* **2004**, *33*, 1072–1077. [CrossRef]

18. Akhmetshina, A.I.; Petukhov, A.N.; Gumerova, O.R.; Vorotyntsev, A.V.; Nyuchev, A.V.; Vorotyntsev, I.V. Solubility of H₂S and CO₂ in imidazolium-based ionic liquids with bis (2-ethylhexyl) sulfosuccinate anion. *J. Chem. Thermodyn.* **2019**, *130*, 173–182. [CrossRef]

19. Jalili, A.H.; Mehrabi, M.; Zoghi, A.T.; Shokouhi, M.; Taheri, S.A. Solubility of carbon dioxide and hydrogen sulfide in the ionic liquid 1-butyl-3-methylimidazolium trifluoromethanesulfonate. *Fluid Phase Equilib.* **2017**, *453*, 1–12. [CrossRef]

20. Yokozeki, A.; Shiflett, M.B.; Junk, C.P.; Grieco, L.M.; Foo, T. Physical and chemical absorptions of carbon dioxide in room-temperature ionic liquids. *J. Phys. Chem. B* **2008**, *112*, 16654–16663. [CrossRef]

21. Kurnia, K.; Harris, F.; Wilfred, C.; Mutalib, M.A.; Murugesan, T. Thermodynamic properties of CO₂ absorption in hydroxyl ammonium ionic liquids at pressures of (100–1600) kPa. *J. Chem. Thermodyn.* **2009**, *41*, 1069–1073. [CrossRef]

22. Carvalho, P.J.; Álvarez, V.H.; Schröder, B.; Gil, A.M.; Marrucho, I.M.; Aznar, M.; Santos, L.M.; Coutinho, J.A. Specific solvation interactions of CO₂ on acetate and trifluoroacetate imidazolium based ionic liquids at high pressures. *J. Phys. Chem. B* **2009**, *113*, 6803–6812. [CrossRef]

23. Shin, E.-K.; Lee, B.-C. High-pressure phase behavior of carbon dioxide with ionic liquids: 1-alkyl-3-methylimidazolium trifluoromethanesulfonate. *J. Chem. Eng. Data* **2008**, *53*, 2728–2734. [CrossRef]

24. Yuan, X.; Zhang, S.; Liu, J.; Lu, X. Solubilities of CO₂ in hydroxyl ammonium ionic liquids at elevated pressures. *Fluid Phase Equilib.* **2007**, *257*, 195–200. [CrossRef]

25. Vega, L.F.; Vilaseca, O.; Llovell, F.; Andreu, J.S. Modeling ionic liquids and the solubility of gases in them: Recent advances and perspectives. *Fluid Phase Equilib.* **2010**, *294*, 15–30. [CrossRef]

26. Rahmati-Rostami, M.; Ghotbi, C.; Hosseini-Jenab, M.; Ahmadi, A.N.; Jalili, A.H. Solubility of H₂S in ionic liquids [hmim][PF6], [hmim][BF4], and [hmim][Tf2N]. *J. Chem. Thermodyn.* **2009**, *41*, 1052–1055. [CrossRef]

27. Shokouhi, M.; Adibi, M.; Jalili, A.H.; Hosseini-Jenab, M.; Mehdizadeh, A. Solubility and diffusion of H₂S and CO₂ in the ionic liquid 1-(2-hydroxyethyl)-3-methylimidazolium tetrafluoroborate. *J. Chem. Eng. Data* **2009**, *55*, 1663–1668. [CrossRef]

28. Jalili, A.H.; Mehdizadeh, A.; Shokouhi, M.; Ahmadi, A.N.; Hosseini-Jenab, M.; Fateminassab, F. Solubility and diffusion of CO₂ and H₂S in the ionic liquid 1-ethyl-3-methylimidazolium ethylsulfate. *J. Chem. Thermodyn.* **2010**, *42*, 1298–1303. [CrossRef]

29. Jou, F.-Y.; Mather, A.E. Solubility of hydrogen sulfide in [bmim][PF6]. *Int. J. Thermophys.* **2007**, *28*, 490. [CrossRef]
30. Yokozeki, A.; Shiflett, M.B. Gas solubilities in ionic liquids using a generic van der Waals equation of state. *J. Supercrit. Fluids* **2010**, *55*, 846–851. [CrossRef]

31. Zhang, S.; Chen, Y.; Ren, R.X.-F.; Zhang, Y.; Zhang, J.; Zhang, X. Solubility of CO₂ in sulfonate ionic liquids at high pressure. *J. Chem. Eng. Data* **2005**, *50*, 230–233. [CrossRef]

32. Kim, Y.; Jang, J.; Lim, B.; Kang, J.W.; Lee, C. Solubility of mixed gases containing carbon dioxide in ionic liquids: Measurements and predictions. *Fluid Phase Equilib.* **2007**, *256*, 70–74. [CrossRef]

33. Al-fnaish, H.; Lue, L. Modelling the solubility of H₂S and CO₂ in ionic liquids using PC-SAFT equation of state. *Fluid Phase Equilib.* **2017**, *450*, 30–41. [CrossRef]

34. Cybenko, G. Approximation by superpositions of a sigmoidal function. *Math. Control Signals Syst.* **1989**, *2*, 303–314. [CrossRef]

35. Ahmadi, M.A.; Ebadi, M.; Hosseini, S.M. Prediction breakthrough time of water coning in the fractured reservoirs by implementing low parameter support vector machine approach. *Fuel* **2014**, *117*, 579–589. [CrossRef]

36. Baghban, A.; Ahmadi, M.A.; Shahraki, B.H. Prediction carbon dioxide solubility in presence of various ionic liquids using computational intelligence approaches. *J. Supercrit. Fluids* **2015**, *98*, 50–64. [CrossRef]

37. Vapnik, V.; Vapnik, V. *Statistical Learning Theory*; Wiley: New York, NY, USA, 1998.

38. Amedi, H.R.; Baghban, A.; Ahmadi, M.A. Evolving machine learning models to predict hydrogen sulfide solubility in the presence of various ionic liquids. *J. Mol. Liquids* **2016**, *216*, 411–422. [CrossRef]

39. Rostami, A.; Arabloo, M.; Kamari, A.; Mohammadi, A.H. Modeling of CO₂ solubility in crude oil during carbon dioxide enhanced oil recovery using gene expression programming. *Fuel* **2017**, *210*, 768–782. [CrossRef]

40. Kim, S.A.; Yim, J.-H.; Lim, J.S. High-pressure phase behavior of binary mixtures containing methylpyrrolidinium derivative ionic liquids and carbon dioxide. *Fluid Phase Equilib.* **2012**, *332*, 28–34. [CrossRef]

41. Tagiuri, A.; Sumon, K.Z.; Henni, A. Solubility of carbon dioxide in three [Tf₂N] ionic liquids. *Fluid Phase Equilib.* **2014**, *380*, 39–47. [CrossRef]

42. Yim, J.-H.; Lim, J.S. CO₂ solubility measurement in 1-hexyl-3-methylimidazolium ([HMIM]) cation based ionic liquids. *Fluid Phase Equilib.* **2013**, *352*, 67–74. [CrossRef]

43. Kim, J.E.; Kim, H.J.; Lim, J.S. Solubility of CO₂ in ionic liquids containing cyanide anions: [c₂mim][SCN], [c₂mim][N(CN)₂], [C₂mim][C(CN)₃]. *Fluid Phase Equilib.* **2014**, *367*, 151–158. [CrossRef]

44. Yim, J.-H.; Song, H.N.; Yoo, K.-P.; Lim, J.S. Measurement of CO₂ solubility in ionic liquids: [BMP][Tf₂N] and [BMP][MeSO₄] by measuring bubble-point pressure. *J. Chem. Eng. Data* **2011**, *56*, 1197–1203. [CrossRef]

45. Hemmati-Sarapardeh, A.; Ameli, F.; Varamesh, A.; Shamshirband, S.; Mohammadi, A.H.; Dabir, B. Toward generalized models for estimating molecular weights and acentric factors of pure chemical compounds. *Int. J. hydrog. Energy* **2018**, *43*, 2699–2717. [CrossRef]

46. Hemmati-Sarapardeh, A.; Varamesh, A.; Husein, M.M.; Karan, K. On the evaluation of the viscosity of nanofluid systems: Modeling and data assessment. *Renew. Sustain. Energy Rev.* **2018**, *81*, 313–329. [CrossRef]

47. Menad, N.A.; Hemmati-Sarapardeh, A.; Varamesh, A.; Shamshirband, S. Predicting solubility of CO₂ in brine by advanced machine learning systems: Application to carbon capture and sequestration. *J. CO₂ Util.* **2019**, *33*, 83–95. [CrossRef]

48. Menad, N.A.; Noureddine, Z.; Hemmati-Sarapardeh, A.; Shamshirband, S. Modeling temperature-based oil-water relative permeability by integrating advanced intelligent models with grey wolf optimization: Application to thermal enhanced oil recovery processes. *Fuel* **2019**, *242*, 649–663. [CrossRef]

49. Fathinasab, M.; Ayatollahi, S.; Hemmati-Sarapardeh, A. A rigorous approach to predict nitrogen-crude oil minimum miscibility pressure of pure and nitrogen mixtures. *Fluid Phase Equilib.* **2015**, *399*, 30–39. [CrossRef]

50. Koza, J.R.; Koza, J.R. *Genetic Programming: On the Programming of Computers by Means of Natural Selection*; MIT Press: Cambridge, MA, USA, 1992; Volume 1.

51. Ferreira, C. Gene expression programming: A new adaptive algorithm for solving problems. *Complex Syst.* **2019**, *3*, 183–199.

52. Teodorescu, L.; Sherwood, D. High energy physics event selection with gene expression programming. *Comput. Phys. Commun.* **2008**, *178*, 409–419. [CrossRef]

53. Rostami, A.; Ebadi, H. Toward gene expression programming for accurate prognostication of the critical oil flow rate through the choke: Correlation development. *Asia-Pac. J. Chem. Eng.* **2017**, *12*, 884–893. [CrossRef]
54. Rostami, A.; Shokrollahi, A. Accurate prediction of water dewpoint temperature in natural gas dehydrators using gene expression programming approach. *J. Mol. Liquids* **2017**, *243*, 196–204. [CrossRef]

55. Ferreira, C. *Gene Expression Programming: Mathematical Modeling by an Artificial Intelligence*; Springer: Berlin/Heidelberg, Germany, 2006; Volume 21.

56. Kamari, A.; Safiri, A.; Mohammadi, A.H. Compositional model for estimating asphaltene precipitation conditions in live reservoir oil systems. *J. Dispers. Sci. Technol.* **2015**, *36*, 301–309. [CrossRef]

57. Jaworska, J.; Nikolova-Jeliazkova, N.; Aldenberg, T. QSAR applicability domain estimation by projection of the training set in descriptor space: A review. *Altern. Lab. Anim.* **2005**, *33*, 445–459. [CrossRef]

58. Haghbakhsh, R.; Soleymani, H.; Raeissi, S. A simple correlation to predict high pressure solubility of carbon dioxide in 27 commonly used ionic liquids. *J. Supercrit. Fluids* **2013**, *77*, 158–166. [CrossRef]

59. Aki, S.N.; Mellein, B.R.; Saurer, E.M.; Brennecke, J.F. High-pressure phase behavior of carbon dioxide with imidazolium-based ionic liquids. *J. Phys. Chem. B* **2004**, *108*, 20355–20365. [CrossRef]

60. Anthony, J.L.; Anderson, J.L.; Maginn, E.J.; Brennecke, J.F. Anion effects on gas solubility in ionic liquids. *J. Phys. Chem. B* **2005**, *109*, 6366–6374. [CrossRef] [PubMed]

<figure><img src="image_12.png" /></figure>

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).