Fabrication of a Cation-Exchange Membrane via the Blending of SPES/N-Phthaloyl Chitosan/MIL-101(Fe) Using Response Surface Methodology for Desalination

In the present work, a novel mixed matrix cation exchange membrane composed of sulfonated polyether sulfone (SPES), N-phthaloyl chitosan (NPHCs) and MIL-101(Fe) was synthesized using response surface methodology (RSM). The electrochemical and physical properties of the membrane, such as ion exchange capacity, water content, morphology, contact angle, fixed ion concentration and thermal stability were investigated. The RSM based on the Box–Behnken design (BBD) model was employed to simulate and evaluate the influence of preparation conditions on the properties of CEMs. The regression model was validated via the analysis of variance (ANOVA) which exhibited a high reliability and accuracy of the results. Moreover, the experimental data have a good fit and high reproducibility with the predicted results according to the regression analysis. The embedding of MIL-101(Fe) nanoparticles contributed to the improvement of ion selective separation by forming hydrogen bonds with the polymer network in the membrane. The optimum synthesis parameters such as degree of sulfonation (DS), the content of SPES and NPHCs and the content of MIL-101(Fe) were acquired to be 30%, 85:15 and 2%, respectively, and the corresponding desalination rate of the CEMs improved to 136% while the energy consumption reduced to 90%. These results revealed that the RSM was a promising strategy for optimizing the preparation factors of CEMs and other similar multi-response optimization studies.


Introduction
RSM integrating mathematics with statistics has been widely utilized to investigate, optimize and model the performance of complex systems, which can determine the significant factors affecting an experiment initiating from the design of experiment (DOE). It aims at reducing the number of experimental runs while maximizing output through the data generated. The data can be used to develop a regression model, which facilitates the acquirement of optimal parameters. There are three experimental design methods including mixed matrix ion exchange membranes have attracted much more attention because of their improvement of IEC, ionic conductivity and interfacial compatibility [6,18,26,27]. The current state of knowledge reveals that the combination of MIL-101(Fe) composite nanoparticles and an acidic functionalized polymer matrix can utilize each of their strengths to improve the performance abilities [28][29][30]. For constructing a MIL-101(Fe)-based composite membrane with flexible and mechanical ability, classic and easy-to-operate methods were adopted to ensure the dispersion of MIL-101(Fe) nanoparticles within the whole of the acidic functionalized polymer matrix. Because of hydrophilic-hydrophilic interactions, MIL-101(Fe) was expected to be confined within the hydrophilic pore/ion conducting ways of the membrane matrix, which is responsible for the pore modification.
Herein, a novel hybrid CEM with multi-dimensional nanocomposite structures by embedding MIL-101(Fe) into the SPES/NPHCs polymetric backbone was synthesized by response surface methodology. The interaction effect of three independent and crucial factors (DS, the ratio of SPES, NPHCs, and MIL-101(Fe) content) on the properties of CEMs were investigated by a Box-Behnken design model. The predicted results were exhibited from three-dimensional response surface plots to validate the interactive effect of synthesized parameters on the properties of the CEMs. The veracity of the regression model between the predicted results and experimental data was confirmed by ANOVA.

Synthesis of MIL-101(Fe)
The MIL-101(Fe) was synthesized according to a hydrothermal method which followed similarly the description of Balu et al. with a corresponding adjustment [31]. The terephthalic acid, FeCl 3 ·6H 2 O and acetic acid were gradually added to 90 mL of N, N-dimethylformamide and treated with an ultrasonic treatment for 30 min to obtain an evenly dispersed mixture. The mixture was transferred into a Teflon-lined stainless-steel bomb and heated at 110 • C for 12-24 h. The brown powder was obtained by centrifugation after the hydrothermal reaction process. The powder was thoroughly washed three times with DMF and methyl alcohol. Afterwards, the mixture was centrifuged and dried at 80 • C for 24 h under vacuum. In order to adjust the suitable size of the MIL-101(Fe) for the membranes and to optimize the operational parameters, the reaction parameters were optimized by BBD methodology and determined to be 1.07 g FeCl 3 ·6H 2 O, 0.49 g terephthalic acid, 3.18 mL acetic acid and 18 h.

Preparation of Sulfonated Polyethersulfone and N-Phthaloyl Chitosan
The drafting of the flexible sulfonic acid group onto the polymer was implemented according to a procedure that has been reported in previous work. In a 200 mL threenecked reaction flask with a nitrogen atmosphere, polyether sulfone (PES, 4 g) was stirred until the solute has dissolved in dichloromethane (50 g) at 25 • C. While the homogeneous solution formed and kept stirring, the chlorosulfonic acid (5-10 mL) was added into the homogeneous solution drip by drip at a constant speed from a constant pressure of a funnel. After full stirring and reaction, a mixture was formed, then the mixture was quickly immersed in ice water with agitation. The solid products were collected and washed by DI water until the pH was approximately 6-7. Then, the obtained sulfonated polyether sulfone was dried at 60 • C for 48 h in a vacuum oven, thus, the synthetic path of the sulfonation polyethersulfone is presented in Figure S1.
The synthesis of the N-phthaloyl chitosan (NPHCs) was performed according to a procedure reported previously. In a 200 mL reaction flask with a nitrogen atmosphere, the chitosan (1 g) and phthalic anhydride (4.48 g) were stirred until the solute was dissolved in the DMF (20 mL). The solution was heated to 130 • C for 6 h. Then, the viscous product was quickly immersed in ice water with agitation. The solid products were collected and purified by ethanol and ethyl ether. Finally, the product was dried at 40 • C for 24 h under vacuum. The synthetic path of the N-phthaloyl chitosan is shown in Figure S2.

Synthesis of Ion Exchange Membranes
The MIL-101(Fe)-sulphonated poly (ether sulfone) membrane was prepared using sulfonated polyether sulfone (SPES), N-phthaloyl chitosan (NPHCs) and MIL-101(Fe) [32][33][34][35]. Firstly, the SPES and NPHCs were added in the NMP with a nitrogen atmosphere. The solution was heated to 60 • C for 6 h and stirred until the solute dissolved. Then, the MIL-101(Fe) nanoparticles were added into the solution and stirring continued for 6 h. In order to ensure the nanoparticles were fully dispersed in the viscous solution, the solution was ultrasonicated for 45 min, then set for 24 h. The MIL-101(Fe)-sulphonated poly (ether sulfone) membrane was formed by casting the viscous solution on a glass plate and spreading it in a homogeneous liquid state. The product was dried at 100 • C for 24 h with air and then dried at 40 • C for 10 h under vacuum. Finally, the MIL-101(Fe)-sulphonated poly (ether sulfone) membrane was immersed in DI water to facilitate separation from the glass plant surface. The synthetic process of the SPES/NPHCs/MIL-101(Fe) CEMs has been presented in Figure 1. immersed in ice water with agitation. The solid products were collected and washed by DI water until the pH was approximately 6-7. Then, the obtained sulfonated polyether sulfone was dried at 60 °C for 48 h in a vacuum oven, thus, the synthetic path of the sulfonation polyethersulfone is presented in Figure S1. The synthesis of the N-phthaloyl chitosan (NPHCs) was performed according to a procedure reported previously. In a 200 mL reaction flask with a nitrogen atmosphere, the chitosan (1 g) and phthalic anhydride (4.48 g) were stirred until the solute was dissolved in the DMF (20 mL). The solution was heated to 130 °C for 6 h. Then, the viscous product was quickly immersed in ice water with agitation. The solid products were collected and purified by ethanol and ethyl ether. Finally, the product was dried at 40 °C for 24 h under vacuum. The synthetic path of the N-phthaloyl chitosan is shown in Figure S2.

Synthesis of Ion Exchange Membranes
The MIL-101(Fe)-sulphonated poly (ether sulfone) membrane was prepared using sulfonated polyether sulfone (SPES), N-phthaloyl chitosan (NPHCs) and MIL-101(Fe) [32][33][34][35]. Firstly, the SPES and NPHCs were added in the NMP with a nitrogen atmosphere. The solution was heated to 60 °C for 6 h and stirred until the solute dissolved. Then, the MIL-101(Fe) nanoparticles were added into the solution and stirring continued for 6 h. In order to ensure the nanoparticles were fully dispersed in the viscous solution, the solution was ultrasonicated for 45 min, then set for 24 h. The MIL-101(Fe)-sulphonated poly (ether sulfone) membrane was formed by casting the viscous solution on a glass plate and spreading it in a homogeneous liquid state. The product was dried at 100 °C for 24 h with air and then dried at 40 °C for 10 h under vacuum. Finally, the MIL-101(Fe)-sulphonated poly (ether sulfone) membrane was immersed in DI water to facilitate separation from the glass plant surface. The synthetic process of the SPES/NPHCs/MIL-101(Fe) CEMs has been presented in Figure 1.

Characterization Method
Powder X-ray diffraction (XRD) patterns of the MIL-101(Fe) and SPES/NPHCs membrane were obtained using a Rigaku Ultima IV X-ray Diffractometer (Tokyo, Japan) operated at 3 kW. The diffractograms were recorded in the scanning angle (2θ) range 5-40 • with a step size of 0.02 • and scanning velocity of 6 • /min. SEM (Hitachi S4800, Tokyo, Japan) was used to investigate the cross-sections and surface morphologies of the CEMs and MIL-101(Fe) at an accelerating voltage of 10 kV and a magnification of 1.00 k [36,37]. An ATR-FTIR Spectrometer (Nicolet 380, Madison, WI, USA) was performed to analyze the functional groups of the CEMs and other materials. The spectra were recorded from 4000 to 400 cm −1 by averaging 32 scans with a resolution of 4 cm −1 . The thermal stability of the membranes was measured by a differential scanning calorimeter (Mettler Toledo DSC 1, Greifensee, Switzerland) from 30 • C to 300 • C under a nitrogen atmosphere at a heating rate of 8 • C·min −1 [38][39][40]. The mechanical strength of the membrane (25 cm 2 rectangular pieces) was measured by using a controlled tensile testing machine model AI-7000-LA supported by Geotech Testing Machines Co., Ltd. (Taichung, China).

Water Content and Hydrophilicity
The membranes were immersed in DI water for 24 h. Then, the wet membranes were quickly wiped by filter paper and weighted (m 1 , g). Subsequently, the wet membranes were dried overnight and re-weighted (m 2 , g) as the dry weight of the membrane. To reduce the experimental error, each sample was measured three times. Water content was obtained from Equation (1): The hydrophilicity of the CEMs was evaluated by a contact angle using the fixed droplet method and a contact angle measuring instrument (Kruss DSA30, Hamburg, Germany). The droplet (DI water, 3 µL) dripping on the membrane surface in the dosing model and the contact angle was captured after 5 s. The measurement was repeated at least 5 times for each sample and then the average value was used. All experiments were performed at room temperature.

Ion Exchange Capacity and Fixed Ion Concentration
Ion exchange capacity (IEC) for the membranes were determined by the back-titration method. The patch of dried membrane was immersed in a HCl solution (1 mol/L) for 48 h and then removed. The remaining HCl solution was titrated by the calibrated NaOH solution until the acid and alkali neutralized. The IEC was determined by the following formula: where C HCL , V HCl , C NaOH and V NaOH are the concentrations and volumes of the HCl and NaOH solution, and m 2 is the weight of the dried membrane sample. The fixed ion concentration (FIC) was calculated by the following formula:

Diffusion Coefficient and Electrochemical Properties of the Membrane
The diffusion coefficient reflects the ability of the electrolyte to pass through the membrane under the effect of concentration, which has a certain correlation to the membrane porosity. The scheme displays the setup used for the measurement of the diffusion coefficient. The ion transport from the chamber II (NaCl solution) to chamber I, which was measured with a conductivity meter ((DDS-307, Shanghai INESA Scientific Instrument Co., Ltd., Shanghai, China), was connected to a computer. The membrane desalination performance was measured in two compartment cells, with an effective area separating the two with vertical membrane of 3.0 cm 2 NaCl solution. The process is shown in Figures S3 and S4. The diffusion coefficient was calculated by the following formula: where K s is the permeability coefficient of sodium chloride; D is the diffusion coefficient of sodium chloride; t m is the thickness of the membrane; v is the volume of sodium chloride solution; θ shows the linear relation between conductivity and content in chamber I; K shows the linear relation between time and content in chamber II; s is the effective area of membrane; and ∆C 0 is the concentration difference between chamber I and chamber II.

Characterization of MIL-101(Fe)
It was necessary to test the stability of the MOFs in a water, weak acid and alkaline environment in order to verify the feasibility of the synthesized MOFs used in the water treatment process. At room temperature, the MIL-101(Fe) materials were immersed in dilute hydrochloric acid (pH = 5) and dilute sodium hydroxide solution (pH = 9), respectively. In deionized water, the MIL-101(Fe) was kept immersed for two days, then rinsed with a large amount of deionized water for 15 min after the immersion. This was repeated three times, then it was dried in a vacuum oven at 40 • C for 24 h, followed by the detection of the MOFs after immersion in the different environments by XRD. As shown in Figure S5, the structure of the MOFs' skeletons exhibited good stability.
The characteristic peaks that appeared at 8.0 • , 8.5 • , 18.5 • , 18.0 • , and 21.3 • are in agreement with the literature [41]. Figure S6 displays the SEM images of the MIL-101(Fe) in optimum conditions. It could be observed that the surface morphology of the MIL-101(Fe) synthesis comprised smooth particles with spindle and irregular polyhedral shapes. Figure S7 exhibits the FTIR spectra of the MIL-101(Fe). The peak of the F-O stretching mode was observed at 533 cm −1 , implying the existence of a metal-oxo bond between the carboxylic group of terephthalic acid and Fe 3+ . The peak of the C-H bending vibrations of benzene in the organic linkers was surveyed at 745 cm −1 . The two intense peaks at 1388 cm −1 and 1535 cm −1 were attributed to the asymmetric and symmetric stretching of the carboxyl groups and illustrated the existence of the dicarboxylic acid anion cross-linking agent in the MIL-101(Fe). Moreover, the peak at 3432 cm −1 belonged to the O-H stretching vibrations of the water molecules adsorbed on the surface.

Characterization of Polymer Structure-Morphology (SPES and NPHCs)
The polymer skeleton is the key component of the membrane. In order to evaluate the polymer structure-morphology, FT-IR spectroscopy was used for the qualitative analysis. Moreover, the chemical nature of the materials, such as the chemical bonds, could be evaluated by identifying the specific absorption peaks for the groups [18,42]. Each compound was analyzed by using FT-IR to determine if the product had been synthesized successfully. Figures 2 and 3 represented the FT-IR of the different SPES/NPHCs-blend membranes. The absorption peak at 1780 cm −1 indicated the presence of NPHCs in the polymeric matrix, while absorption peak at 1026 cm −1 related to the sulfonic groups confirmed the presence of SPES in the membrane. With the increase of the degree of sulfonation and proportion, the characteristic peak shift amplitude of the structure was not changed significantly, indicating that the force between the two components was stable, which validated that the blend membranes had good compatibility.
XRD patterns were used to determine the change of the material crystalline nature. The crystalline nature plays a very important role in understanding the substance's solubility [43,44]. The X-ray diffractograms of the prepared membranes with different DS and SPES/NPHCs ratios are presented in Figures 4 and 5, respectively. The diffraction peaks were offset to different degrees, indicating the interactions between the SPES and NPHCs. With the increase of DS, the peak value decreased (at 2θ = 17 • ) due to the action of the hydrophilic group -SO 3 H. With the increase of the NPHCs proportion, the peak intensity (at 2θ = 22 • ) decreased, which indicated the increase of the amorphous properties of the blended film. When the NPHCs proportion exceeded 20%, the modified chitosan destroyed most of the rigid structure.
hydrophilic group -SO3H. With the increase of the NPHCs proportion, the peak intensity (at 2θ = 22°) decreased, which indicated the increase of the amorphous properties of the blended film. When the NPHCs proportion exceeded 20%, the modified chitosan destroyed most of the rigid structure.  blended film. When the NPHCs proportion exceeded 20%, the modified chitosan destroyed most of the rigid structure.

The Box-Behnken Surface Statistical Design on SPES/NPHCs/MIL-101(Fe) CEMs
The Box-Behnken design (BBD) matrix and the experimental results for the response for the CEMs are presented in Table 1. The matrix was conducted for 17 combinations, consisting of 12 trials and 5 center points for the selected conditions. The response variables and input variables were related by the following second-order polynomial equation:

The Box-Behnken Surface Statistical Design on SPES/NPHCs/MIL-101(Fe) CEMs
The Box-Behnken design (BBD) matrix and the experimental results for the response for the CEMs are presented in Table 1. The matrix was conducted for 17 combinations, consisting of 12 trials and 5 center points for the selected conditions. The response variables and input variables were related by the following second-order polynomial equation: where A, B and C are the values of the SPES content, DS and MIL-101(Fe) content, respectively.
The significance of the three independent variables for response is usually evaluated by an F-value and p-value. Values of "Prob > F" less than 0.0001 indicate the model terms are highly significant. Values of "Prob > F" less than 0.05 indicate the model terms are significant, and vice versa [1]. The predicted R 2 is an index that indicates how well the model predicts the responses to new observations. The ANOVA results for Y 1 of the quadratic model are presented in Table 2. The model "p-value" less than 0.0001 indicated that the matrix response was highly significant. The model "p-value" less than 0.05 verified that the matrix response was significant. The value of "Prob > F" less than 0.0001 exhibited that the model terms were highly significant. In this case, the "p-value" of A, B, A 2 was less than 0.001, demonstrating that these three factors were highly significant for Y 1 . The parameters of C, AC, BC, B 2 were significant for Y 1 . The normal probability plot of studentized residuals is shown in Figure 6a. The data points in this plot were located quite close to the straight line, supporting the significance of the model, and confirming that the assumptions of the analysis were satisfied. The relationship between the actual and predicted values is presented in Figure 6b. Furthermore, a good agreement was observed, indicating that the RSM model was suitable for the data range investigated in this study. The difference between the predicted R p 2 (99.58%) and adjusted R a 2 (99.92%) was 0.0034, demonstrating the high correlation between the observed and the predicted values. The 3D surface plots were graphical diagrams of regression equations showing two factors, while all other factors were maintained at fixed levels. Shown in Figure 7 are the response surface plots showing the influence of the DS, content of SPES and MIL-101(Fe) content for Y 1 . With the increase of MIL-101(Fe) doping amount (from 1% to 3%), the water content of the MMMs increased from 34% to 37%. The addition of the MIL-101(Fe) nanoparticles strengthened the hydrophilicity of the membrane, which was mainly due to the hydrophilic nature of the MIL-101(Fe). In addition, with the addition of hydrophilic MOFs, the hydrophilicity of the membrane was also improved compared with the membrane without MOFs. At the same time, the hydrophilicity of the membrane was also affected by functional groups. There were many hydrophilic functional groups in the membrane that would adsorb water molecules, and the water molecules could have acted as ion transport carriers affecting the separation performance of the membrane. The response surface analysis showed the water content of the membrane increased by about 20% with the increase of DS from 15% to 30%. When the membrane absorbed water and expanded in the water environment, the increase of the water content in the membrane would form ion cluster regions. This was also conducive to reducing the membrane resistance, but it was necessary to avoid an excessive swelling of the membrane due to water absorption, which would reduce the ion selectivity and mechanical strength of the membrane, thus reducing the service life of the membrane. Similarly, the increase of the SPES content meant that the NPHCs content in the membrane decreased, which was more hydrophilic. Comprehensive analysis revealed the WC of the membrane was related to the concentration of the ion-exchange functional groups in the membrane via macromolecular polymers (SPES and NPHCs), and that the impact was greater than from the doping of the MOFs. The increase of the SPES content in the composite membrane also increased the hydrophilic functional groups in the membrane, which in turn affected the separation performance of the membrane.  Table 3 shown the ANOVA results for Y2 of the quadratic model. The value of "P > F" less than 0.0001 indicated that the model terms were highly significant. In this ca  Table 3 shown the ANOVA results for Y2 of the quadratic model. The value of "Prob > F" less than 0.0001 indicated that the model terms were highly significant. In this case, the "p-value" of A and B was less than 0.001, demonstrating that these two factors were highly significant for Y2. The parameters of C, A 2 , and B 2 were significant for Y2. The normal probability plot of studentized residuals is presented in Figure 8a Table 3 shown the ANOVA results for Y 2 of the quadratic model. The value of "Prob > F" less than 0.0001 indicated that the model terms were highly significant. In this case, the "p-value" of A and B was less than 0.001, demonstrating that these two factors were highly significant for Y 2 . The parameters of C, A 2 , and B 2 were significant for Y 2 . The normal probability plot of studentized residuals is presented in Figure 8a. The data points in this plot were located quite close to the straight line, supporting the significance of the model, and confirming that the assumptions of the analysis were satisfied. The relationship between the actual and predicted values is exhibited in Figure 8b. The RSM model was suitable for evaluating this process. The difference between the predicted R p 2 (98.57%) and adjusted R a 2 (99.25%) was 0.0068, demonstrating the high correlation between the observed and the predicted values. The "Lack of Fit F-value" of 0.23 implied that the Lack of Fit was not significant relative to the pure error. Shown in Figure 9 are the response surface plots given the influence of the DS, content of SPES and MIL-101(Fe) content for Y 2 . The ion exchange capacity (IEC) of the membrane was an important characteristic parameter to quantify the concentration of active functional groups contained in the membrane. Through response surface analysis, when the SPES content and DS of the ion exchange membrane were 85% and 30%, the IEC of the membrane dropped slightly from about 1.15 to 1.1, with an increase of the doping amount of MIL-101(Fe) from 1% to 3%. The possible reason for this was that the addition of the MIL-101(Fe) affected the ion sites (sulfonic acid groups). The data indicated that the water content was directly proportional to the IEC, and both increased with the increase of DS (from 0.76 to 1.15). Therefore, when the water content of the membrane was controlled at an appropriate level, the increase of the IEC effectively reduced the resistance of the membrane. The SPES had a significant impact on the IEC. When the content of the MIL-101 (Fe) and DS were fixed, the IEC decreased from 1.37 to 0.78 with an increase of SPES and a decrease of NPHCs.  Table 4 shows the ANOVA results for Y3 of the quadratic model. The value of "Prob > F" less than 0.05 indicated that the model terms were significant. In this case, the parameters of A, B, C, and B 2 were significant for Y3. The normal probability plot of studentized residuals is shown in Figure 10a. The data points in this plot were located quite close to the straight line, supporting the significance of the model, and confirming that the assumptions of the analysis were satisfied. The relationship between the actual and predicted values is exhibited in Figure 10b. These results verified that the RSM model was a promising strategy for optimizing the preparation of CEMs. The difference between the predicted Rp 2 (80.26%) and adjusted Ra 2 (83.46%) was 0.032, which illustrated a good agreement between the experimental results and the predicted values. Shown in Figure  11 are the response surface plots showing the influence of the DS, content of SPES and MIL-101(Fe) content for Y3. With an increase of MIL-101(Fe) doping (from 1% to 3%), the contact angle reduced by 10%. The contact angle of the membrane decreased by 9% with an increase of DS (from 15% to 30%). This variation could be correlated with the water content analysis. Table 4. ANOVA results for Y3 of the RSM of the reduced quadratic model.  Table 4 shows the ANOVA results for Y 3 of the quadratic model. The value of "Prob > F" less than 0.05 indicated that the model terms were significant. In this case, the parameters of A, B, C, and B 2 were significant for Y 3 . The normal probability plot of studentized residuals is shown in Figure 10a. The data points in this plot were located quite close to the straight line, supporting the significance of the model, and confirming that the assumptions of the analysis were satisfied. The relationship between the actual and predicted values is exhibited in Figure 10b. These results verified that the RSM model was a promising strategy for optimizing the preparation of CEMs. The difference between the predicted R p 2 (80.26%) and adjusted R a 2 (83.46%) was 0.032, which illustrated a good agreement between the experimental results and the predicted values. Shown in Figure 11 are the response surface plots showing the influence of the DS, content of SPES and MIL-101(Fe) content for Y 3 . With an increase of MIL-101(Fe) doping (from 1% to 3%), the contact angle reduced by 10%. The contact angle of the membrane decreased by 9% with an increase of DS (from 15% to 30%). This variation could be correlated with the water content analysis.   Table 5 presents the ANOVA results for Y 4 of the quadratic model. The value of "Prob > F" less than 0.0001 indicated that the model terms were highly significant. In this case, the "p-value" of A and A 2 less than 0.001, demonstrated that these two factors were highly significant for Y 4 . The parameter of C was significant for Y 4 . The normal probability plot of studentized residuals is shown in Figure 12a. The data points in this plot were located quite close to the straight line, supporting the significance of the model, and confirming that the assumptions of the analysis were satisfied. The relationship between the actual and predicted values is shown in Figure 12b. It could be observed that the RSM model was a suitable methodology for modeling the synthesizing parameters. The difference between the predicted R p 2 (89.74%) and adjusted R a 2 (95.48%) was 0.0574, demonstrating the high correlation between the observed and the predicted values. The "Lack of Fit F-value" of 0.27 implies the Lack of Fit was not significant relative to the pure error.  Table 5 presents the ANOVA results for Y4 of the quadratic model. The value of "Prob > F" less than 0.0001 indicated that the model terms were highly significant. In this case, the "p-value" of A and A 2 less than 0.001, demonstrated that these two factors were highly significant for Y4. The parameter of C was significant for Y4. The normal probability plot of studentized residuals is shown in Figure 12a-1. The data points in this plot were located quite close to the straight line, supporting the significance of the model, and confirming that the assumptions of the analysis were satisfied. The relationship between the actual and predicted values is shown in Figure 12a-2. It could be observed that the RSM model was a suitable methodology for modeling the synthesizing parameters. The difference between the predicted Rp 2 (89.74%) and adjusted Ra 2 (95.48%) was 0.0574, demonstrating the high correlation between the observed and the predicted values. The "Lack of Fit F-value" of 0.27 implies the Lack of Fit was not significant relative to the pure error.   Figure 13 presents the response surface plots showing the influence of the DS, content of SPES and MIL-101(Fe) content for Y 4 . The fixed ion concentration (FIC) of the membrane was an important indicator of the synergistic effect between the membrane water content and the IEC. When the SPES content and DS of the ion exchange membrane were 85% and 30%, respectively, the FIC decreased from 5 to 4.5. This variation could be correlated with the IEC analysis. From the 3D response surface curve of Y 4 , the FIC decreased as the content of the NPHCs in the MMMs. When the NPHCs content was 15%, the FIC was 4.55, which was caused by the content of the NPHCs.

Effect of MIL-101(Fe) Content on MIL-101(Fe) Hybrid Membrane
With the result in the response surface plot for Y1, Y2, Y3 and Y4, the effects of MIL-101(Fe) concentrations of 0 wt%, 1 wt%, 2 wt% and 3 wt% at casting solution temperatures of 60 °C were studied. Table 6 displays the results of the performance test for these membranes. As mentioned above, the organic linkers provided a good affinity with the MOFs and organic polymers, encouraging the formation of non-covalent bonds between them. Noncovalent interactions enhanced the interface compatibility, promoting the improvement on the properties of the membrane's active layer without any adverse effect on its selectivity. Moreover, the casting solution became relatively viscous with a MIL-101(Fe) concentration from 0.0 wt% to 3.0 wt%, the water content increased from 21.9% to 29.21%, and the contact angle decreased from 67.96° to 60.45° [45][46][47][48]; however, the excessively viscous casting solution had a detrimental influence on the formation of thin membranes, and it also promoted the formation of membranes with uneven thicknesses, causing maldistribution of the membrane pores. As the MIL-101(Fe) concentration increased to 3.0 wt%, part of the MIL-101(Fe) was not well dispersed in the casting solution and led to a blockage of the membrane pore. Therefore, 2.0 wt% of MIL-101(Fe) was determined to be the optimum concentration. Table 6. The physical properties of membranes with different ratios of (SPES+NPHCs)/MOFs (at

Effect of MIL-101(Fe) Content on MIL-101(Fe) Hybrid Membrane
With the result in the response surface plot for Y 1 , Y 2 , Y 3 and Y 4 , the effects of MIL-101(Fe) concentrations of 0 wt%, 1 wt%, 2 wt% and 3 wt% at casting solution temperatures of 60 • C were studied. Table 6 displays the results of the performance test for these membranes. As mentioned above, the organic linkers provided a good affinity with the MOFs and organic polymers, encouraging the formation of non-covalent bonds between them. Noncovalent interactions enhanced the interface compatibility, promoting the improvement on the properties of the membrane's active layer without any adverse effect on its selectivity. Moreover, the casting solution became relatively viscous with a MIL-101(Fe) concentration from 0.0 wt% to 3.0 wt%, the water content increased from 21.9% to 29.21%, and the contact angle decreased from 67.96 • to 60.45 • [45][46][47][48]; however, the excessively viscous casting solution had a detrimental influence on the formation of thin membranes, and it also promoted the formation of membranes with uneven thicknesses, causing maldistribution of the membrane pores. As the MIL-101(Fe) concentration increased to 3.0 wt%, part of the MIL-101(Fe) was not well dispersed in the casting solution and led to a blockage of the membrane pore. Therefore, 2.0 wt% of MIL-101(Fe) was determined to be the optimum concentration.

Membrane Structure-Thermal STABILITY
DSC is a commonly used thermal analysis tool that helps to find the glass transition temperature (Tg) and analyzes the thermal stability of the material. Tg of the polymer is an important criterion for the compatibility of polymer components. The completely miscible polymer blend had a single Tg, while the immiscible polymer blend had a plurality of Tg. As shown in Figure 14, each DSC trace showed a single Tg, indicating good miscibility between the SPES, NPHCs and MIL-101(Fe) due to the formation of van der Waals interactions. In Figure 14, with the increase of MIL-101(Fe) content, the Tg of the blend membrane was decreased, which meant that the thermal stability of the membrane decreased with the increase of MIL-101(Fe) content. The difference in Tg values was observed due to the polymer domain interactions produced by the different forms of various MIL-101(Fe) contents and the differences in the mechanical properties of the blend membranes.

Membrane Structure-Morphology
The microscopic morphology of the membranes helps to determine the importance of membrane in the mechanisms of permeability and selectivity. Thus, the morphology of the prepared membranes was evaluated by scanning electron microscopy. Figure 15 presents SEM images of the blend membranes with surface, section, and partial section, respectively. The surface morphology of the membrane showed it was uniform and smooth without visible flaws, which indicated that the blend membranes had good compatibility. The distribution of tunnels was still relatively uniform. Under the premise of ensuring the mechanical strength and thermal stability of the membrane, the porous ion exchange membrane had a significant effect on the diffusion behavior and exhibited excellent performance in the electrodialysis process. Furthermore, the distribution of MIL-101(Fe) was relatively uniform. Pore size of the membrane is an important parameter for desalination process. Therefore, the pore size distribution of the membrane was measured by Nano Measurer for size distribution analysis to evaluate the permeability and selectivity of the

Membrane Structure-Morphology
The microscopic morphology of the membranes helps to determine the importance of membrane in the mechanisms of permeability and selectivity. Thus, the morphology of the prepared membranes was evaluated by scanning electron microscopy. Figure 15 presents SEM images of the blend membranes with surface, section, and partial section, respectively. The surface morphology of the membrane showed it was uniform and smooth without visible flaws, which indicated that the blend membranes had good compatibility. The distribution of tunnels was still relatively uniform. Under the premise of ensuring the mechanical strength and thermal stability of the membrane, the porous ion exchange membrane had a significant effect on the diffusion behavior and exhibited excellent performance in the electrodialysis process. Furthermore, the distribution of MIL-101(Fe) was relatively uniform. Pore size of the membrane is an important parameter for desalination process. Therefore, the pore size distribution of the membrane was measured by Nano Measurer for size distribution analysis to evaluate the permeability and selectivity of the membrane. Figure S8 presents the average pore size distribution for membranes with different contents (SPES/NPHCs/MIL-101(Fe) at 30% DS. With the increase of NPHCs content, the pore size had a significant increase. When the large amount of NPHCs was added (at 20%), the average pore size distribution was two to three times larger than the 15 wt% NPHCs. This result was attributed to the NPHCs tending to form clusters. It was noteworthy that the pore size did not have a significant increase due to the organic linkers present in the MIL-101(Fe) structures which provided a better affinity with the MIL-101(Fe) and organic polymers. The data of burst strength test was presented in Figure S9.

Water Content, Ion-Exchange Capacity and Hydrophilicity
A large number of hydrophilic functional groups of the ion exchange membr sorb water molecules and water molecules act as ion transport carriers, which affect the membrane separation performance. Hence, it was important to selec branes with the appropriate water content. The water content increased with sulf degree (Table 6) due to a hydrophilic sulfonic acid group. As a chitosan derivat NPHCs contained polar functional groups like hydroxyl, and ether groups that co prove the water content. With the introduction of MIL-101(Fe), the water attracting ity of the membrane increased. These results indicated that the water content was tional to the IEC. Likewise, the NPHCs, sulfonation degree and MIL-101(Fe) all

Water Content, Ion-Exchange Capacity and Hydrophilicity
A large number of hydrophilic functional groups of the ion exchange membrane adsorb water molecules and water molecules act as ion transport carriers, which directly affect the membrane separation performance. Hence, it was important to select membranes with the appropriate water content. The water content increased with sulfonation degree (Table 6) due to a hydrophilic sulfonic acid group. As a chitosan derivative, the NPHCs contained polar functional groups like hydroxyl, and ether groups that could improve the water content. With the introduction of MIL-101(Fe), the water attracting capacity of the membrane increased. These results indicated that the water content was proportional to the IEC. Likewise, the NPHCs, sulfonation degree and MIL-101(Fe) all affected the membrane ion exchange performance. The sulfonic acid groups and NPHCs provided active sites for the proper interaction between ions and the membrane surface, thereby enhancing the feasibility of the ion exchange. The contact angle was used to characterize the outermost changes in membrane modification to evaluate the hydrophilicity of the membrane. The contact angle data of different membranes are presented in Table 6. The contact angle decreased with an increase of MOF content, indicating an increase in hydrophilicity of the blend membrane. This could be ascribed to the agglomeration of the MOF. The excessive MOF content in the membrane (>3 wt%) caused an uneven dispersion of the powder. Accordingly, a MOF content of 2 wt% in the membrane did not affect the hydrophilicity of the membranes.

Membrane Porosity and Desalination
The porous ion exchange membranes had a good property. The diffusion effect depended on the presence of active ion exchange groups and pores throughout the membrane. Therefore, to investigate the change in porosity of the prepared membranes under different reaction conditions, the diffusion coefficient was measured.
As seen in Figure 16, the NaCl diffusion behavior increased with the increase of MOF content, which meant that the porosity of the composite membranes was improved by the effect of the MOF. In Figure 17, the desalination rate of the CEMs increased with the MOF content increasing. The increase in MOF content caused the membrane porosity to accumulate more electrolyte ions, thus, the diffusion of Na + was strengthened, which was consistent with the above-mentioned IEC experiments; however, the mechanical properties of the membrane reduced as the MOF content increased, which resulted in the diffusion coefficient of the membrane being 4 wt% which could not be measured. The membrane M 2 was facilitated, and when comparing it with other membranes, this study gained more satisfactory results. With the effect of the MOF additive on the NaCl diffusion behavior of the prepared CEMs, the MOF could, therefore, improve the desalination performance of the membrane compared with the membrane without a MOF. The energy consumption is presented in Figure 18. It can be observed that with a 2% MIL-101(Fe) loading into the membrane, the desalination rate of the membranes improved to 136% and the energy consumption reduced to 90%.
Membranes 2022, 11, x FOR PEER REVIEW 19 of 23 mulate more electrolyte ions, thus, the diffusion of Na + was strengthened, which was consistent with the above-mentioned IEC experiments; however, the mechanical properties of the membrane reduced as the MOF content increased, which resulted in the diffusion coefficient of the membrane being 4 wt% which could not be measured. The membrane M2 was facilitated, and when comparing it with other membranes, this study gained more satisfactory results. With the effect of the MOF additive on the NaCl diffusion behavior of the prepared CEMs, the MOF could, therefore, improve the desalination performance of the membrane compared with the membrane without a MOF. The energy consumption is presented in Figure 18. It can be observed that with a 2% MIL-101(Fe) loading into the membrane, the desalination rate of the membranes improved to 136% and the energy consumption reduced to 90%.

Conclusions
In this study, we reported the synthesis of hybrid cation exchange membranes (CEMs) by embedding MIL-101(Fe) into a SPES/NPHCs polymetric backbone by the RSM model. The properties of the CEMs observed by SEM, FTIR, XRD and DSC verified the successful reaction and good compatibility between SPES, NPHCs and MIL-101(Fe). As a result, the MIL-101(Fe) was a good pore-filler for the SPES matrix due to a good interaction with the polymer matrix and improvement in the hydrophilic nature of the membrane, which enhanced the membrane's conductivity and stability. According to a regression model and ANOVA, the optimal synthetic factors, such as DS, the ratio of SPES and NPHCs and the content of MIL-101(Fe), were determined to be 30%, 85:15 and 2%, respectively, and the corresponding desalination rate of the CEMs improved to 136%, while the energy consumption reduced to 90%. Therefore, based on the above data, the construction of SPES/NPHCs/MIL101-(Fe) hybrid membranes using a RSM statistical model is a promising approach in the design of high-performance cation exchange membrane for desalination applications.
Supplementary Materials: The following supporting information can be downloaded at: www.mdpi.com/xxx/s1, Figure S1: Schematic representation of sulfonation polyethersulfone; Figure S2: Schematic representation of the prepared N-phthloyl chitosan; Figure S3: The setup used for the measurement of diffusion coefficient; Figure S4: The setup used for the measurement of desali-

Conclusions
In this study, we reported the synthesis of hybrid cation exchange membranes (CEMs) by embedding MIL-101(Fe) into a SPES/NPHCs polymetric backbone by the RSM model. The properties of the CEMs observed by SEM, FTIR, XRD and DSC verified the successful reaction and good compatibility between SPES, NPHCs and MIL-101(Fe). As a result, the MIL-101(Fe) was a good pore-filler for the SPES matrix due to a good interaction with the polymer matrix and improvement in the hydrophilic nature of the membrane, which enhanced the membrane's conductivity and stability. According to a regression model and ANOVA, the optimal synthetic factors, such as DS, the ratio of SPES and NPHCs and the content of MIL-101(Fe), were determined to be 30%, 85:15 and 2%, respectively, and the corresponding desalination rate of the CEMs improved to 136%, while the energy consumption reduced to 90%. Therefore, based on the above data, the construction of SPES/NPHCs/MIL101-(Fe) hybrid membranes using a RSM statistical model is a promising approach in the design of high-performance cation exchange membrane for desalination applications.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/membranes12020144/s1, Figure S1: Schematic representation of sulfonation polyethersulfone; Figure S2: Schematic representation of the prepared N-phthloyl chitosan; Figure S3: The setup used for the measurement of diffusion coefficient; Figure S4: The setup used for the measurement of desalination; Figure