Metal-Organic Frameworks Characterization via Inverse Pulse Gas Chromatography

: The desire to customize the properties of a material through complete control over both its chemical and architectural structure has created a constant and persistent need for efﬁcient and convenient characterization techniques. Inverse gas chromatography (IGC) is considered a useful characterization method for probing the material’s surface properties, like its enthalpies of adsorption, which are the key stimulus components for their adsorption performance. Here, we conclusively review the signiﬁcance of a less common application of the IGC technique for the physicochemical characterization of metal-organic frameworks (MOFs), which are an innovative subclass of porous materials with matchless properties in terms of structure design and properties. This review focuses on the fundamental theory and instrumentation of IGC as well as its most signiﬁcant applications in the ﬁeld of MOF characterization to shed more light on this unique technique. methodology, K.Y. and O.S.; software, K.Y. and O.S.; validation, K.Y., O.S., Z.A. and M.E.; formal analysis, K.Y. and O.S.; investigation, O.S.; resources, K.Y. and O.S.; data curation, K.Y.; writing—original draft preparation, K.Y. and O.S.; writing—review and editing, K.Y. and O.S.; visualization, supervision, and


Introduction
To help identify the potential application of a microporous solid in separation or catalysis, theoretical or experimental determination of its basic adsorption properties and sorption behavior is a beneficial tool. Gathering information about the surface activity of these microporous materials provides a deeper understanding of their performance and continual development. Consequently, it is important to determine the surface properties that indicate the surface's ability to contribute to both dispersive and specific interactions.
The adsorption enthalpy and associated physicochemical parameters of a material may be derived directly from calorimetric measurements or indirectly by the adsorption constant temperature dependence obtained by volumetric, gravimetric, or chromatographic adsorption measurements [1][2][3][4][5][6][7][8][9]. Alkane adsorption molecular modeling via force fields for zeolites has been improved to such an extent over the last decade that great canonical Monte Carlo models can now be applied to forecast zeolites' adsorption behavior at high levels of pore filling [10][11][12].
Inverse gas chromatography (IGC) is recognized as a proper and fast tool for characterizing the surface properties of porous solids, with a minimal adsorbent amount needed for reliable calculations and shorter experimental times compared to other techniques used for thermodynamic characterization. One of the unique features of IGC is that no specific or complicated preparation procedure of samples is necessary. Relative to other surface energy analysis techniques, IGC requires minimal sample preparation [13]. One can easily and accurately use the IGC technique to characterize various crystalline, composite, amorphous, or even fibrous solids.
Without the need for any extrapolation, the enthalpies of adsorption measured with IGC and other methods show excellent consistency [14][15][16]. Since IGC has a superior efficiency at infinite dilution, it can differentiate between characteristic variations that other methods, such as conventional wettability measurements, do not identify [15]. For example, the probability of observing the differences associated with surface group orientations is one advantage of using IGC in researching acid-base specific properties [15].
Compared to analytical gas chromatography (GC), the mobile and stationary phases and roles are reversed in IGC. In IGC, the explored material represents the stationary phase. The injected known molecules within the mobile phase are considered probes used to investigate the physicochemical properties of the stationary phase.
The first attempt to introduce IGC was by Nobel Laureates Martin and Synge in 1941 [17]. In 1960, Kiselev gave the emerged IGC technique its name [18]. Furthermore, he was a key player in developing chromatography and the field of surface chemistry [19]. GC was identified by Smidsrod and Guillet as a versatile "indispensable analytical instrument" for applications other than defining mixture constituents [20]. In the 1970s, IGC became more familiar when it proved to be a promising method in investigating the surface and bulk properties of polymers and related materials [21]. Polymeric studies are still the most significant application of IGC judging by the large number of related publications [22].
In principle, IGC experiments can be conducted by either a pulse or frontal technique. For a pulse technique, a specific portion of the tested molecule is injected. With the help of a carrier gas, this pulse is transferred across the system through a column that is packed with the studied material. A chromatographic peak appears as a result of the adsorption and desorption processes within the column. On the other hand, in the frontal technique, the tested molecule is injected perpetually into the system, and consequently, a breakthrough curve will be the output chromatogram [21]. The advantage of the frontal technique is that, due to its continuous operation, equilibrium can be achieved, whereas the pulse technique relies on the assumption of rapid equilibrium between the probe and the solid surface.
In the case of the continuous flow method, which is a classical form of frontal technique [23], the column can be detached from the gas flow if the concentration is changed and afterwards attached again when the flow is balanced. Then, a positive peak will appear for a higher concentration and adsorption process, while a negative peak will indicate a lower concentration and desorption process. The peak area is directly proportional to the amount of adsorption or desorption. This technique is utilized in many instruments conceived for the dynamic estimation of Brunauer-Emmett-Teller (BET) surface area [15]. However, this approach cannot be considered an IGC technique since retention time is not a parameter in the experimental calculations [15]. In comparison, classical frontal experiments evaluate the nature of retention and, therefore, are considered an IGC technique.
A significant proportion of publications utilized pulse instead of frontal approaches because they are quicker, more convenient, and provide higher accuracy, especially if there are poor interactions between the tested material and probe. The frontal approach is a suitable alternative if a slow equilibrium is the topic of examination.
Some research groups tend to prefer the term "pulse gas chromatography" over the term "inverse gas chromatography," while others used the full expression "inverse pulse gas chromatography." We have noticed that these different terminologies were sometimes misleading, indicated by few cross citations. Furthermore, subsequent studies claim that "there have been almost no studies using the IGC with MOFs" [24].
After the IGC technique was recognized as a sensitive and simple tool for investigating physicochemical properties, it was applied for different types of materials, including non-volatile and solid materials with numerous morphologies. For instance, it was applied to evaluate the impact of grinding on the surface and bulk features of drug molecules [25], which is important for medicinal industry improvements. IGC has also been used to deconvolve the chemical and structural properties of carbon nanotube surfaces [26]. Additionally, IGC was used to characterize polymers, minerals, modified silicas, surfactants, and zeolites [22,27]. However, in the case of metal-organic frameworks (MOFs), which are considered a recently introduced porous material, few studies have used IGC for the characterization of their different properties.
Metal-organic frameworks (MOFs) are a subclass of porous crystalline materials constructed by linking metallic nodes like single-ion or clusters with organic linkers. They have been developed in the last twenty years as potential candidates for solving many of the challenges associated with different applications, such as gas storage, separation, catalysis, and sensing [28][29][30], by means of the anticipated on-demand control of their pore-system (i.e., pore-aperture, shape, size, and functionality). MOFs are considered organic-inorganic hybrid porous materials and have supreme chemical, structural, and topological adaptability [31,32]. Their pore system is tunable via different design and engineering strategies [33]. MOFs can be efficiently designed with the proper properties in terms of pore aperture size, controlled by the organic linkers linked through metallic-nodes and/or shape (channels or cages), to achieve the targeted application. Successful exercising of reticular chemistry ensures the suitable inorganic and organic building blocks, encompassing defined geometrical information that ciphers a designated net, are constructed into a hybrid material framework with the desired topology.
MOFs with a wide range of porosities and chemical functionalities have been synthesized and explored in many applications such as separation and catalysis [30,34]. Separation, for example, can be driven by the difference in affinities and interactions between the MOF framework and the targeted probe molecules, where the structure and chemical functionality of the MOF pores can be varied, for example, via the organic linker length, functionalization, or integration of coordinatively unsaturated metal sites [35][36][37]. MOFs have demonstrated outstanding performances as separating agents either as adsorbents or membranes as a result of their well-defined chemical and structural properties in terms of pore size, functionality, and rigidity/flexibility, which correlate with the targeted adsorbates' intrinsic properties like functionality, size, shape, etc.
Porous materials like MOFs that possess essential features for the selective separation of targeted analytes over others represent a potential alternative to energy-intensive techniques like distillation. An adsorption-based process can afford separation via thermodynamics, where the MOF acts as an adsorbent and shows a comparatively high affinity for a particular molecule due to the enhanced interactions between the adsorbates and the MOF material.
Due to their unique surface properties and separation abilities, much research is focused on using MOFs in applications like sensing, catalysis, or as separating agents for carbon-capture, hydrocarbons separation, xenon/krypton, and many other essential separations.
This study sums up past work concerning the utilization of IGC as a characterization tool to investigate the physicochemical features of MOFs and highlights IGC as a useful, revolutionary, and almost undiscovered technique for investigating and predicting the adsorption properties of MOFs.

Instrumentation and Methods
In analytical GC, the tested sample mixture, whether gas or liquid, is injected into a well-defined analytical column to separate the composition of that mixture for subsequent qualitative and quantitative analysis. The mixture is separated through different interactions of varying strengths between the analytes and the stationary phase. In IGC, a single known probe is usually injected into a packed column to study its retention behavior and correlation with the physicochemical properties of the packing material. Various probes with different well-known properties are used.
The composition of the IGC instrument is similar to that of the analytical GC technique ( Figure 1). An IGC packed column may be glass, fused silica, copper, stainless steel, or Teflon. Wall-coated open-tubular long columns have also been used in several types of research. An inert carrier with a steady flow conducts the probe through the column at a constant rate. In IGC, the main collected parameter is the retention time, which has been effectively used to evaluate the solid surface's dispersive properties and potential specific interactions [15,16,22]. IGC experiments can be carried out under two different conditions: infinite dilution and finite concentration. Infinite dilution has been identified as appropriate for the heat of sorption and surface energy examination of particles [39]. Infinite dilution at low surface coverage indicates a minimal amount of the tested probe molecule is injected. A low concentration of the adsorbate molecules indicates that the interactions occur with the higher energy sites available at the surface of tested materials. Furthermore, probeto-probe interactions are considered negligible due to the limited amount of the probe. Therefore, Henry's Law is applied [40], resulting in a linear isotherm of adsorption and a typical chromatographic peak [15,22]. Infinite dilution is the dominant condition in the IGC technique.
At finite concentration, a large amount of test probe is injected into the system; hence, interaction with all surface sites is expected [41]. This technique can be used to determine solids' heterogeneity of surface energy [42], or to study the role of both mesopores and micropores in adsorption and thermal desorption methods [43]. However, finite concentration IGC is out of the scope of this review.

Theory
The surfaces of solids are normally characterized by their dispersive and specific properties. By measuring the dispersive components of the surface free energy, the dispersive nature is obtained. Specific properties are defined by the factors that calculate the surface propensity to exhibit Lewis acidic or basic features, thus indicating an electron acceptor or donor nature. The total solid surface energy is represented via the summation of both the dispersive and the specific properties.

Dispersive Properties of the Surface
The retention volume (V N ) occupied by a carrier gas from the injection valve to the maximum of the chromatographic peak through the column is estimated using the following equation [44]: Here, t R represents the retention time and t M is the dead time, which is usually experimentally calculated by using an unretained methane gas probe. F a represents the volumetric flow rate, T represents the column temperature, T a represents the ambient temperature, and j represents the James-Martin gas compressibility correction factor, which can be calculated from Equation (2) [45]: where P represents P i /P o , with column pressure (P i ) and outlet pressure (P o ). A specific retention volume (V g ) represents a definite value of the material (as a particular guest molecule-material) affected by several factors, including particle size, packing geometry, and temperature, among others. This is considered equivalent to the point of inflection in breakthrough experiments as an example of frontal techniques [22,46]. The following relationship (3) is used to calculate V g : where w represents the weight of the stationary phase, and P w represents the vapor pressure of water at the used T.
Conventionally, IGC basic equations are derived from V g ; however, the retention volume (V N ) is very often utilized as well. Additionally, V g could be estimated from the net V N through the following equation [47]: Henry constants (K ) represent the slope of isotherm in the linear region. They reflect the proportion of carrier inert gas that is needed to elute the corresponding analytes' maximum concentration within the porous material. Henry constants are usually determined using the retention volume V N and considering ideal analyte behavior [21,48]: where T represents the experimental temperature of the column, and the ideal gas constant is denoted by R.
The ability of the investigated material to discriminate branched alkanes from linear ones is a result of the dissimilarities in the tested probes' physicochemical properties or the material's selective shape properties. To determine the actual reason for the separation, the Henry constants logarithm could be compared to the vapor pressure logarithm plot. For linear alkanes, the relationship is linear; but after the branched alkanes are positioned in that line, the material cannot separate those probes based on their shape and vice versa [49].
The dispersive component of the surface energy (γ D S ) reflects the capability of a solid surface to demonstrate Van der Waals forces or London interactions, and the surface heterogeneity of the probed material. It is a typical factor that estimates the energy of the surfaces of solids and is deemed the equivalent of surface tension in liquids [50]. Two major techniques are employed in order to determine the dispersive interactions between the probes and the tested solid phase: Dorris-Gray [51] and Schultz et al. methods [52]. Dorris-Gray proposed a relationship to estimate γ D S , relying on the Fowkes's approach [50] regarding the used stationary phase and adhesion of methylene group, as shown in Equation (6) [51]: (6) where N represents Avogadro's number, the methylene group cross-sectional area is a CH2 (≈6 Å 2 ), V N,n+1 represent the net retention volumes for n-alkanes with carbon number n + 1 and V N,n for carbon atoms with only n, and γ CH2 represents the surface dispersive free energy of only CH 2 groups (as, for example, in linear polyethylene), in a solid material calculated at a given T ( • C) using the following equation [52]: The Fowkes approach was also the base of the Schultz et al. method, in which the effect of adhesion represents the dispersive free energy between the stationary phase and the probe. Hence, Schultz et al. introduced a relationship between the surface dispersive free energy and V N as follows: where a represents a cross-sectional area and γ d l is the solute probe dispersive free energy. For n-alkanes, a linear relationship is obtained when RT lnV N is plotted against a γ d l . The dispersive free energy for the solute probe can be evaluated from the slope.
The enthalpy change of adsorption (∆H ads ) can be easily obtained by plotting the logarithm of the retention volume, the logarithm-specific retention volume, or Henry's constants gained from isothermal GC measurements versus the inverse temperature. The adsorption enthalpy change ∆H ads (the differential heat of adsorption) represents the interaction strength between the adsorbed molecule and the adsorbent atoms at zero surface coverage: The standard free energy change of adsorption (∆G ads ) is determined from the V N or V g of the vapor probe using Equation (10) [44]: where m represents the adsorbent material mass in the column, π 0 represents the reference two-dimensional surface pressure, S represents the specific surface area of the adsorbent, P • represents the vapor pressure of the adsorbate in the gaseous state, which can be estimated using Antoine's equation: where A, B, and C are the Antoine coefficients from the handbooks and t is the temperature in Celsius. The standard entropy change of adsorption (∆S ads ) is calculated at zero Henry surface coverage using the following equation:

Specific Properties of the Surface
Specific as well as dispersive interactions occur simultaneously for polar probes in the absence of metallic, magnetic, electrostatic, or acid-base surface interactions, which reflect most specific interactions [16]. Polar probes' free energy of adsorption has two types of interactions: dispersive interactions, which are symbolized by ∆G D ads and correspond to London forces, and specific interactions, which are symbolized by ∆G SP ads and are attributed principally to Lewis acid-base contributions: The free adsorption energy-specific component (∆G SP ads ) is estimated as the difference between the polar compound free energy of adsorption (∆G polar ) and the free adsorption energy of hypothetical alkane (∆G re f ) (Equation (14)). The hypothetical alkane must have the same assigned property as a polar probe (e.g., vapor pressure in the Papirer method [53], boiling point in the Brookman and Sawyer method [54], or a γ d l value in Schultz and Lavielle method) ( Figure 2): where The specific component of the adsorption enthalpy (∆H SP ads ) is determined from the ∆G SP ads dependence on temperature, reflecting the Gibbs Helmholtz equation: where the polar probe's specific component of adsorption entropy is represented by ∆S SP ads . Hence, ∆H SP ads can be calculated from the slope of ∆G SP ads T versus 1 T . The specific enthalpy of adsorption of a probe on the tested solid surface is correlated to the solid's acid-base properties. The acceptor constant (K A ) represents the acidity, and the donor constant (K D ) represents the basicity of the tested solid phase; both are determined using Equation (16): where AN* is the modified acceptor number of Riddle and Fowkes [55], and DN is Gutmann's donor number [56]. The value of K A is represented by the slope of the linear relation between − ∆H SP ads AN * vs. DN AN * , and the intercept corresponds to K D .

Applications of IGC for the Characterization of MOFs
IGC is considered a very effective surface energy analysis tool for the characterization of various materials' surface properties. Here, we review recent applications of IGC in analyzing the surface properties of MOF materials. These applications include gases, hydrocarbons, and solvents, which reflect the wide-ranging advantages that can be achieved by applying the IGC technique ( Table 1). The IGC evaluations reviewed here were are carried out using conventional GC at infinite dilution, which is usually equipped with a flame ionization detector (FID) or a thermal conductivity detector (TCD).
The studied MOFs were packed into either metallic or glass columns of 2-80 cm length, with 30 cm being the most commonly used length, and a standard internal diameter (ID) of 3.18 mm with actual utilization of the columns ranging from 0.53 to 5.3 mm ID (Table 1). Whether prepared or purchased, the adsorbent powder was either packed directly into the column or pressed under high pressure (~600 bar) into disks, then broken into fragments and sieved to the proper fraction for a higher homogeneity.
Packing of the tested material was achieved either through vacuum suction of the dry powder or by pushing the suspended powder into the empty tubing. In either case, good packing was usually accomplished using mechanical vibration or an ultrasonic water bath; thereafter, the column ends were plugged with glass wool to avoid the loss of adsorbent.
Silica wall-coated open tubular capillary columns have been used successfully for the IGC calculations for the coated material. Capillary columns coated with MOFs have been used in two forms: either a controlled approach using secondary building units (SBUs) for in situ preparation of the MOF material [24], or as a dynamic coating after preparation [73]. Estimating the physicochemical properties of MOFs incorporated into an organic monolithic capillary column as a composite is an unusual platform for IGC studies; however, it was successful [71].
Generally, high purity nitrogen or helium (≥99.99%) is applied as the carrier gas with flow rates ranging from 20 to 60 mL min −1 ; exceptions include 1 and 5.3 mL min −1 flow rates, which were used in the case of open tubular coated columns [24,73], 5.5 mL min −1 for lower molecular weight N 2 /O 2 gaseous probes [78], and about 0.6 mL min −1 for a monolithic column [71].
Prior to the measurements, the most fulfilling activation procedure of the studied material within the column is performed by heating the column to a suitably high temperature, using a slow ramp heating under a constant flow of the inert carrier gas overnight, that maintained the required thermal stability of the material.
Then, a trace amount of the probe molecules is then injected as a vapor phase. The adsorption data is collected under a wide range of temperatures, starting from room temperature up to 300 • C, depending on the nature of the adsorbed probes, with the retention time of the separated probes being the primary information collected during an IGC experiment.

HKUST-1
HKUST-1 (HKUST = Hong Kong University of Science and Technology), known also as MOF-199, was reported initially by Chui et al. [80]. This MOF entails copper(II)-dimer forming paddle-wheel SBUs, which form coordination bonds with four carboxylate groups from the organic linker trimesic acid (benzene-1,3,5-tricarboxylate (BTC)), and has the formula Cu 3 (BTC) 2 ·xH 2 O. The three-dimensional extension of these paddle-wheel units forms a three-dimensional structure with an underlying tbo topology ( Figure 3). This MOF has a relatively high surface area in the range of 1500-2100 m 2 g −1 .   Other authors have identified this MOF structures with only two types of pores: large 9-12 Å pores connected by 7-8.3 Å windows, and small 5 Å pores accessible by 3.5 Å windows [80,83,84]. In the hydrated HKUST-1 structure, a water molecule is coordinated at the axial position of the paddlewheel unit and can be removed upon an activation process, without any change in MOF crystallinity or loss of porosity. An unoccupied coordination site at the Cu 2+ atoms is called a coordinatively unsaturated site (CUS) and can be accessed by other probe molecules [82]. One metal ion coordination site is unsaturated and therefore open to Lewis basic molecules [80,85]. HKUST-1 possesses Lewis acidic centers from the copper(II)-SBUs, and Lewis basic centers from π electron-rich benzene rings in the BTC molecule.
The Fe-BTC MOF is analogous in structure to HKUST-1 and is synthesized as a porous solid, however, in a gel-like form (surface area is about 1300-1600 m 2 g −1 ) [86]. Fe-BTC MOF, unlike HKUST-1, has not been thoroughly explored or studied because its structure is not highly crystalline that may cause distortion of its pore geometry. Its morphological properties were studied in comparison to MIL-100 (Fe) (MIL = Material of Institute Lavoisier) [87], and Fe-BTC MOF has been found to encompass 21.2% of Fe 3+ as a nodal metal and the BTC linker in a distorted form of MIL-100(Fe) MOF [87,88]. While its composition is similar to that of MIL-100 (Fe), with respect to many applications, Fe-BTC MOF cannot be compared directly with MIL-100 (Fe). HKUST-1 and Fe-BTC MOF were commercialized under the commercial names Basolite C300 and F300, respectively.
Finsy et al. [49] demonstrated a possible use of HKUST-1 MOF for adsorptive separation based on shape selectivity. They found that linear alkanes were adsorbed preferably over methyl-branched alkanes. The high adsorption enthalpy and low Henry constants of the branched alkanes relative to those of their linear isomers suggested shape selectivity, as the Henry constants decreased with branching degree.
Shape-based discrimination was used to examine the logarithm of the vapor pressure versus the logarithm of the Henry constant plot (Figure 4). A linear relationship with 2-methyl-branched alkanes on this trend line was observed for linear alkanes, although not all alkanes branched like 3-methyl-and dimethyl-alkanes behaved similarly. This suggests that HKUST-1 cannot be efficiently used for separating linear alkanes from 2methyl-branched alkanes based on their shape, but it can separate linear alkanes from 3-methyl and dimethyl-branched alkanes. Münch and Mertens [24] applied the IGC method to evaluate the donor properties of the HKUST-1 linker BTC and the acceptor properties of the CUS sites. These experiments were executed using a GC wall-coated silica capillary column of 30 m length and 0.25 mm inner diameter, unlike the majority of IGC studies, which typically utilize a short-packed column (Table 1) [24].
The chromatogram indicates that benzene showed a longer retention time than cyclohexane ( Figure 5), but a lower boiling point. This suggests the existence of specific interactions, for example, π-π interactions between the benzene ring π-electrons of the BTC linker in HKUST-1 and the π-electrons in the benzene moiety of the analyte. Benzene also showed a higher value of ∆H ads than cyclohexane, implying a stronger interaction. A thermodynamic investigation of a group of ethers with different branches, including THF ( Figure 5), revealed a significant effect of the geometry of the analytes on adsorption. It was concluded that there was a combination of two factors: open metal site accessibility, which allowed interaction with oxygen free electron pairs; and an entropic effect caused by the degrees of freedom of the alkyl chains [24].
The Gutmann donor and acceptor numbers of a particular probe analyte permitted a quantifiable assessment of the Lewis acidity and basicity sites in the HKUST-1 structure. The HKUST-1 acceptor properties had a value of 1.6 for K A , and the surface donor properties had a value of 0.4 for K B . From these numbers, it was concluded that HKUST-1 can be categorized as primarily Lewis acidic, which is attributed to the prevalence of CUS [24]. In another study by Autie-Castro et al. [59], n-butylamine thermodesorption analyses were performed using IGC to compare HKUST-1 to its analog Fe-BTC with respect to their acid-base properties. Samples of the proposed materials were saturated with n-butylamine to quantify both sites (acidic and basic) from the desorption curves of n-butylamine using the n-butylamine thermodesorption technique.
The IGC measurements yielded 0.81 and 0.60 as the acidity constants (K A ) and 0.09 and 0.28 as the basicity constants (K B ) for Fe-BTC and HKUST-1, respectively. Moreover, the n-butylamine thermodesorption study indicated a total acidity of 2.06 mmol g −1 for Fe-BTC and 0.49 mmol g −1 for HKUST-1. This behavior may be attributed to the superior polarizing rule of Fe 3+ compared to Cu 2+ due to the higher charge. As a result, both MOFs, HKUST-1 and Fe-BTC, were verified as being predominantly Lewis acidic.
Specific interaction contribution values of 4.1 and 3.2 kJ mol −1 at 200 • C were acquired for the benzene moiety in HKUST-1 and Fe-BTC, respectively. Furthermore, the contributions of the specific interactions regarding THF at 200 • C were found to be 14.5 kJ mol −1 for HKUST-1 and 17.9 kJ mol −1 for Fe-BTC. HKUST-1 has a mildly basic surface character, whereas Fe-BTC has a more prominent acidic character, resulting in more benzene (acidic) interactions on the HKUST-1 material surface and more THF (basic) interactions on the Fe-BTC material surface.
Autie-Castro et al. [57] measured the ∆H ads on both Fe-BTC and HKUST-1 for nalkanes, propane, propylene, benzene, toluene, and xylene. Dispersive interactions were the dominant adsorption mechanism for n-alkanes, and ∆H ads was found to increase with the number of carbon atoms. The electrostatic interaction of propylene and aromatics through its dipole and/or quadrupole moments via the guest-host π-π coupling led to a higher ∆H ads value for propylene than propane.
Higher ∆H ads values were found for HKUST-1 than Fe-BTC, due to the ordered crystal structure and molecular ease of access to the porous structure of HKUST-1. In the case Fe-BTC, which has does not have a well-defined crystal structure, local deformations of its flexible organic ligand may account for the restricted accessibility of pocket-type pores, where small probe molecules were restrained in a small volume and thus contributed more to the adsorption forces of dispersive interactions.
Gutiérrez et al. [58] compared HKUST-1 and Fe-BTC in more depth by comparing the ∆H ads obtained from IGC and the ∆H imm (the enthalpy of immersion) obtained from immersion calorimetry. For the three studied adsorbates (n-heptane (HEP), methylcyclohexane (MCH), and toluene (TOL)), the interaction strength was higher for Fe-BTC than HKUST-1. This was a result of the better accessibility to the active sites in Fe-BTC, as it has a lower fraction of microporous surface area. Although Fe-BTC does not have a specified crystalline structure, the open metal sites can favor electrostatic interactions with the adsorbates, leading to an increase in the overall affinity of the interaction.
The ∆H imm followed the same order as the enthalpy of adsorption, i.e., TOL > HEP > MCH, and was consistent with the ∆H ads measurements for HKUST-1 ( Figure 6). However, unlike ∆H ads , the ∆H imm values for Fe-BTC were in the inverse order with respect to the smallest size of the adsorbates: TOL (4.012) < HEP (4.014) < MCH (4.982) (Figure 7) [89]. This suggests that these electrostatic interactions are less crucial and the accessibility of the adsorbates to the pores is more relevant. The higher cross-section area of the adsorbent n-alkanes on the surface of both MOFs makes the dispersive interaction of the surface energy produced more prominent for HKUST-1 than for Fe-BTC. Additionally, the highest specificity of the interaction was also reported on HKUST-1, which is contrary to what was expected, as indicated by the specific interaction parameter (I sp ), which reflects on the surface acid-base nature. Another research group [38] utilized IGC as a prescreening tool. They effectively examined the affinity of MOFs for targeted probes in the gas phase. For the commercially available MOFs, like HKUST-1 (Basolite C300), Fe-BTC (Basolite F300), MIL-53 (Basolite A100), and ZIF-8 (Basolite Z1200), a preconcentration process for sensitive sensing of nitroalkanes was developed, and the results were associated to those of the Tenax TA porous archetype.
In this study, HKUST-1 showed an enrichment factor 109-fold higher than that of Tenax, as revealed by the high value of ∆H ads for nitromethane (Figure 8), and a dispersive element of surface energy values (γ D S ) above 100 mJ m −2 . As predicted, the enthalpy of adsorption values for nitroalkanes on HKUST-1 (relatively basic) were higher than those on Fe-BTC (Figure 8) since the nitroalkanes were slightly acidic due to the electronwithdrawing effect of nitro groups. To correlate IGC measurements to a frontal technique, the results of IGC experiments were compared to micro-breakthrough ones with downstream thermal desorption experiments. HKUST-1 also displayed impressive improvement factors above 2000 while sampling at a nitromethane source of 1000 ppm.

MIL-Based MOFs
MIL-47 is a rigid MOF composed of infinite chains of VO 6 octahedra connected via the dicarboxylate groups of terephthalate organic (H 2 BDC) ligands. In this structure, a microporous framework is created having one-dimensional diamond-shaped channels with pore apertures of 7.9 Å × 12.0 Å (Figure 9) [90]. The adsorption properties of MIL-47(V) were evaluated by Finsy et al. [60] in a combined IGC experimental and Monte Carlo simulation study. Their report showed that both linear and branched alkanes can adsorb into the pores of MIL-47(V) without any steric hindrance. The dependency of adsorption enthalpy and the Henry constants on the carbon number was significantly correlated between the simulations and experiments. n-Hexane (50.6 kJ mol −1 ) was adsorbed more strongly than benzene (43.4 kJ mol −1 ) and cyclohexane (45.2 kJ mol −1 ) because it has more hydrogen atoms. The adsorption energies of the studied nonpolar molecules were governed by van der Waals interactions, and the adsorption of the benzene probe molecules was further affected by Coulombic interactions.
Similar work on MIL-47(V) to evaluate its physicochemical properties via IGC was done on the adsorption of substituted aromatics in addition to other heterocyclic molecules [61]. Quantitative structure−property relationships (QSPRs) between the adsorption and molecular properties were established using principal component analysis and partial leastsquares regression. With respect to heterofunctionalized probes, the adsorption properties tended to associate with the polarization and/or electronegativity of the substituent. The effect of the resulting dipole moment originating from the substituent placing was observed [61].
The experimental results indicated that the molecular polarizability of the halogen series perfectly fell on the alkyl aromatics series (Figure 10), suggesting that the observed properties were likely due to increased polarization and not to higher electronegativity of the halogen atoms. The oxygen-containing functional groups like nitro, methyl ester, and carbonyl were removed from this pattern. Thus, the expanded π-system enabled further delocalization of π-electrons from the substituent and the aromatic ring. Consequently, the electron-rich substituent attached more tightly to the MOF material than the alkyl fragments. Flexible MOFs are a special class of MOF materials that can reversibly alter their framework when guest molecules are introduced or removed. This leads to unique properties such as the breathing effect and the gate-opening phenomenon, in which pores contract or extend during adsorption [91][92][93]. Examples of breathing materials are MIL-53 [94] and its functional derivatives [91,[95][96][97]. MIL-53 is constructed from MO 4 (OH) 2 octahedra (where M can be Fe 3+ , Cr 3+ , or Al 3+ ) with 1,4-terephthalate (BDC) as the linker. In this way, MIL-53, which is also known under the commercial name Basolite A100, is a crystalline material having one-dimensional diamond-shaped pores with an 8.6 Å free diameter [91,98].
Couck et al. [62] utilized IGC to study the adsorption properties of amino-MIL-53 (Al). Amino-MIL-53 (Al) also has one-dimensional channels, similar to the amino-free MIL-53 analog, with a pore diameter of nearly 7.5 Å. The amino groups are directly inside the pores as a result of using 2-amino terephthalic acid as a linker rather than terephthalic acid [99], which leads to an enhanced affinity towards CO 2 [95].
Couck et al. found that the adsorption enthalpy and Henry's adsorption constants of iso-alkanes were considerably smaller than those of linear alkanes, highlighting the shape selectivity features of amino-MIL-53. Furthermore, the amino groups exposed in the pores of amino-MIL-53 enhance the electrostatic contribution of molecules with double bonds. Although double bonds contribute to total interactions, van der Waals interactions were found to be the dominant element for the adsorption interactions of molecules. Using a simple model, Couck et al. predicted that hydrogen atoms and double bonds contributed to these adsorption properties. The determined confinement factor Z pointed to an increased degree of confinement in amino-MIL-53(Al) compared to MIL-53 and MIL-47. Amino-MIL-53(Al) also exhibited a smaller pore aperture, which is due to either the presence of amino groups or to the framework shrinking as a result of the interaction between amino functional groups and hydroxyl groups.
Sun et al. demonstrated that the enhanced adsorption activities of CO 2 over CH 4 [63] and CH 4 over N 2 [64] over Al-BDC MOFs due to larger pore size. Their XRD patterns were different from those of MIL-53(Al). There was a slight difference compared with the standard MIL-53(Al), which may be as a result of the presence of some free H 2 BDC linkers within the pores and the metastable phase of BDC and Al 3+ as a result of the different conditions of crystal growth. The pore size of Al-BDC MOFs increased from 0.86 (AlBDC Micro ) to 2.58 nm (AlBDC Meso ) via a simple washing and drying process. Subsequently, an analysis of the adsorption equilibrium and thermodynamics of the different gases, CH 4 , CO 2 , and N 2 , was done via volumetric−chromatographic and IGC methods [62].
The initial adsorption heat was estimated to be 11.5 kJ mol −1 for CH 4 and 25.2 for CO 2 on AlBDC Meso MOF, which is approximately 25% less than on AlBDC Micro MOF. These data indicated that enlarged pore size not only increased the CO 2 adsorption selectivity over that of CH 4 , but also decreased the value of ∆H ads , suggesting that it should be the preferred method for selecting a satisfactory absorbent for CO 2 /CH 4 separation [63]. The heats of adsorption of CH 4 , N 2 , and CO 2 were calculated as 15.3, 11.5, and 32.2 kJ mol −1 , respectively, which is aligned with the gravimetric study [64].
MIL-101(Cr) with the formula (Cr 3 O(H 2 O) 2 F(BDC) 3 ) has a very large surface area and cubic cell volume with two kinds of cages: smaller 29 Å cages with pentagonal faces and larger 34 Å cages with both pentagonal and hexagonal faces [100]. The adsorption behavior of MIL-101(Cr) towards p-xylene was investigated using IGC [65]. These data showed that the specific component of the free energy of adsorption (∆G spec ) values of the polar probes were negative, suggesting spontaneous adsorption. The adsorption specific component of the free energy ∆G spec for polar molecules on MIL-101(Cr) were found to be in the order p-xylene > toluene > acetone. The ∆G spec for polar probes was inversely proportional to their dipole moment, although the authors have not provided an explanation for this phenomenon [65].

Isoreticular Metal-Organic Frameworks (IRMOF)
The IRMOF structure consists of Zn 4 O tetranuclear clusters linked by various rigid linear ditopic linkers, giving the obtained frameworks different pore aperture sizes and forming a cubic structure with connected channels (Figure 11) [101].  [103] were studied using IGC, breakthrough adsorption, and static vapor-phase adsorption [66].
IRMOF-1 and MOF-monoclinic have been shown to have different selectivity for the separation of the studied probes. On IRMOF-1, EB was found to be eluted first; however, xylene isomers were eluted similarly. The equilibrium constants controlled the elution order of the xylenes and EB on IRMOF-1. On the other hand, MOF-monoclinic exhibited a preferable adsorption for the p-xylene isomer, which was believed to result from some hindrance due to pore size that governed the diffusion of the isomers in the MOF micropores.
The McReynolds constants of IRMOF-1 showed that the nonpolar nature of that material was advantageous for the thermodynamic separation of the xylene and EB isomers. The intermediate polarity of MOF-monoclinic was attributed to the higher accessibility of its SBUs to isomers, and its pore shape was found to play a more significant role in separation than in IRMOF-1.
The properties of IRMOF-1 can change markedly based on its synthesis, activation method, and contact with atmospheric humidity [104][105][106][107]. Different IRMOF-1 samples prepared with different synthetic schemes were used to examine the impact of structural decomposition on adsorption behavior [67]. The infinite-dilution assumption in the adsorption isotherm linear segment was confirmed as valid through the independence of the chromatographic peak retention times for nonpolar n-alkanes with respect to concentration. On the other hand, the compounds having considerable hydrogen-bond basicity groups, such as diethyl ether, acetone, ethyl acetates and THF, were observed to vary from this ideal behavior.
The independence of the amended retention volume from the flow rate showed that the mass transfer did not affect the retention time. This implies that the amended retention volumes and all related calculated values indicate equilibrium behavior. The calculated values for enthalpy and entropy were more significant than those predicted by Monte Carlo simulations [108]. The presence of a nonporous phase in IRMOF-1 appeared to result in active sites with higher energies than those existing in standard IRMOF-1. The obtained linear relationships amongst the carbon number, the enthalpy of adsorption, and the free energy of adsorption were attributed to the additive nature of dispersive interactions.
The dispersive component of the surface energy (γ D S ) values for IRMOF-1 samples was calculated to be between 48 and 91 mJ m −2 , which is noticeably lower than the values for zeolites and activated carbons, which have been frequently recorded in the range of 200-500 mJ m −2 [51,[109][110][111]. This indicated a comparatively inactive surface for IRMOF-1, which was believed to be due to its large cages. Most of the interactions between the n-alkanes and the organic part of the framework were due to steric constraints near the metallic corners and the size of the n-alkanes. The deformation polarizability was much better suited to measuring negative values for free adsorption energy than for vapor pressure ( Figure 12). Most of the data points outside of the nonpolar line were related to the strong alkali and chlorinated molecules. From this fact and the observation that all the aromatics tended to lie on the nonpolar line, it was concluded that the basicity component of adsorption of the hydrogen bond prevailed for IRMOF-1. Similarly, the influence of the π-bond polarizability seems to be less significant [67].
In another study, Gutiérrez et al. [68] evaluated the interaction of some volatile organic compounds (VOCs) with three IRMOFs (IRMOF-1, IRMOF-8, and IRMOF-10) using IGC in order to assess the adsorption and link it to the morphological structure. The organic linkers of the three MOFs are terephthalic acid (IRMOF-1); 2,6-napthalene dicarboxylic acid (IRMOF-8); and 4,4-diphenyl dicarboxylic acid (IRMOF-10). Thus, they vary only in pore size. Overall, the strength of adsorption was enhanced in the presence of lattice defects and increasing numbers of carbon atoms in the organic linkers; accordingly, IRMOF-10 had the strongest interaction with the VOCs, followed by IRMOF-8 and finally IRMOF-1. Gutiérrez et al. deduced that the dispersive component of the surface free energy and the thermodynamic compensation effect of the n-alkanes exhibited the same centers of adsorption on all three MOFs. Furthermore, they witnessed that the existence of π electronrich aromatic rings or double bonds tended to increase the specific interaction element through a preferred interaction with the aromatic rings. Consequently, the specificity of the interactions was connected to the chemistry of the organic ligands rather than to the structure. Figure 12. Correlation between a VOC deformation polarizability and the free energy of adsorption (retention volume) at 200 • C as measured by IGC shown for (a) IRMOF-1 sample 1; (b) IRMOF-1 sample 2, and (c) IRMOF-1 sample 3. In these plots, the solid circles represent the series of n-alkanes, and the triangles represent other measured analytes. Reprinted with permission from [67]. Copyright (2010) American Chemical Society.

Zeolitic Imidazole Frameworks (ZIFs)
ZIFs have isomorphic topologically with zeolites. They are composed of tetrahedrally linked transition metal ions like zinc, cobalt, etc., that are coordinated by imidazolate-based ligands. The metal-imidazole-metal angle is 145 • , which is analogous to the Si-O-Si angle present in zeolites. ZIF-8 is one of the most utilized members of this family, consisting of 2-methylimidazole ligands connected with Zn 2+ ions to form eight sodalite (sod) cages that construct an overall cubic structure ( Figure 13). ZIF-8 has been commercialized under the commercial name Basolite Z1200. Luebbers et al. [69] applied IGC to investigate the adsorption behavior of many gases and VOCs on ZIF-8. Such highly halogenated compounds, branched alkanes, and aromatic compounds eluted rapidly with small tailing peaks and retention times, demonstrating a molecular sieving behavior that inhibited them from accessing the narrow windows of the MOF structure. On the other hand, the retention behavior of many of the studied compounds, especially n-alkanes, implied an unexpected flexibility of the ZIF-8 structure, even with n-alkanes with kinetic diameters similar or larger than 0.43 nm [112], significantly larger than the recorded ZIF-8 pore window size of 0.34 nm [113,114].
Owing to the mainly organic character of the reachable binding sites, ZIF-8 exhibited a lower surface energy than several other MOFs and microporous materials [113]. H-bond interactions were considered a dominant factor, as in the case of acids, amines, and alcohols that bind too strongly to be adequately studied in this method. Nevertheless, the adsorption of other probes, such as ketones, tends to be far less effective owing to the absence of such strong acid sites. Additionally, electrostatic interactions were essential, as the definite component of the adsorption free energy was observed to be significantly linked to the dipole moments of the adsorbates. In addition, a major increase in the adsorption of alkenes over alkanes was found to originate from the polarizability of the π-electrons, in addition to the interactions resulting from the induced dipole with the ZIF-8 structure.
The molecular sieving behavior and adsorption properties of ZIF-8 were also highlighted by Yusuf et al. [71,115]. They examined the impact of integrating various percentages of ZIF-8 into a BuMA-co-EDMA (butyl methacrylate-co-ethylene dimethacrylate) monolithic capillary column. The adsorption enthalpy values of both the neat monolith and the ZIF-8@(BuMA-co-EDMA) composite monolith did not reach the limit of physical and chemical adsorption (62.8 kJ mol −1 ) [116], suggesting the presence of only physical alkane interactions. The proximity to the heat of liquefaction for n-alkanes [117] and its change in enthalpy of adsorption showed that secondary weak intermolecular forces dominate the interactions through the monolith-n-alkane interface [118]. The addition of ZIF-8 particles to the monolithic organic material increased the dispersive element of the surface energy (γ D S ) values, revealing a more energetic surface. A Lewis basic character with K B /K A > 1 (but with a higher basicity for the neat monolith) was reported based on the results of the K B /K A ratios for both materials and suggests an acidic character for ZIF-8 [71].
In another study [70], the adsorption and separation of ethanol/ethyl acetate (EA)/water mixtures, formed during an industrial process for the synthesis of EA, were studied on ZIF-8. Single-and multicomponent isotherms of adsorption, acquired via gravimetric and breakthrough experiments, showed high adsorption ability and selectivity to EA. The higher interaction strength between EA and ZIF-8 was verified by IGC experiments. At zero coverage, the EA adsorption enthalpy was equivalent to 38.6 kJ mol −1 , while ethanol sorption released less energy (23.4 kJ mol −1 ). The larger polarizability of π-electrons containing ester carbonyl groups than hydroxyl groups may further enhance adsorption [119].
ZIF-68, which has a gme topology, consists of ZnN 4 tetrahedra connected by two derivatives of imidazole linkers; the first is a polar 2-nitroimidazolate linker, and the second is a nonpolar benzimidazolate linker, in a 1:1 ratio. ZIF-68 consists of two diverse channels and three types of cages, which are large, medium, and small [120]. A six-membered channel (diameter approximately 3.4-3.8 Å) forms alternating small and medium cages, where nitro groups are projected within the cages, maintaining the polar character. The larger channel (12-membered with a diameter of approximately 7.5-8.6 Å) is predominantly constructed from imidazole linkers substituted with benzene, resulting in this channel's pronounced nonpolar character. Since the large cages occupy 74% of the total volume, ZIF-68 was assigned an overall nonpolar behavior [121].
An IGC investigation of ZIF-68 [72] demonstrated a lack of shape-selective behavior with the complex pore structure of this material. The logarithm of the Henry constant was plotted as a function of the logarithm of vapor pressure for n-alkanes and isoalkanes to evaluate the shape selectivity ( Figure 14). The same linear correlation for linear and branched alkanes indicated the absence of shape selectivity [122,123]. The zero-coverage adsorption enthalpies of the branched alkanes were approximately 5−8 kJ mol −1 lower than those of the linear chains. This is typical of materials that are not shape selective. The structure demonstrated an overall nonpolar framework with limited polar centers. In contrast to alkanes, alkenes, and aromatics, small polar probes exhibited relatively high adsorption enthalpies. Due to the presence of nitro moieties, the cages in the small channel served as favorable adsorption positions for such polar molecules.

UiO-66 MOF
UiO-66 MOF (UiO = Universitetet i Oslo) is a zirconium-based MOF composed of hexamers of eight coordinated ZrO 6 (OH) 2 polyhedra and BDC linkers [124]. UiO-66 has an fcu topology, with a rigid cubic three-dimensional structure comprised of octahedral and tetrahedral cages with diameters of 11 and 8 Å, respectively. The cages are accessed solely via the triangular pore windows of approximately 5-7 Å (Figure 15). UiO-66 is considered one of the first MOFs that showed reverse shape-selectivity separation properties [125]. In a regular shape-selective material, a molecule having a more compact size with respect to its molecular weight will adsorb preferentially. Probes with higher boiling points and comparable chemical compositions will selectively adsorb on nonselective surfaces over lower-boiling molecules.
Chang et al. [73] demonstrated that a UiO-66-coated capillary column showed both a molecular sieving effect and a reverse shape selectivity. This capillary column favored adsorption of branched alkanes over their linear analogs, thus demonstrating reverse shape selectivity. This column also showed that branched hexane isomers had higher adsorption enthalpies than hexane, indicating stronger interaction with branched alkanes. However, the retention time of the bulkier 1,3,5-trimethylbenzene was smaller than that of n-propyl-and iso-propylbenzene, exhibiting a molecular sieving effect ( Figure 16). Table 2 summarizes the author's attempt to investigate the effect of pore-window aperture and size on adsorption and separation properties of alkane molecules by comparing several MOFs with UiO-66.  Bozbiyik et al. [74] investigated the abovementioned effect of UiO-66 on cycloalkanes. The adsorption and separation of n-hexane and cyclohexane on UiO-66 MOF was studied via IGC. Cyclohexane, which is bulkier and has a higher boiling point than n-hexane, was retained more than n-hexane. The adsorption enthalpy of cyclohexane (−60 kJ mol −1 ) was considerably greater than that of n-hexane (−53 kJ mol −1 ), despite the slightly lower molecular weight. Furthermore, the adsorption entropy of n-hexane was more negative, demonstrating that the bulkier cyclic molecule loses less freedom in the pores of UiO-66 than the linear n-hexane molecule.
In comparison with linear alkanes, the cyclic alkane has both entropic and enthalpic advantages. The composition of the UiO-66 pore system may explain this. In tetrahedral cages of UiO-66, cyclohexane can maintain a relatively high degree of rotational freedom. The enhanced confinement in these small cages results in a great deal of energetic interaction.
Seeking a deeper understanding of the reverse shape-selectivity performance of UiO-66, the same research group employed IGC at zero coverage for the separation of more than 70 molecules, including alkanes, alkenes, aromatics, and cyclic types on UiO-66 and its derived analogs, such as UiO-66-NO 2 , UiO-66-Me, and UiO-66-Me 2 [75].
High adsorption selectivity indicators were found for cycloalkanes in comparison to n-alkanes, such as the ratio of Henry constants (up to 14). Additional functional groups within the linkers in the framework, such as methyl and nitro groups, markedly affected the adsorption properties, primarily by decreasing the effective pore aperture size. Although improved selectivity for UiO-66-Me was observed, the selectivity decreased when an additional CH 3 group was present in UiO-66-Me 2 . This additional CH 3 group enabled confinement factors within the pores, and thus the material's adsorption properties, to be tuned. Selective adsorption was observed as a result of interactions in the smallest octahedral cages. However, intense confinement in the small tetrahedral cage led to stereoselective separation of cis/trans alkenes and disubstituted cycloalkanes. Monte Carlo simulations for pristine UiO-66 MOF and an analysis of the molecular positioning of adsorbates inside the framework supported the experimental results. In addition, the simulations illustrated contributions by entropic factors to the observed selectivity of adsorption.
The effect of chain length on the adsorption of branched and linear alkanes on UiO-66 and its functionalized analogs has also been explored [76]. Characteristic enhancement in Henry constants, adsorption enthalpy, and adsorption entropy occurred at particular carbon chain lengths of n-alkanes for hexane and heptane In case of 2-and 3-methyl alkanes and double branched alkanes, the effect of chain length was even more pronounced. This recognizable change in adsorption took place at a molecular size that matched the pore size of the smallest tetrahedral cages. The resultant selective adsorption emerged from the confinement properties and was a function of the molecule's size and shape.
The compensation plot ( Figure 17) shows the relationship between adsorption entropy and adsorption enthalpy, where the logarithm of pre-exponential factors is presented against the adsorption enthalpies. Generally, a linear connection was witnessed for a group of linear alkanes. The more strongly the probed molecule was adsorbed, the less freedom it had, which was identified as the compensation effect [2,127,128]. This proportional linear connection amongst linear alkanes was perceived for both pristine and SO 2 -functionalized DUT-5 MOFs. The slopes of the compensation lines, which showed an increase in adsorption entropy with respect to the adsorption enthalpy, were similar for both MOFs, which indicated that they had dissimilar adsorption enthalpies for identical molecules. Thus, the loss of freedom for a given interaction was identical in both MOFs. Thus, the shift in pre-exponential factors between the two MOFs resulted from their different adsorption sites and not from their different entropies of adsorption. Based on that, it was determined that the SO 2 functionalization offered better adsorption sites for n-alkane molecules that showed a 4-5 kJ mol −1 higher adsorption enthalpy compared to pristine DUT-5. SO 2 -DUT-5 demonstrated shape-selective activity in the adsorption of branched and linear alkanes at low coverage because of the reduction in pore size. In contrast, due its larger pore size, pristine DUT-5 was not shape selective. The ratio of Henry constants of linear and branched alkanes (i.e., separation factors) for all isoalkane/n-alkane pair separations on SO 2 -DUT-5 were larger than 1.7. DUT-5, in contrast, did not show a marked difference at all for the branched molecules [77].
Alkanes and alkenes adsorbed on SO 2 -DUT-5 to a similar degree. Although DUT-5 did not have a specific interaction with molecules possessing unsaturated bonds, SO 2 -DUT-5 demonstrated an improved interaction with aromatic molecules, although this interaction is not very strong and led to only moderate low-pressure separation. Breakthrough studies have shown that at higher pressures and loading, mixtures of aromatics with alkanes or alkenes can be separated, with both pristine and SO 2 -functionalized DUT-5 MOFs [77]. This distinction is a result of the molecular packing effects rather than specific interactions at a high level of pore filling. Such an effect has already been described several times with other MOFs, for example MIL-47 and MIL-53 [129,130], where separation occurred above a critical level of pore filling, in which the shift occurs from a single alignment of molecules in the one-dimensional pores to a double alignment of molecules.
In conclusion, SO 2 -DUT-5 showed a slight tendency to adsorb π-bond-containing molecules oversaturated molecules, as indicated by the marginally higher Henry constants for unsaturated compounds. Steric effects and specific adsorption played crucial roles and were correlated to the separation selectivity of alkanes/aromatics, also leading to shape selectivity in isoalkane adsorption. Van der Waals interactions dominated the adsorption properties of DUT-5. Consequently, aromatics, cyclo-alkanes, and alkenes exhibited lower values for both Henry constants and adsorption enthalpies than n-alkanes [77].
Due to its ultramicroporous structure with respect to its molecular dimensions, Cu(Qc) 2 [Cu(quinoline-5-carboxyate) 2 ] is a rather attractive MOF [131]. IGC studies were performed on Cu(Qc) 2 to determine its surface free energy of adsorption in order to evaluate the affinity between the O 2 /N 2 mixture [78].
Since O 2 and N 2 are nonpolar, their surface free energies of adsorption are basically reliant on dispersion forces. The chromatographic peaks illustrated that N 2 eluted before O 2 , which can be accredited to the higher affinity of Cu(Qc) 2 towards O 2 than N 2 , or to the high surface free energy or dispersive energy of O 2 in comparison to N 2 under the same conditions ( Figure 18). Using a molecular building block (MBB) approach [132,133], a series of RE-fcu-MOFs emerged, in which hexanuclear rare earth metal clusters served as twelve interconnected MBBs with face-centered cubic (fcu) topology and were connected by fluorinated/nonfluorinated ditopic linkers. The RE-fum-fcu-MOF (fum = fumaric acid) constructed from twelve coordinated RE hexanuclear clusters connected by fum showed a high affinity towards hydrocarbon separations [134].
Y-fum-fcu-MOF demonstrated the optimum window size for size-and/or shapeselective identification of hydrocarbon isomers with kinetic diameters of approximately 4.7 Å [79]. Assen et al. investigated the adsorption of C 4 alkenes with kinetic diameters between linear and branched C 4 alkanes. Different techniques were used, including IGC, to determine adsorption at infinite dilution and confirm the capacity of Y-fum-fcu-MOF for these separations. The retention time for butene isomer GC separation depended on two factors: the accessibility of adsorbates to the pores and the enthalpy of adsorption after the molecules accessed the pores. In spite of the comparable ∆H ads for cis and trans-2-butene, a result of their similar chemical natures, the molecule shape and size were the key factors in determining the adsorption selectivity [79].

Future Outlook
There are many applications of MOFs that depend on their adsorption properties. Hence, we have reviewed the utilization of IGC as a tool to investigate the physicochemical features of MOFs. IGC has helped to evaluate and predict the ability of MOF surfaces to contribute in both dispersive and specific interactions. Unlike an analytical GC technique, in IGC the roles of the mobile and stationary phases are reversed. This means that the material under study becomes the stationary phase and the injected well-known molecules within the mobile phase are probes, which are used to investigate the physicochemical properties of the stationary phase.
An excellent agreement has been revealed between the IGC-measured enthalpies of adsorption, and those estimated using other surface energy analysis techniques, such as calorimetric, volumetric, or gravimetric measurements, and even molecular modeling measurements. Moreover, IGC showed shorter experimental times compared to other techniques without almost any specific preparation of samples for crystalline or amorphous materials. The probability of observing the differences associated with surface group orientations is one of the benefits of using IGC in researching acid-and base-specific properties.
The total solid surface energy is represented by the summation of the dispersive and specific properties. From the thermodynamic properties of a material, its adsorption behavior can be interpreted and predicted. For instance, a plot of the Henry constants logarithm versus the vapor pressure logarithm could be used to investigate whether the material's ability to discriminate between branched and linear alkanes is related to a size-exclusion property. The difference between enthalpy of adsorption and its vaporization analogue is how they demonstrate the strength of adsorption toward the tested probe. Imperfections due to the degradation of MOFs could also be detected by a large value for the entropy of adsorption. Moreover, the activity of the MOF material under examination could be measured and compared with classical microporous materials via the dispersive component of the surface energy as a function of London interactions and surface heterogeneity.
Some considerations should be accounted for while using the IGC technique. A minimal amount of testing probes should be used in order to maintain a monolayer surface coverage. The peak symmetry is should also be carefully controlled, as the precise retention time is a basic source of the subsequent data. Optimum flowrate should be estimated via the Van Deemter curve, as a higher flowrate prevents the system's equilibrium, and a lower flowrate means broader peaks and longer experimental times. The surface area of the examined material must be measured accurately, as it becomes one of the basic properties utilized in further measurements. The molecular area of the polar probes most used in IGC is another issue to be considered as the literature shows a different value for those probes; however, a recent study [135] addressed that problem.
Past research on a variety of materials has proved that the IGC technique can be fruitfully implemented to characterize the surface and bulk properties of solids based on their diverse morphologies and shapes. However, IGC as a tool for MOF characterization is still in its infancy and much of this area remains to be discovered. IGC is adaptable, accurate, and simple, attributes that can cement its position as an essential characterization technique in MOF research. While some MOFs have been investigated using IGC and revealed great potential, many others are still waiting. In addition, many more parameters of MOFs may be evaluated, such as shape, morphology, topology index [136], and surface area measurements [137]. Funding: This project was funded by the National Plan for Science, Technology, and Innovation (MAARIFAH), King Abdul Aziz City for Science and Technology, Kingdom of Saudi Arabia, grant number 14-ADV2447-02.