Printer Friendly

Selection and Configuration of Sorption Isotherm Models in Soils Using Artificial Bees Guided by the Particle Swarm.

1. Introduction

The fate and transport of heavy metals in the soils are a serious concern due to their potential impact on the environment. These reactive substances in soils undergo complex interactions based on the soil properties, type of clay minerals, and pore-fluid parameters. The problem of identifying the fate and transport of these substances in soils requires the understanding of various sorption mechanisms of species in the soil environment [1]. The sorption reactions in soils are, therefore, important in contaminant transport modeling in the subsurface and containment facilities such as landfills. Contaminant sorption process of different solutes in the soils is generally described using isotherms. Sorption isotherms describe the concentration of the sorbed solute to the sorbent concentration at equilibrium. Several linear and nonlinear isotherms are commonly used in various fields of research to understand the sorption mechanism between sorbent (soil solids) and solute [2-4]. Precise selection and configuration of sorption isotherm for a given sorbent-solute combination are essential in the contaminant risk assessment. The selection and configuration depend on chosen isotherm model, experimental representation of the environmental conditions, and parameter estimation method [5]. The sorption model parameters and selection of appropriate isotherm models are accomplished by fitting isotherm models to the experimental data using regression [6]. Nonlinear least-square techniques are widely used for this purpose [4,5,7, 8]. However, it will be shown in this work that classical gradient-based techniques suffer from several limitations.

Nature-inspired algorithms such as algorithms based on swarm intelligence (SI) are global search techniques and are widely used for many optimization problems in engineering. Particle swarm optimization (PSO) [9] and artificial bee colony (ABC) [10] are commonly used SI techniques. These algorithms overcome several limitations of conventional, deterministic techniques. Thus, these global search techniques are widely applied in geotechnical engineering for several optimization problems [11-17]. However, these techniques often prematurely converge to local optimum solution and may not generate the same solution in different runs due to stochastic nature of the algorithm [18]. A hybrid method based on particle swarm optimization and gradient-based algorithm is often used for configuration of sorption isotherms [19]. However, the limitations associated with the slow convergence rate in early stage of the search and lack of diversity in later stage of the search may remain in the basic particle swarm optimization [20, 21]. Therefore, hybrid methods may not balance the limitations that arose from each other. Such hybrid methods are also computationally expensive as they require both gradient information and more number of objective function evaluations.

It will be shown that the performance of conventional optimization techniques and PSO algorithm is greatly hampered by the presence of several local optima on the functional terrain. An inverse model based on artificial bee-colony (ABC) optimization method is proposed in this work for selection and configuration of appropriate sorption isotherms using the experimental sorption data. The proposed solver is further improved by introducing several modifications to the ABC algorithm for better performance. The proposed technique is highly robust; it accurately estimates the sorption model parameters and appropriate isotherms to the experimental data.

2. Theory

Sorption of different chemical constituents in soils is routinely performed in the laboratory under controlled temperature, pressure, and pH. The measured sorption data is generally described by either equilibrium or a kinetic retention process [2]. The present study is limited to only modeling of equilibrium retention process. The measured, equilibrium sorption data is represented as sorbed concentration, [S.sub.i], against the equilibrium concentration in the aqueous phase, [C.sub.i]. Sorption isotherms describe the mathematical dependency of sorbed concentration on the equilibrium concentration. The determination of the sorption model parameters from this sorption data is essential for the design of containment facilities. Further, configuration of the observed data using available models is important in the modeling of contaminant transport and sorption processes for containment applications. Sorption isotherms are classified into S-curve, L-curve, H-curve, and C-curve based on the initial slope of the curve [22-24]. A thorough description of various sorption isotherm models is given elsewhere [24-27]. Some of the commonly used isotherms in geoenvironmental engineering applications are given below.

2.1. Freundlich Isotherm. The Freundlich isotherm is widely used for describing nonideal and reversible adsorption [27, 28]. This isotherm is applied in heterogeneous systems such as soils. The Freundlich isotherm can be expressed as

S = [K.sub.f][C.sup.a], (1)

where S is the amount of solute retained by the soil ([micro]g/g), C is the solute concentration in solution ([micro]/mL), [K.sub.j] is the distribution coefficient (mL/g), and the parameter a is dimensionless constant. The distribution coefficient describes the partitioning of solute species between solid and liquid phases [2]. The linear sorption equation is a particular case of Freundlich expression and can be derived by substituting a = 1.

2.2. Langmuir Isotherm. The Langmuir model is another popularly used sorption model, originally developed by Langmuir [29] to describe the adsorption of gases by solids. This isotherm model is widely used in geoenvironmental applications. For example, the model is recently applied for understanding the phosphorus contamination in soils [30], ammonium removal using zeolites [31], and phenol adsorption in natural soils [32]. The advantage of this model is that a maximum sorptive capacity is used in the model, which is regarded as a measure of the available sorption sites on the solid phase. The Langmuir isotherm is defined as

S/[S.sub.max] = bC/1 + bC, (2)

where b and [S.sub.max] are the fitting parameters. The parameter b (mL/g) is a measure of the bond strength at which the sorbed solutes are held on the soil surface and [S.sub.max] ([micro]g/g) is the maximum sorption capacity. This isotherm is extensively used to describe the sorption of solutes by soils [5]. Often various modifications to (1) and (2) are used [23] as some sorption data do not follow the standard isotherms. These modified forms are presented in the following subsections.

2.3. Freundlich-Langmuir Model. The Freundlich-Langmuir (FL) model is a power function based on the assumption of continuously distributed affinity coefficients [4, 23]. This model can be described as

S/[S.sub.max] = [bC.sup.a]/1 + [bC.sup.a], (3)

where a is the fitting parameter (0 < a [less than or equal to] 1). The Langmuir isotherm is a special case of this model and can be derived by substituting a = 1. The FL model is widely used for simulating the sorption behaviour of heavy metals on clay minerals [33, 34].

2.4. Two-Site Langmuir Model. The two-site Langmuir model is based on the assumption that sorption occurs on two types of surface sites, each with different bonding energies. One site contains a high bonding energy and reacts rapidly, while the other contains a lower bonding energy and reacts more slowly [35]. The modified Langmuir model with two sorbing components is given by

[mathematical expression not reproducible], (4)

where[mathematical expression not reproducible] are the maximum quantities of solute that can be sorbed by the two sites; and [b.sub.1] and [b.sub.2] are constants related to the bonding energies of the components. Several researchers obtained excellent fit of phosphate sorption data with the two-site Langmuir model [36] and metal sorption, that is, Lanthanide on smectite clays [37, 38]. Several other isotherm models such as Farley-Dzombak-Morel [39] and Toth [40] are often used, but the modeling is out of the scope of this work.

3. Estimation of Model Parameters

Accurate estimation of model parameters of (1) through (4) is required for the prediction of fate and transport of contaminants in the environment. These model parameters are frequently estimated by nonlinear, least-square approaches by minimizing the objective function. The objective function is an error measure which is used to determine the accuracy of a given model to describe the measured data. Several error measures, such as coefficient of determination ([R.sup.2]), sum of squared deviations (SSD), mean error (MSE), and root mean square error (RMSE), are commonly used in the inverse analysis [27]. The RMSE is often used to verify the appropriateness of the model to the observed data and is given by

[mathematical expression not reproducible], (5)

where [mathematical expression not reproducible] are the observed and fitted sorption values, respectively, and N is the number of observed data.

As the aforementioned theoretical isotherms are nonlinear, optimization techniques are used to estimate the model parameters of a given model isotherm and to determine the appropriate isotherm for a given data. Nonlinear regression techniques such as Levenberg-Marquardt method are widely used in the literature for parameterization [34,41,42]. A brief description of gradient-based technique is given here.

3.1. Gradient-Based Approach. Gradient-based techniques are the conventional optimization tools for parameter estimation in many engineering applications. The principle behind such techniques is determining an optimum point on the search space where the derivatives of the objective function with respect to the model parameters ([lambda]) become zero. The unknown model parameter vector ([P.sub.i]), in these methods, is determined iteratively. The direction of descent for the gradient-based method is based on

[P.sub.i+i] = [P.sub.i] - [alpha][lambda], (6)

where [alpha] is a scalar which determines the step length in the direction of [lambda]. Several variants of gradient-based algorithms are used that may require the information of second derivative of the objective function (Hessian matrix) for a faster convergence rate.

It was in this work that these techniques suffered from several limitations. These conventional algorithms are highly sensitive to the fitness function (error measure) and on the initial guess. Of late, nature-inspired or swarm intelligence search techniques are used as an alternative to the conventional methods for solving several optimization problems in geoenvironmental engineering [11-13,15-18].

3.2. Swarm Intelligence. Swarm intelligence (SI) is a field of research which develops computational techniques by taking inspiration from the nature for solving the optimization problems. Genetic algorithms (GA), particle swarm optimization (PSO), and artificial bee colony (ABC) are some of the techniques developed by mimicking the nature for the determination of global optimum solutions to the inverse problems. This paper developed several solvers based on some of the SI techniques for predicting the model parameters of isotherm equations and for determining the appropriate isotherms for the measured data. The following SI techniques were used in developing the solvers.

3.2.1. Particle Swarm Optimization. PSO is a class of stochastic, derivative-free, population-based method. The working principle of PSO method is inspired by the social behaviour of animals such as flocking of birds and schooling of fishes [9]. In PSO, each particle (i.e., agent) of the population (i.e., swarm) is thought to be a collision-free bird in an n-dimensional search space, where n is the number of model parameters to be optimized. PSO algorithm starts with initializing predefined number of particles in the search space represented by the objective function. Each particle is updated with a position and velocity at any given iteration in the search space. The position and velocity are updated based on the flying experience (movement) of individual particles and the neighbors to achieve a better position in the search space. Further, the best location of an individual particle in the history (all the previous iterations) and best location experienced in comparison to all the particles in any given iteration are rememorized. Therefore, the movement of the particle is an aggregated acceleration towards its best previously visited location and towards the best individual of a topological neighborhood [43]. The velocity ([v.sub.ij]) and movement ([x.sub.ij]) of the ith particle in the jth dimension are updated using [44, 45]

[v.sub.ij] = [v.sub.ij] + [[phi].sub.1][[gamma].sub.1] {[p.sub.ij] - [x.sub.ij]) + [[phi].sub.2][[gamma].sub.2] {[p.sub.gj] - [x.sub.ij]), (7)

[x.sub.ij] =[x.sub.ij] + [v.sub.ij], (8)

where [p.sub.ij] is the history best position of each individual particle; [p.sub.gj] is the global-best position; [[phi].sub.1] and [[phi].sub.2] are the acceleration constants that accelerate movement of the particles towards individual and global-best positions, respectively; [[gamma].sub.1] and [[gamma].sub.2] are two independent random numbers that are uniformly distributed between 0 and 1 used to stochastically vary the relative pull between [p.sub.ij] and [p.sub.gj]. The randomness is introduced using the constants for providing unpredictable behaviour of the swarm which is useful for preventing from premature convergence. This iterative process continues swarm by swarm until the predefined number of iterations. Individual particles are drawn towards the success of each other during the search process which results in the formation of a cluster of particles at the global optimum region.

PSO is widely used for solving the optimization problems in many engineering disciplines due to the simplicity. Several experiment studies on the application of PSO algorithm on benchmark functions and optimization problems revealed that the algorithm explodes quickly as the first term in the velocity term in (7) increases incessantly [45]. The notable work on improving the update equation are introduction of the inertia factor to control the velocity of the particles [44] and constriction factor which constrains the magnitude of the velocity [45]. The modified form of velocity-update equation is [12, 44, 45]

[v.sub.ij] = [chi]([omega][v.sub.ij] + [[phi].sub.1][[gamma].sub.1]([p.sub.ij] - [x.sub.ij]) + [[phi].sub.2] [[gamma].sub.2]([p.sub.gj] - [x.sub.ij])), (9)

where [chi] is the constriction factor and w is the inertia weight to improve the performance of the PSO algorithm. The first term in (9) is the inertia term to restrict the particles to change the flight direction drastically. Considerable improvements to (9) are available in the literature for improving the performance of the PSO algorithm. The present work uses perturbed PSO algorithm for developing the solver.

3.2.2. Variant of PSO Algorithm. The classical PSO algorithm and many of the variants suffer from few limitations during the search process. The major drawback with these algorithms is that they converge to suboptimal solutions for several problems due to the application of the same movement update equation (10) for gbest particle [18]. Bharat et al. [18] proposed perturbed PSO (PPSO) algorithm which uses the following perturbation equation for the gbest particle:

x(i, j) [left arrow] x(i, j)[l + [[phi].sub.3] ([[gamma].sub.3] - [[gamma].sub.4])], (10)

where [[phi].sub.3] is a constant and [[gamma].sub.3], [[gamma].sub.4] [member of] U(0,1).

Equation (10) is used to randomly perturb the global-best particle in the suboptimal region to avoid convergence. The performance of the PPSO algorithm is further improved using catfish particles [46] into the swarm [14]. The particles are arranged in descending order of their fitness values after few generations (catIter). A known percentage of low-fit solutions are chosen as catfish particles and are repositioned on the D-dimensional search space based on opposition-based learning (OBL). The OBL technique is based on the utilization of opposition numbers of the current positions in the search space [47]. This technique is explained in the following subsections. The catfish particles improve the efficiency by exploring the new feasible regions in the search space. This optimization algorithm is referred to as PCPSO algorithm in this work. PCPSO algorithm is believed to be alleviating premature convergence of the swarm even with small population due to unique combination of perturbation equation for gbest particle and exploration capabilities of the catfish particles [14]. Two other solvers based on SI algorithms are also developed in this work.

3.2.3. Artificial Bee-Colony (ABC) Algorithm. ABC algorithm is a recent addition to the class of global optimization algorithms. ABC algorithm [10] is inspired by the foraging behaviour of bees. It is shown recently on over 50 benchmark functions (unimodal to multimodal, separable to nonseparable functions) that ABC algorithm outperforms other popular SI techniques such as genetic algorithm (GA), differential evolution (DE), and particle swarm optimization (PSO) [48] due to better trade-off between exploration and exploitation.

The algorithm uses several steps inspired by the collective behaviour of bees for the collection and processing of the nectar. The basic ABC algorithm uses three classes of bees, that is, employed bees, onlooker bees, and scout bees. Employed bees are responsible for locating the nectar sources and sharing this information with the onlookers on the dance area of the hive. Onlooker bees select nectar source based on the quality (objective function) with some probability. Therefore, the probability of higher number of onlookers at better quality of food source is higher. The task of employed and onlooker bees is to explore the search space. Scout bees are similar to the catfishes used to exploit new food sources in the search space by random search process. Further, the employed bees with abandoned food sources become the scout bees. Each food source (a set of model parameters) is a feasible solution of the problem and the nectar amount of a food source represents quality of the solution, that is, the objective function. The quality of the solution is represented by fitness value as given by

[mathematical expression not reproducible], (11)

where f is the objective function (5) and x is the model parameter vector of the problem. Half of the colony is assumed to be comprised of employed bees and the other half by the onlookers. Each food source is exploited by only one employed bee. Therefore, the number of employed or the onlooker bees is equal to the number of food sources. The following four steps are used in the algorithm.

(i) Initialization Phase. A predefined number of food sources are randomly generated in the n-dimensional parametric space where "n" is the dimension of the parametric space representing the number of parameters to be optimized. The food sources are generated using

[x.sub.i,j] = [x.sup.min.sub.i,j] + rand (0,1) ([x.sup.max.sub.i,j] - [x.sup.min.sub.i,j]), (12)

where i = 1, 2, ..., [N.sub.food] ([N.sub.food] is the number of food sources and equals half of the colony size); j = 1,2, ..., D (D is the dimension of the problem); [x.sup.min.sub.i,j] and [x.sup.max.sub.ij] are the lower and upper bounds of the jth parameter and ith food source, respectively. The fitness of all the food sources is evaluated using a suitable fitness norm. Further, counters which store the number of trails of solutions are reset to 0 in this phase.

(ii) Employed Bees Phase. Each employed bee produces a modification to the solution in its memory depending on the local (visual) information. The modification to the present food source is determined using the following equation:

[n.sub.i,j] = [x.sub.i,j] + [omega][phi] {[x.sub.i,j] - [x.sub.k,j], (13)

where [omega] represents the inertia factor (usually [omega] = 1 used in canonical ABC), k represents a randomly selected food source different from i, and [phi] is a uniformly distributed random number in the range [-1,1]. The fitness of the newly discovered solution is computed. A greedy selection process is applied between the new ([n.sub.i,j]) and the original ([x.sub.i,j]) solution. The better solution is stored in the memory. The trails count for this food will be reset to 0 if the solution is improved; otherwise, its value is increased by 1. Employed bees share the nectar (solution) information with the onlooker bees on the dance area after this phase.

(iii) Onlooker Bees Phase. The food sources are chosen based on a probability of its nectar amount (quality of the objective function) in this phase. A high number of onlooker bees are swarmed around the same food source when the fitness of the source is high. The probability is calculated according to

[mathematical expression not reproducible], (14)

where the fitness of each agent is calculated according to (11). After this stage, each onlooker bee finds a new food source in the neighborhood based on (13), similar to the employed bees. The new position is replaced with the old one in the memory when the quality of the new position is better than that of the previous one. Otherwise, the bee keeps the old position in memory and the trail count is raised by 1.

(iv) Scout Bees Phase. A food source is abandoned and a new one is selected when the food source does not improve the quality in a predetermined number of iterations, called the "limit." The new food source is selected randomly in the search space using (2).

Some modifications to the original ABC algorithm have been studied on benchmark functions. The modification of the basic ABC involves modifying the position-update equation [49] which was used recently by Zhu and Kwong [49] on six benchmark functions and found favorable results. However, the convergence towards gbest solution cannot be controlled in the Zhu and Kwong [49] work due to lack of acceleration factors. Apart from the modification to position-update equation, the frequency and magnitude of the perturbation are modified [50, 51]. The OBL was used by Gao and Liu [52] for the initialization phase on benchmark functions. Hybrid techniques combining crossover operations of GA with ABC is used recently by Karaboga and Kaya [53].

3.2.4. Modified ABC (MABC) Algorithm. The exploration-exploitation qualities of the classical ABC algorithm are improved in the present work. The following position-update equation is used for the employed bees in the modified algorithm:

[n.sub.i,j] = [x.sub.i,j] + [epsilon][phi]([x.sub.i,j] - [x.sub.k,j]) + r (0,1) [epsilon][C.sub.3]([x.sub.gbest,j] - [x.sub.i,j]), (15)

where r is the random number varying within (0,1), [epsilon] is a constant, and [C.sub.3] is the acceleration coefficient. The perturbation equation in (15) produces perturbed solutions around the present solution in randomly chosen direction of another food source. The information of fittest solutions is not required in this method as the food sources are randomly selected to generate the new solutions.

The movement of the employed bees using the proposed position-update equation is based on the resultant direction of best-fit food source (position) and the randomly chosen food source location. The inertia weight and acceleration coefficient improve the exploration process. The position-update equation for onlooker bees is based on (13) with an adaptive inertia factor. This helps in carrying out the exploitation process for better food sources around the present positions. Further, opposition-based concept is used to generate new food sources for the scout bees.

OBL is based on the utilization of opposition numbers of the current positions in the search space [47]. If [x.sub.ij] = ([x.sub.11], [x.sub.12], [x.sub.1d]) is a food source position in the n-dimensional space, and [x.sub.1j] [member of] [[LL.sub.i], [UL.sub.i]], i = 1, 2, ..., n, the new position based on the OBL is [x.sup.*.sub.ij] = ([x.sup.*.sub.11], [x.sup.*.sub.12], ..., [x.sup.*.sub.1n]) of [x.sub.1j] which is defined as

[x.sup.*.sub.ij] = [LL.sub.i] + [UL.sub.i] - [x.sub.ij], i=1, 2, ..., n, (16)

where LL and UL are the lower and upper limits of the ith dimension.

4. Inverse Analysis

The performance of all the SI techniques is based on the tuning parameters. The tuning parameters of these algorithms are very sensitive to the nature of the parametric space [14,18, 54]. The tuning parameters of all the four algorithms and the measured data are defined in the following subsection.

4.1. Experimental Sorption Data. Observed sorption data of metal substances, that is, Cu, Zn, Ca, and Pb, on various soils were used from the literature for selection and configuration of the described isotherms. Natural soils having various percentages of clay content, different clay minerals, and pH were chosen. The physical properties of the soils as reported in the literature were presented in Table 1. The observed data contains sorbed concentration of the metal species against the corresponding equilibrium concentrations. The details of the experimental data as reported in the literature are given in Table 2.

4.2. Parameter Setting. A number of twenty runs were used to test the performance of the proposed solvers. The number of iterations was fixed to 150 for all the algorithms and for all the isotherms. However, the number of agents was varied based on the number of model parameters, that is, size of the dimensional (parametric) space. The algorithms used 10, 20, and 30 agents for fitting the experimental data with Freundlich isotherm (2-dimensional), Freundlich-Langmuir isotherm (3-dimensional), and two-site Langmuir isotherm (4-dimensional), respectively.

The parameter setting for different algorithms was determined empirically on synthetic test data for which the optima were known. The tuning parameters of both the PSO algorithms, such as acceleration coefficient, inertia factor, and constriction coefficient, are dependent on the given problem [11-14, 18]. The coefficients [C.sub.1] = [C.sub.2] = 2.01 x rand(); [omega] = rand() and [chi] = 1.0 provided the best performance for the present problem. The same parameter setting was used for both the PSO algorithms. The perturbation coefficient of [[phi].sub.3] = 0.1 in (12) provides good convergence for the algorithms [14, 18]. PCPSO algorithm required additional tuning parameters, namely, catIter and number of catfish particles. The 20% worst-fit particles were considered as the catfish particles and repositioned using OBL in every 10 iterations.

The colony size in ABC and MABC algorithms was twice the number of agents considered for other algorithms. In other words, the number of employed bees and onlooker bees was the same and equal to the number of agents for other algorithms. In these algorithms, the function was evaluated three times in each iteration. Therefore, the number of iterations was fixed to 50 to use the same number of function evaluations (150) for each agent in all the algorithms. The scout limit was equal to the colony size multiplied by number of dimensions of the problem. The parameter, [epsilon], was varied from 2.0 to 0.5 from the beginning of the iteration to end of the iteration consistently. The acceleration coefficient, [C.sub.3], in MABC algorithm was fixed to 3.0 for the best performance based on the synthetic tests. The agents were repositioned randomly on the search space when they move out of the boundaries of the n-dimensional parametric space.

The pseudocode of the proposed MABC algorithm was given in the Appendix.

5. Results and Discussion

5.1. Parametric Space. The parametric space for the experimental data of Cr sorption on Alligator soil using nonlinear isotherm equations was computed. The parametric spaces were obtained using two-parameter models of Freundlich and Langmuir isotherms (obtained by substituting a = 1 in Freundlich-Langmuir isotherm) and are shown in Figures 1 and 2, respectively. RMSE (5) between experimentally measured data and theoretically computed data from 22,50,000 sets of model parameters ([K.sub.f], a or [S.sub.max], b) was used. The parametric spaces appeared fairly different in these two cases although the same experimental data was used. Further, these functional terrains contained unique global optimum solutions, several local optima, and a plan surface having the same RMSE values. Further, the parametric space varied with measurement errors in the data and the number of data points. Complexity of the functional terrains increases further with the parametric dimension. The conventional algorithms may converge to one of these local optima when the initial guess was used in such regions of the parametric space. The application of robust algorithms for determining global optimum solutions for the present problem is, therefore, justified.

5.2. Gradient-Based Algorithm. The performance of gradient-based algorithms is dependent on the user-supplied, initial solution (guess), and the type of error measure. The sensitivity of these algorithms on the initial guess was verified qualitatively for the present problem. The "fmincon" MATLAB[R] routine, featuring a constrained optimization of a multivariable function using interior-point technique, was used for finding the best isotherm model parameters that provided best-fit theoretical sorption data to the measured data. The interior-point technique was found to have better convergence abilities than trust-region-reflective or any other available techniques. A finite-difference numerical solution was used to determine the first- and second-order derivatives of the objective function with respect to the model parameters. The tolerance value of the objective function was fixed to 1 x [10.sub.-8] and a very large value was used for maximum number of iterations which served as the stopping criterion. Different initial guess solutions were used to obtain the converged solutions. The performance of the gradient-based algorithm for fitting two-site Langmuir isotherm on the experimental data of Cd sorption on Molokai soil was presented in Table 3. The solutions converged to several suboptimal and local solutions for different initial guess solutions. The predicted model parameters were fairly different from each other, but the RMSEs were very close to each other. Thus, finding a global-best solution was fairly complex in such situations. It is evident from these results that the performance of the gradient-based algorithm was highly sensitive to the initial guess and the convergence was impeded by the presence of local minima.

5.3. Parameter Estimation by SI Based Solvers. Model parameters of Freundlich, Freundlich-Langmuir, and two-site Langmuir isotherms on the experimental data were predicted using the developed solvers. The performance of the solvers was determined based on the mean fitness value, best-fit solution, and standard deviation (STD) computed from 20 independent runs as the SI techniques are stochastic in nature and may produce different solutions in different runs. The robustness of the algorithms is evaluated by comparing the mean fitness value with the best fitness (i.e., minimum error between theoretical and experimental data) and by determining the standard deviation. The robust solver predicts best-fit solution very close to the mean fitness solution and predicts smallest standard deviation, that is, close to zero. The model parameters obtained by different algorithms for all the 21 experiments using Freundlich, Freundlich-Langmuir, and two-site Langmuir isotherms were given separately in Tables 4(a), 5(a), and 6(a), respectively. The best and worst solutions obtained in 20 independent runs based on the objective function values were provided. The Langmuir isotherm was not used in the comparison as it was the special case of Fredlund-Langmuir isotherm with a= 1.

The mean fitness, best-fit solution, and standard deviation by different solvers for fitting the experimental data using Freundlich isotherm were presented in Table 4(b). The achieved best results were highlighted in bold. Performance comparison of PSO and PCPSO algorithms showed that PCPSO algorithm improved the mean fitness and standard deviations adequately without compromising on the best solution. The modification to PSO algorithm by perturbation equations and catfish swarms helped PCPSO algorithm to achieve repeatability of the best-fit solution in different runs. The performance comparison of PCPSO and ABC algorithms showed that the performance was fairly the same in finding the best (global) solutions on several different experimental data. However, ABC algorithm improved the performance of the solver with better predicted mean fitness values and standard deviation that suggested that the ABC algorithm outperformed the two PSO algorithms for fitting the sorption isotherms. Performance comparison of ABC and MABC showed that the performance of MABC algorithm was superior to other algorithms; the optimizer produced the same values of mean fit and best-fit solutions on all the experimental data. Thus, zero standard deviation is achieved in all the tests. The performance comparison of the algorithms for the 2-parameter Freundlich isotherm was in the order of MABC > ABC > PCPSO > PSO.

The computed mean fit, best-fit, and standard deviation by different algorithms on the experimental data using the Freundlich-Langmuir isotherm was presented in Table 5(b). The best solution in each case was highlighted in bold. The PCPSO algorithm improved the mean fitness and outperformed the classical PSO algorithm. However, the basic PSO algorithm converges to global-best solutions for many experimental data when compared to the PCPSO algorithm. The convergence of the PCPSO algorithm is relatively slow for fitting the data to Freundlich-Langmuir isotherm due to the additional exploitation and exploration steps by using perturbation equation and catfish swarm, respectively. Similarly, ABC algorithm improved the mean fitness but failed to converge to global solution in many experiments. On the other hand, MABC algorithm predicted the global-best solution for all the experimental data without compromising on the mean fit value. The computed standard deviations for all the test cases were close to zero. The MABC algorithm, therefore, balanced the exploration-exploitation components efficiently. The performance order of the algorithms in improving the mean fitness was MABC > ABC > PCPSO > PSO for the Freundlich-Langmuir isotherm. On the other hand, performance order in obtaining global-best solutions was MABC > PSO > ABC > PCPSO. The relative convergence behaviour of different solvers for fitting the experimental sorption data of Cu on McLaren soil was presented in Figure 3. The PSO and PCPSO algorithms converged quickly in the first few iteration and got trapped at the local optima. The convergence profiles of ABC and MABC algorithms were fairly similar and continued the search until the termination.

The performance measures by different solvers on the experimental data using the two-site Langmuir isotherm using 20 independent runs were presented in Table 6(b). The predicted best solutions were highlighted in bold. The results showed that the performance in terms of mean fitness value by PCPSO and ABC algorithms was improved slightly compared to the PSO algorithm. However, this improvement by the PCPSO and ABC algorithms was not significant for fitting the two-site Langmuir isotherm on the experimental data. Conversely, the performance of PCPSO and ABC algorithms was inferior to PSO algorithm in finding the best-fit solutions. The convergence behaviour of gbest agents against number of iterations for PSO, PCPSO, and ABC algorithms for 10 independent runs were given in Figures 4-6, respectively. The observed sorption data of Zn on Kula soil was used. PSO poorly established the repeatability of the best-fit solutions (presence of very thick band of fitness solutions) and converged to multiple solutions in different runs. In some runs, the optimizer prematurely converged to the local optima. On the other hand, the convergence of PCPSO and ABC algorithms was similar but predicted an improved mean fitness solution. However, these algorithms were not robust in converging to global-best solutions. Performance of MABC algorithm showed that the algorithm obtained global-best solutions on all the data without compromising on the predicted mean best fitness. In case of Pb sorption on Cecil and Lafitte soil and Cd sorption on Calciorthid soils, the PSO outperformed the MABC algorithm in predicting the best-fit solutions. However, the predicted model parameters and RMSEs were practically identical. The convergence behaviour of gbest agent with respect to the number of iterations for MABC algorithm for Zn sorption on Kula soil in 10 independent runs was presented in Figure 7. The figure showed that gbest agent quickly converged to the global-best solutions in all runs by establishing a small mean fitness value. Further, the computed standard deviations for all the test cases were nearly zero. The algorithms performance order in improving the mean fitness was MABC > ABC = PCPSO = PSO. On the other hand, performance order of algorithms in locating global solutions was MABC > PSO > ABC > PCPSO. The relative convergence behaviour of different optimizers on the experimental sorption data of Zn on Alligator soil was presented in Figure 8. The PSO and PCPSO algorithms quickly converged in the beginning but were trapped in the local optimum solutions. The convergence profiles of ABC and MABC optimizers were similar and continued the search process until the termination. However, MABC solver outperformed all the other solvers and provided excellent convergence. Figure 9 compares theoretical sorption data obtained using the isotherm parameters predicted by different solvers for the observed sorption data of Zn on Alligator soil. The comparison clearly showed that the MABC optimizer provided an excellent theoretical fit to the observed data. Theoretical solutions based on the predicted model parameters by MABC algorithm were used to fit the experimental data. Theoretical and experimental data for Lead sorption on Cecil soil and Chromium sorption on Alligator soil were given in Figures 10-11, respectively. Thus, the proposed MABC solver provided an excellent performance in fitting all three isotherms on all the experimental data.

6. Isotherm Selection

Performance comparison of different solvers for the estimation of isotherm model parameters from the experimental data was presented earlier. The suitable isotherms are often selected in the contaminant transport studies in soils by comparing optimized fitness values obtained by fitting different isotherms to the sorption data. The best-fit isotherms for the observed data based on the computed mean fitness values by different solvers were presented in Table 7. The suitable isotherms were highlighted in bold that provided lowest RMSE values by different solvers. It was clear that all the solvers except MABC algorithms failed to predict select appropriate isotherm due to premature convergence of the solvers while estimating the model parameters. The MABC algorithm selected suitable isotherms for all the experimental data due to excellent performance in the parameter estimation. This clearly explained the robustness of the MABC algorithm in the selection of appropriate isotherm. Similarly, the selection of appropriate isotherm based on best-fit values, obtained over 20 independent runs, was presented in Table 8. The results were similar to the former case, where the algorithms other than MABC failed to predict the appropriate isotherm for some experimental data. The results indicated that the MABC algorithm outperformed all the other solvers for finding the accurate model parameters and in the selection of appropriate isotherm for the observed data. Further, the solver is robust on various problems that consist of different parametric dimensions.

7. Summary and Conclusions

In this paper, we analyzed conventional and stochastic optimization algorithms for fitting the nonlinear isotherm models to the experimental sorption data. The disadvantage of using conventional optimization techniques is demonstrated for the present problem. These techniques were found to be highly sensitive to the initial solution and the performance is impeded by the presence of local minima. Inverse models based on swarm intelligence are developed for fitting Freundlich, Freundlich-Langmuir, and two-site Langmuir isotherm models to the observed sorption data after tuning the model parameters. However, the estimated solutions vary in different runs due to the stochastic nature of these algorithms. It was observed that all the tested optimization algorithms used from the literature are relatively slow in nature as they do not balance the exploitation and exploration components effectively. A modified bee-colony optimization algorithm is proposed in this work. Comparative analysis reveals that the proposed solver, MABC, predicts excellent results on all the experiments. The results showed that the proposed solver has superior convergence capabilities due to good exploration-exploitation components in estimating global best, mean fitness, and standard deviations. Further, the proposed solver is proved to be robust on various functional terrains. The best-fit isotherms for observed data based on the computed mean fitness values by different solvers are also obtained. The improvement in the performance of the proposed solver is significantly better and reliable for studying the fate of contaminants in soils. Therefore, a great caution should be exercised when the existing solvers are used for the selection and configuration of nonlinear isotherms.

Appendix

See Pseudocode 1.
Pseudocode 1: Pseudocode of the proposed algorithm.

Set Iteration = 0;
INITIALIZE
Initialize food sources (equation (10));
WHILE (the termination conditions are not met)
    % Employed Bee phase
    Evaluate nectar amount (fitness) of food sources (equation (5));
    FOR (each employed bee)
      Produce a new food source using (15);
      Evaluate the fitness of the new food source (equation (5));
      Greedy selection to choose best food source;
    ENDFOR
    Compute probability for each food source (equation (14));
    % Onlooker Bee's phase
    FOR (each onlooker bee)
      Send each onlooker bees to food sources based on P;
      Produce a new food source using (13);
    Evaluate the fitness of the new food source (equation (5));
    Greedy selection to choose best food source;
    END FOR
    % Scout Bee's phase
    IF (any employed bee becomes scout bee)
        Send scout bee to a new food source based on OBL (equation
          (16));
    ENDIF
    Memorize the best solution achieved so far;
    Iteration = Iteration + 1;
END WHILE


http://dx.doi.org/10.1155/2017/3497652

Competing Interests

The author declares that they have no competing interests.

References

[1] T. L. Theis, "Reactions and transport of trace metals in groundwater," in Metal Speciation: Theory, Analysis, and Application, J. R. Kramer and H. E. Allen, Eds., CRC/Lewis, Boca Raton, Fla, USA, 1988.

[2] H. M. Selim and M. C. Amacher, Reactivity and Transport of Heavy Metals in Soils, CRC Press, Boca Raton, Fla, USA, 1997.

[3] K. V. Kumar and S. Sivanesan, "Prediction ofoptimum sorption isotherm: comparison of linear and non-linear method," Journal of Hazardous Materials, vol. 126, no. 1-3, pp. 198-201, 2005.

[4] G. P. Jeppu and T. P. Clement, "A modified Langmuir-Freundlich isotherm model for simulating pH-dependent adsorption effects," Journal of Contaminant Hydrology, vol. 129-130, pp. 46-53, 2012.

[5] C. H. Bolster and G. M. Hornberger, "On the use of linearized langmuir equations," Soil Science Society of America Journal, vol. 71, no. 6, pp. 1796-1806, 2007

[6] W. R. Roy, I. G. Krapac, S. F. J. Chou, and R. A. Griffin, Batch-Type Procedures for Estimating Soil Adsorption of Chemicals, EPA/530/SW-87/006-F, U.S. Environmental Protection Agency, Washington, DC, USA, 1992.

[7] R. M. Allen-King, P. Grathwohl, and W. P Ball, "New modeling paradigms for the sorption of hydrophobic organic chemicals to heterogeneous carbonaceous matter in soils, sediments, and rocks," Advances in Water Resources, vol. 25, no. 8-12, pp. 985-1016, 2002.

[8] Y.-S. Ho, "Selection of optimum sorption isotherm," Carbon, vol. 42, no. 10, pp. 2115-2116, 2004.

[9] J. Kennedy and R. C. Eberhart, "Particle swarm optimization," in Proceedings of the IEEE International Conference on Neural Networks, pp. 1942-1948, Institute of Electrical and Electronics Engineers, Perth, Australia, 1995.

[10] D. Karaboga, "An idea based on honeybee swarm for numerical optimization," Tech. Rep. TR06, Erciyes University, Engineering Faculty, Computer Engineering Department, Kayseri, Turkey, 2005.

[11] T. V. Bharat, "Agents based algorithms for design parameter estimation in contaminant transport inverse problems," in Proceedings of the IEEE Swarm Intelligence Symposium (SIS '08), pp. 1-7, Missouri, USA, September 2008.

[12] T. V. Bharat, P V. Sivapullaiah, and M. M. Allam, "Swarm intelligence-based solver for parameter estimation of laboratory through-diffusion transport of contaminants," Computers and Geotechnics, vol. 36, no. 6, pp. 984-992, 2009.

[13] T. V. Bharat, P. V. Sivapullaiah, and M. M. Allam, "Swarm intelligence based inverse model for characterization of groundwater contaminant source," Electronic Journal of Geotechnical Engineering, vol. 14, pp. 1-14, 2009.

[14] T. V. Bharat, P V. Sivapullaiah, and M. M. Allam, "Robust solver based on modified particle swarm optimization for improved solution of diffusion transport through containment facilities," Expert Systems with Applications, vol. 39, no. 12, pp. 10812-10820, 2012.

[15] Y. M. Cheng, L. Li, S.-C. Chi, and W. B. Wei, "Particle swarm optimization algorithm for the location of the critical non-circular failure surface in two-dimensional slope stability analysis," Computers and Geotechnics, vol. 34, no. 2, pp. 92-103, 2007.

[16] R. Kalatehjari, N. Ali, M. Hajihassani, and M. Kholghi Fard, "The application of particle swarm optimization in slope stability analysis of homogeneous soil slopes," International Review on Modelling and Simulations, vol. 5, no. 1, pp. 458-465, 2012.

[17] W. Gao, "Back analysis algorithm in geotechnical engineering based on particle swarm optimization," Yantu Lixue/Rock and Soil Mechanics, vol. 27, no. 5, pp. 795-798, 2006.

[18] T. V. Bharat, P. V. Sivapullaiah, and M. M. Allam, "Accurate parameter estimation of contaminant transport inverse problem using particle swarm optimization," in Proceedings of the IEEE Swarm Intelligence Symposium (SIS '08), pp. 1-7, IEEE, St. Louis, Mo, USA, September 2008.

[19] L. S. Matott and A. J. Rabideau, "ISOFIT--a program for fitting sorption isotherms to experimental data," Environmental Modelling and Software, vol. 23, no. 5, pp. 670-676, 2008.

[20] S. Katare, A. Kalos, and D. West, "A hybrid swarm optimizer for efficient parameter estimation," in Proceedings of the IEEE Congress on Evolutionary Computation, pp. 309-315, Portland, Ore, USA, June 2004.

[21] T. Ying, Y. P. Yang, and J. C. Zeng, "An enhanced hybrid quadratic particle swarm optimization," in Proceedings of the 6th International Conference on Intelligent Systems Design and Applications (ISDA '06), pp. 980-985, Jinan, China, October 2006.

[22] C. H. Giles, D. Smith, and A. Huitson, "A general treatment and classification of the solute adsorption isotherm. I. Theoretical," Journal of Colloid And Interface Science, vol. 47, no. 3, pp. 755-765, 1974.

[23] G. Sposito, The Surface Chemistry of Soils, Oxford University Press, New York, NY, USA, 1984.

[24] C. Hinz, "Description of sorption data with isotherm equations," Geoderma, vol. 99, no. 3-4, pp. 225-243, 2001.

[25] V. Murali and L. A. G. Aylmore, "Competitive adsorption during solute transport in soils: 1. Mathematical models," Soil Science, vol. 135, no. 3, pp. 143-150, 1983.

[26] H. M. Selim, M. C. Amacher, and I. K. Iskandar, "Modeling the transport of heavy metals in soils," CRREL-Monograph 90-2, U.S. Army Corps of Engineers, Hanover, NH, USA, 1990.

[27] K. Y. Foo and B. H. Hameed, "Insights into the modeling of adsorption isotherm systems," Chemical Engineering Journal, vol. 156, no. 1, pp. 2-10, 2010.

[28] O. S. Pokrovsky, A. Probst, E. Leviel, and B. Liao, "Interactions between cadmium and lead with acidic soils: experimental evidence of similar adsorption patterns for a wide range of metal concentrations and the implications of metal migration," Journal of Hazardous Materials, vol. 199-200, pp. 358-366, 2012.

[29] I. Langmuir, "The adsorption of gases on plane surfaces of glass, mica and platinum," The Journal of the American Chemical Society, vol. 40, no. 9, pp. 1361-1403, 1918.

[30] C. G. Rossi, D. M. Heil, N. B. Bonuma, and J. R. Williams, "Evaluation of the Langmuir model in the Soil and Water Assessment Tool for a high soil phosphorus condition," Environmental Modelling and Software, vol. 38, pp. 40-49, 2012.

[31] R. Malekian, J. Abedi-Koupai, S. S. Eslamian, S. F. Mousavi, K. C. Abbaspour, and M. Afyuni, "Ion-exchange process for ammonium removal and release using natural Iranian zeolite," Applied Clay Science, vol. 51, no. 3, pp. 323-329, 2011.

[32] B. Subramanyam and D. Ashutosh, "Adsorption isotherm modeling of phenol onto natural soils--applicability of various isotherm models," International Journal of Environmental Research, vol. 6, no. 1, pp. 265-276, 2012.

[33] O. Altin, H. O. Ozbelge, and T. Dogu, "Use of general purpose adsorption isotherms for heavy metal-clay mineral interactions," Journal of Colloid and Interface Science, vol. 198, no. 1, pp. 130-140, 1998.

[34] H. B. Bradl, "Adsorption of heavy metal ions on soils and soils constituents," Journal of Colloid and Interface Science, vol. 277, no. 1, pp. 1-18, 2004.

[35] C. C. Travis and E. L. Etnier, "A survey of sorption relationships for reactive solutes in soil," Journal of Environmental Quality, vol. 10, no. 1, pp. 8-17, 1981.

[36] I. C. R. Holford, R. W. M. Wedderburn, and G. E. G. Mattingly, "A Langmuir two-surface equation as a model for phosphate adsorption by soils," Journal of Soil Science, vol. 25, no. 2, pp. 242-255, 1974.

[37] M. Vidal, M. J. Santos, T. Abrao, J. Rodriguez, and A. Rigol, "Modeling competitive metal sorption in a mineral soil," Geoderma, vol. 149, no. 3-4, pp. 189-198, 2009.

[38] E. Galunin, M. D. Alba, M. J. Santos, and M. Vidal, "Effects of the presence of Fe(0) on the sorption of lanthanum and lutetium mixtures in smectites," Applied Clay Science, vol. 65-66, pp. 162-172, 2012.

[39] K. J. Farley, D. A. Dzombak, and F. M. M. Morel, "A surface precipitation model for the sorption of cations on metal oxides," Journal of Colloid and Interface Science, vol. 106, no. 1, pp. 226-242, 1985.

[40] J. Toth, "Uniform interpretation of gas/solid adsorption," Advances in Colloid and Interface Science, vol. 55, pp. 1-239, 1995.

[41] M. I. El-Khaiary, "Least-squares regression of adsorption equilibrium data: comparing the options," Journal of Hazardous Materials, vol. 158, no. 1, pp. 73-87, 2008.

[42] B. Subramanyam and A. Das, "Linearized and non-linearized isotherm models comparative study on adsorption of aqueous phenol solution in soil," International Journal of Environmental Science & Technology, vol. 6, no. 4, pp. 633-640, 2009.

[43] R. Eberhart and J. Kennedy, "New optimizer using particle swarm theory," in Proceedings of the 6th International Symposium on Micro Machine and Human Science, pp. 39-43, Nagoya, Japan, October 1995.

[44] Y. Shi and R. Eberhart, "Modified particle swarm optimizer," in Proceedings of the IEEE International Conference on Evolutionary Computation (ICEC '98), pp. 69-73, Piscataway, NJ, USA, May 1998.

[45] M. Clerc and J. Kennedy, "The particle swarm--explosion, stability, and convergence in a multidimensional complex space," IEEE Transactions on Evolutionary Computation, vol. 6, no. 1, pp. 58-73, 2002.

[46] L.-Y. Chuang, S.-W. Tsai, and C.-H. Yang, "Catfish particle swarm optimization," in Proceedings of the IEEE Swarm Intelligence Symposium (SIS '08), St. Louis, Mo, USA, September 2008.

[47] X. Z. Gao, X. Wang, S. J. Ovaska, and K. Zenger, "A hybrid optimization method of harmony search and opposition-based learning," Engineering Optimization, vol. 44, no. 8, pp. 895-914, 2012.

[48] D. Karaboga and B. Akay, "A comparative study of artificial Bee colony algorithm," Applied Mathematics and Computation, vol. 214, no. 1, pp. 108-132, 2009.

[49] G. Zhu and S. Kwong, "Gbest-guided artificial bee colony algorithm for numerical function optimization," Applied Mathematics and Computation, vol. 217, no. 7, pp. 3166-3173, 2010.

[50] B. Akay and D. Karaboga, "A modified Artificial Bee Colony algorithm for real-parameter optimization," Information Sciences, vol. 192, pp. 120-142, 2012.

[51] W.-F. Gao and S.-Y. Liu, "A modified artificial bee colony algorithm," Computers & Operations Research, vol. 39, no. 3, pp. 687-697, 2012.

[52] W. Gao and S. Liu, "Improved artificial bee colony algorithm for global optimization," Information ProcessingLetters, vol. 111, no. 17, pp. 871-882, 2011.

[53] D. Karaboga and E. Kaya, "An adaptive and hybrid artificial bee colony algorithm (aABC) for ANFIS training," Applied Soft Computing, vol. 49, pp. 423-436, 2016.

[54] A. Chatterjee and P. Siarry, "Nonlinear inertia weight variation for dynamic adaptation in particle swarm optimization," Computers & Operations Research, vol. 33, no. 3, pp. 859-871, 2006.

Tadikonda Venkata Bharat

Department of Civil Engineering, Indian Institute of Technology Guwahati, Guwahati, Assam 781039, India

Correspondence should be addressed to Tadikonda Venkata Bharat; tvb@iitg.ernet.in

Received 23 March 2016; Revised 10 October 2016; Accepted 22 November 2016; Published 18 January 2017

Academic Editor: David Glass

Caption: Figure 1: Parametric space of Freundlich isotherm for Cr sorption on Alligator soil.

Caption: Figure 2: Parametric space of Langmuir isotherm for Cr sorption on Alligator soil.

Caption: Figure 3: Convergence of gbest agent in PSO (*), PCPSO (o), ABC (**), and MABC (***) optimizers for fitting Freundlich-Langmuir isotherm on the experimental data of Cu sorption on McLaren soil (iterations = Iter x function evaluations for gbest agent, for ABC and MABC optimizers).

Caption: Figure 4: Convergence of gbest particle in PSO optimizer in different independent runs for experimental sorption of Zn sorption on Kula soil.

Caption: Figure 5: Convergence of gbest particle in PCPSO optimizer in different independent runs for experimental sorption of Zn sorption on Kula soil.

Caption: Figure 6: Convergence of gbest agent in ABC optimizer in different independent runs for experimental sorption of Zn sorption on Kula soil (here, iterations = Iter x function evaluations for gbest agent).

Caption: Figure 7: Convergence of gbest agent in MABC optimizer in different independent runs for experimental sorption of Zn sorption on Kula soil (here, iterations = Iter x function evaluations for gbest agent).

Caption: Figure 8: Convergence of gbest agent in PSO (*), PCPSO (o), ABC (**), and MABC (***) optimizers for fitting two-site Langmuir isotherm on the experimental data of Zn sorption on Alligator soil (iterations = Iter x function evaluations for gbest agent, for ABC and MABC optimizers).

Caption: Figure 9: Theoretical fits of two-site Langmuir isotherm to experimental data of Zn sorption on Alligator soil by different solvers.

Caption: Figure 10: Theoretical fit of Freundlich-Langmuir isotherm to the experimental data of Pb sorption on Cecil soil by MABC solver.

Caption: Figure 11: Theoretical fit of Freundlich-Langmuir isotherm to the experimental data of Cr sorption on Alligator soil by MABC solver.
Table 1: Description of the experimental data.

#    Soil          Reference                      pH

1    Alligator     Buchter et al., 1989           4.8
2    McLaren       Selim and Amacher, 1997        NA
3    Cecil         Selim and Amacher, 1997        5.4
4    Kula          Selim and Amacher, 1997        5.9
5    Webster       Selim and Amacher, 1997        7.6
6    Lafitte       Selim and Amacher, 1997        3.9
7    Windsor       Buchter et al., 1989           5.3
8    Molokai       Buchter et al., 1989           6.0
9    Spodosol      Buchter et al., 1989           4.3
10   Calciorthid   Buchter et al., 1989           8.5
11   Hartsells     Bolster and Hornberger, 2007   NA
12   Pembroke      Bolster and Hornberger, 2007   NA
13   Loring        Bolster and Hornberger, 2007   NA

#    Soil type                CEC        Grain size (%)

                                     Sand   Silt   Clay

1    Clay, Montmorillo-nitic  30.3   5.9    39.4   54.7
2    NA                       NA     NA     NA     NA
3    Clay, Kaolinitic         2.4    30     18.8   51.2
4    Sandy Loam               27     66.6   32.9   0.5
5    Clay Loam                14.1   27.5   48.6   23.9
6    Sandy Loam               26.9   60.7   21.7   17.6
7    Sandy Loam               2.0    76.8   20.5   2.8
8    Clay Loam, Kaolinitc     11     25.7   46.2   28.2
9    Sand                     2.7    90.2   6.0    3.8
10   Sandy Loam               14.7   70     19.3   10.7
11   Fine Loamy               NA     NA     NA     NA
12   Fine Silty               NA     NA     NA     NA

13   Fine Silty               NA     NA     NA     NA

Table 2: Details of the input data.

#    Reactive    Soil          Number     Reference
     metal                     of data
     substance                 points

1       Cu        Alligator       7         Selim and Amacher, 1997
2       Cu         McLaren        5         Selim and Amacher, 1997
3       Cu          Cecil         5         Selim and Amacher, 1997
4       Zn        Alligator       10        Selim and Amacher, 1997
5       Zn          Kula          7         Selim and Amacher, 1997
6       Zn         Webster        8         Selim and Amacher, 1997
7       Pb          Cecil         6         Selim and Amacher, 1997
8       Pb        Spodosol        5         Selim and Amacher, 1997
9       Pb         Lafitte        7         Selim and Amacher, 1997
10      Pb        Alligator       8         Selim and Amacher, 1997
11      Cr        Alligator       15          Buchter et al., 1989
12      Cr         Windsor        12          Buchter et al., 1989
13      Cr          Kula          10          Buchter et al., 1989
14      Cd         Molokai        13          Buchter et al., 1989
15      Cd          Kula          10          Buchter et al., 1989
16      Cd         Windsor        14          Buchter et al., 1989
17      Cd        Spodosol        14          Buchter et al., 1989
18      Cd       Calciorthid      7           Buchter et al., 1989
19      Ca        Hartsells       6       Bolster and Hornberger, 2007
20      Ca        Pembroke        6       Bolster and Hornberger, 2007
21      Ca         Loring         6       Bolster and Hornberger, 2007

Table 3: Performance of the gradient-based solver on the experimental
data of Cd sorption on Molokai soil using two-site Langmuir isotherm.

                      Initial guess

[mathematical    [b.sub.1]   [mathematical    [b.sub.2]
expression not               expression not
reproducible]                reproducible]

920              0.117       23.7             4.4
880              0.147       11.5             7
100              0.1         10               7
100              0.1         10               0.5
100              0.1         100              0.5

                      Predicted solution                  Error

[mathematical    [b.sub.1]   [mathematical    [b.sub.2]
expression not               expression not
reproducible]                reproducible]

923.84           0.119       21.1746          4.884       0.071872
880.29           0.130       19.5152          5.054       0.072161
490.211          0.320       356.701          0.075       0.089178
515.777          0.308       324.213          0.075       0.089366
429.705          0.350       429.915          0.079       0.088687

Table 4: (a) Best and Worst solutions (set of model parameters)
achieved by different solvers in 20 independent runs f[O.sub.R]
fitting the experimental data on Freundlich isotherm model. (b)
Performance of different solvers for fitting Freundlich isotherm to
the experimental data.

                         (a)

#    Parameters         PSO            PCPSO

                  Best    Worst   Best    Worst

1    [K.sub.f]    252.6   228.5   252.6   244.2
     a            0.533   0.518   0.533   0.501
     [O.sub.R]    0.145   0.151   0.145   0.148
2    [K.sub.f]    28.40   635.4   28.51   22.21
     a            0.729   1.91    0.729   0.796
     [O.sub.R]    0.111   2.499   0.111   0.126
3    [K.sub.f]    22.25   4512    22.25   14.62
     a            0.559   0.006   0.559   0.748
     [O.sub.R]    0.023   1.753   0.023   0.114
4    [K.sub.f]    28.50   4.45    28.50   42.81
     a            0.994   4.19    0.994   0.973
     [O.sub.R]    0.080   2.49    0.080   0.188
5    [K.sub.f]    231.1   226.9   231.1   201.7
     a            0.728   1.601   0.728   0.739
     [O.sub.R]    0.040   0.640   0.040   0.073
6    [K.sub.f]    736.8   309.8   736.8   871.9
     a            0.679   1.155   0.679   0.725
     [O.sub.R]    0.097   1.191   0.097   0.111
7    [K.sub.f]    265.7   234.9   265.7   1480
     a            0.695   1.409   0.695   1.437
     [O.sub.R]    0.059   0.544   0.059   0.819
8    [K.sub.f]    126.6   64.57   126.6   126.9
     a            0.750   1.693   0.750   0.750
     [O.sub.R]    0.085   0.665   0.085   0.085
9    [K.sub.f]    994.7   1176    994.7   977.2
     a            1.645   3.581   1.645   1.638
     [O.sub.R]    0.133   1.080   0.133   0.133
10   [K.sub.f]    1848    4872    1848    6063
     a            0.861   3.152   0.861   1.216
     [O.sub.R]    0.111   2.583   0.111   0.283
11   [K.sub.f]    3.493   37.66   3.493   3.027
     a            0.501   0.515   0.501   0.541
     [O.sub.R]    0.091   1.041   0.091   0.116
12   [K.sub.f]    8.519   166.7   8.519   8.451
     a            0.503   0.945   0.503   0.447
     [O.sub.R]    0.120   1.400   0.120   0.134
13   [K.sub.f]    63.61   32.06   63.61   44.67
     a            0.605   2.680   0.605   0.639
     [O.sub.R]    0.083   2.127   0.083   0.169
14   [K.sub.f]    86.28   832.5   86.28   83.87
     a            0.744   1.090   0.744   0.731
     [O.sub.R]    0.136   0.997   0.136   0.137
15   [K.sub.f]    194.4   1365    194.4   176.9
     a            0.738   1.124   0.738   0.725
     [O.sub.R]    0.108   0.681   0.108   0.113
16   [K.sub.f]    14.95   4873    14.95   18.51
     a            0.787   5.339   0.787   0.783
     [O.sub.R]    0.094   5.300   0.094   0.134
17   [K.sub.f]    5.301   64.29   5.301   5.595
     a            0.835   1.063   0.835   0.859
     [O.sub.R]    0.068   1.083   0.068   0.075
18   [K.sub.f]    407.7   979.9   407.7   459.0
     a            0.602   0.940   0.602   0.669
     [O.sub.R]    0.079   0.347   0.079   0.103
19   [K.sub.f]    81.69   99.29   81.69   92.05
     a            0.264   1.434   0.264   0.283
     [O.sub.R]    0.051   1.052   0.051   0.080
20   [K.sub.f]    43.70   1.428   43.70   11.57
     a            0.419   2.678   0.419   1.204
     [O.sub.R]    0.029   1.121   0.029   0.379
21   [K.sub.f]    83.26   33.30   83.26   86.52
     a            0.312   1.010   0.312   0.305
     [O.sub.R]    0.031   0.468   0.031   0.034

                         (a)

#    Parameters         ABC            MABC

                  Best    Worst   Best    Worst

1    [K.sub.f]    252.6   262.2   252.6   252.6
     a            0.533   0.598   0.533   0.533
     [O.sub.R]    0.145   0.158   0.145   0.145
2    [K.sub.f]    28.41   29.95   28.41   28.37
     a            0.729   0.706   0.729   0.729
     [O.sub.R]    0.111   0.112   0.111   0.111
3    [K.sub.f]    22.25   26.06   22.25   22.41
     a            0.559   0.50    0.559   0.556
     [O.sub.R]    0.023   0.042   0.023   0.023
4    [K.sub.f]    28.50   28.70   28.50   28.50
     a            0.994   0.985   0.994   0.994
     [O.sub.R]    0.080   0.080   0.080   0.080
5    [K.sub.f]    231.1   231.1   231.1   231.1
     a            0.728   0.728   0.728   0.728
     [O.sub.R]    0.040   0.040   0.040   0.040
6    [K.sub.f]    736.8   788.9   736.8   736.8
     a            0.679   0.701   0.679   0.679
     [O.sub.R]    0.097   0.100   0.097   0.097
7    [K.sub.f]    265.7   265.6   265.7   265.7
     a            0.695   0.695   0.695   0.695
     [O.sub.R]    0.059   0.059   0.059   0.059
8    [K.sub.f]    126.6   126.6   126.6   126.6
     a            0.750   0.750   0.750   0.750
     [O.sub.R]    0.085   0.085   0.085   0.085
9    [K.sub.f]    994.7   1414    994.7   994.7
     a            1.645   1.859   1.645   1.645
     [O.sub.R]    0.133   0.156   0.133   0.133
10   [K.sub.f]    1848    2993    1848    1846
     a            0.861   1.049   0.861   0.860
     [O.sub.R]    0.111   0.164   0.111   0.111
11   [K.sub.f]    3.493   3.606   3.493   3.493
     a            0.501   0.460   0.501   0.501
     [O.sub.R]    0.091   0.104   0.091   0.091
12   [K.sub.f]    8.519   8.518   8.519   8.519
     a            0.503   0.502   0.503   0.503
     [O.sub.R]    0.120   0.120   0.120   0.120
13   [K.sub.f]    63.61   63.66   63.61   63.61
     a            0.605   0.605   0.605   0.605
     [O.sub.R]    0.083   0.083   0.083   0.083
14   [K.sub.f]    86.28   96.99   86.28   86.28
     a            0.744   0.767   0.744   0.744
     [O.sub.R]    0.136   0.145   0.136   0.136
15   [K.sub.f]    194.4   194.4   194.4   194.4
     a            0.738   0.738   0.738   0.738
     [O.sub.R]    0.108   0.108   0.108   0.108
16   [K.sub.f]    14.95   31.20   14.95   14.95
     a            0.787   1.348   0.787   0.787
     [O.sub.R]    0.094   0.661   0.094   0.094
17   [K.sub.f]    5.301   5.289   5.301   5.301
     a            0.835   0.835   0.835   0.835
     [O.sub.R]    0.068   0.068   0.068   0.068
18   [K.sub.f]    407.7   483.14  407.7   407.7
     a            0.602   0.665   0.602   0.602
     [O.sub.R]    0.079   0.101   0.079   0.079
19   [K.sub.f]    81.69   69.36   81.69   81.69
     a            0.264   0.375   0.264   0.264
     [O.sub.R]    0.051   0.098   0.051   0.051
20   [K.sub.f]    43.70   44.22   43.70   43.73
     a            0.419   0.413   0.419   0.418
     [O.sub.R]    0.029   0.029   0.029   0.029
21   [K.sub.f]    83.26   86.20   83.26   83.26
     a            0.312   0.294   0.312   0.312
     [O.sub.R]    0.031   0.034   0.031   0.031

                             (b)

#            PSO algorithm                   PCPSO algorithm

     Mean      Best      Std           Mean      Best      Std

1    0.3555    0.1449    3.91E - 01    0.1452    0.1449    4.11E - 05
2    0.4190    0.1112    6.79E - 01    0.1119    0.1112    1.30E - 03
3    0.3079    0.0233    4.57E - 01    0.0311    0.0233    9.85E - 03
4    0.2045    0.0798    3.08E - 01    0.0819    0.0798    6.05E - 03
5    0.0614    0.0404    6.30E - 02    0.0404    0.0404    1.41E - 05
6    0.1143    0.0973    5.11E - 02    0.0973    0.0973    1.85E - 05
7    0.0809    0.0589    6.62E - 02    0.0589    0.0589    1.13E - 04
8    0.1476    0.0852    1.53E - 01    0.0853    0.0852    4.40E - 04
9    0.1537    0.1330    4.59E - 02    0.1330    0.1330    1.69E - 05
10   0.1465    0.1111    8.61E - 02    0.1119    0.1111    2.42E - 03
11   0.4355    0.0911    7.65E - 01    0.1335    0.0912    3.72E - 02
12   0.2852    0.1198    4.96E - 01    0.1450    0.1198    2.74E - 02
13   0.6264    0.0826    1.11E + 00    0.0881    0.0825    5.61E - 03
14   0.4260    0.1355    5.63E - 01    0.1397    0.1355    6.97E - 03
15   0.1828    0.1081    2.24E - 01    0.1142    0.1081    5.65E - 03
16   0.1155    0.0944    6.32E - 02    0.1185    0.0944    2.00E - 02
17   0.6646    0.0676    9.28E - 01    0.1019    0.0676    2.36E - 02
18   0.1821    0.0790    1.74E - 01    0.0891    0.0790    1.51E - 02
19   0.2400    0.0509    3.83E - 01    0.0658    0.0509    1.60E - 02
20   0.1677    0.0292    2.62E - 01    0.0437    0.0292    2.05E - 02
21   0.4273    0.0310    6.15E - 01    0.0444    0.0310    1.48E - 02

                             (b)

#             ABC algorithm                   MABC algorithm

     Mean      Best      Std           Mean      Best      Std

1    0.1457    0.1449    2.43E - 03    0.1449    0.1449    0.00E + 00
2    0.1226    0.1112    2.60E - 02    0.1112    0.1112    0.00E + 00
3    0.0233    0.0233    6.44E - 05    0.0233    0.0233    2.10E - 07
4    0.0827    0.0798    7.70E - 03    0.0798    0.0798    0.00E + 00
5    0.0404    0.0404    0.00E + 00    0.0404    0.0404    0.00E + 00
6    0.0979    0.0973    1.89E - 03    0.0973    0.0973    0.00E + 00
7    0.0589    0.0589    0.00E + 00    0.0589    0.0589    0.00E + 00
8    0.0852    0.0852    0.00E + 00    0.0852    0.0852    0.00E + 00
9    0.1330    0.1330    5.09E - 06    0.1330    0.1330    1.51E - 06
10   0.1113    0.1111    1.93E - 04    0.1010    0.1010    6.61E - 06
11   0.0915    0.0911    1.13E - 03    0.0911    0.0911    0.00E + 00
12   0.1198    0.1198    1.00E - 07    0.1198    0.1198    0.00E + 00
13   0.0826    0.0826    4.20E - 07    0.0825    0.0825    0.00E + 00
14   0.1355    0.1355    1.00E - 07    0.1355    0.1355    0.00E + 00
15   0.1081    0.1081    1.12E - 05    0.1081    0.1081    0.00E + 00
16   0.0944    0.0944    1.86E - 06    0.0944    0.0944    0.00E + 00
17   0.0677    0.0677    1.70E - 07    0.0676    0.0676    0.00E + 00
18   0.0810    0.0790    5.88E - 03    0.0790    0.0790    0.00E + 00
19   0.0531    0.0509    6.04E - 03    0.0509    0.0509    0.00E + 00
20   0.0378    0.0292    2.51E - 02    0.0292    0.0292    2.58E - 06
21   0.0591    0.0310    8.42E - 02    0.0176    0.0176    0.00E + 00

Table 5: (a) Best and worst solutions (set of model parameters)
achieved by different solvers in 20 independent runs for fitting the
experimental data on Freundlich-Langmuir isotherm model. (b)
Performance of different solvers for fitting Langmuir-Freundlich
Isotherm to the experimental data.

                        (a)

#    Parameters         PSO             PCPSO

                   Best    Worst   Best    Worst

1    [S.sub.max]   1174    358.2   1166    574.8
     b             0.403   5.403   0.408   2.195
     a             0.779   0.998   0.781   0.987
     [O.sub.R]     0.128   0.306   0.128   0.202
2    [S.sub.max]   711.9   143.3   668.0   403.2
     b             0.034   4.400   0.036   0.133
     a             0.988   0.329   0.972   0.521
     [O.sub.R]     0.093   0.406   0.095   0.243
3    [S.sub.max]   645.0   91.33   791.3   93.47
     b             0.031   7.348   0.029   4.705
     a             0.688   0.418   0.615   0.520
     [O.sub.R]     0.023   0.315   0.024   0.309
4    [S.sub.max]   845.8   82.25   1261    79.96
     b             0.041   2.705   0.024   4.236
     a             0.964   0.639   0.958   0.995
     [O.sub.R]     0.152   0.674   0.131   0.667
5    [S.sub.max]   2244    283.2   2334    771.3
     b             0.132   2.437   0.125   0.570
     a             0.868   0.639   0.860   0.720
     [O.sub.R]     0.031   0.399   0.031   0.217
6    [S.sub.max]   1420    380.1   1460    3503
     b             1.280   5.470   1.201   0.242
     a             0.950   0.561   0.940   0.667
     [O.sub.R]     0.074   0.381   0.074   0.107
7    [S.sub.max]   1530    317.3   1639    3547
     b             0.253   5.047   0.230   0.077
     a             0.871   0.839   0.856   0.501
     [O.sub.R]     0.046   0.358   0.046   0.173
8    [S.sub.max]   2458    210.7   1753    213.7
     b             0.058   3.177   0.084   5.192
     a             0.851   0.686   0.909   0.993
     [O.sub.R]     0.082   0.424   0.084   0.425
9    [S.sub.max]   2288    174.3   4857    237.4
     b             0.223   5.257   0.105   6.490
     a             0.999   0.401   0.973   1.000
     [O.sub.R]     0.266   0.514   0.263   0.444
10   [S.sub.max]   2327    649.8   2352    752.4
     b             1.402   1.949   1.389   6.703
     a             1.000   0.608   0.999   1.000
     [O.sub.R]     0.101   0.313   0.101   0.191
11   [S.sub.max]   45.41   8.379   40.78   7.229
     b             0.098   2.172   0.115   6.333
     a             0.617   0.639   0.634   0.999
     [O.sub.R]     0.070   0.395   0.071   0.405
12   [S.sub.max]   45.91   12.61   45.86   15.79
     b             0.342   1.831   0.343   4.944
     a             0.763   0.528   0.764   0.994
     [O.sub.R]     0.031   0.362   0.031   0.283
13   [S.sub.max]   748.9   131.3   758.3   1480
     b             0.114   6.554   0.112   0.051
     a             0.795   1.000   0.795   0.467
     [O.sub.R]     0.037   0.467   0.037   0.194
14   [S.sub.max]   951.2   96.21   895.3   125.1
     b             0.150   9.307   0.166   6.732
     a             0.912   0.702   0.922   1.000
     [O.sub.R]     0.079   0.745   0.079   0.549
15   [S.sub.max]   931.6   155.3   943.3   157.1
     b             0.403   2.035   0.395   6.723
     a             0.918   0.647   0.916   1.000
     [O.sub.R]     0.040   0.426   0.040   0.357
16   [S.sub.max]   275.9   18.57   651.8   88.24
     b             0.076   3.857   0.028   0.309
     a             0.916   0.896   0.869   0.955
     [O.sub.R]     0.082   0.573   0.069   0.206
17   [S.sub.max]   380.1   5.862   275.9   7.135
     b             0.016   7.608   0.022   1.255
     a             0.905   0.701   0.896   0.979
     [O.sub.R]     0.068   0.797   0.071   0.581
18   [S.sub.max]   1288    291.0   1288    451.1
     b             0.609   4.950   0.608   5.687
     a             0.766   0.746   0.766   1.000
     [O.sub.R]     0.032   0.278   0.032   0.159
19   [S.sub.max]   181.1   121.4   181.1   127.2
     b             1.058   6.737   1.059   6.168
     a             0.660   0.241   0.661   0.997
     [O.sub.R]     0.024   0.179   0.024   0.099
20   [S.sub.max]   244.8   100.5   239.3   101.3
     b             0.205   5.056   0.212   5.464
     a             0.661   0.741   0.666   0.435
     [O.sub.R]     0.015   0.177   0.016   0.184
21   [S.sub.max]   266.2   132.2   273.1   167.2
     b             0.501   5.629   0.480   1.656
     a             0.540   0.374   0.529   0.935
     [O.sub.R]     0.011   0.175   0.011   0.046

                        (a)

#    Parameters         ABC             MABC

                   Best    Worst   Best    Worst

1    [S.sub.max]   1176    887.7   1167    1361
     b             0.399   0.730   0.407   0.301
     a             0.779   0.905   0.780   0.725
     [O.sub.R]     0.128   0.132   0.128   0.129
2    [S.sub.max]   744.1   695.6   674.6   905.3
     b             0.032   0.041   0.036   0.026
     a             0.986   0.904   1.000   0.958
     [O.sub.R]     0.093   0.102   0.092   0.096
3    [S.sub.max]   1072    276.5   904.7   402.7
     b             0.021   0.080   0.023   0.049
     a             0.618   0.790   0.633   0.758
     [O.sub.R]     0.022   0.051   0.020   0.027
4    [S.sub.max]   3426    1586    5878    8999
     b             0.009   0.019   0.005   0.003
     a             1.000   0.999   1.000   0.999
     [O.sub.R]     0.078   0.096   0.076   0.077
5    [S.sub.max]   2259    1664    2429    2280
     b             0.131   0.201   0.119   0.129
     a             0.869   0.943   0.854   0.864
     [O.sub.R]     0.031   0.035   0.031   0.031
6    [S.sub.max]   1424    1382    1420    1434
     b             1.273   1.356   1.281   1.258
     a             0.949   0.959   0.950   0.947
     [O.sub.R]     0.074   0.074   0.074   0.074
7    [S.sub.max]   1678    1259    1690    1780
     b             0.222   0.346   0.219   0.204
     a             0.853   0.927   0.851   0.842
     [O.sub.R]     0.046   0.048   0.046   0.046
8    [S.sub.max]   2234    1640    2715    2260
     b             0.064   0.093   0.052   0.064
     a             0.856   0.940   0.841   0.865
     [O.sub.R]     0.082   0.085   0.082   0.082
9    [S.sub.max]   13333   3801    13342   13268
     b             0.036   0.136   0.036   0.036
     a             0.992   0.999   1.000   0.995
     [O.sub.R]     0.249   0.258   0.248   0.249
10   [S.sub.max]   2340    2326    2328    2390
     b             1.399   1.401   1.402   1.350
     a             1.000   0.998   1.000   1.000
     [O.sub.R]     0.101   0.101   0.101   0.101
11   [S.sub.max]   45.48   38.90   45.47   44.64
     b             0.098   0.121   0.098   0.099
     a             0.616   0.650   0.617   0.615
     [O.sub.R]     0.070   0.072   0.070   0.070
12   [S.sub.max]   45.97   43.91   45.91   46.06
     b             0.341   0.371   0.342   0.340
     a             0.762   0.775   0.763   0.761
     [O.sub.R]     0.031   0.032   0.031   0.031
13   [S.sub.max]   755.0   639.9   752.7   756.5
     b             0.113   0.141   0.113   0.112
     a             0.793   0.833   0.794   0.792
     [O.sub.R]     0.037   0.041   0.037   0.037
14   [S.sub.max]   949.0   924.0   951.3   954.9
     b             0.150   0.157   0.150   0.149
     a             0.912   0.920   0.912   0.911
     [O.sub.R]     0.079   0.079   0.079   0.079
15   [S.sub.max]   929.1   959.3   931.4   933.6
     b             0.405   0.377   0.403   0.402
     a             0.919   0.911   0.918   0.918
     [O.sub.R]     0.040   0.040   0.040   0.040
16   [S.sub.max]   503.8   294.5   497.5   523.3
     b             0.037   0.072   0.037   0.035
     a             0.870   0.917   0.869   0.860
     [O.sub.R]     0.065   0.079   0.065   0.065
17   [S.sub.max]   580.9   250.3   635.9   541.8
     b             0.010   0.025   0.009   0.011
     a             0.866   0.936   0.865   0.869
     [O.sub.R]     0.063   0.079   0.063   0.063
18   [S.sub.max]   1296    1111    1287    1277
     b             0.601   0.796   0.609   0.623
     a             0.765   0.797   0.766   0.770
     [O.sub.R]     0.032   0.035   0.032   0.032
19   [S.sub.max]   181.3   177.3   181.2   180.5
     b             1.055   1.122   1.057   1.069
     a             0.660   0.684   0.660   0.664
     [O.sub.R]     0.024   0.024   0.024   0.024
20   [S.sub.max]   244.5   201.6   247.8   267.0
     b             0.205   0.255   0.202   0.185
     a             0.662   0.763   0.657   0.632
     [O.sub.R]     0.015   0.019   0.015   0.016
21   [S.sub.max]   267.6   222.6   265.7   252.1
     b             0.499   0.701   0.502   0.553
     a             0.539   0.641   0.540   0.562
     [O.sub.R]     0.011   0.016   0.011   0.011

                             (b)

#     PSO algorithm                  PCPSO algorithm

      Mean     Best     Std          Mean     Best     Std

1     0.2008   0.1282   1.13E - 01   0.1333   0.1283   7.13E - 03
2     0.2156   0.0955   1.03E - 01   0.1613   0.0940   5.93E - 02
3     0.1168   0.0247   1.11E - 01   0.0686   0.0239   7.20E - 02
4     0.3748   0.1252   1.86E - 01   0.2916   0.1386   9.92E - 02
5     0.0890   0.0308   1.00E - 01   0.0444   0.0316   1.13E - 02
6     0.0852   0.0742   1.92E - 02   0.0819   0.0742   2.20E - 02
7     0.1303   0.0457   1.56E - 01   0.0509   0.0458   7.64E - 03
8     0.1234   0.0824   9.87E - 02   0.1109   0.0827   1.79E - 02
9     0.3188   0.2591   4.23E - 02   0.2923   0.2689   2.44E - 02
10    0.1800   0.1010   8.93E - 02   0.1070   0.1010   1.63E - 02
11    0.1699   0.0704   1.27E - 01   0.0873   0.0724   1.49E - 02
12    0.1354   0.0309   1.38E - 01   0.0598   0.0310   4.64E - 02
13    0.1514   0.0367   1.38E - 01   0.0674   0.0369   4.01E - 02
14    0.2191   0.0788   1.83E - 01   0.1939   0.0788   1.25E - 01
15    0.1535   0.0396   1.65E - 01   0.0491   0.0396   1.16E - 02
16    0.1775   0.0651   1.64E - 01   0.1219   0.0829   5.20E - 02
17    0.2818   0.0726   2.34E - 01   0.1616   0.0834   4.53E - 02
18    0.0578   0.0320   3.03E - 02   0.0483   0.0320   3.19E - 02
19    0.0514   0.0235   3.96E - 02   0.0360   0.0235   2.47E - 02
20    0.0466   0.0154   5.14E - 02   0.0376   0.0156   5.14E - 02
      0.0428   0.0110   6.00E - 02   0.0282   0.0113   3.27E - 02

#     ABC algorithm                  MABC algorithm

      Mean     Best     Std          Mean     Best     Std

1     0.1288   0.1282   7.48E - 04   0.1282   0.1282   1.32E - 05
2     0.0975   0.0943   4.08E - 03   0.0928   0.0924   1.02E - 03
3     0.0329   0.0237   9.39E - 03   0.0211   0.0201   1.16E - 03
4     0.1214   0.1063   1.42E - 02   0.0796   0.0771   2.08E - 03
5     0.0358   0.0309   1.92E - 03   0.0315   0.0308   6.33E - 04
6     0.0743   0.0742   8.81E - 05   0.0742   0.0742   8.01E - 05
7     0.0477   0.0460   1.75E - 03   0.0459   0.0457   2.59E - 04
8     0.0864   0.0835   3.88E - 03   0.0819   0.0817   1.67E - 04
9     0.2681   0.2607   5.12E - 03   0.2500   0.2490   1.78E - 03
10    0.1014   0.1010   2.72E - 04   0.1011   0.1010   7.26E - 05
11    0.0791   0.0715   7.33E - 03   0.0706   0.0704   1.64E - 04
12    0.0331   0.0312   2.61E - 03   0.0310   0.0309   3.78E - 05
13    0.0410   0.0373   3.96E - 03   0.0367   0.0367   5.56E - 05
14    0.0811   0.0791   2.73E - 03   0.0788   0.0788   3.11E - 06
15    0.0405   0.0396   1.33E - 03   0.0396   0.0396   3.72E - 05
16    0.0756   0.0660   8.24E - 03   0.0652   0.0651   6.61E - 05
17    0.0971   0.0731   2.31E - 02   0.0638   0.0625   1.22E - 03
18    0.0344   0.0323   1.65E - 03   0.0321   0.0320   9.92E - 05
19    0.0245   0.0236   1.09E - 03   0.0236   0.0235   3.75E - 05
20    0.0175   0.0159   1.83E - 03   0.0155   0.0154   9.48E - 05
      0.0140   0.0113   2.62E - 03   0.0117   0.0110   8.64E - 04

Table 6: (a) Best and worst solutions (set of model parameters)
achieved by different solvers in 20 independent runs for fitting the
experimental data on two-site Langmuir isotherm model. (b) Performance
of different solvers for fitting two-site Langmuir isotherm to the
experimental data.

                                 (a)

#     Parameters                      PSO             PCPSO

                                 Best    Worst   Best    Worst

1     [mathematical expression   78.86   502.3   87.57   103.1
      not reproducible]
      [b.sub.1]                  9.983   2.119   8.895   2.838
      [mathematical expression   944.4   9.705   956.9   785.6
      not reproducible]
      [b.sub.2]                  0.380   4.070   0.355   0.701
      [O.sub.R]                  0.127   0.195   0.128   0.134
2     [mathematical expression   737.2   109.5   484.1   1679
      not reproducible]
      [b.sub.1]                  0.031   6.837   0.056   0.433
      [mathematical expression   0.434   22.23   3.097   8.473
      not reproducible]
      [b.sub.2]                  3.091   8.363   3.898   4.842
      [O.sub.R]                  0.094   0.412   0.111   0.302
3     [mathematical expression   324.4   92.02   295.0   149.6
      not reproducible]
      [b.sub.1]                  0.016   0.558   0.039   0.091
      [mathematical expression   42.26   17.06   13.12   12.81
      not reproducible]
      [b.sub.2]                  0.722   5.077   5.458   6.979
      [O.sub.R]                  0.036   0.245   0.045   0.114
4     [mathematical expression   953.9   124.5   3152    150.3
      not reproducible]
      [b.sub.1]                  0.035   0.999   0.009   0.153
      [mathematical expression   0.021   0.151   2.076   128.6
      not reproducible]
      [b.sub.2]                  1.318   6.773   1.065   0.144
      [O.sub.R]                  0.126   0.521   0.083   0.283
5     [mathematical expression   1701    344.3   1697    1257
      not reproducible]
      [b.sub.1]                  0.161   0.545   0.161   0.290
      [mathematical expression   27.52   131.3   30.87   5.442
      not reproducible]
      [b.sub.2]                  9.245   1.002   7.133   6.675
      [O.sub.R]                  0.029   0.249   0.029   0.046
6     [mathematical expression   1327    938.6   1320    816.3
      not reproducible]
      [b.sub.1]                  1.319   1.722   1.331   2.105
      [mathematical expression   49.07   164.2   48.82   565.2
      not reproducible]
      [b.sub.2]                  9.972   4.804   9.851   0.763
      [O.sub.R]                  0.074   0.084   0.074   0.075
7     [mathematical expression   1468    564.7   1307    1086
      not reproducible]
      [b.sub.1]                  0.171   1.296   0.237   0.400
      [mathematical expression   122.2   0.0002  66.06   21.15
      not reproducible]
      [b.sub.2]                  2.649   5.624   4.290   5.282
      [O.sub.R]                  0.042   0.145   0.045   0.051
8     [mathematical expression   1812    213.5   1627    3296
      not reproducible]
      [b.sub.1]                  0.064   1.652   0.074   0.031
8     [mathematical expression   22.14   8.328   20.94   13.98
      not reproducible]
      [b.sub.2]                  9.999   3.118   7.233   7.006
      [O.sub.R]                  0.077   0.3447  0.078   0.109
9     [mathematical expression   3504    480.8   7609    1415
      not reproducible]
      [b.sub.1]                  0.049   1.503   0.065   0.366
      [mathematical expression   2823    5.276   1.465   1.004
      not reproducible]
      [b.sub.2]                  0.116   3.801   6.330   5.246
      [O.sub.R]                  0.253   0.347   0.252   0.279
10    [mathematical expression   2299    1452    1243    2318
      not reproducible]
      [b.sub.1]                  1.399   2.151   1.399   1.385
      [mathematical expression   31.51   151.7   1088    7.399
      not reproducible]
      [b.sub.2]                  1.399   2.949   1.399   6.203
      [O.sub.R]                  0.101   0.108   0.101   0.101
11    [mathematical expression   29.81   15.23   25.70   8.095
      not reproducible]
      [b.sub.1]                  0.053   0.823   0.074   0.334
      [mathematical expression   2.836   2.167   2.609   1.953
      not reproducible]
      [b.sub.2]                  5.293   0.019   6.006   3.362
      [O.sub.R]                  0.075   0.182   0.078   8.095
12    [mathematical expression   37.94   21.89   36.58   19.27
      not reproducible]
      [b.sub.1]                  0.276   1.282   0.215   0.607
      [mathematical expression   3.505   20.95   6.099   7.660
      not reproducible]
      [b.sub.2]                  8.607   0.072   4.588   3.475
      [O.sub.R]                  0.034   0.058   0.037   0.102
13    [mathematical expression   664.7   128.3   608.9   1004
      not reproducible]
      [b.sub.1]                  0.053   2.270   0.059   0.041
      [mathematical expression   66.76   13.99   79.34   28.35
      not reproducible]
      [b.sub.2]                  1.784   6.641   1.375   7.985

      [O.sub.R]                  0.031   0.375   0.034   0.091
14    [mathematical expression   922.4   222.5   804.5   672.6
      not reproducible]
      [b.sub.1]                  0.119   1.379   0.161   0.164
      [mathematical expression   21.01   8.495   10.72   59.10
      not reproducible]
      [b.sub.2]                  4.916   3.209   7.793   1.574
      [O.sub.R]                  0.072   0.309   0.075   0.086
15    [mathematical expression   897.9   265.4   877.5   220.8
      not reproducible]
      [b.sub.1]                  0.333   2.947   0.355   1.430
      [mathematical expression   37.85   1.458   34.51   666.4
      not reproducible]
      [b.sub.2]                  6.720   4.238   6.991   0.290

      [O.sub.R]                  0.039   0.232   0.039   0.046
16    [mathematical expression   383.3   31.57   429.5   133.8
      not reproducible]
      [b.sub.1]                  0.038   0.876   0.034   0.200
      [mathematical expression   4.521   1.361   2.832   1.973
      not reproducible]
      [b.sub.2]                  3.485   2.994   6.396   5.194
      [O.sub.R]                  0.070   0.366   0.072   0.167
17    [mathematical expression   246.9   4.863   166.9   8.144
      not reproducible]
      [b.sub.1]                  0.019   7.005   0.034   5.133
      [mathematical expression   0.747   2.472   1.131   0.611
      not reproducible]
      [b.sub.2]                  7.408   8.274   4.310   6.367

      [O.sub.R]                  0.063   0.705   0.097   0.666
18    [mathematical expression   1532    670.5   1161    706.8
      not reproducible]
      [b.sub.1]                  0.171   2.061   0.277   1.726
      [mathematical expression   205.3   88.30   198.0   157.2
      not reproducible]
      [b.sub.2]                  9.999   6.937   9.995   4.266
      [O.sub.R]                  0.041   0.088   0.043   0.089
19    [mathematical expression   165.1   77.05   110.5   135.2
      not reproducible]
      [b.sub.1]                  0.023   2.275   3.815   2.047
      [mathematical expression   117.8   72.91   80.92   16.50
      not reproducible]
      [b.sub.2]                  3.519   2.438   0.111   5.067

      [O.sub.R]                  0.010   0.039   0.015   0.038
20    [mathematical expression   167.6   84.16   165.5   152.5
      not reproducible]
      [b.sub.1]                  0.135   2.507   0.135   0.236
      [mathematical expression   24.21   12.32   24.87   10.73
      not reproducible]
      [b.sub.2]                  7.752   6.932   8.009   1.760
      [O.sub.R]                  0.014   0.157   0.015   0.041
21    [mathematical expression   144.2   105.3   144.4   117.0
      not reproducible]
      [b.sub.1]                  0.087   6.227   0.184   3.267
      [mathematical expression   103.3   29.01   77.79   29.79
      not reproducible]
      [b.sub.2]                  3.563   4.744   5.675   4.661

      [O.sub.R]                  0.004   0.122   0.011   0.089

                                 (a)

#     Parameters                       ABC           MABC

                                 Best    Worst   Best    Worst

1     [mathematical expression   83.03   146.9   80.01   82.78
      not reproducible]
      [b.sub.1]                  9.669   4.322   10.00   8.945
      [mathematical expression   942.9   872.8   969.4   959.4
      not reproducible]
      [b.sub.2]                  0.369   0.357   0.360   0.367
      [O.sub.R]                  0.127   0.130   0.127   0.128
2     [mathematical expression   353.8   648.5   685.9   694.3
      not reproducible]
      [b.sub.1]                  0.031   0.038   0.035   0.034
      [mathematical expression   334.4   2.301   2.541   0.725
      not reproducible]
      [b.sub.2]                  0.039   0.370   0.036   0.121
      [O.sub.R]                  0.092   0.094   0.092   0.092
3     [mathematical expression   347.3   269.7   387.9   318.4
      not reproducible]
      [b.sub.1]                  0.021   0.018   0.012   0.024
      [mathematical expression   24.72   59.81   57.98   24.70
      not reproducible]
      [b.sub.2]                  1.470   0.257   0.322   1.026
      [O.sub.R]                  0.026   0.034   0.018   0.026
4     [mathematical expression   2934    860.1   5356    2199
      not reproducible]
      [b.sub.1]                  0.010   0.025   0.005   0.014
      [mathematical expression   5.534   690.5   0.012   40.61
      not reproducible]
      [b.sub.2]                  0.040   0.020   1.435   0.026
      [O.sub.R]                  0.079   0.106   0.076   0.084
5     [mathematical expression   1824    1615    1839    1991
      not reproducible]
      [b.sub.1]                  0.138   0.157   0.140   0.118
      [mathematical expression   37.83   90.92   31.16   46.78
      not reproducible]
      [b.sub.2]                  7.241   1.663   9.996   5.910
      [O.sub.R]                  0.029   0.035   0.028   0.030
6     [mathematical expression   1270    1257    1352    1208
      not reproducible]
      [b.sub.1]                  1.335   1.560   1.248   1.538
      [mathematical expression   81.96   33.59   58.03   94.98
      not reproducible]
      [b.sub.2]                  5.946   6.702   9.901   2.984
      [O.sub.R]                  0.075   0.075   0.074   0.075
7     [mathematical expression   1218    1017    2029    1327
      not reproducible]
      [b.sub.1]                  0.159   0.242   0.076   0.131
      [mathematical expression   253.9   287.8   264.7   244.5
      not reproducible]
      [b.sub.2]                  1.345   0.939   1.516   1.673
      [O.sub.R]                  0.042   0.048   0.039   0.043
8     [mathematical expression   1495    1168    1799    1490
      not reproducible]
      [b.sub.1]                  0.080   0.141   0.064   0.081
8     [mathematical expression   28.30   2.741   22.07   29.66
      not reproducible]
      [b.sub.2]                  4.511   6.117   9.820   3.233
      [O.sub.R]                  0.080   0.091   0.077   0.081
9     [mathematical expression   5915    1985    13067   7344
      not reproducible]
      [b.sub.1]                  0.080   0.266   0.036   0.066
      [mathematical expression   15.16   1.310   1114    4.637
      not reproducible]
      [b.sub.2]                  0.826   0.126   0.010   2.912
      [O.sub.R]                  0.253   0.266   0.247   0.252
10    [mathematical expression   2270    2398    2315    1717
      not reproducible]
      [b.sub.1]                  1.442   1.362   1.404   1.496
      [mathematical expression   9.769   4.237   12.33   559.0
      not reproducible]
      [b.sub.2]                  0.729   3.670   1.289   1.274
      [O.sub.R]                  0.101   0.101   0.101   0.101
11    [mathematical expression   28.04   20.63   30.08   26.57
      not reproducible]
      [b.sub.1]                  0.061   0.127   0.051   0.060
      [mathematical expression   2.696   2.222   2.908   2.632
      not reproducible]
      [b.sub.2]                  5.611   5.403   5.195   6.527
      [O.sub.R]                  0.075   0.097   0.075   0.079
12    [mathematical expression   39.82   40.01   38.18   36.05
      not reproducible]
      [b.sub.1]                  0.236   0.150   0.265   0.199
      [mathematical expression   3.828   7.602   3.610   6.909
      not reproducible]
      [b.sub.2]                  8.296   3.785   8.471   4.168
      [O.sub.R]                  0.035   0.041   0.034   0.038
13    [mathematical expression   647.9   561.3   660.9   613.6
      not reproducible]
      [b.sub.1]                  0.068   0.111   0.052   0.088
      [mathematical expression   53.71   27.61   71.29   34.55
      not reproducible]
      [b.sub.2]                  2.114   4.196   1.652   3.133
      [O.sub.R]                  0.034   0.042   0.031   0.037
14    [mathematical expression   876.0   907.8   919.8   821.1
      not reproducible]
      [b.sub.1]                  0.136   0.115   0.120   0.139
      [mathematical expression   17.56   39.02   19.75   14.65
      not reproducible]
      [b.sub.2]                  5.244   3.319   5.278   6.665
      [O.sub.R]                  0.072   0.080   0.072   0.075
15    [mathematical expression   860.1   993.8   895.7   797.8
      not reproducible]
      [b.sub.1]                  0.323   0.216   0.325   0.375
      [mathematical expression   49.64   104.8   41.51   73.83
      not reproducible]
      [b.sub.2]                  5.345   3.371   6.266   3.133
      [O.sub.R]                  0.040   0.045   0.039   0.043
16    [mathematical expression   382.3   240.8   653.4   515.8
      not reproducible]
      [b.sub.1]                  0.029   0.069   0.011   0.014
      [mathematical expression   11.02   5.644   34.18   21.70
      not reproducible]
      [b.sub.2]                  1.587   2.380   0.595   1.060
      [O.sub.R]                  0.071   0.089   0.069   0.074
17    [mathematical expression   248.7   106.9   329.6   100.6
      not reproducible]
      [b.sub.1]                  0.018   0.045   0.014   0.069
      [mathematical expression   2.265   2.451   0.701   47.52
      not reproducible]
      [b.sub.2]                  1.709   1.468   8.615   0.018
      [O.sub.R]                  0.075   0.127   0.060   0.119
18    [mathematical expression   1117    990.5   1411    1397
      not reproducible]
      [b.sub.1]                  0.259   0.241   0.200   0.197
      [mathematical expression   216.2   293.5   202.7   207.3
      not reproducible]
      [b.sub.2]                  9.364   7.183   9.997   9.775
      [O.sub.R]                  0.046   0.062   0.041   0.042
19    [mathematical expression   114.4   114.2   127.5   98.13
      not reproducible]
      [b.sub.1]                  3.763   0.687   0.034   4.354
      [mathematical expression   79.45   51.84   117.4   87.77
      not reproducible]
      [b.sub.2]                  0.088   9.713   3.506   0.170
      [O.sub.R]                  0.015   0.027   0.010   0.019
20    [mathematical expression   153.0   103.2   166.4   145.1
      not reproducible]
      [b.sub.1]                  0.125   0.133   0.131   0.093
      [mathematical expression   36.22   65.74   26.07   55.86
      not reproducible]
      [b.sub.2]                  2.292   0.750   5.570   1.244
      [O.sub.R]                  0.015   0.027   0.014   0.019
21    [mathematical expression   137.0   140.39  147.4   139.7
      not reproducible]
      [b.sub.1]                  0.117   0.269   0.094   0.283
      [mathematical expression   96.69   68.81   99.24   66.28
      not reproducible]
      [b.sub.2]                  3.978   6.682   3.827   7.002
      [O.sub.R]                  0.006   0.014   0.004   0.015

                                 (b)

      PSO algorithm                   PCPSO algorithm

      Mean     Best     Std           Mean     Best     Std

1     0.1375   0.1273   2.24E - 02    0.1294   0.1275   1.29E - 03
2     0.1870   0.0959   1.15E - 01    0.1406   0.0996   3.41E - 02
3     0.1562   0.0358   1.13E - 01    0.0569   0.0346   1.76E - 02
4     0.2318   0.1334   1.14E - 01    0.1962   0.1191   7.12E - 02
5     0.0670   0.0277   7.99E - 02    0.0376   0.0294   4.21E - 03
6     0.0775   0.0743   9.28E - 03    0.0746   0.0744   1.19E - 04
7     0.0432   0.0381   3.46E - 03    0.0471   0.0421   3.10E - 03
8     0.1035   0.0776   4.41E - 02    0.0883   0.0783   6.50E - 03
9     0.2751   0.2473   2.75E - 02    0.2642   0.2489   8.20E - 03
10    0.1010   0.1010   7.04E - 05    0.1010   0.1010   1.34E - 05
11    0.0888   0.0748   1.47E - 02    0.1093   0.0831   1.65E - 02
12    0.0513   0.0341   2.80E - 02    0.0536   0.0353   1.82E - 02
13    0.0776   0.0322   9.46E - 02    0.0535   0.0318   1.23E - 02
14    0.1240   0.0719   1.13E - 01    0.0760   0.0721   3.89E - 03
15    0.0487   0.0393   1.65E - 02    0.0440   0.0396   5.70E - 03
16    0.1185   0.0706   1.15E - 01    0.1156   0.0695   2.93E - 02
17    0.1318   0.0791   3.39E - 02    0.1418   0.0686   5.39E - 02
18    0.0629   0.0419   1.85E - 02    0.0557   0.0423   9.23E - 03
19    0.0251   0.0140   7.53E - 03    0.0290   0.0216   4.60E - 03
20    0.0162   0.0143   3.31E - 03    0.0270   0.0152   2.30E - 02
21    0.0220   0.0030   2.95E - 02    0.0140   0.0100   1.86E - 03

      ABC algorithm                  MABC algorithm

      Mean     Best     Std          Mean      Best     Std

1     0.1291   0.1278   1.17E - 03   0.1278    0.1273   7.05E - 04
2     0.1017   0.0938   7.20E - 03   0.0954    0.0924   4.59E - 03
3     0.0408   0.0323   8.71E - 03   0.0271    0.0200   3.62E - 03
4     0.1296   0.1011   2.56E - 02   0.0820    0.0772   5.01E - 03
5     0.0357   0.0315   2.24E - 03   0.0287    0.0277   8.36E - 04
6     0.0747   0.0744   1.12E - 04   0.0746    0.0743   1.03E - 04
7     0.0438   0.0407   2.63E - 03   0.0418    0.0385   2.60E - 03
8     0.0843   0.0786   4.03E - 03   0.0777    0.0772   6.28E - 04
9     0.2578   0.2523   2.53E - 01   0.2504    0.2475   3.90E - 03
10    0.1010   0.1010   3.49E - 05   0.1010    0.1010   2.34E - 05
11    0.0802   0.0751   5.69E - 03   0.0751    0.0747   3.90E - 04
12    0.0402   0.0344   4.30E - 03   0.0343    0.0341   1.77E - 04
13    0.0379   0.0332   2.40E - 03   0.0351    0.0316   2.20E - 03
14    0.0746   0.0724   1.73E - 03   0.0721    0.0719   1.12E - 04
15    0.0495   0.0398   2.61E - 03   0.0396    0.0393   3.43E - 04
16    0.0838   0.0712   1.30E - 02   0.0707    0.0688   1.50E - 03
17    0.1091   0.0876   1.35E - 02   0.0680    0.0596   1.05E - 02
18    0.0524   0.0435   4.49E - 03   0.0456    0.0423   2.45E - 03
19    0.0190   0.0129   4.32E - 03   0.0160    0.0111   3.51E - 03
20    0.0173   0.0149   1.76E - 03   0.0148    0.0143   5.62E - 04
21    0.0136   0.0097   2.93E - 03   0.0079    0.0030   2.41E - 03

Table 7: Performance of different solvers for the selection of
Appropriate algorithm based on mean solution.

#             PSO                  PCPSO

1      Two-site Langmuir     Two-site Langmuir
           (0.1375)              (0.1294)
2      Two-site Langmuir        Freundlich
           (0.1870)              (0.1119)
3     Freundlich-Langmuir       Freundlich
           (0.1168)              (0.0311)
4         Freundlich            Freundlich
           (0.2045)              (0.0819)
5         Freundlich         Two-site Langmuir
           (0.0614)              (0.0376)
6      Two-site Langmuir     Two-site Langmuir
           (0.0775)              (0.0746)
7      Two-site Langmuir     Two-site Langmuir
           (0.0432)              (0.0471)
8      Two-site Langmuir        Freundlich
           (0.1035)              (0.0853)
9         Freundlich            Freundlich
           (0.1537)              (0.1330)
10     Two-site Langmuir     Two-site Langmuir
           (0.1010)              (0.1010)
11     Two-site Langmuir    Freundlich-Langmuir
           (0.0888)              (0.0873)
12     Two-site Langmuir     Two-site Langmuir
           (0.0513)              (0.0536)
13     Two-site Langmuir     Two-site Langmuir
           (0.0776)              (0.0535)
14     Two-site Langmuir     Two-site Langmuir
           (0.1240)              (0.0760)
15     Two-site Langmuir     Two-site Langmuir
           (0.0487)              (0.0440)
16        Freundlich         Two-site Langmuir
           (0.1155)              (0.1156)
17     Two-site Langmuir        Freundlich
           (0.1318)              (0.1019)
18    Freundlich-Langmuir   Freundlich-Langmuir
           (0.0578)              (0.0483)
19     Two-site Langmuir     Two-site Langmuir
           (0.0251)              (0.0290)
20     Two-site Langmuir     Two-site Langmuir
           (0.0162)              (0.0270)
21     Two-site Langmuir     Two-site Langmuir
           (0.0220)              (0.0140)

#              ABC                    MABC

1      Freundlich-Langmuir     Freundlich-Langmuir
            (0.1288)                (0.1282)
2      Freundlich-Langmuir     Freundlich-Langmuir
            (0.0975)                (0.0928)
3          Freundlich          Freundlich-Langmuir
            (0.0233)                (0.0211)
4          Freundlich          Freundlich-Langmuir
            (0.0827)                (0.0796)
5       Two-site Langmuir       Two-site Langmuir
            (0.0357)                (0.0290)
6      Freundlich-Langmuir     Freundlich-Langmuir
            (0.0743)                (0.0742)
7       Two-site Langmuir       Two-site Langmuir
            (0.0438)                (0.0418)
8       Two-site Langmuir       Two-site Langmuir
            (0.0843)                (0.0777)
9          Freundlich              Freundlich
            (0.1330)                (0.1330)
10      Two-site Langmuir      Freundlich/Two-site
                                    Langmuir
            (0.1010)                (0.1010)
11     Freundlich-Langmuir     Freundlich-Langmuir
            (0.0791)                (0.0706)
12     Freundlich-Langmuir     Freundlich-Langmuir
            (0.0331)                (0.0310)
13      Two-site Langmuir       Two-site Langmuir
            (0.0379)                (0.0351)
14      Two-site Langmuir       Two-site Langmuir
            (0.0746)                (0.0721)
                                   Freundlich-
15     Freundlich-Langmuir      Langmuir/Two-site
                                    Langmuir
            (0.0405)                (0.0396)
16     Freundlich-Langmuir     Freundlich-Langmuir
            (0.0756)                (0.00652)
17         Freundlich          Freundlich-Langmuir
            (0.0677)                (0.0638)
18     Freundlich-Langmuir     Freundlich-Langmuir
            (0.0344)                (0.0321)
19      Two-site Langmuir       Two-site Langmuir
            (0.0190)                (0.0160)
20      Two-site Langmuir       Two-site Langmuir
            (0.0173)                (0.0148)
21      Two-site Langmuir       Two-site Langmuir
            (0.0136)                (0.0079)

Table 8: Performance of different solvers for the selection of
appropriate algorithm based on best solution.

#            PSO                    PCPSO

1     Two-site Langmuir       Two-site Langmuir
          (0.1273)                 (0.1275)

2    Freundlich-Langmuir     Freundlich-Langmuir
          (0.0955)                 (0.0940)
3        Freundlich               Freundlich
          (0.0233)                 (0.0233)
4        Freundlich               Freundlich
          (0.0798)                 (0.0798)
5     Two-site Langmuir       Two-site Langmuir
          (0.0277)                 (0.0294)
6    Freundlich-Langmuir     Freundlich-Langmuir
          (0.0742)                 (0.0742)
7     Two-site Langmuir       Two-site Langmuir
          (0.0381)                 (0.0421)
8     Two-site Langmuir       Two-site Langmuir
          (0.0776)                 (0.0783)
9        Freundlich               Freundlich
          (0.1330)                 (0.1330)
         Freundlich-             Freundlich-
10    Langmuir/Two-site       Langmuir/two-site
          Langmuir                 Langmuir
          (0.1010)                 (0.1010)
11   Freundlich-Langmuir     Freundlich-Langmuir
          (0.0704)                 (0.0724)
12   Freundlich-Langmuir     Freundlich-Langmuir
          (0.0309)                 (0.0310)
13    Two-site Langmuir       Two-site Langmuir
          (0.0322)                 (0.0318)
14    Two-site Langmuir       Two-site Langmuir
          (0.0719)                 (0.0721)
15    Two-site Langmuir    Freundlich-Langmuir/two-
                                site Langmuir
          (0.0393)                 (0.0396)
16   Freundlich-Langmuir      Two-site Langmuir
          (0.0651)                 (0.0695)
17       Freundlich               Freundlich
          (0.0676)                 (0.0676)
18   Freundlich-Langmuir     Freundlich-Langmuir
          (0.0320)                 (0.0320)
19    Two-site Langmuir      Freundlich-Langmuir
          (0.0140)                 (0.0235)
20    Two-site Langmuir       Two-site Langmuir
          (0.0143)                 (0.0152)
21   Freundlich-Langmuir      Two-site Langmuir
          (0.0110)                 (0.0100)

#            ABC                    MABC

1     Two-site Langmuir      Two-site Langmuir
          (0.1278)                (0.1273)
                                Freundlich-
2     Two-site Langmuir      Langmuir/two-site
                                  Langmuir
          (0.0938)                (0.0924)
3        Freundlich          Two-site Langmuir
          (0.0233)                (0.0200)
4        Freundlich         Freundlich-Langmuir
          (0.0798)                (0.0771)
5    Freundlich-Langmuir     Two-site Langmuir
          (0.0309)                (0.0279)
6    Freundlich-Langmuir    Freundlich-Langmuir
          (0.0742)                (0.0742)
7     Two-site Langmuir      Two-site Langmuir
          (0.0407)                (0.0390)
8     Two-site Langmuir      Two-site Langmuir
          (0.0786)                (0.0772)
9        Freundlich              Freundlich
          (0.1330)                (0.1330)
         Freundlich-       Freundlich/Freundlich-
10    Langmuir/Two-site      Langmuir/two-site
          Langmuir                Langmuir
          (0.1010)                (0.1010)
11   Freundlich-Langmuir    Freundlich-Langmuir
          (0.0715)                (0.0704)
12   Freundlich-Langmuir    Freundlich-Langmuir
          (0.0312)                (0.0309)
13    Two-site Langmuir      Two-site Langmuir
          (0.0332)                (0.0316)
14    Two-site Langmuir      Two-site Langmuir
          (0.0724)                (0.0719)
15   Freundlich-Langmuir     Two-site Langmuir
          (0.0396)                (0.0393)
16   Freundlich-Langmuir    Freundlich-Langmuir
          (0.0660)                (0.0651)
17       Freundlich          Two-site Langmuir
          (0.0677)                (0.0596)
18   Freundlich-Langmuir    Freundlich-Langmuir
          (0.0323)                (0.0320)
19    Two-site Langmuir      Two-site Langmuir
          (0.0129)                (0.0111)
20    Two-site Langmuir      Two-site Langmuir
          (0.0149)                (0.0143)
21    Two-site Langmuir      Two-site Langmuir
          (0.0097)                (0.0030)
COPYRIGHT 2017 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2017 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Bharat, Tadikonda Venkata
Publication:Advances in Artificial Intelligence
Article Type:Report
Date:Jan 1, 2017
Words:17521
Previous Article:iWordNet: A New Approach to Cognitive Science and Artificial Intelligence.
Next Article:Natural Language Processing and Fuzzy Tools for Business Processes in a Geolocation Context.
Topics:

Terms of use | Privacy policy | Copyright © 2019 Farlex, Inc. | Feedback | For webmasters