Research Article  Open Access
Study of Isothermal, Kinetic, and Thermodynamic Parameters for Adsorption of Cadmium: An Overview of Linear and Nonlinear Approach and Error Analysis
Abstract
Reports about presence and toxicity of Cd^{2+} in different chemical industrial effluents prompted the researchers to explore some economical, rapid, sensitive, and accurate methods for its determination and removal from aqueous systems. In continuation of series of investigations, adsorption of Cd^{2+} onto the stem of Saccharum arundinaceum is proposed in the present work. Optimization of parameters affecting sorption potential of Cd^{2+} including pH, contact time, temperature, sorbent dose, and concentration of sorbate was carried out to determine best suited conditions for maximum removal of sorbate. To understand the nature of sorption process, linear and nonlinear forms of five sorption isotherms including Freundlich and Langmuir models were employed. Feasibility and viability of sorption process were evaluated by calculating kinetics and thermodynamics of the process, while error analysis suggested best fitted sorption model on sorption data. Thermodynamic studies demonstrated exothermic nature of reaction, while kinetic studies suggested pseudosecond order of reaction.
1. Introduction
Environmental pollution should be taken into special consideration because it is a very serious matter affecting every type of organism at every level. The most adversely affected environmental resource is water [1]. As water is an essential element for the survival of living beings, it is very necessary to keep it pure and clean [2]. Quality of drinking water is of prime importance for mankind because waterborne diseases can decimate population of the whole area. These diseases arise due to toxic release of chemicals from industrial zones [3]. Particularly in industrial areas, these waterborne diseases are a great threat towards safety of water supplies. Other sources which may pollute water include domestic waste, pesticides run off from agricultural land, metal plating operations, and so on. Key contaminants present in water include heavy metals, chlorinated hydrocarbons, pathogens, detergents, pesticides, algal nutrients, trace organic compounds, dyes, and so on. These hazardous substances are of concern because of their ultimate effect on survival of human life [4]. Heavy metals like Cd, Zn, Ni, and Pb are present in relatively major amounts in industrial effluents and enter in rivers and oceans and ultimately pollute groundwater leading to adverse effects on aquatic life. Metals resist the process of biodegradability and hence remain in ecosystem, affecting food chain and human health [5]. In order to provide clean environment and healthy lifestyle to our coming generations, it is necessary to remove hazardous pollutants from the environment. In environmental restoration areas, conventional techniques are practiced to eradicate those pollutants from the environment which include chemical precipitation, evaporative method, electrolytic extraction, reverse osmosis, ion exchange, and electrochemical and membrane processes [6]. All these methods are costly and produce large amount of sludge which is difficult to dispose of. Use of biological materials for the removal of pollutants from aqueous media is considered superior to other methods in terms of cost effectiveness and simple design. It is a surface phenomenon, in which pollutants get accumulated on the surface of the adsorbent material. Binding nature is based on type of sorbent and sorbate, but mostly physisorption or chemisorption takes place [7]. Materials with ease in availability and low cost are preferred for the purpose. In this context, agrowastes are considered a significant material for adsorption. Binding capacity of these materials can be intensified by physical and chemical treatments and heat therapy [8].
To explore the appropriate adsorbent, it is necessary to establish equilibrium correlation of sorbent to predict behavior of sorbent under different experimental conditions. This equilibrium correlation is developed by using equilibrium isotherms. These isotherms express way of sorbent interaction with the surface of adsorbent, that is, whether it is monolayer or multilayer sorption [9]. Similarly, thermodynamic studies are of prime importance to predict whether the adsorption is spontaneous or not. Furthermore, it provides information about suitable temperature range for sorption and nature of sorbent and sorbate at equilibrium [10].
The aim of the present research was to explore Saccharum arundinaceum for adsorption of cadmium under different operating conditions including pH, contact time, initial concentration, and temperature. The application of linear and nonlinear forms of equilibrium isotherms was to determine appropriate isotherm for the purpose, and thermodynamic and kinetic studies were performed to check reaction nature of the adsorption phenomenon. Error analysis based on five different error functions was also performed.
2. Materials and Methods
2.1. Preparation of Adsorbent
On the basis of literature survey and indigenous availability of agrowaste materials, the stem of Saccharum arundinaceum (hardy sugar cane) was collected from different regions of Sargodha District, Pakistan. After collection, the sample was properly washed with deionized water to remove dust and surface impurities. The sorbent was initially dried in an open container at room temperature and later in an electric oven (Model, LEB120) at 105°C for 24 h to remove all the moisture contents. The dried sorbent was ground, and appropriate particle size was separated by sieves and was stored for further analyses.
2.2. Chemicals
All the chemicals, reagents, and solvents used in the present work were of analytical reagent grade and purchased from Merck (Germany) or SigmaAldrich (Germany). Standard solutions were prepared, and successive dilutions were made with doubledistilled water to make working solutions.
2.3. Pretreatment of Sorbents
Saccharum arundinaceum was pretreated with HCl (0.1 M) and NaOH (0.1 M) to evaluate the effects of acid and base treatments on pore size, that is, pore area, pore volume, and sorption capacity. For chemical treatment, the sorbent (20 g) was stirred for 4 h in 1 L solution of 0.1 M NaOH or HCl followed by filtration and extensive washing with distilled water to remove any traces of acid/base. After that, the treated sorbent material was dried at 110°C and stored in airtight zipper bags at −4°C before further use.
2.4. Characterization of Sorbent
To determine different physical and chemical parameters affecting adsorption, it is necessary to characterize the sorbent. Therefore, physical and chemical characterization was done by scanning electron microscopy (SEM) and Fourier transform infrared spectroscopy (FTIR).
2.4.1. Scanning Electron Microscopy
Surface analysis was performed using scanning electron microscope JEOL model 2300. SEM provides information about surface area available for adsorption and morphology of sorbent [11]. Analysis of each sorbent was carried out in optimized conditions under argon atmosphere.
2.4.2. Fourier Transform Infrared Spectroscopy
Functional groups present in structure of sorbent were determined by Fourier transform infrared spectrophotometer (Model Shimadzu AIM8800). These functional groups are responsible for adsorption of sorbate on the surface of sorbent, and their detection helps in determining the nature of binding interactions between sorbate and sorbent surface [12]. Diffused reflectance infrared technique (DRIFT) was used for analysis taking KBr as a background reagent.
2.5. Equilibrium Isotherms
In order to study adsorption pathway and equilibrium relationship between sorbent and sorbate, it is necessary to design proper adsorption isotherms. Isotherms predict the appropriate parameters and behavior of sorbent towards different sorption systems [13]. In this context, linear and nonlinear models are utilized using Microsoft Excel®2007 (equilibrium isotherms applied on the present work are given in Table 1s of supplementary data).
2.6. Error Functions
In order to determine best fitting of linear or nonlinear models onto adsorption data, it is necessary to calculate the error function [14]. These error functions include sum square error, hybrid functional error, average relative error, sum of absolute error, nonlinear chisquare, and so on (calculated error functions and their equations are present in Table 2s of supplementary data).
2.7. Thermodynamic Investigations
Thermodynamic investigations are another important parameter of adsorption studies. For thermodynamic studies, the adsorption experiment was carried out at different temperature conditions and calculated parameters included enthalpy (), entropy (), and Gibbs free energy ().
For this purpose, (1)–(3) were appliedwhere is the natural gas constant and is the constant at equilibrium and is calculated aswhere is the concentration of sorbent at equilibrium condition.
2.8. Adsorption Kinetics
In batch adsorption process, kinetic studies provide information about optimum conditions, mechanism of sorption, and possible rate controlling step. For this purpose, linear and nonlinear form of pseudofirst and pseudosecondorder kinetics is applied on adsorption data [15]. In order to check the effect of contact time (10–70 min) on adsorption, initial concentration of 100 mg/L for cadmium was prepared and 100 ml of this sample was used for study. Sorbent (0.5 g) was added in this cadmium solution and applied for shaking at 150 rpm speed. After fixed interval of time, the sample was removed from flask and analyzed for cadmium concentration by atomic absorption spectrophotometer. The amount of cadmium adsorbed at different time intervals was calculated by employing the following formula:where is the amount of cadmium adsorbed at any time , and are initial and equilibrium concentrations, respectively. The volume of cadmium solution taken is represented by V(L), and is the amount of sorbent in g.
2.8.1. PseudoFirstOrder Kinetics
In order to calculate pseudofirstorder kinetics for adsorption system, following equations were used:where is the amount adsorbed at time , is the equilibrium amount, t is time in minutes, and is the rate constant.
2.8.2. PseudoSecondOrder Kinetics
For pseudosecondorder kinetics, linear and nonlinear forms were applied as follows:
3. Results and Discussion
3.1. Effect of Pretreatment
Pretreatment has promising effect on adsorption potential of Saccharum arundinaceum. Results reveal that basetreated (97.5%) sorbent shows good efficiency for cadmium sorption as compared to raw (91.15%) and acidtreated (57.6) sorbent as shown in Figure 1. Adsorption capacity depends upon functional groups present on the surface of sorbent and its microporous structure [16]. Increase in sorption capacity by base treatment can be attributed to hydroxyl groups created on the surface of adsorbent by base treatment or modification of cell wall components by base [17]. Decrease in adsorption after acid treatments was found as the binding sites available on the surface of biosorbent got destructed due to acid [18]. So, basetreated Saccharum arundinaceum was used for adsorption analysis.
3.2. Characterization of Sorbents
Saccharum arundinaceum was characterized in terms of surface morphology and functional group analysis by scanning electron microscopy and Fourier transform infrared spectroscopy.
3.2.1. Scanning Electron Microscopy
Three native and two treated sorbents (acid and basetreated Saccharum arundinaceum) were analyzed through scanning electron microscope to study surface morphology. Large pore size available on the surfaces of native and basetreated sorbent was responsible for enhanced adsorption on the surface of these agrowaste materials. Results for SEM analysis are given in Figure 2. Hollow cavities appear in the structure of raw adsorbent, which were responsible for binding of sorbate onto the surface of sorbent. Acid treatment decreases these cavities by deforming surface of the sorbent, so adsorption decreases after acid treatment because surface becomes smooth and thin adsorption layer is formed on the sorbent surface. Raw and basetreated sorbent surface is found rough and cylindrical, which make possible multilayer and thick adsorption on the sorbent surface as compared to smooth surface. Results obtained in the SEM micrograph are in good agreement with reported data [19].
(a)
(b)
(c)
3.2.2. Fourier Transform Infrared Spectroscopy
Fourier transform infrared spectrometer provides information about functional groups present on the surface of sorbent and makes possible attachment of sorbate [20]. FTIR spectra of sorbent obtained in the range of 4000–450 cm^{−1} wavenumber and major functional groups present in adsorbent are listed in Table 1. A broadband appears in the range of 3000–3700 cm^{−1} which was due to –OH stretching vibration of hydroxyl functional groups including hydrogen bonding because broadband of –OH group in this range is indication of hydrogen bonding present in the compound. This peak appears in raw and basetreated sorbent but disappears in case of acid treated due to reaction of –OH group with acid hydrogen. Stretching band of –CH appears in 2900–3000 cm^{−1} wavenumber range for all sorbents. The peak at 1750 cm^{−1} appears due to C=O and 1200 cm^{−1} due to CO functional group. In some cases, –CN also appears at 1049 cm^{−1} value. Vibration due to secondary amide appears at 1645 cm^{−1} [21]. For adsorption purpose, significant role is played by –OH group and heteroatoms to attach sorbate on the surface.

3.3. Adsorption Study
Adsorption study was performed by the batch adsorption method by varying different parameters including pH, contact time, and initial concentration of sorbate to find best suited conditions for the removal of cadmium from aqueous media.
3.3.1. Effect of Contact Time
Contact time was varied from 10 to 100 minutes under neutral conditions with constant amount of sorbent (1 g), initial concentration (60 ppm), and shaking speed (150 rpm), and results obtained are shown in Figure 1s (Supplementary material).
Maximum adsorption was achieved at 60 minutes time interval and no significant increase found by further increase in time. Initially, excess of vacant places are available on the surface of sorbent, and uptake of metal ions was more, so there was continuous increase in adsorption capacity by increasing time slot from zero to 60 minutes. But, further increase could not cause sufficient change in adsorption of metals as vacant spaces are already filled, and equilibrium is achieved [22].
3.3.2. Effect of pH
Initial pH of the adsorption system has significant role in adsorption of sorbate, as it affects the surface morphology of the sorbent and binding nature of sorbate. The range of pH selected was 2–10 with 1 g sorbent and 60 ppm initial concentration of sorbate. The result given in Figure 2s (Supplementary material) reveals the fact that adsorption capacity is quite low under acidic conditions. When pH is increased, it causes increase in adsorbed amount of sorbate on the surface of sorbent. At low pH value, metals have to compete with H^{+} ions for adsorption on sorbent surface since H^{+} ions are present in excess at that pH value. But when pH value is raised, it causes significant increase in adsorption due to attraction developed between negatively charged surfaces of sorbent by –OH groups and positively charged metal ions [23]. So, for cadmium, optimum pH range was found from 6 to 8; in this range, cadmium shows best adsorption behavior. When pH is further increased, there is decline in adsorption capacity due to formation of metal hydrides.
3.3.3. Effect of Initial Concentration
Initial concentration of sorbate is another important parameter, which affects the adsorption phenomenon. For this purpose, initial concentration of cadmium was varied in range of 10–100 ppm by keeping all other parameters constant (Figure 3s Supplementary material).
Rapid increase in adsorption capacity was observed initially for adsorption of cadmium on Saccharum arundinaceum as vacant spaces were available on the surface of sorbent. So, rise in concentration also raised adsorption of sorbate on available sites [24]. Adsorbent readily occupies these adsorption sites, and adsorption capacity has positive influence of concentration in this range. Further increase in concentration from 60 to 100 ppm has no significant effect on adsorption phenomenon. Surface of adsorbent becomes saturated with sorbate, and after establishment of equilibria, increase in concentration has no significant influence on adsorption phenomenon. Previous studies also report that accommodation for sorbate decreases as concentration is very high due to unavailability of resident sites [25].
3.3.4. Effect of Temperature
The effect of temperature on adsorption was studied by using temperature range 20, 30, 40, and 50°C with pH 6 (Figure 4s) (Supplementary material). Adsorption of cadmium onto Saccharum arundinaceum agrowaste was found to increase with increase in temperature. At high temperature, intraparticle diffusion increases and more adsorption sites are created which boost up adsorption phenomenon [26].
3.4. Equilibrium Isotherms
Adsorption system can be designed by adsorption isotherms commonly known as equilibrium isotherms which represent the amount of solute adsorbed per unit weight of sorbent [27]. These isotherms use equilibrium concentration of sorbent at constant temperature. In order to remove effluents from the system, particular design is optimized to generate proper correlation for experimental data which is called adsorption isotherm. Researches proposed many isotherms in this regard which are based on the adsorption system including Langmuir, Freundlich, Redlich–Peterson, Temkin, and Elovich [28, 29]. Sorption was carried out by employing linear as well as nonlinear adsorption models by varying initial concentration from 10 to 100 ppm.
3.4.1. Freundlich Isotherm
Freundlich adsorption isotherm was developed for the heterogeneous system, and it gives concept of multilayer adsorption on the surface of sorbent (Figure 5s Supplementary material).
Parameters calculated for Freundlich isotherm by employing its linear and nonlinear form are given in Table 4. Freundlich isotherm was obtained by plotting log Cad versus log Ce. K_{F} and n are constants obtained from intercept and slope, respectively. Freundlich adsorption capacity (K_{F}) is an indicator of a system, whether it is favorable for adsorption or not. Adsorption is considered promising if value of K_{F} is found in range of 1–20, and results reveal that in the present study, K_{F} was 9.5 and 12.2, respectively, for linear and nonlinear approaches of Freundlich adsorption isotherm. Similarly, adsorption intensity represented by n indicates fitness of model for adsorption purposes if value of n is above 1. Value of R^{2} obtained from the plot is significant (0.9446) representing good fitness of this model for adsorption of cadmium onto Saccharum arundinaceum.
3.4.2. Langmuir Isotherm
Langmuir adsorption isotherm is based on monolayer adsorption of metal ions on the surface of agrowastes, and energy of adsorption system is considered constant. In order to calculate Langmuir model, initial concentrations were changed from 10 to 100 ppm with 1 g sorbent amount and 1 h shaking time. Distribution of metal ions between liquid and solid surface was calculated by equations given in Table 1, employing linear and nonlinear forms of the model. Langmuir adsorption isotherm was obtained by plotting Ce/Cad versus Ce as shown in Figure 6s Supplementary material. R^{2} value obtained for plot was found satisfactory showing fitness of model on the adsorption experiment. Q_{o} represents metal ion uptake per unit mass of adsorbent (mg/g) and b is Langmuir constant [15]. A dimensionless constant R_{L} is calculated by using Langmuir constant, and initial concentration represents model fitness for a particular system. If value of R_{L} falls between 0 and 1, the system is considered appropriate for adsorption purpose and Table 2 shows results which are in this range. Furthermore, experimental data and predicted results obtained for the present work were found in close correlation with low value of residual sum of square (0.006) making this model applicable for the present work.

3.4.3. Dubinin–Radushkevich Isotherm
Dubinin–Radushkevich isotherm was designed as an empirical model for adsorption of vapors onto solid surface. It is successfully applied for adsorption of heterogeneous system including solid and liquid. This model is considered more general than Langmuir because in its derivation homogenous surface and constant sorption potential are not assumed [30]. The relationship given in Table 1s (supplementary data) was employed to relate ln Cad with ε^{2}, where ε^{2} is Polanyi potential which is based on temperature, natural gas constant, and equilibrium concentration as given in the following equation:
Slope of the plot gives value of and the intercept is q_{s}. The model showed good applicability on the adsorption system in nonlinear form with high value of R^{2}.
Dubinin–Radushkevich isotherm has found very promising applications for determination of nature of sorption, whether it is physical or chemical. For this purpose, obtained from the slope of the plot is used in the following equation:
The value of calculated for the present research was 0.764 and suggests physical nature of sorption. Because the value of below 8 kj/mol reflects physical sorption and 8–16 kj/mol reflects chemical sorption (Figure 7s Supplementary material).
3.4.4. Temkin Isotherm
Temkin adsorption isotherm discusses interaction of sorbent and sorbate, and the model is based on assumption that heat of adsorption will not remain constant. It decreases due to interaction between sorbent and sorbate during adsorption phenomenon [31]. Linear and nonlinear forms of Temkin model are given in Table 1s (Supplementary data). Equilibrium constant of binding K_{T} provides information about binding energy, and β expresses heat of adsorption for a particular adsorption experiment (Figure 8s Supplementary material). Linear form of Temkin model is found more suitable with high value of binding constant as given in Table 2. The model indicates the exothermic nature of adsorption reaction as B > 0 which is an indicator of heat release during the process [32].
3.4.5. Elovich Isotherm
According to Elovich model, mechanism of adsorption is based on chemical reactions which are responsible for adsorption. Plot of ln C_{ad}/C_{e} versus C_{ad} gives R^{2} value close to unity. K_{E} and Q_{m} are obtained from intercept and slope of plot, respectively. K_{E} shows initial sorption rate and Q_{m} is adsorption constant. Initial sorption rate obtained from linear form of Elovich model is quite high (35,100.411) as compared to nonlinear form (11.7891), so making linear form adequate to describe adsorption of cadmium onto Saccharum arundinaceum. Furthermore, R^{2} value (0.9033) for linear form is also high than nonlinear form (0.835) (Figure 9sSupplementary material).
3.5. Error Analysis for Equilibrium Isotherms
In order to check the fit of adsorption model to experimental data, error functions are used [33]. In the present work, six error functions were applied on linear and nonlinear form of data by minimizing the error function in a range of concentration used for analysis by employing solver addin with Microsoft Excel 2010. Results for optimization of equilibrium isotherms by error analysis are given in Table 3. For meaningful results, a comparison of each error function for linear and nonlinear forms was made (Figure 10s Supplementary material).

For linear form of adsorption isotherms, a comparison of error functions reflects that Langmuir, Freundlich, and Elovich isotherms have good correlation with experimental values for the present adsorption study. These isotherms give low values for most of error functions. Applicability of these models for removal of cadmium ions from aqueous media is also studied by many researchers [34, 35]. R^{2}: Temkin > Langmuir > Elovich > Freundlich > Dubinin–Radushkevich RSS: Temkin > Dubinin–Radushkevich > Elovich > Freundlich > Langmuir ARE: Dubinin–Radushkevich > Temkin > Freundlich > Langmuir > Elovich EABS: Dubinin–Radushkevich > Langmuir > Temkin > Freundlich > Elovich Chisquare (χ^{2}): Elovich > Dubinin–Radushkevich > Freundlich > Temkin > Langmuir
Similar study was carried out by employing nonlinear form of adsorption models, and results are summarized below. Nonlinear form of Temkin isotherm was not found suitable for adsorption of cadmium onto Saccharum arundinaceum agrowaste because of high value of error functions. Freundlich and Elovich isotherms have been proved to be suitable models for this study with low value for error functions. R^{2}: Langmuir > Freundlich > Dubinin–Radushkevich > Temkin > Elovich RSS: Temkin > Dubinin–Radushkevich > Freundlich > Langmuir > Elovich ARE: Temkin > Dubinin–Radushkevich > Langmuir > Elovich > Freundlich EABS: Temkin > Dubinin–Radushkevich > Elovich > Langmuir > Freundlich Chisquare (χ^{2}): Temkin > Freundlich > Dubinin–Radushkevich > Elovich > Langmuir
A comparison between linear and nonlinear approaches of each adsorption isotherm was also made to select the most appropriate form for adsorption study. Linear form of Freundlich adsorption isotherm was found superior over nonlinear form with small error functions in most of the cases. R^{2}: Freundlich (linear approach) < Freundlich (nonlinear) RSS: Freundlich (linear approach) < Freundlich (nonlinear) ARE: Freundlich (linear approach) > Freundlich (nonlinear) EABS: Freundlich (linear approach) < Freundlich (nonlinear) Chisquare (χ^{2}): Freundlich (linear approach) < Freundlich (nonlinear)
For Langmuir adsorption isotherm, error function for nonlinear form was obtained high as compared with linear to exception of R^{2}. R^{2}: Langmuir (linear approach) > Langmuir (nonlinear approach) RSS: Langmuir (linear approach) < Langmuir (nonlinear approach) ARE: Langmuir (linear approach) < Langmuir (nonlinear approach) EABS: Langmuir (linear approach) < Langmuir (nonlinear approach) Chisquare (χ^{2}): Langmuir (linear approach) < Langmuir (nonlinear approach)
Applicability of linear Dubinin–Radushkevich model was found better than nonlinear due to small value of error function (except ARE). R^{2}: Dubinin–Radushkevich (linear approach) < Dubinin–Radushkevich (nonlinear) RSS: Dubinin–Radushkevich (linear approach) < Dubinin–Radushkevich (nonlinear) ARE: Dubinin–Radushkevich (linear approach) > Dubinin–Radushkevich (nonlinear) EABS: Dubinin–Radushkevich (linear approach) < Dubinin–Radushkevich (nonlinear) Chisq/χ^{2}: Dubinin–Radushkevich (linear approach) < Dubinin–Radushkevich (nonlinear)
Linear form of Elovich isotherm, which is also based on multilayer sorption on the surface of sorbent, shows small value of ARE and EABS, but other error functions were found lower for nonlinear form. R^{2}: Elovich (linear approach) > Elovich (nonlinear approach) RSS: Elovich (linear approach) > Elovich (nonlinear approach) ARE: Elovich (linear approach) < Elovich (nonlinear approach) EABS: Elovich (linear approach) < Elovich (nonlinear approach) Chisquare (χ^{2}): Elovich (linear approach) > Elovich (nonlinear approach)
Linear approach for Temkin isotherm was found favorable for adsorption of cadmium ions onto Saccharum arundinaceum with low error function. R^{2}: Temkin (linear approach) > Temkin (nonlinear approach) RSS: Temkin (linear approach) < Temkin (nonlinear approach) ARE: Temkin (linear approach) < Temkin (nonlinear approach) EABS: Temkin (linear approach) < Temkin (nonlinear approach) Chisquare (χ^{2}): Temkin (linear approach) < Temkin (nonlinear approach)
3.6. Thermodynamic Studies
Effect of temperature on adsorption was studied by using temperature range 20, 30, 40, and 50°C with pH 6 and variable initial concentration (30–120 ppm). Adsorption of cadmium onto Saccharum arundinaceum agrowaste was found to increase with increase in temperature. At high temperature, intraparticle diffusion increases and more adsorption sites are created which boost up adsorption phenomenon. Results for thermodynamic study are given in Table 4:

Plot of log versus 1/T was obtained with Rsquared value 0.927. Slop and intercept provide value of ∆H° and ∆S°, respectively, as shown in (10).
∆G° was calculated by employing (1) given in Section 2 in temperature range 292–328 K. Results show a negative value for Gibb’s free energy at all temperature ranges, and ∆G° increases with the increase in temperature. These negative values represent spontaneous nature as well as feasibility of adsorption reaction [36]. Decrease in ∆G° with the increase in temperature reflects better sorption at elevated temperature. The positive value for change in enthalpy is due to endothermic nature of adsorption of cadmium. Enthalpy was also found positive because randomness in system increases due to solidliquid interaction during adsorption phenomenon. Sorption energy calculated by Dubinin–Radushkevich model was found below 1 which is an indication of physical nature of cadmium sorption on the surface of sorbent. E < 8 kJ/mol is representative of physical sorption, and E > 8–16 kJ/mol is due to chemical sorption [37]. For adsorption of cadmium value of E is found below 8, so adsorption of cadmium occurred on the surface and no chemical bonding took place between sorbent and sorbate. Similar results for adsorption of cadmium onto agrowaste were reported in the literature [38].
3.7. Adsorption Kinetics
Adsorption kinetics has prime importance in describing solute uptake rate and time required for adsorption process. In the present work, kinetic study was performed at different time intervals for cadmium adsorption by employing linear and nonlinear forms of pseudofirst and secondorder kinetics. Results indicate that amount of cadmium adsorbed increases with the increase in time interval; however, this increase was sharp in the start of reaction and gradually magnitude of adsorption decreases down. Initially, plenty of active sites were available on the surface of sorbent, so sharp rise in adsorption occurred, but these sites got occupied with the passage of time, so magnitude of adsorption gradually decreases [39].
3.7.1. PseudoFirstOrder Kinetics
For pseudofirstorder kinetic model, log (Q_{e} − Q_{t}) was plotted against time interval and value of k was obtained from slope of the line and Q_{e} from intercept. Initial sorption rate, h, was calculated by the following equation:
Poor correlation was obtained for linear form of the model with low value of R^{2} (0.0918). Results indicate that adsorption of cadmium onto Saccharum arundinaceum does not follow pseudofirstorder kinetics. Nonlinear form of pseudofirstorder kinetics was obtained by using Microsoft Excel 2010 [40].
3.7.2. PseudoSecondOrder Kinetics
Secondorder kinetics is applicable on small amount of initial concentration for determination of initial sorption rate. Different linear and nonlinear forms of pseudosecondorder kinetics are given in Table 3s (Supplementary material). These four linear forms of pseudosecondorder kinetic models were applied on experimental data, and Figures 11s–15s (Supplementary material) show results for these models. Coefficient of determination (R^{2}), found for type 1, was quite high indicating the best fitting of this model on adsorption data of cadmium. Results obtained for pseudosecondorder kinetic are given in Table 4s (Supplementary material). Theoretical results obtained for amount of cadmium adsorbed at equilibrium are found best fitted with experimental data for pseudosecondorder kinetics. For nonlinear form of pseudosecondorder, a computerbased procedure was used in Microsoft Excel 2010 using solver addin method as reported in the literature [41]. For sorption of cadmium onto Saccharum arundinaceum pseudosecondorder kinetic may describe the method of adsorption in quite appropriate way as compared to pseudofirstorder approach. Furthermore, nonlinear form of pseudosecondorder gives close results to experimental data.
3.8. Effect of Interfering Ions
Process of adsorption becomes complicated in case of multicomponent adsorption as many interactions of sorbent and sorbate are involved. Effect of interfering ions was measured by observing adsorption of one metal ion, and then, by addition of interfering ion change in adsorption capacity was noted. Interfering effect was calculated by employing the following equation:where is the % adsorption of mixture of two metal ions and is the % adsorption of pure metal on selected sorbent. If the value of is found equal to 1, then there is no effect of interfering ions on adsorption phenomenon. However, if it is found less than 1, then adsorption capacity is found to be affected by addition of these interfering ions. To study interference, ions were divided into three categories as monovalent, bivalent, and trivalent ions based on valences of ions. One ion was selected from each class of ions to check the interfering effect on adsorption of cadmium ions. Metals are attached on the surface of sorbent through electrostatic forces, and competition of metals for sorbent place is mainly based on metal ion charge and its attraction towards functional groups present on the surface on adsorbent. Results for interference of metal ions are summarized in Table 5. Metals with high charge value were found to have maximum effect on adsorption of cadmium as compared with those with low charge. Anions have also found to affect the adsorption phenomenon of metal ions but their interference is comparatively very less than cations. Adsorption of anions on sorbent is dependent on charge present on the surface of sorbent. Since negatively charged hydroxyl groups are present on the surface of adsorbent, adsorption of anions is not as favored as cations [42].

4. Conclusion
Removal of cadmium was performed by employing the stem powder of Saccharum arundinaceum. In order to generate proper correlation for the removal of cadmium, five adsorption isotherms were applied on experimental data including Freundlich, Langmuir, Dubinin–Radushkevich, Elovich, and Temkin. Error analysis provides information about fitness of these models on experimental data. The model with minimum error was selected best for adsorption data. Order of equilibrium isotherms according to increasing RSS value was found as Temkin > Dubinin–Radushkevich > Elovich > Freundlich > Langmuir.
Linear form of Freundlich and Langmuir models was found best fitted with minimum value of error. Effect of temperature on cadmium adsorption was investigated by thermodynamic analysis, and it was found to increase with increase in temperature. Gibb’s free energy (∆G° = −612.34 at 303 K) revealed spontaneous nature of sorbentsorbate binding reaction, and it followed pseudosecondorder kinetics.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Supplementary Materials
Figures 1s–4s: effect of different parameters (contact time, pH, initial concentration of cadmium, and temperature) which were adjusted during adsorption for maximum removal of cadmium from aqueous media. Figures 5s–9s: graphical results of isothermal study. Data derived from these figures are present in Table 4 of the main manuscript. Figure 10s: thermodynamic studies. Results of Figures 11s–15s of kinetic studies are described in these figures. Linear and nonlinear forms of equilibrium isotherms are given in Tables 1s and 2s. Table 2s provides information about equations of error functions applied on the results. Tables 3s and 4s describe kinetic studies and formulas applied on work for calculation of pseudofirst and secondorder kinetics. (Supplementary Materials)
References
 E. Szabó, Z. Pap, G. Simon, A. Dombi, L. Baia, and K. Hernádi, “New insights on the simultaneous removal by adsorption on organoclays of humic acid and phenol,” Water, vol. 8, no. 1, p. 21, 2016. View at: Publisher Site  Google Scholar
 J. N. Halder and M. N. Islam, “Water pollution and its impact on the human health,” Journal of Environment and Human, vol. 2, no. 1, pp. 36–46, 2015. View at: Publisher Site  Google Scholar
 T. Vaughan, C. W. Seo, and W. E. Marshall, “Removal of selected metal ions from aqueous solution using modified corncobs,” Bioresource Technology, vol. 78, no. 2, pp. 133–139, 2001. View at: Publisher Site  Google Scholar
 P. Húmpola, H. Odetti, J. C. MorenoPiraján, and L. Giraldo, “Activated carbons obtained from agroindustrial waste: textural analysis and adsorption environmental pollutants,” Adsorption, vol. 22, no. 1, pp. 23–31, 2016. View at: Publisher Site  Google Scholar
 V. K. Gupta, M. Gupta, and S. Sharma, “Process development for the removal of lead and chromium from aqueous solutions using red mud—an aluminium industry waste,” Water Research, vol. 35, no. 5, pp. 1125–1134, 2001. View at: Publisher Site  Google Scholar
 Z. Carmen and S. Daniela, “Textile organic dyes–characteristics, polluting effects and separation/elimination procedures from industrial effluents–a critical overview,” in Organic Pollutants Ten Years after the Stockholm Convention–Environmental and Analytical Update, pp. 55–81, InTech, Croatia, 2012. View at: Google Scholar
 G. Crini, “Recent developments in polysaccharidebased materials used as adsorbents in wastewater treatment,” Progress in Polymer Science, vol. 30, no. 1, pp. 38–70, 2005. View at: Publisher Site  Google Scholar
 M. Vidali, “Bioremediation. An overview,” Pure and Applied Chemistry, vol. 73, no. 7, pp. 1163–1172, 2001. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 Y. Salameh, N. AlLagtah, M. N. M. Ahmad, S. J. Allen, and G. M. Walker, “Kinetic and thermodynamic investigations on arsenic adsorption onto dolomitic sorbents,” Chemical Engineering Journal, vol. 160, no. 2, pp. 440–446, 2010. View at: Publisher Site  Google Scholar
 O. Carmody, R. Frost, Y. Xi, and S. Kokot, “Surface characterisation of selected sorbent materials for common hydrocarbon fuels,” Surface Science, vol. 601, no. 9, pp. 2066–2076, 2007. View at: Publisher Site  Google Scholar
 M. Iqbal, A. Saeed, and S. I. Zafar, “FTIR spectrophotometry, kinetics and adsorption isotherms modeling, ion exchange, and EDX analysis for understanding the mechanism of Cd^{2+} and Pb^{2+} removal by mango peel waste,” Journal of Hazardous Materials, vol. 164, no. 1, pp. 161–171, 2009. View at: Publisher Site  Google Scholar
 Y. S. Ho, C. T. Huang, and H. W. Huang, “Equilibrium sorption isotherm for metal ions on tree fern,” Process Biochemistry, vol. 37, pp. 1421–1430, 2002. View at: Publisher Site  Google Scholar
 L. S. Chan, W. H. Cheung, S. J. Allen, and G. McKay, “Error analysis of adsorption isotherm models for acid dyes onto bamboo derived activated carbon,” Chinese Journal of Chemical Engineering, vol. 20, no. 3, pp. 535–542, 2012. View at: Publisher Site  Google Scholar
 J. Febrianto, A. N. Kosasih, J. Sunarso, Y.H. Ju, N. Indraswati, and S. Ismadji, “Equilibrium and kinetic studies in adsorption of heavy metals using biosorbent: a summary of recent studies,” Journal of Hazardous Materials, vol. 162, no. 12, pp. 616–645, 2009. View at: Publisher Site  Google Scholar
 S.J. Park and Y.S. Jang, “Pore structure and surface properties of chemically modified activated carbons for adsorption mechanism and rate of Cr(VI),” Journal of Colloid and Interface Science, vol. 249, no. 2, pp. 458–463, 2002. View at: Publisher Site  Google Scholar
 M. Salah Azab and P. J. Peterson, “The removal of cadmium from water by the use of biological sorbents,” Water Science and Technology, vol. 21, no. 12, pp. 17051706, 1989. View at: Publisher Site  Google Scholar
 A. Kapoor and T. Viraraghavan, “Fungal biosorption—an alternative treatment option for heavy metal bearing wastewaters: a review,” Bioresource Technology, vol. 53, no. 3, pp. 195–206, 1995. View at: Publisher Site  Google Scholar
 B. Tansel and P. Nagarajan, “SEM study of phenolphthalein adsorption on granular activated carbon,” Advances in Environmental Research, vol. 8, no. 34, pp. 411–415, 2004. View at: Publisher Site  Google Scholar
 O. B. Belskaya, A. V. Lavrenov, I. G. Danilova, M. O. Kazakov, R. M. Mironenko, and V. A. Likholobov, “FTIR spectroscopy of adsorbed probe molecules for analyzing the surface properties of supported Pt (Pd) catalysts,” in Infrared Spectroscopy–Materials Science, Engineering and Technology, INTECH Open Access Publisher, London, UK, 2012. View at: Google Scholar
 B. Stuart, Infrared Spectroscopy, Wiley, Hoboken, NJ, USA, 2005.
 S. Yang, J. Li, D. Shao, J. Hu, and X. Wang, “Adsorption of Ni (II) on oxidized multiwalled carbon nanotubes: effect of contact time, pH, foreign ions and PAA,” Journal of Hazardous Materials, vol. 166, no. 1, pp. 109–116, 2009. View at: Publisher Site  Google Scholar
 M. Malandrino, O. Abollino, A. Giacomino, M. Aceto, and E. Mentasti, “Adsorption of heavy metals on vermiculite: influence of pH and organic ligands,” Journal of Colloid and Interface Science, vol. 299, no. 2, pp. 537–546, 2006. View at: Publisher Site  Google Scholar
 B. B. Johnson, “Effect of pH, temperature, and concentration on the adsorption of cadmium on goethite,” Environmental Science and Technology, vol. 24, no. 1, pp. 112–118, 1990. View at: Publisher Site  Google Scholar
 T. Teka and S. Enyew, “Study on effect of different parameters on adsorption efficiency of low cost activated orange peels for the removal of methylene blue dye,” International Journal of Innovation and Scientific Research, vol. 8, pp. 106–111, 2014. View at: Google Scholar
 G. Ks and S. L. Belagali, “Removal of heavy metals and dyes using low cost adsorbents from aqueous medium, a review,” IOSR Journal of Environmental Science, Toxicology and Food Technology, vol. 4, no. 3, pp. 56–68, 2013. View at: Publisher Site  Google Scholar
 I. D. Mall, V. C. Srivastava, and N. K. Agarwal, “Removal of OrangeG and Methyl Violet dyes by adsorption onto bagasse fly ash—kinetic study and equilibrium isotherm analyses,” Dyes and pigments, vol. 69, no. 3, pp. 210–223, 2006. View at: Publisher Site  Google Scholar
 B. H. Hameed, I. A. W. Tan, and A. L. Ahmad, “Adsorption isotherm, kinetic modeling and mechanism of 2, 4, 6trichlorophenol on coconut huskbased activated carbon,” Chemical Engineering Journal, vol. 144, no. 2, pp. 235–244, 2008. View at: Publisher Site  Google Scholar
 Y. Liu and Y.J. Liu, “Biosorption isotherms, kinetics and thermodynamics,” Separation and Purification Technology, vol. 61, no. 3, pp. 229–242, 2008. View at: Publisher Site  Google Scholar
 R. M. AllenKing, 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. View at: Publisher Site  Google Scholar
 C. P. Schulthess and D. L. Sparks, “Equilibriumbased modeling of chemical sorption on soils and soil constituents,” in Advances in Soil Science, pp. 121–163, Springer, Berlin, Germany, 1991. View at: Google Scholar
 V. Patrulea, A. Negrulescu, M. M. Mincea, L. D. Pitulice, O. B. Spiridon, and V. Ostafe, “Optimization of the removal of copper(II) ions from aqueous solution on chitosan and crosslinked chitosan beads,” BioResources, vol. 8, no. 1, 2013. View at: Publisher Site  Google Scholar
 R. Rostamian, M. Najafi, and A. A. Rafati, “Synthesis and characterization of thiolfunctionalized silica nano hollow sphere as a novel adsorbent for removal of poisonous heavy metal ions from water: Kinetics, isotherms and error analysis,” Chemical Engineering Journal, vol. 171, no. 3, pp. 1004–1011, 2011. View at: Publisher Site  Google Scholar
 E. Demirbas, M. Kobya, and A. E. Konukman, “Error analysis of equilibrium studies for the almond shell activated carbon adsorption of Cr(VI) from aqueous solutions,” Journal of Hazardous Materials, vol. 154, no. 1–3, pp. 787–794, 2008. View at: Publisher Site  Google Scholar
 D. Gusain, V. Srivastava, M. Sillanpää, and Y. C. Sharma, “Kinetics and isotherm study on adsorption of chromium on nano crystalline iron oxide/hydroxide: linear and nonlinear analysis of isotherm and kinetic parameters,” Research on Chemical Intermediates, vol. 42, no. 9, pp. 7133–7151, 2016. View at: Publisher Site  Google Scholar
 B. Singha and S. K. Das, “Biosorption of Cr(VI) ions from aqueous solutions: Kinetics, equilibrium, thermodynamics and desorption studies,” Colloids and Surfaces B: Biointerfaces, vol. 84, no. 1, pp. 221–232, 2011. View at: Publisher Site  Google Scholar
 J. RomeroGonzález, J. R. PeraltaVidea, E. Rodrı́guez, S. L. Ramirez, and J. L. GardeaTorresdey, “Determination of thermodynamic parameters of Cr(VI) adsorption from aqueous solution onto Agave lechuguilla biomass,” Journal of Chemical Thermodynamics, vol. 37, no. 4, pp. 343–347, 2005. View at: Publisher Site  Google Scholar
 E. Mekonnen, M. Yitbarek, and T. R. Soreta, “Kinetic and thermodynamic studies of the adsorption of Cr(VI) onto some selected local adsorbents South African,” Journal of Chemistry, vol. 68, pp. 45–52, 2015. View at: Publisher Site  Google Scholar
 Y. Bulut and H. Aydın, “A kinetics and thermodynamics study of methylene blue adsorption on wheat shells,” Desalination, vol. 194, no. 1–3, pp. 259–267, 2006. View at: Publisher Site  Google Scholar
 J. Lin and L. Wang, “Comparison between linear and nonlinear forms of pseudofirstorder and pseudosecondorder adsorption kinetic models for the removal of methylene blue by activated carbon,” Frontiers of Environmental Science and Engineering in China, vol. 3, no. 3, pp. 320–324, 2009. View at: Publisher Site  Google Scholar
 Y.S. Ho, “Secondorder kinetic model for the sorption of cadmium onto tree fern: a comparison of linear and nonlinear methods,” Water Research, vol. 40, no. 1, pp. 119–125, 2006. View at: Publisher Site  Google Scholar
 D. Mohan and K. P. Singh, “Singleand multicomponent adsorption of cadmium and zinc using activated carbon derived from bagasse—an agricultural waste,” Water Research, vol. 36, no. 9, pp. 2304–2318, 2002. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Fozia Batool et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.