Next Article in Journal
The Impact of Filler Geometry on Polylactic Acid-Based Sustainable Polymer Composites
Next Article in Special Issue
Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes
Previous Article in Journal
BODIPY-Based Molecules, a Platform for Photonic and Solar Cells
Previous Article in Special Issue
The Proton Dissociation of Bio-Protic Ionic Liquids: [AAE]X Amino Acid Ionic Liquids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Viscosity of Ionic Liquids: Application of the Eyring’s Theory and a Committee Machine Intelligent System

by
Seyed Pezhman Mousavi
1,
Saeid Atashrouz
2,
Menad Nait Amar
3,
Abdolhossein Hemmati-Sarapardeh
1,4,*,
Ahmad Mohaddespour
5,* and
Amir Mosavi
6,7,8,*
1
Department of Petroleum Engineering, Shahid Bahonar University of Kerman, Kerman 76169-13439, Iran
2
Department of Chemical Engineering, Amirkabir University of Technology (Tehran Polytechnic), Hafez 424, Tehran 15875-4413, Iran
3
Département Etudes Thermodynamiques, Division Laboratoires, Sonatrach, Boumerdes 35000, Algeria
4
College of Construction Engineering, Jilin University, Changchun 130600, China
5
College of Engineering and Technology, American University of the Middle East, Dasman, Kuwait
6
Faculty of Civil Engineering, Technische Universität Dresden, 01069 Dresden, Germany
7
School of Economics and Business, Norwegian University of Life Sciences, 1430 Ås, Norway
8
John von Neumann Faculty of Informatics, Obuda University, 1034 Budapest, Hungary
*
Authors to whom correspondence should be addressed.
Molecules 2021, 26(1), 156; https://doi.org/10.3390/molecules26010156
Submission received: 27 October 2020 / Revised: 5 December 2020 / Accepted: 15 December 2020 / Published: 31 December 2020
(This article belongs to the Special Issue Ion Solvation)

Abstract

:
Accurate determination of the physicochemical characteristics of ionic liquids (ILs), especially viscosity, at widespread operating conditions is of a vital role for various fields. In this study, the viscosity of pure ILs is modeled using three approaches: (I) a simple group contribution method based on temperature, pressure, boiling temperature, acentric factor, molecular weight, critical temperature, critical pressure, and critical volume; (II) a model based on thermodynamic properties, pressure, and temperature; and (III) a model based on chemical structure, pressure, and temperature. Furthermore, Eyring’s absolute rate theory is used to predict viscosity based on boiling temperature and temperature. To develop Model (I), a simple correlation was applied, while for Models (II) and (III), smart approaches such as multilayer perceptron networks optimized by a Levenberg–Marquardt algorithm (MLP-LMA) and Bayesian Regularization (MLP-BR), decision tree (DT), and least square support vector machine optimized by bat algorithm (BAT-LSSVM) were utilized to establish robust and accurate predictive paradigms. These approaches were implemented using a large database consisting of 2813 experimental viscosity points from 45 different ILs under an extensive range of pressure and temperature. Afterward, the four most accurate models were selected to construct a committee machine intelligent system (CMIS). Eyring’s theory’s results to predict the viscosity demonstrated that although the theory is not precise, its simplicity is still beneficial. The proposed CMIS model provides the most precise responses with an absolute average relative deviation (AARD) of less than 4% for predicting the viscosity of ILs based on Model (II) and (III). Lastly, the applicability domain of the CMIS model and the quality of experimental data were assessed through the Leverage statistical method. It is concluded that intelligent-based predictive models are powerful alternatives for time-consuming and expensive experimental processes of the ILs viscosity measurement.

1. Introduction

The attention in green chemical technologies has resulted in the growth of a new class of highly tunable and special compounds named ionic liquids (ILs) [1]. Ionic liquids (ILs) were first introduced in 1914 by having to report the physical characteristics of ethyl-ammonium nitrate ( [ NHHH 2 ] [ NO 3 ] ) [2]. Basically, ILs are formed by the combination of organic cations and organic or inorganic anions [3]. The wide types of cations are pyridinium, phosphonium, imidazolium, and ammonium. Moreover, the commonly applied class of anions includes phosphates, halides, and sulfates [2,4]. These types of liquids can preserve their state at room temperature—which results in molten salts—and, thus, are very useful in a number of temperature-sensitive processes such as biocatalysts, etc. [5].
The features of ILs depend on a number of factors, such as the extent of the cation and anion parts, cation and anion types, and the number of branches in the elements [6]. One of the most critical properties of ILs is tuneability that allows them to acquire many of their desired features by a proper combination of anions and/or cations, side chains, and task-specific groups [6]. Ionic liquids also exhibit incredibly low volatility, hence reducing air pollution once implemented [7]. Some of the other individual features of ILs are thermal and chemical stability, nonflammability, high heat capacity, unique permittivity, high ion conductivity, and good electrical, high viscosity, electrochemical stability, ease of recycling, high solubility, and high solvability capacity for both polar and nonpolar compounds [4,7,8,9,10]. Some of the primary industrial utilization of ILs include applications in lithium batteries as eco-friendly electrolytes and dissolving various compounds of organometallic. Other broad usages of ILs touch different fields such as chemical synthesis, absorption, nanomaterials synthesis, catalytic reactions, enhanced oil recovery, bioremediation, biotechnology, and electrolytes in batteries, bio-enzyme/catalysis stabilization, and membrane separation technology [11,12,13,14,15,16,17,18,19].
The study of the chemical and physical characteristics of ILs is crucial due to their high sensitivity to the addition of small quantities of impurities such as halides or water [20,21]. These characteristics include density, speed of sound, surface tension, refractive index, electrical conductivity, and viscosity. To judge their purity, the exact measurements of these properties can always be applied [20,21]. Among the properties described above is viscosity as one of the most crucial physicochemical properties that help measure the purity of ILs [21,22,23] and help understand fluid dynamics or measuring intermolecular forces [24,25,26,27,28]. Acquiring high viscosity is required in applications such as lubricants and supported membrane separation processes, while low viscosity values are typically acceptable to increase mass transfer rates and minimize pumping costs [29,30].
In general, ILs show a wide range of viscosity from 10-105 CP, while the viscosity of ILs is two to three times higher than that of traditional organic compounds [31]. For instance, the viscosity of 1-hexyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide ( [ C 6 mim ] [ Tf 2 N ] ) is 70 CP, while the toluene viscosity at room temperature is 0.6 CP [32]. Therefore, high viscosity ILs are suitable to be used as the stationary phase of gas and liquid chromatography, lubricants, and so on. [32]. Overall, experimental measurements of thermodynamic characteristics, chemical, and physical properties of ILs, such as viscosity in an extensive range of temperature and pressure, is more important. It is also not always feasible to calculate the proprieties of ILs because they are usually costly and time-consuming [33]. Therefore, it is necessary to develop modern and reliable predictive approaches to predict the physicochemical properties of ILs [34].
Various computational methods such as group contribution methods (GCM), quantitative structure−property relationships (QSPR), and intelligent approaches (IA) can be used to predict the viscosity of ILs [35,36]. To this end, Gardas and Coutinho [37] performed a modeling investigation of viscosity of ILs by applying GCM for 500 data points from 29 ILs (based on imidazolium, pyrrolidinium, and pyridinium) in a wide range of temperature (293–393 K). The result denotes that the absolute average relative deviation (AARD) for predicting the viscosity of ILs is equal to 7.7%. Gharagheizi et al. [38] performed another study to estimate the ILs viscosity by creating a group contribution model for 443 different ILs (1672 data points) at atmospheric and a wide range of temperature from 253.15 to 433.15 K where they obtained AARD of about 6.32%. Lazzús et al. have obtained a linear model to forecast the ILs viscosity based on GCM at the temperature range of 253–395 K where AARD of the model obtained was about 4.5% [39]. A group contribution method (GCM) based on feed-forward neural network (FF-NN) to estimate the viscosity of ILs has been proposed by Paduszynski et al. [40] for 13,000 data points (1484 ILs) where pressure and temperature ranged from 0.06–350 MPa and 253–573 K, respectively, and AARD was about 11.4%. Zhao et al. [41] proposed nonlinear (support vector machine) and linear (multiple linear regression) QSPR models to model 1502 experimental data points (89 ILs) in a wide range of temperature (253.15–395.2 K) and pressure (0.1–300 MPa), where AARD for linear and nonlinear models were obtained 10.68% and 6.58%, respectively. Therefore, proposing new, accurate models to predict the viscosity of ILs yet to be accomplished.
In this paper, 2813 experimentally obtained viscosity values for 45 ILs were gathered to establish three kinds of predictive models. The first approach includes obtaining a simple correlation (Model (I)) based on temperature, pressure, molecular weight, critical volume, boiling temperature, critical temperature, critical pressure, and acentric factor. On the other hand, five advanced models with higher accuracies were developed to forecast the viscosity at different pressures and temperatures based on Model (II) (thermodynamic properties) and Model (III) (chemical structure) by using intelligent models comprising multilayer perceptron networks optimized by Levenberg–Marquardt algorithm (MLP-LMA) and Bayesian Regularization (MLP-BR), decision tree (DT), least-square support vector machine optimized by bat algorithm (BAT-LSSVM), and committee machine intelligent system (CMIS). Furthermore, Eyring’s theory (ET) was used for estimating the viscosity of ILs containing pure systems based on temperature and boiling temperature. It should be noted that various graphical and statistical criteria were considered to investigate the reliability of the proposed approaches in order to obtain the most accurate approach.

2. Viscosity Data of Ionic Liquids

The generalization and accuracy of a model highly depend on the variety and the number of data points involved in its development. For this aim, a databank—including 2813 experimental viscosity data from 45 ILs of different bases in a broad range of pressures, 0.06–298.9 (MPa), temperatures, 253.15–573 (K), and viscosities of pure ILs, 1.13–9667.6 (MPa.s)—was used to build the models [11,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61]. The ILs cations consist of pyrrolidinium [ Pyr ] + , imidazolium [ Im ] + , ammonium [ Am ] + , phosphonium [ Ph ] + , and pyridinium [ Py ] + . Moreover, anions include hexafluorophosphate [ PF 6 ] , tetrafluoroborate [ BF 4 ] , bis[(trifluoromethyl) sulfonyl]imide [ Tf 2 N ] , ethyl sulfate [ EtSO 2 ] , trifluoromethanesulfonate [ CF 3   SO 3 ] , tris(pentafluoroethyl)trifluorophosphate [ FAP ] , hydrogensulfate [ HSO 4 ] , trifluoromethanesulfonate [ TFO ] , iodid [ I ] , nitrate [ NO 3 ] , diethylphosphate [ DEP ] , dimethylphosphate [ DPO 4 ] , and methanesulfonate [ Mesy ] . The chemical structures of designated ILs are presented in Figure 1. The ionic liquids (cations and anions), pressure, temperature, and abbreviations of ILs are provided in Table 1. In addition, the full characteristics (name, unit, min, max, and mean) of the databank used for modeling in this study are shown in Table 2. In this study, the viscosity is predicted by three approaches: Model I) a simple correlation model, Model II) intelligent models based on thermodynamic properties, pressure, and temperature (Equation (1)), and Model III) intelligent models based on chemical structure, pressure, and temperature (Equation (2)).
η = f ( T ,   P , M w ,   T c , T b ,   P c ,   ω ,   V c )
η = f (T, P, Chemical structure)

3. Model Development

3.1. Calculation of Pure Viscosity Based on Eyring’s Theory-ET

Kirkwood and co-workers [62] have established a robust kinetic theory about monatomic liquids’ transport characteristics. However, this theory does not lead to quick and easy results to apply. Henry Eyring and co-workers proposed the absolute rate theory [63,64,65]. The individual molecules are in constant motion in a pure liquid at rest. However, the motion is largely confined to vibration of each molecule formed by its nearest neighbors because of the close packing inside a “cage”. This “cage” is demonstrated by an energy barrier of height Δ G 0 + ^ N A . Where N A stands for the Avogadro number (molecules/g-mol). To “escape” from the cage in the stationary fluid, a molar free energy of activation is needed that is denoted by Δ G 0 + ^ here (see Figure 2).
Based on Eyring’s theory (ET), a molecule escapes from its “cage” into an attached “hole” in a resting liquid. Thus, the molecules move in each of the direction in jumps of length “ά” at a frequency “ f ” per molecule. The rate expression determines the frequency:
f = K T p   exp ( Δ G 0 + ^ T R )
where K denotes the Boltzmann constant (J/K), P is the Planck constant (J·s), while R is the gas constant (J/mole·k). T and Δ G 0 + ^ are the absolute temperature (K) and the molar activation energy in the fluid at rest, respectively. The frequency of molecular reconfigurations is increased in a fluid flowing in the x-direction with a gradient of velocity ( d v x d y ). Figure 2 shows the potential energy barrier as distorted under the applied stress τ y x that represents as following equation:
Δ G + ^ = Δ G 0 + ^ ± ( γ / ά ) ( τ y x   Q ˜ 2 )
in which ± ( γ / ά ) ( τ y x   Q ˜ 2 ) shows the estimation of the work done on the molecules. Q ˜ denotes the volume of a mole of liquid. The frequency of forward jumps and backward jumps are “ f + ” and “ f ”, respectively. Then, the combination of Equations (3) and (4), is given as follows:
f = K T p exp ( Δ G 0 + ^ T R   ) exp ( ± γ τ y x   Q ˜ 2 ά   T R )
The net velocity that determines how far the molecules in layer “A” goes of those in layer “B” (Figure 2) is the distance traveled in each jump ( ά ) times the net frequency of advancing jumps ( f + f ). The following equation is applied:
f x A f x B = ( f + f )   ω
A linear velocity profile can be seen over a very short distance “ ά ” between the two layers, so that:
d v x d y = ( γ / ά ) ( f + f + )
Finally, Equations (5) and (7) are combined as follows:
d v x d y = ( γ / ά ) ( K T p exp ( Δ G 0 + ^ T R   ) ) ( exp ( + γ τ y x   Q ˜ 2 ά   T R ) exp ( γ τ y x   Q ˜ 2 ά   T R )                                                                       = = ( γ ά ) ( K T p exp ( Δ G 0 + ^ T R   ) ) ( 2 sinh   γ τ y x   Q ˜ 2 ά   T R )
Additionally, if γ τ y x   Q ˜ 2 ά   T R 1 , then the Taylor series can be used. Eventually, the viscosity can be obtained using the following equation:
η = ( γ ά ) 2   N A h / Q ¯ exp ( Δ G 0 + ^ T R   )
The unity factor, ( ά γ ), is a simplification that involves no loss of accuracy, as Δ G 0 + ^ is obtained empirically to make sure that the equation provides consistent values with experimental values. On the other hand, the calculated Δ G 0 + ^ (free energies of activation) through fitting Equation (9) to experimental viscosity data versus temperature are found to be nearly constant for a specific fluid and are corresponded to the internal energy of vaporization ( Δ U ^ v a p = Δ H vap-RT Δ Z vap) at the normal boiling point, as follows [63]:
Δ G 0 + ^ 0.408   Δ U ^ v a p
By using this empiricism and setting ά γ = 1, Equation (9) becomes:
η = N A p / Q ¯ exp ( 0.408 Δ U ^ v a p   T R )
The Trouton’s rule provides an accurate estimation of the energy of vaporization at the normal boiling point, as follows:
Δ U ^ v a p Δ H ^ v a p T b R 9.4 T b R
According to this approximation, Equation (11) becomes:
η = N A p / Q ¯ exp ( λ   T b   T )
In this work, in Equation (13), a “ λ ” term was added for each ILs based on GRG in Excel software. This term is not constant, while changing with each ionic liquid. Equations (11) and (13) show good consistency with the apparently successful and long-used empiricism η = A exp ( B / T ) . The theory shows a decrease in viscosity with temperature.

3.2. Generalized Reduced Gradient

The generalized reduced gradient (GRG) approach is frequently used as a solver for multivariable problems. This scheme is designed to integrate and solve linear and non-linear problems based on the concept of the reduced gradients [66]. It controls the component in such a way that when the process transitions from one stage to another, the active constraints remain to be satisfied. In fact, at a given point x, GRG provides a linear estimation for the gradient. The constraints and objective gradient are simultaneously solved, and the gradient of the objective function can be represented as constraint gradients. Then, the search area becomes smaller by moving in a practical path. Following expressions denote an objective function f(z) that is subjected to h(z) [67]:
Minimize: f(z) = z
Subjected   to :   h k ( z ) = 0
In the following form, GRG can be adjusted [67]:
d f d z k = z k t f z i t f ( h z i ) 1   h z k
Basically, f(z) will be minimum under two simple conditions that are df(z) = 0 or d f ( z ) d z k [68,69].

3.3. Decision Tree-DT

The first decision tree (DT) was the model first developed by Morgan and Sonquist [70], Automatic Interaction Detection (AID) that was introduced by Morgan and Sonquist [70]. This approach is a non-parametric supervised learning method that is applicable to both classification and regression problems. The first algorithm for the tree classification was THAID which Messenger and Mandell have suggested [71]. Learning and classification are two steps in the DT approaches. During the learning phase, the algorithm generates a tree from a set of training samples that have been classified. In the following step, unclassified data are classified using the tree developed in the learning phase [72]. The decision tree (DT) is successfully applied in many different fields such as speech recognition, remote sensing, radar signal classification, expert system, character recognition, and medical diagnosis. They are relatively inexpensive in terms of computing and appropriate accuracy. The decision tree is capable of breaking down a complicated decision-making process into a collection of simple decisions, thereby simplifying the decision [73]. A flow chart-like structure for the decision tree consists of branches, internal nodes, and root nodes. The whole of the sample space was shown via the top node with no income branch that is called the root node. The nodes with one incoming branch and more outgoing edges are classified as the internal nodes or test. The leaves or terminal nodes are identified as the other nodes that indicate the final results. The decision tree has made with three parts: pruning, stopping and splitting [74]. Splitting means that the data are divided into a number of subsets based on the most significant attribute testing that is also applicable to training instances. For the variance reduction, standard deviation reduction, and classification tree, the various criteria may be handled including the Gini index, information gain, classification error, gain ratio, and towing [75].

3.4. Multilayer Perceptron Neural Network—MLPNN

Artificial Neural Network (ANN) is modeling technique which is inspired from the human brain network as a smart computing plan. The simplest element of ANN processing is known as neurons in which connections interrelate and are organized into various layers. Neural networks are used in many aspects, including identification, estimation of functions, recognition of patterns, clustering etc. [76]. ANN has been applied in various fields, including electronics, medical, aerospace, petroleum and chemical industries [77,78]. The most popular ANN model is the multilayer perceptron (MLP). In an MLP model, there are several layers between the input and output layers, which are called hidden layers [79]. The hidden layers are specific connections between the inputs and outputs of the models. In addition, the number of neurons should be specified through a trial and error process in the hidden layer. The number of neurons in the first and last layers is regulated with input and output parameters. The use of two hidden layers in more complex problems is more appropriate than one hidden layer, but one hidden layer in the MLP model should usually be used in simple problems [80]. The value of each hidden/output neuron is calculated through the multiply and summation of the previous neurons to the neurons’ weights, and ultimately, a bias term is added to this summing-up [79]. This result is passed through the activation functions that are defined as follows:
Tangsig = Tan h : f ( n ) = e n e n e n + e n = 2 1 + e 2 n 1
ArcTan : f ( a ) = tan 1 ( a )
Logsig   =   Sigmoid : f ( y ) = 1 1 + e y
Linear = Pureline :   f ( z ) = z
Sinusid :   f ( z ) = sin ( z )
In hidden layers, logsig and tansig are the frequently applied activation functions, while pureline is usually considered the output layer’s activation function. To illustrate how the outcomes are gained from an MLP model, consider an illustrative model with two hidden layers having logsig and tansig as transfer functions, respectively, and an output layer with pureline as the transfer function. The output of the models is achieved by:
Output = pureline ( w 3 × ( logsig ( w 2 × ( tansig ( x ) + b 1 ) ) + b 2 ) + b 3 )
In the above-mentioned equation, pureline, logsig, and tansig are activation functions, respectively, b 1 denotes bias vectors of the first hidden layer and b 2 , b3 are bias vectors of the second hidden layer and the output layer, respectively. Additionally, w 1 , w 2 , and w 3 are the matrix weight of the first layer, the matrix weight of the second layer, and the matrix weight of the output layer, respectively [80].
The optimization algorithm is one of the most key roles in the performance of the MLP model which are used for training the model. Thus, two main optimization algorithms have been used in this study, including LMA, Levenberg-Marquardt, and BR, Bayesian Regularization. Additional information on developing LMA and BR in the MLP training phase can be found elsewhere [81,82,83,84]. A scheme of the MLP network that was used in this study is represented in Figure S1 (in the Supplementary File).

3.5. Least Square Support Vector Machine—LSSVM

A support vector machine (SVM) is a supervised and powerful, intelligent tool applied on known input/output data for various purposes such as pattern recognition, problem classification, and regression analysis [85,86]. The least-square SVM (LSSVM) was proposed by Vandewalle and Suykens [87] as a newer version of the support vector machine. Further, the LSSVM has emerged to improve and prevail the typical shortcomings of the SVM approach and simplify its solution [85]. The LSSVM updates the optimization constraints and mathematically determines the regression error. In reality, regression error in SVM algorithms during the learning process is optimized, and it is numerically defined and resolved in LSSVM approaches [85]. The penalized function of the least square SVM (LSSVM) method is described as shown below [85,88,89,90]:
Q L S S V M = λ × k = 1 N e k 2 + w T 2
where λ is the summation of regression errors, while T shows the transpose matrix. Equation (23) is subjected to the following constraint [85]:
y k = e k + b + w T φ ( x k ) ,   k = 1 ,   2     N
where e k denotes the regression error of N training objects, y shows the output vector of the model, b is the intercept or the bias of linear regression. Additionally, T expresses the transposed matrix, and w represents linear regression slope. The weight function (w) is expressed as follows [85]:
w = k = 1 N x k × a k ,   where ,   a k = 2 × e k × λ
Where a k denotes the Lagrangian multiplier. λ corresponds to the relation of specific weights and also the weight of all regression errors, while e k represents regression error related to whole databank. Using least square SVM (LSSVM) method, another form of Equation (25) is represented as follows [85,87,88,89,90]:
y = b + k = 1 N x k T × x × a k
Consequently, the Lagrange multipliers are described as the following equation [85]:
a k = ( y k b ) ( 2 λ ) + x × x k T
Using the following Kernel function, the aforementioned equation for linear regression can be rewritten as [85]:
f ( x ) = b + k = 1 N a k × k ( x , x k )
In which k ( x , x k ) shows the Kernel function. k ( x , x k ) denotes the dot product of Φ ( x k ) and Φ ( x ) T as follows [85]:
K ( x , x k ) = Φ ( x k ) × Φ ( x ) T
In the present study, one of the most well-known Kernel functions, i.e., the radial basis Kernel function, was applied. This latter is defined as follows:
K ( x , x k ) = exp ( | | x k x | | 2 σ 2 )
In the above equation, σ 2 is a regulation parameter that is to be obtained by an optimization algorithm. Briefly, there are two adjustable parameters in the LSSVM technique, namely λ and σ 2 , which should be optimized during training step. The parameters for tuning were optimized using the bat (BAT) algorithm. This approach was named BAT-LSSVM. In this approach, σ 2 and λ are optimized using BAT algorithm for determining viscosity of ILs. The schematic representation of this approach for predicting viscosity of ILs is showed in Figure S2 (in the Supplementary File).

3.6. Committee Machine Intelligent System (CMIS)

In an extensive number of studies, various artificial intelligence modelling technique are employed, and after selecting the most accurate model, the other approaches are left out. The better choice is to utilize these models to build a more precise model. Many years ago, Nilsson presented a Committee Machine Intelligent System (CMIS) as a kind of artificial neural network. Through this process, the basic goal is to divide and conquer until resolving a problem. In addition, the outputs of each approach have been mixed and the benefits of all work with little extra calculation have been achieved. Afterward, the model can outperform the best single network [91,92]. As a matter of fact, in this technique, different predictive models are combined in order to form a more efficient and accurate model. Different methods have emerged for merging several models in a single model. The categories of committee machines can be divided into the two following types [93]:
(1)
Static structure
(2)
Dynamic structure
An appropriate approach based on simple averaging or weighted averaging is employed to linearly combine the resolutions and getting the best model [94]. Herein, an improved weighted averaging approach was used and a bias factor was added to the linear equation. Each model’s contribution extent in a CMIS corresponds to the absolute coefficient of that model in the CMIS linear equation. Figure 3 shows a scheme of the CMIS approach in this study.

3.7. Optimization Technique

Bat Algorithm (BAT)

In 2010, Yang Xinshe introduced the bat algorithm (BA) by analyzing the behaviors and the features of the microbat [95]. In the rest of this section, we shed light on the echolocation and details of this nature-inspired algorithm. Yang and Deb [96] used the idealized and approximate rules in the BA approaches with the following steps:
  • All the species of the bat utilize echolocation to sense distance, and bats ‘know’ the discrepancy among food/prey and background obstacles in some magical techniques.
  • In order to search prey, the bats can fly fortuitously with the velocity v i at position x i with a frequency f min , loudness A 0 , and a variable wavelength λ . Bats can spontaneously adjust the wavelength and/or frequency of their generated pulses and regulate the level of pulse emission r in the range of [0,1], reliant on the nearness of their goal.
  • Although there are various methods to regulate the loudness, it is usually assumed that the loudness is between a positive A 0 and a minimum constant amount, which is represented by A min .
The following equations show the motion of a virtual bat by Yang’s method [97]:
f i = f m i n ( f m i n f m a x ) . β
v i t = v i t 1 ( x b e s t x i t ) . f i
x i t = v i t + x i t 1
where β is between 0 and 1 and denotes a random obtained from monotonic, min and max are the suffixes that are shown minimum and maximum value, respectively. The velocity of the bat is v i , the current iteration is indicated as “ t ”. The global close-best solution discovered so far over the whole population is denoted x b e s t . The location of the ith bat in the solution space is x i and f denotes the frequency used by the bat to seek for its hunt. Furthermore, one of the roles in the system is also assumed to be the rate of pulse emission from the bat. The symbol r i [ 0 , 1 ] is the pulse emission rate, also, the suffix i denotes the ith bat. A random number is created in each iteration and compared with r i . The random walk is denoted a local search strategy, if r i is smaller than a random walk. The below-shown equation represents a new solution for the bat:
x n e w = r a n d . A t + x o l d
where the average loudness of all bats at the current time step is A t A i t , while rand [ 1 , 1 ] . The loudness A i and pulse emission rate r i are also updated after updating the positions of the bat, only when the close-best global solution is updated and A i is bigger than the random number. Equations (35) and (36) state the update of r i and A i [97]:
A i t + 1 = A i t . β
r i t + 1 = [ exp ( α t ) + 1 ] r i 0
where β and α are constants. As a result, 0 < β < 1 and α > 0 [97].
The above-described steps are reiterated until reaching a stopping condition.

4. Model Assessment

4.1. Statistical Criteria

In order to evaluate the validity of the obtained models, the mathematical formula for statistical assessment parameters including average relative deviation (ARD%), determination coefficient ( R 2 ), standard deviation (SD), root mean square error (RMSE), and average absolute relative deviation (AARD%) were used. These statistical parameters are detailed as follows:

4.1.1. Determination Coefficient ( R 2 )

Regression determination is a measure of fitting quality revealing the accuracy of the model. Accordingly, if its value is close to 1, the model matches the data more accurately. The mathematical formula of R 2 is given below:
R 2 = i = 1 N P ( η i e x p η ¯ ) 2 i = 1 N P ( η i c a l η i e x p ) 2 i = 1 N P ( η i e x p η ¯ ) 2

4.1.2. Average Relative Deviation (ARD%)

The relative deviation of the estimated data from the experimental is measured by the percentage of ARD with the following equation:
ARD % = 100 N P j = 1 ( η j e x p η j e s t η j e x p )
Positive and negative ARD (%) denote underestimation and overestimation of a model, respectively.

4.1.3. Standard Deviation (SD)

It is a measure of the scattering of data around the mean. It is defined as:
SD = ( 1 N 1 j = 1 N P ( η j e x p η j e s t λ j e x p ) 2 ) 1 2

4.1.4. Average Absolute Relative Deviation (AARD%)

The AARD value is a measure of the relative absolute deviation of the predicted/represented data from the actual/real data. It is represented by:
AARD ( % ) : 100 × j = 1 N P | η j e x p η j e s t η j e x p | N P

4.1.5. Root Mean Square Error (RMSE)

The root-mean-square error (RMSE) is a widely used calculation of the differences between values (sample or population values) expected by a model or estimator and the observed values. It is denoted by:
RMSD = 1 N P j = 1 N P ( η j e x p η j e s t ) 2
where N P points out the numbers of data points and η is the experimental/real values of viscosity of ILs, while η ¯ , superscripts “est”, and “exp” are the average of viscosity of ILs obtained by experiments/real data, the estimated value, and the experimental/real value, respectively.

4.2. Graphical Evaluation of the Models

To better assess the proposed models and analyze their predictive accuracy, several graphical plots were employed in this study. The plots include cross-plot, cumulative frequency diagram, and error distribution diagram. In the error distribution, the percent relative deviation is plotted versus the target/real values to measure the distribution of error around the zero-error line and to show whether the model has an error trend or not. Cross-plot includes sketching the estimated/represented value obtained by the model versus the experimental data. Then, between the experimental data and predicted/represented values, a straight line of a 45° (unit slope line) is drawn. Eventually, the closer the plotted data to this line, the higher is the model accuracy. For a cumulative frequency diagram, the larger part of the approximations will be in a standard error range where the cumulative frequency is calculated from the absolute relative error.

5. Result and Discussion

5.1. Development of Models

In this study, models were developed based on 2813 data points (45 ionic liquids) that were gathered from the literature. The random databank division was conducted to attain two subsets, the test (20% of the whole dataset) and training sets (80% of the whole dataset). In fact, the validity of the constructed models and over-fitting problems were monitored based on the “testing” subset (563 data points). The models’ structure and the adjustment of their tuning parameters were performed based on the “training” subset (2250 data points) [98].
As it was mentioned early, T, P, M w , V c , T b , T c , P c , ω , and chemical structure were considered as the input parameters, while the output was the viscosity of ILs (MPa·s). Firstly, a new empirical correlation was established based on the GRG approach. The following relation was obtained for viscosity of ILs:
log ( η ) = ( a T i + b P j + c M w k + d V C l + e T b m + f T c n + g P c o + h ω p ) T q
where T c and P c denote the critical temperature and the critical pressure of ionic liquids, respectively. Additionally, ω , T, and P are the acentric factor, temperature, and pressure, respectively, while M w is the molecular weight of ILs, V c is the critical volume of ILs, and T b is the boiling temperature of the ILs. Moreover, the other parameters are the adjustable coefficients of the correlation which are shown in Table S1 (in the Supplementary File). The ARD%, AARD%, SD, RMSE, and R 2 for the proposed correlation are 13.44%, 28.34%, 0.368, 394.12, and 0.207, respectively. The cross-plots between the logarithm of predicted data were plotted against the logarithm of experimental viscosity data for the presented correlation (Equation (42)) is shown in Figure S3 (in the Supplementary File). This plot depicts a medium-uniform distribution of the predictions around the unit-slope line. On the other hand, Eyring’s theory (ET) was applied to estimate the viscosity of ILs based on boiling temperature ( T b ) and temperature (T) for the gathered database. The achieved correlation is expressed as follows:
η = N A p / Q ¯ exp ( λ   T b   T )
where η is the predicted viscosity of pure ILs (MPa·s). N A and p are the Avogadro number ( mole 1 ) and the Plank constant (J·s), respectively. Q ¯ denotes the volume of a mole of liquid (m3/mole), T b and T represent the boiling temperature (K) and temperature (K), respectively, while λ is a factor that was obtained based on GRG and for each ILs is not constant. The average absolute relative deviation (AARD) is 21.86%. Additionally, Figure 4 represents a cross-plot between the logarithm of predicted of ILs viscosity versus the logarithm of experimental of ILs viscosity. As seen, a moderate conformity was noticed as the data points were not very close to the diagonal line. The theory of Eyring is not sufficient in this study, because Arrhenius dependency does not match the experimental transport features of ILs. Viscosity values reduce significantly with temperature rises, as generally reported for all ILs measured. Furthermore, the thermally enabled transport features of ILs are generally defined by the Vogel–Tamman–Fulcher (VTF) development because of the development of an underlying complex energy landscape with a multiplicity of local potential energy minimums and a broad dispersion of energy barriers [99,100,101,102].
Further, DT, MLP and LSSVM-BAT modelling techniques were employed to predict the viscosity of ILs based on Model (II) and Model (III). In the first stage based on Model (II), the MLP training phase was performed utilizing two learning techniques, namely, BR and LMA. For the best architecture of the MLP network with BR and LMA optimizers, using three hidden layers were determined to be the most appropriate for predicting the viscosity of ILs. Therefore, each layer for MLP-LMA and MLP-BR encompasses 11-11-9 and 13-11-9 neurons, respectively. The best transfer functions for all of three layers of MLP-LMA and MLP-BR models were found as tansig. The optimum values of the main parameters of the developed LSSVM model, i.e., σ 2 and λ, were investigated using the BAT algorithm. For Model (II), the obtained values for σ 2 and λ are 25.4817 and 5,793,328.591, respectively. Table 3 reports the statistical assessment of the constructed schemes.
According to the results presented in the Table 3, the MLP-BR model is the most reliable technique compared to the other models since it provides accurate predictions for the whole dataset. The proposed models including MLP-LMA, MLP-BR, LSSVM-BAT, and DT were combined into a CMIS. The coefficients of the constructed CMIS model were optimized through a multiple linear regression and the following equation was obtained for this latter:
η I L s i = a + b y D T i + c y L S S V M - B A T i + d y M L P - B R i + e y M L P - L M A i
where a, b, c, d, and e were determined to be 0, 8.77 × 10 5 , 0.049344, 0.619175, and 0.330275, respectively.
In the other procedure based on Model (III) (chemical structure, pressure, and temperature) the MLP models with two optimization algorithms, including LMA and BR were developed. The hidden layers for MLP-LMA and MLP-BR were chosen 12-11-9 and 11-11-9, respectively, where the first, second, and third numbers correspond to the number of neurons of the first, second, and third hidden layer, respectively. Additionally, tansig was found to be the best activation function for the three hidden layers of MLP-LMA and MLP-BR models. The tuning parameters (namely λ and σ 2 ) are the most important in the LSSVM-BAT model. The parameters of the LSSVM method were determined as 6,635,314.0686 and 27.3672 for λ and σ 2 , respectively, using BAT optimization algorithm. Table 4 shows the statistical parameters for predicting the viscosity of ILs based on Model (III). Accordingly, MLP with BR optimization is a better choice for viscosity estimation based on Model (III). Therefore, a CMIS method that combines the proposed paradigms, i.e., MLP-LMA, MLP-BR, LSSVM -BAT, and DT, was implemented. The appropriate coefficients of the CMIS model were specified employing a multiple linear regression. For this goal, the following correlation is provided:
η I L s i = q + m y D T i + n y L S S V M - B A T i + d y M L P - B R i + e y M L P - L M A i
where the coefficients q to e are as shown below:
q = 0.052; m = 0.0001; n = 0.151; d = 0.466; e = 0.379.

5.2. Statistical Evaluation

In order to indicate the error margin, various statistical indexes, including ARD, AARD, RMSE, SD, and R 2 were computed for each of the developed model. As reported in Table 3, the computed AARD, ARD, SD, and RMSE for the CMIS are small values for both test and training sets, as well as the whole dataset. The value of AARD for CMIS model based on Model (II) is 3.24%, which shows the higher accuracy of this model compared to the four ANN-based models. Briefly, Table 3 depicted that the models are ranked as follows in terms of accuracy:
CMIS > MLP-BR > MLP-LM > LSSVM-BAT > DT
In the other hand, based on Model (III), the same statistical parameters comprising ARD, AARD, RMSE, SD, and R 2 are also presented for MLP-BR, MLP-LMA, LSSVM-BAT, DT, and CMIS which are summarized in Table 4. The models are categorized according to their precision as follows: DT < LSSVM-BAT < MLP-LMA < MLP-BR < CMIS. As can be deduced, the CMIS model could provide more accurate predictions compared to the other established models.

5.3. Graphical Error Analysis

Various graphical error evaluations based on Models (II) and (III) were investigated to give a more transparent view of the performance of the models. To assess the accuracy of Model (II), the logarithm of predicted viscosity data was plotted against the logarithm of corresponding experimental values in Figure 5. Therefore, Figure 5 shows comparative cross-plots for test and training sets for CMIS model (cross-plot for other models such as: MLP-BR, MLP-LMA, DT, and LSSVM-BAT are presented in Figure S4 in the supplementary File). A high-dense accumulation of the points can be seen around the unit slop for both test and training datasets. On the other hand, higher accuracy of the constructed CMIS compared to the other models for approach (II) is clearly shown in these figures. Accordingly, the CMIS paradigm is introduced as the best predictive method. Figure 6 represents the relative deviation (RD) against the experimental viscosity of ILs for the CMIS model. Additionally, Figure S5 (in the Supplementary File) reveals the experimental viscosity of ILs against the relative deviation for MLP-BR, MLP-LMA, LSSVM-BAT, and DT. As shown in these figures, for LSSVM-BAT and DT, the data distribution near the zero line is larger than the developed CMIS model. This observation confirms the results reported in Table 3. The implementation of the CMIS results in the minimum total relative deviation guiding to the narrowest error margin.
Furthermore, Figure 7 is plotted to assess the capability of the models with respect to the AARD (%) of the proposed empirical correlation (Equation (42)), smart models developed based on Model (II), Eyring’s theory in this study, and also other models form the literature. Clearly, the CMIS model and then the MLP-BR approach are more accurate, have a higher flexibility, and are more suitable to be used for the viscosity prediction of ILs approach. Figure 8 depicts the data cumulative frequency versus the ARD for the MLP-LMA, MLP-BR, LSSVM-BAT, DT, and CMIS models. As shown in this figure, the CMIS approach is more appropriate than other models and has a high accuracy in predicting viscosity of ILs. The figure implies that employing the CMIS model about 80% of the ILs viscosity data were predicted with an ARD of less than 4.5% and more than 72% of them were predicted with an ARD of less than 2.3%.
Furthermore, the so-called Taylor diagram [103] is represented in Figure 9 and Figure 10 based on Model (II) to provide a more comprehensive presentation. This diagram considers various statistical parameters. As shown in Figure 9 and Figure 10, standard deviations and correlation coefficients of each model (DT, LSSVM-BAT, MLP-LMA, MLP-BR, and CMIS) are employed to quantify the variance between the predicted and measure viscosities of the ILs. RMSE is another measure considered as red circles in this diagram. An ideal predictive model is a model that has an R 2 equal to 1 in the Taylor diagram. The CMIS model exhibited the outperformance with RMSE values of 9.553 and 9.035 for the training and testing phases, respectively, while the other models showed RMSE values of more than 10. Thus, the Taylor diagram approves the outperformance of the CMIS once again as its predictions are the closest to the experimental measurements. In conclusion, the proposed model is very advantageous since it is developed based on a large dataset and the employment of various models enhances its exactness and consistency.
On the other hand, graphical evaluation was conducted for the proposed paradigms based on Model (III) to investigate their accuracy and effectiveness. Figure 11 and Figure S6 (in the supplementary File) illustrate the logarithm of measured viscosity values versus the logarithm of estimated ones for CMIS approach and other smart models, respectively. These figures reflect that data points have formed a compact zone in the proximity of the unit slope line (45° line). In other words, the figures reveal that all intelligent models could provide consistent responses with target values. Furthermore, despite the high performance of all developed models, the constructed CMIS model could provide more reliable results thank to the satisfactory statistical criteria. Figure 12 and Figure S7 (in the Supplementary file) represent a comparison between the developed models in the regard of the relative deviations of models. Obviously, the dense concentration of points around the zero line reflects that the CMIS model is able to predict the viscosity of ILs with the lowest possible relative errors. Cumulative frequency plots are comparative plots, which are usually represented to simply clarify the competency of the models. This convenient plot has been sketched in Figure 13. As shown in this diagram, the CMIS shows the largest cumulative frequency at a specific absolute relative deviation (ARD). To clarify this, about 95% of data points were predicted by the CMIS model when we approach to the ARD of 10%, while the corresponding values for MLP-BR, MLP-LMA, LSSVM-BAT, and DT models were 88%, 83%, 80%, and 60%, respectively. Therefore, the superiority of the CMIS is once again demonstrated. In addition to the aforementioned point of Figure 13, this figure also demonstrates that the CMIS model could predict about 100% of data with lower ARD compared to the other models. Figure 14 shows the visual comparison between AARD (%) for different approaches including DT, MLP-LMA, MLP-BR, LSSVM-BAT, and CMIS for estimating the viscosity of ILs based on Model (III). As can be seen, the accuracy of CMIS-based models is generally more than other smart models.
In order to investigate the strength and reliability of the proposed models with regard to temperature alterations, Model (II) and Model (III) were used. Figure 15a,b show the effect of temperature variations at a constant pressure (75 MPa) on 1-ethyl-3-methylimidazolium ethylsulfate based Model (II) and (III), respectively. Obviously, the predicted viscosity of ILs by all of the smart models (MLP-BR, MLP-LMA, LSSVM-BAT, and CMIS) except the DT model is in agreement to the experimental dataset. According to Figure 15, the DT approaches could not provide consistent estimations for the viscosity of ILs and the predicted values are associated with significant errors. When the temperature varies, the CMIS, MLP-BR, MLP-LMA, and LSSVM-BAT paradigms show similar physical trends. Figure 15a,b indicate a reduction in the viscosity of ILs when the temperature rises denoting an inverse relationship between the viscosity and temperature. The discussed results confirm the expected trend as well as the common knowledge of chemical thermodynamics.
As shown in Table 3 and Table 4, among the best CMIS models for Model (II) and Model (III), the CMIS Model (II) has the lowest AARD value of 3.293%; hence, it was kept for further evaluations.

5.4. Identifying Outliers in Experimental Data and Applicability Domain of CMIS Model

The objective of outlier (or aberrant) detection is found for groups of data (or individual data) which deviate a lot from the bulk data in a database [104]. For this aim, the Leverage approach is one of the well-established strategies [35,104,105]. In this technique, the Hat matrix (H) and standardized residuals (R) are two main concepts [105]. The standardized residual (R) is calculated for each data point based on following equation:
R i = z i ( M S E ( 1 H i i ) ) 1 2
where MSE denotes the model mean square error, while z i and H i i represent the error and the Hat indices (Leverage) of the ith data point [106]. In addition, Hat index (or Leverage) can be determined as below:
H = X ( X t X ) 1   X t
where X t points out the transpose matrix, X represents a two-dimensional q × w matrix (where, “q” and “w” are the quantity of data points and the dimension of the model, respectively).
The Williams plot was applied in order to investigate the outliers after the measurement of the R and Hat index (H) [107,108]. In this diagram, the standardized residuals are sketched against Hats. The Leverage limit ( H * ) is a parameter defined as 3 a b , where b stands for the quantity of data points and a is the number of model parameters plus one [109]. Standardized residual (R = 3) is an acceptable data point within the range of [−3,+3] of standard deviations from the average to cover 99% of normally distributed data [35,105]. The model is statistically applicable if the large percentage of data points are within the H * H 0 and 3 R 3 [35]. Usually, data points within the range of −3 ≤ R ≤ 3 and H * ≤ H are regarded as “Good High Leverage” points, and are not in the applicability domain but are predicted well. By contrast, the data points with R values of higher than 3 or lower than −3 are considered as “Bad High Leverage”. These are the points predicted with large uncertainty and are located outside of the model applicability domain. These are the experimentally suspected data points.
In the present study, it is clear that the performance of the CMIS as the best model based on Model (II) (thermodynamic properties, pressure, and temperature) can be significantly influenced by the consistency of the employed data. The proposed model showed a H * value of 0.0095. The Williams plot of the CMIS model is depicted in Figure 16. Accordingly, nearly all data points appear to be within 0 ≤ H ≤ 0.0095 and −3 ≤ R ≤ 3 reflecting that the developed CMIS system is statistically applicable. Commonly, the less standardized residual values of experimental data, the more reliable they are. Nevertheless, with respect to Figure 16, 0.63% of all data (18 suspected data points) are in the ranges of R < −3 or R > 3 and consequently, they are regarded as outliers that are associated with large uncertainty. Furthermore, 2.73% of all data (77 data points) have H > 0.0095. These are the points that are out of Leverage, no matter what their Hat’s (Leverage) value is, all of them are located in the range of −3 ≤ R ≤ 3, and consequently, they are Good High Leverage.

5.5. Relative Importance of Input Variables

Sensitivity analysis was performed to determine the magnitude of the impact of all input variables on the viscosity of ILs using the CMIS model. The input parameters are based on Model (II). Figure 17 shows the relative importance of the inputs. The relevancy factor (r) is a measure that determines the extent of the impact of each input parameter on the pure viscosity of ILs as the model output. Positive values in this technique demonstrate a direct relationship between the output and the corresponding input parameter, while negative values reveal an inverse one. The values of the relevancy factor (or r) are evaluated as the following equation:
r ( I i , η ) = j = 1 n ( I i j I ¯ i ) ( η j η ¯ ) ( j = 1 n ( I i , j I ¯ j ) 2 j = 1 n ( η j η ¯ ) 2 ) 0.5
where n, I i , j , and I ¯ i are the number of datasets, the jth value, and mean of the i -th input variable, respectively. η ¯ denotes the average value of the estimated/represented viscosity of ILs, while η j expresses the jth value of the represented/predicted viscosity of ILs. Figure 17 depicts each parameter’s relative effect based on Model (II) on the pure viscosity of ionic liquids (as the output). As can be seen, the input parameters, namely, acentric factor and pressure, positively affect the output model. Any increase in pressure and acentric factor cause to an increase in the viscosity of pure ionic liquids. As is obvious in Figure 17, other parameters such as T, M w , V c , T b , T c , and P c show a decreasing trend due to their negative relevancy factors, which means that increasing these parameters would decrease in viscosity of ILs. Obviously, the viscosity of ILs is mainly influenced by temperature.

6. Conclusions

In this study, 2813 experimental data from 45 ILs were initially modeled with the aim of obtaining a simple correlation based on temperature, pressure, molecular weight, critical volume, boiling temperature, critical temperature, critical pressure, and acentric factor of ILs. In addition, Eyring’s theory was applied based on temperature and boiling temperature to estimate the viscosity of ILs. In another scheme, the proposed models were linked to a univalent model as CMIS. To study the efficiency of the model, a comparison between the results of CMIS and experimental data was made using both statistical and graphical methods. The model showed a stable performance and high accuracy based on R2, RMSE, AARD, and ARD definitions. The effects of input parameters on model outputs were also analyzed using the relevancy factor. Finally, the leverage statistical approach was used to assess the reliability and validity of the employed dataset. In this regard, the Williams’s plot was employed to investigate the applicability domain of the constructed paradigm and find the recorded data. The results indicate that just a few data points were outside the applicability domain. It can be concluded that the developed CMIS can be used as a valuable tool to predict the viscosity of ILs with high accuracy as well as acquiring accurate knowledge of IL’s physicochemical properties in various chemical engineering processes where needed.

Supplementary Materials

The following are available online, Figure S1: The structure of MLPNN used in this study, Figure S2: Schematic of procedure for development of BAT-LSSVM model, Figure S3: Cross plot of the proposed correlation for viscosity of ILs, Figure S4: Logarithm of experimental viscosity data against predicted values based on Model (II), Figure S5: Relative deviation of predictions of various models versus logarithm of experimental data based on Model (II), Figure S6: Cross plot of the proposed developed models based on Artificial Neural Network based on Model (III), Figure S7: Relative deviation distribution for developed models in this study for estimation of the viscosity of ILs based on Model (III). Table S1: The adjustable parameters to proposed correlation in this study.

Author Contributions

Conceptualization, S.P.M. and S.A.; Methodology, S.P.M. and S.A.; Soft-ware, S.P.M., S.A., M.N.A., A.H.-S., A.M. (Ahmad Mohaddespour) and A.M. (Amir Mosavi); Validation, S.P.M., S.A., M.N.A., A.H.-S., A.M. (Ahmad Mohaddespour) and A.M. (Amir Mosavi); Formal analysis, S.P.M., S.A., M.N.A. and A.H.-S.; Investigation, S.P.M., S.A. and A.H.-S.; Resources, A.H.-S., A.M. (Ahmad Mohaddespour) and A.M. (Amir Mosavi); Data curation, S.P.M. and S.A.; Writing—original draft preparation, S.P.M.; Writing—review and editing, S.P.M., S.A., M.N.A., A.H.-S., A.M. (Ahmad Mohaddespour) and A.M. (Amir Mosavi); Visualization, S.P.M.; Supervision, A.H.-S. and A.M. (Ahmad Mohaddespour); Project administration, A.H.S.; Funding acquisition, A.M. (Amir Mosavi). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available in Supplementary Material.

Acknowledgments

The support of Alexander von Humboldt Foundation is acknowledged.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Salgado, J.; Regueira, T.; Lugo, L.; Vijande, J.; Fernández, J.; Garcia, J. Density and viscosity of three (2, 2, 2-trifluoroethanol+ 1-butyl-3-methylimidazolium) ionic liquid binary systems. J. Chem. Thermodyn. 2014, 70, 101–110. [Google Scholar] [CrossRef]
  2. Zafarani-Moattar, M.T.; Majdan-Cegincara, R. Viscosity, Density, Speed of Sound, and Refractive Index of Binary Mixtures of Organic Solvent + Ionic Liquid, 1-Butyl-3-methylimidazolium Hexafluorophosphate at 298.15 K. J. Chem. Eng. Data 2007, 52, 2359–2364. [Google Scholar] [CrossRef]
  3. Atashrouz, S.; Zarghampour, M.; Abdolrahimi, S.; Pazuki, G.; Nasernejad, B. Estimation of the Viscosity of Ionic Liquids Containing Binary Mixtures Based on the Eyring’s Theory and a Modified Gibbs Energy Model. J. Chem. Eng. Data 2014, 59, 3691–3704. [Google Scholar] [CrossRef]
  4. Schmidt, H.R.; Stephan, M.; Safarov, J.; Kul, I.; Nocke, J.; Abdulagatov, I.; Hassel, E. Experimental study of the density and viscosity of 1-ethyl-3-methylimidazolium ethyl sulfate. J. Chem. Thermodyn. 2012, 47, 68–75. [Google Scholar] [CrossRef]
  5. Van Rantwijk, F.; Sheldon, R.A. Biocatalysis in Ionic Liquids. Chem. Rev. 2007, 107, 2757–2785. [Google Scholar] [CrossRef]
  6. Freemantle, M. An Introduction to Ionic Liquids; Royal Society of Chemistry: London, UK, 2010. [Google Scholar]
  7. Plechkova, N.V.; Seddon, K.R. Applications of ionic liquids in the chemical industry. Chem. Soc. Rev. 2008, 37, 123–150. [Google Scholar] [CrossRef]
  8. Wasserscheid, P.; van Hal, R.; Bösmann, A. 1-n-Butyl-3-methylimidazolium ([bmim]) octylsulfate—an even ‘greener’ ionic liquid. Green Chem. 2002, 4, 400–404. [Google Scholar] [CrossRef]
  9. Ghatee, M.H.; Zare, M.; Zolghadr, A.R.R.; Moosavi, F. Temperature dependence of viscosity and relation with the surface tension of ionic liquids. Fluid Phase Equilibria 2010, 291, 188–194. [Google Scholar] [CrossRef]
  10. Chiappe, C.; Pieraccini, D. Ionic liquids: Solvent properties and organic reactivity. J. Phys. Org. Chem. 2005, 18, 275–297. [Google Scholar] [CrossRef]
  11. Wu, T.-Y.; Chen, B.-K.; Hao, L.; Kuo, C.-W.; Sun, I.-W. Thermophysical properties of binary mixtures $1−methyl−3−pentylimidazoliumtetrafluoroborate+polyethyleneglycolmethylether$. J. Taiwan Inst. Chem. Eng. 2012, 43, 313–321. [Google Scholar] [CrossRef]
  12. Jiang, H.; Wang, J.; Zhao, F.; Qi, G.; Hu, Y. Volumetric and surface properties of pure ionic liquid n-octyl-pyridinium nitrate and its binary mixture with alcohol. J. Chem. Thermodyn. 2012, 47, 203–208. [Google Scholar] [CrossRef]
  13. Wu, T.-Y.; Chen, B.-K.; Hao, L.; Lin, K.-F.; Sun, I.-W. Thermophysical properties of a room temperature ionic liquid (1-methyl-3-pentyl-imidazolium hexafluorophosphate) with poly(ethylene glycol). J. Taiwan Inst. Chem. Eng. 2011, 42, 914–921. [Google Scholar] [CrossRef]
  14. Lopes, J.N.C.; Gomes, M.F.C.; Husson, P.; Pádua, A.A.H.; Rebelo, L.P.N.; Sarraute, S.; Tariq, M. Polarity, Viscosity, and Ionic Conductivity of Liquid Mixtures Containing [C4C1im][Ntf2] and a Molecular Component. J. Phys. Chem. B 2011, 115, 6088–6099. [Google Scholar] [CrossRef] [PubMed]
  15. Hezave, A.Z.; Dorostkar, S.; Ayatollahi, S.; Nabipour, M.; Hemmateenejad, B. Dynamic interfacial tension behavior between heavy crude oil and ionic liquid solution (1-dodecyl-3-methylimidazolium chloride ([C12mim][Cl]+distilled or saline water/heavy crude oil)) as a new surfactant. J. Mol. Liq. 2013, 187, 83–89. [Google Scholar] [CrossRef]
  16. Volkl, J.; Muller, K.J.; Mokrushina, L.; Arlt, W. A Priori Property Estimation of Physical and Reactive CO2 Absorbents. Chem. Eng. Technol. 2012, 35, 579–583. [Google Scholar] [CrossRef]
  17. Królikowska, M.; Hofman, T. Densities, isobaric expansivities and isothermal compressibilities of the thiocyanate-based ionic liquids at temperatures (298.15–338.15K) and pressures up to 10MPa. Thermochim. Acta 2012, 530, 1–6. [Google Scholar] [CrossRef]
  18. Hezave, A.Z.; Dorostkar, S.; Ayatollahi, S.; Nabipour, M.; Hemmateenejad, B. Investigating the effect of ionic liquid (1-dodecyl-3-methylimidazolium chloride ([C12mim] [Cl])) on the water/oil interfacial tension as a novel surfactant. Colloids Surf. A Physicochem. Eng. Asp. 2013, 421, 63–71. [Google Scholar] [CrossRef]
  19. Ciocirlan, O.; Croitoru, O.; Iulian, O. Densities and Viscosities for Binary Mixtures of 1-Butyl-3-Methylimidazolium Tetrafluoroborate Ionic Liquid with Molecular Solvents. J. Chem. Eng. Data 2011, 56, 1526–1534. [Google Scholar] [CrossRef]
  20. Torrecilla, J.S.; Tortuero, C.; Cancilla, J.C.; Díaz-Rodríguez, P. Neural networks to estimate the water content of imidazolium-based ionic liquids using their refractive indices. Talanta 2013, 116, 122–126. [Google Scholar] [CrossRef]
  21. Torrecilla, J.S.; Tortuero, C.; Cancilla, J.C.; Díaz-Rodríguez, P. Estimation with neural networks of the water content in imidazolium-based ionic liquids using their experimental density and viscosity values. Talanta 2013, 113, 93–98. [Google Scholar] [CrossRef]
  22. Fernández, A.; Garcia, J.; Torrecilla, J.S.; Oliet, M.; Rodriguez, F. Volumetric, transport and surface properties of [bmim][MeSO4] and [emim][EtSO4] ionic liquids as a function of temperature. J. Chem. Eng. Data 2008, 53, 1518–1522. [Google Scholar] [CrossRef]
  23. Nieuwenhuyzen, M.; Seddon, K.R.; Teixidor, F.; Puga, A.V.; Viñas, C. Ionic Liquids Containing Boron Cluster Anions. Inorg. Chem. 2009, 48, 889–901. [Google Scholar] [CrossRef] [PubMed]
  24. Zhou, Q.; Wang, L.-S.; Chen, H.-P. Densities and Viscosities of 1-Butyl-3-methylimidazolium Tetrafluoroborate + H2O Binary Mixtures from (303.15 to 353.15) K. J. Chem. Eng. Data 2006, 51, 905–908. [Google Scholar] [CrossRef]
  25. Zhu, A.; Wang, J.; Liu, R. A volumetric and viscosity study for the binary mixtures of 1-hexyl-3-methylimidazolium tetrafluoroborate with some molecular solvents. J. Chem. Thermodyn. 2011, 43, 796–799. [Google Scholar] [CrossRef]
  26. Gong, Y.-H.; Shen, C.; Lu, Y.-Z.; Meng, H.; Li, C.-X. Viscosity and Density Measurements for Six Binary Mixtures of Water (Methanol or Ethanol) with an Ionic Liquid ([BMIM][DMP] or [EMIM][DMP]) at Atmospheric Pressure in the Temperature Range of (293.15 to 333.15) K. J. Chem. Eng. Data 2011, 57, 33–39. [Google Scholar] [CrossRef]
  27. Lashkarblooki, M.; Hezave, A.Z.; Al-Ajmi, A.M.; Ayatollahi, S. Viscosity prediction of ternary mixtures containing ILs using multi-layer perceptron artificial neural network. Fluid Phase Equilibria 2012, 326, 15–20. [Google Scholar] [CrossRef]
  28. Rooney, D.W.; Jacquemin, J.; Gardas, R.L. Thermophysical Properties of Ionic Liquids; Springer: Berlin/Heidelberg, Germany, 2009; Volume 290, pp. 185–212. [Google Scholar]
  29. Wu, W.; Li, W.; Han, B.; Zhang, Z.; Jiang, T.; Liu, Z. A green and effective method to synthesize ionic liquids: Supercritical CO2 route. Green Chem. 2005, 7, 701–704. [Google Scholar] [CrossRef]
  30. Lu, Y.; Liu, G.L.; Lee, L.P. High-Density Silver Nanoparticle Film with Temperature-Controllable Interparticle Spacing for a Tunable Surface Enhanced Raman Scattering Substrate. Nano Lett. 2005, 5, 5–9. [Google Scholar] [CrossRef]
  31. Yu, G.; Zhao, D.; Wen, L.; Yang, S.; Chen, X. Viscosity of ionic liquids: Database, observation, and quantitative structure-property relationship analysis. AIChE J. 2011, 58, 2885–2899. [Google Scholar] [CrossRef]
  32. Widegren, J.A.; Magee, J.W. Density, Viscosity, Speed of Sound, and Electrolytic Conductivity for the Ionic Liquid 1-Hexyl-3-methylimidazolium Bis(trifluoromethylsulfonyl)imide and Its Mixtures with Water†. J. Chem. Eng. Data 2007, 52, 2331–2338. [Google Scholar] [CrossRef]
  33. Burrell, G.L.; Burgar, I.M.; Separovic, F.; Dunlop, N.F. Preparation of protic ionic liquids with minimal water content and 15N NMR study of proton transfer. Phys. Chem. Chem. Phys. 2010, 12, 1571–1577. [Google Scholar] [CrossRef] [PubMed]
  34. Atashrouz, S.; Mirshekar, H.; Hemmati-Sarapardeh, A.; Moraveji, M.K.; Nasernejad, B. Implementation of soft computing approaches for prediction of physicochemical properties of ionic liquid mixtures. Korean J. Chem. Eng. 2016, 34, 425–439. [Google Scholar] [CrossRef]
  35. Hosseinzadeh, M.; Hemmati-Sarapardeh, A. Toward a predictive model for estimating viscosity of ternary mixtures containing ionic liquids. J. Mol. Liq. 2014, 200, 340–348. [Google Scholar] [CrossRef]
  36. Barycki, M.; Sosnowska, A.; Gajewicz, A.; Bobrowski, M.; Wileńska, D.; Skurski, P.; Giełdoń, A.; Czaplewski, C.; Uhl, S.; Laux, E.; et al. Temperature-dependent structure-property modeling of viscosity for ionic liquids. Fluid Phase Equilibria 2016, 427, 9–17. [Google Scholar] [CrossRef]
  37. Gardas, R.L.; Coutinho, J.A. A group contribution method for viscosity estimation of ionic liquids. Fluid Phase Equilibria 2008, 266, 195–201. [Google Scholar] [CrossRef]
  38. Gharagheizi, F.; Ilani-Kashkouli, P.; Mohammadi, A.H.; Ramjugernath, D.; Richon, D. Development of a group contribution method for determination of viscosity of ionic liquids at atmospheric pressure. Chem. Eng. Sci. 2012, 80, 326–333. [Google Scholar] [CrossRef]
  39. Lazzús, J.A.; Pulgar-Villarroel, G. A group contribution method to estimate the viscosity of ionic liquids at different temperatures. J. Mol. Liq. 2015, 209, 161–168. [Google Scholar] [CrossRef]
  40. Paduszyński, K.; Domańska, U. Viscosity of Ionic Liquids: An Extensive Database and a New Group Contribution Model Based on a Feed-Forward Artificial Neural Network. J. Chem. Inf. Model. 2014, 54, 1311–1324. [Google Scholar] [CrossRef]
  41. Zhao, Y.; Huang, Y.; Zhang, X.; Zhang, S. A quantitative prediction of the viscosity of ionic liquids using S σ-profile molecular descriptors. Phys. Chem. Chem. Phys. 2015, 17, 3761–3767. [Google Scholar] [CrossRef]
  42. Gaciño, F.M.; Paredes, X.; Comuñas, M.J.P.; Fernández, J. Effect of the pressure on the viscosities of ionic liquids: Experimental values for 1-ethyl-3-methylimidazolium ethylsulfate and two bis (trifluoromethyl-sulfonyl) imide salts. J. Chem. Thermodyn. 2012, 54, 302–309. [Google Scholar] [CrossRef]
  43. Gaciño, F.M.; Paredes, X.; Comuñas, M.J.P.; Fernández, J. Pressure dependence on the viscosities of 1-butyl-2, 3-dimethylimidazolium bis (trifluoromethylsulfonyl) imide and two tris (pentafluoroethyl) trifluorophosphate based ionic liquids: New measurements and modelling. J. Chem. Thermodyn. 2013, 62, 162–169. [Google Scholar] [CrossRef]
  44. Xu, Y.; Chen, B.; Qian, W.; Li, H. Properties of pure n-butylammonium nitrate ionic liquid and its binary mixtures of with alcohols at T=(293.15 to 313.15)K. J. Chem. Thermodyn. 2013, 58, 449–459. [Google Scholar] [CrossRef]
  45. Yu, Z.; Gao, H.; Wang, H.; Chen, L. Densities, Viscosities, and Refractive Properties of the Binary Mixtures of the Amino Acid Ionic Liquid [bmim][Ala] with Methanol or Benzylalcohol atT=(298.15 to 313.15) K. J. Chem. Eng. Data 2011, 56, 2877–2883. [Google Scholar] [CrossRef]
  46. Xu, Y.; Yao, J.; Wang, C.; Li, H. Density, Viscosity, and Refractive Index Properties for the Binary Mixtures of n-Butylammonium Acetate Ionic Liquid + Alkanols at Several Temperatures. J. Chem. Eng. Data 2012, 57, 298–308. [Google Scholar] [CrossRef]
  47. Koller, T.; Rausch, M.H.; Schulz, P.S.; Berger, M.; Wasserscheid, P.; Economou, I.G.; Leipertz, A.; Fröba, A.P. Viscosity, Interfacial Tension, Self-Diffusion Coefficient, Density, and Refractive Index of the Ionic Liquid 1-Ethyl-3-methylimidazolium Tetracyanoborate as a Function of Temperature at Atmospheric Pressure. J. Chem. Eng. Data 2012, 57, 828–835. [Google Scholar] [CrossRef]
  48. Domańska, U.; Zawadzki, M.; Lewandrowska, A. Effect of temperature and composition on the density, viscosity, surface tension, and thermodynamic properties of binary mixtures of N-octylisoquinolinium bis $(trifluoromethyl)sulfonyl$ imide with alcohols. J. Chem. Thermodyn. 2012, 48, 101–111. [Google Scholar] [CrossRef]
  49. Kanakubo, M.; Nanjo, H.; Nishida, T.; Takano, J. Density, viscosity, and electrical conductivity of N-methoxymethyl-N-methylpyrrolidinium bis(trifluoromethanesulfonyl)amide. Fluid Phase Equilibria 2011, 302, 10–13. [Google Scholar] [CrossRef]
  50. Fendt, S.; Padmanabhan, S.; Blanch, H.W.; Prausnitz, J.M. Viscosities of Acetate or Chloride-Based Ionic Liquids and Some of Their Mixtures with Water or Other Common Solvents. J. Chem. Eng. Data 2011, 56, 31–34. [Google Scholar] [CrossRef]
  51. Mokhtarani, B.; Mojtahedi, M.M.; Mortaheb, H.R.; Mafi, M.; Yazdani, F.; Sadeghian, F. Densities, Refractive Indices, and Viscosities of the Ionic Liquids 1-Methyl-3-octylimidazolium Tetrafluoroborate and 1-Methyl-3-butylimidazolium Perchlorate and Their Binary Mixtures with Ethanol at Several Temperatures. J. Chem. Eng. Data 2008, 53, 677–682. [Google Scholar] [CrossRef]
  52. Domańska, U.; Skiba, K.; Zawadzki, M.; Paduszyński, K.; Królikowski, M. Synthesis, physical, and thermodynamic properties of 1-alkyl-cyanopyridinium bis $(trifluoromethyl)sulfonyl$ imide ionic liquids. J. Chem. Thermodyn. 2013, 56, 153–161. [Google Scholar] [CrossRef]
  53. Da Rocha, M.M.A.; Ribeiro, F.M.S.; Ferreira, A.I.; Coutinho, J.A.; Santos, L.M.N.B.F. Thermophysical properties of [CN−1C1im][PF6] ionic liquids. J. Mol. Liq. 2013, 188, 196–202. [Google Scholar] [CrossRef]
  54. Diogo, J.C.F.; Caetano, F.J.P.; Fareleira, J.M.N.A.; Wakeham, W.A. Viscosity measurements of three ionic liquids using the vibrating wire technique. Fluid Phase Equilibria 2013, 353, 76–86. [Google Scholar] [CrossRef]
  55. Pires, J.; Timperman, L.; Jacquemin, J.; Balducci, A.; Anouti, M. Density, conductivity, viscosity, and excess properties of (pyrrolidinium nitrate-based Protic Ionic Liquid+propylene carbonate) binary mixture. J. Chem. Thermodyn. 2013, 59, 10–19. [Google Scholar] [CrossRef] [Green Version]
  56. Liu, Q.-S.; Li, P.-P.; Welz-Biermann, U.; Chen, J.; Liu, X.-X. Density, dynamic viscosity, and electrical conductivity of pyridinium-based hydrophobic ionic liquids. J. Chem. Thermodyn. 2013, 66, 88–94. [Google Scholar] [CrossRef]
  57. Liu, X.; Afzal, W.; Prausnitz, J.M. Unusual trend of viscosities and densities for four ionic liquids containing a tetraalkyl phosphonium cation and the anion bis(2,4,4-trimethylpentyl) phosphinate. J. Chem. Thermodyn. 2014, 70, 122–126. [Google Scholar] [CrossRef]
  58. Qian, W.; Xu, Y.; Zhu, H.; Yu, C. Properties of pure 1-methylimidazolium acetate ionic liquid and its binary mixtures with alcohols. J. Chem. Thermodyn. 2012, 49, 87–94. [Google Scholar] [CrossRef]
  59. Ochkedzan-Siodłak, W.; Dziubek, K.; Siodłak, D. Densities and viscosities of imidazolium and pyridinium chloroaluminate ionic liquids. J. Mol. Liq. 2013, 177, 85–93. [Google Scholar] [CrossRef]
  60. González, E.J.; González, B.; Calvar, N.; Dominguez, Á. Physical properties of binary mixtures of the ionic liquid 1-ethyl-3-methylimidazolium ethyl sulfate with several alcohols at T=(298.15, 313.15, and 328.15) K and atmospheric pressure. J. Chem. Eng. Data. 2007, 52, 1641–1648. [Google Scholar] [CrossRef]
  61. Yan, F.; He, W.; Jia, Q.; Wang, Q.; Xia, S.; Ma, P. Prediction of ionic liquids viscosity at variable temperatures and pressures. Chem. Eng. Sci. 2018, 184, 134–140. [Google Scholar] [CrossRef]
  62. Bird, R.B.; Stewart, W.E.; Lightfoot, E.N. Transport Phenomena; John Wiley & Sons: New York, NY, USA, 1960; p. 413. [Google Scholar]
  63. Kincaid, J.F.; Eyring, H.; Stearn, A.E. The Theory of Absolute Reaction Rates and its Application to Viscosity and Diffusion in the Liquid State. Chem. Rev. 1941, 28, 301–365. [Google Scholar] [CrossRef]
  64. Glasstone, S.; Laidler, K.J.; Eyring, H. The Theory of Rate Processes; MacGraw-Hill Book Co: New York, NY, USA, 1941; p. 477. [Google Scholar]
  65. Eyring, H. Viscosity, Plasticity, and Diffusion as Examples of Absolute Reaction Rates. J. Chem. Phys. 1936, 4, 283–291. [Google Scholar] [CrossRef]
  66. Gill, P.E.; Murray, W.; Wright, M.H. Practical Optimization; Academic Press: London, UK, 1981. [Google Scholar]
  67. Ameli, F.; Hemmati-Sarapardeh, A.; Dabir, B.; Mohammadi, A.H. Determination of asphaltene precipitation conditions during natural depletion of oil reservoirs: A robust compositional approach. Fluid Phase Equilibria 2016, 412, 235–248. [Google Scholar] [CrossRef]
  68. Wilde, D.J.; Beightler, C.S. Foundations of Optimization; Prentice Hall Inc.: Upper Saddle River, NJ, USA, 1967. [Google Scholar]
  69. Sharma, R.; Glemmestad, B. On Generalized Reduced Gradient method with multi-start and self-optimizing control structure for gas lift allocation optimization. J. Process. Control. 2013, 23, 1129–1140. [Google Scholar] [CrossRef]
  70. Morgan, J.N.; Sonquist, J.A. Some results from a non-symmetrical branching process that looks for interaction effects. Young 1963, 8, 5. [Google Scholar]
  71. Loh, W.-Y. Fifty Years of Classification and Regression Trees. Int. Stat. Rev. 2014, 82, 329–348. [Google Scholar] [CrossRef] [Green Version]
  72. Breiman, L.; Friedman, J.H.; Olshen, R.A.; Stone, C.J. Classification and Regression Trees; The Wadsworth & Brooks/Cole Mathematics Series; Springer: Berlin/Heidelberg, Germany, 1984. [Google Scholar]
  73. Safavian, S.R.; Landgrebe, D. A survey of decision tree classifier methodology. IEEE Trans. Syst. Man Cybern. 1991, 21, 660–674. [Google Scholar] [CrossRef] [Green Version]
  74. Song, Y.-Y.; Lu, Y. Decision tree methods: Applications for classification and prediction. Shanghai Arch. Psychiatry 2015, 27, 130–135. [Google Scholar]
  75. Patel, N.; Upadhyay, S. Study of Various Decision Tree Pruning Methods with their Empirical Comparison in WEKA. Int. J. Comput. Appl. 2012, 60, 20–25. [Google Scholar] [CrossRef]
  76. Tatar, A.; Shokrollahi, A.; Mesbah, M.; Rashid, S.; Arabloo, M.; Bahadori, A. Implementing Radial Basis Function Networks for modeling CO2-reservoir oil minimum miscibility pressure. J. Nat. Gas. Sci. Eng. 2013, 15, 82–92. [Google Scholar] [CrossRef]
  77. Tahmasebi, P.; Hezarkhani, A. A fast and independent architecture of artificial neural network for permeability prediction. J. Pet. Sci. Eng. 2012, 86, 118–126. [Google Scholar] [CrossRef]
  78. Atashrouz, S.; Pazuki, G.; Alimoradi, Y. Estimation of the viscosity of nine nanofluids using a hybrid GMDH-type neural network system. Fluid Phase Equilibria 2014, 372, 43–48. [Google Scholar] [CrossRef]
  79. Lashkarbolooki, M.; Hezave, A.Z.; Ayatollahi, S. Artificial neural network as an applicable tool to predict the binary heat capacity of mixtures containing ionic liquids. Fluid Phase Equilibria 2012, 324, 102–107. [Google Scholar] [CrossRef]
  80. Hemmati-Sarapardeh, A.; Ghazanfari, M.-H.; Ayatollahi, S.; Masihi, M. Accurate determination of the CO2-crude oil minimum miscibility pressure of pure and impure CO2streams: A robust modelling approach. Can. J. Chem. Eng. 2016, 94, 253–261. [Google Scholar] [CrossRef]
  81. Hagan, M.; Menhaj, M. Training feedforward networks with the Marquardt algorithm. IEEE Trans. Neural Netw. 1994, 5, 989–993. [Google Scholar] [CrossRef]
  82. Foresee, F.D.; Hagan, M. Gauss-Newton approximation to Bayesian learning. In Proceedings of the International Conference on Neural Networks, Houston, TX, USA, 9–12 June 1997; pp. 1930–1935. [Google Scholar]
  83. Mackay, D.J.C. Bayesian Interpolation. Neutral Comput. 1992, 447, 415–447. [Google Scholar] [CrossRef]
  84. Amar, M.N.; Noureddine, Z. An efficient methodology for multi-objective optimization of water alternating CO2 EOR process. J. Taiwan Inst. Chem. Eng. 2019, 99, 154–165. [Google Scholar] [CrossRef]
  85. Suykens, J.A.K.; Vandewalle, J. Least Squares Support Vector Machine Classifiers. Neural Process. Lett. 1999, 9, 293–300. [Google Scholar] [CrossRef]
  86. Amar, M.N.; Zeraibi, N. Application of hybrid support vector regression artificial bee colony for prediction of MMP in CO2-EOR process. Petroleum 2018. [Google Scholar] [CrossRef]
  87. Qin, L.-T.; Liu, S.-S.; Liu, H.-L.; Zhang, Y.-H. Support vector regression and least squares support vector regression for hormetic dose–response curves fitting. Chemosphere 2010, 78, 327–334. [Google Scholar] [CrossRef]
  88. Eslamimanesh, A.; Gharagheizi, F.; Mohammadi, A.H.; Richon, D. Phase Equilibrium Modeling of Structure H Clathrate Hydrates of Methane + Water “Insoluble” Hydrocarbon Promoter Using QSPR Molecular Approach. J. Chem. Eng. Data 2011, 56, 3775–3793. [Google Scholar] [CrossRef]
  89. Gharagheizi, F.; Eslamimanesh, A.; Farjood, F.; Mohammadi, A.H.; Richon, D. Solubility Parameters of Nonelectrolyte Organic Compounds: Determination Using Quantitative Structure–Property Relationship Strategy. Ind. Eng. Chem. Res. 2011, 50, 11382–11395. [Google Scholar] [CrossRef]
  90. Eslamimanesh, A.; Gharagheizi, F.; Illbeigi, M.; Mohammadi, A.H.; Fazlali, A.; Richon, D. Phase equilibrium modeling of clathrate hydrates of methane, carbon dioxide, nitrogen, and hydrogen+water soluble organic promoters using Support Vector Machine algorithm. Fluid Phase Equilibria 2012, 316, 34–45. [Google Scholar] [CrossRef]
  91. Chen, C.-H.; Lin, Z.-S. A committee machine with empirical formulas for permeability prediction. Comput. Geosci. 2006, 32, 485–496. [Google Scholar] [CrossRef]
  92. Naftaly, U.; Intrator, N.; Horn, D. Optimal ensemble averaging of neural networks. Netw. Comput. Neural Syst. 1997, 8, 283–296. [Google Scholar] [CrossRef] [Green Version]
  93. Haykin, S.; Network, N. A comprehensive foundation. Neural Netw. 2004, 2, 41. [Google Scholar]
  94. Hashem, S.; Schmeiser, B. Approximating a Function and Its Derivatives Using MSE-Optimal Linear Combinations of Trained Feedforward Neural Networks; Purdue University, Department of Statistics: West Lafayette, IN, USA, 1993. [Google Scholar]
  95. Yang, X.-S. A New Metaheuristic Bat-Inspired Algorithm. In Recent Advances in Computational Optimization; Springer: Berlin/Heidelberg, Germany, 2010; pp. 65–74. [Google Scholar]
  96. Yang, X.-S.; Deb, S. Eagle Strategy Using Lévy Walk and Firefly Algorithms for Stochastic Optimization. In Recent Advances in Computational Optimization; Springer: Berlin/Heidelberg, Germany, 2010; pp. 101–111. [Google Scholar]
  97. Yang, X.S.; He, X.-S. Bat algorithm: Literature review and applications. Int. J. Bio-Inspired Comput. 2013, 5, 141. [Google Scholar] [CrossRef] [Green Version]
  98. Shokrollahi, A.; Arabloo, M.; Gharagheizi, F.; Mohammadi, A.H. Intelligent model for prediction of CO2—Reservoir oil minimum miscibility pressure. Fuel 2013, 112, 375–384. [Google Scholar] [CrossRef]
  99. Vogel, H. The law of relationbetween the viscosity of liquids and the temperature. Phys. Z 1921, 22, 645–646. [Google Scholar]
  100. Garcia-Garabal, S.; Vila, J.; Rilo, E.; Dominguez-Pérez, M.; Segade, L.; Tojo, E.; Verdia, P.; Varela, L.M.; Cabeza, O. Transport properties for 1-ethyl-3-methylimidazolium n-alkyl sulfates: Possible evidence of grotthuss mechanism. Electrochim. Acta 2017, 231, 94–102. [Google Scholar] [CrossRef]
  101. Tammann, G.; Hesse, W. Die Abhängigkeit der Viscosität von der Temperatur bie unterkühlten Flüssigkeiten. Zeitschrift für Anorganische und Allgemeine Chemie 1926, 156, 245–257. [Google Scholar] [CrossRef]
  102. Fulcher, G.S. Analysis of Recent Measurements of the Viscosity of Glasses. J. Am. Ceram. Soc. 1925, 8, 339–355. [Google Scholar] [CrossRef]
  103. Taylor, K.E. Summarizing multiple aspects of model performance in a single diagram. J. Geophys. Res. Atmos. 2001, 106, 7183–7192. [Google Scholar] [CrossRef]
  104. Shateri, M.; Ghorbani, S.; Hemmati-Sarapardeh, A.; Mohammadi, A.H. Application of Wilcoxon generalized radial basis function network for prediction of natural gas compressibility factor. J. Taiwan Inst. Chem. Eng. 2015, 50, 131–141. [Google Scholar] [CrossRef]
  105. Hemmati-Sarapardeh, A.; Aminshahidy, B.; Pajouhandeh, A.; Yousefi, S.H.; Hosseini-Kaldozakh, S.A. A soft computing approach for the determination of crude oil viscosity: Light and intermediate crude oil systems. J. Taiwan Inst. Chem. Eng. 2016, 59, 1–10. [Google Scholar] [CrossRef]
  106. Atashrouz, S.; Mirshekar, H.; Mohaddespour, A. A robust modeling approach to predict the surface tension of ionic liquids. J. Mol. Liq. 2017, 236, 344–357. [Google Scholar] [CrossRef]
  107. Mokarizadeh, H.; Atashrouz, S.; Mirshekar, H.; Hemmati-Sarapardeh, A.; Mohaddespour, A. Comparison of LSSVM model results with artificial neural network model for determination of the solubility of SO2 in ionic liquids. J. Mol. Liq. 2020, 304, 112771. [Google Scholar] [CrossRef]
  108. Atashrouz, S.; Mirshekar, H.; Hemmati-Sarapardeh, A. A soft-computing technique for prediction of water activity in PEG solutions. Colloid Polym. Sci. 2017, 295, 421–432. [Google Scholar] [CrossRef]
  109. Atashrouz, S.; Hemmati-Sarapardeh, A.; Mirshekar, H.; Nasernejad, B.; Moraveji, M.K. On the evaluation of thermal conductivity of ionic liquids: Modeling and data assessment. J. Mol. Liq. 2016, 224, 648–656. [Google Scholar] [CrossRef]
Figure 1. Chemical structure of ionic liquids (ILs) in this study.
Figure 1. Chemical structure of ionic liquids (ILs) in this study.
Molecules 26 00156 g001aMolecules 26 00156 g001bMolecules 26 00156 g001c
Figure 2. Illustration of an escape process in the flow of a liquid. Molecule 1 must pass through a “bottleneck” to reach the vacant site.
Figure 2. Illustration of an escape process in the flow of a liquid. Molecule 1 must pass through a “bottleneck” to reach the vacant site.
Molecules 26 00156 g002
Figure 3. Committee Machine Intelligent System (CMIS) model offered in this study.
Figure 3. Committee Machine Intelligent System (CMIS) model offered in this study.
Molecules 26 00156 g003
Figure 4. Cross plot of the proposed Eyring’s theory for viscosity of ILs.
Figure 4. Cross plot of the proposed Eyring’s theory for viscosity of ILs.
Molecules 26 00156 g004
Figure 5. Logarithm of experimental viscosity data against predicted value based on Model (II) for CMIS model in this study.
Figure 5. Logarithm of experimental viscosity data against predicted value based on Model (II) for CMIS model in this study.
Molecules 26 00156 g005
Figure 6. Relative deviation of prediction of CMIS model versus logarithm of experimental data based on Model (II).
Figure 6. Relative deviation of prediction of CMIS model versus logarithm of experimental data based on Model (II).
Molecules 26 00156 g006
Figure 7. Comparison between absolute average relative deviation (AARD, %) of different models for estimation of viscosity of ILs based on Model (II).
Figure 7. Comparison between absolute average relative deviation (AARD, %) of different models for estimation of viscosity of ILs based on Model (II).
Molecules 26 00156 g007
Figure 8. Cumulative frequency of absolute relative deviation for different models based on Model (II).
Figure 8. Cumulative frequency of absolute relative deviation for different models based on Model (II).
Molecules 26 00156 g008
Figure 9. Taylor diagram for the offered methods based on Model (II).
Figure 9. Taylor diagram for the offered methods based on Model (II).
Molecules 26 00156 g009
Figure 10. Zonation of the plotted points in the Taylor diagram based on Model (II).
Figure 10. Zonation of the plotted points in the Taylor diagram based on Model (II).
Molecules 26 00156 g010
Figure 11. Cross plot of the proposed developed model based on Model (III) for CMIS approach.
Figure 11. Cross plot of the proposed developed model based on Model (III) for CMIS approach.
Molecules 26 00156 g011
Figure 12. Relative deviation distribution for developed CMIS model in this study for estimation of the viscosity of ILs based on Model (III).
Figure 12. Relative deviation distribution for developed CMIS model in this study for estimation of the viscosity of ILs based on Model (III).
Molecules 26 00156 g012
Figure 13. Cumulative frequency of absolute relative deviation in different models based on Model (III).
Figure 13. Cumulative frequency of absolute relative deviation in different models based on Model (III).
Molecules 26 00156 g013
Figure 14. Comparison among AARD (%) of different models for prediction of viscosity of ILs based on Model (III).
Figure 14. Comparison among AARD (%) of different models for prediction of viscosity of ILs based on Model (III).
Molecules 26 00156 g014
Figure 15. Effect of temperature on 1-ethyl-3-methylimidazolium ethylsulfate. (a): Model (II); (b): Model (III).
Figure 15. Effect of temperature on 1-ethyl-3-methylimidazolium ethylsulfate. (a): Model (II); (b): Model (III).
Molecules 26 00156 g015
Figure 16. The Williams plot of the whole dataset for the CMIS model based on Model (II).
Figure 16. The Williams plot of the whole dataset for the CMIS model based on Model (II).
Molecules 26 00156 g016
Figure 17. Relevancy factor of the CMIS model inputs based on Model (II).
Figure 17. Relevancy factor of the CMIS model inputs based on Model (II).
Molecules 26 00156 g017
Table 1. The selected ionic liquids in the present study.
Table 1. The selected ionic liquids in the present study.
Component of ionic liquidAbbreviationnT
(K)
P (MPa)
1-butyl-3-methylimidazolium hexafluorophosphate[C4mim] [PF6]238273.15–413.150.1–249.3
1-octyl-3-methylimidazolium hexafluorophosphate[C8mim] [PF6]132273.15–363.150.1–175.9
1-hexyl-3-methylimidazolium hexafluorophosphate[HMIM] [PF6]179273.15–238.50.1–238.5
1-octyl-3-methylimidazolium tetrafluoroborate[C8mim] [BF4]141273.15–363.150.1–224.2
1-hexyl-3-methylimidazolium tetrafluoroborate[C6mim] [BF4]183283.15–368.150.1–121.8
1-butyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl] imide[C4mim] [Tf2N]344273.15–5730.1–298.9
1-ethyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl] imide[C2mim] [Tf2N]225263.15–388.190.1–125.5
1-octyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl] imide[C8mim] [Tf2N]25278–363.150.1
1-hexyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl] imide[C6mim] [Tf2N]236258.15–433.150.1–124
1-butyl-3-methylimidazolium trifluoromethanesulfonate[C4mim] [CF3SO3]25283.15–363.150.1
1-ethyl-3-methylimidazolium ethylsulfate[C2mim] [EtSO4]137253.15–388.190.1–75
1-hexylpyridinium bis[(trifluoromethyl)sulfonyl] imide[HPy] [Tf2N]8283–3430–1
1-butylpyridinium bis[(trifluoromethyl)sulfonyl] imide[BPy] [Tf2N]9283.15–353.150.1
1-butyl-1-methylpyrrolidinium bis[(trifluoromethyl) sulfonyl]imide[C4MPyr] [Tf2N]148273.15–5730.1–102.9
1-ethylpyridinium ethylsulfate[EPy] [ESO4]8283–3430.1
trimethylhexylammonium bis[(trifluoromethyl)sulfonyl]imide[N1116] [Tf2N]1293.150.1
Trimethylbuthlammonium bis[(trifluoromethyl)sulfonyl]imide[N1114] [Tf2N]17293.15–388.510.1
1-butyl-3-methylimidazolium tris(pentafluoroethyl) trifluorophosphate[C4mim] [FAP]1293.150.1
1,2-dimethylimidasolium bis[(trifluoromethyl)sulfonyl] imide[DMIM] [Tf2N]1298.150.1
trihexyl(tetradecyl)phosphonium tris(pentafluoroethyl) trifluorophosphate[P6,6,6,14] [FAP]181268.15–373.150.1
1-butyl-1-methylpyrrolidinium tris(pentafluoroethyl) trifluorophosphate[C4mpyrr] [FAP]67283.15–373.150.1–150
1-butyl-1-methylpyrrolidinium trifluoromethanesulfonate[BMPyr] [TfO]67293.15–373.150.1–150
1-ethyl-3-methylimidazolium hydrogensulfate[C2mim] [HSO4]22268.15–373.150.1
trimethylpropylammonium bis[(trifluoromethyl)sulfonyl] imide[N1113] [Tf2N]6293–3180.1
1-heptyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl] imide[C7mim] [Tf2N]12930.1
1-undecyl-3-methylimidazolium tetrafluoroborate[C11MIM] [BF4]8293–3630.1
1-butyl-3-methylimidazolium iodid[C4mim] [I]35289.15–388.150.1
1-butyl-3-methylimidazolium nitrate[C4mim] [NO3]27283.15–363.150.1
1-dodecyl-3-methylimidazolium hexafluorophosphate[C12MIM] [PF6]4333.15–363.150.1
1-octyl-3-methylimidazolium nitrate[C8mim] [NO3]16283.15–363.150.1
1-hexyl-3-methylimidazolium nitrate[C6mim] [NO3]14283.15–363.150.1
1-butylpyridinum tetrafluoroborate[BPy] [BF4]70278.15–338.150.1–65.9
1-hexylpyridinium bis[(trifluoromethyl)sulfonyl] imide[C6Py] [Tf2N]9298.15–398.150.1
1-heptyl-3-methylimidazolium hexafluorophosphate[C7mim] [PF6]13293.15–263.150.1
1-ethyl-3-methylimidazolium diethylphosphate[C2mim] [DEP]17292.15–373.150.1
1-pentyl-3-methylimidazolium hexafluorophosphate[C5mim] [PF6]13293.15–263.150.1
1-nonyl-3-methylimidazolium hexafluorophosphate[C9mim] [PF6]12303.15–363.150.1
1,2-dimethyl-3-propylimidazolium tetrafluoroborate[M1,2P3im] [BF4]8289.15–343.150.1
1-butyl-4-methylpyridinium tetrafluoroborate[mbpy] [BF4]48283.15–333.150.1–65
1,3-dimethylimidazolium dimethylphosphate[C1mim] [DPO4]7293.15–323.150.1
1,2-dimethyl-3-propylimidazolium bis[(trifluoromethyl)sulfonyl] imide[M1,2P3im] [Tf2N]16290–3650.1
1-ethyl-3-methylimidazolium methylsulfate[C2mim] [MSO4]27283.15–373.150.1
1-ethyl-3-methylimidazolium methanesulfonate[C2mim] [mesy]45278.15–363.150.1
1-butyl-3-methylimidazolium perchlorate[C4mim] [CLO4]15283.15–383.150.1
1-butyl-2,3-dimethylimidazolium tetrafluoroborate[BDmim] [BF4]7298.15–353.150.1
Table 2. Statistical details about gathered databank in the present work.
Table 2. Statistical details about gathered databank in the present work.
ParameterSymbolUnitMinMaxMean
TemperatureTK253.15573325.63
PressurePMPa0.06298.9024.45
Molecular WeightMwg/mole201.22515.13346.65
Critical TemperatureTcK520.061534.631005.87
Critical PressurePcbar2.6357.6022.29
Critical VolumeVccm3/mol550.652573.60992.83
Acentric factorΩ-0.211.100.59
Boiling TemperatureTbK410.771130.30723.93
Experimental viscosityη expMPa.s1.139667.62191.91
Table 3. Statistical parameters of the developed models for prediction of the viscosity of ILs based on Model (II).
Table 3. Statistical parameters of the developed models for prediction of the viscosity of ILs based on Model (II).
DTLSSVM–BATMLP–LMAMLP–BRCMIS
Training set
ARD%−3.273−0.402−0.326−0.1700.011
AARD%13.3667.8994.6473.5533.256
RMSE152.76711.5348.4658.5179.533
SD0.2190.1620.1020.0640.084
R20.8800.9990.9990.9990.999
Number of Data point22502250225022502250
Test set
ARD%−4.139−0.190−0.171−1.1400.258
AARD%17.4548.1514.9295.0043.117
RMSE223.79325.35622.36820.0449.035
SD0.2440.1400.0890.2320.050
R20.7510.9940.9970.9970.999
Number of Data point563563563563563
Total
ARD%−3.447−0.345−0.260−0.348−0.207
AARD%14.1847.9414.7073.8413.293
RMSE169.38315.33112.54811.77011.812
SD0.2250.1580.0990.1180.083
R20.8530.9980.9990.9990.999
Number of Data point28132813281328132813
Table 4. Statistical parameters of the developed models for prediction of the viscosity of ILs based on Model (III).
Table 4. Statistical parameters of the developed models for prediction of the viscosity of ILs based on Model (III).
DTLSSVM−BATMLP−LMAMLP−BRCMIS
Training set
ARD%−1.1080.3570.020−0.219−0.161
AARD%13.5895.5524.5223.7683.422
RMSE15.4449.0279.4417.9238.894
SD0.2920.1290.1140.0760.073
R20.9980.9990.9990.9990.999
Number of Data point22502250225022502250
Test set
ARD%5.384−0.991−0.464−0.977−0.412
AARD%24.3456.6044.6244.9493.454
RMSE24.30928.98422.67310.0566.844
SD1.9600.1610.0920.1500.064
R20.9950.9910.9980.9990.999
Number of Data point563563563563563
Total
ARD%0.1900.087−0.076−0.371−0.172
AARD%15.7425.7634.4534.0043.426
RMSE17.58015.27513.1978.3939.505
SD0.9160.1360.1100.0950.073
R20.9980.9980.9990.9990.999
Number of Data point28132813281328132813
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mousavi, S.P.; Atashrouz, S.; Nait Amar, M.; Hemmati-Sarapardeh, A.; Mohaddespour, A.; Mosavi, A. Viscosity of Ionic Liquids: Application of the Eyring’s Theory and a Committee Machine Intelligent System. Molecules 2021, 26, 156. https://doi.org/10.3390/molecules26010156

AMA Style

Mousavi SP, Atashrouz S, Nait Amar M, Hemmati-Sarapardeh A, Mohaddespour A, Mosavi A. Viscosity of Ionic Liquids: Application of the Eyring’s Theory and a Committee Machine Intelligent System. Molecules. 2021; 26(1):156. https://doi.org/10.3390/molecules26010156

Chicago/Turabian Style

Mousavi, Seyed Pezhman, Saeid Atashrouz, Menad Nait Amar, Abdolhossein Hemmati-Sarapardeh, Ahmad Mohaddespour, and Amir Mosavi. 2021. "Viscosity of Ionic Liquids: Application of the Eyring’s Theory and a Committee Machine Intelligent System" Molecules 26, no. 1: 156. https://doi.org/10.3390/molecules26010156

APA Style

Mousavi, S. P., Atashrouz, S., Nait Amar, M., Hemmati-Sarapardeh, A., Mohaddespour, A., & Mosavi, A. (2021). Viscosity of Ionic Liquids: Application of the Eyring’s Theory and a Committee Machine Intelligent System. Molecules, 26(1), 156. https://doi.org/10.3390/molecules26010156

Article Metrics

Back to TopTop