 Open Access
 Total Downloads : 34
 Authors : P. Pallavi , P. King , Y. Prasanna Kumar
 Paper ID : IJERTV8IS080046
 Volume & Issue : Volume 08, Issue 08 (August 2019)
 Published (First Online): 16082019
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Biosorption of Malachite Green by Raphanus Raphanistrum LeavesExperimentation and Process Optimization
P. Pallavi 1*
1Ph.D Scholar, JNTUK
P. King 2, Y. Prasanna Kumar
2Department of Chemical Engineering, Andhra University, Visakhapatnam530003, (Andhra Pradesh) India
3Visakha Institute of Engg. & Techn. Visakhapatnam, India
Abstract The study was focused on identification and describing the binding mechanism of the dye by Raphanus Raphanistrum in a biosorption process optimized by Response Surface Methodology based on the Central Composite Design (CCD). The parameters influencing the biosorption capacity, such as contact time, solution pH, temperature, biosorbent dosage, average biosorbent size and initial dye concentration were investigated. Langmuir, Freundlich and Temkin sorption isotherms were applied to the biosorption equilibrium data. The biosorption isotherm data was found to fit the Freundlich model. Pseudofirst order, Pseudo second order, Elovchi models were used to describe the biosorption kinetic studies. Pseudosecond order kinetics gives the rate of biosorption and mass transfer kinetics. Thermodynamic analysis revealed that the biosorption process is spontaneous and exothermic in nature. Using Central Composite Design (CCD) the process parameters like contact time (min), solution pH, initial dye concentration (mg/L) and temperature (K) were analyzed. The secondorder polynomial equation, response surface plots were used to quantitatively determine the relationships between dependent and independent variables.
Keywords Malachite green, Raphanus Raphanistrum Leaves (RRL), Biosorption Isotherms, Kinetics, thermodynamic studies, Central Composite Design (CCD).
INTRODUCTION
Out of various industries, the dyeing industry is an important sector of the chemical industry. Dyes are substances capable of coloring fabrics to such an extent that the color cannot be removed by rubbing or washing. Dyes are widely used for the manufacture of textiles, paper, leather, fur, hair, food, drugs, cosmetics, wax, grease, petroleum products, and plastics [1]. India's dye industry produces many types of dyes and pigments. Production of dyestuff and pigments in India is close to 80,000 tones. India is the second largest exporter of dyestuffs and intermediates among the developing countries, after China. These dyes are being used in the textile industry at nearly 80%. However, a recent study conducted under the National Biodiversity strategy and Action Plan (BSAP) has revealed that chemical colors have all but wiped out India's wonderful vegetable dyes.
Color removal has become a major scientific interest, and hence significant attention has been directed towards the exploration of various biosorbent materials [2].
The Indian textile industries now predominantly use synthetic organic dyes like direct dyes; disperse dyes, vat dyes, acid dyes, reactive dyes, etc. The large variety of dyes and chemicals used in an attempt to make more attractive popular shades of fabrics for a competitive market render the effluent very complex. During the last decade, environmental issues associated with dyestuff production and their applications have grown significantly [3].
Malachite green (MG) dye has been used as a food coloring agent, food additive, medical disinfectant as well as a dye in silk, wool, jute, leather, cotton, paper and acrylic industries [4]. It is also extensively used as a topical fungicide [5] and echo parasiticide in aquaculture industries throughout the world [6]. The dye has been reported to cause serious health effects including effect on immune system and reproductive system as well as its genotoxic and carcinogenic potentials [6]. Researchers have been working on different ways of removing dyes from wastewater for many years. Dye wastewater is very difficult to treat, because of its synthetic origin and complex aromatic structure. Therefore, different procedures have been developed: adsorption onto materials such as activated carbon [3, 6] or waste products [7, 8], physical and chemical degradation [9, 10] in addition to a large number of other techniques; Fentons oxidation, electro coagulation [11], ozonization . Biosorption is the fast, economic, feasible and ecofriendly technique which also utilized the waste generated from industrial processing. In the present work, Raphanus Raphanistrum powder was tested for removal of Malachite green (MG) dye from aqueous solutions. The equilibrium, kinetics and thermodynamic data of the biosorption process were investigated. Information has been reported on the statistical optimization of experiment conditions for the biosorption of MG dye by RR biosorbent.
MATERIALS AND METHODS
Preparation of Raphanus Raphanistrum Leaf powder preparation procedure
Raphanus Raphanistrum leaves were collected from Vignans Engineering College, Vadlamudi.
The collected leaves were washed several times with deionized water until the wash water contains no dirt. These washed leaves were completely dried under sunlight for 30 days and were powdered using domestic mixer. The dried powder was sieved for 63212 Âµm.
Preparation of stock solution
Malachite green dye solution was prepared. A known weight of 0.1 g of dye was weighed and a standard stock solution of concentration 1000 mg/L was prepared in double distilled water and further working solutions of concentrations (25, 50, 75 and 100 mg/L) were prepared in stoppered bottles and stored.
EXPERIMENTATION
Batch biosorption studies were carried out at 303K on a rotary shaker at 180 rpm with time intervals of 5135 min. The biosorption studies were completed under effective process parameters of solution pH 210, dye concentrations 25100 mgL1, weight of the biosorbent 0.020.2 g and particle size of the biosorbent vary from 63 (242 mesh) 212 (72 mesh) m. The mixed biosorbent solutions were taken out, filtered and analyzed for Malachite green dye concentration in an UV Spectrophotometer at 618nm. By varying the initial concentration of the dye solution from 25100mg/L at equilibrium biosorption time=90 min, pH=6, dosage=0.1g and temperature=303K along with thermodynamic parameters which gives the feasibility of the process, 26 experiments were designed by the Design Expert Software with the optimized parameters.
The % biosorption is given as
% Biosorption = (CiCf)/Ci*100 (1)
The dye uptake onto RR is obtained by:
V (Ci C f )
surface. At 25 mg/L of initial dye concentration, the biosorption efficiency was increased from 82.64% to 97.76% with a contact time of 90 min.
Fig.1: Effect of contact time on the biosorption efficiency of MG dye onto RR biosorbent
Effect of solution pH
Biosorption Capacity is highly influenced by the solution pH. The degree of biosorption onto biosorbent surface is due to the surface charge on the biosorbent, which in turn influenced by the solution pH. Dissociation of functional groups present on the surface and subsequently shift in reaction kinetics and equilibrium characteristics takes place with change of solution pH. The effect of solution pH on biosorption efficiency of MG dye is shown in Fig.2. The biosorption efficiency increases with an increase in pH gradually and attains equilibrium at pH 6. At pH 6, the maximum % biosorption obtained was 97.76% for an initial dye concentration of 25 mg/L.
qt
1000w
(2)
Fig.2: Effect of solution pH on the biosorption efficiency of MG dye onto
Where qt is the amount of dye adsorbed on the RRL biosorbent surface (mg/g),
Ci is the initial concentration of olute in the solution before biosorption (mg/L),
Cf is the final concentration of solute in the solution after biosorption (mg/L),
V is the volume of the dye solution (ml) and
w is the weight of the Biosorbent (g).
RESULTS AND DISCUSSION
Effect of contact time:
The biosorption efficiency for different initial dye concentrations from 25 mg/L – 100 mg/L, with respect to contact time is represented in Fig.1. More active sites are available for biosorption at the initial stages of contact period. At 90 min, the %biosorption was constant even after increase in contact time, which indicates availability of less number of accessible vacant sites on the biosorbent
RR biosorbent
Effect of Initial dye concentration:
Dye removal is highly influenced by initial dye concentration.Fig.3. shows that with an increase in initial dye concentration the percentage biosorption was decreased and dye uptake was increased. At low concentration, most of the dye molecules adsorb on active sites and hence % biosorption was high at the initial dye concentration. The percentage of biosorption was decreased and dye uptake was increased with an increase in initial dye concentrations, increased number of dye molecules could not accommodate number of active sites for a fixed dosage of biosorbent. The dye uptake was increased from 12.22 to 47.34 mg/g with an increase in initial dye concentration from 25 to 100 mg/L at 303K.
Fig.3: Effect of initial dye concentration on the biosorption efficiency of MG dye onto RR biosorbent
Effect of Biosorbent dosage
The effect of biosorbent dosage on biosorption capacity of RR biosorbent is shown in Fig.4. With an increase of biosorbent dosage from 0.02 to 0.1g, the % biosorption increased gradually and attained a maximum % biosorption at 0.1g. The percentage biosorption increases proportionally with an increase of weight of the biosorbent dosage. The increase in biosorption with the increasing dosage of biosorbent is basically due to binding of almost all dye molecules and establishment of dynamic equilibrium between dye molecules on the surface and in aqueous solution. The % biosorption was increased from
81.84 to 97.76% with an increase of biosorbent dosage from 0.02 to 0.2 g for an initial concentration of 25 mg/L.
Fig.4: Effect of biosorbent dosage on the biosorption efficiency of MG dye onto RR biosorbent
Effect of average biosorbent size
The influence of average biosorbent size on biosorption with different initial dye concentrations is shown in Fig.5. The percentage biosorption decreased from 98.8% to 79.6% with an increase of particle size from 63Âµm to 212
Âµm for an initial concentration of 25mg/L. When compared with larger particles, the smaller biosorbent particles will have more specific surface area, where rapid biosorption process takes place with high mass transfer.
Fig.5: Effect of average biosorbent size on the biosorption efficiency of MG dye onto RR biosorbent
Effect of temperature
The effect of temperature on biosorption for different initial dye concentrations is shown in Fig.6. As the temperature increases the percentage biosorption decreases gradually, due to physicochemical bonds present between dye molecules and surface of biosorbent may be weakened resulting in the low percentage biosorption. Moreover, at high temperatures the solubility of dyes in aqueous solution increases, thus attractive forces between dye molecules and solution were stronger than attractive forces between dye molecules and surface. At a contact time of 90 min, the maximum biosorption was observed as 97.76% at 303K for an initial concentration of 25mg/L.
Fig.6: Effect of temperature on the biosorption efficiency of MG dye onto RR biosorbent
Biosorption Isotherms
The applicability of Langmuir isotherm model indicates the homogeneity with finite number of active sites on the surface of biosorbent. The Langmuir model also suggests the possibility of uniform coverage of biosorbent surface with dye molecules. The equation applied to analyze the Langmuir constants are:
(3)
where Ce is the concentration of the dye in the solution at equilibrium (mg/L), qe is the quantity of dye adsorbed at equilibrium (mg/g), qmax is the maximum biosorption capacity of the dye (mg/g) and KL is Langmuir constant (L/mg). A plot of Ce vs. Ce/qe was plotted. The values of constant and maximum biosorption are evaluated from equation (3). The high value of biosorption capacity (qmax) indicates the availability of high surface area.
Equation (4) gives the separation factor (RL) values represented in Table1. Ci is the initial concentration
(mg/L). The values of RL which are between 0 and 1 indicate that biosorption is favourable with an increase in the initial dye concentration.
0.5
0.4
0.5
0.4
RL=1/ (1+KLCi) (4)
y (303K) = 0.014x + 0.0412; RÂ² = 0.9807 y(313K)= 0.0139x + 0.0654; RÂ² = 0.9964 y(323 K)= 0.0107x + 0.1253; RÂ² = 0.9343
y (333K)= 0.0073x + 0.1925; RÂ² = 0.978
y (303K) = 0.5962x + 2.8729; RÂ² = 0.9968 y (313K) = 0.628x + 2.5704; RÂ² = 0.9757 y (323K)= 0.7101x + 2.1436; RÂ² = 0.9982 y(333K)= 0.8194x + 1.7379; RÂ² = 0.9984 y (343K) = 0.8729x + 1.3044; RÂ² = 0.9996
4.0
3.8
3.6
3.4
4.0
3.8
3.6
3.4
ln qe
ln qe
3.2
3.0
2.8
2.6
2.4
2.2
y(343K)= 0.0055x + 0.3031; RÂ² = 0.9923
1.0 0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0
ln C
e
Ce/qe (g/L)
Ce/qe (g/L)
Fig.8: Freundlich isotherms at different temperatures for biosorption of
0.3 MG dye onto RR biosorbent
0.2
0.1
0.0
0 2 4 6 8 10 12 14 16 18
Ce (mg/L)
Temkin isotherm equation is given as
qe = BT ln (AT) + BT ln (Ce) (6)
The biosorption potential and the biosorption energy variations in biosorption process are studied under Temkin isotherm. qe and Ce are the quantity of dye adsorbed
Fig.7: Langmuir isotherm at different temperatures for biosorption of MG dye onto RR biosorbent
For a heterogeneous surface with a nonuniform multilayer coverage of biosorbent surface with dye molecules, the Freundlich isotherm is derived. The Freundlich isotherm model is given by
(mg/g) and concentration of the dye in the solution at equilibrium (mg/L) respectively. Equation (6) gives the values of temkin constants BT and AT. A graph of ln Ce vs. qe was plotted and the values are tabulated in Table1. Furthermore, Freundlich isotherm was best fitted with data than Temkin isotherm.
50
40
50
40
y (303K) = 15.169x + 19.503; RÂ² = 0.976
ln qe ln K f

1 ln C n f
(5)
e
e
y (313K) = 15.878x + 11.771; RÂ² = 0.9971 y (323K) = 17.006x + 2.3373;RÂ² = 0.9535 y (333K)= 19.245x – 6.7595; RÂ² = 0.9709
qe and Ce are the quantity of dye adsorbed at equilibrium (mg/g) and concentration of the dye in the solution at
equilibrium (mg/L) respectively. A plot of lnCe vs. lnqe is plotted in Fig.8. Freundlich constant values of Kf and 1/nf were calculated by using equation (5). With an increase of temperature, the values of Kf decreased which indicates the exothermic nature of biosorption. The nf values are greater than one indicates the process was favourable and due to chemical interactions between biosorbate and biosorbent. The constant values are indicated in Table1. The R2 value represents Freundlich isotherm best suits for the data.
y(343K) = 19.516x – 15.01; RÂ² = 0.9688 30
qe (mg/L)
qe (mg/L)
20
10
0
1.0 0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0
ln Ce
Fig.9: Temkin isotherms at different temperatures for biosorption of MG dye onto RR biosorbent
Isotherm Model
Parameters
Temperature (K)
303
313
323
333 343
Langmuir
qmax (mg/g)
71.42
71.9
93.4
136.9
181.8
KL (L/g)
0.33
0.21
0.08
0.03
0.018
R2
0.98
0.99
0.93
0.97
0.9923
RL=1/(1+KLCo)
0.10
0.15
0.32
0.519
0.689
Freundlich
Kf (mg/g)/(L/g)n
17.68
13.07
8.53
5.68
3.68
nf (L/g)
1.67
1.59
1.40
1.22
1.14
R2
0.9
0.97
0.99
0.998
0.99
Temkin
AT (L/g)
3.61
2.09
1.14
0.703
0.463
BT(L/mg)
166.07
163.8
157.9
143.8
146.12
R2
0.97
0.99
0.95
0.970
0.968
Kinetic Modeling
Studies of biosorption kinetics provide the necessary information to evaluate the biosorption rate and to identify other factors that influence the biosorption process which in turn are used to optimize the efficiency of the biosorption process. The kinetic studies were conducted using pseudofirstorder, pseudosecondorder and Elovich kinetic models to best fit the data. The pseudofirstorder model generally expressed as
that for higher concentrations the Elovich model was not suitable to fit the data.
ln(qe qt ) k f t ln qe
(7)
Where kf (min1) is the firstorder biosorption process rate constant and qt (mg/g) is the weight of MG dye adsorbed on biosorbent surface at time t. Fig.10. shows a graph of t verses ln (qe qt) at different concentrations (25 100 mg/L) using the dye at optimum conditions of pH 6 and temperature 303K. The values of Pseudofirstorder kinetic model are shown in Table2. The values of R2 are low which suggest the nonapplicability of pseudofirst order kinetic model.
Pseudo secondorder kinetic model is expressed as
Fig.10: Pseudo First order Kinetic model for biosorption of MG dye onto RR biosorbent
k q
k q
2
2
t 1 t 1
qt qe s e
(8)
where ks (g mg1 min1) is the secondorder biosorption process rate constant. A plot of t vs. t/qt was plotted as shown in Fig.11 represents the secondorder kinetic model. The constants are evaluated from equation (8). The increase of R2 value suggested the applicability of pseudo second order kinetics model.
The applicability of second order model represents the presence of both physisorption and chemisorption process. The overall biosorption process can be limited by chemisorption. The pseudo second order was fitted with the data well.
Elovich model gives the information about the physisorption or chemisorption mechanism. Elovich model is used to describe the second order kinetics of chemisorptions and it is based on the assumption that the adsorption surface is energetically heterogeneous. The Elovich model is expressed as
Fig: 11. Pseudo Second order Kinetic model for biosorption of MG dye onto RR biosorbent
q 1 ln 1 ln t (9)
Where is initial rate and representing activation energy (g/mg). A graph of ln t vs. q is drawn where the constants are evaluated from equation (9). The constants and correlation coefficients are given in Table2, indicate
Fig: 12. Elovich kinetic model for the biosorption of MG dye onto RR biosorbent
Table.3: Kinetic rate constants for MG dye biosorption onto RR
Ci mg/L
Pseudofirstorder
PseudoSecondorder
Elovich Model
kf
qe
R2
ks
qe
R2
R2
25
1.1
2.2
0.9
0.12
11.7
0.9
2.5
2.525×10
10
0.8
50
1.0
4.8
0.9
0.05
23.0
0.9
1.1
1.089×10
9
0.8
75
0.9
9.5
0.9
0.02
34.3
0.9
0.5
2.046×10
6
0.9
100
1.1
14.5
0.9
0.01
44.4
0.9
0.4
2.863×10
6
0.8
Thermodynamic studies
Thermodynamic parameters give the feasibility of biosorption process. Thermodynamic studies provide information about the effect of temperature on biosorption process, which is essential for the estimation of energy requirement. The thermodynamic parameters are represented in Vant Hoff equation as
ln KD= So/R Ho/RT (10)
Where Ho is standard enthalpy, So is standard entropy, and Go is standard Gibbs free energy, R is the universal gas constant (8.314 J/mol K), T is the absolute temperature, K. Go represents the feasibility and spontaneous nature of the biosorption process. The enthalpy (Ho) and entropy (S0) are calculated from the equation (10) Go and KD are related by the following equation,
negative value of So indicates an increase in the randomness at the solidsolute interface during biosorption. With an increase of temperature, the feasibility of biosorption decreased with simultaneous decrease in negative values of Go.
Go RT ln K D (11)
A plot was plotted between 1/T vs. ln KD, shown in Fig.13. for different initial dye concentrations. The values of Ho, So and Go are represented in Table3. The negative values of Ho represents the reaction is exothermic. The
Fig.13: Thermodynamic studies at different concentrations for biosorption of MG dye onto RR biosorbent
Table.3. Thermodynamic energy parameters for the biosorption of MG dye using RR biosorbent
Ci
(mg/L)
Ho
(kJ/m ol)
So
(kJ/mol.
K)
Go(kJ/mol)
303K
313K
323K
333K
343K
25
41.66
0.11
8.33
7.23
6.13
5.03
3.93
50
36.85
0.09
9.58
8.68
7.78
6.88
5.98
75
30.06
0.07
8.85
8.15
7.45
6.75
6.05
100
25.32
0.06
7.14
6.54
5.94
5.34
4.74
Optimization using RESPONSE SURFACE METHODOLOGY (RSM)
A secondorder polynomial equation was developed from Table6 as regression equation:
2
2
% Removal (Y) = – 5266.07+11.12X1+13.64X2 4.17X3+35.52X40.05X 20.31X 20.7X 2
1 2 3
Response surface methodology (RSM) is a collection of
mathematical and statistical technique that is useful to model and analyze problems in which a response of interest is influenced by several variables. The independent and dependent variables in biosorption process are optimized using the RSM modeling technique. The most
0.05X4 +0.13X1X2+0.02X1X30.01X1X40.01X2X3
0.07X2X4+0.02X3X4
Where X1, X2, X3 and X4 are the code values for the independent variables (contact time, solution pH, initial dye concentration, temperature), X1X2, X1X3, X1X4, X2X3,
X2X4, X3X4, X 2,X 2, X 2 and X 2 are the significant model
1 2 3 4
important effects, interactions and analysis of variance (ANOVA) are performed using Design Expert software.
The four parameters opted for the process optimization: contact time (min), solution pH, initial dye concentration (mg/L), temperature (K).
k
k
By using the least square regression methodology, a second order multiple regression analysis equation along with process parameters are obtained.
terms for the biosorption of MG dye. A plot of experimental values versus predicted values is plotted in Fig.14. The interaction between the variables is represented by the response surface plots. (Fig.15 (ivi)) The surface plot (Fig. 15 (i) to 15(vi)) had a clear peak, where the optimum conditions falls inside the design boundary.
Prediction of optimal values
By using RSM, the maximum biosorption efficiency
Y o
k
k

i xi
i 1
ii i 1
k
k
ij
ij
x 2 i
i j
xi x j
obtained was 98% under the optimal conditions of t=88.8
min, pH= 5.9, Co =23.7mg/L, and T=301.2K. The biosorption efficiency was 97.76% under experimental process conditions.
Where Y is the % removal, o is the constant, i is the incline or straight impact of the data element xi, ii is the quadratic effect of input factor xi, ij is the linear by linear interaction effect between the input factor xi and is the residual term
Table 4. Experimental variables and their levels.
Variable 
Process parameter 
Range of Process parameters 

2 
1 
0 
1 
2 

X1 
Contact Time (min) 
80 
85 
90 
95 
100 
X2 
Solution pH 
2 
4 
6 
8 
10 
X3 
Initial dye concentration (mg/L) 
15 
20 
25 
30 
35 
X4 
Temperature (K) 
293 
298 
303 
308 
313 
Process parameters 
Initial observed 
Critical values 
Final observed limits 
Contact time (min) 
80 
88.81 
100 
Solution pH 
2 
5.96 
10 
Initial effluent concentration (mg/L) 
15 
23.77 
35 
Temperature (K) 
293 
301.26 
313 
Fig 14. values for the biosorption of MG dye.
Surface Plots
CONCLUSIONS

The present study confirms that Raphanus Raphanistrum as an effective biosorbent for the removal of Malachite green from an aqueous solution.

The batch experiments show that biosorption of dye onto Raphanus Raphanistrum is highly dependent on contact time, solution pH, initial dye concentration, biosorbent dosage, average biosorbent size and temperature.

The percentage biosorption increases with contact time, pH and biosorbent dosage and decreases with initial dye concentration, average biosorbent size and temperature.

The equilibrium data were analyzed by nonlinear Freundlich isotherm model with an increase in R2 values.

The biosorption kinetic uptake for Malachite green by Raphanus Raphanistrum at various initial dye concentrations was analyzed by nonlinear curve fitting analysis method to fit the pseudosecond order equations.

The negative values of enthalpy indicate that the process is exothermic, and the negative values of entropy indicates that the process is random and the negative value of Gibbs free energy represents that the process is feasible and spontaneous in nature.

The experimental optimum values of biosorption are t=90 min, pH=6, Ci=25mg/L and a T= 303 K having a percentage biosorption efficiency of 97.76%.

RSM results showed the optimum process parameter values as t=88.8 min, pH=5.9, Ci= 23.7 mg/L and T=301.2 K having a percentage biosorption efficiency of 98%.
REFERENCES

Indu Shekhar Thakur. "Environmental Biotechnology" I.K International Private Limited, New Delhi. 1 (2006) 422

Mathur. N.,P.Bhatnagar and P. Bakre, Assessing Mutagenicity of Textile dyes from Pali (Rajasthan) using Ames Bioassay. Applied ecology
and Environmental Research, 4(1) (2005) I2I.

Equilibrium and thermodynamic studies for dye removal using biosorption, Narayan Saibaba KV & P.King, International Journal of Research in Engineering & technology. International journal of research in engineering & technology 1 (3) (2014) 24.

Culp SJ, Beland FA. Malachite green: a toxicological review. Int J Toxicol15(1996) 238.

Hussein MMA, Wada S, Hatai K, Yamamoto A. Antimicotic activity of eugenol
against some water molts. Symp Dis Asian Aqua Philippines 1999.

Srivastava S, Sinha R, Roy D. Toxicological effects of malachite green. Aquatic Toxicology 66(2004) 329.

Removal of dyes from aqueous solution by adsorption, biosorption and phytoremediation techniques by G I.A.W. Tan, A.L. Ahmad, B.H. Hameed. 4(3) (2006) 216.

T. Santhi, S. Manonman, T.Smita and K. Mahalakshmi, Rasayan Journal of Chem., 2 (4) (2009)813.

R. Sanghi, B. Bhattacharya, Review on decolorisation of aqueous dye solutions by low cost adsorbents, Color. Technol. 118 (5) (2002)269.

A. Mittal, L. Krishnan, V.K. Gupta, Removal and recovery of malachite green from wastewater using an agricultural waste material, deoiled soya, Sep. Purif. Technol. 43 (2) (2005) 133

E. Malkoc, Y. Nuhoglu, Potential of tea factory waste for chromium
(VI) removal from aqueous solutions: thermodynamic and kinetic studies, Sep. Purif. Technol. 54(3)(2007)29.