Boosting CO2 adsorption and selectivity in metal–organic frameworks of MIL-96(Al) via second metal Ca coordination

Aluminum trimesate-based MOF (MIL-96-(Al)) has attracted intense attention due to its high chemical stability and strong CO2 adsorption capacity. In this study, CO2 capture and selectivity of MIL-96-Al was further improved by the coordination of the second metal Ca. To this end, a series of MIL-96(Al)–Ca were hydrothermally synthesised by a one-pot method, varying the molar ratio of Ca2+/Al3+. It is shown that the variation of Ca2+/Al3+ ratio results in significant changes in crystal shape and size. The shape varies from the hexagonal rods capped in the ends by a hexagonal pyramid in MIL-96(Al) without Ca to the thin hexagonal disks in MIL-96(Al)–Ca4 (the highest Ca content). Adsorption studies reveal that the CO2 adsorption on MIL-96(Al)–Ca1 and MIL-96(Al)–Ca2 at pressures up to 950 kPa is vastly improved due to the enhanced pore volumes compared to MIL-96(Al). The CO2 uptake on these materials measured in the above sequence is 10.22, 9.38 and 8.09 mmol g−1, respectively. However, the CO2 uptake reduces to 5.26 mmol g−1 on MIL-96(Al)–Ca4. Compared with MIL-96(Al)–Ca1, the N2 adsorption in MIL-96(Al)–Ca4 is significantly reduced by 90% at similar operational conditions. At 100 and 28.8 kPa, the selectivity of MIL-96(Al)–Ca4 to CO2/N2 reaches up to 67 and 841.42, respectively, which is equivalent to 5 and 26 times the selectivity of MIL-96(Al). The present findings highlight that MIL-96(Al) with second metal Ca coordination is a potential candidate as an alternative CO2 adsorbent for practical applications.


Introduction
Global warming is considered a serious disaster facing our planet. MOFs have been studied extensively for capturing greenhouse gases from gaseous mixtures. 1 Capture of these gases by adsorption techniques is practically used in industries. Conventional porous materials such as activated carbon, 2 zeolites 3 and metal oxides 4 are reliable adsorbents to control emission of greenhouse gases into the atmosphere. Recently, metal-organic frameworks (MOFs) have been a hot topic due to their great potential in different industrial applications 5 including greenhouse gas storage and separation. 6,7 In comparison with conventional materials, these MOF materials have a limited thermal and chemical stability; but possess larger surface areas with tunable pore sizes and pore volumes and facile functionalisation, making them attractive alternative adsorbents to be applied in various environmental elds. MOFs are synthesised by coordinating multidentate organic linkers with transition metal ions (or their clusters) into periodic porous frameworks. 8 Popularly, single metal-organic frameworks are based on a transition metal such as Al 3+ , 9 Fe 3+ , 10 Cr 3+ , 11 V 3+ , 12 Zn 2+ , 13,14 Mg 2+ , 15 Mn 2+ , 16 Co 3+ , 7 Cu 2+ , 17 Zr 4+ 18 or Ti 4+ . 19 MOFs are distinguished as porous materials, which can be easily modied to suit various applications. [20][21][22] One of these important modications is to synthesise MOFs with a diversity of metals accommodated in their structure. 23 This modication can be accomplished using two methods called direct synthesis modication (DSM) and post synthesis modication (PSM) for the synthesis of mixed metal-organic frameworks (M-MOFs). 24 DSM aims to enhance the textural characteristics and active functionalities of the parent MOFs by adding the second metal and other chemicals in one pot during the synthetic process. 25 However, DSM is very limited in preparation of M-MOFs because of the unfavourable existence of the central and second metal in the same reaction pot, which may result in unstable frameworks with weak topologies and functional groups. 26 To overcome this problem, suitable metals of similar ionic radii and coordination geometry are considered to maintain the integrity of the structures in M-MOFs during the DSM synthesis. On the other hand, PSM is mostly used in the preparation of M-MOFs for different applications. 27 This modication can be achieved by controlling the activation procedure, which can build M-MOFs with well-constructed features. 28 The main limitation of this method lies in unmatched physical or chemical properties between the main and the second metals. 29 Furthermore, due to the small aperture size and the high steric hindrance within the structures of MOFs, it is difficult to achieve the desirable metal molar ratios by PSM or to attain a uniform distribution of the two metal ions in the nal M-MOF, which may lead to undesirable alternation or collapse of the nal M-MOF structure. 30 Thus DSM can be more reliable to synthesise M-MOF if the selection of a second metal is carefully controlled according to the research target. Likewise, it has been conrmed that the preparation of M-MOFs by using incompatible metals in the direct synthesis results in the formation of a similar structure of the parent MOFs with a very low content of the second metal, but the textural properties and vacant metal sites are enriched. 31 Recently, M-MOFs have been synthesised to improve the characteristics of MOFs for use in specic applications such as catalysis, 32 sensing, illumines, 33 and gas storage. 34 In addition, CO 2 has been a major component of greenhouse gases that have been hugely emitted from fossil fuel-red power plants. 35 CO 2 can be effectively adsorbed and separated from other gas (such as N 2 ) by MOFs or M-MOFs with high selectivity. 36 In particular, the direct synthesis of M-MOFs leads to an increase in the concentration of open metal sites or defects to further increase the CO 2 uptake and selectivity. 37 Furthermore, the selectivity of CO 2 /N 2 is also governed by textural properties and functionalities of MOFs. 38 Nowadays, many researchers attempt to modify the recently developed MOFs to enhance the storage capacity of CO 2 and simultaneously increase its separation factor from other gases. 39 However, there are few studies dedicated to using the second metal to modify MOFs to improve their CO 2 adsorption. This is the major motivation of the present study. Recently, a microporous Al trimesate-based MOF, denoted MIL-96-(Al), has attracted intense attention because of its high hydrothermal stability and strong CO 2 adsorption capacity caused by the good affinity for CO 2 due to the presence of Al Lewis acid sites and bridging -OH groups in the framework. 40 In this work, we further improved the performance of MIL-96(Al) for CO 2 uptake with enhanced selectivity of CO 2 /N 2 by modifying its textural properties via the introduction of the incompatible metal Ca into the framework. For this purpose, DSM was applied to synthesise MIL-96(Al)-Ca samples using Al 3+ as the main metal and Ca 2+ as the second one. A series of MIL-96(Al)-Ca samples with different ratios of Ca 2+ /Al 3+ were synthesised and characterised. Aer Ca coordination, the MIL-96(Al)-Ca samples were activated by methanol to exchange most coordinated Ca 2+ and to enhance the concentration of vacant metal sites and defects in the molecular structure of MIL-96(Al)-Ca samples. Subsequently, these specially tailored bimetallic MOFs were tested for CO 2 and N 2 adsorptions. The present results indicate that MIL-96(Al) with a low Ca content can be used as an excellent adsorbent for capturing CO 2 ; by contrast, MIL-96(Al) with high Ca content can be applied as a novel adsorbent for separating CO 2 from N 2 .

Synthesis procedure
MIL-96(Al) and MIL-96(Al)-CaN (N ¼ 1, 2, 3, 4) were synthesised according to a previously reported typical procedure. 41 Al (NO 3 ) 3 $9H 2 O (19.69 mmol, 7.39 g) and benzene-1,3,5-tricarboxylic acid (BTC, C 9 H 6 O 6 ) (6.32 mmol, 1.33 g) were mixed with 26.93 mL of deionised water inside a 125 mL Teon container for around 40 min. Aer that, different amount of calcium carbonate (CaCO 3 ) was added inside the above mixture according to the ratios given in Table S1. † Next, the Teon container was capped and assembled into the steel case with tightly sealing and heated in a preheating oven at 493 K for 48 h. Aer that, the autoclave was le to cool at the laboratory temperature. Then, the vacuum ltration was done, and the crystalline collected product was washed thoroughly with deionised water. The splashed product was dried under vacuum at laboratory temperature and further dried at 373 K in an oven. The as-synthesised samples were activated by solvent exchange activation as described in the ESI. † This method is very active with using methanol as the exchangeable solvent.

Characterisation
Powder X-ray diffraction (XRD) measurements were taken using a Bruker D8 advance X-ray diffractometer with a Cu-Ka radiation source. A Fourier-transform infrared spectroscopy (FTIR) spectrometer (PerkinElmer) was used to nd out Fourier transformer infrared spectra in a wave number ranged from 650 to 4000 cm À1 . Thermal behaviour of the samples was determined by a thermogravimetric analysis (TGA) instrument (TGA/DSC1 STARe system-METTLER TOLEDO). The morphological description was achieved by using Zeiss Neon 40EsB FESEM. Elemental analysis of the samples was achieved via the elemental analysis instrument (ICP-OES). More characterisation details can be referred from the ESI section. †

Adsorption study
The CO 2 (Laser Grade, 99.995%), He (UHP Grade, 99.999%) and N 2 (UHP Grade, 99.999%) gases were supplied from BOC limited in Australia. Tristar-Plus-3020 and ASAP2050 (Micrometrics instruments, USA) were used for CO 2 adsorption analysis at atmospheric pressure and high pressure respectively. While N 2 adsorption was measured by ASAP2050 at 273 K and 290 kPa. An activated sample (0.09-0.15 g) was rstly dried in the oven for at least 1 h then the empty sample tube was weighed. Aer that, the dried sample was transferred into a weighed tube. Next, the lled tube was heated at 473 K under a vacuum for 8 h by a sample preparation system (VacPrep 061). Then, the net weight of the degassed sample was calculated. Finally, the degassed sample was analysed by the above instruments using a suitable analysing method. In addition, crushed ice was used to adjust the bath temperature at 273 K while the heating mantle of ASAP2050 was used to adjust other temperatures.  Table S1. † The trace Ca content in the sample is possibly due to the Ca leaching out from the MOF structure during the methanol exchange activation process. The X-ray powder diffraction (XRPD) patterns of MIL-96(Al) and MIL-96(Al)-Ca samples are shown in Fig. 1a. The patterns demonstrate a high similarity in the MIL-96(Al) structure as reported previously. 42 However, as the molar ratios of Ca 2+ /Al 3+ in MIL-96(Al)-Ca3 and MIL-96(Al)-Ca4 increase, some defects are observed. As shown in the pattern of MIL-96(Al)-Ca3, the positions of the peaks overlap with the peaks in other samples but their intensity is signicantly reduced, which is the prominent evidence of reducing the crystal size to nanometer range. 43 Further justication could be that the second metal would cause the defects if it is incompatible with the main metal as these two metals have different coulombic charges and ionic sizes. 44 FTIR spectra of MIL-96(Al)-Ca samples and MIL-96(Al) are shown in Fig. 1b. It can be seen that the infrared band at 1630 and 1332 cm À1 is allocated to the OH group bending in the plane of MIL-96(Al), MIL-96(Al)-Ca1 and MIL-96(Al)-Ca2. While this band disappears in the spectra of MIL-96(Al)-Ca3 and MIL-96(Al)-Ca4 of higher Ca 2+ content. This observation might be attributed to the further interaction of OH À bending with Ca 2+ in the metal centre instead of its interaction with OH group of non-coordinated BTC molecules inside the pores. 45 The spectrum of carboxylic groups in the BTC is usually seen in the spectral range from 1200 to 1700 cm À1 . Specically, protonated BTC restricted in the structure is detected at the band of 1684 cm À1 in the MIL-96(Al) and MIL-96(Al)-Ca1 while this is merged with the band of carboxylic group in the deprotonated BTC in other samples of higher Ca 2+ content. Likewise, the vibrations of asymmetric and symmetric stretching COO of Al-carboxylate groups in the deprotonated BTC at the activated MIL-96(Al)-Ca and the parent sample are displayed in two pairs of the peaks. The rst one is at 1660 and 1456 cm À1 as shown in the spectra of MIL-96(Al)-Ca except that in the spectrum of MIL-96(Al)-Ca3 the rst peak in this pair is seen at 1624 cm À1 ; and the second pair is at 1598 and 1396 cm À1 in parent sample while these peaks shi to 1600 and 1400 cm À1 in MIL-96(Al)-Ca1 and MIL-96(Al)-Ca2; and shis to 1567 and 1399 cm À1 , and 1578 and 1395 cm À1 in MIL-96(Al)-Ca3 and MIL-96(Al)-Ca4 respectively. 41 The spectrum of MIL-96(Al)-Ca3 is more affected because its crystal sizes are very ne. Consequently, the intensities of the peaks are reduced and shied. Moreover, the peaks at 760 and 735 cm À1 are assigned to C-H out the plane of the BTC ring via the structures of the MOFs. 46 The vibration band at 687 cm À1 is assigned to O-Al-O or O-(Ca) Al-O. This peak is attenuated as the Ca 2+ content increased. It seems that the vibrational motion of Al centre is changed because Ca 2+ occupies some active sites of aluminium cluster.

Results and discussions
Thermal gravimetric analysis (TGA) proles presented in Fig. 1c show that all samples have similar thermal stability. More specically, the thermal proles have exposed three steps of weight loss, the rst step being around 373 K, which can be attributed to evaporation of the moisture during heating to the boiling point of water. 47 The second step is caused by the burning of restricted molecules of the protonated BTC inside the pores. Also, this step might be affected by a coordinated Ca 2+ in the metal centres, which is clearly seen at around 600 K in the parent sample as shown in the previous studies. 41 As the Ca 2+ content in the main structures increase, this step shis to a higher temperature until it disappears in the thermal prole of MIL-96(Al)-Ca4 that conrms, higher content of Ca in the structure displaces the restricted free BTC out of the pores. Eventually, the third step represents the weight loss due to collapsing the whole structure at around 830 K because the connections between the linkers and the metal centres were broken.  Fig. 2a shows a hexagonal rod capped by two hexagonal bipyramids likes what the parent sample (MIL-96(Al)) displays in Fig. 2e and the length of its crystals have signicantly reduced. More specically, crystals of the parent sample had lengths ranged from 5 to 12 mm, and diameters ranged from 2 to 5 mm whereas the lengths of the crystals of MIL-96(Al)-Ca1 appear much shorter (in the range of 2 to 8 mm). Moreover, the general morphological description of MIL-96(Al)-Ca2 remains as the morphology of the parent sample, whereas the diameter of the rod signicantly reduces to 100 nm as shown in Fig. 2b. In addition, Fig. 2c shows that the shape of the crystals in MIL-96(Al)-Ca3 is a long rod with a very thin diameter (from 25 to 75 nm) lacking the hexagonal bipyramid on its ends. In contrast, Fig. 2d demonstrates a different morphology in MIL-96(Al)-Ca4 that is a thin hexagonal platelet (disk) with a diameter ranged from 10 to 20 mm and a thickness in a few nm. Remarkably, DSM can produce various morphological views and crystal sizes when the conditions of the synthesis are changed such as temperature, time and solvent. 48 In this study, using Ca 2+ as the second metal in the single pot synthesis caused the signicant change in the morphology because it is directly affected the crystal anisotropy. 49 This study conrms that the solvent exchange activation plays a main role in enhancing the textural properties by opening the metal sites of the prepared MIL-96(Al)-Ca samples. There are two reasons to justify this statement. One is the variety of the ionic sizes and valence numbers of Al 3+ and Ca 2+ and another is the efficiency of the activation process. Therefore, the second metal (Ca 2+ ) does not replace the main metal (Al 3+ ) as schematically shown in Fig. 3a and b, but it might coordinate with an active site of Al 3+ metal centre (Fig. 3b). The methanol exchange activation is a very reliable method to leach out most coordinated Ca 2+ by methanol leaving the opened metal sites for gas capture (Fig. 3c). However, the activated MIL-96(Al)-Ca samples still have traces of Ca 2+ in the nal products as presented in Table S1. † This journal is © The Royal Society of Chemistry 2020 RSC Adv., 2020, 10, 8130-8139 | 8133

CO 2 and N 2 adsorption behavior
The N 2 adsorption/desorption isotherms at 77 K, micropore distribution and mesopore distribution of the MIL-96(Al)-Ca samples compared to the parent sample are shown in Fig. 4ac. For each sample shown in Fig. 4a, a different hysteresis is obtained due to the presence of different pore sizes, including micropores and mesopores or their combined pores. Thus, it is expected that the slow desorption mechanism (percolation theory) might happen. 50 Moreover, as shown in Table 1, the macropore content is remarkably increased as the content of Ca 2+ is enhanced by more than 0.09%, and it is approved by a sharp increase of N 2 adsorption when the pressure approached the atmospheric pressure. This is similar behaviour to that of MIL-100(Fe)-Ca in our previous work. 31 In addition, a signicant change is observed in N 2 adsorption of MIL-96(Al)-Ca4 which exhibits a larger hysteresis loop than other samples due to the majority of the mesopore and macropore in its structure. Fig. 4b shows the micropore distribution of MIL-96(Al)-Ca 1, 2, 3 and 4 compared to the parent sample. The availability of the smallest pore is seen in the samples of lower Ca 2+ content, the pore diameter is 0.   The mesopore size distribution curves presented in Fig. 4c clearly show that both of the parent sample and MIL-96(Al)-Ca2 exposed very low cumulative volume (<0.0005 cm 3 g À1 ) in a single unique peak of pore size centred at 4 and 13 nm respectively. MIL-96(Al)-Ca2 demonstrates the largest mesopore diameter with the lowest in the cumulative volume. Whereas in MIL-96(Al)-Ca1, the maximum cumulative volume is doubled (0.001 cm 3 g À1 ) in smaller mesopore diameter of  3 nm, and that volume increases in MIL-96(Al)-Ca3 (0.0027 cm 3 g À1 ) and MIL-96(Al)-Ca4 (0.0013 cm 3 g À1 ) at mesopore diameter of 4 and 5 nm respectively. As a result, both population and the length of mesopore are signicantly enhanced in the samples of the higher Ca 2+ content. 51,52 Table 1 illustrates the BET surface areas, average pore sizes, pore volumes, micropore contents. Recognisably, the BET surface area increases at the lowest loading of Ca 2+ (0.09%) in MIL-96(Al)-Ca1, but it is declines at the highest Ca The pore diameter in MIL-96(Al)-Ca4 is signicantly higher than other samples. This can be interpreted by the expansion of a cage size, which is governed by the number of vertices and their connectivity. Therefore, the larger pore is possible to be dominated when the connectivity of vertices is reduced. 53 It means when some of the active Al 3+ sites are occupied by Ca 2+ , some of the vertices may lose some of their connectivity. Consequently, further enhancement in the pore size may occur when the Ca 2+ are desolvated via the activation process. However, the pore volume is increased to 0.3, 0.32 and 0.53 cm 3 g À1 in MIL-96(Al)-Ca1, MIL-96(Al)-Ca2 and MIL-96(Al)-Ca3 respectively, which is higher than the pore volume in MIL-96(Al) and reduced to 0.151 cm 3 g À1 in MIL-96(Al)-Ca4 due to its lowest surface area. Finally, the microporosity is signicantly decreased when the concentration of Ca 2+ is raised to 1.4%. Hence, the micropore content is dropped from 94% in the parent sample to 8.4% in MIL-96(Al)-Ca4. Fig. 5a shows the CO 2 adsorption at 273 K and high pressure up to 950 kPa. It is noticed that the CO 2 adsorption capacity is 8.09 mmol g À1 in the parent sample. Whereas this noticeably enhanced in MIL-96(Al)-Ca1, MIL-96(Al)-Ca2 and MIL-96(Al)-Ca3 as it is 10.22, 9.38 and 8.64 mmol g À1 respectively. However, CO 2 adsorption decreases to about 5.26 mmol g À1 in MIL-96(Al)-Ca4 because further increasing in Ca 2+ content leads to decrease in the BET surface area, pore-volume, micropore content. 54 It seems that the pore volume expands when the high pressure is applied. Therefore, CO 2 uptake is dramatically increased in the MIL-96(Al)-Ca samples with the Ca 2+ content lower than 0.74%. Systematically, the CO 2 adsorption in microporous materials increases when the pore volume and BET surface area are augmented and vice versa. 55 It can be conrmed from Table 2 that the adsorption of CO 2 is signicantly boosted in most of the MOFs in this work and typically in MIL-96(Al)-Ca1 and MIL-96(Al)-Ca2. Although, almost all samples in the tabulated references have higher surface areas, they present lower CO 2 adsorption capacity than the present work listed samples. Fig. 5b shows that the isosteric heat of CO 2 adsorption changes depending on the CO 2 coverages in MIL-96(Al)-Ca samples. The heat of adsorption is reduced with increasing the Ca 2+ content. Generally, in the sample of the lowest Ca 2+ content, the heat of adsorption slightly changes when the coverage of CO 2 increases. More specically, MIL-96(Al)-Ca1 and MIL-96(Al) have exposed similar heat of CO 2 adsorption as it is 23 kJ mol À1 at CO 2 coverage of 60 mmol g À1 . In contrast, the heat of CO 2 adsorption in MIL-96(Al)-Ca2 and MIL-96(Al)-Ca3 obviously decrease to 20 and 19 kJ mol À1 respectively. However, MIL-96(Al)-Ca4 shows a different behavior, the heat of adsorption is decreased further than that in MIL-96(Al)-Ca3 when the coverage of CO 2 ranges from 10 to 38 mmol g À1 , then it is increased over than both of MIL-96(Al)-Ca2 and MIL-96(Al)-Ca3 at a higher coverage of CO 2 . This heat of CO 2 adsorption is 18 kJ mol À1 at CO 2 coverage of 10 mmol g À1 and is 22 kJ mol À1 at CO 2 coverage of 55 mmol g À1 .
The adsorption of CO 2 and heat of adsorption are profoundly affected by functionalities and sizes of the pores. Furthermore, heterogeneity of the pores can be dominated in mesoporous materials. 56 As a result, the homogenous adsorption can be anticipated to occur in the microporous adsorbents such as MIL-96(Al), MIL-96(Al)-Ca1 and MIL-96(Al)-Ca2 due to the homogeneity of the ne pores. Accordingly, the heat of CO 2 adsorption is approximately uniform in those adsorbents.
The N 2 adsorption capacities of MIL-96(Al)-Ca samples and the parent sample at extended pressure up to about 290 kPa are illustrated in Fig. 5c. Similarly, to the CO 2 adsorption, the highest N 2 adsorption is seen in MIL-96(Al)-Ca1. This material adsorbs 1.2 mmol g À1 of N 2 , which is higher than that of the parent sample and other MIL-96(Al)-Ca samples. Remarkably, N 2 adsorption on MIL-96(Al)-Ca4 is very low (0.12 mmol g À1 ) due to the high average pore size. 57 However, when the materials are ultrane microporous as in the parent sample and MIL-96(Al)-Ca1 of smaller pores, the adsorption affinity of N 2 molecules towards the surfaces of the pores is dramatically enhanced. In other words, the molecules of N 2 can be actively interacted with the surface of the pores, as well as with each other due to the very limited space in the micropore.
In large pores, the interactions between the N 2 molecules themselves are not likely to happen at normal temperatures. Also, in these conditions, the N 2 molecules are weakly attached to the adsorption sites on the adsorbents because N 2 molecules have a low quadrupole moment and low polarisability in addition to absence of the intensive electric eld inside the large pores. 58 Therefore, the adsorption capacities of N 2 are extremely lower than those of CO 2 , which has much higher quadrupole moment and polarisability and lower kinetic diameter. The selectivity of CO 2 over N 2 at 100 kPa and 273 K is shown in Fig. 5d. In general, MIL-96(Al)-Ca4 presents signicant separation factor among other samples in this study. The selectivity of CO 2 /N 2 at 1 bar and 273 K is 13.46 in MIL-96(Al) while it is 11.68, 12, 13.08, and 60. It is clearly observed that a further increase in the pore size has signicantly participated in enhancing the selectivity of CO 2 /N 2 on MIL-96(Al)-Ca4. The improvement of this selectivity is because CO 2 has higher polarisability, a higher quadrupole moment and a lower kinetic diameter. These parameters can highly enhance the adsorption of CO 2 in the microporous and mesoporous materials with high superiority of adsorption in the ultrane-microporous adsorbents. 59 On the other hand, as it  was mentioned above, the N 2 molecule at STP conditions has very low affinity to be adsorbed in the porous materials with a majority of larger mesopore. 38,58 However, although the microporous content, BET and CO 2 adsorption is very low in MIL-96(Al)-Ca4, the affinity of this adsorbent for CO 2 is much higher than that of N 2 due to presence the Lewis base-calcium oxide coordinated on the metal center 60 . Thus, it produces a very high selectivity compared to other adsorbents reported in the literature as shown in Table 2.

Conclusions
Four samples of MIL-96(Al)-Ca and MIL-96(Al) were successfully synthesised, namely MIL-96(Al)-Ca1, MIL-96(Al)-Ca2, MIL-96(Al)-Ca3 and MIL-96(Al)-Ca4 with different Ca content. They were well characterised and used for CO 2 and N 2 adsorption. It was found that Ca 2+ content in the nal product is very low due to the incompatibility in the ionic size of Al 3+ and Ca 2+ and the high performance of the methanol exchange activation, but inuencing the CO 2 capture and selectivity signicantly. The results reveal that CO 2 adsorption is enhanced as the Ca 2+ content increases up to 0.09% while it signicantly dropped when the Ca 2+ content reaches 1.4%. Moreover, the N 2 adsorption sharply reduces when the Ca 2+ content increases in MIL-96(Al)-Ca3 and MIL-96(Al)-Ca4. Consequently, the selectivity of CO 2 /N 2 in MIL-96(Al)-Ca4 increases by a factor of 5 compared to MIL-96(Al). The present results suggest that MIL-96(Al)-Ca1 is a potential candidate for adsorption of CO 2 while MIL-96(Al)-Ca4 is an excellent adsorbent for separation of CO 2 from N 2 . These adsorbents may also have high potential usage for water treatment, deserving further investigations.

Conflicts of interest
The authors declare no conict of interest.