Experimental Volumetric Hydrogen Uptake Determination at 77 K of Commercially Available Metal-Organic Framework Materials

: Storage is still limiting the implementation of hydrogen as an energy carrier to integrate the intermittent operation of renewable energy sources. Among different solutions to the currently used compressed or liquiﬁed hydrogen systems, physical adsorption at cryogenic temperature in porous materials is an attractive alternative due to its fast and reversible operation and the resulting reduction in storage pressure. The feasibility of cryoadsorption for hydrogen storage depends mainly on the performance of the used materials for the speciﬁc application, where metal-organic frameworks or MOFs are remarkable candidates. In this work, gravimetric and volumetric hydrogen uptakes at 77 K and up to 100 bar of commercially available MOFs were measured since these materials are made from relatively cheap and accessible building blocks. These materials also show relatively high porous properties and are currently near to large-scale production. The measuring device was calibrated at different room temperatures to calculate an average correction factor and standard deviation so that the correction deviation is included in the measurement error for better comparability with different measurements. The inﬂuence of measurement conditions was also studied, concluding that the available adsorbing area of material and the occupied volume of the sample are the most critical factors for a reproducible measurement, apart from the samples’ preparation before measurement. Finally, the actual volumetric storage density of the used powders was calculated by directly measuring their volume in the analysis cell, comparing that value with the maximum volumetric uptake considering the measured density of crystals. From this selection of commercial MOFs, the materials HKUST-1, PCN-250(Fe), MOF-177, and MOF-5 show true potential to fulﬁll a volumetric requirement of 40 g · L − 1 on a material basis for hydrogen storage systems without further packing of the powders. slope in the adsorption plateau increases. Microscopy images were taken with a Leica TL5000 optic microscope coupled to a Leica DMC5400 camera. MOF-177, ZIF-8, HKUST-1, MIL-100(Fe), PCN-205(Fe), CAU-10, UIO-66, UIO-67, and Zr-FA. The crystal phase of all the materials is stable during the hydrogen adsorption measurement, but special care was necessary for materials MOF-5 and UIO-67 to avoid contact with moisture during preparation. All the materials show permanent porosity during adsorption measurements, and their gravimetric uptakes of hydrogen at 77 K and high pressure depend on their porous properties BET area and pore volume. The measurement accuracy depends on the total available area of analyzed materials, with at least 1000 m 2 of sorbent area needed to reduce the correction error below 10%. However, the sample must ﬁt in the cold volume of the analysis cell. Material MOF-177 shows the highest gravimetric uptake at high pressure and higher deliverable capacity for isothermal operation. Materials MOF-5 and PCN-250(Fe) show the highest potential for volumetric hydrogen storage considering the measured density of the MOF crystals, over 50 g · L − 1 . However, the actual volumetric uptake considering the density of the powders is lower, mainly due to the presence of death volumes between MOF particles. The difference between calculated and measured density depends on the size and shape of crystals, the difference being smaller for powders made out of homogeneous and large crystals. The best crystals ´ self-packing is for material HKUST-1. The powders of materials MOF-177, PCN-250(Fe), and HKUST-1 store hydrogen with densities over 40 g · L − 1 without further packing. These results indicate that materials HKUST-1, PCN-250(Fe), and MOF-177 are the most promising commercially available MOFs for volumetric hydrogen storage without an intensive packing process of as-synthesized powders.


Introduction
Storage is still a limiting factor in implementing the large-scale use of hydrogen as an energy carrier to assist the world in the transition to using 100% green energy [1]. Despite the acceptance and use of compressed gas and the liquified state for hydrogen storage and transportation, more energetically efficient and operationally safe solutions will be requested. These solutions would preferably need milder temperature and pressure conditions, although still involving fast loading and release kinetics similar to the compressed gas [2]. Among the different suggested mechanisms as alternative solutions for hydrogen storage, physical adsorption or physisorption can be advantageous because of its almost complete reversibility, fast adsorption/desorption kinetics, relatively high storage density [3,4], and lower pressure or higher temperature than those involved during compression or liquefaction of hydrogen, respectively, for a similar storage density [1,5].
The uptake of hydrogen by physisorption was measured for a wide range of porous materials like zeolites [6,7], with a high affinity but a usually limited capacity because of their low porosity, and carbon-based materials [4,8], with a relatively high adsorption
PXRD patterns were collected with Cu-Kα radiation (λ = 1.50406 nm) on a D8 Advanced diffractometer (Bruker AXS, Germany) equipped with a Lynxeye detector. Samples were measured in reflection geometry in a 2Ɵ range from 3° to 50° with a step size of 0.009°. N2 adsorption/desorption isotherms were measured at 77 K on a volumetric ASAP 2020 device (Micromeritics) in the range of relative pressures p/p0 from 10 −7 to 0.996. For analysis, 0.2-0.3 g of each material were loaded in the analysis probe, followed by thermal activation (see Table 1 for activation temperature) under dynamic vacuum to reach a pressure value below 10 −6 mbar. The warm and cold free spaces were calculated by He expansion on the activated sample after the N2 experiment. The BET area (ABET) was calculated by following Rouquerol and Llewellyns' rule to select the optimal BET range [26]. The total pore volume (Vp) was calculated from the N2 uptake at p/p0 ≈ 0.99, and the micropore volume (Vµp) from the corresponding N2 uptake at p/p0 value immediately before the mesoporous or the interparticle condensation, when the slope in the adsorption plateau increases. Microscopy images were taken with a Leica TL5000 optic microscope coupled to a Leica DMC5400 camera.
PXRD patterns were collected with Cu-K α radiation (λ = 1.50406 nm) on a D8 Advanced diffractometer (Bruker AXS, Germany) equipped with a Lynxeye detector. Samples were measured in reflection geometry in a 2Θ range from 3 • to 50 • with a step size of 0.009 • . N 2 adsorption/desorption isotherms were measured at 77 K on a volumetric ASAP 2020 device (Micromeritics) in the range of relative pressures p/p 0 from 10 −7 to 0.996. For analysis, 0.2-0.3 g of each material were loaded in the analysis probe, followed by thermal activation (see Table 1 for activation temperature) under dynamic vacuum to reach a pressure value below 10 −6 mbar. The warm and cold free spaces were calculated by He expansion on the activated sample after the N 2 experiment. The BET area (A BET ) was calculated by following Rouquerol and Llewellyns' rule to select the optimal BET range [26]. The total pore volume (V p ) was calculated from the N 2 uptake at p/p 0 ≈ 0.99, and the micropore volume (V µp ) from the corresponding N 2 uptake at p/p 0 value immediately before the mesoporous or the interparticle condensation, when the slope in the adsorption plateau increases. Microscopy images were taken with a Leica TL5000 optic microscope coupled to a Leica DMC5400 camera. Hydrogen adsorption experiments were performed in a HPVA II (high-pressure volumetric analyzer) from Micromeritics. For a typical experiment, 0.3-0.9 g of material were loaded in the sample tube. Activation was performed in the HPVA in the same conditions that those used for N 2 adsorption/desorption experiments. Densities for He and hydrogen gases were calculated from p and T values by MBWR EOS of McCarty and Arp and Helmholtz EOS, respectively, as recommended by NIST [27,28]. The controlled room temperature is stable within ±1 • C, and the manifold temperature is kept at 32-33 • C. The device uses two pressure transducers for reading manifold's pressure up to 1000 mbar (accuracy 0.15% of the reading) and up to 200 bar (accuracy 0.04% of the full-scale). The figure of merit is always more than 100% during a measurement compared to the accuracy limit [29]. The experiments were performed without a liquid nitrogen level control, and the sample tube was covered by a cryojacket which kept the cryogenic liquid level approximately constant during the experiment. Adsorption equilibration was assured by using 20 min of equilibration time in all the adsorption steps. The reported excess uptake was averaged from different measurements with an intermediate evacuation for at least six hours at room temperature. All the isotherms are plotted from Figures S5-S15 in Supplementary Materials. More details about raw data correction and calibration are provided in Appendix A. The apparent or skeleton volume of materials was measured in the HPVA by subtracting the free space of the sample tube occupied with materials to the free space of the empty tube, averaged from 100 He-expansion measurements at different room temperatures. Finally, the skeleton density (ρ sk ) was calculated by dividing the measured mass of outgassed material by the measured sample's volume. The crystal density of each MOF (ρ MOF ) was calculated from Equation (1) from ρ sk and the pore volume (corresponding to micropores or total pore volume). The bed density (ρ bed ) was calculated from the dry mass of material and the measured sample's volume using the shown tool in Figure S17 in Supplementary Materials and Equation (2) [30]. The total gravimetric uptake of materials was calculated from the measured excess uptake, the pore volume, and the gas density at the specific pressure and temperature by Equation (3). The volumetric uptake of MOFs was calculated by multiplying the gravimetric total uptake times the corresponding density value.

Results
The crystalline phase of all materials is stable after their preparation and hydrogen adsorption experiments (Figure 2 and Figure S1 in Supplementary Materials), except material UIO-67, which loses the crystalline phase after hydrogen adsorption experiments despite its previously reported thermal and chemical stability [31]. However, a slight reduction in the pore volume after methanol and ethanol adsorption experiments previously indicated some instability issues [32]. Since the sample's volume contracted in the sample holder during waiting for the PXRD analysis, material UIO-67 probably lost the crystal phase between the adsorption analysis and the PXRD measurement due to contact with ambient, similar to the behavior of other MOF compounds [33]. An evolving diffraction signal at 8.94 • after hydrogen adsorption in material MOF-5 is not present in the simulated pattern, indicating structural changes because of the known material's instability to moisture [24]. Material MOF-177 is also unstable to ambient moisture, but there is no difference in the PXRD pattern after hydrogen adsorption experiments and X-ray measurements, indicating higher stability to ambient moisture than MOF-5. The PXRD pattern for material MIL-100(Fe) starts at a lower diffraction angle of 2θ because this structure's most relevant diffraction signals appear at 2-3 • as a consequence of its mesoporous cavities.

Results
The crystalline phase of all materials is stable after their preparation and hydrogen adsorption experiments (Figures 2 and S1 in Supplementary Materials), except material UIO-67, which loses the crystalline phase after hydrogen adsorption experiments despite its previously reported thermal and chemical stability [31]. However, a slight reduction in the pore volume after methanol and ethanol adsorption experiments previously indicated some instability issues [32]. Since the sample's volume contracted in the sample holder during waiting for the PXRD analysis, material UIO-67 probably lost the crystal phase between the adsorption analysis and the PXRD measurement due to contact with ambient, similar to the behavior of other MOF compounds [33]. An evolving diffraction signal at 8.94° after hydrogen adsorption in material MOF-5 is not present in the simulated pattern, indicating structural changes because of the known material's instability to moisture [24]. Material MOF-177 is also unstable to ambient moisture, but there is no difference in the PXRD pattern after hydrogen adsorption experiments and X-ray measurements, indicating higher stability to ambient moisture than MOF-5. The PXRD pattern for material MIL-100(Fe) starts at a lower diffraction angle of 2θ because this structure's most relevant diffraction signals appear at 2°-3° as a consequence of its mesoporous cavities. Studied materials show type I N2 adsorption/desorption isotherms at 77 K, according to the IUPAC classification, corresponding to microporous materials ( Figure 3) [26,34]. MOF-177 shows the highest N2 uptake and a calculated BET area of 4449 m²·g −1 , with a pore volume of 1.80 cm³·g −1 . MOF-5 shows a lower surface area than expected from previous measurements, around 3500-3600 m²·g −1 [16]. A N2 adsorption/desorption experiment of this material performed in the hydrogen analyzer gave a value ca. 3100 m²·g −1 of BET area, indicating that some structural damage could be produced during the transference of material from hydrogen to nitrogen experiments, as PXRD suggested. The BET area of material MIL-100(Fe) is ca. 200 m²·g −1 higher than that for PCN-250 despite being Studied materials show type I N 2 adsorption/desorption isotherms at 77 K, according to the IUPAC classification, corresponding to microporous materials ( Figure 3) [26,34]. MOF-177 shows the highest N 2 uptake and a calculated BET area of 4449 m 2 ·g −1 , with a pore volume of 1.80 cm 3 ·g −1 . MOF-5 shows a lower surface area than expected from previous measurements, around 3500-3600 m 2 ·g −1 [16]. A N 2 adsorption/desorption experiment of this material performed in the hydrogen analyzer gave a value ca. 3100 m 2 ·g −1 of BET area, indicating that some structural damage could be produced during the transference of material from hydrogen to nitrogen experiments, as PXRD suggested. The BET area of material MIL-100(Fe) is ca. 200 m 2 ·g −1 higher than that for PCN-250 despite being made out of the same trinuclear iron oxide SBU. The first uses tritopic trimesic acid as linker, while the latest uses tetratopic 3,3',5,5'-azobenzenetetracarboxylic acid yielding a structure with lower pore volume and higher crystal density. Besides, material MIL-100(Fe) reaches the adsorption plateau at higher p/p 0 as a consequence of its larger cavities size in the range of mesopores (24-29 Å) [35,36]. Materials UIO-66, MOF-5, MIL-100(Fe), ZIF-8, and, slightly, MOF-177 show an increase in the adsorption uptake at high p/p 0 values, indicating interparticle condensation [26]. The more important this contribution is, the more significant is the difference between total pore volume and micropore volume (Table 1). For UIO-66, this adsorption step starts at lower p/p 0 than in the other MOFs, near 0.9, resulting in a type IV isotherm corresponding to mesoporous materials despite the microporous structure of UIO-66 [37]. This isotherm also shows an H3 hysteresis loop between sorption and desorption, indicating slit-like mesopores, probably between UIO-66 crystals, similar to previous reported UIO-66 samples [38]. Material UIO-66 was not distributed as a powder but as big blocks which needed to be ground before characterization (see fragments in Figure S2 in Supplementary Materials). These slit-like mesopores might be space between the intergrown microporous crystals of the MOF. The framework of UIO-67 is isorreticular to UIO-66 but uses a larger linker, 4,4'-biphenyldicarboxylic acid, instead of terephthalic acid, which expands the pore and increases both pore volume and surface area. The calculated BET area of this material is similar to that previously measured [31] and that calculated by simulation [39] since there was no contact with ambient moisture during the material's pre-treatment to the adsorption analysis. Materials Zr-FA (UIO-66FA) and CAU-10 use fumaric acid as linker, showing the lowest surface areas. made out of the same trinuclear iron oxide SBU. The first uses tritopic trimesic acid as linker, while the latest uses tetratopic 3,3',5,5'-azobenzenetetracarboxylic acid yielding a structure with lower pore volume and higher crystal density. Besides, material MIL-100(Fe) reaches the adsorption plateau at higher p/p0 as a consequence of its larger cavities size in the range of mesopores (24-29 Å) [35,36]. Materials UIO-66, MOF-5, MIL-100(Fe), ZIF-8, and, slightly, MOF-177 show an increase in the adsorption uptake at high p/p0 values, indicating interparticle condensation [26]. The more important this contribution is, the more significant is the difference between total pore volume and micropore volume ( Table 1). For UIO-66, this adsorption step starts at lower p/p0 than in the other MOFs, near 0.9, resulting in a type IV isotherm corresponding to mesoporous materials despite the microporous structure of UIO-66 [37]. This isotherm also shows an H3 hysteresis loop between sorption and desorption, indicating slit-like mesopores, probably between UIO-66 crystals, similar to previous reported UIO-66 samples [38]. Material UIO-66 was not distributed as a powder but as big blocks which needed to be ground before characterization (see fragments in Figure S2 in Supplementary Materials). These slit-like mesopores might be space between the intergrown microporous crystals of the MOF. The framework of UIO-67 is isorreticular to UIO-66 but uses a larger linker, 4,4'-biphenyldicarboxylic acid, instead of terephthalic acid, which expands the pore and increases both pore volume and surface area. The calculated BET area of this material is similar to that previously measured [31] and that calculated by simulation [39] since there was no contact with ambient moisture during the material's pre-treatment to the adsorption analysis. Materials Zr-FA (UIO-66FA) and CAU-10 use fumaric acid as linker, showing the lowest surface areas.   [39]. Measured and calculated crystal densities are reasonably similar considering the measurement limitations in the skeleton density and the pore volume. The most significant differences appear for materials UIO-66, Zr-FA, and CAU-10, for which the measured density values are   [39]. Measured and calculated crystal densities are reasonably similar considering the measurement limitations in the skeleton density and the pore volume. The most significant differences appear for materials UIO-66, Zr-FA, and CAU-10, for which the measured density values are lower than the theoretical ones. A lower density could indicate an overestimation of the pore volume or even a low efficiency in removing solvent molecules during outgassing. For material CAU-10, the difference between measured and theoretical density is big enough to consider this last as overestimated. Indeed, the measured BET area herein is almost three times higher than the collected value in the CoRE MOF database. The deviation for material MIL-100(Fe) could be due to a wrong calculation of the micropore volume since it is a mesoporous structure. Indeed, the theoretical density value is among calculated crystal and bed densities.  (1) from ρ sk and V p . 3 Calculated powder density by Equation (1) from ρ sk and V µp . 4 Calculated crystal density from crystallographic information. 5 Calculated bed density from the measured volume of powders.
Unlike nitrogen adsorption, there are neither standard procedures nor reference materials to set a comparison basis for hydrogen cryoadsorption measurements, and it is possible to find irreproducible reported hydrogen uptakes in the literature [40,41] unless special considerations and good experimental practices are adequately observed. Most concerning aspects to reproducibility are related to the preparation of materials before adsorption, technical specifications of analyzers, and, in most cases, calibration procedures to correct recorded data [21]. For this reason, a calibration was carried out as described in the experimental section considering multiple analysis conditions, including room temperature and analysis duration, and the carbon material NORIT RAW 0.8 mm, which took part in an interlaboratory analysis for hydrogen uptake at 77 K [42], was used to ensure comparability in the hydrogen uptake measurements ( Figure S3 in Supplementary Materials).
Among the analyzed commercial MOFs, material MOF-177 shows the highest hydrogen excess uptake at 77 K due to its highest surface area (see Figure 4 and Table 1). All the isotherms show a maximum in the excess uptake at a pressure value between 20 bar and 50 bar because the adsorbed phase behaves as an incompressible fluid after its superficial saturation, reducing the excess uptake for higher pressure where the gas-phase density continues increasing [12]. Generally, the maximum gravimetric excess of materials at 77 K correlates with their BET area as described by Chahine's rule (ca. 1 wt. % of excess hydrogen each 500 m 2 /g of BET area) [4]. There are minor deviations to this trend due to differences in affinities to the sorbed gas, but they could also be owed to differences in the activation degree or even structure degradation before N 2 or hydrogen adsorption experiments, as seems to happen in material MOF-5. Materials MOF-5 and PCN-250 (Fe) have similar surface area values but a hydrogen uptake difference of 0.5 wt. %. The necessary pressure to reach the adsorption maximum is lower for the Fe-based material since it has coordinatively unsaturated metallic sites or open metal sites (OMS), which interact more strongly with hydrogen molecules [43]. The uptake of this material is 1 wt. % higher than that predicted by Chahine's rule, probably due to these OMS oriented to cavities [44]. Similar behavior is observed in OMS-rich material HKUST-1, with similar uptake than UIO-67 despite its lower surface area. The excess uptake of this last material is higher than that in isorreticular UIO-66 and Zr-FA materials, made out of shorter linkers and smaller cavities. Material ZIF-8, without OMS, adsorbs a higher amount of hydrogen than MIL-100(Fe) with a higher surface area and OMS [45]. The larger cavities size of structure MIL-100 could be the reason for this lower uptake, driving to a weaker interaction between gas molecules and adsorption surface due to the lower confinement effect of the gas in pores [11]. For this material, the more-separated OMS in the structure might dilute the effect of these preferential adsorption sites in increasing the superficial packing density of the adsorbed phase [44,46]. Finally, ZIF-8 and UIO-67 materials show hysteresis in desorption, which could be related, in the case of ZIF-8 material, to its previously-described gate-opening effect [47,48]. For material UIO-67, this could also indicate some flexibility in the network during desorption. The total uptake of gas adds the compressed-gas phase in the pore volume of materials to the superficial excess capacity (Equation (3)). For this reason, material MOF-177, first, and MOF-5, later, show the highest total stored amounts on a gravimetric basis ( Figure S16 in Supplementary Materials) since they have the highest excess uptakes and pore volumes. The effect of a larger pore volume is evident when comparing material MIL-100(Fe) with PCN-250(Fe), UIO-67, and HKUST-1: despite the lower excess uptake of the former, the total uptakes are similar due to the larger amount of compressed gas in pores out of the excess phase. At low pressure, materials HKUST-1 and PCN-250(Fe) show higher gravimetric uptake than MOF-177, indicating that the deliverable capacity of the former material will be higher due to a smaller amount of sorbed gas at low pressure.
that predicted by Chahine's rule, probably due to these OMS oriented to cavities [44]. Similar behavior is observed in OMS-rich material HKUST-1, with similar uptake than UIO-67 despite its lower surface area. The excess uptake of this last material is higher than that in isorreticular UIO-66 and Zr-FA materials, made out of shorter linkers and smaller cavities. Material ZIF-8, without OMS, adsorbs a higher amount of hydrogen than MIL-100(Fe) with a higher surface area and OMS [45]. The larger cavities size of structure MIL-100 could be the reason for this lower uptake, driving to a weaker interaction between gas molecules and adsorption surface due to the lower confinement effect of the gas in pores [11]. For this material, the more-separated OMS in the structure might dilute the effect of these preferential adsorption sites in increasing the superficial packing density of the adsorbed phase [44,46]. Finally, ZIF-8 and UIO-67 materials show hysteresis in desorption, which could be related, in the case of ZIF-8 material, to its previously-described gateopening effect [47,48]. For material UIO-67, this could also indicate some flexibility in the network during desorption. The total uptake of gas adds the compressed-gas phase in the pore volume of materials to the superficial excess capacity (Equation (3)). For this reason, material MOF-177, first, and MOF-5, later, show the highest total stored amounts on a gravimetric basis ( Figure S16 in Supplementary Materials) since they have the highest excess uptakes and pore volumes. The effect of a larger pore volume is evident when comparing material MIL-100(Fe) with PCN-250(Fe), UIO-67, and HKUST-1: despite the lower excess uptake of the former, the total uptakes are similar due to the larger amount of compressed gas in pores out of the excess phase. At low pressure, materials HKUST-1 and PCN-250(Fe) show higher gravimetric uptake than MOF-177, indicating that the deliverable capacity of the former material will be higher due to a smaller amount of sorbed gas at low pressure. The precision of measurements is higher for materials MOF-177, MOF-5, and HKUST-1 compared to materials CAU-10, ZIF-8, and UIO-67 (see Figures S5-S15 in Supplementary Materials). In addition to the reproducibility of the adsorption itself, because of differences in the analysis temperature, the material's activation degree between experiments, or even degradation of the adsorbent [21], the observed deviation herein is more related to the correction procedure of the raw data (see Figure S4 in Supplementary Materials). The ambient factors affecting the blank analyses are also affecting the sorption experiments. For this reason, the used correction factor was averaged for several blank analyses during a long operating time; therefore, it considers the effect of the possible room-temperature's variation during or between experiments. Figure 5 represents the difference at 100 bar between corrected data with the highest and the lowest correction factors versus the available area of the analyzed material (A BET ·m mat for each sample). The larger the available area, the smaller the correction error because of the more significant contribution of the adsorption itself. Instead of an average correction factor, a specific blank analysis reproducing the same analysis conditions and room temperature would be ideal, slowing the data acquisition down and reducing the experimental flexibility. In this experimental setup, ca. 1000 m 2 of total sorbent area are necessary to reduce the correction error below 10%. Nevertheless, a larger surface area available for analysis does not always lead to greater accuracy. The hydrogen uptake of material MOF-177 is lower when using a 3.6 times higher amount of adsorbent (see Figure S15 in Supplementary Materials) because the occupied volume of the sample is 2.91 cm 3 , which is larger than the effective cold analysis volume (2.46 ± 0.05 cm 3 ). In this case, the fraction of material out of the temperature-controlled volume reduces the overall sample's uptake. However, ZIF-8's powder occupies 1.95 cm 3 and adsorbs the same amount of hydrogen after a 36% packing of its volume (see Figures S9 and S10 in Supplementary Materials). Therefore, the sample's volume must always be contained where the analysis temperature is known and controlled during the experiment.
experiments. For this reason, the used correction factor was averaged for several blank analyses during a long operating time; therefore, it considers the effect of the possible room-temperature's variation during or between experiments. Figure 5 represents the difference at 100 bar between corrected data with the highest and the lowest correction factors versus the available area of the analyzed material (ABET·mmat for each sample). The larger the available area, the smaller the correction error because of the more significant contribution of the adsorption itself. Instead of an average correction factor, a specific blank analysis reproducing the same analysis conditions and room temperature would be ideal, slowing the data acquisition down and reducing the experimental flexibility. In this experimental setup, ca. 1000 m² of total sorbent area are necessary to reduce the correction error below 10%. Nevertheless, a larger surface area available for analysis does not always lead to greater accuracy. The hydrogen uptake of material MOF-177 is lower when using a 3.6 times higher amount of adsorbent (see Figure S15 in Supplementary Materials) because the occupied volume of the sample is 2.91 cm³, which is larger than the effective cold analysis volume (2.46 ± 0.05 cm³). In this case, the fraction of material out of the temperature-controlled volume reduces the overall sample's uptake. However, ZIF-8's powder occupies 1.95 cm³ and adsorbs the same amount of hydrogen after a 36% packing of its volume (see Figures S9 and S10 in Supplementary Materials). Therefore, the sample's volume must always be contained where the analysis temperature is known and controlled during the experiment. The volumetric uptake is not only proportional to the porous properties of materials but also the adsorbent density [49,50], and it represents the adsorbate packing efficiency in the occupied space by the material [51,52]. All the analyzed materials show higher volumetric storage density than the compressed gas at 77 K and the same pressure ( Figure 6 and Figure S18 in Supplementary Materials) except material CAU-10. The volumetric gain of the latest is negative from 70 bar regardless of the used MOF density value. From this selection, only materials MOF-5 and PCN-250(Fe) show potential to store 50 g·L −1 , the technical target from US DOE for usable hydrogen for light-duty vehicles [53], when their uptakes are calculated from the measured crystal density. The usable amount is smaller than the uptake, but these two amounts are nearer using a temperature-pressure swing adsorption system 100 bar 77 K-5 bar 160 K [54]. In comparison, the Fe-based material needs a lower pressure to reach this storage density because of the higher excess uptake at low pressure. These data indicate that the deliverable capacity of this material will be lower unless a higher desorption temperature is used to release the same amount of the stored hydrogen than in MOF-5. It will also mean higher liquid nitrogen consumption during the adsorption process to remove the resulting heat, a higher heat consumption (or lower thermal insulation) during the hydrogen release, and lower liquid nitrogen usage during long-term storage. Materials MOF-177, UIO-67, and HKUST-1 could potentially fulfill the target of 40 g·L −1 for rechargeable portable power equipment for 2025 [55]. The  (1) is relevant for calculating this storage density for materials MIL-100(Fe), MOF-5, ZIF-8, and, more importantly, UIO-66. For materials HKUST-1, CAU-10, MOF-177, and PCN-250, the calculated crystal and bed densities are similar because of the minor contribution of interparticle condensation in the N 2 -adsorption measurements to the total pore volume (see Figure 3).  The measured bed density is used to calculate the actual hydrogen uptake in the occupied volume by the powder (Figure 6). The actual uptake is always lower than that calculated by Equation (1) and needs a direct volume s measurement for a specific powder to include micro/mesopore volumes, interparticle space, or even further dead volumes in the storage system [52]. Only the bed storage density of HKUST-1 is close to the maximum value, over 40 g·L −1 , indicating a better self-packing degree, probably owed to relatively big polyhedral crystals with a regular size. Also, materials MOF-177 and PCN-250(Fe) are distributed as relatively big crystals ( Figure S2 from Supplementary Materials), reaching higher volumetric uptake than 40 g·L −1 . Crystals of PCN-250(Fe) are also polyhedral, but their sizes are not as homogeneous as in the HKUST-1 sample, which may increase void spaces between crystals [18,19]. Crystals of CAU-10 are tiny, with irregular shapes but similar in size. The significant difference between measured and calculated bed density can be attributed to the higher concentration of morphology defects in the tiny crystals' external surfaces, creating more interparticle space [19]. Crystal morphology is not distinguished in MOF-5, UIO-67, ZIF-8, and Zr-FA materials, and ultrafine powders are observed instead. For this reason, their volumetric uptakes could probably be enhanced by packing since powders are easier to pack than big crystals [18,19]. Materials MOF-177 and MOF-5 show similar packing density values than previously reported ones, and it is in principle possible to further compact them up to 0.34 and 0.28 g/cm 3 , respectively, without reducing their porous properties and excess gravimetric hydrogen uptake [30,56], finally storing up to 41 g·L −1 for material MOF-5 and 43 g·L −1 for MOF-177.

Conclusions
The gravimetric and volumetric hydrogen uptake at the cryogenic temperature of 77 K was measured in commercially available MOF materials MOF-5, MOF-177, ZIF-8, HKUST-1, MIL-100(Fe), PCN-205(Fe), CAU-10, UIO-66, UIO-67, and Zr-FA. The crystal phase of all the materials is stable during the hydrogen adsorption measurement, but special care was necessary for materials MOF-5 and UIO-67 to avoid contact with moisture during preparation. All the materials show permanent porosity during adsorption measurements, and their gravimetric uptakes of hydrogen at 77 K and high pressure depend on their porous properties BET area and pore volume. The measurement accuracy depends on the total available area of analyzed materials, with at least 1000 m 2 of sorbent area needed to reduce the correction error below 10%. However, the sample must fit in the cold volume of the analysis cell. Material MOF-177 shows the highest gravimetric uptake at high pressure and higher deliverable capacity for isothermal operation. Materials MOF-5 and PCN-250(Fe) show the highest potential for volumetric hydrogen storage considering the measured density of the MOF crystals, over 50 g·L −1 . However, the actual volumetric uptake considering the density of the powders is lower, mainly due to the presence of death volumes between MOF particles. The difference between calculated and measured density depends on the size and shape of crystals, the difference being smaller for powders made out of homogeneous and large crystals. The best crystals self-packing is for material HKUST-1. The powders of materials MOF-177, PCN-250(Fe), and HKUST-1 store hydrogen with densities over 40 g·L −1 without further packing. These results indicate that materials HKUST-1, PCN-250(Fe), and MOF-177 are the most promising commercially available MOFs for volumetric hydrogen storage without an intensive packing process of as-synthesized powders.

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