Quantitative Structure–Property Relationships from Experiments for CH4 Storage and Delivery by Metal–Organic Frameworks

Quantitative structure–property relationships (QSPRs) can be applied to metal–organic frameworks (MOFs) to allow for reasonable estimates to be made of the CH4 storage performance. QSPRs are available for CH4 storage of MOFs, but these were obtained from Grand Canonical Monte Carlo (GCMC) simulations which have come under scrutiny and of which the accuracy has been questioned. Here, QSPRs were developed from experimental data and insights are provided on how to improve storage and deliverable CH4 storage capacity based on material properties. Physical properties of MOFs, such as density, pore volume, and largest cavity diameter (LCD), and their significance for CH4 storage capacity were assessed. One relationship that was found is that CH4 gravimetric storage capacity is directly proportional to Brunauer–Emmett–Teller (BET) surface area (r2 > 90%). The QSPRs demonstrated the effect of van der Waals forces involved in CH4 adsorption. An assessment was made of the accuracy of QSPRs made by GCMC as compared to QSPRs derived from experimental data. Guidelines are provided for optimal design of MOFs, including density and pore volume. With the recent achievement of the gravimetric 2012 DOE CH4 storage target, the QSPRs presented here may allow for the prediction of structural descriptors for CH4 storage capacity and delivery.


Introduction
Metal-organic frameworks (MOFs) are compounds comprised of metal nodes coordinated to organic ligands [1][2][3][4][5][6]. The metal ions form a cluster known as the secondary building unit (SBU), which can be linked by organic compounds through reticular synthesis to form designed three-dimensional structures for CH 4 storage and delivery ( Figure 1). Reticular synthesis allows for the molecular assembly of MOFs with a structure tailored for CH 4 storage. HKUST-1 (Hong Kong University of Science and Technology) is among the MOF compounds that can be used for methane storage [7]. HKUST-1 is a MOF used as a benchmark for medium-pressure ranges (35-65 bar) [8]. HKUST-1 is composed of copper centers and benzene-1,3,5-tricarboxylate (BTC) organic linkers [9,10]. This MOF is notable and has a surface area of 1800 m 2 g −1 [9,11,12]. Adapted with permission from Reference [9]. Copyright 2013 American Chemical Society.
The design of MOFs for high CH 4 storage and deliverable capacity requires an understanding of how MOFs' structure and properties affect storage and deliverable capacity. This is a highly effective strategy because MOFs can be designed using a molecular building-block approach to possess certain geometrical and topological features, which lead to the measured density, pore volume, void fraction, etc. of these materials. The relationships between these properties and CH 4 storage and deliverable capacity developed through quantitative structure-property relationships (QSPRs) allow these materials to be designed, synthesized, and tested to meet CH 4 storage targets more efficiently. One storage target was set by the 2012 United States (US) Department of Energy (DOE) program for methane-storage systems, called Methane Opportunities for Vehicular Energy (MOVE) [13]. The target set by the program was to achieve a gravimetric capacity of 0.5 g (methane) g −1 and a volumetric capacity of 350 cm 3 (STP = standard temperature and pressure equivalent: T = 273.15 K, P = 101.325 kPa) cm −3 , which corresponds to 250 g (methane) L −1 [14][15][16]. The critical issue for this approach is how to develop structure-property relations which are effective and accurate based on key structural features, in order to be able to make a quick estimate of whether the MOF is likely to meet the storage target.
The development of QSPRs provides models that may be used to yield estimates of CH 4 storage performance by provided key material properties, without the theory required for first-principles calculations. QSPRs developed based on three-dimensional descriptors have been wrong [17]. QSPR models have been developed using Grand Canonical Monte Carlo simulations [18]. For computational methods, there is always a question regarding the necessary degree of detail associated with the molecular structure of the MOF. Different computational approaches are available, including molecular mechanics and ab initio programs. Structures can be complex and heterogeneous, Crystals 2020, 10, 700 3 of 15 and various 3D molecular descriptors may therefore be used. In addition, material structures are not static but are in fact dynamic, and have various molecular conformations. This may lead to misleading information from QSPRs developed using simulations. This is a fundamental problem that needs to be addressed for QSPRs to become more effective and reliable.
In this study, QSPR analysis was used to predict CH 4 adsorption capacity and, more importantly, delivery, based on experimental properties and the underlying molecular structure of MOFs. For practical applications, deliverable capacity is an important variable to assess the performance of these materials. The deliverable capacity is a measure of the amount of available gas released between the upper storage pressure and lower pressure (5 bar) required at the engine inlet of the adsorption-desorption process. The deliverable capacity is affected by properties of the MOF material such as pore architecture [19], pore spacing [20], linker [21,22], functional groups [23,24], and hydrophobicity [25,26]. For this QSPR study, experimental descriptors were selected and QSPR analysis was conducted via regression analysis. BET surface area, pore volume, and density were experimental descriptors used in the analysis, as well as the largest cavity diameter (LCD) and pore limiting diameter (PLD) which were used to molecularly represent the pores. The accuracy of computational QSPR models was assessed as compared to the developed experimental model. In addition, a relationship between experimental properties of MOFs and CH 4 storage and delivery was developed and investigated.

Materials and Methods
Over 150 MOFs with known experimentally determined crystal structures, densities, pore volumes, isosteric heats of adsorption (Q st ), CH 4 uptakes, and CH 4 deliveries were selected for QSPR analysis (Table 1 and Supplementary Materials). In addition, the computation-ready experimental (CoRE) MOF database was used to supplement the experimental data to obtain additional physical properties of the MOFs [27]. The CoRE database does not contain any disordered structures. Because the QSPR models were determined from experimental data, the reliability of the work was not affected by the disorder, group position disorder, or lattice deformation of the MOF, or the properties that are not included in the CoRE database. These factors were already taken into account during measurement of the properties obtained through N 2 and CH 4 adsorption. In the future, artificial intelligence (AI) and machine-learning algorithms may be able to play a role in incorporating and modeling these features in the structures of these real materials to maximize CH 4 storage [28]. The CH 4 uptake includes excess CH 4 adsorbed. The framework topologies, LCD, PLD, accessible surface area, metal sites, and whether or not the MOFs had open metal sites were noted. LCD was used as an important descriptor to determine whether the pores were efficiently utilized for adsorption and packing of CH 4 molecules without dead or free space. LCD was selected based on its previously demonstrated success in improving CH 4 storage [15]. The objective was to increase the adsorption over that of a pressurized CH 4 tank that does not contain a MOF. Organic linkers and functional groups were also selected to better analyze the relationships between structure and performance. This information allowed the MOFs to be classified and grouped. This classification was used to partition MOFs into subsets based on their common properties to obtain additional QSPRs.  [29].  a Vp: pore volume. b Total volumetric uptake, cm 3 (STP) cm −3 . c The deliverable amount is defined as the difference in total uptake between 5 bar and the specified upper limiting working pressure under isothermal conditions. Q st is the isosteric heat of adsorption of Data were obtained at 298 K and at 65 and 35 bar. Deliverable capacities were determined based on the difference between the available gas released between the upper storage pressure and lower pressure (5 bar) required at the engine inlet of the adsorption-desorption process. Gravimetric uptakes were calculated when not available, based on the density of the MOF and the density of CH 4 gas at STP. LCD, PLD, and metal sites were obtained from the CoRE MOF database [27]. Adsorption isotherms were used when available and needed to determine deliverable capacity when not reported.

MOF
The electronic environment of the MOF was also considered as a descriptor, including the electronegativity of functional groups in MOFs and their effects on van der Waals interactions between the MOF and CH 4 . Correlational analysis was used to establish regressions. Various structural descriptors were used to describe the molecular dimensions of the MOFs, as well as intensive experimental variables as opposed to extensive variables.
The results were plotted using a 5D visualization platform named the CH 4 Adsorption MOF Data Explorer, which was generated using experimental data to establish structure-property relationships. This allowed for the generation of structure-property figures at various pressures and conditions. The explorer allowed for the assessment of uptake and delivery of CH 4 . In the explorer, MOFs that were considered outliers could be monitored. Various structural parameters were assessed, including the pore volume, Vp, BET surface area, void fraction, Q st , LCD, LPD, etc., and snapshots were captured of their behavior for QSPR analysis.
Molecular and experimental descriptors of the MOFs were analyzed and selected for QSPR model building. The model was validated by calculating the coefficient of determination or r 2 . r 2 is the proportion of the variance in the dependent variable that is predictable from the independent variable. r 2 was determined using the following equation: where SSR is the sum of squares of regression and SST is the total sum of squares. r 2 measures how effective a model is at predicting CH 4 uptake or delivery.

Results
Important variables influencing CH 4 adsorption capacity and CH 4 delivery were identified through QSPR analysis based on experiments. These structural descriptors were used to predict CH 4 uptake and delivery in order to provide insights into MOF design and to improve material design and discovery. Important structural parameters which have a significant effect on CH 4 storage and delivery were identified. The density of the MOF and the LCD were found to be significant.
The relationship between the properties of MOFs is important in the design of better CH 4 adsorbents. An increase in pore volume, V p , was associated with an increase in BET surface area ( Figure 2). Increasing BET surface area was associated with decreasing density and increasing LCD. For a BET surface area greater than 6000 m 2 /g, the pore volume was estimated to be greater than 3 cm 3 /g. Larger BET surface areas were associated with lower density and larger LCD. The relationship between pore volume and BET surface for MOFs from experiments was provided by the following equation.
BET sur f ace area For this equation, r 2 = 0.881. The fact that there is an estimated relationship between pore volume and BET may support the selection one structural descriptor to estimate CH4 uptake and delivery to simplify QSPR analysis. Because of the minute details and the wide design freedom accessible for MOFs, not every structure is expected to exactly fit this model, but most structures experimentally investigated to date do, as indicated by the r 2 of nearly 0.90. CH4 volumetric and gravimetric uptakes were compared at 35 bar and 298 K as functions of density and LCD ( Figure 3). Increasing volumetric uptake was associated with decreasing gravimetric uptake. The results showed that as the density of the MOF increased from 0.3 g/cm 3 to greater than 1 g/cm 3 , the gravimetric uptake decreased from around 0.24 g/g to around 0.125 g/g. As the mass of MOF per unit volume increased, the mass of CH4 stored decreased. As LCD decreased, gravimetric uptake also decreased. As LCD decreased, volumetric uptake increased. Principally, CH4 For this equation, r 2 = 0.881. The fact that there is an estimated relationship between pore volume and BET may support the selection one structural descriptor to estimate CH 4 uptake and delivery to simplify QSPR analysis. Because of the minute details and the wide design freedom accessible for MOFs, not every structure is expected to exactly fit this model, but most structures experimentally investigated to date do, as indicated by the r 2 of nearly 0.90. CH 4 volumetric and gravimetric uptakes were compared at 35 bar and 298 K as functions of density and LCD (Figure 3). Increasing volumetric uptake was associated with decreasing gravimetric uptake. The results showed that as the density of the MOF increased from 0.3 g/cm 3 to greater than 1 g/cm 3 , the gravimetric uptake decreased from around 0.24 g/g to around 0.125 g/g. As the mass of MOF per unit volume increased, the mass of CH 4 stored decreased. As LCD decreased, gravimetric uptake also decreased. As LCD decreased, volumetric uptake increased. Principally, CH 4 adsorption takes place in the pores. The diameters of these pores are four or five times the molecular diameter of CH 4 . Fundamentally, the interaction between CH 4 and the pore plays a significant role in adsorption. With increasing pore size, the force associated with the interaction between CH 4 and the pore decreases rapidly, such that the equilibrium adsorption becomes like that of a plane surface. The consequence of this is that any pores with widths greater than 2 nm, which excludes mesopores, do not significantly improve of CH 4 storage as compared to a system with a compressed gas. To design high-performing adsorbents for CH 4 storage, the volume of the micropores per unit volume of adsorbent should be maximized to improve CH 4 storage as compared to a compressed gas system. At the same time, void volume in the system packed with adsorbent should be minimized. As shown in Figure 4, gravimetric uptake at 35 bar increased with gravimetric uptake at 65 bar. Increasing gravimetric uptake was associated with decreasing density and increasing LCD (Figure 4). For the relationship shown in Figure 4: Gravimetric uptake at 65 bar g g = 1.6997 × Gravimetric uptake at 35 bar g g − 0.0801 where r 2 = 0.9331.
Crystals 2020, 10, 700 7 of 15 gravimetric uptake was associated with decreasing density and increasing LCD (Figure 4). For the relationship shown in Figure 4: where r 2 = 0.9331.   CH4 gravimetric uptake as a function of BET surface area was assessed at 35 bar and 298 K for different MOF densities and LCDs ( Figure 5). Increasing the BET surface area increased gravimetric uptake. As the BET surface area increased, there was usually an increase in the LCD and a decrease in the density. A QSPR model was built, yielding two different equations for two different pressures at 298 K. The model was validated by calculating the coefficient of determination or r 2 .
where r 2 = 0.7973, indicating that this equation can be used to predict CH4 uptake. At 65 bar, where r 2 = 0.7973, indicating that this equation can be used to predict CH 4 uptake. At 65 bar, Gravimetric Uptake at 65 bar g g = 5 × 10 −5 × BET Sur f ace Area m 2 g + 0.0916 (5) where r 2 = 0.9086m indicating that this equation can be used to predict CH 4 uptake. The predictive power of the QSPR model proposed in Equation (5) is shown in Figure 6 for a test set of MOFs. The model was found to have excellent predictive power over the range of data tested.  In practice, gravimetric deliverable capacity is a more important parameter than gravimetric uptake for energy storage applications. Structural features of MOFs should be designed for low CH4 uptake at low pressure (5 bar). As shown in Figure 7, two main structural features which showed an apparent trend were the experimental density and LCD of the MOF. Increasing BET surface area increased gravimetric delivery at 35 bar. This was associated with a decrease in density. Density and LCD are intricately linked to void space. Another material property to consider is the Qst of CH4. Decreasing the Qst from around 22 kJ/mol to around 15 kJ/mol was associated with increasing gravimetric delivery from 0.06 g/g to around 0.  In practice, gravimetric deliverable capacity is a more important parameter than gravimetric uptake for energy storage applications. Structural features of MOFs should be designed for low CH4 uptake at low pressure (5 bar). As shown in Figure 7, two main structural features which showed an apparent trend were the experimental density and LCD of the MOF. Increasing BET surface area increased gravimetric delivery at 35 bar. This was associated with a decrease in density. Density and LCD are intricately linked to void space. Another material property to consider is the Qst of CH4. Decreasing the Qst from around 22 kJ/mol to around 15 kJ/mol was associated with increasing gravimetric delivery from 0.06 g/g to around 0.2 g/g, based on experimental results at 35 bar (SI). At 65 bar, gravimetric delivery increased when decreasing Qst from 22 kJ/mol to around 15 kJ/mol, and decreased when further decreased to around 10 kJ/mol. These structural features must be optimized A B To meet the DOE's target of 0.5 g/g at 65 bar and 298 K, it was estimated that the BET surface area of the MOF should be greater than 7000 m 2 /g, based on the QSPR analysis presented here. At 100 bar and 270 K, a CH 4 uptake of 0.66 g/g was met by NU-1501-Al, which has a BET area of 7310 m 2 /g [40]. As indicated by Figure 3, it is difficult to simultaneously achieve both the gravimetric and volumetric DOE CH 4 storage targets due to the observation that with increasing gravimetric uptake, volumetric uptake decreases.
In practice, gravimetric deliverable capacity is a more important parameter than gravimetric uptake for energy storage applications. Structural features of MOFs should be designed for low CH 4 uptake at low pressure (5 bar). As shown in Figure 7, two main structural features which showed an apparent trend were the experimental density and LCD of the MOF. Increasing BET surface area increased gravimetric delivery at 35 bar. This was associated with a decrease in density. Density and LCD are intricately linked to void space. Another material property to consider is the Q st of CH 4 . Decreasing the Q st from around 22 kJ/mol to around 15 kJ/mol was associated with increasing gravimetric delivery from 0.06 g/g to around 0.2 g/g, based on experimental results at 35 bar (SI). At 65 bar, gravimetric delivery increased when decreasing Q st from 22 kJ/mol to around 15 kJ/mol, and decreased when further decreased to around 10 kJ/mol. These structural features must be optimized to maximize gravimetric delivery. From a holistic point of view where both the gravimetric and volumetric delivery targets are met, a particular strategy may work to achieve this, as shown in Figure 8. This strategy involves the use of medium-sized LCDs and densities of around 0.8 g/cm 3 . For this design challenge, using structures with large LCDs of 25 Å or more is not sufficient anymore. From a holistic point of view where both the gravimetric and volumetric delivery targets are met, a particular strategy may work to achieve this, as shown in Figure 8. This strategy involves the use of medium-sized LCDs and densities of around 0.8 g/cm 3 . For this design challenge, using structures with large LCDs of 25 Å or more is not sufficient anymore. density, and the diameter of each circle represents the largest cavity diameter (LCD) of each MOF.
From a holistic point of view where both the gravimetric and volumetric delivery targets are met, a particular strategy may work to achieve this, as shown in Figure 8. This strategy involves the use of medium-sized LCDs and densities of around 0.8 g/cm 3 . For this design challenge, using structures with large LCDs of 25 Å or more is not sufficient anymore.  In the case of volumetric CH 4 delivery at 35 bar and 298 K, volumetric delivery decreased with increasing BET surface area. This was accompanied by a decrease in density and increase in LCD (Supplementary Figure S1). For a BET surface area of around 6000 m 2 /g, the volumetric delivery was around 75 cm 3 (STP)/cm 3 . Volumetric uptake at 65 bar increased linearly with volumetric uptake at 35 bar (Supplementary Figure S2). Comparing volumetric delivery at 65 bar to that at 35 bar, it was noted that there was an optimum at a density of around 1 g/cm 3 (Figure 9). Therefore, this indicated that the optimum structural properties for CH 4 delivery were nearly the same at various pressures. For volumetric delivery at 35 bar, the optimum experimental result was obtained at lower BET surface areas and smaller LCDs. This coincided with Q st values of around 17.5 kJ/mol to 19 kJ/mol. In the case of volumetric CH4 delivery at 35 bar and 298 K, volumetric delivery decreased with increasing BET surface area. This was accompanied by a decrease in density and increase in LCD (Supplementary Figure S1). For a BET surface area of around 6000 m 2 /g, the volumetric delivery was around 75 cm 3 (STP)/cm 3 . Volumetric uptake at 65 bar increased linearly with volumetric uptake at 35 bar (Supplementary Figure S2). Comparing volumetric delivery at 65 bar to that at 35 bar, it was noted that there was an optimum at a density of around 1 g/cm 3 (Figure 9). Therefore, this indicated that the optimum structural properties for CH4 delivery were nearly the same at various pressures. For volumetric delivery at 35 bar, the optimum experimental result was obtained at lower BET surface areas and smaller LCDs. This coincided with Qst values of around 17.5 kJ/mol to 19 kJ/mol.

Discussion
The results shown here provide several QSPR models that can be used to design MOFs for CH4 storage at 298 K. Observations were made on how to optimize gravimetric CH4 delivery and volumetric CH4 delivery through correlations. Some trends were also observed regarding the structures of these materials. One of the key questions asked during this work was whether

Discussion
The results shown here provide several QSPR models that can be used to design MOFs for CH 4 storage at 298 K. Observations were made on how to optimize gravimetric CH 4 delivery and volumetric CH 4 delivery through correlations. Some trends were also observed regarding the structures of these materials. One of the key questions asked during this work was whether computational QSPR analysis was accurate. Furthermore, can advanced material design (mesopores, low temperature) be incorporated into QSPR analysis? Computational results are good and provide insights, but sometimes they miss some important or minute details from the real systems.
Regarding previous studies, QSPRs were established using 130,000 hypothetical MOF structures [41]. However, because the structures were hypothetical and the results were computationally obtained, this work sought to verify and complete QSPR analysis using experimental results and synthesized MOF structures. The results from the QSPR models can aid future MOF development and be used to double-check how accurate and realistic GCMC simulation predictions are. Critical to both endeavors is the identification of key structural parameters that can be controlled through reticular synthesis, which significantly affect CH 4 uptake and delivery. Computationally, Fernandez et al. [41] focused on the following structural properties in their QSPR analysis: dominant pore diameter, surface area, and void fraction. It was concluded that for CH 4 storage at 35 bar, the density of the MOF should be greater than 0.43 g/cm 3 and the void fraction should be greater than 0.52. One drawback of this statement is that the simulated MOFs were hypothetical and no upper bound is given to narrow the range. In another study, 122,835 hypothetical MOFs were simulated for methane storage [32]. At 65 bar, the optimal structural parameters for these hypothetical MOFs was a Q st between 10.5 and 13.0 kJ/mol, an LCD between 10 and 12 Å, and a volumetric surface area between 2100 and 2300 m 2 /cm 3 . In another study, 45 MOFs were simulated for CH 4 storage at 65 bar for an investigation on structural parameters, of which void fraction and Q st were found to be most significant [18]. The optimum Q st was 30 kJ/mol with a void fraction of 0.9 at 35 bar. More recently, it was demonstrated that these hypothetical MOF structures can be synthesized [42]. However, most hypothetical MOFs have not actually been tested for CH 4 adsorption.
One certainty from simulated CH 4 adsorption on MOFs is that there does not seem to be a clear consensus across all reports on what the best structural parameters are. In this work, more than 150 MOFs were analyzed and their CH 4 adsorption results from experiments were used to create a previously unreported QSPR model. This gave realistic results that can be referenced and will allow better design of MOFs for CH 4 adsorption. QSPR correlations were created for simplicity and were modeled as linear equations which possessed r 2 values of around 0.9. From the analysis of CH 4 data from 150 MOFs, it was observed that estimates of an optimum Q st of more than 30 kJ/mol are too high at 65 bar and 298 K. As far as density is concerned, the value depends on whether high volumetric or gravimetric delivery or uptake is being targeted. Starkly different optima were obtained in these cases. The recommendation given here is to design MOFs with an optimum density that is determined based on whether the storage target is volumetric or gravimetric, the pressure, and the temperature of operation for this process. For example, if the objective is to maximize volumetric delivery, then a density of around 1 g/cm 3 seems to be sufficient. In the case of gravimetric delivery optimization, lower densities seem best.
Ultimately, the complexity and the number of structural properties of MOFs make it difficult to perform optimizations as observed in GCMC simulations. In this work, these structural parameters were related. The relationship and effect of one structural property on another allows the analysis to essentially be simplified, and for fewer structural properties to be focused on for QSPR analysis and the proposition of simplified correlations. Again, as shown in Figure 2, these properties are related and clear trends were observed from the experiments. It was observed that pore volume and surface area are related. Increases in either parameter were accompanied by a decrease in density and increase in LCD. The specific values of LCD and density can be obtained from the 5D interactive data explorer assembled for CH 4 adsorption. The relationship between surface area and pore volume were obtained by a linear relation to help simplify analysis. While it is true that not all MOFs fit exactly to the proposed model, the majority did. The fit depends on whether advanced material-design strategies that have recently swept this area of research are used, including ultraporosity. Nevertheless, a clear framework and QSPR analysis equations are provided in this work to estimate the properties required with some confidence, and extrapolate required surface areas to meet current and future storage targets. For future work, it is proposed that the predictions from the models provided here be implemented and that experimental data and machine learning be applied to generate more complex models for the efficient design and synthesis of CH 4 adsorbents.

Conclusions
Quantitative structure-property relationship (QSPR) approaches, when applied to MOFs, allow for reasonable estimates to be made of CH 4 storage performance. Here, QSPRs were developed from experimental data and insights are provided on how to improve storage and deliverable CH 4 storage capacity based on material properties. Geometrical features of MOFs such as density, pore volume, and LCD, and their significance for CH 4 storage capacity, were assessed. One relationship identified is that CH 4 gravimetric storage capacity is proportional to BET surface area (r 2 > 90%). QSPR demonstrated the effect of open metal sites required for van der Waals forces involved in CH 4 adsorption. Guidelines are provided for optimal design of MOFs, including density and pore volume. With the recent achievement of the 2012 DOE target, the QSPRs presented here may allow for the prediction of structural descriptors for CH 4 storage capacity and delivery [28].
These predictions can be made at various pressures and a range of densities, and associated void fractions can be estimated based on the storage target desired. This was done using large-scale, experimental data, which has not previously been available to allow for comparison with GCMC model predictions. Linear relationships were found between gravimetric CH 4 delivery, CH 4 uptake, and BET surface area. Based on these relationships, surface area should be around 7000 m 2 /g or higher to achieve the 0.5 g/g gravimetric uptake target set by the DOE in 2012 at 65 bar and 298 K. Topology has an important effect, and the effect of functionalization was taken into account in the model. Pore volume and diameter play important roles, and the structural parameters were related to each other to simply the QSPR analysis and to allow for predictions to be made easily. Overprediction and underprediction of experimental data were observed in the GCMC simulations.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4352/10/8/700/s1, Figure S1: QSPR for CH 4 storage in MOFs. CH 4 volumetric uptake at 298 K and 35 bar is plotted against density and the largest cavity diameter (LCD) denoted by the diameter of the circles, Figure S2: QSPR for CH 4 storage in MOFs. CH 4 volumetric uptake at 298 K and 35 bar and 65 bar is plotted against density and the largest cavity diameter (LCD) denoted by the diameter of the circles, Figure S3. QSPR for CH 4 storage in MOFs. CH 4 gravimetric uptake at 298 K and 65 bar on the y-axis is plotted against BET surface area in m 2 /g on the x-axis, Figure S4. QSPR for CH 4 storage in MOFs. CH 4 volumetric delivery at 298 K and 35 bar on y-axis is plotted against CH 4 volumetric delivery at 298 K and 65 bar on the x-axis, Figure S5. QSPR for CH 4 storage in MOFs. CH 4 gravimetric uptake at 298 K and 65 bar (y-axis) is plotted against CH 4 gravimetric uptake at 298 K and 35 bar on the x-axis, Figure S6. QSPR for CH 4 storage in MOFs. CH 4 gravimetric uptake at 298 K and 35 bar (g/g) is plotted against BET surface area in m 2 /g, Figure S7. BET surface area on the y-axis (m 2 /g) plotted against pore volume of MOF on x-axis (cm 3 /g), Spreadsheet 1. Pore volume, uptake and deliverable capacity of selected MOFs.
Funding: This research was funded by United Arab Emirates University, grant no G00002618.