Background

Dyes are widely used in textile, leather, paper, rubber, plastics, cosmetics, pharmaceutical, food plants, printing and dyeing industry in China, which consumed large quantity of wastewater from different steps in the dyeing and finishing processes. It is reported that the discharge volumes of these wastewater reaches 14.13 × 108 tons per year [1]. These dyes from such industrial effluents are usually released into the environment which can be carcinogenic, mutagenic, allergenic, toxic to the aquatic life, the food web and can be resistant to natural biological degradation due to presence of these heat, light, chemical and other exposures stable [2], thus upsetting aquatic life [3]. Hence, the removal of such colored agents from aqueous effluents becomes significantly environmental, technical, and commercial importance [4].

MG is a cationic dye, an N-methylated diaminotriphenylmethane dye and has been widely used for the dyeing of leather, wool and silk, distilleries, jute, paper, cotton, leather, etc. [5], which is generally used as a strong anti-fungal, anti-bacterial and anti-parasitical agent in aquaculture [6]. The powerful antimicrobial activity of MG has been attributed to inhibition of intracellular enzymes, intercalation into DNA, and/or interaction with cellular membranes. For its high toxic, high residual, teratogenic, carcinogenic, mutagenic and other side effects to fresh water fish [7, 8] and the detection of MG though foodstuff destined to human consumption, alarmed the health hazards against human being. Though the use of this dye has been banned in several countries, it is still being used in many areas due to its low cost, ready availability and efficacy [9, 10]. In 2002, MG has been prohibited to be used as raw material for all food animals by Ministry of Agriculture of China, but the results of domestic market survey showed that: malachite green was still commonly used in aquiculture and aquatic products trafficking.

There are numbers of methods, like decolourisation, coagulation method, chemical oxidation, electrolytic process to eliminate dyes from wastewater. However, studies confirmed that the products formed after degradation of MG were also not safe and had carcinogenic potential [9, 11]. Adsorption is widely studied and applied for its low cost, high efficiency and easy operation in recent years [12]. The removal of dyes from wastewater is considered to be an important application of adsorption process using suitable adsorbent [13]. Low cost materials such as Prosopis cineraria[14], Hen feathers [15], Jute stick [16], Tea leaves [17], Tricholoma lobayense[1], Bagasse [5] and Peanut husk [18] have been reported to be used as biosorbent for the efficient removal of dyes from aqueous solution.

In this paper, a novel biomaterial P. ostreatus was used to remove MG from aqueous solutions. The macro-fungus biomasses (dead, living or their derivatives) had been studied of their adsorptive properties. Previous studies have reported that they have the capacity to bind metals [19, 20], dyes [21, 22] or other organic pollutant [23]. The macro-fungus P. ostreatus, commonly known as a mushroom, is an edible basidomycete. No matter in temperate or tropical climates, P. ostreatus can be found easily. It is now the second most important cultivated mushrooms in the world [24], and is one of the most common edible mushrooms in China. Its biomass is cheap and easily available. Therefore, the purpose of this study was to evaluate the potential of P. ostreatus as a biosorbent to remove MG from aqueous solutions. The operational parameters such as biosorbent dose, pH and ionic strength were studied. Biosorption isotherm models and kinetics models were studied. Taguchi method was used to determine the significant factors and the optimum levels of experimental factors for MG biosorption. FTIR was used to research the characteristics of biosorbent surface.

Methods

Adsorbate

MG was used as a representative cationic dye in this study. It was of analytical grade chemicals (C.I. = 42,000, chemical formula = C23H25ClN2, FW = 364.92, λmax = 617 nm) and purchased from Kelong Chemical Reagent Factory, Chengdu, China. The stock solution (1000 mg/l) of MG was prepared by dissolving accurately weighed amount of the dye in distilled water. Experimental solutions of the desired concentrations were obtained by successive diluting the stock solution by distilled water.

Biosorbent

The macro-fungus P. ostreatus was purchased from a mushroom production site near Chengdu, Sichuan. The samples were repeatedly cleaned with deionized water to remove adhering dirt and soluble impurities, and boiled with distilled water three times, each for 3–5 minutes filter with the distilled water until the filtered water was cleared. They were dried at 50°C until keep constant weight, then crushed and sieved through 200 mesh (< 75 μm) sieve to get the smaller particles. The sieved biosorbent were kept in a plastic sealable bag and stored in a desiccator without any further chemical or physical treatment before adsorption experiments.

Characterization of biosorbent

Fourier transform infrared spectrometer (FTIR) (NEXUS-650, America) was used to analyse the functional groups on the biosorbent surface, which may be responsible for MG biosorption.

The pH at the point of zero charge of the biosorbent (pHPZC), namely the pH value required to give zero net surface charge, was determined by mass titration. Three solutions with different initial pH values (3, 6 and 11, respectively) were prepared and with the addition of 0.05 M aqueous solution of HNO3 or NaOH, while NaNO3 were used as the background electrolyte. For each initial pH, six Erlenmeyer flasks were filled with 30 mL of the solution and different amounts of biosorbent were added (0.05%, 0.1%, 0.5%, 1%, 5% and 10% by weight). Then the mixture was shaken for 24 h at a thermostatic swing shaker. The pH was measured after the biosorbent was separated. A plot of the equilibrium pH versus mass fraction yielded a curve showing a plateau and the pHPZC was identified as the point at which the change of pH is zero. The pHPZC was then taken as the average of the three asymptotic pH values [22].

Experimental methods and measurements

The concentration of MG in the solution before/after equilibrium biosorption was determined by using a double beam UV/vis spectrophotometer to find out the absorbance at λmax of 617 nm. The professional data processing and mapping software Origin 8.5 was used for data analysis and graphics production. Calibration curve for MG was done before all the experiments, by the regression equation, the concentration of MG of every process could calculate precisely. The dye concentrations were measured at time t = 0 and at equilibrium, the biosorption capacity at equilibrium, q e (mg/g), was calculated by using the following equation [25]:

q e = C 0 - C e v w
(1)
P = C 0 - C e C 0 100 %
(2)

Where C 0 is the concentrations of MG at initial time and C e is the concentrations of MG at equilibrium, respectively (mg/L). V is the volume of the solution (L) and W is the biosorbent dose (g), q e is the biosorption capacity per unit mass of biosorbent (mg/g) and P for the removal efficiency (%).

Batch biosorption studies

A series of experiments were carried out in 150 ml stoppered glass Erlenmeyers flasks with 50 ml MG solution to investigate the operational parameters such as the effect of biosorbent dose, initial MG concentration, contact time, pH and ionic strength. Erlenmeyers flasks were placed at a shaker (SUKUN) with a shaking of 180 rpm and at 25°C for 6 h. After biosorption, all the samples were centrifuged with a Centrifuge (Thermo) at 10,000 g for 5 min to separate the solid phase from the liquid phase, and then analysis of the residual concentration of MG by spectrophotometer. To investigate the effect of biosorbent dose on MG removal, different weight of fungal biosorbent (varying from 0.05 g to 0.35 g) was added in each MG solution (40 mg/L). The effects of the medium pH on biosorption capacities of the biosorbent was studied from 2 to 10 (i.e. 2, 4, 6, 8 and 10), which was adjusted by adding a few drops of diluted aqueous solutions of HCl or NaOH (0.10 M). Ionic strengths of MG solutions was studied at four different NaCl (i.e. 0.05, 0.1, 0.15 and 0.2 M) and CaCl2 (i.e. 0.05, 0.1, 0.15 and 0.2 M). All experiments were conducted in duplicates and two controls were performed. Dye control (dye solution without biosorbent) were also shaken using the same procedure. Biosorbent control (biosorbent in double-distilled water without MG) with the same procedure.

Biosorption isotherms and kinetics studies

Biosorption isotherm models are used to indicate the interactions between adsorbate and biosorbent. The biosorption isotherm models were analyzed using the Langmuir [26] and Freundlich [27] isotherm models to determine the mechanism of biosorption in this study. For the biosorption isotherms studies, biosorption experiments were carried out at different initial MG concentration (i.e. 20, 40, 60, 80 and 120 mg/L).

Three common kinetics models, pseudo-first-order [28], pseudo-second-order [16] and intraparticle diffusion [29] were used to study the mechanism of adsorption. For the kinetic studies, the experimental samples were taken at predetermined time intervals for analyzing the residual dye concentration. Table 1 represents the equation of all the above-mentioned models.

Table 1 Models equations of biosorption isotherms and kinetics studies

Taguchi method

The Taguchi method was used to study the effect of various control factors on the biosorption efficient of MG biosorption from aqueous solution and to determine the significant factors and the optimum levels of experimental factors.

The Experiment design of using the Taguchi method is as follows:

  1. 1.

    The percent removal and biosorption capacity were determined as the quality characteristics to be optimized.

  2. 2.

    Selecting the controllable factors that can be set and maintained and determining their alternative levels. In this paper, the design factors include pH, contact time, biosorbent dose, initial MG concentration and three levels of each factor are represented in Table 2.

  3. 3.

    Determining and designing the orthogonal array. L9 (34) orthogonal array experiment was selected and designed as shown in Table 2.

  4. 4.

    A series of experiments were conducted from one to another as the designed orthogonal experiment.

  5. 5.

    Statistical Package for the Social Sciences (SPSS) 18.0 statistical software was used for the range analysis and the analysis of variance (ANOVA) of the data and then to determine the most important controllable factors, which can maximize the percent removal and biosorption capacity, then to select the optimal levels for those factors.

Table 2 The results of Taguchi method experiment (L 9 (3 4 ) orthogonal array experiment design)

Results and discussion

Discussion

Characterization of biosorbent

The FTIR spectra of biosorbent that before-biosorption and after-biosorption to confirm the functional groups are responsible for the biosorption process which were shown in Figure 1.

Figure 1
figure 1

FTIR spectrum of biosorbent (before-biosorption and after-biosorption).

The biosorbent had much more numbers of biosorption peaks before MG biosorption, indicating that a sufficient number of adsorption sites were existed on the surface of biosorbent. For biological biosorbent materials, it always had broad adsorption peak from 3500 to 3200 cm-1. The intense peaks at 3500 to 3200 cm-1 and 1557 cm-1 representing the stretching vibrations of amino groups [30]. The broad peak, which was observed at the level of 3500 to 3200 cm-1 for both before-biosorption and after-biosorption, indicated the presence of –OH groups and –N-H groups. The peaks at 2924, 2869 and 2852 cm-1 represent the existence of the methyl and methylene groups, and the 2869 cm-1 peak disappeared in the FTIR spectrum of after-biosorption biosorbent, indicating that these groups were responsible for the biosorption process. The peaks at 1736 and 1653 cm-1 represent the characteristics of carbonyl group stretching from aldehydes and ketones in the FTIR spectrum of biosorbent, and the peak shifted to 1652 cm-1 after biosorption [31]. The peak at 1417 cm-1 was the characteristics of bending vibration of –OH. The peak at 1375 cm-1 was for –CH3 stretching vibrations. –C-N stretching vibration of amino acid was observed at the peaks of 1321, 1252 cm-1. The peak at around 1153 cm-1 was indicative of –P = O stretching vibrations [32]. The peaks at 1083, 1042 cm-1 were the characteristics of bending vibration of –C-OH. The peak at 898 cm-1 proved that the sulfur in the protein or amino acid also existed. The band between 617 and 536 cm-1 for the biosorbent indicated that the existence of –C-N-C scissoring, which can be only found in protein structure. The FTIR spectra of biosorbent indicated that the functional groups such as, carboxyl, hydroxyl, amino and phosphonate groups on the biosorbent surface could be the potential adsorption sites for MG biosorption.

The pHPZC of biosorbent was 6.52, when the pH was below the pHPZC, the biosorbent surface was charged positively to limit the adsorption of MG. Nevertheless, when the pH was above the pHPZC of biosorbent, the surface was charged negative ion to improve the adsorption of MG.

Effect of biosorbent dose on MG removal

Figure 2 shows that the removal percentage increased with increasing in the biosorbent dose. At equilibrium time of 6 h, the percent removal increased from 58.98% to 99.58% for increasing biosorbent dose from 0.05 to 0.3 g. Above 4 g/L of biosorbent dose, the percent removal of dye MG did not bring about a significant improvement in biosorption at biosorption equilibriums. For economic considerations, optimum biosorbent dose was found to be 4 g/L. The increase of MG removal is owing to the increase of the available biosorption surface areas and availability of more possible biosorption sites with increasing biomass. Similar observations were previously reported for removal of methylene blue from aqueous solution by spent tea leaves [17].

Figure 2
figure 2

Effect of biosorbent dose on the MG biosorption.

Effect of pH on MG removal

Dyes in textile industry are complex organic compounds. It is well known that the initial pH of the solution can affect the degree of ionization of the dyes. The surface of the adsorbents is composed of polysaccharides, proteins and lipids with different functional groups (such as amino, carboxyl, thiol and phosphate groups) could bind of the dye molecules. The initial pH affects the structural stability of MG and its color intensity [33]. In the preliminary experiment, it observed by naked eyes that the dye color reduction gradually when adjustment solution pH more than 10 and lower than 2. It may be due to the structural of MG changes being effected at higher pH. MG (pKa = 10.3), the cationic dye, gets protonated in the acidic medium and deprotonated at higher pH. Consequently, the initial pH ranges from 2 to 10 in the experiment.

Figure 3 shows that the percent removal increased gradually from 35.98% to 95.92% with increasing of the pH from 2 to 10, which an increase of nearly 60% of the percent removal. For MG, a cationic dye, existed in aqueous solution in the form of positively charged ion. In addition, the pHPZC of biosorbent was 6.52, when the pH from 2 to 6 which was below the pHPZC, the biosorbent surface was charged positively to limit the adsorption of MG. Furthermore, the biosorbent surface was negatively charged when the pH was above the pHPZC of biosorbent [34]. Nevertheless, when the pH value was greater than 6, the removal rate didn’t rise as expected. The pH dependence of dye adsorption was mainly influenced by two factors: firstly, distribution of dye in the solution phase and secondly, overall charge of the adsorbent. A similar result was reported for the biosorption of Neutral Red onto peanut husk [18].

Figure 3
figure 3

Effect of solution pH on the biosorption of MG.

Effect of ionic strength

There are various impurities, which could be acids, alkalis, salts or metal ions in the wastewaters from textile-manufacturing or dye-producing industries. The presence of cations such as Na+, K+, Cu2+, Ca2+ and Cr3+, which are the most common metal ions existing in dye wastewater that may affect the performance of the biosorption process.

Figure 4 shows the effect of monovalent Na+ and divalent Ca2+ on the removal percentage of dye. It was observed that biosorptive capacity for biosorbent decreased with increasing in the concentration of Na+ and Ca2+. As concentration of Na+ changed from 0.05 M to 0.2 M, the percent removal of MG changed from 87.22% to 49.01%, with a decreased of 38.21%. As concentration of Ca2+ changed from 0.05 M to 0.2 M, the percent removal of MG changed from 83.73% to 43.57%, with a decreased of 40.16%. From 0.05 to 0.2 M, all the concentration, divalent ions Ca2+ had a greater inhibition than monovalent ions Na+ on the percent removal, which ranged from 3.49% to 5.44%. With increasing in the concentration, the gap was bigger. It was an adverse effect of ionic strength on MG removal. It may be owed to the possibility of ion exchange mechanisms in the biosorption process. The competition exist between Na+, Ca2+ and positively charged MG molecules for the same binding sites on the biosorbent surface. Ca2+ had greater positive charge than Na+ and had a stronger competitive ability. Similar observations were previously reported in aqueous solution for removal of chromium (VI) by Dunaliella species [35] and dye removal by Poly (propylene imine) Dendrimer [36].

Figure 4
figure 4

Effect of Na + and Ca 2+ on percent removal.

Biosorption isotherms

As seen of Table 3, the regression coefficient R2 (0.996-0.998) of the linear Freundlich model were higher than the determination coefficient R2 (0.906-0.972) of Langmuir model and R2 (0.951-0.956) of Temkin model at all the temperatures, which indicated that Freundlich model was well fitted to data than the Langmuir and Temkin isotherm for this study. Moreover, the maximum MG biosorption capacity (q m ) 111.11 mg/g, 125 mg/g, 125 mg/g which got from Langmuir model at 15°C, 25°C, 35°C were found to be much higher compared with the experimental result 24.437 mg/g. These demonstrated that the process of MG biosorption was in terms of Freundlich model, not Langmuir and Temkin model. For Freundlich model, KF (l/g) and n are Freundlich constants which indicated the ease of biosorption and adsorption intensity, respectively. The value of n was bigger than 1, which indicated the biosorption of MG actively and was a heterogeneous nature biosorption and had rapid adsorption process at studied conditions [37]. KF and R2 at 25°C were higher than 15°C, 35°C. The reason could be that with the temperature increased, the molecular of MG movement became faster, more collision frequency between biosorbent and sorbate. It may also be the higher affinity of binding sites on the biosorbent surfaces.

Table 3 Isotherm constants for MG biosorption

Adsorption kinetics

As seen of Figure 5, the MG biosorption occurred in two stages. In the first stage, it was a rapid uptake within 60 min, then in the second stage, it was slow uptake from 60 to 240 min and reached equilibriumat at 240 min, in this stage, the percent removal of MG did not bring about a significant change.

Figure 5
figure 5

Biosorption at different contact time and different initial MG concentrations.

As seen of Table 4, the R2 values for pseudo-first-order model was range from 0.912 to 0.989, which was lower than the R2 values for pseudo-second-order model which ranged from 0.999 to 1.000. Moreover, the theoretical qe,cal values for all the studied MG concentration in terms of pseudo-second-order model were more close to the experimental qe,exp values compared with the pseudo-first-order model.

Table 4 Kinetic parameters for the removal of MG (temperature = 25°C, stirring rate = 180 rpm, pH = 8)

The initial concentration of MG in the solution remarkably influenced the biosorption. It was noted that initial concentration increased the biosorption of MG (Figure 5). This increase in uptake capacity of the biosorbent with the increase in initial dye concentration is due to higher availability of dye ions for the biosorption [38].

Figure 6 was the linearized plots of the pseudo-second-order kinetic model for MG biosorption. The data of adsorption kinetics were fitted well to the pseudo-second-order linear model for different MG concentration.

Figure 6
figure 6

The pseudo-second-order model for MG biosorption at different initial concentrations.

As showed in Figure 7, the biosorption process was a multi-steps process which included a fast step and a slow step. In the first step, the adsorption active sites on external surface were enough for a fast biosorption. After the first step, the biosorption would get to a steady phase which the intraparticle diffusion was diffusion-controlled.

Figure 7
figure 7

Intraparticle diffusion model for MG biosorption at different initial concentrations.

As seen in Table 5, Kid,1 was larger than Kid,2 which demonstrated the biosorption in the first step got a larger rate of transporting the adsorbed MG from exterior to interior of biosorbent particle. The Ai,2 was bigger than Ai,1 the larger intercept for the greater boundary layer effect, which demonstrated that after a fast biosorption process, most of MG was on the biosorbent, and most of the adsorption active sites on external surface were occupied and made a thick layer. Therefore, the biosorption rate was controlled by the intraparticle diffusion rate. As seen in Figure 7, there were no plots pass through the origin, thus, it demonstrated that intraparticle diffusion was not the only rate control step. Maybe there were other effect such as complexation or ion exchange mechanisms took effect in biosorption process, and the R2 values for intraparticle diffusion model was ranged from 0.8735 and 0.9957 which were lower than the R2 values for pseudo-second-order kinetic model.

Table 5 Intraparticle diffusion kinetic model (temperature = 25°C, stirring rate = 180 rpm, pH = 8)

Many factors influence the adsorption behavior such as structure of dye and adsorbent surface properties. There are many forces which control the adsorption such as van der Waals, hydrogen bonding, steric effect, ion exchange etc. It can be concluded by the results of kinetics that the biosorption process was well in terms of the pseudo-second-order model, which demonstrated that chemisorption was the rate-controlling mechanism through sharing or exchange ions between biosorbent and MG in the biosorption process [39].

Taguchi method

Range analysis (R) of the results of experiment

As seen of Table 2, the trial 5 with the combination of pH 8, contact time 6 h, biosorbent 0.25 g, and initial MG concentration 80 mg/L had the best percent removal and the trial 8 with the pH 9, contact time 6 h, biosorbent 0.15 g, and initial MG concentration 120 mg/L had the best biosorption capacity. Range analysis (R) was used to study the importance of the four factors and get the optimal combination.

The results of R are shown in Table 6. The value of R indicates which factor is more important for the percent removal. K1, K2, and K3 are the average values of the results of the three different levels of each control factor and can determine the optimal level of the specific control factor in Taguchi method experiment. The R was used to estimate the effect of the specific control factor. No matter the value of ki (K1, K2, and K3) or the value of R, all of them comply with the principle of “The larger the better”. As seen, for percent removal, AK2,BK2,CK3,DK1 were larger than others results, therefore, the combination of A2B2C3D1 which pH 8, contact time 6 h, biosorbent 0.25 g, initial MG concentration 80 mg/L could reach the best percent removal. Through the experimental test, the percent removal reached 89.58%. From Table 6, it can be seen for RC>RD>RA>RB, indicating that the biosorbent does was the most important factor on percent removal. For the reason that sufficient biosorbent may afford available sorption surface and availability of more adsorption sites. For biosorption capacity, AK3,BK2,CK1,DK3 were larger than others results, therefore, the combination of A3B2C1D3 which pH 9, contact time 6 h, biosorbent 0.15 g, initial MG concentration 120 mg/L could have the best biosorption capacity. From Table 6, it can be seen for RC>RD>RB>RA, indicating that the biosorbent does was the most important factor on biosorption capacity. For the reason that excessive MG molecule could be fully biosorbed on the surface of the biosorbent. Through the experimental test, the biosorption capacity reached 32.33 mg/g.

Table 6 Range analysis of the experiment

Analysis of variance

Analysis of Variance (ANOVA) is a mathematical method which distinguishes the differences between the test results of different experimental conditions and the differences between the test results of causal factors, and the result of ANOVA could obtain which control factor had significantly effect on percent removal and the biosorption capacity [40]. In order to find which factor significantly effect on remove rate and the biosorption capacity, the conditions differences and random differences influenced by any special factor need to be compared. Square sum of deviations were divided by the Corresponding degree of freedom, the reason of the sum of squares were related to the number of observation and the number of groups, therefore, conditions differences and random differences could not compared directly.

In the ANOVA, the Fischer ratio (F test) was used to quantify the whole process. The F values of each factor were the ratio of the mean of the squared deviations to the mean of the squared error. The critical values for F test of confidence level at 0.1 and 0.05 were got in experimental design handbooks. Then the F values compared with the critical value at different confidence levels. If the F values was larger than F0.05, indicating that the experiment factor was particularly significant effect on the results and marked with “**”. If F value was larger than F0.10 and smaller than F0.05, indicating that the experiment factor was significant effect on the results and marked with “*”. At last, if the F value was smaller than F0.10, indicating that the experiment factor had no significant effect on the results.

As seen of Table 7, biosorbent does and initial MG concentration had significant influences on the percent removal. Biosorbent does and initial MG concentration had particularly significant influences on biosorption capacity. The values of F were 47.716 and 27.279 for biosorbent does and initial MG concentration, respectively. From the results of ANOVA, it indicated that the biosorption process was not influenced by the time and the environment pH under the condition of this study. In practice, it will show stability in the biosorption process.

Table 7 ANOVA for percent removal and biosorption capacity in the L 9 (3 4 ) orthogonal array experiment

Conclusions

This study confirmed that Pleurotus ostreatus has the potential to be an economical biosorbent for the removal of MG, a hazardous cationic dye, from wastewater. The biosorption effect in alkaline condition was better than in acidic condition, but pH affected the structural stability of MG and its color intensity when pH was more than 10. MG was almost adsorbed by biosorbent when the initial MG concentration was low in solution. Biosorbent dose, time and pH were important for biosorption in the single factor analysis. Monovalent Na+ and divalent Ca2+ were found to have an adverse effect on MG removal. Equilibrium time was nearly 240 min. The biosorption equilibrium data were followed the Freundlich isotherm model. The biosorption process was found to follow the pseudo-second-order kinetic model, indicating that chemical biosorption was the rate-controlling mechanism through sharing or exchange ions between biosorbent and MG in the biosorption process.

For the Taguchi method experiment, the optimum condition for percent removal and biosorption capacity were found out. The combination of A2B2C3D1 which pH 8, contact time 6 h, biosorbent 0.25 g, initial MG concentration 80 mg/L could reach the highest percent removal 89.58%. And the combination of A3B2C1D3 which pH 9, contact time 6 h, biosorbent 0.15 g, initial MG concentration 120 mg/L could reach the largest biosorption capacity 32.33 mg/g. The results of FTIR spectroscopy showed the existence of the functional groups such as, carboxyl, hydroxyl, amino and phosphonate groups on the biosorbent surface.