RSM Modeling and Optimization of CO2 Separation from High CO2 Feed Concentration over Functionalized Membrane

The challenges in developing high CO2 gas fields are governed by several factors such as reservoir condition, feed gas composition, operational pressure and temperature, and selection of appropriate technologies for bulk CO2 separation. Thus, in this work, we report an optimization study on the separation of CO2 from CH4 at high CO2 feed concentration over a functionalized mixed matrix membrane using a statistical tool, response surface methodology (RSM) statistical coupled with central composite design (CCD). The functionalized mixed matrix membrane containing NH2-MIL-125 (Ti) and 6FDA-durene, fabricated in our previous study, was used to perform the separation performance under three operational parameters, namely, feed pressure, temperature, and CO2 feed concentration, ranging from 3.5–12.5 bar, 30.0–50.0 °C and 15–70 mol%, respectively. The CO2 permeability and CO2/CH4 separation factor obtained from the experimental work were varied from 293.2–794.4 Barrer and 5.3–13.0, respectively. In addition, the optimum operational parameters were found at a feed pressure of 12.5 bar, a temperature of 34.7 °C, and a CO2 feed concentration of 70 mol%, which yielded the highest CO2 permeability of 609.3 Barrer and a CO2/CH4 separation factor of 11.6. The average errors between the experimental data and data predicted by the model for CO2 permeability and CO2/CH4 separation factor were 5.1% and 3.3%, respectively, confirming the validity of the proposed model. Overall, the findings of this work provide insights into the future utilization of NH2-MIL-125 (Ti)/6FDA-based mixed matrix membranes in real natural gas purification applications.


Introduction
The East Natuna Block, located in the Greater Sarawak Basin near Natuna Island, Indonesia, is the largest natural gas field in Southeast Asia. It contains 46 trillion cubic feet (Tcf) of natural gas, and is identified as a high carbon dioxide (CO 2 ) content gas field, with CO 2 concentrations up to 70 mol%. In addition, about 13 trillion cubic feet (Tcf) of high CO 2 gas fields in Malaysia remain unexploited, with CO 2 concentrations of more than 40 mol%. According to the requirements of the gas processing plant, the CO 2 content should be reduced to 8 mol% in order to comply with the criteria of the downstream sales gas process [1]. High CO 2 contents can reduce the heating value of natural gas [2] and cause corrosion issues in pipelines in the presence of water [3]. Thus, it is crucial to investigate technologies for separating CO 2 from methane at high CO 2 feed concentrations. Among conventional technologies such as absorption [4], adsorption [5], cryogenic distillation, and membrane separation [6], membrane separation has gained the most attention due to its small footprint, simple operation and energy efficiency [7].
The evolution of membrane materials in separation technology started with polymeric, followed by inorganic, mixed matrix and metal organic frameworks (MOFs), etc. Mixed matrix membranes comprised of polymeric material and inorganic filler were developed to address the trade-off in performance experienced by polymeric membranes [8] and reproducibility concerns regarding the fabrication of inorganic membranes [9]. Lately, incorporating MOF fillers in polymeric membranes for the formation of mixed matrix membranes has been the focus of many researchers due to the fact that the high porosity of MOF fillers can improve gas separation performance [10].
On top of that, numerous efforts have been made by researchers to further enhance the separation performance of mixed matrix membranes by functionalizing metalorganic frameworks (MOFs) fillers with -NH 2 functional groups. The presence of -NH 2 in the filler induces strong affinity toward CO 2 via hydrogen bonding, which enhances the interaction between CO 2 and the filler framework [11,12], and subsequently, improves the performance of the membrane in terms of CO 2 separation. Owing to its tunable pore size and affinity to CO 2 gas [13][14][15], NH 2 -MIL-125 (Ti) has recently been studied as a filler for the formation of membranes for gas separation [3,16,17]. Previous studies reported that the separation performance of NH 2 -MIL-125 (Ti)/polysulfone mixed matrix membranes was superior to that of neat polysulfone membranes [17]. This finding was attributed to the existence of porous fillers, which created an extra pathway for the penetration of CO 2 gas and promoted the diffusion of gas molecules through the membrane. Furthermore, in another work, a NH 2 -MIL-125 (Ti)/Matrimid mixed matrix membrane was fabricated for CO 2 /CH 4 separation and a CO 2 /CH 4 gas pair selectivity of 50 was obtained when 15 wt% filler was loaded into the membrane [16]. This was mainly due to the small aperture size (6 Å) of NH 2 -MIL-125 (Ti), which allowed CO 2 molecules to interact with the -OH groups in Ti clusters, as well as interacting with -NH 2 groups on the linker. Also, our previous study discovered that the inclusion of NH 2 -MIL-125 (Ti) filler in 6FDA-durene polymer matrices significantly improved the CO 2 permeability and CO 2 /CH 4 gas pair selectivity, i.e., by 119% and 331%, respectively, compared to neat 6FDA-durene membranes [3]. The improvement in performance was mainly attributed to the high porosity of NH 2 -MIL-125 (Ti) fillers and the attraction of amine functional groups to CO 2 molecules. Although several studies have been published on the utilization of functionalized MIL-125 (Ti) as a filler for the preparation of mixed matrix membranes in CO 2 /CH 4 separation, these studies mostly focused on single gas permeation, binary gas separation (50:50) and the effect of feed pressure on the membrane performance [16,17]. Hence, further research evaluating the separation performance of the membrane at high CO 2 feed concentrations and various temperatures and CO 2 feed pressures in binary gas mixtures is still crucial.
According to the findings in the literature, operating temperature influences the CO 2 permeability of the membrane by affecting the mobility of the polymer chain and free volume [18] in the polymer, thus resulted in a change of CO 2 diffusivity [19]. In contrast, feed pressure may affect the gas solubility in the membrane, which affects the CO 2 permeability [11]. Furthermore, the CO 2 concentration in the feed can cause significant changes in CO 2 permeability through the membrane, which is mainly due to the higher sorption capability of CO 2 compared to CH 4 . Hence, study of the correlation between the operational parameters and membrane separation performance is essential. In fact, several approaches, including response surface methodology (RSM), full factorial design and Taguchi methods, are available to investigate these questions.
For an experimental assessment and the optimization of the separation parameters, including significant responses and independent variables, comprehensive analytical tools and efficient optimization tools are practical. The multivariant method enabled the identification of the interaction between the variable parameters and responses. Furthermore, for our optimization study, the response surface methodology (RSM) was identified as the optimal tool. This technique combines the mathematical and statistical analyses by reducing the number of experiments. As a result, various researchers have reported the application of the RSM method in their studies.
Jusoh et al. [20] investigated the effect of pressure, temperature and CO 2 feed concentration on CO 2 /CH 4 binary gas separation over a ZIF-8/6FDA-based mixed matrix membrane. They reported that the optimum parameter for CO 2 /CH 4 separation performance was obtained at a feed pressure of 4.76 bar, a temperature of 30 • C, and a CO 2 feed concentration of 90 mol%. The latest study reported by Afarani et al. [18] utilized a polyurethane-zeolite 3A mixed matrix membrane to assess the impact of three independent variables, namely, zeolite content (0-24 wt%), operating temperature (25-45 • C), and operating pressure (0.2-0.1 MPa), on the single gas permeation performance. The optimum gas permeation performance of polyurethane-3A zeolite membrane was identified at a zeolite loading of 18 wt%, a temperature of 30 • C, and a pressure of 0.8 MPa. It can be observed that the RSM method has great potential for analyses of the relationship between the operational parameters and responses. Thus, the process of optimization of the separation performance over NH 2 -MIL-125 (Ti)/6FDAdurene mixed matrix membranes by varying the operational parameters including feed pressure, temperature and CO 2 feed concentration using RSM method is presented in this study.
The 7.0 wt% NH 2 -MIL-125 (Ti)/6FDA-durene mixed matrix membrane reported in our previous work [3] was used in this study, since the membrane demonstrated the best performance in CO 2 and CH 4 single gas permeation testing. First, the interaction of operational parameters, including feed pressure (A), temperature (B), and CO 2 feed concentration (C), toward the CO 2 and CH 4 permeability and CO 2 /CH 4 separation factor were analyzed in detail using the RSM statistical method paired with the CCD optimization tool. Then, the significance of models was further evaluated by using analysis of variance (ANOVA). Next, 3D surface plots were used to visualize the linear, quadratic, and interaction between the operating parameters and significant responses. Subsequently, the optimum condition for the separation performance was generated and validated by experimental work.
The premixed gases used for the gas separation test were purchased from Air Products Sdn Bhd, (Kuala Lumpur, KL, Malaysia) with CO 2 compositions ranging from 15 mol% to 70 mol%.

Synthesis of 6FDA-Durene
In this study, 6FDA-durene polyimide was synthesized using a dual-step procedure reported in the literature [21]. First, equimolar amounts of durene diamine and 6FDA dianhydride monomers were dissolved in NMP for 24 h at room temperature under a nitrogen atmosphere. Next, polyimide was formed after adding propionic anhydride and triethylamine to the mixture solution. Then, the solution was precipitated in methanol and washed several times with methanol. The resulting polymer was dried under a vacuum for 24 h at 150 • C.

Synthesis of NH 2 -MIL-125(Ti)
The NH 2 -MIL-125 (Ti) filler was synthesized by following the method reported in our previous work [3]. Firstly, 17.5 g of pivalic acid was dissolved in a mixture solution containing 125 mL of acetonitrile and 5 mL of tetraisopropyl orthotitanate. Next, the mixture solution was placed in a Teflon-lined autoclave reactor and heated at 100 • C for 84 h to obtain white crystals. The white crystals were filtered and dried in an oven at 80 • C for 24 h. Then, 2.4 g of white crystals were dispersed in a mixture solution containing 50 mL methanol and 50 mL DMF. After that, the solution was added to a mixture solution containing 3.3 g of 2-aminoterephthalic acid and 75 mL of DMF. The final mixture was heated in an oven at 100 • C for 24 h. The resulting suspension was centrifuged at 7800 rpm for 5 min. Afterwards, the NH 2 -MIL-125 (Ti) fillers were washed several times with methanol and DMF before being dried in an oven at 60 • C for 24 h.

Fabrication of Dense Membrane
In our previous work [3], the 7.0 wt% NH 2 -MIL-125 (Ti) /6FDA-durene mixed matrix membrane showed the highest performance in terms of CO 2 and CH 4 single gas permeation; therefore, this membrane was utilized to obtain the optimum condition in CO 2 /CH 4 binary gas separation. The mixed matrix membrane was prepared by dispersing fillers and dissolving the polymer in DCM separately. After stirring, the suspension was sonicated to disperse the fillers in DCM. Then, the priming step was introduced by adding 10 wt% of the polymer solution into the filler suspension under stirring to induce the polymer-filler interface [22] before the suspension was sonicated again. The remaining polymer solution was added to the suspension, which was then sonicated. Subsequently, the suspension was vigorously stirred for 1 h before being cast onto a petri dish. The petri dish was then covered with perforated aluminum foil for solvent evaporation at room temperature for 24 h. After complete evaporation, the membrane film was carefully removed from the petri dish and dried in an oven at 60 • C. Then, the membrane was further dried under vacuum at 60 • C for 24 h, followed by thermal annealing at 250 • C for 24 h.

Binary Gas Separation Measurement
Binary gas separation testing was carried out by using a custom-made permeation test rig as described in our previous work [3]. The compositions of gases in the retentate and permeate gas streams were determined using a gas analyzer (Fuji, NDIR Gas Analyzer ZPAJ). Then, the permeability of the gases was calculated using Equations (1) and (2), as follows [20]: where P CO2 and P CH4 are the permeability of CO 2 and CH 4 (Barrer, 1 Barrer = 1 × 10 −10 cm 3 (STP).cm/s·cm 2 ·cmHg), V p is the volumetric flowrate (cm 3 (STP)/s), A m and t are the operational area (cm 2 ) and the thickness (cm) of the mixed matrix membrane, respectively, p h and p l are the pressure of feed and permeate side, respectively (cmHg), and x and y represent the volume fraction of the component in retentate and permeate streams by assuming ideal gas behaviors. The CO 2 /CH 4 separation factor was calculated using Equation (3), as follows [20]: where α CO2/CH4 indicates the CO 2 /CH 4 separation factor, and x and y represent the volume fraction of the component in the retentate and permeate streams, respectively.

Optimization of Membrane Separation Performance
The optimization of CO 2 separation from CH 4 was conducted based on the experimental runs generated by Design-Expert software V9.0 (Stat-Ease Inc., Minneapolis, MN, USA). The central composite design (CCD) available in the software was selected due to its various advantages, including its flexibility, efficiency, and continuous run. In the current study, three operational parameters were selected, namely, feed pressure (A), temperature (B), and CO 2 feed concentration (C). Meanwhile, three significant responses were determined, i.e., CO 2 permeability, CH 4 permeability and CO 2 /CH 4 separation factor. All parameters were tested at three measurement levels from low to high range. Then, the alpha (α) was fixed at 1, which is considered to be a face-centered design. The range and level of the operational parameters are shown in Table 1. Table 1. Experimental range and levels of the operational parameters.

Operational Parameters Units
Coded Measurement Levels The interaction between the operational parameters, i.e., feed pressure, temperature and CO 2 feed concentration, and responses, i.e., the permeability of the gases and the CO 2 /CH 4 separation factor, were further analyzed via analysis of variance (ANOVA). The experimental results were fitted into empirical models (second-order polynomial function) expressed in Equation (4) in order to correlate the corresponding responses over the operational parameters [18].
where Y is the corresponding response, B 0 is a fixed term, and B i , B ij and B ii are linear, interaction and quadratic terms, respectively. Meanwhile, x i and x j represent the coded terms for operational parameters.
The fitted model for each corresponding response was analyzed in terms of the statistical significance of the operational parameters and interactions using the F-and p-values. The F-value is defined as the ratio of the mean square model to the mean square error. The mean square model is the sum of squares divided by the number of degrees of freedom. The p-value indicates the significance of the model; therefore, p < 0.05 shows that the regression model is relevant, and the null hypothesis (H o ): p ≤ 0.5 is rejected. In addition, the consistency of the models was quantified by the coefficient of determination (R 2 ). Residual is the unexplained variation by the fitted model. 'Lack of fit' denotes that the lack of fit is insignificant in comparison to the pure error.

Characteristics of Membrane
The physicochemical properties, morphology and filler distribution, thermal stability, fractional free volume (FFV) analysis as well as CO 2 and CH 4 single gas performance of the NH 2 -MIL-125 (Ti)/6FDA-durene membranes at various filler loadings were discussed in detail in our previous work [3]. As noted in our previous work, the 7.0 wt% NH 2 -MIL-125 (Ti)/6FDA-durene mixed matrix membrane exhibited the highest CO 2 permeability and CO 2 /CH 4 gas pair selectivity, i.e., 1115.7 Barrer and 37.1, respectively, at a feed pressure of 3.5 bar and at room temperature. In this work, the performance of this membrane in mixed gas separation with various operating parameters, particularly at high CO 2 feed concentration, was investigated using a statistical approach consisting of a central composite design (CCD) paired with a response surface methodology (RSM). Figure S1 and Table S1 display the characterization results comprising the XRD, FESEM, EDX and FFV for the 7.0 wt% NH 2 -MIL-125 (Ti)/6FDA-durene mixed matrix membrane. In reference to Figure S1a, the significant peaks observed at 6.9 • and 9.8 • corresponded to the X-ray reflection planes (011) and (002) of the NH 2 -MIL-125 (Ti) fillers. Meanwhile, the fillers were encapsulated and uniformly dispersed in the 6FDA-durene polymer matrix, as shown in Figure S1b. On the other hand, Table S1 compares the CO 2 and CH 4 single gas permeation performance and FFV values of pure and 7.0 wt% NH 2 -MIL-125 (Ti)/6FDA-durene mixed matrix membrane.

Central Composite Design (CCD)
From CCD, 20 trial experiments, including eight factorial points, six axial points, and six central point replicates in a randomized sequence, were generated by the DoE software. The replicated center points were used to estimate the pure error for the lack of fit test. The randomized sequence for the experimental run was intended to limit the impact of uncontrollable factors. Therefore, experimental 20 runs were conducted and three significant responses, namely, CO 2 permeability, CH 4 permeability, and CO 2 /CH 4 separation factor, were measured. Table 2 shows the condition of the operational parameters and significant responses of the membrane separation performance observed in the experiments. From Table 2, it can be observed that the permeability of CO 2 and CH 4 were in the range of 293.2 to 794.4 Barrer and 28.7 to 147.5 Barrer, respectively. The percentage errors between the actual and predicted values of membrane separation performances are listed in Table S2. Meanwhile, the CO 2 /CH 4 separation factor obtained ranged from 5.6 to 13.0. The margin of error of the experimental results was ±5%. Furthermore, compared to single gas permeation testing [3], the performance of the membrane in binary gas separation demonstrated lower values. These difference could be explained by the competitive sorption between CO 2 and CH 4 gas molecules [23]. The presence of a competitive gas may substantially affect the gas penetration over the membrane [24].

CO 2 Permeability
The quadratic polynomial model suggested by the DoE software for CO 2 permeability is shown in Equation (5) where A, B and C denote the feed pressure (bar), temperature ( • C) and CO 2 feed concentration (mol%), respectively. Table 3 presents the findings of the ANOVA and regression analysis for CO 2 permeability over the membrane. As shown in Table 3, model F and p values of 27.11 and <0.05 were achieved, respectively, suggesting that the model is statistically significant. In this case, A, B, C, AB, AC, B 2 , C 2 are statistically significant model terms. It can be observed from Table 3 that an R 2 value of 0.96 was achieved, which validated the accuracy of the model for CO 2 permeability.  Figure 1a-c display a 3D plot of the effect of different operational parameters on CO 2 permeability. The increment of temperature at a feed pressure of 3.5 bar and a fixed CO 2 feed concentration of 42.5 mol% resulted in higher CO 2 permeability, as observed in Figure 1a. This could be due to the higher CO 2 diffusivity [19] caused by the increment of polymer chain mobility and free volume [18]. Meanwhile, the increment of feed pressure at temperatures between 30 to 40 • C resulted in a slight improvement in CO 2 permeability due to the increase in gas solubility as a consequence of thermodynamic promotion [11]. On the other hand, a decreasing trend of CO 2 permeability can be seen for temperatures ranging from 40 to 50 • C, mainly because of the CO 2 sorption isotherm, which follows the dual-mode sorption mechanism [25]. Generally, higher temperatures lead to higher CO 2 permeability; however, in this study, we observed a different trend. This might have been related to the lower solubility of gas molecules with increasing temperature. The adverse impact of temperature on sorption enthalpy would have affected the transport and sorption behavior of gases over the membrane [26]. From the obtained results, a maximum CO 2 permeability of 609.1 Barrer was achieved at a feed pressure of 3.5 bar and a temperature of 40 • C. Meanwhile, the lowest CO 2 permeability, i.e., 321.5 Barrer, was achieved at a feed pressure of 8 bar and a temperature of 30 • C. Figure 1b shows the effect of CO 2 feed concentration and feed pressure on CO 2 permeability at a constant temperature of 40 • C. It can be seen in Figure 1b that the lowest CO 2 permeability, i.e., 470.7 Barrer, was achieved at a CO 2 feed concentration of 15 mol% and a feed pressure of 8 bar. In contrast, the highest CO 2 permeability, 792.2 Barrer, was obtained at a feed pressure of 3.5 bar and a CO 2 feed concentration of 70 mol%. Additionally, at constant temperature, increasing the CO 2 feed concentration in the feed caused significant improvement in CO 2 permeability; this was mainly due to the higher sorption of CO 2 than CH 4 over the membrane. Furthermore, the interaction between CO 2 molecules and -NH 2 elements also contributed to the improvement of CO 2 diffusion through the membrane [27]. On the other hand, the CO 2 permeability decreased slightly with increasing feed pressure. This result was attributed to the saturation of available sorption sites, which resulted in a lower solubility coefficient [28]. Furthermore, this trend was also related to the dual-mode sorption and diffusion mechanism on gas transport behavior over the membrane. Figure 1c illustrates the effect of CO 2 feed concentration and temperature on CO 2 permeability at a fixed feed pressure of 8 bar. As shown in Figure 1c, increasing the CO 2 feed concentration led to an increase in CO 2 permeability due to the higher CO 2 sorption effect relative to CH 4 , which, subsequently, enhanced the sorption of CO 2 over CH 4 gas molecules in the membrane. A similar incremental trend of CO 2 permeability was observed under fixed feed pressure conditions, as shown in Figure 1b. When the temperature increased from 30 • C to 40 • C at a constant feed pressure of 8 bar, the CO 2 permeability improved due to the increase in movement and flexibility of polymer chains, as well as the kinetic energy of gas molecules [26]. Furthermore, it could be that the impact of temperature on the gas transport and the sorption behavior compensated for the adverse effect of temperature on gas solubility. However, a slight decrease in CO 2 permeability was discovered after increasing the temperature from 40 • C to 50 • C. This was mainly due to the sorption competition between the CO 2 and CH 4 gases [25]. The reduction in gas solubility caused by the reduced interaction of CO 2 gas with the polymer matrix was attributed to the increment of adsorption energy with temperature; consequently, the thermodynamic effect overcame the kinetic effect during gas penetration [29]. From Figure 1c, it may be seen that the maximum CO 2 permeability, i.e., 650.2 Barrer, was achieved at a feed pressure of 8 bar, a CO 2 feed concentration of 70 mol% and a temperature of 40 • C. A parity plot for CO 2 permeability is shown in Figure S2a. It can be seen that the actual and predicted values of the responses were scattered near to the 95% prediction limits.   where A, B and C denote the feed pressure (bar), temperature (°C) and CO2 feed tration (mol%), respectively. Table 4 demonstrates the ANOVA and regression analysis of CH4 permeab shown, F and p values of 29.95 and < 0.05, respectively, were achieved, indicating model terms were statistically significant. Additionally, A, B, C, AB, AC, BC, B were statistically significant model terms in this context. In addition, an R 2 valu was obtained, confirming the accuracy of the model for CH4 permeability.

CH 4 Permeability
Equation (6) where A, B and C denote the feed pressure (bar), temperature ( • C) and CO 2 feed concentration (mol%), respectively. Table 4 demonstrates the ANOVA and regression analysis of CH 4 permeability. As shown, F and p values of 29.95 and < 0.05, respectively, were achieved, indicating that the model terms were statistically significant. Additionally, A, B, C, AB, AC, BC, B 2 and C 2 were statistically significant model terms in this context. In addition, an R 2 value of 0.96 was obtained, confirming the accuracy of the model for CH 4 permeability.    Figure 2a. The increase in chain mobility and free volume resulted in higher diffusivity, and thus, increased the penetration of CH 4 molecules through the membrane [30]. The relationship between temperature and permeability can be defined using Arrhenius equation, where CH 4 molecules exhibit higher activation energy than CO 2 molecules, and therefore, promote the diffusion of nonpolar gases over the glassy polymer. Meanwhile, a reduction of CH 4 permeability was found with the increment of temperature from 40 • C to 50 • C. This was ascribed to the reduction in solubility, which subsequently enhanced the restriction of CH 4 gas permeation through the membrane [31]. Meanwhile, increasing the feed pressure caused a slight reduction in CH 4 permeability, mainly due to the decrease in the dual-sorption properties for CH 4 gas [19], as well as the high compressibility of CH 4 gas molecules [32] in the membrane. Figure 2b shows the effect of CO 2 feed concentration and feed pressure on CH 4 permeability at a constant temperature of 40 • C. Referring to Figure 2b, the CH 4 permeability obtained at feed pressures ranging from 3.5 bar to 12.5 bar demonstrated a similar trend to the results presented in Figure 2a. On the other hand, an increase in CO 2 feed concentration from 15 mol% to 42.5 mol% led to a reduction of CH 4 permeability, due to the lower sorption of CH 4 compared to CO 2 , as well as the slower diffusion of CH 4 gas over the membrane [20]. However, increasing the CO 2 feed concentration from 42.5 mol% to 70 mol% enhanced the permeability of CH 4 . The increasing trend in CH 4 permeability was because of the swelling of polymer chain packing at higher CO 2 feed concentrations, and therefore, increased the segmental mobility [19]. On the other hand, the lowest CH 4 permeability, i.e., 33.9 Barrer, was seen with a CO 2 feed concentration of 42.5 mol% and a feed pressure of 12.5 bar. In contrast, the maximum CH 4 permeability, 133.3 Barrer, was achieved at a CO 2 feed concentration of 70 mol%, a temperature of 40 • C and a feed pressure of 3.5 bar (Figure 2b). Figure 2c demonstrates the effect of CO 2 feed concentration and temperature on CH 4 permeability at a constant feed pressure of 8 bar. As shown, the CH 4 permeability gradually improved with increasing the temperature, owing to the enhanced mobility and flexibility of the polymer chains [26], as well as to the higher activation energy of CH 4 [33]. Meanwhile, a reduction in CH 4 permeability was observed at CO 2 feed concentrations ranging from 15 mol% to 42.5 mol%, and an improvement in CH 4 permeability was found when the CO 2 feed concentration increased from 42.5 mol% to 70 mol%. This trend was identical to that shown in Figure 2b. Apart from that, the highest CH 4 permeability, i.e., 106.1 Barrer, was obtained with a CO 2 feed concentration of 70 mol%, a feed pressure of 8 bar and a temperature of 40 • C (Figure 2c); this was associated with concentration polarization effects [34]. Meanwhile, the lowest CH 4 permeability, 28.7 Barrer, was observed at a CO 2 feed concentration of 42.5 mol% and a temperature of 30 • C. The increase in CH 4 permeability with the increase in temperature and CO 2 feed concentration was mainly due to the increase in polymer free volume caused by the alteration of polymer chain packing and intersegmental motion [33], as well as the swelling effect induced by CO 2 gas. The presence of high CO 2 feed concentrations caused the rapid diffusion of CH 4 gas through the membrane [35], because CH 4 gas molecules are more accessible to a swollen polymer matrix. A parity plot for CH 4 permeability is shown in Figure S2b. As shown, very few data points fall outside the 95% prediction limits.
Polymers 2022, 14, x FOR PEER REVIEW polymer matrix. A parity plot for CH4 permeability is shown in Figure S2b. As very few data points fall outside the 95% prediction limits.

CO2/CH4 Separation Factor
The quadratic polynomial model for CO2/CH4 separation factor is defined tion 7 in the form of coded values. where A, B and C denote the feed pressure (bar), temperature (°C) and CO2 feed tration (mol%), respectively.
The ANOVA analysis of the CO2/CH4 separation factor is presented in Table be seen that model F and p values of 18.04 and <0.05 were achieved, respectively ing that the model terms are statistically significant, while A, B, C, BC and C 2 are cally significant model terms. Furthermore, an R 2 value of 0.94 was achieved, sho accuracy of the model for CO2/CH4 separation factor.

CO 2 /CH 4 Separation Factor
The quadratic polynomial model for CO 2 /CH 4 separation factor is defined in Equation (7) in the form of coded values.
where A, B and C denote the feed pressure (bar), temperature ( • C) and CO 2 feed concentration (mol%), respectively. The ANOVA analysis of the CO 2 /CH 4 separation factor is presented in Table 5. It can be seen that model F and p values of 18.04 and <0.05 were achieved, respectively, indicating that the model terms are statistically significant, while A, B, C, BC and C 2 are statistically significant model terms. Furthermore, an R 2 value of 0.94 was achieved, showing the accuracy of the model for CO 2 /CH 4 separation factor.   Figure 3 displays 3D plots of the effect of feed pressure, temperature, and CO 2 feed concentration on the CO 2 /CH 4 separation factor over the membrane. As shown, increasing the feed pressure resulted in a higher CO 2 /CH 4 separation factor owing to the higher CO 2 adsorption capacity compared to CH 4 , as well as the good compatibility of the polymerfiller, which was induced by the presence of -NH 2 groups in the membrane [27]. The enhanced CO 2 /CH 4 separation factor was related to the improvement of CO 2 permeability and the reduction of CH 4 permeability observed in Figures 1a and 2a. In contrast, increasing the temperature revealed a negligible impact on CO 2 /CH 4 separation factor, since the contribution of solubility and diffusivity is interchangeable. Additionally, as shown in Figure 3a, the highest CO 2 /CH 4 separation factor, i.e., 13.0, was obtained at a temperature of 30 • C and a feed pressure of 12.5 bar. In comparison, the lowest CO 2 /CH 4 separation factor, 7.3, was observed at a feed temperature of 50 • C and a feed pressure of 3.5 bar. Figure 3b illustrates the effect of CO 2 feed concentration and feed pressure on the CO 2 /CH 4 separation factor at a constant temperature of 40 • C. It can be observed that the CO 2 /CH 4 separation factor increased with increasing feed pressure. This was mainly due to the higher sorption and permeation of CO 2 molecules over the membrane compared to CH 4 [33], which resulted in higher CO 2 permeability and lower CH 4 permeability, as observed in Figures 1b and 2b, respectively. The higher CO 2 permeability also contributed to the increase in CO 2 solubility and the reduction of CH 4 solubility, and thus, improved the CO 2 /CH 4 solubility selectivity [36]. Furthermore, the CO 2 /CH 4 separation factor displayed an increasing trend with increasing CO 2 feed concentration up to 42.5 mol%, and then decreased when the CO 2 feed concentration was increased further to 70 mol%. The increment of CO 2 feed concentration caused a reduction of CH 4 permeability and led to an increase in CO 2 /CH 4 separation factor, owing to the smaller kinetic diameter of CO 2 compared to CH 4 . Meanwhile, increasing the CO 2 feed concentration further from 42.5 mol% to 70 mol% caused an improvement of CH 4 permeability, but a reduction in the CO 2 /CH 4 separation factor was observed due to the early stage of CO 2 -induced plasticization behavior [33,37]. The highest separation factor, i.e., 13.0, was achieved at a feed pressure of 12.5 bar and CO 2 a feed concentration of 42.5 mol%. Meanwhile, the lowest CO 2 /CH 4 separation factor, 6.0, was observed at a feed pressure of 3.5 bar and a CO 2 feed concentration of 15 mol%, mainly resulting from the presence of a larger amount of CH 4 gas (85 mol%) in the feed mixture. Figure 3c shows the effect of CO 2 feed concentration and temperature on the CO 2 /CH 4 separation factor at a pressure of 8 bar. From Figure 3c, increasing the temperature at a fixed feed pressure of 8 bar caused a slightly drop in the CO 2 /CH 4 separation factor. This reduction was mainly due to the improvement of CO 2 and CH 4 permeability, as explained in Sections 3.3 and 3.4. Furthermore, at a fixed feed pressure of 8 bar, the increase of CO 2 feed concentration from 15 mol% to 42.5 mol% resulted in the enhancement of the CO 2 /CH 4 separation factor, which was consistent with the increase of CO 2 permeability, as shown in Figure 1c. On the other hand, increasing the CO 2 feed concentration from 42.5 mol% to 70 mol% led to a slight drop in the CO 2 /CH 4 separation factor due to the saturation of CO 2 gas molecules inside the polymer voids [11]. A parity plot for CO 2 /CH 4 separation factor is shown in Figure S2c. Overall, the actual and predicted values of the responses were scattered near the 95% prediction limits. As shown in Figure S2c, very few data points fall outside the 95% prediction limits.

Optimization of CO2/CH4 Separation Performance
The primary objective of this work is to identify the optimal operational co for membrane separation in terms of CO2 permeability and CO2/CH4 separatio Therefore, the optimization conditions included maximizing the CO2 feed conc as an operational parameter, as well as CO2 permeability and CO2/CH4 separatio Additionally, the desirability function (DF) was used in RSM to optimize a seq quadratic models. The geometric mean of individual desirability (d) was used to the total desirability (D), as expressed in Equation 8 [38]: where Dd is the total desirability and dn is the nth desirability, n = 1,2, …, n. T desirability was measured from 0 to 1, where 0 represents the most undesirable and 1 the most desirable. The optimum solution and its desirability generated by the DOE software is Table 6 in term of actual values. The optimum operational parameters suggested feed concentration, feed pressure, and temperature were 70 mol%, 12.5 bar and respectively, which yielded optimal CO2 permeability and separation factor va 571.9 Barrer and 11.9, respectively.

Optimization of CO 2 /CH 4 Separation Performance
The primary objective of this work is to identify the optimal operational conditions for membrane separation in terms of CO 2 permeability and CO 2 /CH 4 separation factor. Therefore, the optimization conditions included maximizing the CO 2 feed concentration as an operational parameter, as well as CO 2 permeability and CO 2 /CH 4 separation factor. Additionally, the desirability function (DF) was used in RSM to optimize a sequence of quadratic models. The geometric mean of individual desirability (d) was used to calculate the total desirability (D), as expressed in Equation (8) [38]: where D d is the total desirability and d n is the nth desirability, n = 1, 2, . . . , n. The total desirability was measured from 0 to 1, where 0 represents the most undesirable response and 1 the most desirable.
The optimum solution and its desirability generated by the DOE software is listed in Table 6 in term of actual values. The optimum operational parameters suggested for CO 2 feed concentration, feed pressure, and temperature were 70 mol%, 12.5 bar and 34.7 • C, respectively, which yielded optimal CO 2 permeability and separation factor values, i.e., 571.9 Barrer and 11.9, respectively.

Validation of the Optimum Condition
In order to validate the optimum conditions and predicted responses shown in Table 6, three repeated experiments were conducted based on the suggested conditions; the results are presented in Table 7. The percentage (%) error for CO 2 permeability ranged from 3.6% to 6.5% relative to the standard deviation of 1.5%. On the other hand, the percentage (%) error for CO 2 /CH 4 separation factor ranged from 0.8% to 5.0% with a standard deviation of 2.1%. Meanwhile, the average errors obtained for CO 2 permeability and CO 2 /CH 4 separation factor were 5.3% and 2.8%, respectively. Overall, the average error for the experimental and predicted values were within 5%, indicating that the model validity reached 95% of the prediction interval. Consequently, the model was successfully validated, and as such, the optimization of operational parameters was achieved utilizing the RSM approach.

Conclusions
In conclusion, the effects of operational parameters, i.e., feed pressure, temperature and CO 2 feed concentration, were identified as dominant factors influencing the separation performance of a NH 2 -MIL-125 (Ti)-6FDA/durene membrane. CO 2 feed concentration demonstrated a significant effect on CO 2 permeability, whereas feed pressure was the primary parameter influencing CH 4 permeability and CO 2 /CH 4 separation factor, with F-values of 102.5 and 59.3, respectively. The R 2 values obtained ranged from 0.94 to 0.96, indicating that the regression models were statistically significant. The optimum operational parameters, i.e., a feed pressure of 12.5 bar, a temperature of 34.7 • C and a CO 2 feed concentration of 70 mol%, yielded the maximum CO 2 permeability, i.e., 609.3 Barrer, and a CO 2 /CH 4 separation factor of 11.6. The average errors for CO 2 permeability and CO 2 /CH 4 separation factor were 5.3% and 2.8%, respectively, suggesting that the model was 95% reliable. Overall, the experimental findings show that the RSM paired with the CCD method is a better strategy to obtain optimal CO 2 /CH 4 separation performance with a NH 2 -MIL-125 (Ti)/6FDA-based mixed matrix membrane. The present research provides an experimental reference for further improvements and scale-ups of membranes in CO 2 /CH 4 separation with high CO 2 feed concentrations.

Supplementary Materials:
The following supporting information can be downloaded at https: //www.mdpi.com/article/10.3390/polym14071371/s1; Figure S1: The structural properties of the 7.0 wt NH 2 -MIL-125 (Ti)/6FDA-durene membrane including XRD, FESEM and EDX; Table S1: The FFV values and single gas permeation performance of membranes; Figure S2: The parity plot of predicted and actual data for the model of (a) CO 2 permeability, (b) CH 4 permeability, and (c) CO 2 /CH 4 separation factor with a 95% prediction interval.; Table S2: The actual and predicted values of the membrane separation performances with percentage of error.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.