Translate this page into:
Optimization and modeling of synthesis parameters of neodymium(III) bromide by dry method using full factorial design analysis
⁎Corresponding author. berkanima@yahoo.fr (Madjid Berkani)
-
Received: ,
Accepted: ,
This article was originally published by Elsevier and was migrated to Scientific Scholar after the change of Publisher.
Peer review under responsibility of King Saud University.
Abstract
The synthesis of neodymium(III) bromide (NdBr3) by sintering brominating of neodymium oxide (Nd2O3) with ammonium bromide (NH4Br) was investigated. The influence of various synthesis parameters (temperature, contact time and stoichiometry) on the reaction yield was studied and optimized. The main interaction effects of the synthesis parameters on the reaction yield were also determined by a full 23 factorial designs with six replicates at the center point.
This study showed that the optimum conditions for the synthesis of NdBr3 are following: contact time t = 60 min, stoichiometry in moles Nd2O3:NH4Br = 1:24 and temperature T = 400 °C. The reaction yield for these parameters was equal to 97.80%. The first order model was obtained to predict the reaction yield as a function of these three parameters. It was shown that all parameters have a significant positive influence on reaction yield. In addition it was pointed out also that the interaction effects between them are significant.
Keywords
Neodymium(III) bromide
Lanthanide halides
Reaction yield
Full factorial design
1 Introduction
In the past 50 years, lanthanide halides were the subject of several studies. The majority of these studies were focused on their chemical and physical properties. The physicochemical properties of the compounds they form, for instance with the alkali halides, can be examined and correlated in a way which is unique since the lanthanides regularly vary along the series they form in the periodic table. This is due to the presence of f electrons in the lanthanide ions, which do not participate in the formation of the chemical bonding (Gschneidner et al., 2003). These compounds are used in various fields such as: reprocessing of nuclear wastes (Kobayashi et al., 1998), recycling of spent nuclear fuel (Ogawa et al., 1995; Ogawa and Igarashi, 1998), catalysis and chemical synthesis (Doty et al., 2007), medical imaging (Beckert et al., 2016), and in the lighting industry (Muklejohn et al., 1993; Markus et al., 2005). These different applications require knowledge of their physicochemical, structural and thermodynamic properties. The study of these properties requires high purity of these salts in order to have accurate data (Rycerz and Gaune-Escard, 2008; Kolodziej et al., 2009; Chojnacka et al., 2014). These anhydrous salts that are commercially available still contain small amounts of water and cannot be used for such a study. Therefore, good control of the different stages of synthesis and optimization of parameters that can affect the synthesis reaction is the preliminary step to cross (Gaune-Escard et al., 1994; Peringer et al., 2008; Mendil et al., 2013). The determination of these optimum conditions requires a large number of experiments and takes a lot of time. These drawbacks can be avoided by applying the experimental design method, which allows studying the effects of several factors on one or more responses, and to find a mathematical model that connects the response to the factors (Kincl et al., 2005). This method has several advantages such as: providing maximum information with minimum experiments, the possibility of increasing the number of factors studied or their levels, taking into account possible interactions between factors, simple modeling results and a good search of the optimal response (Lazic, 2006; Tinsson, 2010). The aims of the present work were as follows: firstly, optimization of synthesis parameters (temperature, contact time and molar ratio neodymium oxide: ammonium bromide) of neodymium(III) bromide (NdBr3) by sintering brominating of neodymium oxide (Nd2O3) with ammonium bromide (NH4Br), secondly, application of a 23 full factorial design in order to determine a mathematical model connecting these three factors with reaction yield and to gain insight into how the various factors interact and influence the response.
2 Experimental
NdBr3 was synthesized from neodymium oxide Nd2O3 (Sigma-Aldrich, 99.9%) by sintering brominating with NH4Br (Sigma-Aldrich, 99%). A mixture of Nd2O3 and NH4Br in well-defined proportions was ground and homogenized in a mortar, and placed into alumina crucible. Then, the crucible was inserted into a quartz reactor and placed in an electric furnace. The furnace temperature was controlled by the programmer EUROTHERM 2404. This resulted in a process where volatile products were absorbed in a scrubber filled with sodium hydroxide. The entire experiments were carried out in an argon atmosphere. The yield of the reaction was determined by chemical analysis. Neodymium(III) bromide is soluble in water whereas unreacted neodymium oxide is water insoluble. Thus dissolution of the reaction product transfers neodymium bromide into solution, whereas neodymium oxide stays in insoluble solid residue. The resulting product of the reaction was cooled down, and then dissolved in 150 cm3 of 0.05 M HCl solution using magnetic stirrer for 10 min. In this step obtained neodymium(III) bromide was transferred to the solution. The solution of 0.05 M HCl was used to prevent hydrolysis of neodymium ions formed during the dissolution of the neodymium bromide. Buchner funnel with hard paper filter was used to separate the leach solution from the sediments. Filtrates (leach solutions) were collected in a 200 cm3 volumetric flask, and filled by the solution of 0.05 M HCl to the mark. Then a 20 cm3 of the resulting solution was diluted 10 times in a 200 cm3 volumetric flask. Next, 10 cm3 of the diluted solution was transferred to a conical flask and 10 cm3 of buffer solution (CH3COOH/CH3COONa), pH = 5.6 was added. The prepared solution was heated almost to boil (80 °C), and then four drops of xylenol orange solution were added. The chemical analysis of the content of Nd was performed by a complexometric titration with a standard EDTA solution until the color changed from purple–red to light yellow. The titration was repeated 3 times.
3 Full factorial designs
A full factorial design is an experimental design with a number of levels limited to two (low and high) for each factor; it takes into account all the combinations between the levels of factors during the experiment. The full factorial design allows modeling of the first degree with or without interactions and it quickly finds the application limits when the number of factors increases. This design also allows us to determine the effects of the factors on the response and the effects of interactions between different factors (Montgomery, 2012).
4 Results and discussion
4.1 Thermodynamic aspect of the reaction
Before synthesis of NdBr3 by sintering brominating with NH4Br, we performed a thermodynamic study of the reaction to determine the standard thermodynamic parameters such as enthalpy, entropy and the standard Gibbs free energy of this reaction. These thermodynamic parameters provided information about the exothermic or endothermic synthesis process, the degree of disorder, and clarified the possibility of spontaneous synthesis process.
The synthesis of NdBr3, from Nd2O3 and NH4Br by dry route, can be described by the following reaction:

The results of Fig. 1a show that the reaction is endothermic between room temperature and T = 606 K (sublimation temperature of NH4Br). The sudden increase of enthalpy observed at 373 K is associated with the passage of water from liquid to a gaseous state.
Fig. 1b shows that the entropic contribution is positive over the entire temperature range, this is due to the increase in disorder, which can be explained by the increase in the number of moles of gas.
The free enthalpy of reaction as a function of the temperature was calculated according to the equation:
4.2 Thermal analysis of reagents by TG–DTG
The two reagents (Nd2O3 and NH4Br) used for the synthesis of NdBr3, have been analyzed with TG–DTG in order to verify their stability in function of the temperature. The TG–DTG curves obtained with heating rate 10 °C/min under a nitrogen atmosphere are shown in Figs. 2 and 3.

Fig. 2 shows that NH4Br begins to lose the weight at 179 °C and greater weight loss, corresponding to decomposition temperature according to the reaction (3), was observed at 333 °C.
The curves TG–DTG of Nd2O3 (Fig. 3) shows that the oxide does not lose any weight and that during this process the oxide remained stable. This observation confirms the literature information concerning the lanthanide oxides (Sarbak et al., 2005).
4.3 Optimization of the synthesis parameters
In this work, we studied the influence of three main synthesis parameters (contact time, stoichiometry and temperature) on the reaction yield.
4.3.1 Influence of the contact time
The study of contact time influence on the reaction yield was realized by varying it from 10 to 120 min while maintaining the two other parameters constant (the molar ratio Nd2O3: NH4Br = 1:6 and experimental temperature T = 250 °C).
The reaction yield, noted R, was calculated with the following equation:

This figure shows that the yield of the reaction R (%) increases with the increase in reaction time up to 60 min. Beyond this value, the yield remains almost constant and stabilizes at an average of 30%. Therefore, the optimum contact time of this reaction is 60 min, and for this reason, we set the contact time to 60 min in order to study the effect of the stoichiometry and temperature.
4.3.2 Influence of the stoichiometry
The study of this parameter was carried out in order to see the effect of excess NH4Br on the reaction yield. It is based on the variation of the stoichiometry in moles of Nd2O3:NH4Br from 1:6 to 1:30. The two other parameters were kept constant (t = 60 min and T = 250 °C). The results obtained are represented in Fig. 5.
Increase in the reaction yield R (%) was observed with excess of NH4Br up to the molar ratio Nd2O3:NH4Br = 1:24. Beyond this value, the yield was stabilized at a level of about 46%. The optimal molar ratio Nd2O3:NH4Br was found to be 1:24.
4.3.3 Influence of the temperature
The thermodynamic study of the synthesis reaction of NdBr3 showed that the synthesis is promoted by an increase in temperature. According to this result, we have chosen to vary the temperature from 250 to 450 °C, and kept the two other parameters at their optimum level (t = 60 min and molar ratio Nd2O3:NH4Br = 1:24). Fig. 6 represents the experimental results of this study.
It is evident that an increase in temperature up to 400 °C has a positive effect on the reaction yield. Beyond this value, the yield decreases slightly. This is probably due to the loss of NH4Br related to its decomposition at high temperatures (Mendil et al., 2013). As a result, the optimum temperature for synthesis of NdBr3 by sintering brominating with NH4Br is 400 °C for which yield reaches 97.80%.
4.4 Factorial design analysis
After optimization of the three factors influencing the reaction yield of the synthesis of NdBr3, we have applied the experiment design method to model the reaction yield depending on the three factors. This method allows us to find a mathematical model that links the reaction yield to the three factors, and allows us to examine the main and interaction effects of factors on the reaction yield.
In this study, a full factorial design 23 with six experiments at the central point and therefore a total of 14 experiments were used. The replication of experiments in central point is in order to evaluate the experimental error (Mendonça et al., 2011). The high, center and low levels defined for the 23 factorial designs were listed in Table 1. The low and high levels of the factors were selected according to the optimization results.
| Factors | Low level (−1) | Center level (0) | High level (+1) |
|---|---|---|---|
| ( ) Contact time (min) | 10 | 35 | 60 |
| ( ) Stoichiometry in moles (Nd2O3:NH4Br) | 1:6 | 1:15 | 1:24 |
| ( ) Temperature (°C) | 250 | 325 | 400 |
The mathematical model in the coded form for the full factorial design 23 is given as follows:
The results obtained from the 14 experiments are shown in Table 2, with the low (−1), center (0), and high (+1) levels as specified in Table 1.
| Experiments | R (%) | |||
|---|---|---|---|---|
| 1 | −1 | −1 | −1 | 7.40 |
| 2 | 1 | −1 | −1 | 30.84 |
| 3 | −1 | 1 | −1 | 14.99 |
| 4 | 1 | 1 | −1 | 46.04 |
| 5 | −1 | −1 | 1 | 21.84 |
| 6 | 1 | −1 | 1 | 79.76 |
| 7 | −1 | 1 | 1 | 30.07 |
| 8 | 1 | 1 | 1 | 97.80 |
| 9 | 0 | 0 | 0 | 47.40 |
| 10 | 0 | 0 | 0 | 48.55 |
| 11 | 0 | 0 | 0 | 48.17 |
| 12 | 0 | 0 | 0 | 49.32 |
| 13 | 0 | 0 | 0 | 47.01 |
| 14 | 0 | 0 | 0 | 50.86 |
4.4.1 Statistical analysis
After realization of the 14 experiments, we have calculated the various coefficients of the model (
), and then we have determined the main effects and interactions between different factors. This effect which represents the difference between the average value of the response (R) at the high level (+) and the average value of the response at low level (−) was determined using the formula (6) (Cavalitto and Mignone, 2007).
The significance of the different calculated coefficients has been checked by the Student's t-test, which allows us to determine whether among the coefficients of the model, there are no significant coefficients, which will be eliminated from the regression equation because their influence on the reaction yield is negligible (Erto et al., 2010).
The ti value which represents the ratio of the coefficient to the estimated parameter and standard deviation residual are calculated using the following equation:
The regression coefficients, the effects and t values are shown in Table 3.
| Term | Effects | Coefficients | Standard error | ti | Significance |
|---|---|---|---|---|---|
| b0 | 41.09 | 41.09 | 0.9506 | 83.15 | S |
| b1 | 45.03 | 22.51 | 0.9506 | 45.55 | S |
| b2 | 12.27 | 06.13 | 0.9506 | 12.41 | S |
| b3 | 32.55 | 16.28 | 0.9506 | 32.95 | S |
| b12 | 04.36 | 02.18 | 0.9506 | 4.41 | S |
| b13 | 17.79 | 08.90 | 0.9506 | 18.01 | S |
| b23 | 00.87 | 00.44 | 0.9506 | 0.89 | NS |
| b123 | 00.55 | 00.28 | 0.9506 | 0.57 | NS |
S: Significant, NS: Not Significant.
For a 95% confidence level and 5 degrees of freedom (m-1), the tabulated value of STUDENT (tcrit value) is 2.57.
The results are illustrated by means of Pareto charts (Fig. 7). The vertical line in Pareto charts indicates minimum statistically significant effect magnitude for a 95% confidence level. The effects which have ti less than 2.57 are statistically insignificant (Montgomery, 2012; Carmona et al., 2005).
According to the Pareto chart (Fig. 7), the interaction, stoichiometry - temperature (coded ) and time-stoichiometry-temperature (coded ) are insignificant for a 95% confidence level.
Thus, the reduced model, represented by (Eq. (5)) in terms of coded parameters after excluding the insignificant terms for the predicted reaction yield
(%) is as follows:
4.4.2 Variance analysis (ANOVA)
The significance and adequacy of the regression model after excluding insignificant coefficients was evaluated by the analysis of variance (ANOVA) (Sen and Swaminathan, 2004).
The ANOVA is based on the comparison of the variation associated with the model that represents the sum of the quadratic variation of the experimental yield (R) and the predicted yield ( ) to the experimental average yield ( ) and variation associated with experimental error which represents the sum of the squared residual between the experimental and predicted yields (Kasiri and Khataee, 2011).
The sum of squared total (SST), sum of squares for regression (SSR) and sum of squared errors (SSE) were computed from the reaction yield as follows (Tinsson, 2010):
| Experiments | |||||
|---|---|---|---|---|---|
| 1 | 7.40 | 7.25 | 0.022 | 1145.32 | 1135.19 |
| 2 | 30.84 | 30.11 | 0.533 | 120.62 | 105.12 |
| 3 | 14.99 | 15.15 | 0.026 | 673.02 | 681.35 |
| 4 | 46.04 | 46.73 | 0.476 | 31.78 | 24.48 |
| 5 | 21.84 | 22.01 | 0.029 | 364.15 | 370.66 |
| 6 | 79.76 | 80.47 | 0.504 | 1550.58 | 1495.17 |
| 7 | 30.07 | 29.91 | 0.026 | 125.05 | 121.50 |
| 8 | 97.80 | 97.09 | 0.506 | 3135.71 | 3215.84 |
| SUM | – | – | SSE = 2.121 | SSR = 7146.222 | SST = 7149.300 |
The validity of the regression equation is based on the statistical of Fisher, which consists in comparing of F-value defined as the ratio between the mean-square-effect and the mean-square-error, and the tabulated value of the F-distribution at a level of significance and a certain number degrees of freedom.
The value of F is calculated by the following equation:
The model is a good predictor of the experimental results if the F-value is greater than the tabulated value of the F-distribution (Yetilmezsoy et al., 2009; Liu et al., 2004).
The calculated F-value for reaction yield “F = 1347.57” is much greater than the tabulated ( at the 95% significance), the regression equation is thus very significant on the response, and the adjustment of the model is adequate.
4.4.3 Correlation coefficients
The correlation coefficient or determination coefficient (R2) is the ratio between the variation due to the regression and the total variation, this coefficient enable us to check the correlation between the experimental and the predicted responses (Liu and Chiou, 2005; Santos and Boaventura, 2008; Deming and Morgan, 1993), it is calculated as follows:

The goodness of fit of the model was also checked by adjusted coefficient
(Eq. (14)), which corrects the correlation coefficient value R2 for the sample size and the number of terms in the model.
In our study, the value of (0.993) was found very close to the corresponding value R2.
4.4.4 Residual analysis
The residuals analysis, which represents the difference between the experimental and predicted responses, allows us to represent the unexplained part of the model and check whether there is a failure in the latter. This test consists of analyzing the plot of observed residues depending on the predicted responses (Fig. 9) (Kasiri and Khataee, 2011; Sado and Sado, 2000).
The look of the graph shows a random distribution of points, confirming the nonexistence of a relationship between predicted responses and residue. So the model of first level that we have established well explains the experimental results (Kasiri and Khataee, 2011; Goupy and Creighton, 2006).
4.5 Main and interaction effects
A main effect can be defined as the effect of one factor on the response, ignoring the effects of all other factors, it represents the change in the response when the factor changes from low to high levels. The effect of a factor is positive if the response R (%) increases when the factor changes from low to high levels. Otherwise (the response R (%) decreases when the factor changes from low to high levels) the effect is negative (Ponnusami et al., 2007).
The main effects of the experimental results for the respective low and high levels of contact time, molar ratio and temperature were illustrated in Fig. 10.
From the graph (Fig. 10), we can see that an increase in contact time (a), stoichiometry (b), temperature (c), from low to high levels resulted in an increase in the reaction yield R (%). This indicates that all factors have a positive effect on the response (Palanikumar et al., 2009; Lorenzen and Anderson, 1993).
Contrast to the main effect, the interaction effect represents the effect of a factor on the response as a function of other factors. An interaction between two factors is significant when the effect of a factor on the response depends on the level of the other factor, i.e. whether the lines of the two factors are not parallel (Lazić, 2004; Mathialagan and Viraraghavan, 2005).
The interactions between factors were illustrated in Fig. 11.
Fig. 11 shows the significant interactions between contact time and stoichiometry (Fig. 11a), as well as between the contact time and the temperature (Fig. 11b). However, the interaction between stoichiometry and temperature (Fig. 11c) is neglected (the lines are almost parallel). These plots also indicate that the interaction between the contact time and the temperature ( ) is stronger than between the contact time and stoichiometry ( ). These plots also show that the effect of the contact time is larger when the temperature and the stoichiometry are at high (+1) levels.
5 Conclusion
In these study, the synthesis parameters (contact time, stoichiometry and temperature) of neodymium(III) bromide by dry method are optimized. According to the results, the optimum conditions for the synthesis of NdBr3 are as follows: contact time (t = 60 min), stoichiometry in moles (Nd2O3:NH4Br = 1:24) and temperature (T = 400 °C), for which the reaction yield reached 97.80%.
The full factorial design was used to allow the evaluation of the most important factors for synthesis of neodymium(III) bromide. This modeling leads to a first degree mathematical model, which relates the yield to the three operating parameters. It has also shown that an increase from lower to higher levels in factors has a significantly positive influence on reaction yield R (%), with presents of significant interaction between contact time (t) and stoichiometry, and between contact time (t) and temperature. The validity of this study was limited to, contact time between 10 and 60 min, stoichiometry between 1:6 and 1:24 and temperatures between 250 and 400 °C.
Acknowledgments
The work was finances by a statutory activity subsidy from the Polish Ministry of Science and Higher Education for the Faculty of Chemistry of Wroclaw University of Technology.
Financial support by the Thematic Agency for Research in Science and Technology (ATRST) from the Algerian Ministry of Higher Education and Scientific Research is gratefully acknowledged.
M. B and Y. B wish to thank the Department of Chemistry of Wroclaw University of Technology for hospitality and support during this work.
References
- Thermochemical Data of Pure Substances (third ed.). Germany: VHC Verlagsgesellschaft mbH; 1995.
- Medical imaging scintillators from glass-ceramics using mixed rare-earth halides. Opt. Mater.. 2016;60:513-520.
- [Google Scholar]
- Biosorption of chromium using factorial experimental design. Process Biochem.. 2005;40(2):779-788.
- [Google Scholar]
- Application of factorial and Doehlert designs for optimization of protopectinase production by a Geotrichum klebahnii strain. Process Biochem.. 2007;42(2):175-179.
- [Google Scholar]
- Direct preparation of anhydrous lanthanide chlorides from lanthanide oxides chlorinated by NH4Cl. Chin. Rare Earths. 2008;2:015.
- [Google Scholar]
- Phase diagram and specific conductivity of the DyBr 3–CsBr binary system. J. Alloy Compd.. 2014;582:505-510.
- [Google Scholar]
- Experimental Design: a Chemometric Approach. Elsevier; 1993.
- Structure and properties of lanthanide halides. Int. Soc. Opt. Photon. 2007 670705-670705-11
- [Google Scholar]
- Experimental and statistical analysis of trichloroethylene adsorption onto activated carbon. Chem. Eng. J.. 2010;156(2):353-359.
- [Google Scholar]
- Calorimetric investigation of NdCl3- MCl liquid mixtures (where M is Na, K, Rb, Cs) Thermochim. Acta. 1994;236:67-80.
- [Google Scholar]
- Introduction aux plans d’expériences (3ème édition). Paris: Dunod; 2006.
- Gschneidner. K.A., Bünzli. J.C., Pecharsky V.K., 2003. Handbook on the Physics and Chemistry of Rare Earths, Vol. 33, first ed., Elsevier, North Holland.
- Photooxidative decolorization of two organic dyes with different chemical structures by UV/H2O2 process: experimental design. Desalination. 2011;270:151-159.
- [Google Scholar]
- Application of experimental design methodology in development and optimization of drug release method. Int. J. Pharma.. 2005;291:39-49.
- [Google Scholar]
- Phase diagram and electrical conductivity of the AgBr–NdBr 3 binary system. J. Alloy. Compd.. 2009;480(2):834-838.
- [Google Scholar]
- Design of Experiments in Chemical Engineering: a Practical Guide. New York: John Wiley & Sons; 2006.
- Design and Analysis of Experiments: Section 2.3. Design of Experiments in Chemical Engineering: a Practical Guide. 2004. p. :262-367.
- Handbook of Chemistry and Physics (87th ed.). Boca Raton, USA: CRC Press (Taylor and Francis Group); 2007.
- Optimal decolorization efficiency of Reactive Red 239 by UV/TiO 2 photocatalytic process coupled with response surface methodology. Chem. Eng. J.. 2005;112:173-179.
- [Google Scholar]
- Optimal production of sulphuric acid by Thiobacillus thiooxidans using response surface methodology. Process Biochem.. 2004;39(12):1953-1961.
- [Google Scholar]
- Design of Experiments: a No-name Approach. New York: CRC Press; 1993.
- High temperature gas phase chemistry for the development of advanced ceramic discharge lamps. J. Phys. Chem. Solids. 2005;66(2):372-375.
- [Google Scholar]
- Biosorption of pentachlorophenol by fungal biomass from aqueous solutions: a factorial design analysis. Environ. Technol.. 2005;26(5):571-580.
- [Google Scholar]
- Determination of the optimum conditions for the synthesis of praseodymium (III) chloride. C. R. Chim.. 2013;16(9):795-798.
- [Google Scholar]
- Application of full factorial design and Doehlert matrix for the optimisation of beef tallow methanolysis via homogeneous catalysis. Fuel Process. Technol.. 2011;92(3):342-348.
- [Google Scholar]
- Design and Analysis of Experiments (eighth ed.). New York: John Wiley & Sons; 2012.
- Muklejohn, S.A., Devonshire, R., Trindel, D.L., 1993. High Temperature Lamp Chemistry III, New Orleans.
- Pyrochemical processes in advanced nuclear programs-with emphasis on management of long-lived radionuclides. High Temp. Mater. Processes.. 1998;2:587-596.
- [Google Scholar]
- Ogawa, T., Yamagushi, S.Y., Itoh, A., Mukayama, T., Hanada, M., Haire, R.G., Proc. Int; 1995. Conference on Evaluation of Emerging Nuclear Fuel Dense Systemes, France.
- Assessment of some factors influencing tool wear on the machining of glass fibre-reinforced plastics by coated cemented carbide tools. J. Mater. Process. Technol.. 2009;209:511-519.
- [Google Scholar]
- On the synthesis of LaCl3 catalysts for oxidative chlorination of methane. J. Appl. Catal. A: Gen.. 2008;350:178-185.
- [Google Scholar]
- Biosorption of reactive dye using acid-treated rice husk: factorial design analysis. J. Hazard. Mater.. 2007;142:397-403.
- [Google Scholar]
- Lanthanide (III) halides: thermodynamic properties and their correlation with crystal structure. J. Alloy. Compd.. 2008;450:167-174.
- [Google Scholar]
- Plans d’expérience, de l’expérience à l’assurance de qualité. France: Afnor; 2000.
- Adsorption modelling of textile dyes by sepiolite. Appl. Clay Sci.. 2008;42:137-145.
- [Google Scholar]
- Binary Rare Earth Oxides. Dordrecht, The Netherlands: Kluwer Academic Publishers; 2005.
- Response surface modeling and optimization to elucidate and analyze the effects of inoculum age and size on surfactin production. Biochem. Eng. J.. 2004;21(2):141-148.
- [Google Scholar]
- Plans d'expérience: constructions et analyses statistiques. Berlin: Springer Science & Business Media; 2010.
- Response surface modeling of Pb (II) removal from aqueous solution by Pistacia vera L.: Box-Behnken experimental design. J. Hazard. Mater.. 2009;171:551-562.
- [Google Scholar]
