High quantities of pollutants such as synthetic dyes, metal ions, phenols, and fertilizers discharge into surface and groundwater. The productive sources of these pollutants are mining, chemical manufacturing, metallurgical, tannery, battery manufacturing industries, textile, leather, paper, and plastics industries etc [1, 2]. Wastewaters containing dyes, even at very low concentration, can disturb the ecosystem both the environment and the human. Colored wastewater causes a reduction of sunlight penetration and affects the photosynthetic activity of aquatic organisms. In addition, these kinds of pollutants are toxic, carcinogenic, mutagenic, or teratogenic for various microbiological and fish species . Methylene blue (MB), as one of the commonly known dye, belongs to a thiazine cationic dye used in various applications like coloring paper cotton, wools and temporary hair colorant [4, 5]. Therefore, it is necessary to eliminate the risk of trapping the colors in wastewater and environmental hazards .
Various processes such as electrochemical techniques, membrane separation, ion exchange, adsorption, photo-catalysis, and photo-oxidation had used to eliminate dyes from the dyeing wastewater before its discharge to freshwater [7-9]. Among these processes, adsorption is the most widely used process due to high efficiency to treat dyes in more concentrated forms and low cost of the application and ease of operation [10, 11]. Various kinds of adsorbents have been used to remove MB dye from wastewater such as peanut hulls , walnut wood , Albizia lebbeck seed pods , peanut sticks wood , almond shell , corn straw , Harmal seeds residue  and papaya leaves . In order to improve the reactivity and adsorption capacity of adsorbent for the removal of dye, various modification materials and methods have been used; such as cationic surfactants , magnetic particles , metal oxides , chemical activating agents , organic reagents  and amino functionalities . Among these materials, magnetic nanoparticles have been considered to be an ideal candidate for modification of surface because of their unique large surface areas, well-defined pore sizes, high pore volume, and easy separation by applying an external magnetic field [24-26]. Therefore, studies on the synthesis of magnetic nanoparticles have increased during the past decade using a variety of techniques such as hydrothermal, co-precipitation, sol-gel, chemical vapor deposition, microwave assisted, ultrasonic irradiation and anodic arc plasma . Microwave radiations possess both electrical and magnetic properties. Microwave radiation is able to interact directly with the particles inside the material and changes electromagnetic energy into heat transfer inside the dielectric materials. Using Microwave radiation has the advantages such as lower activation temperature, immediate startup and shutdown, smaller steps, improved safety and improved efficiency [28. 29]. Further, Microwave radiation gives nanoparticles with smaller particle size, narrow size distribution, and high purity . Thus, these features make the microwave-assisted process as a superior method over the others as magnetic stirring applied. Myrtus Communis is a wild plant and belongs to the family of Myrtaceae that grows rapidly in different regions. Different parts of this plant were used in the food and medicine industries, such as leaves, white star-like delicate flowers and dark-blue berry fruit containing the seeds . Even though there are studies on the removal of the dye with magnetic nanoparticles covered activated carbon (AC), there is no any study for the removal of MB by using AC obtained from Myrtus Communis.
In order to find out the optimum removal efficiency in a batch system, it is necessary to identify effective parameters such as pH, contact time, dose, initial concentration, temperature, etc. Taguchi method, as a statistical tool and useful design of experimental approach, is a multivariate optimization method that reduces the cost of experiments and time of experimental investigation and presents more information about the whole process [32, 33]. However, there are few studies reported in the literature regarding the utilization of Taguchi statistical method in the field of dye removal from aqueous solutions and wastewaters.
In the present research, Fe3O4/AC nanocomposite was successfully synthesized by microwave assisted in situ co-precipitation technique at room temperature. The Taguchi statistical method with L27 orthogonal array was applied to find the optimum parameters for MB removal using new magnetite activated carbon adsorbent. The analysis of variance (ANOVA) and signal-to-noise (S/N) ratio techniques were used for identification of the most effective parameters on the removal efficiency of MB dye from aqueous solutions using magnetite activated carbon adsorbent. In addition, adsorption isotherm, kinetic and thermodynamic studies were performed to find the adsorption mechanism and efficiency of Fe3O4/AC nanocomposite.
MATERIALS AND METHODS
All chemicals and reagents used in this study were of analytical reagent grade. FeCl3.6H2O, FeCl2.4H2O, H2SO4, NaHCO3, NaOH, and MB dye were purchased from Merck. The Chemical structure of MB (chemical formula of C16H18N3SCl, λmax of 663 nm) was shown in Fig. 1. The stock solution (100 mg/L) of MB was prepared by dissolving 10 mg of MB dye in 100 ml deionized water and was diluted with deionized water to the desired MB concentrations.
Preparation of Myrtus Communis activated carbon
Myrtus Communis leaves (Raw M) were washed with tap water and then boiled in deionized water for 2 h to remove the water-soluble phenolic and other organic compounds. Subsequently, they were washed several times with deionized water and dried at 80 °C for 12 h. The dried samples were crushed and sieved (100–120 mesh). The activated carbon was prepared according to Demirbas. et al : 50 g Sieved sample was soaked in the concentrated H2SO4 solution and was activated in the air oven at 105 °C for 24 h. In order to remove the free acid, the mixture was washed with distilled water and then immersed in a 1% NaHCO3 solution. The carbonized material was placed in a ceramic crucible and activated in the furnace at 700 °C with a retention time of 3 h. Finally, the activated carbon was filtered and washed with deionized water and dried in the air oven at 105 °C.
Preparation of Fe3O4/AC
Fe3O4/AC nanocomposite was prepared by co-precipitation of Fe(II) and Fe(III) ions in aqueous solution by NaOH using microwave assisted.
At first, the activated carbon was added to a 100 mL solution containing FeCl3.6H2O, FeCl2.4H2O (molar ratio of 2:1), at room temperature and stirred for 10 min. Then, 50 ml of NaOH solution (0.5 mol/L) was added drop-wise to raise the suspension pH to around 9. The prepared sorbent was heated under microwave irradiation with an output power of 500 W for 20 min. The prepared nanocomposite was washed with deionized water until the pH of the solution become neutral and then dried in an oven at 80°C for 12 h.
The morphology of the Myrtus Communis leaves, activated carbon and Fe3O4/AC nanocomposite was examined using the field emission scanning electron microscope (FESEM) (Hitachi-S4160, Japan) under an acceleration voltage of 20 kV. X-ray diffraction studies of Fe3O4/AC nanocomposite sample was done by using X-ray diffractometer (APD 2000, Italy). The organic functional groups of Myrtus Communis leaves activated carbon and Fe3O4/AC nanocomposite was determined by FTIR spectroscopy (Shimadzu-8400S, Japan) in the wavelength of 4000–500 cm−1. The pH of MB dye solution was adjusted by using 0.1 mol/L H2SO4 and NaOH and measured with a pH meter (Metrohm-827, USA). The samples were dried in program controller JEIO TECH (CF-02G KOREA). Also, a microwave oven (MOLINEX, MVV 200130 2450 MHZ) was used for synthesized of the Fe3O4/AC nanocomposite. The residual MB concentration in the aqueous solution before and after the adsorption process was measured by a UV/Vis spectrophotometer (Perkin-Elmar-lambada25).
MB adsorption experiments
In this research, the batch experiments were studied by Taguchi statistical method. The percentage MB removal was estimated for different experimental conditions including pH, adsorbent dose, contact time, initial dye concentration and temperature. The amount of MB adsorbed per unit mass of Fe3O4/AC nanocomposite (adsorption capacity) was calculated using Eq. (1):
The percentage MB removal was calculated using the following equation:
in these formulas, Co (mg/L) and Ce (mg/L) are the concentration of MB dye before and after adsorption experiments, respectively. qe (mg/g) is the adsorption capacity of MB dye onto Fe3O4/AC nanocomposite, V (L) is the solution volume of MB dye, and m (g) is the amount of Fe3O4/AC nanocomposite adsorbent.
In this research, five controllable factors including pH, adsorbent dose, contact time, initial dye concentration and temperature were designed in three levels as shown in Table 1. A standard L27 (35) orthogonal array design (OAD) was used to determine the optimum experimental conditions for maximum removal of MB dye and the designed experimental runs are given in Table 2. The signal/noise (S/N) ratio is used to measure the deviation of the output characterization from the desired values [35, 36]. In this method, the quality characteristics of data are categorized into “the larger is the better”, “the normal is the best”, and “the smaller is the better” types. According to the goal of this research, “the larger is the better” was selected to define the optimum conditions. The S/N ratio is defined as :
In this formula, n is the replication number of the experiment and yi is the experimental response. After calculating values of S/N ratio for each experimental test, in order to obtain the optimum conditions of the removal process, the highest value of S/N ratio was selected as the optimum condition for removal process .
RESULTS AND DISCUSSION
Fig. 2 shows the FESEM images of Myrtus Communis leaves, activated carbon and Fe3O4/AC nanocomposite at different magnification, respectively. As seen from the image, the surface of Raw M leaf is soft, smooth and without any porosity. Activation process creates some changes on the surface of Raw M. These changes including a porous and rough surface with deep pores that is favorable for removal dye on the sorbent. According to Fig. 2(c), the Fe3O4 nanoparticles cover the surface of activated carbon with good homogeneity in shape and size. As seen the Fig. 2c, microwave irradiation provided nanoparticles with higher purity and uniformly in a shorter time. These observations are similar to the morphology of nanoparticles by Osouli-Bostanabad et al. and Hernández-Hernández et al. [27, 28].
The XRD pattern of Fe3O4/AC nanocomposite (Fig. 3) showed the diffraction peaks at 2θ= 29.4, 35.4, 43.2, 52.9, 57.2 and 62.6. The existence of these peaks in Fe3O4/AC nanocomposite confirms the incorporation and uniform distribution of Fe3O4 nanoparticles in the nanocomposite. Furthermore, the narrow sharp peaks indicate that the Fe3O4 nanoparticles are well crystallized.
The FTIR spectrums of Raw M activated carbon and Fe3O4/AC nanocomposite are presented in Fig. 4, respectively. The broad adsorption band at 3463 cm-1 is related to the bonded –OH stretching vibrations of OH and COOH functionalities . The band at 2360 cm−1 is corresponding to the carboxyl group. The band at approximately 1596 cm-1 is attributed to the C=C vibrations . The peaks at 1157-1118 cm-1 indicate the C-O-C groups . After activation process, several changes or shifts happened in the spectrum of activated carbon. The new peak at 2923 cm-1 is related to aldehydic C-H and alkane C-H stretching vibrations . As can be seen in Fig. 4(c), the band at 667 cm–1 (belong to the Fe3O4) is confirmation of the successful bonding of Fe coupling onto the Fe3O4/AC nanocomposite surface .
Taguchi results and effects of parameters
Results of experiments designed by the Taguchi method and the calculated S/N (signal to noise) ratio for each experiment are presented in Table 2. The obtained results of Table 2 show that the percentage of MB removal varied from 80.6% to 96.2%, and S/N ratios from 38.12 to 39.6. The results are given in Table 3 and Fig. 5 indicates that the best level for each controllable factor is level 3 for pH, level 3 for adsorbent dose, level 2 for contact time, level 1 for initial dye concentration and level 1 for temperature. On the other hand, the best levels for the MB removal are pH 7, Fe3O4/AC nanocomposite dose 0.03 g, contact time 30 min, initial MB concentration 25 mg/L and temperature 298 K. It can be observed from the above results that the order of importance of controllable factors for the removal of MB using Fe3O4/AC nanocomposite is pH, adsorbent dose, temperature, initial dye concentration and contact time, respectively. It shows that the pH factor is the most significant controllable factor of MB removal. Another ability of the Taguchi method is to predict the removal percent of MB dye in optimum condition. The amount predicted by the Taguchi method is 98.8% that is very close to this achieved experimentally in the same condition with the closeness of 98.0%.
Effect of pH
The pH of solution is one of the important parameters for the removal of dyes and it can effect on surface charge of adsorbent . According to the above results, the highest S/N ratio was achieved at the highest level of pH. On the other hands, the S/N values increases by increasing pH from 3 to 7 and maximum S/N occurred at pH 7. At low pH and acidic conditions of the solution, the dimethylamine groups in MB ions are protonated. Moreover, the surface of Fe3O4/AC nanocomposite is surrounded by the hydrogen ions that compete with MB ions for interaction with the active sites of the adsorbent. Therefore, there is a lower removal at lower pH values and adsorption process is weak due to the competition of H3O+ and MB cations. At higher pH and neutral or basic conditions of the solution, hydrogen bonding and electrostatic interactions between MB ions and functional groups of Fe3O4/AC nanocomposite during the adsorption process cause an efficient removal of MB ions from aqueous solution (Scheme 1) [44, 45]. There are abundant of hydroxyl and groups on Fe3O4/AC nanocomposite surface. As can be seen in Scheme 1, the adsorption process occurred through hydrogen bonding between the hydroxyl groups on the Fe3O4/AC nanocomposite and the electronegative groups on the MB dye molecule. Furthermore, the presence negatively charged functional groups on Fe3O4/AC nanocomposite (such as O- and COO-) and their interaction with the electronegative groups on the MB cause an efficient removal of MB ions.
Effect of adsorbent dose
The adsorbent dose is an important parameter for the removal of dyes due to presence vacant active sites and pores on the surface of adsorbent . Experiments were designed with 3 levels of adsorbent dose from 0.005 to 0.03 g. The highest S/N ratio was achieved at the highest level of adsorbent dose and maximum removal occurred at an adsorbent dose of 0.03 g. As can be seen from Fig. 4, the removal efficiency of MB increased by increasing adsorbent dose because of raising number vacant sites toward dye molecule at constant dye concentration. On the other hand, by increasing the amount of adsorbent in aqueous solution, the surface area increased, and consequently, the number of active sites available for adsorption process increased. The obtained result from this study is in accordance with the similar studies [5, 47, 48].
Effect of contact time
According to obtained results from Table 3, the highest S/N ratio was achieved at the second level of contact time. On the other words, the maximum percentage MB removal occurred at contact time 30 min. It shows that the variation of S/N ratios of contact time is the smallest and contact time is the least important variable influencing the dye removal. At the beginning of the adsorption process, the fast adsorption occurs on the surface of Fe3O4/AC nanocomposite. This may be due to the availability of active sites and functional groups on the surface and large surface area of adsorbent. By passing time, reducing vacant sites and creating repulsive forces between the dye molecules on the solid and bulk phase might be responsible for decreasing the MB removal [49, 50].
Effect of initial concentration
The effect of MB initial concentration in the range of 25 to 65 mg/L was investigated and the highest S/N ratio was achieved at the first level of initial concentration. On the other hands, the S/N value for MB concentration decreases from 39.04 to 38.93 by increasing initial concentration from 25 to 65 mg/L and maximum S/N occurred at MB initial concentration of 25 mg/L. According to these results, the removal efficiency of dye decreases with the increasing MB initial concentration. This phenomenon might be due to the saturation of the surface of Fe3O4/AC nanocomposite adsorbent or possible repulsive force between adsorbed layers on the surface of adsorbent and remaining bulk dye molecules .
Effect of temperature
Experiments of related to temperature were designed with three levels of temperature from 298 to 338 K and the highest S/N ratio was achieved at the lowest level of temperature. In the other hands, after the temperature of 298 K of the solution, the removal of dye decreased. The continuous decrease in the removal efficiency of MB dye may be due to the weak attraction between MB and adsorbent  and or increasing temperature of solution leads to creation the tendency of dye molecules to escape from the solid phase to bulk phase .
Analysis of variance (ANOVA)
In order to investigate significant parameters and percent contribution of each controllable factor to the response, the analysis of variance (ANOVA) was used. According to the results of Table 4, the degree of freedom (DOF) for each factor was 2 with a total DOF of 26. A p-value less than 0.05 indicate that the estimated model is significant . The obtained results of this indicator show that pH, adsorbent dose, initial dye concentration, and temperature are significant parameters and contact time is an insignificant parameter. Based on the ANOVA results, the rank order of the percentage contributions of each factor is as follows: (1) the pH of solution (66.81%), (2) adsorbent dose (25.54%), (3) temperature (4.83%), (4) initial dye concentration (1.23%) and (5) contact time (0.32%). According to Table 4, the parameter of contact time is the least important variable influencing the dye removal. Similar reports have been observed in other studies [33, 35]. The obtained results of Tables 3 and 4 show that the order of influence of controllable factors for the removal is same by comparison of the highest S/N ratio and percentage contributions of ANOVA.
In this part of the study, the adsorption kinetics of MB onto Fe3O4/AC nanocomposite were investigated by using the pseudo-first-order, the pseudo-second-order, and the intra-particle diffusion model.
The pseudo-first-order  model is based on the rate of adsorbate seizing the adsorption sites is proportional to the number of untaken adsorption sites. This model is:
In this formula, qt and qe are the amount of MB (mg/g) adsorbed at time t and equilibrium, respectively, and k1 (1/min) is the rate constant of the pseudo-first-order equation of the adsorption. The values of constants (qe and k1) are determined by plotting log (qe−qt) versus t.
The pseudo-second-order  model is based on the sorbent adsorbed the adsorbate chemically. This model is:
In this formula, qt and qe are the amount of MB (mg/g) adsorbed at time t and equilibrium, respectively, and k2 (g/mg.min) is the rate constant of the pseudo-second-order equation of the adsorption. The values of constants (qe and k2) are determined by plotting t/qt versus t.
The intra-particle diffusion  model is:
In this formula, ki is the initial rate of intra-particle diffusion (mg/g.min 1/2), and y (mg/g) is the intercept. The values of constants (ki and y) are determined by plotting qt versus t1/2.
The best model of kinetic adsorption is confirmed by correlation coefficients (R2). Moreover, the applicability of the best kinetic model to describe the adsorption process is confirmed by nonlinear Chi-square test (χ2) and average relative error (ARE) . These error functions are:
In these formulas, qe,exp is the observation value from the batch experiment in an equilibrium state, qe,cal is calculated value from the kinetic equations (equations 4-6) and N is the number of observations in the experimental kinetic.
The result of the effect of contact time on the adsorption of MB onto Fe3O4/AC nanocomposite is shown in Fig. 6a. The plots of pseudo-first-order and pseudo-second-order are shown in Fig. 6(b and c) and obtained results from these plots are listed in Table 5. The results show that the value of qe,cal (222.2 mg/g) is near to value of qe,exp (223.9 mg/g) for the pseudo-second-order model. Moreover, the value of the correlation coefficient is very high and close to unity. In addition, the χ2 and ARE values for the pseudo-second-order are 0.012 and 0.084, which are lower than those for the pseudo-first-order model. According to the above results, the pseudo-second-order model best fit the adsorption kinetic of MB dye onto the Fe3O4/AC nanocomposite and this indicates that the adsorption process was enhanced by chemisorption [5, 52].
In a solid/liquid adsorption process, solute diffusion into sorbent is done by external mass transfer, intra-particle diffusion, or both of them . Therefore, the intra-particle diffusion model was used to analyze the mechanism of MB adsorption. As shown in Fig. 6d, the plot of qt versus t1/2 does not pass through the origin and it separated into two linear regions. This kind of plot shows that the adsorption process is multiple stages. Therefore, external mass transfer and intra-particle diffusion control the adsorption process and intra-particle diffusion is not the sole rate-controlling step.
The aim of this part is to study of equilibrium relationship between adsorbate and adsorbent. In the other way, isotherm models can describe the correlation between the amount of adsorbate adsorbed on the surface of adsorbent and the amount of adsorbate remained in solution . This relationship was investigated by Langmuir, Freundlich, and Dubinin-Radushkevich isotherm models.
The Langmuir isotherm  is based on monolayer adsorption on the surface of the adsorbent. According to this model, there is no interaction between adsorbate ions and all adsorption active sites are energetic uniform. The linear form of the model is:
In this formula, qm (mg/g) is the maximum amount of MB adsorbed at equilibrium and b (L/mg) is the Langmuir constant related to adsorption energy. The values of the Langmuir constants are determined by plotting Ce/qe versus Ce.
The feasibility of the Langmuir isotherm model was investigated by the RL factor (the dimensional factor) . Based on this factor, the adsorption process is favorable since the RL values in the range of 0 < RL < 1.
The Freundlich isotherm  is based on multilayer adsorption on the surface of the adsorbent. In this model, the adsorption process occurs by using the non-uniform distribution of adsorption heat and affinities over the heterogeneous surface. The linear form of the model is:
In this formula, kf ((mg/g).(L/mg)1/n) indicates the adsorption capacity and n indicates the adsorption intensity. The values of the Freundlich constants are determined by plotting log qe versus log Ce.
The Dubinin–Radushkevich isotherm  applies to determine the adsorption type onto a heterogeneous surface. The linear form of the model is:
In this formula, Qm (mmol/g) is the theoretical adsorption capacity, kD (mol2/J2) is the constant related to the adsorption free energy, ε is the polanyi potential (ε =RT ln(1+1/Ce)), R (8.314 J/mol.K) is the gas constant and T (K) is the absolute temperature. The values of the Dubinin–Radushkevich constants are determined by plotting ln qe versus ε2.
The mean free energy of adsorption, as a suitable factor for estimating the type of adsorption process, was calculated using the following equation:
The adsorption capacity of Fe3O4/AC nanocomposite for MB dye was increased with increasing metal dye concentration (Fig. 7a). Fig. 7 (b-d) shows the plots of the Langmuir, Freundlich and Dubinin-Radushkevich isotherms and obtained results are listed in Table 6. By comparing the values of correlation coefficient R2 of three isotherms, it is seen that the Langmuir isotherm has the highest value of correlation coefficient (0.982). This means that the equilibrium adsorption data for MB removal fit well to the Langmuir isotherm. In addition, the values of RL are in the range of 0.087-0.035, which mentioned that the adsorption of MB onto Fe3O4/AC nanocomposite is favorable under the conditions selected in this research. Although the value of the correlation coefficient for the Freundlich isotherm is lower as compared to the Langmuir isotherm, but the value of n is higher than 1 (2.193) representing a favorable adsorption. The calculated value of E for the Dubinin–Radushkevich isotherm falls between 8.0 and 16.0 kJ/mol. Thus MB is chemically adsorbed by Fe3O4/AC nanocomposite . The maximum adsorption capacity (qm) of MB adsorption was found to be 384.6 mg/g according to the Langmuir isotherm. According to these results, the Langmuir isotherm best represents the experimental data. In the other words, adsorption of MB onto the surface of adsorbent is monolayer and homogeneous. The Langmuir maximum adsorption capacities of the various adsorbents along with Fe3O4/AC nanocomposite are presented in Table 7 [4, 5, 7, 13, 59-64]. The difference in maximum adsorption capacity of the various adsorbents is due to the difference in adsorbent type, activation method and an activation agent. Another factor is the difference in experimental conditions such as pH, temperature, and initial concentration, etc .
In this part of the study, the thermodynamic parameters including the change in free energy (ΔG°), enthalpy (ΔH°) and entropy (∆S°) were determined with equations as follow:
In these formulas, T (K) is the absolute temperature and R (8.314 J/mol.K) is the gas constant and K is the equilibrium adsorption constants of the isotherm fits .
The values of ΔH° and ∆S° are determined by plotting ln K against 1/T (Fig. 8). Based on the linear plot of Fig. 8 and obtained results in Table 8, it is obvious that the adsorption process onto Fe3O4/AC nanocomposite is spontaneous and thermodynamically favorable. The decrease of negative values of ΔG° with the temperature raising means that the adsorption process is more energetically favorable at lower temperatures . The negative value of ΔHo indicates that the adsorption process is exothermic and the positive value of ΔSo indicates the increased randomness at the solid–solution interface.
Adsorption-desorption experiments were done three times to check the reusability of Fe3O4/AC nanocomposite. Firstly, a solution of HCl (0.1 M) was used to desorb MB ions from Fe3O4/AC nanocomposite, and then the nanocomposite was washed with deionized water to reach neutral pH. After three cycles of adsorption-desorption, the removal efficiency of MB dye is shown in Fig. 9. As can be seen in Fig. 9, Fe3O4/AC nanocomposite has a potential to be used repeatedly at least for three adsorption-desorption cycles with an insignificant decrease in its removal efficiency.
In this study Fe3O4/AC nanocomposite, a new activated carbon modified iron oxide magnetic nanoparticle was synthesized by microwave assisted and characterized. We used the Taguchi design of experiment to find the best conditions for the removal of cationic dye from solution. Based on S/N ratios, the optimum process conditions were pH 7, Fe3O4/AC nanocomposite dose 0.03 g, contact time 30 min, initial MB concentration 25 mg/L and temperature 298 K. According to ANOVA results, the influence of the controllable factors in descending order is pH > adsorbent dose > temperature > initial dye concentration > contact time. The kinetics of MB adsorption on Fe3O4/AC nanocomposite follows the pseudo-second-order model. The equilibrium data followed the Langmuir isotherm with an adsorption maximum of 384.6 mg/g. This nanocomposite was effectively used for the removal of methylene blue (MB) dye from the aqueous solution.
The authors greatly acknowledgment the financial support from Young Researchers and Elite Club, Arak Branch, Islamic Azad University, Arak, Iran (No. 96066).
CONFLICTS OF INTEREST
The authors declare that there are no conflicts of interest regarding the publication of this manuscript.