PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2010; 5(11): e13948.
Published online 2010 November 11. doi:  10.1371/journal.pone.0013948
PMCID: PMC2978715

The Role of Lactic Acid Adsorption by Ion Exchange Chromatography

Richard C. Willson, Editor

Abstract

Background

The polyacrylic resin Amberlite IRA-67 is a promising adsorbent for lactic acid extraction from aqueous solution, but little systematic research has been devoted to the separation efficiency of lactic acid under different operating conditions.

Methodology/Principal Findings

In this paper, we investigated the effects of temperature, resin dose and lactic acid loading concentration on the adsorption of lactic acid by Amberlite IRA-67 in batch kinetic experiments. The obtained kinetic data followed the pseudo-second order model well and both the equilibrium and ultimate adsorption slightly decreased with the increase of the temperature at 293–323K and 42.5 g/liter lactic acid loading concentration. The adsorption was a chemically heterogeneous process with a mean free energy value of 12.18 kJ/mol. According to the Boyd_plot, the lactic acid uptake process was primarily found to be an intraparticle diffusion at a lower concentration (<50 g/liter) but a film diffusion at a higher concentration (>70 g/liter). The values of effective diffusion coefficient Di increased with temperature. By using our Equation (21), the negative values of ΔG° and ΔH° revealed that the adsorption process was spontaneous and exothermic. Moreover, the negative value of ΔS° reflected the decrease of solid-liquid interface randomness at the solid-liquid interface when adsorbing lactic acid on IRA-67.

Conclusions/Significance

With the weakly basic resin IRA-67, in situ product removal of lactic acid can be accomplished especially from an open and thermophilic fermentation system without sterilization.

Introduction

Lactic acid is a very important organic acid with a wide range of applications in the food, pharmaceutical, leather and textile industries. It is classified as GRAS for general-purpose food additives by the FDA of USA. Because of the properties such as mild acidic taste, nonvolatility, lack of odor and bacteriostasis, lactic acid is used as an acidulant, taste enhancer, pH regulator and preservative in the food industry. Sodium lactate, calcium lactate, ethyl lactate and the biodegradable polymer, polylactic acid (PLA), have already been used in pharmaceutical and cosmetic applications, especially PLA for medical applications such as controlled-release drugs, surgical sutures, and prostheses [1]. Lactic acid can be conventionally produced by chemical or enzymatic synthesis, but is mainly by fermentation using lactic acid bacteria, such as some Lactobacillus and thermophilic Bacillus species, or the fungus Rhizopus oryzae [2], [3].

To date, many attempts have been made to improve the fermentation productivity which suffers from serious product inhibition. CaCO3 is one of the traditional additives for product precipitation and neutralization. It is added into lactic acid fermentation broths to neutralize the free lactic acid and thus to minimize product inhibition. Although in downstream processing, the raw calcium lactate is concentrated, crystallized, separated from dissolved impurities by filtration, and acidified with sulfuric acid, some impurities are still in the cake and quite a lot of calcium lactate is washed away and becomes a pollutant to our environment.

Recently, several novel approaches have been investigated for lactic acid recovery, such as solvent extraction, electrodialysis, etc. However, due to the hydrophilic nature, lactic acid is hardly extractable by the common organic solvents. Although reactive extraction has been considered to be an interesting alternative to the conventional process, it needs high amounts of solvents and the toxic effects of the extractants and diluents restrict its application. With the development of the membrane process, electrodialysis fermentation is promising since it can remove the lactic acid from the broth and maintain the pH in a proper value. However, the approach also encounters many problems such as membrane fouling, deionization of the fermentation broth and a higher operating cost [4].

Ion exchange technique has been widely used in bio-separations. As for lactic acid adsorbents, they must possess the important characteristics of high capacity and selectivity for lactic acid over water and substrate, regenerability, and biocompatibility with microorganisms. Recently, many different adsorbents have been investigated for lactic acid removal from fermentation broth such as poly(4-vinylpyridine) resin (PVP), Reillex 425, MWA-1, VI-15, Amberlite IRA-400, IRA-420, IRA-900, IRA-92, IRA-35 and IRA-96, etc [5][8]. A few papers [8], [9], [10] have reported that IRA-67 resin is a useful adsorbent for lactic acid or some other organic acid extraction. However, little systematic study has been done on the mechanism and adsorption efficiency of lactic acid under different operating conditions.

The purpose of this work is to investigate the adsorption kinetics of Amberlite IRA-67 by examining temperature, resin dose and lactic acid concentration. Also, we examine the adsorption mechanisms at different temperatures and lactic acid loading concentrations. These data could be very useful for the in situ removal of lactic acid especially by thermophilic Bacillus producers in unsterilized media. Adsorption thermodynamics is also studied in this condition.

Results

Adsorption kinetics

Information on the kinetics of lactic acid uptake by IRA-67 resin was needed for selecting optimum operating conditions for in situ removal of lactic acid from the aqueous environment. Figure 1 showed the changes in adsorption of lactic acid at different contact times at 293, 303, 313 and 323K, respectively. As seen from the chart, the value of qt (g lactic acid/g wet resin) decreased somewhat with the increase of operating temperature, but the effect was not significant. Considering the possible bacterial contamination during the process of in situ removal of the product lactic acid, 323K was thought to be a promising operating temperature which would sharply decrease such contamination. As to the reversible reaction of An external file that holds a picture, illustration, etc.
Object name is pone.0013948.e001.jpg, the change of lactic acid concentration (An external file that holds a picture, illustration, etc.
Object name is pone.0013948.e002.jpg) is relative with (CLA−Ce)n and n is 1 or 2 [11].

Figure 1
Effects of contact time and temperature on adsorption of lactic acid by IRA-67.

When n  = 1, then:

equation image
(1)

After calculation, the following equation is obtained:

equation image
(2)

Substituting Equation (2) into Equation (26), then:

equation image
(3)

where R(s) stands for raw resin (R) in solid state (s), LA (aq) represents the free lactic acid (LA) in liquid state (aq), R_LA stands for resin occupied by lactic acid, C0 is the initial concentration of lactic acid (g/liter), CLA and Ce are the concentrations of lactic acid at time t (min) and equilibrium (g/liter) respectively, qe is the amount of lactic acid adsorbed at equilibrium (g/g wet resin), qt is the amount of lactic acid adsorbed at time t (min), and k1 is the equilibrium rate constant of the pseudo-first order adsorption.

The rate constants were obtained via the straight line plots of ln(qe−qt) vs t under different experimental conditions. The values of k1 and correlation coefficient An external file that holds a picture, illustration, etc.
Object name is pone.0013948.e006.jpg were calculated from these plots (Table 1).

Table 1
Kinetic parameters of pseudo-first order and pseudo-second order reaction at various temperatures.

When n  = 2, then:

equation image
(4)

After calculation, the following Equation (5) is derived:

equation image
(5)

Substituting Equation (5) into Equation (26), then:

equation image
(6)

Equation (6) can be rearranged as a linear form:

equation image
(7)

and

equation image
(8)

where k2 is the rate constant of the pseudo-second order equation, and h is the initial sorption rate (g/g wet resin per min) [12].

The k2, An external file that holds a picture, illustration, etc.
Object name is pone.0013948.e014.jpg and qe values were calculated by a plot of t/qt vs t and were given in Table 1. The linear plots of t/qt vs t were shown in Figure 2 for the pseudo-second order model for adsorption of lactic acid by IRA-67 at 293–323K.

Figure 2
Linear plots of t/qt vs t for the pseudo-second order model.

The results in Table 1 revealed that the pseudo-second order kinetic model provided a better correlation for adsorption of lactic acid by IRA-67 compared to the pseudo-first order model, since the An external file that holds a picture, illustration, etc.
Object name is pone.0013948.e015.jpg values (0.994–0.999) were much higher than the An external file that holds a picture, illustration, etc.
Object name is pone.0013948.e016.jpg (0.900–0.954) and the experimental qe,exp values were much closer to the theoretical qe2,cal values than to the theoretical qe1,cal values.

Adsorption isotherm models

Adsorption equilibrium is achieved when the solution concentration does not vary with the contact time. These equilibrium data can be applied to predict the adsorption models and related theories linked to the adsorption equilibrium. In the present study, the Langmuir, Freundlich and Dubinin-Radushkevich (D-R) isotherm models are employed and their expressions are given by Equations (9), (10) and (11), respectively [13], [14], [15].

The Langmuir theory assumes that adsorption takes place at specific homogeneous adsorbent sites and once a lactic acid molecule occupies a site, no further adsorption can occur at the same site. The linear form is given by:

equation image
(9)

where Ce is the lactic acid concentration at equilibrium (g/liter), qm is the maximum adsorption capacity of IRA-67 (g lactic acid/g wet resin), KD is equilibrium constant of the Langmuir isotherm model (g/liter). KD and qm were calculated to be 1.148 and 0.1923 (R2  = 0.950).

The Freundlich adsorption isotherm model assumes that the adsorbent is composed of a heterogeneous adsorption surface with different classes of adsorption sites and the linear form is expressed by:

equation image
(10)

where Kf is a constant related to the adsorption capacity and n is an empirical parameter linked to the adsorption intensity varying with the heterogeneity of the adsorbent. The values of Kf and n were 0.14 and 13.89 according to the plot of lnqe vs lnCe presented in Figure 3 (R2  = 0.996).

Figure 3
The Freundlich isotherm linear plot for adsorption of lactic acid.

The D-R isotherm model is more frequently used to determine whether the nature of adsorption process is physical or chemical, when the equilibrium data are submitted. The linear equation of the D-R isotherm model is represented as:

equation image
(11)

where qE is the equilibrium concentration of lactic acid adsorbed on IRA-67 (mol/g), qM is the maximum lactic acid adsorption capacity (mol/g), E is the adsorption mean free energy (J/mol) and ε represents the Polanyi potential (ε  = RTln(1 + 1/Ce)). qM and E were calculated to be 2.140*10−3 mol/g and 12.18 kJ/mol (R2  = 0.940) by the lnqE vs ε2 plot.

The results from the D-R isotherm model revealed that the adsorption of lactic acid was a chemical reaction process, since the obtained value of adsorption mean free energy E lay within 8–16 kJ/mol in this case [16]. Although both the Langmuir and D-R isotherm models were to some extent able to describe the relationship between the amount of equilibrium adsorption and equilibrium concentration of lactic acid, via contrasting the value of correlation coefficient R2, the Freundlich adsorption isotherm model was found to perfectly simulate the relationship between qe and Ce and the n value (n[dbl greater-than sign]1) indicated that the adsorption of lactic acid by IRA-67 was favorable under the studied condition.

Effects of resin dose and lactic acid concentration

The effects of IRA-67 dose and lactic acid concentration were studied. Figures 4 and and55 showed a series of contact time curves with wet resin dose varying from 0.5 to 1.75 g/10 ml. Similar plots were also obtained with lactic acid concentrations from 21.5 to 99.0 g/liter. Table 2 exhibited the results as a plot of t/qt against t for adsorption of lactic acid for the pseudo-second order model with the correlation parameters of pseudo-second order rate, as well as the adsorption rate constant k2, initial adsorption rate h, and the equilibrium adsorption capacity q2 as functions of resin dose and lactic acid concentration. Similar results of pseudo-first order model were also obtained (Table S1). k2 and h for IRA-67 showed an increase but q2 decreased with an increase in the above resin dose. On the other hand, k2, h and q2 showed an increase with the lactic acid loading concentrations from 21.5 to 99.0 g/liter.

Figure 4
Effect of resin dose on adsorption of lactic acid.
Figure 5
Pseudo-second order kinetics for adsorption of lactic acid at various resin doses.
Table 2
Kinetic parameters of pseudo-second order reaction at various resin doses and lactic acid concentrations.

The corresponding linear curves of the values of q2 and h can be respectively depicted as a function of ms or C0 for IRA-67 by [12]:

equation image
(12)
equation image
(13)
equation image
(14)
equation image
(15)

substituting the values of q2 and h from Table 2 respectively into Equations (12), (13), (14) and (15) to obtain the empirical parameters listed in Table 3 and then into Equations (7) and (8), respectively. Finally, the relationships of qt against ms, t and qt against C0, t are expressed as follows:

equation image
(16)
equation image
(17)
Table 3
Empirical parameters for predicted q2 and h from ms and C0.

These equations can be used to predict the amount of lactic acid adsorption at any reaction time at the given resin dose and lactic acid concentration. The three dimensional plots of Equations (16) and (17) were shown in Figures 6 and and77.

Figure 6
Effect of resin dose on adsorption of lactic acid onto resin at different contact times.
Figure 7
Effect of lactic acid concentration on its adsorption onto resin at different contact times.

Adsorption mechanisms

In the adsorption process, either the film diffusion or the intraparticle diffusion will control the overall rate of adsorption. It is very useful to predict the rate-limiting step of the adsorption process by studying the resin application. The dynamic data have been analyzed by the model given by Boyd [17] as expressed by the following equations:

equation image
(18)
equation image
(19)
equation image
(20)

where F is the fractional attainment of equilibrium at time t (min), Di (cm2/s) is the effective diffusion coefficient of lactic acid in the resin phase, and r (cm) is the radius of the resin particle.

The values of F can be calculated using Equation (20). A value of Bt corresponds with each value of F by Equation (18) and the values of Di can be obtained by Equation (19). The control mechanism of either film diffusion or intraparticle diffusion for lactic acid adsorption by IRA-67 can be determined by the linear plots of Bt vs t. The results in Table 4 demonstrated that at the lower lactic acid concentrations (<50 g/liter), the adsorption was an intraparticle diffusion, since the linear plots passed through the origin. However, at the higher concentrations (>70 g/liter), the linear plots did not pass through the origin, which indicated that the adsorption was a film diffusion. Moreover, the value of F had nothing with the lactic acid concentration and was only dependent on the Di/r2 ratio. In addition, the values of Di changed very little with the increase of lactic acid concentration but increased with the rise of temperature. The reason is that the increased mobility of the ingoing resin at the higher temperature surpassed the effluence of retarding forces [18].

Table 4
Parameters of Boyd model at various lactic acid concentrations and temperatures.

Adsorption thermodynamics

Figure 8 showed the change of adsorption at different temperatures. According to the chart, the amount of adsorption decreased with the increase of temperature, but the variation was not large. In order to study the thermodynamic behavior of lactic acid adsorption by IRA-67, thermodynamic parameters such as changes in free energy (ΔG°), enthalpy (ΔH°) and entropy (ΔS°) were applied. For the reaction R↔RLA, based on the theory described by Beard and Qian [19], we deduced for the first time the formula for calculating the change of Gibbs free energy (ΔG) as follows:

equation image
(21)

where Na is the total number of effective binding sites at different temperature and lactic acid loading concentration, and Nb is the number of effective binding sites occupied by lactic acid at different contact times. In this equation, the values of Na and Nb are affected by both the operating temperature and lactic acid concentration. However, in the following calculations, it is found that the value of Nb/(Na−Nb) is only linked to temperature. When the reaction reaches the equilibrium state at ΔG  = 0, then ΔG°, ΔH° and ΔS° are calculated by the following equations:

equation image
(22)
equation image
(23)
equation image
(24)
equation image
(25)

where Kd is the distribution coefficient for the adsorption, Ne is the number of effective binding sites occupied by lactic acid in equilibrium state, T is the solution temperature (K) and R is universal gas constant of 8.314 J/(mol•K).

Figure 8
Effect of temperature on adsorptive capacity of IRA-67.

The distribution coefficient Kd could be also obtained by plotting ln(Ce/qe) vs Ce and elongating to zero Ce as suggested by Khan and Singh [20]. However, in the following calculation with such Kd values, the R2 value of lnKd vs 1/T plot is lower than when using Equation (21). Perhaps the main reason is that it is very hard to obtain a high value of R2 from plots of ln(Ce/qe) vs Ce. Since this indicates the capability of the resin to retain lactic acid, Kd would be a useful parameter to compare the adsorptive capacities of different resins for lactic acid, if measured under same experimental conditions [21].

ΔG° can be determined by Equation (22), ΔH° and ΔS° were respectively calculated from the slope and intercept of van't Hoff plots of lnKd vs 1/T as presented in Figure 9. The change in Gibbs free energy (ΔG°) for adsorption of lactic acid onto IRA-67 was available with an initial concentration of 42.5 g/liter lactic acid and the results were given in Table 5. The obtained negative values of ΔG° indicated that the adsorption of lactic acid onto IRA-67 was spontaneous. Furthermore, the negative ΔG° descending with the increase of temperature indicated a decrease in adsorption at higher temperatures [22]. The negative value of enthalpy ΔH°, calculated to be −32.82 kJ/mol, indicated that the adsorption processes was exothermic in nature at 293–323K. The obtained negative value of ΔS° exhibited the decreased randomness at the solid-solution interface during the adsorption of lactic acid on the active sites of IRA-67. It was also noted that, although decreasing against the rising temperature, the values of Kd were still large at a relatively high temperature, which supported the possibility of industrial lactic acid adsorption by IRA-67.

Figure 9
Plot of lnKd versus 1/T at lactic acid concentration of 42.5 g/liter.
Table 5
Thermodynamic parameters for adsorption of lactic acid onto Amberlite IRA-67 at different temperatures.

Discussion

We have described the dynamic equilibrium and thermodynamic evaluation on adsorption of lactic acid by the weakly basic resin IRA-67. The pseudo-second order model fits well with the kinetic data, in which the equilibrium and ultimate adsorption are 0.182 g and 0.199 g at 323K, respectively. Moreover, the Freundlich isotherm model fits the equilibrium data well, which indicates that the adsorption process is heterogeneous. Therefore, the lactic acid molecule is absorbed first on the higher energy site prior to the lower energy site. The resin dose and lactic acid loading concentration respectively somehow exhibit negative and positive effects on unit resin adsorption, and it is possible to predict the adsorption amount of lactic acid at any time with various resin doses and lactic acid loading concentrations.

It is first found that the overall rate of adsorption is mainly controlled by the intraparticle diffusion at a lower lactic acid concentration (<50 g/liter) but is controlled by the film diffusion at a relatively higher concentration (>70 g/liter) via a slightly increased Di with the increase of temperature. Thus, the effect of temperature on lactic acid adsorption is very limited. In this research, Equation (21) could be favorably matched the formula in reference [20] and provides an easier alternative way to determine resin adsorption thermodynamics. The adsorption process for lactic acid is spontaneous and exothermic. Also, the values of the distribution coefficient Kd changes inversely with temperature resulting in a slightly decreased amount of lactic acid adsorbed on the resin. However, the overall value is still large at a relatively high temperature due to the fact that the distribution coefficient demonstrates the capability of the resin to retain lactic acid, and the higher the value of Di, the lower is the adsorption of lactic acid.

To date, the simultaneous removal of lactic acid in the fermentation has not been carried out as an industrial application. The fatal limitation is that during the in situ product removal process of lactic acid, it is fairly difficult to avoid the risk of contamination caused by other bacteria. The reason is that the optimal fermentation temperature for most industrial strains is in the range of 303K to 315K. To achieve continued lactic acid recovery on the industrial scale, one of the most likely ways is to use thermophilic producers and to efficiently adsorb the lactic acid by resin at high temperature, thus dramatically reducing the possibility of contamination. Finally, the results of this work indicate that, due to the characteristics of quick and high adsorption especially at a high temperature, use of IRA-67 is possible for in situ product removal of lactic acid. This, coupled with the use of thermophilic producer, would allow low-carbon clean manufacture without contamination by other bacteria. Thus, the production cost would sharply decrease and environmental pollution would be avoided.

Materials and Methods

Materials

Amberlite IRA-67 (Rohm and Hass Company, USA), a weakly basic gel-type polyacrylic resin with a functional tertiary amine group, was selected for this study. L-lactic acid (purity>98%, Sigma-Aldrich Corporation, USA) was used as HPLC standard and a conventional lactic acid reagent of AR grade (purity  = 85–90%, Tianjin Guangfu Fine Chemical Institute, China) was used for sorption experiments.

Adsorption kinetics and equilibrium studies

Due to the resin's nature, Amberlite IRA-67 hardly occurs in the OH form. Before utilization, the resin was treated with 4% NaOH for 24 h, then washed with double-distilled water. One gram of wet Amberlite IRA-67 was mixed with 10 ml of 42.5 g/liter lactic acid solution and shaken at 150 rpm in a 100-ml conical flask. Duplicate samples were taken at various time intervals at 293, 303, 313 and 323K. The equilibrium was achieved at 24 h. The amount of lactic acid adsorbed onto the resin at equilibrium was calculated from the mass balance of the equation as given below:

equation image
(26)

where C0 and Ce are the initial and equilibrium concentrations of lactic acid solution (g/liter), qe is the amount of lactic acid adsorbed at equilibrium (g/g wet resin), V is the volume of the lactic acid solution (liter), and W is the mass of the resin used (g).

Desorption process

Different characteristics of resins require different elution protocols for desorption. For IRA-67 resin (free base form), NaOH (2–4%), NH4OH (2–4%) and Na2CO3 (4–8%) are the commonly used regenerants. In this study, 4% NaOH was used to completely desorb the adsorbed lactic acid with a contact time of 24 h.

Analytical methods

Lactic acid concentration was determined using an Agilent 1100 series high-performance liquid chromatography (HPLC) system (Hewlett-Packard Corporation, USA), a VWD UV detector, and a Bio-Rad Aminex HPX-87H column maintained with a thermostat at 323K. The mobile phase was 5 mM sulfuric acid at a flow rate of 0.60 ml/min.

Supporting Information

Table S1

Kinetic parameters of pseudo-first-order reaction at various resin doses and lactic acid concentrations.

(0.03 MB DOC)

Acknowledgments

The authors thank Dr. Arnold L. Demain (R.I.S.E., Drew University, USA) for critically reading of the manuscript, and thank Dexin Han and Binzhe Wang for helping with the lactic acid adsorption experiment.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: The work for this manuscript is financially supported by National Basic Research Program of China (2007CB714305 and 2011CB707401) (973 Program: http://www.most.gov.cn/eng/programmes1/200610/t20061009_36223.htm), National High-Tech R & D Program of China (2008AA10Z336) (863 Program: http://www.most.gov.cn/eng/programmes1/200610/t20061009_36225.htm) and Key Project of Tianjin Natural Science Foundation of China (08JCZDJC15100) (http://www.tstc.gov.cn/). The funders had no role in study design, data collection and analysis, decision to publish, or preparation for this manuscript.

References

1. El-Beyrouty C, Huang V, Darnold CJ, Clay PG. Poly-L-lactic acid for facial lipoatrophy in HIV. Ann Pharmacother. 2006;9:1602–1606. [PubMed]
2. Kasuga T, Ota Y, Nogami M, Abe Y. Preparation and mechanical properties of poly lactic acid composites containing hydroxyapatite fibers. Biomaterials. 2001;22:19–23. [PubMed]
3. Qin J, Zhao B, Wang X, Wang L, Yu B, et al. Non-sterilized fermentative production of polymer-grade L-lactic acid by a newly isolated thermophilic strain Bacillus sp. 2–6. PLoS ONE. 2009;4:1–7. [PMC free article] [PubMed]
4. Wang E, Hatanaka H, Iijima S, Tokebayashi T, Shi Z, et al. Control of cell and lactate concentration in a hollow fiber bioreactor for lactic acid fermentation. J Chem Eng Jpn. 1988;21:36–40.
5. Chen CC, Ju LK. Adsorption characteristics of polyvinylpyridine and activated carbon for lactic acid recovery from fermentation of Lactobacillus delbrueckii. Sep Sci Technol. 1998;33:1423–1437.
6. Chabani M, Amrane A, Bensmaili A. Equilibrium sorption isotherms for nitrate on resin Amberlite IRA 400. J Hazard Mater. 2009;165:27–33. [PubMed]
7. Srivastava A, Roychoudhury K, Sahai V. Extractive lactic acid fermentation using ion exchange resin. Biotechnol Bioeng. 1992;39:607–613. [PubMed]
8. Moldes AB, Alonso JL, Parajo JC. Recovery of lactic acid from simultaneous saccharification and fermentation media using anion exchange resins. Bioproc Biosyst Eng. 2003;25:357–363. [PubMed]
9. Uslu H, İnici İ, Bayazit ŞS. Adsorption equilibrium data for acetic acid and glycolic acid onto Amberlite IRA-67. J Chem Eng Data. 2010;55:1295–1299.
10. Patel M, Bassi A, Zhu JJX, Gomaa H. Investigation of a dual-particle liquid-solid circulating fluidized bed bioreactor for extractive fermentation of lactic acid. Biotechnol Progr. 2008;24:821–831. [PubMed]
11. Ho YS, McKay G. Pseudo-second order model for sorption process. Process Biochem. 1999;34:451–465.
12. Ho YS, McKay G. A kinetic study of dye sorption by biosorbent waste product pith. Resour Conserv Recy. 1999;25:171–193.
13. Langmuir I. The adsorption of gases on plane surfaces of glass, mica and platinum. J Am Chem Soc. 1918;40:1361–1403.
14. Freundlich H. Adsorption in solution. Z Phys Chem. 1906;40:1361–1368.
15. Dubinin MM, Zaverina ED, Radushkevich LV. Sorption and structure of active carbons. I. Adsorption of organic vapors. Zhurnal Fizicheskoi Khimii. 1947;21:1351–1362.
16. Helfferich F. Ion Exchange. New York, USA: McGraw Hill; 1962.
17. Boyd GE, Adamson AW, Myers LS., Jr The exchange adsorption of ions from aqueous solutions by organic zeolites. II. Kinetics. J Am Chem Soc. 1947;69:2836–2848. [PubMed]
18. Gupta VK, Ali I. Removal of DDD and DDE from wastewater using bagasse fly ash, a sugar industry waste. Wat Res. 2001;35:33–40. [PubMed]
19. Beard DA, Qian H. Relationship between thermodynamic driving force and one-way fluxes in reversible processes. PLoS ONE. 2007;2:1–4. [PMC free article] [PubMed]
20. Khan AA, Singh RP. Adsorption thermodynamics of carbofuran on Sn (IV) arsenosilicate in H+, Na+, and Ca2+ forms. Colloids Surf. 1987;24:33–42.
21. Fontes MPF, Gomes PC. Simultaneous competitive adsorption of heavy metals by the mineral matrix of tropical soils. Appl Geochem. 2003;18:795–804.
22. Sarl A, Tuzen M. Equilibrium, thermodynamic and kinetic studies on aluminum biosorption from aqueous solution by brown algae (Padina pavonica) biomass. J Hazard Mater. 2009;171:973–979. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science