Research Article  Open Access
Zvezdelina Lyubenova Yaneva, Bogdana Koumanova Koumanova, Nedyalka Valkanova Georgieva, "Linear and Nonlinear Regression Methods for Equilibrium Modelling of pNitrophenol Biosorption by Rhizopus oryzae: Comparison of Error Analysis Criteria", Journal of Chemistry, vol. 2013, Article ID 517631, 10 pages, 2013. https://doi.org/10.1155/2013/517631
Linear and Nonlinear Regression Methods for Equilibrium Modelling of pNitrophenol Biosorption by Rhizopus oryzae: Comparison of Error Analysis Criteria
Abstract
The study assessed the applicability of Rhizopus oryzae dead fungi as a biosorbent medium for pnitrophenol (pNP) removal from aqueous phase. The extent of biosorption was measured through five equilibrium sorption isotherms represented by the Langmuir, Freundlich, RedlichPeterson, multilayer and FritzSchlunder models. Linear and nonlinear regression methods were compared to determine the bestfitting equilibrium model to the experimental data. A detailed error analysis was undertaken to investigate the effect of applying seven error criteria for the determination of the singlecomponent isotherm parameters. According to the comparison of the error functions and to the estimation of the corrected Akaike information criterion (), the Freundlich equation was ranked as the first and the FritzSchlunder as the second bestfitting models describing the experimental data. The present investigations proved the high efficiency (94%) of Rhizopus Oryzae as an alternative adsorbent for pNP removal from aqueous phase and revealed the mechanism of the separation process.
1. Introduction
Nitrophenols are classified as moderately to highly toxic even at low concentrations when present in wastewaters. Thus, they are a source of serious social and hygienic problems as important occupational and environmental pollutants [1]. The undesirable impact of these recalcitrant organics on organisms, as well as their considerable release through industrial wastes, provoked the researchers to search for innovative effective technologies for the remediation of nitrophenolsrich wastes before their discharge into natural water streams [2]. Nitroaromatic compounds are widely used as pesticides, explosives, solvents, and intermediates in the synthesis of dyes and other chemicals. Many of these compounds and their transform products are of significant toxicological concern [3].
The monosubstituted pnitrophenol (pNP) could enter the environment during its production and use in hightemperature coal conversion, as a byproduct in the enzymatic hydrolysis of parathion, methyl parathion, and Nacetylpaminophenol. It is also found in suspended particulate matter in the atmosphere, originating mostly from secondary photochemical reactions in the air [3]. The toxicity of pnitrophenol on biological systems has led to its classification as a priority pollutant by the United States Environmental Protection Agency (EPA) [4].
The literature reports a number of studies on the efficiency of various adsorbents for phenols and substituted phenols removal from wastewaters [5]. The high adsorption activity of various biopolymers (unmodified/modified chitosan, lignin, etc) towards phenols, chlorophenols, and dyes was scientifically proven [6–8]. Fungi belonging to the genera Rhizopus and Penicillium have already been studied as a potential biomass for the removal of heavy metals (Cu, Zn, Cr, Cd) [9–11], nonylphenol [12], polycyclic aromatic hydrocarbons (PACs) [13], cyanide [14], and dyes from aqueous solutions [8]. The saprophytic microorganisms Rhizopus oryzae were proved as a good adsorbent for the removal of the organochloropesticide lindane [15], and rhodamine B from water [16].
The walls of R. oryzae consist of chitin and chitosan. The mechanism of biosorption is quite complicated. The capacity of dead cells might be higher, equal or even lower when compared to that of live cells. In most cases, however, the application of dead biomass is more effective [17]. Chitin is the second most abundant natural carbohydrate polymer next to cellulose. Many applications of R. oryzae are due to the secondary amino groups of chitosan which show polycationic chelating and filmforming properties along with high solubility in dilute acetic acid [18].
Over the past few decades, linear regression has been developed as a major option in designing adsorption systems [19–23]. However, recent investigations have indicated the growing discrepancy (between the predictions and experimental data) and disability of the model, propagating towards a different outcome [24]. A number of error functions (sum of the squares of errors, ERRSQ; hybrid fractional error function, HYBRID; Marquardt’s percent standard deviation, MPSD; average relative error, ARE; Sum of the absolute errors, EABS; chisquare function, ; Akaike information criterion, AIC, etc.) were highlighted and discussed for the comparisons of linear and nonlinear isotherm models [24–26]. According to recent studies, the expanding of the nonlinear isotherms represents a potentially viable and powerful tool, leading to the superior improvement in the area of adsorption science [27]. However, according to Foo and Hameed [24], further explorations in this area are recommended. Moreover, these error estimating functions do not take into account the number of parameters in the models. Hence, they cannot be used for selection of model(s) appropriateness and determination of its/their magnitude. The literature reports studies applying the Akaike information criterion (AIC) for model selection and ranking in the field of veterinary medicine [28], genetics [29], aquatic biogeochemical modeling [30], FTIR spectroscopy [31], pharmacy [32], and so forth. The number of studies subjected to the applicability of AIC for statistical modelling of toxic organics sorption on alternative sorbents, however, is limited.
In this context, the scope of the present study was to investigate the applicability of Rhizopus oryzae as an alternative “ecofriendly” adsorbent for pnitrophenol removal from aqueous phase at equilibrium conditions. Linear and nonlinear regression analysis was performed to determine model parameters. Seven error functions were applied to evaluate, compare, and rank the feasibility of the five applied isotherm models (Langmuir, Freundlich, RedlichPeterson, multilayer, and FritzSchlunder).
2. Materials and Methods
2.1. Adsorbate
pNitrophenol (pNP) (Merck, 99%), without further purification, was used as a sorbate in the investigations. Its physicochemical characteristics are presented elsewhere [33].
2.2. Adsorbent
Rhizopus oryzae used in the recent study as a sorbent was supplied by the International Mycological Institute in Surrey, UK, in the form of IMI strain 266680. The microorganisms were isolated from a soil in Sri Lanka. The spores dried at low temperature were reactivated and cultivated in malt extract (17 g dm^{−3} malt extract and 3 g dm^{−3} mycological peptone, dissolved in distilled water at pH ). The malt extract was inoculated by a standard sterile method and incubated at 32°C for 3 days in a platform shaker at 175 rpm. Three ceramic granules were added to each of the batch reactors with malt extract to limit the microorganisms threadlike growth. The biomass obtained was precisely washed out consecutively with maternal lye and distilled water and then dried in an oven at 50°C. The dried biomass was ground in a hammer mill and screened. The fraction used in the recent investigations was 0.15−0.50 mm.
2.3. Equilibrium Studies
The equilibrium experiments were accomplished using model solutions of pNP in distilled water. The investigations were carried out at temperature °C and pH 6.1. Solutions were in the concentration range 5–25 mg dm^{−3}. Known amounts of Rhizopus oryzae, 0.3 g, were added to 100 cm^{3} of the model solutions in screw cap jars. The jars were shaken on a platform shaker. The solute from each jar was then filtered. The residual pNP concentrations in the liquid phase () were determined spectrophotometrically. All investigations were performed in triplicate. The corresponding values of pNP solid phase concentrations () were calculated by the mass balance: where , mg dm^{−3}, is the initial adsorbate concentration in the liquid phase, , and , , is the adsorbent mass.
SPECORD UVVIS, Carl Zeiss Jena, spectrophotometer was used for concentration determinations at maximum absorption wavelength 228 nm. The pH was measured using an LHP 403T TACUSSEL pHmeter. Blanks containing no adsorbate and replicates of each adsorption point were used for each series of experiments.
2.4. Isotherm Modeling
The design and efficient operation of adsorption processes require equilibrium adsorption data for use in kinetic, dynamic, and mass transfer models [34]. In order to optimize the design of a specific sorbate/sorbent system, it is important to establish the most appropriate correlation for the experimental equilibrium data [35].
The biosorption behavior of pNP on R. oryzae in the present research was modelled by the Langmuir, Freundlich, RedlichPeterson, multilayer, and FritzSchlunder (fourparameter) isotherm equations (Table 1).
2.5. Error Analysis
In the present study, linear and nonlinear regression analysis was performed to determine the values of the isotherm model parameters. Six different error functions (, ERRSQ, HYBRID, MPSD, ARE, and EABS) were examined to evaluate the applicability of each model isotherm equation to the experimental data (Table 2) using the solver addin functions of Microsoft Excel software.

3. Results and Discussion
In order to assess the fate of nitrophenols in wastewater and to control their mobility and reactivity during remediation processes, the sorption behaviour and mechanism of these toxic contaminants must be understood and revealed. The knowledge of sorbate/sorbent adsorption behaviour at equilibrium is essential for environmental engineering and science as the derived isotherms reveal the specific relation between the pollutant concentration and its uptake degree by the solid phase at constant temperature.
In this section, among the five studied isotherms models, the bestfitting one was determined by the use of seven wellknown error functions to calculate the error deviation between experimental and predicted equilibrium adsorption data, after both linear and nonlinear analysis. In all of the error methods it was assumed that both the liquid phase concentration and the solid phase concentration contribute equally to weighting the error criterion for the model solution procedure. The experimental investigations were conducted at biomass concentration 3.00 g dm^{−3}.
The values of the model parameters and the isotherm error deviation data for the Langmuir, Freundlich, and RedlichPeterson equations determined by linear regression analysis and for the multilayer model obtained through the secondorder polynomial form of (7) (see Table 1), are presented in Table 3.
 
Note. The highest and the lowest ERRSQ, HYBRID, MPSD, ARE, and EABS are in bold. 
Alternative isotherm parameters were also determined by nonlinear regression using five error functions (ERRSQ, HYBRID, MPSD, ARE, and EABS) and the corrected Akaike information criterion (). The values of the Langmuir, Freundlich, RedlichPeterson, multilayer and FritzSchlunder model constants and the isotherm error deviation data are presented in Table 4.
 
Note. The lowest ERRSQ, HYBRID, MPSD, ARE, EABS, and AIC_{c} are in bold. 
The experimental data points of pNP adsorption on Rhizopus oryzae dead biomass (Figures 1, 2, 4, and 5) outlined a distinct steep vertical section in the lowconcentration region, where a sharp leap of nitrophenol solid phase concentration from 0 to 2.63 mg g^{−1} was distinguished. Consequently, higher extent of adsorption during the initial stages of the process could be expected. Besides, a clearly marked plateau, that is, a horizontal section in the higher concentration range, was not observed. The latter could be interpreted by means of the presence of a number of vacant active sites, as well as a larger continuance of the adsorption process in the later stages. The abovestated assumptions were consistent with the higher values of the equilibrium constant for the first layer adsorption ( 1.935: linear regression; 2.7971−3.2225: nonlinear regression) than that for multilayer sorption ( 0.026: linear regression; 0.0265−0.0336: nonlinear regression), calculated by the multilayer model (Tables 3 and 4). The experimentally established maximum equilibrium capacity of the biomass was 4.6 mg g^{−1}, the monolayer capacity according to the linear Langmuir model— 5.01 mg g^{−1} (Table 3), and according to the nonlinear model—4.28−7.91 mg g^{−1} (Table 4). The second polynomial form of the multilayer equation predicted a maximum monolayer adsorption capacity 3.399 mg g^{−1} (Table 3), while the values of this parameter calculated by applying the nonlinear approach ranged between 3.0834 and 3.401 mg g^{−1} (Table 4).
The comparative analysis between the values of the error functions, obtained through the linear approach, outlined the threeparameter RedlichPeterson model as the one with the highest 0.9997 (Table 3). However, the Freundlich equation characterized with the lowest ERRSQ, HYBRID, ARE, MPSD and EABS error values. The data from Table 4 showed analogous tendency. Hence, is seemed that the Freundlich equation was the most suitable model presenting the most satisfactory description of the studied biosorption phenomenon.
Statistically, it is expected that the higher the number of parameters in a model equation, the closer the theoretical estimates should be to the empirical data. Moreover the error functions HYBRID and MPSD could be accepted as the most indicative, adequate and essentially meaningful when determining the best fit isotherm model, as the number of the isotherm parameters is accounted only by them [35]. The data in Tables 3 and 4, however, displayed the contrary as the threeparameter (RedlichPeterson and multilayer) and the fourparameter (FritzSchlunder) models characterized with lower extend of suitability to the experimental data points due to the higher error values.
To prove the latter observations, the modes of the experimental and model isotherms obtained on the basis of the linear and nonlinear approach were also compared (Figure 1, 2, 4, 5). Obviously, the Langmuir isotherms presented in Figure 1 did not correlate the experimental equilibrium data satisfactorily. It was observed that the HYBRID and EABS methods yielded the best fit in the low concentration range, while the LTFM isotherm could be applied for modeling the high concentration region.
Among the five applied equations, the Freundlich (Figure 2) and FritzSchlunder (Figure 5) isotherms could be identified as the most suitable for modelling the equilibrium sorption behaviour of pNP on R. oryzae covering the entire concentration range. Besides, the individual model curves, derived on the basis of the studied error functions, practically coincided in both cases. Considering the theoretical bases of the Freundlich model, the studied separation process could be described either as nonideal and reversible sorption, not restricted to monolayer formation on heterogeneous surface, or as multilayer sorption with nonuniform distribution of adsorption heat and affinities over the heterogeneous surface [24]. Besides, as the value of the parameter is below unity according to all implemented error analysis procedures, it is indicative of a chemisorption process.
As stated earlier (Section Introduction), chitin and chitosan are the main components of R. oryzae cell walls. Due to the better hydrophilic properties of chitosan, it expands in an aqueous medium, the active sites which predetermine hydrogen bonding, become more pronounced and easily accessible. The amino (–NH_{2}) and hydroxyl (–OH) groups serve as the coordination and reaction sites. Besides, the pKa value of amino group (R–NH_{2}) in the structure of chitosan is 6.3, and amino group dissociates partly into even at pH 6.1 (Figure 3). The pKa value of pnitrophenol is around 7.1, whereby at certain pH below the pKa of the dissociating solutes, pnitrophenol exists as neutral form and above the pKa value, pnitrophenol exists in ionic form [4]. O–H bond can be broken off easily, and nitro group causes to earn the resonance stability to structure by helping to the delocalization of negative charge [7]. Hence, there is a possibility of chemical interaction between this positive charge and negative charge existent and delocalized in the anionic structure of the organic molecule (Figure 3). Chitin is carrying on linear amino group per glucose unit and thus exhibitis much higher uptake capacity. The amino group has an electron pair available for coordination and behaves like a strong Lewis base. The amine nitrogen on each chitin monomer has been suggested as the active site for 4NP adsorption [51].
Hence, it is supposed that the adsorption of pnitrophenol molecules was due mainly to the formation of Hbonds between chitin and chitosan functional groups and NO_{2}– and OH– groups in the organic molecules (Figure 3).
The relatively high value of the regression coefficient obtained through the second polynomial form of the multilayer model ( 0.9967) (Table 3) confirmed its applicability for describing the experimental data. The plots presented in Figure 4, however, demonstrated a deviation of all model isotherms for the middle concentration region. Probably, the biosorption of the mononitrophenol was accomplished predominantly on the first biosorbent layers, but the possibility of multilayer sorption, although at a lower extend, could not be neglected.
The comparison between the values of the six error estimating functions and the modes of experimental versus model isotherms outlined the Freundlich (Figure 2) and/or FritzSchlunder (Figure 5) equations as the most adequate models describing pNP biosorption on R. oryzae. In the context of the current work, the corrected Akaike information criterion was used as to verify the conclusions for the best fitting model withdrawn on bases of error analysis, so to rank the five isotherm models.
The AIC developed by Akaike is a methodology for model selection in a situation where more than one model has been fitted to experimental data and screening of the candidate models is crucial to the objectives of the research work. Akaike’s general approach not only allows the best model to be identified, but also allows the ranking of the rest of the models under consideration [46].
The data in Table 4 show that the for the Freundlich model has the minimum value ( −14.44). This implies that the theoretical data obtained from this model fits the experimental results better than the other four equilibrium isotherm equations. Consequently, this statistical tool confirmed the results from the error analysis and proved the better correlation of the experimental data with the nonlinear forms of the conventional twoparameter Freundlich equation. The ranking of the isotherm equations was made on the basis of the Relative Akaike Weight (RAW), , (Table 2). The values of (Table 4) clearly distinguished the appropriateness of the five isotherm equations ranking the twoparameter Freundlich model as the first ( 1.00) and the fourparameter FritzSchlunder model as the second () best fitted ones to the experimental data. Similar results had been obtained by Chatterjee et al. [49]. They observed that experimental data gave better fitting to Langmuir than to RedlichPeterson model using AIC to rank them. The investigations of Mutua [46] outlined that the Freundlich model ranked first among seven isotherm equations used to model the adsorption of Cu^{2+} and Cd^{2+} onto bentonite and modified bentonite clay.
Comparative analyses between the presented experimental results and the literature cited biosorption capacity of various R. oryzaebased biosorbents were accomplished and presented in Table 5.

Considering the values of the individual system parameters (initial sorbate concentration, biosorbent concentration, etc.), it could be concluded that the maximum biosorption capacity of R. oryzae dead biomass towards pNP ( 4.6 mg g^{−1}) was satisfactory. Besides, the present investigations proved the high efficiency (94%) of Rhizopus oryzae, obtained at 5 mg dm^{−3} and 3 g dm^{−3}, as an alternative adsorbent for pNP removal from aqueous phase and revealed the mechanism of the separation process at laboratory scale. Due to the harmful effects of these organic compounds and the urge for sustainable development for the natural water resources at global scale, however, the wastewaters containing them must be treated before being discharged to receiving natural water bodies. Consequently, the future perspectives of the applicability of the proposed “ecofriendly” biosorbent for industrial purposes include largescale investigations with real wastewaters containing mononitrophenols and/or contaminated natural waters.
4. Conclusion
The singlecomponent biosorption of pNP on Rhizopus oryzae dead biomass from aqueous medium was investigated on the basis of equilibrium studies. The experimental data was modelled and evaluated using five isotherm models and seven optimization and error functions, including the corrected Akaike information criterion () as a statistical estimating and ranking tool. The linear transform model provided the highest regression coefficient for the case of RedlichPeterson isotherm. The analysis of the other error functions for linear/nonlinear optimization and outlined the conventional twoparameter Freundlich equation as generally the bestfitting model among the other two, three, and fourparameter isotherms applied. The laboratory studies proved the high efficiency (94% uptake extend) of Rhizopus oryzae dead biomass as an alternative adsorbent for removal of the toxic and recalcitrant organic contaminant pNP from aqueous phase and outlined its future applicability in largescale wastewater treatment technologies.
Nomenclature
:  FritzSchlunder constant, dm^{3} g^{−1} 
:  Langmuir isotherm constant, dm^{3} g^{−1} 
:  RedlichPeterson isotherm constant, dm^{3} g^{−1} 
AIC:  Akaike information criterion 
:  Corrected Akaike information criterion 
ARE:  The average relative error 
:  FritzSchlunder constant, dm^{3} g^{−1} 
:  RedlichPeterson isotherm constant 
:  Equilibrium adsorbate concentration in the liquid phase, mg dm^{−3} 
:  Initial adsorbate concentration in the liquid phase, mg dm^{−3} 
:  Particle diameter, mm 
EABS:  The sum of the absolute errors 
ERRSQ:  The sum of the squares of errors 
FSch:  FritzSchlunder model 
HYBRID:  The hybrid fractional error function 
:  Langmuir isotherm constant, dm^{3} g^{−1} 
:  Freundlich isotherm constant, dm^{3} g^{−1} 
:  RedlichPeterson isotherm constant, dm^{3} g^{−1} 
:  equilibrium constant for the first layer adsorption in the multilayer isotherm model 
:  Equilibrium constant for multilayer adsorption in the multilayer isotherm model 
LTFM:  The linear transform model 
:  Adsorbent concentration, g dm^{−3} 
MPSD:  Marquardt’s percent standard deviation 
:  number of isotherm parameters 
:  the heterogeneity factor in the Freundlich model 
:  Number of experimental data points 
:  Maximum monolayer adsorption capacity in the multilayer isotherm model, mg g^{−1} 
:  The equilibrium sorbate concentration in the solid phase, mg g^{−1} 
:  Experimental pNP equilibrium solid phase concentration, mg g^{−1} 
:  Model calculated pNP equilibrium solid phase concentration, mg g^{−1} 
:  The sorption capacity at equilibrium and at time , mg g^{−1} 
:  Correlation coefficient 
SSR:  The sum of the squares for the residual 
:  Temperature, K 
:  Time, min 
:  Solution volume, dm^{3} 
:  Adsorbent mass, g. 
:  FritzSchlunder exponent 
:  FritzSchlunder exponent 
:  Maximum absorbance wavelength, nm 
:  Akaike weight for the th model 
:  The difference between the of the best fitting model and that of model . 
References
 A. Ayar, S. Gürsal, A. A. Gürten, and O. Gezici, “On the removal of some phenolic compounds from aqueous solutions by using a sporopolleninbased ligandexchange fixed bed—isotherm analysis,” Desalination, vol. 219, no. 1–3, pp. 160–170, 2008. View at: Publisher Site  Google Scholar
 K. Lin, J. Pan, Y. Chen, R. Cheng, and X. Xu, “Study the adsorption of phenol from aqueous solution on hydroxyapatite nanopowders,” Journal of Hazardous Materials, vol. 161, no. 1, pp. 231–240, 2009. View at: Publisher Site  Google Scholar
 W. Huang, C. Yao, S. Jin, S. Ying, and X. Shen, “Sorption of pnitrophenol onto sediment in the presence of cetylpyridinium chloride and Pb(NO_{3})_{2}: influence of pH,” Journal of Hazardous Materials, vol. 155, no. 12, pp. 225–229, 2008. View at: Publisher Site  Google Scholar
 W. S. W. Ngah and S. Fatinathan, “Chitosan flakes and chitosanGLA beads for adsorption of pnitrophenol in aqueous solution,” Colloids and Surfaces A, vol. 277, no. 1–3, pp. 214–222, 2006. View at: Publisher Site  Google Scholar
 P. Barkakati, A. Begum, M. L. Das, and P. G. Rao, “Adsorptive separation of Ginsenoside from aqueous solution by polymeric resins: equilibrium, kinetic and thermodynamic studies,” Chemical Engineering Journal, vol. 161, no. 12, pp. 34–45, 2010. View at: Publisher Site  Google Scholar
 I. Uzun and F. Güzel, “Kinetics and thermodynamics of the adsorption of some dyestuffs and pnitrophenol by chitosan and MCMchitosan from aqueous solution,” Journal of Colloid and Interface Science, vol. 274, no. 2, pp. 398–412, 2004. View at: Publisher Site  Google Scholar
 I. Uzun and F. Güzel, “Rate studies on the adsorption of some dyestuffs and pnitrophenol by chitosan and monocarboxymethylated(mcm)chitosan from aqueous solution,” Journal of Hazardous Materials, vol. 118, no. 1–3, pp. 141–154, 2005. View at: Publisher Site  Google Scholar
 S. J. Allen and B. Koumanova, “Decolourization of water/wastewater using adsorption (Review),” Journal of the University of Chemical Technology and Metallurgy, vol. 40, no. 3, pp. 175–192, 2005. View at: Google Scholar
 B. J. Mcafee, W. D. Gould, J. C. Nadeau, and A. C. A. Da Costa, “Biosorption of metal ions using chitosan, chitin, and biomass of Rhizopus oryzae,” Separation Science and Technology, vol. 36, no. 14, pp. 3207–3222, 2001. View at: Publisher Site  Google Scholar
 Q. I. Yang, J. L. Wang, and Z. Xing, “Biosorption of cadmium by fungal biomass of Aspergillus niger,” Biomedical and Environmental Sciences, vol. 18, no. 3, pp. 141–145, 2005. View at: Google Scholar
 K. C. Bhainsa and S. F. D'Souza, “Removal of copper ions by the filamentous fungus, Rhizopus oryzae from aqueous solution,” Bioresource Technology, vol. 99, no. 9, pp. 3829–3835, 2008. View at: Publisher Site  Google Scholar
 W. Lang, C. Dejma, S. Sirisansaneeyakul, and N. Sakairi, “Biosorption of nonylphenol on dead biomass of Rhizopus arrhizus encapsulated in chitosan beads,” Bioresource Technology, vol. 100, no. 23, pp. 5616–5623, 2009. View at: Publisher Site  Google Scholar
 B. Ma, H. Chen, M. Xu et al., “Adsorption of polycyclic aromatic hydrocarbons (PAHs) on Rhyzopus Oryzae cell wall: application of cosolvent models for validating the cell wallwater partition coefficient, Approaches for Measuring Kinetic and Thermodynamic Properties of Processes at the Mineral/Water Interface,” in Proceedings of the International Annual Meetings ASA, CSSA and SSSA, pp. 303–304, Long Beach, Calif, USA, November 2010. View at: Google Scholar
 R. R. Dash, C. Balomajumder, and A. Kumar, “Cyanide removal by combined adsorption and biodegradation process,” Iranian Journal of Environmental Health Science & Engineering, vol. 3, no. 2, pp. 91–96, 2006. View at: Google Scholar
 S. Ghosh, S. K. Das, A. K. Guha, and A. K. Sanyal, “Adsorption behavior of lindane on Rhizopus oryzae biomass: physicochemical studies,” Journal of Hazardous Materials, vol. 172, no. 1, pp. 485–490, 2009. View at: Publisher Site  Google Scholar
 S. K. Das, J. Bhowal, A. R. Das, and A. K. Guha, “Adsorption behavior of rhodamine B on Rhizopus oryzae biomass,” Langmuir, vol. 22, no. 17, pp. 7265–7272, 2006. View at: Publisher Site  Google Scholar
 B. Koumanova, P. Peeva, S. J. Allen et al., “Biosorption from aqueous solutions by eggshell membranes and Rhizopus oryzae: equilibrium and kinetic studies,” Journal of Chemical Technology & Biotechnology, vol. 77, pp. 539–545, 2002. View at: Google Scholar
 J. S. Rhee, M. W. Jung, and K. J. Paeng, “Evaluation of chitin and chitosan as a sorbent for the preconcentration of phenol and chlorophenols in water,” Analytical Sciences, vol. 14, no. 6, pp. 1089–1092, 1998. View at: Google Scholar
 O. Hamdaoui and E. Naffrechoux, “Modeling of adsorption isotherms of phenol and chlorophenols onto granular activated carbon. Part II. Models with more than two parameters,” Journal of Hazardous Materials, vol. 147, no. 12, pp. 401–411, 2007. View at: Publisher Site  Google Scholar
 Q. S. Liu, T. Zheng, P. Wang, J. P. Jiang, and N. Li, “Adsorption isotherm, kinetic and mechanism studies of some substituted phenols on activated carbon fibers,” Chemical Engineering Journal, vol. 157, no. 23, pp. 348–356, 2010. View at: Publisher Site  Google Scholar
 Z. Zhang, Z. Zhang, Y. Fernández et al., “Adsorption isotherms and kinetics of methylene blue on a lowcost adsorbent recovered from a spent catalyst of vinyl acetate synthesis,” Applied Surface Science, vol. 256, no. 8, pp. 2569–2576, 2010. View at: Publisher Site  Google Scholar
 J. Rahchamani, H. Z. Mousavi, and M. Behzad, “Adsorption of methyl violet from aqueous solution by polyacrylamide as an adsorbent: isotherm and kinetic studies,” Desalination, vol. 267, no. 23, pp. 256–260, 2011. View at: Publisher Site  Google Scholar
 V. M. Vučurović, R. N. Razmovski, and M. N. Tekić, “Methylene blue (cationic dye) adsorption onto sugar beet pulp: equilibrium isotherm and kinetic studies,” Journal of the Taiwan Institute of Chemical Engineers, vol. 43, no. 1, pp. 108–111, 2012. 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
 K. V. Kumar, K. Porkodi, and F. Rocha, “Isotherms and thermodynamics by linear and nonlinear regression analysis for the sorption of methylene blue onto activated carbon: comparison of various error functions,” Journal of Hazardous Materials, vol. 151, no. 23, pp. 794–804, 2008. View at: Publisher Site  Google Scholar
 H. LeónSantiestebán, M. Meraz, K. Wrobel, and A. Tomasini, “Pentachlorophenol sorption in nylon fiber and removal by immobilized Rhizopus oryzae ENHE,” Journal of Hazardous Materials, vol. 190, no. 1–3, pp. 707–712, 2011. View at: Publisher Site  Google Scholar
 M. A. Abdullah, L. Chiang, and M. Nadeem, “Comparative evaluation of adsorption kinetics and isotherms of a natural product removal by Amberlite polymeric adsorbents,” Chemical Engineering Journal, vol. 146, no. 3, pp. 370–376, 2009. View at: Publisher Site  Google Scholar
 R. Lancelot, M. Lesnoff, and J. J. McDermott, “Use of Akaike information criteria for model selection and inference. An application to assess prevention of gastrointestinal parasitism and respiratory mortality of Guinean goats in Kolda, Senegal,” Preventive Veterinary Medicine, vol. 55, no. 4, pp. 217–240, 2002. View at: Publisher Site  Google Scholar
 J. Ødegård, J. Jensen, G. Klemetsdal, P. Madsen, and B. Heringstad, “Genetic analysis of somatic cell score in Norwegian cattle using random regression testday models,” Journal of Dairy Science, vol. 86, no. 12, pp. 4103–4114, 2003. View at: Google Scholar
 C. P. McDonald and N. R. Urban, “Using a model selection criterion to identify appropriate complexity in aquatic biogeochemical models,” Ecological Modelling, vol. 221, no. 3, pp. 428–432, 2010. View at: Publisher Site  Google Scholar
 A. Vaia and N. V. Sahinidis, “Simultaneous parameter estimation and model structure determination in FTIR spectroscopy by global MINLP optimization,” Computers and Chemical Engineering, vol. 27, no. 6, pp. 763–779, 2003. View at: Publisher Site  Google Scholar
 Y. Li, A. M. Rauth, and X. Y. Wu, “Prediction of kinetics of doxorubicin release from sulfopropyl dextran ionexchange microspheres using artificial neural networks,” European Journal of Pharmaceutical Sciences, vol. 24, no. 5, pp. 401–410, 2005. View at: Publisher Site  Google Scholar
 A. A. M. Daifullah and B. S. Girgis, “Removal of some substituted phenols by activated carbon obtained from agricultural waste,” Water Research, vol. 32, no. 4, pp. 1169–1177, 1998. View at: Publisher Site  Google Scholar
 S. J. Allen, Q. Gan, R. Matthews, and P. A. Johnson, “Comparison of optimised isotherm models for basic dye adsorption by kudzu,” Bioresource Technology, vol. 88, no. 2, pp. 143–152, 2003. View at: Publisher Site  Google Scholar
 M. C. Ncibi, “Applicability of some statistical tools to predict optimum adsorption isotherm after linear and nonlinear regression analysis,” Journal of Hazardous Materials, vol. 153, no. 12, pp. 207–212, 2008. View at: Publisher Site  Google Scholar
 I. Langmuir, “The constitution and fundamental properties of solids and liquids. Part I. Solids,” The Journal of the American Chemical Society, vol. 38, no. 2, pp. 2221–2295, 1916. View at: Google Scholar
 H. Freundlich, “Über die adsorption in iösungen,” The Journal of Physical Chemistry, vol. 57, pp. 385–470, 1907. View at: Google Scholar
 O. Redlich and D. L. Peterson, “A useful adsorption isotherm,” The Journal of Physical Chemistry, vol. 63, pp. 1024–1026, 1959. View at: Google Scholar
 J. Wang, C. P. Huang, H. E. Allen, D. K. Cha, and D. W. Kim, “Adsorption characteristics of dye onto sludge particulates,” Journal of Colloid and Interface Science, vol. 208, no. 2, pp. 518–528, 1998. View at: Publisher Site  Google Scholar
 W. Fritz and E. U. Schluender, “Simultaneous adsorption equilibria of organic solutes in dilute aqueous solutions on activated carbon,” Chemical Engineering Science, vol. 29, no. 5, pp. 1279–1282, 1974. View at: Google Scholar
 M. I. ElKhaiary and G. F. Malash, “Common data analysis errors in batch adsorption studies,” Hydrometallurgy, vol. 105, no. 34, pp. 314–320, 2011. View at: Publisher Site  Google Scholar
 A. Gunay, “Application of nonlinear regression analysis for ammonium exchange by natural (Bigadiç) clinoptilolite,” Journal of Hazardous Materials, vol. 148, no. 3, pp. 708–713, 2007. View at: Publisher Site  Google Scholar
 M. Hadi, M. R. Samarghandi, and G. McKay, “Equilibrium twoparameter isotherms of acid dyes sorption by activated carbons: study of residual errors,” Chemical Engineering Journal, vol. 160, no. 2, pp. 408–416, 2010. View at: Publisher Site  Google Scholar
 A. Kapoor and R. T. Yang, “Correlation of equilibrium adsorption data of condensible vapours on porous adsorbents,” Gas Separation and Purification, vol. 3, no. 4, pp. 187–192, 1989. View at: Google Scholar
 D. W. Marquardt, “An algorithm for least squares estimation of nonlinear parameters,” Journal of the Society for Industrial and Applied Mathematics, vol. 11, pp. 431–441, 1963. View at: Google Scholar
 F. M. Mutua, “The use of the Akaike information criterion in the identification of an optimum flood frequency model,” Hydrological Sciences Journal, vol. 39, no. 3, pp. 235–244, 1994. View at: Google Scholar
 J. F. Porter, G. McKay, and K. H. Choy, “The prediction of sorption from a binary mixture of acidic dyes using single and mixedisotherm variants of the ideal adsorbed solute theory,” Chemical Engineering Science, vol. 54, no. 24, pp. 5863–5885, 1999. View at: Publisher Site  Google Scholar
 A. Seidel and D. Gelbin, “On applying the ideal adsorbed solution theory to multicomponent adsorption equilibria of dissolved organic components on activated carbon,” Chemical Engineering Science, vol. 43, no. 1, pp. 79–88, 1988. View at: Google Scholar
 S. Chatterjee, S. K. Das, R. Chakravarty, A. Chakrabarti, S. Ghosh, and A. K. Guha, “Interaction of malathion, an organophosphorus pesticide with Rhizopus oryzae biomass,” Journal of Hazardous Materials, vol. 174, no. 1–3, pp. 47–53, 2010. View at: Publisher Site  Google Scholar
 Y. Yang, Z. Li, G. Wang et al., “Computational identification and analysis of the key biosorbent characteristics for the biosorption process of Reactive Black 5 onto fungal biomass,” PloS ONE, vol. 7, no. 2, pp. 1–8, 2012. View at: Google Scholar
 T. Şişmanoǧlu, “Removal of some fungicides from aqueous solution by the biopolymer chitin,” Colloids and Surfaces A, vol. 297, no. 1–3, pp. 38–45, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Zvezdelina Lyubenova Yaneva 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.