Full Factorial Experimental Design Analysis of Reactive Dye Removal by Carbon Adsorption

The objective of this study was to investigate the removal of Remazol Yellow dye from aqueous solutions by adsorption on activated carbon prepared by chemical activation of sunflower seed cake. It was found that the carbon content of biomass increases up to 65.12% after activation and carbonization processes. The maximum percentage dye removal was obtained as 82.12% with 0.4 g/50mL adsorbent dosage at 313 K. The Langmuir model showed the best fit with equilibrium isotherm data. The interactions were evaluated with respect to both pseudo-first-order and pseudo-second-order reaction kinetics. The adsorption process was found to follow the pseudo-second-order model. To optimize the operating conditions, the effects of pH, adsorbent dosage, and initial dye concentration were investigated by full factorial experimental design method; adsorbent dosage was found as the most significant factor with P = 0.02 lower than 95% confidence level. The obtained results are very promising since (i) the utilization of sunflower seed cake activated carbon (SSCAC) played a critical role in the adsorption of dye; (ii) sunflower seed cake was an intriguing, low-cost, and easily available material. It can be an alternative adsorbent precursor for more expensive adsorbents used for Remazol Yellow (RY) removal.


Introduction
Rapidly changing technologies, industrial products, and applications are causing worldwide waste problem and contaminating the environment.If waste products are improperly managed, public health and the environment could be threatened [1,2].Among all industrial sectors, textile industries are rated as high polluters, taking into consideration the volume of discharge and effluent composition [1,3].Globally, thousands of the dye stuffs are being synthesized daily and also being released in the environment in the form of effluents during synthesis and dyeing processes [4].Most of the dyes decompose to give out hazardous products, such as carbon monoxide (CO), carbon dioxide (CO 2 ), nitrogen oxides, and hydrogen chlorides, and can reduce light penetration and photosynthesis [5].
Currently, synthetic dyes are largely used in many fields, for example, in various branches of the textile industry, in leather tanning industry, in paper production, in food technology, in agricultural research, in light-harvesting arrays, in photoelectrochemical cells, and in hair colorings [6][7][8][9][10][11][12][13][14][15].The wastewater from textile and dyestuff industries is characterized by high alkalinity, biological oxidation demand, chemical oxidation demand, and total dissolved solids with dye concentrations generally below 1 g/dm 3 [16,17].Reactive dyes are the most common dyes used due to their favorable characteristics of bright color, water fastness, simple application techniques, and low energy consumption.They exhibit a wide range of different chemical structures, primarily based on substituted aromatic and heterocyclic groups [18,19].They are not easily biodegradable; thus, the color may remain in the effluent even after extensive treatment [18,20].
Dye wastewater is usually treated by physical or chemical treatment processes.These include flocculation combined with flotation, electroflocculation, membrane filtration, electrokinetic coagulation, electrochemical destruction, ionexchange, irradiation, adsorption, precipitation, ozonation, and katox treatment methods involving the use of activated carbon and air mixtures.However, these technologies are generally unsuccessful in color removal, expensive, and less adaptable to a wide range of dye wastewaters [16,21].The low cost, simple design, easy handling, and sludge-free cleaning operations have established the adsorption technique as more effective and convenient in comparison to other techniques [1,22].
Adsorption now plays a key role in modern industries, mainly in the field of environmental protection engineering, with the increasing environmental consciousness of people all over the world.Adsorption processes are being employed widely for large-scale biochemical, chemical, and environmental recovery and purification applications [23].Over the last few decades, adsorption has been recognized as an influential separation process and has become an attractive option for the removal of dyes from industrial effluents [24].Adsorbents are used in various processes; activated carbons are among the most effective adsorbents because of their excellent adsorption capacity for organic targets [24].Many studies have been undertaken to investigate the use of low-cost adsorbents such as peat, bentonite, steel-plant slag, fly ash, china clay, maize cob, wood shavings, silica [16,[25][26][27], bacterial biomass and biopolymers [22,28], coir pith [20], sugar beet pulp [29], sugarcane bagasse pith [30], jute fiber [19,31], hen feathers [32], soybean [33,34], and wheat husk [35] for color removal.Unfortunately, these low-cost adsorbents have generally low adsorption capacities and require large amounts of adsorbents.Therefore, a need arises to find new, economical, easily available, and highly effective adsorbents [16].
Sunflower (Helianthus annuus L.) is one of the foremost oil seed crops cultivated for the production of cooking oil [36].The total world production of 38 million tons of sunflower seed, grown on about 25 million hectares [37], go almost exclusively to oil extraction, providing 8.2% of total world volume, estimated at around 107 million tones [38,39].In Turkey, sunflower is grown in the Thrace, the Aegean, and the Black Sea regions for oil production; thus, Turkey is one of the major sunflower producers and exporter countries.Oil extraction from the seeds is accompanied by the coproduction of lignocellulosic biomass, in the form of sunflower seed hulls, which comprise 30% of the sunflower seeds.These residues are usually disposed of by burning or by deposition in landfills, but conversion to higher-value products would be preferable.One such product could be activated carbon [36].
No information is available on Remazol Yellow removal from aqueous solution by sunflower biomass in the literature.This low-cost material may be especially suitable for various applications in developing countries and small-scale industries.Consequently, the scope of this study was not only to produce chemically activated sunflower seeds as an alternative low-cost adsorbent for removal of Remazol Yellow dye from aqueous solutions, but also to characterize it and to investigate the effects of operating parameters such as solution pH, adsorbent dosage, initial dye concentration, temperature, contact time, and electrolytes used.Equilibrium isotherm data were fitted to Langmuir and Freundlich equations and constants of isotherm equations were determined.Furthermore, pseudo-first and second-order kinetic models were also used to analyze adsorption kinetics.Finally, the most effective parameter (pH, adsorbent dosage, and initial dye concentration) in the adsorption process for dye removal was defined by using Full Factorial Design method and analysis of variance (ANOVA) statistical approach.

Preparation of Sunflower
Seed Cake Adsorbent and Aqueous Dye Solutions.Sunflower seed cake (SSC) obtained from Central Anatolia region of Turkey was air dried, crushed, and sieved to obtain mean sizes.Bulk density of chosen size was calculated about 0.47 g/cm 3 .Sunflower seed cake with reduced particle size was modified by impregnation with 50% wt NH 4 Cl solution and dried at 100 ∘ C in an oven.Chemical treatment with NH 4 Cl was applied to help break down the lignin complex in order to improve the performance of the adsorbent.Modified biomass material was then carbonized at 350 ∘ C to produce activated carbon (green carbon) at low temperature, washed with distilled water until the pH of the washing solution reached 6-7, and dried at ambient temperature.Finally, sieved biomass and dried sunflower seed cake activated carbon (SSCAC) were stored in plastic bottles for further use.
Reactive textile dye Remazol Yellow (RY) was used without purification.Dye solutions were prepared using distilled water to prevent and minimize possible interferences.Stock solution of dye was prepared by dissolving 1.0 g of RY dye in 1000 mL of distilled water.By diluting the standard solution of dye, different working concentrations (50-250 mg/L) were arranged.

Methodologies Used for Low
Temperature Activated Carbon Characterization.Different characterization techniques have been used to identify the produced activated carbon and the main constituents of the sunflower seed cake activated carbon.This will now be discussed.The ultimate analysis of SSC and SSCAC was performed using Elemental Analyzer (Leco CNH628 S628) to find carbon, hydrogen, nitrogen, and oxygen contents of materials by using helium, dry air, and oxygen gases.The complete combustion of all organic samples was carried out by operating Elemental Analyzer's furnace at 950 ∘ C. Structure and preliminary analyses were carried out to complete the proximate analysis of SSC.Functional groups of SSC and SSCAC were estimated by Fourier transform infrared (FTIR) spectroscopy (Perkin Elmer Spectrum 100).The FT-IR spectrums of SSC and SSCAC were obtained using the ATR technique (with a diamond protected Attenuated Total Reflectance crystal unit) with a resolution of 4 cm −1 after 100 scans.Surface morphologies of SSC and SSCAC were observed by employing a scanning electron microscope (Zeiss Supra VP 40) with an accelerating voltage of 5 kV.The samples were sputter-coated with platinum (Quorum Q 150 R ES DC Sputter).

Batch Adsorption Studies.
The effects of important parameters such as pH, adsorbent dosage, initial dye concentration, temperature, contact time, and electrolyte on the adsorption of Remazol Yellow were studied.Batch adsorption experiments were conducted in a set of conical flasks containing 50 mL dye solution of different operating conditions until the equilibrium was reached.The suspensions were then filtered and dye concentrations in the supernatant solutions were measured.Standard curves were developed at  max , 422 nm for Remazol Yellow dye, through the measurement of the dye solution absorbance by UV/Visible Spectrophotometer (V-530 Jasco UV/VIS).
The amount of dye adsorbed per unit mass was calculated using the following equation [18]: where   is the amount of dye adsorbed on the adsorbent (mg/g),  0 and   (mg/L) are the liquid phase concentrations of dye at initial and equilibrium, respectively. (L) is the volume of dye solution and  (g) is the amount of the green carbon used.The dye removal efficiency was defined as [18]  = (  0 −    0 ) 100. (2) 2.4.Kinetic and Isotherm Models.Interaction of adsorbate with the adsorbent materials is indicated by the adsorption isotherm.It is important to establish the most appropriate correlations for the equilibrium data of the system for optimizing the design of adsorption process to remove RY dye.The experimental data of RY dye was analyzed using linear Langmuir and Freundlich isotherm models to investigate the effect of temperature on the equilibrium capacity of SSCAC.Linear regression is generally used to find which isotherm is best fitted, and also the correlation coefficients of  2 are evaluated to compare the practicality of isotherm equation.The Langmuir isotherm assumes monolayer adsorption onto a surface containing a finite number of identical adsorption sites [40].Langmuir theory is based on two assumptions that the adsorbed layer is unimolecular and the forces of interaction between sorbed molecules are negligible [36].Adsorption is assumed to take place at specific homogeneous sites with the adsorbent.Once a dye molecule occupies a site, no further adsorption can take place at that site.The linear form of Langmuir isotherm equation is represented as [6] where   (mg/g) and   (L/mg) are Langmuir constants related to maximum adsorption capacity and rate of adsorption, respectively.Also,   (mg/g) is the amount of dye adsorbed per unit weight of adsorbent,   (mg/L) is the concentration of the dye solution at equilibrium.  and   values are obtained from the slope and intercept of the linear plot (  /  is plotted against   ).The Freundlich isotherm assumes heterogeneous surface energies on the adsorbent surface [41], in which the energy term in Langmuir equation varies as a function of the surface coverage.A linear form of the Freundlich equation is derived by [42] log where   (mg/g(L/mg) 1/ ) and  are the Freundlich adsorption constant and heterogeneity factor, respectively.  is related to the bonding energy and 1/ value is related to the adsorption intensity.The intercept   and the slope 1/ are obtained by plotting ln   versus ln   .
To illuminate the possible mechanism of RY dye adsorption process, pseudo-first-order and pseudo-second-order kinetic models were tested by using the data obtained from adsorption kinetic experiments.
The pseudo-first-order kinetic model of Lagergren is extensively studied and defined as [43] log where  1 is the rate constant of pseudo-first-order adsorption (min −1 ),  is the contact time (min), and   and   (mg/g) are the amounts of RY adsorbed at equilibrium and at time  (min), respectively [18].By using this linear equation and drawing a plot of log(  −  ) versus , the slope of (− 1 /2.303) and the intercept of (log   ) are obtained.
The pseudo-second-order model estimates the behavior over the whole range adsorption [18] and this model is also based on the adsorption capacity of the solid phase.The pseudo-second-order model has a supposition that the sorption process involves chemisorption mechanism and is defined as [38] where  2 is the rate constant of pseudo-second-order adsorption (g/mg ⋅ min),  is the contact time (min), and   and   (mg/g) are the amounts of RY adsorbed at equilibrium and at time  (min), respectively [38].If pseudo-second-order kinetic equation is suitable,   and  2 can be determined experimentally from the slope and intercept of plot /  versus  [44].The factorial design describes which factor shows more impact and influences the variation of one factor on the other factors [48].The three factors pH, adsorbent dosage, and initial concentration of dye were varied at two levels as given in Table 1 to investigate their effects on RY dye removal (%).Experiments for full factorial design were conducted in a set of conical flasks containing 50 mL dye solution of known pH, concentration, and adsorbent dose for 1 h at 293 K until the equilibrium was reached.After one hour of contact time, the suspensions were filtered and dye concentrations in the supernatant solutions were measured using a UV-vis spectrophotometer.

Ultimate and Proximate Analysis.
The results for ultimate and proximate analyses of sunflower seed cake and ultimate analysis of activated carbon are presented in Table 2.
Low ash and high carbon content showed that the sunflower seed cake was suitable for activated carbon production.After undergoing activation process, the carbon content and calorific value increased significantly whereas the oxygen content decreased in the sample, which indicates that activated carbon was more carbonaceous material than sunflower seed cake.

Fourier Transforms Infrared Spectroscopy (FT-IR).
The FT-IR spectra can provide valuable information about the chemical compositions of the materials.FT-IR analysis of sunflower seed cake and activated carbon was performed in the range of 4000-380 cm −1 in order to explore surface characteristics (see Figure 1).The FT-IR vibrational spectra of the sunflower seed cake and chemically activated low temperature green carbon are shown in Figures 1(a) and 1(b), respectively.The obtained FT-IR spectra revealed that the sunflower seed cake contained more bands than the prepared activated carbon.Note that the surface functional groups of sunflower seed cake experienced chemical changes during NH  at 1235 cm −1 was assigned to CH=CH stretching and also small bands ranging from 1100 to 1000 cm −1 were assigned to C-O stretching vibrations of lignin [49].Peaks located at bandwidths of 700-600 cm −1 could be assigned to C-H outof-plane bending in benzene derivatives.The C-H out-ofplane bending in benzene derivatives was commonly found on the surface of various activated carbons [18].

Morphological Characterization (SEM).
The physical morphologies and surface properties of the sunflower seed cake and activated carbon prepared by chemically activation and low temperature carbonization methods were examined by using scanning electron microscopy (SEM) technique.
Figure 2 shows the textural structure investigation of SSC and SSCAC particles from SEM micrograph at 500x and 5000x magnification, respectively.It can be concluded that SEM images of SSC and SSCAC showed severe differences.SSC had a rough surface, and thick wall structure caused a little porosity.Chemical activation and low temperature carbonization methods led opening of thick wall; thus, surface of SSCAC contained more pores than SSC surface and this results in higher surface area.

Equilibrium Adsorption Edges of Reactive Dye. Various factors that affect the adsorption performance of Remazol
Yellow on the activated carbon were examined at equilibrium.

Effect of Solution pH on SSCAC Adsorption.
The pH of the dye solution plays an important role in the adsorption capacity, where it affects both the degree of ionization of the dye as well as the surface properties of the adsorbent [6].
In this work, the influence of pH on the dye adsorption was  studied while the initial dye concentration, contact time, temperature, and amount of SSCAC were fixed at 100 mg/L, 60 min, 293 K, and 0.10 g/50 mL, respectively.The effect of pH on adsorption of dye was studied within pH range 2-8.The pH was adjusted by adding a small amount of 0.1 M HCl and/or 0.1 M NaOH; adjusted pH was measured using pH meter (Thermo Scientific Orion 3 Star).After 60 min equilibration at various pH values illustrated that the percentage adsorption decreased from 22.70% to 7.09% for RY removal when pH increased from 2 to 8 (see Figure 3); thus, the pH value of 2 was selected as the optimum for performing the adsorption studies.
The favorable adsorption of the dye at relatively low pHs in the range of 1.0-5.0 may be attributed to the electrostatic interactions, since at lower solution pH, the SSCAC may get positively charged, which enhances the negatively charged reactive dye anions through electrostatic forces of attraction and hence increase the adsorption capacity [38].Also, this behavior of dye is in accordance with the known fact that the adsorption decreases with increasing pH for anionic dyes (as the Remazol Yellow in the present case), while it increases with increasing pH for cationic dyes [6].

Effect of Adsorbent Dosage on SSCAC Adsorption.
Adsorbent dosage is an important parameter for the adsorption process as it determines the capacity of an adsorbent for a given initial concentration of the adsorbate.The effect of adsorbent dosage for optimum solution pH value of 2 was studied by varying adsorbent amount from 0.1 to 0.5 g/50 mL in increments of 0.05 g.Effects of SSCAC amount on both the percentage dye removal and the adsorption capacity were shown in Figure 4.It was observed that dye removal increased with increasing adsorbent dosage.The percentage dye removals of 61.0% and 69.10% were quite close to each other for 0.4 RY removal (%) q e (mg/g) q e (mg/g) Figure 4: The adsorbent dosage effect for Remazol Yellow adsorption (initial dye concentration: 100 mg/L, contact time: 60 min, temperature: 293 K, pH: 2).and 0.5 g/50 mL adsorbent dosages, respectively.When lowcost material was considered, 0.4 g SSCAC/50 mL solution was chosen as an optimum dosage.Increase in adsorption with adsorbent dosage can be attributed to increased adsorbent surface area and availability of more adsorption sites [6].However, the increase in the adsorbent dosage provided a remarkable decrease in the amount of dye uptake per gram of adsorbent (  ).This could be explained with two factors.First, the increase in adsorbent mass at fixed dye concentration and volume will lead to unsaturation of adsorption sites through the adsorption process, and, second, the reduction in adsorbent capacity is likely to be due to particle aggregation, resulting from high adsorbent mass.Such aggregation would lead to a decrease in total surface area of the adsorbent and an increase in diffusional path length [49].Hence, the remaining part of the experiments was carried out with adsorbent dosages of 0.1 g/50 mL and 0.4 g/50 mL (w/v) to compare lower and higher adsorbent dosages for RY dye.A-q e (mg/g) B-q e (mg/g) q e (mg/g)

Effect of Initial Dye Concentration on SSCAC Adsorption.
The effect of initial dye concentration was studied by varying initial dye concentration from 50 to 250 mg/L for 0.4 g/50 mL adsorbent dosage.It has been known that increasing the initial dye concentration results in an increase in the adsorption capacity because the initial dye concentration provides a driving force to overcome mass transfer resistances between the adsorbent and adsorption media [44].Initially, adsorbate molecules should encounter the boundary layer effect.Then, it should diffuse from boundary layer film onto adsorbent surface and, finally, it should diffuse into the porous structure of the adsorbent [18].Therefore, at higher initial dye concentration, the number of molecules competing for the available sites on the surface of activated carbon was high, hence, resulting in higher RY adsorption capacities.In the case of lower concentrations, the ratio of dye to the available sorption sites was low and higher adsorption yields were obtained.The adsorption percentage was found to decrease with increase dye concentration.This may be due to the saturation of surface area and active sites of adsorbent [41].When 0.1 g adsorbent dosage was used, the adsorption uptakes of RY at equilibrium increase from 7.42 to 11.29 mg/g as the initial dye concentration increases from 50 to 250 mg/L.Moreover, when 0.4 g/50 mL adsorbent dosage was used, the adsorption uptakes of RY at equilibrium increase from 4.60 to 10.91 mg/g as the initial dye concentration increases from 50 to 250 mg/L (Figure 5).A similar phenomenon was observed for the adsorption of reactive dyes from aqueous solution on activated carbons [18,41,44].The maximum percentage dye removal of RY was determined as 29.67% at an initial concentration of 50 mg/L for 0.1 g/50 mL adsorbent dosage and 73.61% at an initial concentration of 50 mg/L for 0.4 g/50 mL adsorbent dosage at 293 K.This interesting result shows that high concentration of RY dye attained maximum percentage dye removal even at low adsorbent dosage.Hence, sunflower seed cake activated carbon can be considered an efficient adsorbent for the removal of reactive dye from aqueous solution.A-303 K-removal (%) A-313 K-removal (%) A-293 K-q e B-293 K-q e B-303 K-q e B-313 K-q e A-303 K-q e A-313 K-q e

Effect of Contact Time and Temperature on SSCAC
Adsorption.The adsorption of RY on SSCAC at different temperatures was studied as a function of contact time (20-120 min) in order to determine the equilibrium time; experiments were conducted at 293, 303, and 313 K in an isothermal water bath shaker with initial dye concentration of 100 mg/L and rotation speed of 120 rpm.A rapid adsorption can be seen at the initial stage of the contact period.This is most likely due to a large number of surface sites for adsorption.However, it gradually slowed down until it reached equilibrium.After a lapse of time, the remaining surface sites are difficult to be occupied because of the repulsion between the solute molecules of the solid and bulk phases.Thus, adsorption took long time to reach equilibrium [18].At the point of dye removal reached a constant value where no more dye was removed from the solution.The amount of dye being adsorbed onto the adsorbent was in a state of dynamic equilibrium with the amount of dye desorbed from the adsorbent.The time required for this state of equilibrium was called the equilibrium time.The amount of dye adsorbed at the equilibrium time reflected the maximum dye adsorption capacity of the adsorbent under these particular conditions [42].It was indicated that the contact time is needed for RY solution with initial concentration of 100 mg/L to reach equilibrium was 100 min (see Figure 6).The maximum percentage removal of RY was found as 84.12% obtained at 313 K with 0.4 g/ 50 mL adsorbent dosage.The improvement in adsorption with temperature may be related to an increase in the number of active surface sites available for adsorption on adsorbent, in the porosity, and in the total pore volume of the adsorbent.It could also be due to the decrease in the thickness of the boundary layer surrounding the sorbent with temperature.Therefore, mass transfer resistance of adsorbate in the boundary layer decreases.In addition, this can be a result of an increase in the mobility of the dye molecule with an increase in their kinetic energy and the enhanced rate of intraparticle diffusion of sorbate with the rise of temperature [44].The results indicated that the adsorption reaction of RY adsorbed by SSCAC was an endothermic process in nature.

Effect of Electrolytes on SSCAC Adsorption.
Reactive dyes are the major cause for complaint.Exhaust reactive dyeing required high salt concentration (up to 80 g/L of Na 2 SO 4 /NaCl); salt is added to shift the equilibrium of dye from the aqueous phase to the solid (fiber) phase [50].Thus, effects of different electrolytes (NaNO 3 , NaCl, NH 4 NO 3 , and CaCl 2 ⋅ 2H 2 O) on the adsorption kinetics were investigated and shown in Figure 7.It was also convenient to study the effect of various salt concentrations (0.01, 0.05, 0.1, and 0.2 M) on the adsorption behavior.Electrolyte effect experiments were carried out with 100 mg/L RY solutions at 293 K and pH of 2. It was found that the addition of these electrolytes to dye solution increased the percentage dye removal and adsorption capacity.The color removal was also dependent on the concentration of added electrolytes.When CaCl 2 ⋅ 2H 2 O content increased from 0.01 M to 0.2 M, the maximum adsorption capacities for RY increased from 16.09 to 16.97 mg/g.Guelli Ulson de Souza et al. [51] found that the salt addition positively influences the adsorption compared to the experiments achieved without salts.It can be preferred that the salt cations neutralize the negative charge of the carbon surface enabling the adsorption of more molecules or the cations to act directly on the negative adsorbate ions [52].

Equilibrium Adsorption Isotherms.
As described previously, the RY adsorption isotherms carried out at 293 K values obtained by both models are summarized in Table 3.
According to Table 3, the highest correlation coefficients ( 2 values) were obtained from Langmuir equation; thus, the Langmuir model displayed better fit to the adsorption data than the Freundlich model.This result indicated that the surface of SSCAC for adsorption of RY was made up of homogeneous adsorption patches.

Adsorption Kinetic Models.
Reactive RY dye was adsorbed onto SSCAC with different adsorbent dosages (0.1 and 0.4 g/50 mL) for different time intervals and obtained experimental data were investigated to have a good agreement between pseudo-first-order and pseudo-second-order models.The adsorption rate constants ( 1 and  2 ),   and correlation coefficient ( 2 ) values at 293, 303, and 313 K were listed in Table 4.In accordance with Table 4, all of the correlation coefficients obtained for pseudo-second-order kinetic model were higher than pseudo-first-order kinetic model  2 values and were close to unity.In the case of pseudo-second order kinetics, the calculated  ,cal values were very close to the experimental  ,exp values, while pseudo-first-order model did not give reasonable values.That shows that RY adsorption on SSCAC did not match with the pseudo-firstorder model.This result exposed that a chemisorption mechanism most likely controlled the adsorption.Adsorption capacities were increased by raising temperature.This tendency proved that the adsorption occurred endothermically.

Full Factorial Design of Experiments. Dye removal by
an adsorbent in a batch system usually depends on several factors.These are acidity of medium (pH), adsorbent dosage, initial dye concentration, contact time, temperature, and so forth.Only one factor is varied by time and the others are fixed when any factor is optimized; subsequently, the best value obtained by this procedure is fixed and other factors will be varied by the time; thus, using the univariate procedure to optimize all variables is very time consuming.
The interactions among all the factors are neglected in univariate procedure, so the best condition could not be achieved [53].Full factorial design of experiments was performed by practicing with different levels of factors' all probable combinations [54].The results of the experimental design were analyzed using MINITAB 14 statistical software to evaluate the effects as well as the statistical parameters, the statistical plots (Pareto, normal probability of the standardized effects, main effects, and interaction plots). 2 3 full factorial design having 8 experiments for RY removal was studied, and a matrix was established according to their high and low levels, represented by +1 and −1, respectively.The coded values of variables with the responses (% removal efficiency) were illustrated in Table 5.The interactions between independent factors were determined with analysis of variance (ANOVA) and the main effects of dye adsorption were identified based on the  value with >95% of confidence level.
The following codified equation was used to explain the 2 3 factorial designs of RY removal by SSCAC: where  is the predicted response (removal efficiency %),  0 represents the global mean,   is the regression coefficient corresponding to the main factor effects and interactions,  is the initial solution pH,  is the adsorbent dose (g/L), and  is the initial dye concentration (mg/L).The main and interaction effects, coefficients of the model, standard deviation of each coefficient, regression coefficients, standard errors, and  and  values were shown in Table 6.All of the main factors (pH, adsorbent dosage, and initial dye concentration) and their interactions () were significant at a 5% of probability level ( < 0.05).When the factor effect is negative, removal efficiency decreases as the factor is changed from low to high levels (as seen from pH and initial dye concentration).In contrast, if the effects are positive, removal efficiency increases for high level of the same factor (as seen from adsorbent dosage).Furthermore, the fit models, submitted square correlation coefficient ( 2 ) of 0.9996, were in good agreement with the statistical model.Figure 8(a) shows the main effects of the three factors (, , and ) on removal efficiency (%).The effect of a factor is the change in response produced by the change in level of factor.This is frequently called a main effect as it refers to the primary factor of interest in the experiment [55].It was concluded that the larger the vertical line is (see Figure 8(a)), the larger the change in removal efficiency (%) when it is changing from level −1 to level +1.Please note that the statistical significance of a factor is directly related to the length of the vertical line [56].The  effect was insignificant, when it was compared with other effects.Thus the  effect was neglected and did not included in the model equation.Then, model equation, related to the level of parameters and removal efficiency, was obtained by substituting the regression coefficients in (7): Equation ( 8) indicated that two-variable interactions were significant.Proof of large positive () and negative ( and ) interactions was quite intense, thus, could not be ignored from the model.Although the main effects gave a clear idea, the interaction between those two parameters would favor a better statement of the process.obtained from ANOVA illustrated the possible positive and negative two-variable interactions among the variables , , and  for removal efficiency (%).It was observed that the effect of pH was more noticeable when the adsorbent dosage was high (see Figure 8(b)), but at lower adsorbent dosage, effect of pH was not so high.On the contrary, pH effect and adsorbent dosage effect were high at lower initial dye concentration.
Figure 9(a) shows the normal probability plot of the standardized effects with  = 0.05 to evaluate the significance of each factor and its interactions on removal efficiency (%).Normal probability plot could be separated into two regions, the region with percent above 50% where the factors were indicated positive coefficients ( and ) and the region with percent below where the factors were indicated negative coefficients (, , , and ).All these factors and interactions denoted as a circle were not significant and the effects shown as a square were significant.
Analysis of variance is a statistical method that partitions the total variation into its component parts each of which is associated with a different source of variation [48].The interaction effects are easily estimated and tested by using the usual ANOVA.The ANOVA results of RY were shown in Table 7.The sum of the squares used to estimate factors affect and Fisher's  ratios (defined as the ratio of mean square effect and the mean square error) and  values (defined as the level of significance leading to the rejection of the null hypothesis) were also represented.
Relative importance of the individual and interaction effects was given by the Pareto chart of the standardized effects in Figure 9(b).In order to identify whether the calculated effects were significantly different from zero, Student's  -test was performed and horizontal columns in Pareto chart showed these values for each effect.For a 95% confidence level and seven degrees of freedom  value was equal to 12.71.The minimum statistically significant effect magnitude for 95% confidence level is represented by the vertical line in the chart.Four values higher than 12.71 ( = 0.05) were located at right of the dash line and were significant.
It can be concluded that the adsorbent dosage was the strongest effect of the overall adsorption procedure of SSCAC.The 2 coefficient was the largest negative coefficient for model equation (8).The positive value of its coefficient meant that RY removal by SSCAC was preferred at adsorbent dosage of 8 g/L.An increase of adsorbent dosage caused an uptrend in the percentage removal of RY dye, since rising adsorbent dosage provided a suitable surface area for adsorption.The second important factor to the overall optimization of the batch adsorption process was pH.The increase in the pH involved a striking decrease of RY uptake by SSCAC.The third important factor was the initial dye concentration of dye.Surface area and active sites of SSCAC may be saturated at higher concentrations; thus, removal efficiency (%) decreased with raising the initial dye concentration.The interaction of two factors  was the fourth significant factor.The coefficient of this interaction had the negative value; thus, a decrease in pH of the solution with a deduction of the adsorbent dosage caused an increase in removal efficiency (%).This antagonistic effect would not be distinguished in the univariate optimization of the dye removal process.

Conclusions
Green carbon produced by chemical activation of sunflower seed cake with an activation agent 50% NH 4 Cl was capable of removing reactive Remazol Yellow dye molecules from aqueous solutions.The adsorption was favored at acidic medium with pH value 2 and adsorption efficiency (%) was also found to increase with the increase in adsorbent dosage, contact time, temperature, and addition of electrolyte.The kinetic study showed that dye-activated carbon adsorption systems followed by pseudo-second-order model with high correlation coefficients and the process was endothermic.The equilibrium data was in good agreement with the Langmuir model and dimensionless separation factors (  values) within the range of zero to one showed that the adsorption was favorable.Additionally, the influence of pH (2-8), adsorbent dosage (2-8 g/L), and initial dye concentration (100-250 mg/L) on removal efficiency (%) was designed by using 2 3 full factorial design.It was examined by using analysis of variance (ANOVA), -test, and -test.According to Pareto Chart, normal probability plot, main effects, and interaction plots in variance analysis, the most significant factors on removal efficiency (%) were found to be adsorbent dosage (), pH (), initial dye concentration (), and the interaction between pH and adsorbent dosage (), respectively.Conversely, the other interactions were not effective on RY adsorption.Due to the obtained results, green carbon which is produced from chemical activation of sunflower seed cake could be employed as an effective and low-cost adsorbent.Therefore, this successful adsorbent could be considered as an alternate to commercial activated carbons for the removal of reactive remazol yellow dye from aqueous solutions.

Figure 8 :
Figure 8: Plots of (a) main effects and (b) interaction effects for the removal efficiency (%).

Figure 9 :
Figure 9: (a) Normal probability plot of the standardized effects at  = 0.05, (b) Pareto chart of standardized effects on the removal efficiency (%) for RY.

Table 1 :
High and low levels of factors.

Table 2 :
The results for ultimate and proximate analyses.
a Estimated by difference.

Table 3 :
Langmuir and Freundlich coefficients and regression correlation coefficients for adsorption of RY dyes onto SSCAC at 293 K.

Table 5 :
Design matrix for Remazol Yellow.