In Silico Screening of Metal−Organic Frameworks and Zeolites for He/N2 Separation

In silico screening of 10,143 metal−organic frameworks (MOFs) and 218 all-silica zeolites for adsorption-based and membrane-based He and N2 separation was performed. As a result of geometry-based prescreening, structures having zero accessible surface area (ASA) and pore limiting diameter (PLD) less than 3.75 Å were eliminated. So, both gases can be adsorbed and pass-through MOF and zeolite pores. The Grand canonical Monte Carlo (GCMC) and equilibrium molecular dynamics (EMD) methods were used to estimate the Henry’s constants and self-diffusion coefficients at infinite dilution conditions, as well as the adsorption capacity of an equimolar mixture of helium and nitrogen at various pressures. Based on the obtained results, adsorption, diffusion and membrane selectivities as well as membrane permeabilities were calculated. The separation potential of zeolites and MOFs was evaluated in the vacuum and pressure swing adsorption processes. In the case of membrane-based separation, we focused on the screening of nitrogen-selective membranes. MOFs were demonstrated to be more efficient than zeolites for both adsorption-based and membrane-based separation. The analysis of structure–performance relationships for using these materials for adsorption-based and membrane-based separation of He and N2 made it possible to determine the ranges of structural parameters, such as pore-limiting diameter, largest cavity diameter, surface area, porosity, accessible surface area and pore volume corresponding to the most promising MOFs for each separation model discussed in this study. The top 10 most promising MOFs were determined for membrane-based, vacuum swing adsorption and pressure swing adsorption separation methods. The effect of the electrostatic interaction between the quadrupole moment of nitrogen molecules and MOF atoms on the main adsorption and diffusion characteristics was studied. The obtained results can be used as a guide for selection of frameworks for He/N2 separation.


Introduction
Helium is widely used in aviation, space and electronic industries, science and medicine due to its unique physical properties [1][2][3]. Today helium is mostly recovered from natural gas by the cryogenic distillation method. The traditional process used for helium separation [3][4][5] consists of the following stages: preliminary purification (separation from C 2+ hydrocarbons, water, H 2 S and CO 2 ), cryogenic separation of methane followed by cryogenic separation of nitrogen in a nitrogen rejection unit (NRU) yielding feedstock with helium concentration 1-3%, separation of the N 2 /He mixture in a helium recovery unit yielding 50-70% crude helium vapor stream, purification of crude helium in a helium upgrade unit producing outcoming gas with helium concentration about 90%, deep purification in a helium purification unit producing helium with 99.995% purity and liquefaction of commodity helium. Membrane-based separation, vacuum swing adsorption (VSA) and pressure swing adsorption (PSA) processes are alternatives for the cryogenic distillation method [3,6,7]. Specific feature of these methods is that they can be used both IRMOF-1, which was equal to 2.4 [29]. Another representative of the IRMOF family IRMOF-3 demonstrated the selectivity of 2.5 [29]. In addition, the He/N 2 selectivity of membranes made of MIL-53 [30] and [Cu 2 (bza) 4 (pyz)] n [31] was reported to be equal to 2.4 and 3.9, respectively.
If the experimental measurement of the selectivity, adsorption capacity and permeability is not possible, simulation of adsorption and diffusion characteristics can give valuable information for estimation of the suitability of different materials for the He/N 2 separation. The development of zeolite database IZA [32] and MOF databases CSD MOF [33] and CoRE MOF [34] containing the unit cell parameters and atomic coordinates initiated in silico screening of efficient materials for adsorption-based [35][36][37][38][39] and membrane-based [40][41][42] separation. For example, Zarabadi-Poor et al. carried out in silico screening of almost 500 MOF structures from the DFT-optimized CoRE MOFs database [43] for helium separation from natural gas via N 2 /He separation [44]. The API (Adsorbent Performance Indicator) metric [45] and membrane selectivity were used to determine the top performing MOF structures in the case of adsorption-based and membrane-based separation, respectively. We used the results obtained in this work as a reference for comparison. Most of the studies are focused on the screening of materials for membrane-based helium recovery from natural gas [10,46]. Kadioglu et al. screened 139 MOF structures for He/CH 4 membrane-based separation [47]. It was shown that the top performing MOF structures had a pore limiting diameter in the range of 3.8-4 Å. Qiao and coworkers predicted gas permeability and selectivity of MOFs for membrane-based separation of He/CH 4 and He/N 2 at infinite dilution condition [48]. Based on the obtained results, the Top 5 MOF structures for He/N 2 separation show selectivity in the range of 3.02 to 3.53 and permeability in the range of 8.72 × 10 2 to 8.72 × 10 3 Barrer. Daglar et al. studied incorporation of MOF fillers into polymers to obtain mixed matrix membranes (MMM) for 11 different gas separations processes [49]. It was shown that all MOF fillers improve He permeability with slight changes in He/N 2 membrane selectivity. Helium permeability and He/N 2 selectivity of such MMM varied from 39.2 to 1.13 × 10 4 Barrers and from 0.8 to 622, respectively. Gas permeabilities and selectivities for covalent organic frameworks (COFs) from the CURATED COF database [50] and COF/polymer MMMs for helium separation were predicted in two studies by different groups (Aydin et al. [51] and Feng et al. [52]). It was shown that COFs had a linear correlation between He permeability and He/N 2 selectivity. For He/N 2 separation, the selectivity of COF structures varies from 5.8×10 -2 to 8.7. In both studies it was shown that the addition of COF fillers improved the He permeability of MMMs without significantly changing the He/N 2 selectivity. Based on the results of these studies, we think that MOF/COF membranes have low He selectivity at room temperature, while polymer membranes with high He/N 2 selectivity have low He permeability. Thus, two strategies are possible for creating efficient membranes for the He/N 2 separation: incorporation of porous fillers into the polymers to obtain MMMs with high He/N 2 selectivity and He permeability, or the use of nitrogen-selective porous medium membranes. Since the first approach has been extensively studied in the literature, we have focused on the screening of nitrogen-selective membranes. The goal of this study was to perform screening of zeolites and MOFs for search of relationships between their structural parameters and performance characteristics for adsorption-based and membrane-based He/N 2 separation.

Adsorption-Based Gas Separation
To test the force field models used in this study, calculated N 2 adsorption isotherms for a number of well-known MOFs at 295-298 K and pressure up to 20 bar were compared with the corresponding literature data [53][54][55][56][57][58][59]. Figure 1 (left) demonstrates that the suggested force field model can be used to predict the nitrogen adsorption in various MOFs with good precision. Several force fields were used to simulate the nitrogen adsorption in zeolites: Dreiding [60], TraPPE [61] and the force field developed by Vujić and Lyubartsev [62]. For each force field model, N 2 uptake values at 303 K and 1 bar pressure were calculated for five different pure silica frameworks, and the obtained results were compared with the corresponding experimental data [63] ( Figure S1). The model by Vujić and Lyubartsev was used hereafter because it predicts the experimental data with the best precision. Figure 1 (right) demonstrates that this force field makes it possible to predict nitrogen adsorption isotherms with good precision both on pure silica zeolites and on several aluminophosphates [64,65]. The experimental conditions used to measure the nitrogen adsorption isotherms and the corresponding references are reported in Tables S2 and S3. Testing the used force field models for helium adsorption simulation is very difficult because the value measured in the adsorption experiment is excess adsorption. However, before the experiment, dead space is measured using the same helium. So, it is not possible to measure the helium adsorption isotherm by traditional methods. suggested force field model can be used to predict the nitrogen adsorption in var MOFs with good precision. Several force fields were used to simulate the nitrogen ads tion in zeolites: Dreiding [60], TraPPE [61] and the force field developed by Vujić and bartsev [62]. For each force field model, N2 uptake values at 303 K and 1 bar pressure w calculated for five different pure silica frameworks, and the obtained results were c pared with the corresponding experimental data [63] ( Figure S1). The model by Vujić Lyubartsev was used hereafter because it predicts the experimental data with the precision. Figure 1 (right) demonstrates that this force field makes it possible to pre nitrogen adsorption isotherms with good precision both on pure silica zeolites an several aluminophosphates [64,65]. The experimental conditions used to measure th trogen adsorption isotherms and the corresponding references are reported in Table  and S3. Testing the used force field models for helium adsorption simulation is very ficult because the value measured in the adsorption experiment is excess adsorp However, before the experiment, dead space is measured using the same helium. So not possible to measure the helium adsorption isotherm by traditional methods. The performed screening of MOFs demonstrated that Henry's constants for He sorption are in the range from 1.25 × 10 −8 to 6.89 × 10 −6 mol/kg/Pa, and for N2 adsorp they vary from 7.08 × 10 −8 to 1.21 × 10 −3 mol/kg/Pa. The N2/He adsorption selectivi infinite dilution exceeds one for all the studied frameworks reaching a maximum valu 13,829. The highest Henry's constants for N2 adsorption ( Figure S4) and nitrogen ads tion selectivity ( Figure S5) were observed at PLD about 4.5 Å and LCD about 6 Å. result can be explained by the fact that the intermolecular interaction potentials ove in small cavities leading to an increase in the heats of adsorption and Henry's const The analysis of the relationship between structural parameters and Henry's constant He adsorption ( Figure S4) revealed a trend to the 0 He K growth with an increase in a sible surface area, porosity and accessible pore volume. Meanwhile, a trend to the decrease with an increase in the framework density was observed. In the case of K these trends were less evident due to the presence of an additional contribution from trostatic interaction between the quadrupole moment of the nitrogen molecule and pa atomic charges of the MOF atoms. The analysis of the correlation between the struc parameters and N2 adsorption selectivity ( Figure S5) showed a trend of the decrea 0 , 2/ ads N He S with an increase in accessible surface area and porosity, and a decrease in sity. Similar relationships between the structural parameters and adsorption selectivi nitrogen were also observed for zeolites ( Figure S5). The highest adsorption selectivi nitrogen calculated at infinite dilution for zeolites was as high as 26. A large variety o The performed screening of MOFs demonstrated that Henry's constants for He adsorption are in the range from 1.25 × 10 −8 to 6.89 × 10 −6 mol/kg/Pa, and for N 2 adsorption they vary from 7.08 × 10 −8 to 1.21 × 10 −3 mol/kg/Pa. The N 2 /He adsorption selectivity at infinite dilution exceeds one for all the studied frameworks reaching a maximum value of 13,829. The highest Henry's constants for N 2 adsorption ( Figure S4) and nitrogen adsorption selectivity ( Figure S5) were observed at PLD about 4.5 Å and LCD about 6 Å. This result can be explained by the fact that the intermolecular interaction potentials overlap in small cavities leading to an increase in the heats of adsorption and Henry's constants. The analysis of the relationship between structural parameters and Henry's constants for He adsorption ( Figure S4) revealed a trend to the K 0 He growth with an increase in accessible surface area, porosity and accessible pore volume. Meanwhile, a trend to the K 0 He decrease with an increase in the framework density was observed. In the case of K 0 N2 , these trends were less evident due to the presence of an additional contribution from electrostatic interaction between the quadrupole moment of the nitrogen molecule and partial atomic charges of the MOF atoms. The analysis of the correlation between the structural parameters and N 2 adsorption selectivity ( Figure S5) showed a trend of the decrease in S 0 ads,N2/He with an increase in accessible surface area and porosity, and a decrease in density. Similar relationships between the structural parameters and adsorption selectivity to nitrogen were also observed for zeolites ( Figure S5). The highest adsorption selectivity to nitrogen calculated at infinite dilution for zeolites was as high as 26. A large variety of the MOFs chemical composition leads to a greater variety of their structural parameters compared to zeolites ( Figure S3). As a result, the selectivity of the most promising MOFs in adsorption-based separation of nitrogen-helium mixtures exceeds the best of the zeolite frameworks by several orders of magnitude.
To move from the estimation of ideal adsorption properties of MOFs and zeolites at infinite dilution to gas mixtures, adsorption of an equimolar mixture of helium and nitrogen was simulated at 0.01, 0.1, 0.3 and 1 Mpa and 298 K. The obtained results ( Figure S6) demonstrate that at low adsorption selectivity S 0 ads,N 2 /He and S mix ads,N 2 /He are almost equal. However, when the adsorption selectivity increases, the S mix ads,N 2 /He /S 0 ads,N 2 /He ratio becomes less than one.
When the pressure increases, the difference between the two selectivity values grows as well. So, the higher the pressure in the system, the greater the effect of competitive adsorption. Two models of He/N 2 adsorption-based separation at room temperature were considered in this study. In the first model, corresponding to the conditions of vacuum swing adsorption, the adsorption pressure was equal to 0.1 Mpa and the desorption pressure was equal to 0.01 Mpa. In the second model, corresponding to conditions of pressure swing adsorption, the adsorption pressure was equal to 1 Mpa and the desorption pressure was equal to 0.1 Mpa. As regenerability tended to decrease with the APS increase, only frameworks with regenerability above 80% were considered. The highest APS values were observed for frameworks with high nitrogen ∆N and S mix ads,N 2 /He . Top 50 best MOFs have APS > 57 mol/kg and APS > 70 mol/kg for the VSA and PSA gas separation models, respectively. (Figure 2a,b). In the case of PSA, the APS metric of the most promising MOFs is higher by more than an order of magnitude than that of the most efficient zeolites (Figure 2d). In the case of VSA, this difference reaches two orders of magnitude ( Figure 2c). This result clearly demonstrates how promising MOFs are for adsorption-based separation of helium and nitrogen. In addition to the APS metric, sometimes a more complex API (Adsorbent Performance Indicator) metric [45], which additionally takes into account enthalpy of adsorption, is used in the literature: Here, ∆H ads,N2 is enthalpy of nitrogen adsorption; constants A, B and C are equal to 0.5, 2 and 1 [45]. As adsorption is an exothermal process, for large-scale industrial installations for separation based on vacuum swing adsorption, the use of an adsorbent with high adsorption enthalpy of results in an increase of the adsorber temperature leading to the decrease of the target component adsorption. On the other hand, heat is absorbed during the adsorbent regeneration leading to the decrease in the adsorber temperature, which makes the adsorbent regeneration more difficult. So, the higher the enthalpy of adsorption, the greater the difference of regenerability from the ideal model value. Therefore, an efficient adsorbent should have high adsorption capacity, high selectivity and low enthalpy of adsorption. In this study the API metric was calculated for all the studied MOFs and zeolites with regenerability above 80%. Figure S7 demonstrates a linear correlation between API and APS for both VSA and PSA. So, the use of both metrics leads to the same set of the most promising frameworks.
In this study, structure-adsorption performance relationships were investigated. Two databases with Top 50 MOFs based on the APS metric were constructed for the VSA and PSA separation processes. Then, smoothed probability density distributions (PDF) were built for several structural parameters for the database of all studied MOFs and for the databases of Top 50 MOFs for the VSA and PSA separation processes. These structural parameters included pore limiting diameter, largest cavity diameter, accessible surface area, accessible pore volume, density and porosity. The range of effective structural parameters was determined using a criterion that this range included more than 90% of Top 50 MOFs. The importance of each structural parameter was estimated by comparing PDFs for all studied MOFs and Top 50 MOFs. If PDFs were about the same for both databases, this structural parameter was considered to have little effect on the adsorption performance. Meanwhile, if PDFs were very different and PDF for TOP 50 MOFs had a narrow distribution, it was possible to claim that certain range of optimal parameters existed. For instance, Figure 3 demonstrates that for the VSA separation process there is a narrow range 3.75 Å < PLD < 4.8 Å corresponding to more than 90% of Top 50 MOFs, and PDF for them is different from PDF for all MOFs. So, the following optimal structural parameters were determined for VSA ( Figure S8): 3.75 Å < PLD < 4.8 Å, 4.4 Å < LCD < 6 Å, 100 m 2 /g < ASA < 700 m 2 /g, 0.02 cm 3 /g < AV < 0.09 cm 3 /g, 1300 kg/m 3 < density < 2500 kg/m 3 , 0.04 < VF < 0.15. For PSA the optimal ranges are wider ( Figure S9): 3.75 Å < PLD < 6.3 Å, 4.2 Å < LCD < 7.1 Å, 300 m 2 /g < ASA < 1400 m 2 /g, 0.03 cm 3 /g < AV < 0.18 cm 3 /g, 900 kg/m 3 < density < 2500 kg/m 3 , 0.05 < VF < 0.2. To determine the chemical composition-adsorption performance relationships, the probabilities of finding certain metal atoms in the MOF structure were calculated for all the studied MOFs and the database of Top 50 MOFs for the VSA and PSA separation processes ( Figure 4). Ga, Ru, V, Er, Gd, La, U and Ca were found to be the most suitable metals for the VSA separation model based on the probabilities of their presence among Top 50 MOFs and all studied MOFs. Meanwhile, the most widespread metal atoms Zn, Cu, Cd, Co and Mn (present in 54.7% of all MOFs) were scarcely present among Top 50 MOFs (less than 10%). In the case of the PSA model, U, Al, Er, Be and Mg were found to be the most suitable metals. Characteristics of the most promising MOFs (Top 10) for the VSA and PSA separation processes are reported in Tables S6 and S7. As the pressure increase leads to a much more significant selectivity decrease than the adsorption capacity increase, the best MOFs for VSA have higher APS than the most promising MOFs for PSA. In this study, structure-adsorption performance relationships were investigated. Two databases with Top 50 MOFs based on the APS metric were constructed for the VSA and PSA separation processes. Then, smoothed probability density distributions (PDF) were built for several structural parameters for the database of all studied MOFs and for the databases of Top 50 MOFs for the VSA and PSA separation processes. These structural parameters included pore limiting diameter, largest cavity diameter, accessible surface area, accessible pore volume, density and porosity. The range of effective structural pa- spread metal atoms Zn, Cu, Cd, Co and Mn (present in 54.7% of all MOFs) were scarcely present among Top 50 MOFs (less than 10%). In the case of the PSA model, U, Al, Er, Be and Mg were found to be the most suitable metals. Characteristics of the most promising MOFs (Top 10) for the VSA and PSA separation processes are reported in Tables S6 and S7. As the pressure increase leads to a much more significant selectivity decrease than the adsorption capacity increase, the best MOFs for VSA have higher APS than the most promising MOFs for PSA. It is interesting to compare the results of the MOF screening obtained in this study with earlier literature data. A total of 213 MOFs were studied for VSA gas separation with very similar operational conditions (equimolar mixture of helium and nitrogen, gas pressure during the adsorption cycle 1.2 bar, pressure during the adsorbent regeneration 0.1 bar) [44]. UVEXAV was found to be the most promising MOF with API = 680 and adsorption selectivity 2 , / mix ads N He S = 222.7. In our study we discovered 24 MOFs with superior API metric and 34 MOFs with higher nitrogen selectivity. Additionally, in the same earlier study [44] it was demonstrated that the electrostatic contribution for most MOFs was negligible. Despite the fact that a nitrogen molecule has a relatively a low quadrupole moment, we believe that it is important to consider the electrostatic interactions between nitrogen molecules and the framework atoms during the screening. To estimate their contributions, Henry's constants and nitrogen enthalpies of adsorption were calculated with the account of electrostatic interactions and without them ( Figure S10). The account of electrostatic interactions results in the growth of the median value of Henry's constants and nitrogen adsorption enthalpies by 10.5% and 3.4%, respectively. Due to the great variety of the MOFs chemical composition, larger spread in Henry's constants and nitrogen adsorption enthalpies was observed for them in comparison with zeolites. In this study, it was found that the account of electrostatic interactions results in the increase of the Henry's constants by at least a factor of 1.5 for 15% of all considered MOFs. Meanwhile, the analysis of Top 50 MOFs for VSA and PSA revealed a similar increase in Henry's constants for 80% and 58% of MOFs, respectively. So, the contribution of electrostatic interactions between quadrupole moments of nitrogen molecules and the MOF atoms is significant for the most promising MOFs. Since the chemical composition of both Top 50 MOFs includes a large number of different metal atoms with different atomic weights and partial charges, it can be assumed that metal atoms affect the adsorption performance primarily through the structure topology rather than through their contribution to the intermolecular interaction. Therefore, MOF screening for adsorption-based separation of N2/He should be based primarily on the search for MOFs with optimal structural parameters. Thus, an ideal adsorbent should have uniformly narrow pores without pockets or cavities and a significant density in order to provide high nitrogen adsorption enthalpy and N2/He selectivity. As a result, such an ideal structure will have low porosity, pore volume and accessible surface area. It is interesting to compare the results of the MOF screening obtained in this study with earlier literature data. A total of 213 MOFs were studied for VSA gas separation with very similar operational conditions (equimolar mixture of helium and nitrogen, gas pressure during the adsorption cycle 1.2 bar, pressure during the adsorbent regeneration 0.1 bar) [44]. UVEXAV was found to be the most promising MOF with API = 680 and adsorption selectivity S mix ads,N 2 /He = 222.7. In our study we discovered 24 MOFs with superior API metric and 34 MOFs with higher nitrogen selectivity. Additionally, in the same earlier study [44] it was demonstrated that the electrostatic contribution for most MOFs was negligible. Despite the fact that a nitrogen molecule has a relatively a low quadrupole moment, we believe that it is important to consider the electrostatic interactions between nitrogen molecules and the framework atoms during the screening. To estimate their contributions, Henry's constants and nitrogen enthalpies of adsorption were calculated with the account of electrostatic interactions and without them ( Figure S10). The account of electrostatic interactions results in the growth of the median value of Henry's constants and nitrogen adsorption enthalpies by 10.5% and 3.4%, respectively. Due to the great variety of the MOFs chemical composition, larger spread in Henry's constants and nitrogen adsorption enthalpies was observed for them in comparison with zeolites. In this study, it was found that the account of electrostatic interactions results in the increase of the Henry's constants by at least a factor of 1.5 for 15% of all considered MOFs. Meanwhile, the analysis of Top 50 MOFs for VSA and PSA revealed a similar increase in Henry's constants for 80% and 58% of MOFs, respectively. So, the contribution of electrostatic interactions between quadrupole moments of nitrogen molecules and the MOF atoms is significant for the most promising MOFs.
Since the chemical composition of both Top 50 MOFs includes a large number of different metal atoms with different atomic weights and partial charges, it can be assumed that metal atoms affect the adsorption performance primarily through the structure topology rather than through their contribution to the intermolecular interaction. Therefore, MOF screening for adsorption-based separation of N 2 /He should be based primarily on the search for MOFs with optimal structural parameters. Thus, an ideal adsorbent should have uniformly narrow pores without pockets or cavities and a significant density in order to provide high nitrogen adsorption enthalpy and N 2 /He selectivity. As a result, such an ideal structure will have low porosity, pore volume and accessible surface area.

Membrane-Based Gas Separation
To test the force fields used in this study for prediction of diffusion properties of MOFs and zeolites, calculated He and N 2 permeances in the temperature range of 298-301 K were compared with the corresponding literature data [22,26,29,30,[66][67][68][69][70][71]. Figure S2 demonstrates that the models used in this study can be applied only for estimation of the membrane permeance. In addition to the lack of ideality in the used force field models, there are several additional reasons leading to deviation of the calculated permeance values from the experimental ones. First, most polycrystalline membranes based on MOFs and zeolites have low selectivity and high permeance due to the presence of defects both inside their framework and in the packing of crystals. The majority of membranes demonstrate selectivity close to that of the Knudsen diffusion model, which indicates that the size of pores between the crystallites can exceed 2 nm. Second, only permeability can be calculated by simulations. To connect this value with permeance, it is necessary to know the membrane thickness. Except for single-crystal membranes [27], the thickness of the MOF or zeolite layer on the support can be substantially varied resulting in deviation between the calculated and experimental values. Third, an ideal model imitating the adsorbate behavior at infinite dilution was used to calculate permeability in this study. Such a model does not take into account the effects of competitive adsorption and gas mixture composition on the diffusion properties. As helium is weakly adsorbed, and its concentration determined by GCMC is substantially lower than the N 2 concentration, it is necessary to consider much larger framework fragments than for simulation of adsorption to obtain correct values of self-diffusion coefficients. In turn, this leads to a significant increase in the simulation time. Therefore, it is rational to use the ideal model at the first screening stage with the following correction of the membrane permeability and selectivity for Top 10 frameworks using Equations (10) and (11).
At the first stage of screening for the most promising materials for membrane-based separation of He and N 2 , diffusion coefficients were calculated at infinite dilution conditions. The analysis of the dependence of self-diffusion coefficients on the structural properties of the studied MOFs ( Figure S12) demonstrates a natural trend towards an increase in self-diffusion coefficients with increasing pore size (PLD and LCD) and porosity. The greatest difference between D 0 N 2 and D 0 He was observed in the PLD range under 6 Å and LCD range under 8 Å. The diffusion selectivity S 0 di f He/N 2 in such small pores can reach several orders of magnitude ( Figure S13a,b). This result can be explained by the difference in the kinetic diameters of helium and nitrogen molecules. Similar relationships between structural parameters and diffusion selectivity were also observed for zeolites ( Figure S13). Similar to the adsorption-based separation, the screening demonstrated that great variety of MOFs makes it possible to find MOFs with much higher diffusion selectivity than that of zeolites (maximum S 0 di f He/N 2 for MOFs is equal to 1197 compared to 29 for zeolites).
The analysis of the screening results showed that the membrane selectivity S 0  In addition, the effect of electrostatic interaction on self-diffusion coefficients and nitrogen permeability as well as membrane selectivity to helium was studied ( Figure S10).  In addition, the effect of electrostatic interaction on self-diffusion coefficients and nitrogen permeability as well as membrane selectivity to helium was studied ( Figure S10). The account of electrostatic interaction results in a minor decrease in D 0 N 2 (charged) relative to D 0 N 2 (non−charged) . As electrostatic interaction has a more significant effect on Henry's adsorption constants compared to self-diffusion coefficients, on the average the nitrogen permeability grows after the account of electrostatic interaction. As a result, the membrane selectivity to helium decreases. So, electrostatic interactions between the quadrupole moment of nitrogen molecules and the framework atoms should be taken into account during both adsorption-based and membrane-based screening.

The account of electrostatic interaction results in a minor decrease in
The perspectives of using MOFs and zeolites as membrane materials for separation of He and N 2 can be estimated by comparing calculated selectivity and permeability with the upper bonds obtained in different years for polymeric membranes. Upper bond is an empirical relation P i = kS mem i/j n where constants k and n are determined from the selectivity vs. permeability graph for experimental data on He/N 2 separation over various polymeric membranes. Figure 6 presents upper bonds obtained by Robeson in 1991 [72] and 2008 [12] and by Wu et al. in 2019 [73]. The upper bonds for polymeric membranes demonstrate that the best MOFs and zeolites trail polymeric membranes in selectivity S 0 mem He/N 2 but have higher permeability P 0 He . Based on this fact, several variants of using MOFs and zeolites for separation of He and N 2 are possible: creation of materials from a combination of polymers and MOFs/zeolites or using MOF/zeolite materials selective to nitrogen. In the former case, combination of materials can result in the growth in permeability (in comparison with polymeric membranes) and increase in helium selectivity (in comparison with MOFs/zeolites). In several publications [49][50][51] it was shown that such an approach makes it possible to obtain more efficient materials for separation of different gas mixtures than using MOFs along as the membrane materials. However, due to the low membrane selectivity of MOFs and zeolites to helium, apparently, any other mesoporous materials with the membrane selectivity close to the Knudsen diffusion selectivity can be used for separation of helium and nitrogen. In the second case, if MOFs or zeolites selective to nitrogen are used as the membrane, the retentate will be enriched with helium and the permeate will be enriched with nitrogen. As using MOFs and zeolites as membrane materials is more promising in the second case, Top 50 frameworks most selective to nitrogen (S 0 mem N 2 /He > 22.6) were identified. Not a single zeolite was included in this list. selectivity of MOFs and zeolites to helium, apparently, any other mesoporous materials with the membrane selectivity close to the Knudsen diffusion selectivity can be used for separation of helium and nitrogen. In the second case, if MOFs or zeolites selective to nitrogen are used as the membrane, the retentate will be enriched with helium and the permeate will be enriched with nitrogen. As using MOFs and zeolites as membrane materials is more promising in the second case, Top 50 frameworks most selective to nitrogen Structure-membrane separation performance relationships were also studied. Their analysis was performed using the same technique that was used for adsorption-based separation and was described above in Section 3.1. The following optimal structural parameters were determined for membrane-based separation ( Figure S14): 3.75 Å < PLD < 4.65 Å, 4 Å < LCD < 5.8 Å, 100 m 2 /g < ASA < 450 m 2 /g, 0.01 cm 3 /g < AV < 0.06 cm 3 /g, 1100 kg/m 3 < density < 2400 kg/m 3 , 0.02 < VF < 0.1. Note that the PDF functions for Top 50 MOFs and all the studied MOFs are substantially different, and the range of optimal structural parameters is much narrower than for adsorption-based separation. Similarly, chemical composition-membrane separation performance relationships were analyzed. Top 50 MOFs include 22 metals. It means that the membrane properties are primarily determined by its structural parameters rather than by chemical composition. Still, based on the ratio of probabilities of the presence among Top 50 MOFs and in all the studied MOFs, Ga and Al were found to be the most promising metals for membrane-based separation.
The characteristics of the most promising MOFs (Top 10) for membrane-based separation are reported in Table S8. All MOFs from Top 10 demonstrate high adsorption selectivity predominating over their diffusion selectivity and high N2 permeability exceeding the He permeability by more than an order of magnitude. At the second stage of Structure-membrane separation performance relationships were also studied. Their analysis was performed using the same technique that was used for adsorption-based separation and was described above in Section 3.1. The following optimal structural parameters were determined for membrane-based separation ( Figure S14): 3.75 Å < PLD < 4.65 Å, 4 Å < LCD < 5.8 Å, 100 m 2 /g < ASA < 450 m 2 /g, 0.01 cm 3 /g < AV < 0.06 cm 3 /g, 1100 kg/m 3 < density < 2400 kg/m 3 , 0.02 < VF < 0.1. Note that the PDF functions for Top 50 MOFs and all the studied MOFs are substantially different, and the range of optimal structural parameters is much narrower than for adsorption-based separation. Similarly, chemical composition-membrane separation performance relationships were analyzed. Top 50 MOFs include 22 metals. It means that the membrane properties are primarily determined by its structural parameters rather than by chemical composition. Still, based on the ratio of probabilities of the presence among Top 50 MOFs and in all the studied MOFs, Ga and Al were found to be the most promising metals for membrane-based separation.
The characteristics of the most promising MOFs (Top 10) for membrane-based separation are reported in Table S8. All MOFs from Top 10 demonstrate high adsorption selectivity predominating over their diffusion selectivity and high N 2 permeability exceeding the He permeability by more than an order of magnitude. At the second stage of screening, these MOFs were used in simulation by GCMC and MD using an equimolar mixture of He and N 2 at 298 K and 3 bar, which corresponds to the conditions of real membrane-based separation of gases. The simulation results (Table S8) indicate that the switch to the equimolar mixture leads to simultaneous decrease of the nitrogen adsorption selectivity, He/N 2 diffusion selectivity and permeability. Nevertheless, the values of S 0 mem N 2 /He and S mix mem N 2 /He are very similar indicating that the suggested technique used for the initial screening makes it possible to predict promising frameworks while using minimum computer power. Earlier [44], based on the results of in silico screening of 500 MOFs, the highest membrane selectivity S 0 mem N 2 /He = 44.91 was found for LIFWOO. In this study, we determined 16 MOFs with even higher membrane selectivity. Thus, highly nitrogen-selective membranes are characterized by the predominance of the N 2 /He adsorption selectivity over N 2 /He diffusion selectivity. Therefore, as in the case of adsorption-based separation, an ideal membrane should have narrow uniform pores and a significant density. The membrane properties are primarily determined by its structural parameters rather than by chemical composition.

MOF and Zeolite Databases
"Computation Ready, Experimental Metal−Organic Framework Database" (CoRE MOF 2019) [34] was used as the parent MOF database. At the first step only ordered frameworks where all solvent molecules were removed were selected from this database. The resulting database consisted of 10143 MOFs. Then, the Zeo++ algorithm [74] was used to calculate density, porosity (VF), accessible surface area (ASA), pore limiting diameter (PLD) and largest cavity diameter (LCD) for each framework. Additionally, Zeo++ was used to identify MOFs with pockets not accessible both for helium and nitrogen molecules. In the following computations, such pockets were blocked. The accessible pore volume (AV) was determined using RASPA-2.0 software package [75] with helium molecule used as a probe. Further, only MOFs with non-zero ASA and PLD > 3.75 Å were selected from the database. These conditions would allow both gases to penetrate through the membrane (the kinetic diameter is equal 3.64 Å for N 2 and 2.56 Å for He). So, the number of analyzed MOFs was shortened to 5156.
IZA database [32] implemented in iRASPA visualization package [76] was used as the parent database of all-silica zeolites and zeolite-like materials. Structural parameters were determined for all frameworks from this database using a procedure similar to the one used earlier for MOFs. After removing frameworks with PLD below 3.75 Å, the number of analyzed zeolite frameworks was reduced to 110.

Simulation Details
The adsorption and diffusion properties of the studied frameworks were simulated using equilibrium molecular dynamics (EMD) and Grand canonical Monte Carlo (GCMC) methods in the RASPA-2.0 package. The interactions were described by the sum of Lennard-Jones 6-12 (LJ) and Coulomb interaction potentials. The nitrogen molecule was simulated according to the TraPPE [77] force field as a dumbbell with a rigid bond between the atoms with the length of 1.1 Å. The LJ interaction parameters for each atom were ε N /k B = 36.0 K and σ N = 3.31 Å. The N 2 quadrupole moment was described by three charges: two -q charges where q = 0.482e that were located at the centers of the nitrogen atoms and +2q charge located at the molecule center of masses. The helium molecule was simulated by a single-charge model [78] with parameters σ He = 2.64 Å and ε He /k B = 10.9 K. The constants of the LJ 6-12 potential for atoms in MOFs were simulated according to the Dreiding [60] force field. If the necessary parameters were missing, the required data were retrieved from the UFF [79] force field. The force field developed by Vujić and Lyubartsev [62] was used for simulation of adsorption and diffusion in zeolites. Cross constants of the LJ interaction were calculated using the Lorentz-Berthelot mixing rule. The LJ interaction was cut at the distance R cutoff = 12.8 Å. Its potential was shifted to zero starting from the distance 0.9 R cutoff . Determination of partial atomic charges by REPEAT [80] or DDEC [81] methods that have high precision requires periodic DFT calculations for each studied structure. A detailed and systematic analysis of the effect of the choice of framework partial atomic charges on CO 2 adsorption in six different widely studied MOFs predicted by molecular simulations was performed in the study [82]. It was shown that the partial charges obtained by both DDEC and REPEAT methods yielded isotherms that were in good agreement with each other. In our previous studies, it was shown that partial charges obtained by DDEC and REPEAT methods lead to similar results for modeling hydrogen adsorption in SAPO-11 (zeolite-like material) at 77 K [83] and water adsorption in CAU-10-H (MOF) at 298 K [84]. Thus, it can be expected that the use of the DDEC or REPEAT methods for calculating atomic partial charges will lead to similar results in the case of modeling nitrogen adsorption in MOFs and zeolites. A reasonable alternative is to use pre-trained models obtained by machine learning based on the CoRE MOF DDEC [43,85] containing 2932 optimized structures with partial atomic charges calculated by the DDEC method. So, partial atomic charges of the framework atoms were determined using pre-trained Random Forest model in the PACMOF code [86]. Long-range Coulomb interactions were calculated using the Ewald summation technique. The size of the MOF and zeolite structure fragments was selected to ensure the minimum distance in each direction exceeded 2R cutoff . In each simulation, we assumed that the MOF or zeolite structure was rigid and did not contain any defects.
Henry's constants were calculated using the Widom particle insertion method. Henry's constants were calculated at infinite dilution conditions for 10 5 cycles at 298 K. The adsorption selectivity at infinite dilution S 0 ads,i/j was determined as the ratio of Henry's constants K 0 i [10]: The adsorption selectivity for separation of gas mixtures S mix ads,i/j was calculated for 0.01, 0.1, 0.3 and 1 MPa pressures at 298 K as follows [10]: Here, N i is the adsorption of the i-th component determined from GCMC, x i is the molar fraction of the i-th component. For a gas mixture, the adsorption values N i were simulated using the following GCMC moves: translation, rotation, insertion, deletion, reinsertion, identity exchange of He and N 2 molecules. The fugacity coefficient was calculated from the Peng-Robinson equation of state.
Adsorbent Performance Score (APS) used to estimate the adsorbent efficiency was calculated as follows [87]: The sorbent regenerability was calculated as [10]: Self-diffusion coefficient was determined by equilibrium molecular dynamics from the root-mean-square particle displacement using the following formula [88]: Here, N is the number of molecules, → r j (t) is the position of the j-th particle at the moment t, d is the dimension of the system. Self-diffusion coefficient at infinite dilution D 0 i was simulated by positioning 30 adsorbate molecules in the MOF fragment with all interactions between the gas molecules switched off. The simulation was performed in the NVT ensemble (at constant number of particles, volume and temperature) using the Nosé-Hoover thermostat. These conditions simulate the properties of the adsorbate at infinite dilution. After the initial position of the adsorbate in the framework was generated, for the first 50 ps the system was subjected to equilibration before the data collection. The MD simulation time was 20 ns. The diffusion selectivity at infinite dilution S 0 di f ,i/j was calculated as follows [42]: Permeability P 0 i was determined using Equation (7) [42]: Membrane selectivity S 0 mem,i/j was estimated as follows [42]: For selected Top 10 frameworks, self-diffusion coefficient D self,i was calculated for each component of the gas mixture at the adsorbate concentrations determined by GCMC. As helium is weakly adsorbed, and its concentration determined by GCMC is significantly lower than that of N 2 , the framework fragment used in the EMD method was 27 times larger than during the initial screening at infinite dilution. Five independent EMD runs were performed to collect statistics. Diffusion selectivity for a mixture of gases S mix di f ,i/j was determined as the ratio of self-diffusion coefficients [42]: Permeability of the i-th mixture component P mix i was calculated according to Equation (10) [42]: Here, ϕ is the adsorbent porosity, c i is the concentration of the i-th mixture component at the feed side of the membrane determined by GCMC, f i is the fugacity of the i-th mixture component before the membrane. This model assumes vacuum permeating pressure. The membrane selectivity for the mixture of gases S mix mem,i/j was determined as follows [42]:

Conclusions
In silico screening of 10143 MOFs and 218 all-silica zeolites for adsorption-based and membrane-based separation of He and N 2 was performed. GCMC and EMD methods were used to calculate Henry's constants, adsorption at different pressures and self-diffusion coefficients for He and N 2 . These values were used to calculate major metrics, such as adsorption, diffusion and membrane selectivity, regenerability and permeability. Dependence of adsorption selectivity S mix ads,N 2 /He in an equimolar mixture of He and N 2 on the pressure in the system was studied. The effect of electrostatic interaction between the quadrupole moment of nitrogen molecules and framework atoms on the main adsorption and diffusion metrics was studied. MOFs were demonstrated to be more efficient than zeolites in both adsorption-based and membrane-based separation. Top 10 most promising MOFs for membrane-based, VSA and PSA separation methods were determined. The analysis of structure-adsorption and membrane performance relationships made it possible to determine the ranges of structural parameters, such as pore limiting diameter, largest cavity diameter, surface area, porosity, accessible surface area and pore volume, corresponding to the most promising MOFs for each separation model considered in this study. A similar analysis was performed to find out the optimal MOF chemical composition. The screening performed in this study can be called topological because the structural flexibility, possible presence of structural defects or modification of linkers with functional groups were not taken into account. Similarly, in the case of zeolites possible changes of the Si/Al ratio, variation of the cation composition or modification of their structure by isomorphous substitutions were not considered. Despite these limitations, one can expect that the ranges of optimal structural parameters and trends of adsorption and membranes metrics determined in this study will be correct even if all these factors are taken into account. The obtained results can be used as a guide for selection of frameworks for He/N 2 separation.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/molecules28010020/s1, Section S1: Comparison of molecular simulations and experimental data; Section S2: Structure-adsorption performance relationships; Section S3: Comparison of ideal and mixture adsorption selectivities; Section S4: Relationship between API and APS metrics; Section S5: Comparison of the structural properties PDF for the best 50 MOF structures for VSA and PSA with all MOF structures; Section S6: Characteristics of Top 10 MOFs for adsorption-based separation; Section S7: Effect of electrostatic interaction on He/N 2 separation; Section S8: Structurediffusion performance relationships; Section S9: Comparison of the structural properties PDF for the best 50 MOF structures for membrane separation process with all MOF structures; Section S10: Chemical composition-membrane separation performance relationships; Section S11: Characteristics of Top 10 MOFs for membrane-based separation; References.