ZrBDC-Based Functional Adsorbents for Small-Scale Methane Storage Systems

Metal-organic frameworks (MOF), potentially porous coordination structures, are envisioned for adsorption-based natural gas (ANG) storage, including mobile applications. The factors affecting the performance of the ANG system with a zirconium-based MOF with benzene dicarboxylic acid as a linker (ZrBDC) as an adsorbent were considered: textural properties of the adsorbent and thermal effect arising upon adsorption. The high-density ZrBDC-based pellets were prepared by mechanical compaction of the as-synthesized MOF powder at different pressures from 30 to 240 MPa at 298 K without a binder and mixed with polymer binders: polyvinyl alcohol (PVA) and carboxyl methylcellulose (CMC). The structural investigations revealed that the compaction of ZrBDC with PVA under 30 MPa was optimal to produce the ZrBDC-PVA adsorbent with more than a twofold increase in the packing density and the lowest degradation of the porous structure. The specific total and deliverable volumetric methane storage capacities of the ZrBDC-based adsorbents were evaluated from the experimental data on methane adsorption measured up to 10 MPa and within a temperature range from 253 to 333 K. It was measured experimentally that at 253 K, an 100 mL adsorption tank loaded with the ZrBDC-PVA pellets exhibited the deliverable methane storage capacity of 172 m3(NTP)/m3 when the pressure dropped from 10 to 0.1 MPa. The methane adsorption data for the ZrBDC powder and ZrBDC-PVA pellets were used to calculate the important thermodynamic characteristic of the ZrBDC/CH4 adsorption system—the differential molar isosteric heat of adsorption, which was used to evaluate the state thermodynamic functions: entropy, enthalpy, and heat capacity. The initial heats of methane adsorption in powdered ZrBDC evaluated from the experimental adsorption isosteres were found to be ~19.3 kJ/mol, and then these values in the ZrBDC/CH4 system decreased at different rates during adsorption. In contrast, the heat of methane adsorption onto the ZrBDC-PVA pellets increased from 19.4 kJ/mol to a maximum with a magnitude, width, and position depended on temperature, and then it fell. The behaviors of the thermodynamic state functions of the ZrBDC/CH4 adsorption system were interpreted as a variation in the state of adsorbed molecules determined by a ratio of CH4-CH4 and CH4-ZrBDC interactions. The heat of adsorption was used to calculate the temperature changes of the ANG systems loaded with ZrBDC powder and ZrBDC pellets during methane adsorption under adiabatic conditions; the maximum integrated heat of adsorption was found at 273 K. The maximum temperature changes of the ANG system with the ZrBDC materials during the adsorption (charging) process did not exceed 14 K that are much lower than those reported for the systems loaded with activated carbons. The results obtained are of direct relevance for designing the adsorption-based methane storage systems for the automotive industry, developing new gas-power robotics systems and uncrewed aerial vehicles.


Introduction
A wealth of modern studies is devoted to finding solutions to reduce carbon dioxide emissions from power plants, which increase the amount of greenhouse gases in the atmosphere and contribute to global climate change [1].One of the ways to solve this problem is to replace liquid hydrocarbon fuels, which are refined petroleum products, with natural gas (NG), which contains more than 70 percent methane and produces the lowest carbon dioxide emissions.There are proven reserves of 47,805 billion cubic meters of NG in Russia [2].Therefore, there is considerable scope for expanding NG use, including energy systems of heating and electricity generating plants and onboard vehicular applications (motor vehicles and aircraft).However, at normal temperature and pressure (NTP), NG has a lower energy density than that of gasoline and diesel, and this is a major setback to bringing NG into widespread use as a fuel [3][4][5].At the same time, in many respects, the available technologies of compressed (CNG) and liquified natural gas (LNG) fail to meet efficiency and safety requirements [6][7][8].Therefore, new technological solutions are required to increase the energy density of NG, and an adsorbed NG (ANG) method seems to be the most effective one for accumulation and storage [9][10][11].Based on a rich set of experimental data reported over the past twenty years [11][12][13][14][15][16][17][18], activated carbons and metal-organic framework structures are recognized as the most promising adsorbents for employing in ANG systems.
For example, the textural properties of the most efficient activated carbons as MAXSORB-III [19] can provide the specific methane storage capacity of over 200 m 3 (NTP)/m 3 at the pressure of 3.4 MPa and 293 K.However, a further increase in the efficiency of carbon adsorbents is limited mainly by the standard methods of their synthesis, including carbonization and activation.It is known that an increase in the activation degree of a carbon material induces an increase both in the volume of micropores and the specific surface area [20].However, it inevitably leads to the expansion of pores and reduces the packing density [21], which negatively affects its specific total methane capacity.
Meanwhile, the synthesis of metal-organic frameworks (MOF) proceeds according to a fundamentally different mechanism, and the final sorption behaviors depend mainly on the structural parameters of a porous cluster formed by an organic ligand-metal ion pair, which makes it possible to form a variety of MOF structures with various properties.Wilmer et al. [22] and later Simon et al. [23] performed the screening analysis of more than a hundred thousand hypothetical MOF compounds in order to establish the relationship between the structural parameters and adsorption properties for methane.They found an optimal MOF structure that ensures the DOE 2012 target value of total methane adsorption capacity for a monolith crystal of 263-266 m 3 (STP)•m -3 at 3.5 MPa and 298 K [15].
Many studies of note reported on the methane adsorption in MOF, which attained or even went ahead of the DOE target achievement by changing the thermodynamic conditions: either by increasing pressure [24] or decreasing temperature [25,26].Moreover, there is a discrepancy in the theoretical calculations and experimental data caused by incorrect use of the crystal density of as-prepared powdered MOF in calculating a specific volumetric methane adsorption capacity [27].Indeed, most MOF compounds are synthesized as crystals ranging in size from tens to hundreds of microns, and the resultant sample is a polydisperse microcrystalline powder.However, the use of MOF powders as adsorbents does not allow achieving the high productivity of the ANG system due to their unsatisfactory performance: low values of packing density and mechanical resistance, insufficient rates of mass and heat transfer, adsorption cyclical instability, formation of dust.The low packing density does not provide the high adsorption capacity of an ANG tank (or adsorber) loaded with the MOF powder.Therefore, for the practical application of the MOF-based materials in the ANG systems, it is necessary to develop methods to increase their packing density via shaping [3,[28][29][30].The advantages of the shaped MOF materials are that they display better operational features such as high packing density and mechanical strength, humidity resistance, thermal conductivity, and resistance to physical stress, compared to the initial powders [30][31][32].
One of the approaches to MOF shaping implies in situ growth of MOF films on a porous substrate via the layer by layer method [33], atomic layer deposition [34], or their combination [35].Thakkar et al. [36] reported on the successful application of the 3D printing technology to fabricate the MOF-74(Ni) and UTSA-16(Co) monoliths with MOF contents of up to 80 and 85 wt.%, respectively.Posttreatment strategies employing the mixing of presynthesized MOFs with additives followed by microwave-assisted polymerization [37] or freeze-drying [38,39] are proposed to fabricate shaped MOF-based composites with different MOFloadings and densities.Tian with coworkers enabled to prepare a monolithic porous Cu-based MOF (HKUST-1) with an increased density by sol-gel process [40].The monolithic MOF showed robust mechanical properties and increased methane uptake.
Another approach implies mechanical shaping of MOF powders: dry or wet granulation [41], spray-drying [42], and pressing or extrusion under high pressures [43].Mechanical pressing (compaction) of MOF powder is the simplest, flexible, and affordable technique used to prepare shaped MOF materials with higher density, thermal conductivity, and improved resistance to abrasion [30].Two options for preparing MOF-based pellets or granules are possible.First, a binder is used to fabricate MOF pellets, which facilitates the shaping process and improves the mechanical strength and wear resistance of MOF-based material [44,45].Another approach is to treat an assynthesized MOF powder without any binders [46,47].However, the compaction of MOF powders, especially under high pressures, the granulation, or extrusion shaping can cause irreversible degradation of the porous structure of MOF, and as a result, deterioration of the sorption performance [46].Nandasiri with coworkers reported [48] that the shaping of different MOFs such as MIL-101 (Cr), MIL-53 (Al), MOF-5, and HKUST-1 reduced the specific surface area, S BET , by 50% and the volume of pores by ~43%.

2
Adsorption Science & Technology the decomposition occurs at temperatures above 450 °С, and it is resistant to boiling water and can withstand low concentrations of H 2 S, SO x , and NO x [52].It was found that the shaping of the ZrBDC powder under 665 MPa did not significantly affect the structure; its specific surface and micropore volume decreased by no more than 2% [53].
It is worthwhile noting that the major impediments to large-scale industrial use of MOFs are the complexity of the synthesis process, the high cost of reagents, and, accordingly, the final product [54].In this regard, it is advisable to link the prospects of using MOFs with their application as an adsorbent to store alternative fuels in small-scale energy systems as uncrewed air vehicles [55,56].Optimal shaping of MOF powders into a tailorable form is a significant step towards the implementation of this engineering solution.
In the previous studies [28,57,58], we tested the conditions for preparing high-density functional materials from Cu-based MOF (CuBTC, where BTC is 1,3,5-benzene tricarboxylic acid) and ZrBDC.The primary purpose of the studies was to develop an optimal technique for fabricating the shaped materials from CuBTC and ZrBDC with reproducible properties.In particular, the shaping of the CuBTC and ZrBDC powders by applying pressure up to 160 and 30 MPa, respectively, both in binder-free form and mixed with polyvinyl alcohol (PVA) as a binder, allowed us to prepare the materials with sufficiently high and reproducible performance characteristics.The methane adsorption capacities of thus prepared materials were measured at an experimental ANG bench at 3.5 and 10 MPa and compared with that calculated using the Dubinin theory of volume filling of micropores (TVFM) [59,60].The experimental values occurred to be close to the calculation results at high pressure, whereas at 3.5 MPa, a noticeable discrepancy between them was found.The obtained results made it possible to consider the MOF-based materials as promising adsorbents for ANG technology.
Further development of these investigations requires working out optimal conditions for preparing the functional materials from ZrBDC intended for ANG systems.For this purpose, it is necessary to justify the choice of a binder and regime of compaction pressures based on how much the textural properties of the resulting materials change.
In the present study, we examine the factors, which determine the performance of such an ANG system.First, we synthesize chemically and thermally stable ZrBDC as a basis for preparing functional materials intended for loading the ANG tank (adsorber).Second, we focus on identifying the optimal conditions for producing the MOF-based adsorbent by the simplest way of shaping-mechanical pressing.For this purpose, we compare the results of comprehensive experimental research of textural and methane adsorption properties of the as-prepared ZrBDC powder and ZrBDCbased materials shaped without any binder and with two kinds of the binder: polyvinyl alcohol (PVA) and carboxymethylcellulose (CMC).The changes in the textural properties can be detected by structural investigations, including those performed by such methods as low-temperature nitrogen adsorption, X-ray diffraction (XRD), and scanning electron microscopy (SEM).
Direct measurements of the amount of methane delivered from a mini (100 cmcm 3 ) ANG storage tank loaded with the ZrBDC pellets with PVA over a wide temperature range and pressures up to 10 MPa provide a possibility to evaluate the real efficiency of the adsorbent and compare it with the data calculated from the experimental methane adsorption isotherms for the MOF powder and binder-free shaped material.In our work, we measure adsorption of the main NG component-methane, although there are many other impurities (C 2 -C 5 hydrocarbons, carbon dioxide, water, sulfur compounds), which affect the performance of adsorbent in the ANG system [61].For example, Zhang et al. found that HKUST-1 lost the adsorption capacity by 32% over 200 sorption/desorption cycles due to irreversible blocking of pores by propane molecules [62].Although similar cycling investigations in the presence of all NG components are determinative for designing ZrBDC-based materials with the long operational life necessary for their use in an NG-powered vehicle, these questions remained outside the scope of our work.
In the rest of the paper, we concentrated on thermal effects arising during the exothermic/endothermic adsorption/desorption processes, which are also essential for ANG storage performance [63][64][65][66].Indeed, the heat release/absorption effects during the charge/discharge processes change the temperature of the adsorbent, thereby lowering the gas storage performance below the isothermal performance [63,64] and increasing the duration of the charge/discharge processes [65].On the other hand, the heat of adsorption depends on the strength of the adsorbateadsorbent interactions and is used to get deep insight into the adsorption process [67][68][69].Therefore, keeping two goals in mind, we calculate the thermodynamic functions of the CH 4 /ZrBDC material adsorption system from the methane adsorption data in wide ranges of temperatures and pressures for the powdered ZrBDC and ZrBDC pellets.First, by evaluating the differential isosteric heat of adsorption, on the one hand, we calculate temperature changes in the ANG storage tank system during the charging process [65,69] and thereby solve a problem of practical relevance.Second, the thermodynamic characteristics of adsorption equilibrium as the entropy, enthalpy, and heat capacity calculated as a function of temperature and amount of adsorbed methane enable us to solve a purely scientific problem, namely, to examine the course of a complex adsorption process, taking into account the changes in the state of adsorbed molecules in pores.

Materials and Methods
2.1.Synthesis of ZrBDC.ZrBDC was synthesized by the solvothermal method [30] by mixing 0.34 g BDC (C 6 H 4 (COOH) 2 ) (99 + %, Acros Organics) in 30 mL N,Ndimethylformamide (DMF) (chemically pure according to the Russian State Industry Standards or GOST 20289-74, Spectrokhim) with 0.53 g ZrCl 4 (98%, Acros Organics) in a glass flask.The solution was stirred for 15 min using a magnetic stirrer without heating.Then, the mixture was placed into an autoclave, and the synthesis was carried out for from the initial solution by vacuum filtration and repeatedly washed by hot (~60 °С) DMF.The sample was dried at room temperature for 24 h and afterward in a ventilated thermostat at 110-160 °С for 2-5 h.

Preparation of the ZrBDC-Based Shaped Materials.
A hydraulic press was used to prepare the shaped MOFbased materials by compressing the ZrBDC powder in a cylindrical press die mold at the compaction pressure from 30 to 240 MPa for 2 min.The Zr-BDC powder was shaped with no binders (ZrBDC-F) and with polymer binders: a 5% PVA aqueous solution (ZrBDC-PVA) and a 1% carboxymethylcellulose sodium salt solution (ZrBDC-CMC).The choice of these binders as optimal was based on the results of the previous study [57].The MOF:binder mass ratio of the shaped mixture was 1 : 1, which, as shown in [53], allowed one to increase the density and strength of the MOF material with minimal changes in its porous structure (specific surface and micropore volume).
The shaped samples were dried in an oven at 110 °С for 12 h.The duration and temperature were determined as appropriate from a prior experiment, including periodic weight measurements of a probe sample.The resulting ZrBDC-based samples were cylindrical pellets with a diameter of 38 mm and a height of 12 mm.These dimensions provided the minimal void space in the ANG vessel (adsorber).The linear dimensions of the pellets were determined by averaging over five measurements according to the standard requirements [70].
Figure 1(a) shows the MOF samples shaped with PVA as a binder-ZrBDC-PVA before loading to the ANG adsorber (b).
Prior to loading the adsorber, the ZrBDC-PVA pellets were regenerated at 160 °C down to the pressure of 10 Pa.
2.3.Adsorptive.The adsorptive used was high purity (99.999%) methane purchased from the Moscow Gas Processing Plant (Russia).According to the data reported in [71], methane has the following physicochemical properties: molar mass μ = 16:0426 g/mol; critical temperature T cr = 190:56 K; critical pressure p cr = 4:599 MPa.The compressibility and enthalpy of gaseous methane were calculated from the state equations with the virial coefficients provided in the handbook [72].
2.4.Characterization.Elemental composition and surface morphology of the adsorbents were examined by SEM using a Quanta 650 FEG microscope (FEI Company, USA) equipped with an Oxford Energy Dispersive X-ray (EDX) detector.The elemental composition of the sample was determined by averaging ten measurements.
Porous characteristics of the as-prepared and shaped ZrBDC materials were estimated from the adsorption/ desorption isotherm of standard nitrogen vapors at 77 K measured by a Quantachrome Autosorb iQ multifunctional surface area analyzer.The structural and energy characteristics (SEC), such as the specific micropore volume W 0 [cm 3 / g], standard characteristic energy of nitrogen adsorption E [kJ/mol], and average effective radius of micropores x 0 [nm], were evaluated by using the theory of volume filling of micropores (TVFM) [59,60].The specific surface area (S BET ) of the adsorbents was estimated using the Brunauer-Emmet-Teller (BET) equation [73].In addition, the specific surface area of mesopores, S ME , was also calculated from the nitrogen adsorption/desorption data at 77 K using the well-known Kelvin equation [74].The specific mesopore volume was calculated as W ME = W S − W 0 , where W S is the total pore volume calculated from the nitrogen adsorption at the relative pressure p/p s = 0:99.Here, we have provided the data for the ZrBDC sample activated at 160 °C under conditions of thermal-vacuum regeneration.This sample demonstrated the highest SEC values.
XRD patterns of as-prepared ZrBDC were acquired on an Empyrean (Panalytical BV) diffractometer in Bragg-Brentano geometry using a nickel-filtered CuKα (λ Cu = 0:1542 nm) radiation in the 2θ angular range of 0 The packing density of as-prepared ZrBDC powder, d, was measured following the standard method [75].
2.5.Methane Adsorption Measurements.High-pressure methane adsorption was measured by the combined volumetric-gravimetric method using the original adsorption setup designed in IPCE RAS [76].The scheme is provided below in Figure 2.
Before the measurements, the adsorbent was regenerated by thermal evaporation at 160 °C to the residual pressure of no more than 1 Pa.
The value of absolute methane adsorption in the MOFbased adsorbents, а, was evaluated as follows: Here, N is the amount of methane introduced into the volume of the measuring unit, g; V is the entire volume of the system, cm 3 ; V a is the volume of the adsorbent with micropores, cm 3 ; ρ g is the density of a gaseous phase at specified pressure P and temperature T, g/cm 3 ; m 0 is the mass of regenerated adsorbent, g; μ is the molar mass of methane, g/ mmol.The value of V а was calculated as a sum of the adsorbent skeleton volume determined from helium picnometry measurements, V He , and the product m 0 × W 0 .The value of W 0 was calculated from the standard low-temperature nitrogen adsorption data by TVFM equations [59,60].
The accuracy of methane adsorption measurements was ±5% at a confidence level of 0.95 for most experimental values.The uncertainty of adsorption measurements was calculated following the standard requirements [70].
The deliverable volumetric methane capacity of the shaped MOF, i.e., the amount of methane delivered from an ANG vessel when the pressure dropped from 10 to 0.1 MPa, was measured on an ANG bench according to the recently reported route [77] (see Figure 3).
The main element of the scheme is a stainless steel ANG tank or adsorber (1), see Figure 1  compressor (5), and a high-pressure receiver ( 6), supplied the adsorber with methane.The ANG bench involves several gas-lock valves V1-4 connected by the high-pressure pipelines.With valve V4 closed, valves V1-V3 were opened, and methane was pumped to the adsorber up to a required pressure (1, 3.5, 7, and 10 MPa).The pressure regulators PR1-2 ensure the required level of pressure with an accuracy of ±0.15%.After achieving the equilibrium in the system, valves V1-V3 were closed, V4 was opened, and the delivery was started.The process was carried out at a rate of about 10 L/min until the excess pressure of 10, 7, 3.5, and 1 MPa decreased to 0.1 MPa.The FQI gas counter readings were taken only after reaching equilibrium pressure and temperature.The relative error in measuring the amount of gas by FQI was ±1%.
To take into account a dead volume, before the measurements, the system was calibrated at the same thermodynamic conditions using a mockup (MO) mounted inside the ANG vessel.The mockup volume was equal to that of the stack of the ZrBDC pellets.The correction to the deliverable storage capacity was calculated as the difference between the FQI readings: V FQI−ANG and V FQI−MO .
Before the measurements, three cycles of methane displacement adsorption up to 25 MPa were carried out to prepare the bench.

Results and Discussion
3.1.SEM Analysis of as-Synthesized ZrBDC.As follows from the SEM images at various magnifications (Figures 4(a) and 4(b)), the surface of as-prepared ZrBDC is formed by the interlocking of cubic micron-sized particles.
As follows from Table 1, which summarizes the SEM-EDX data, the chemical composition of the ZrBDC sample is determined by the reagents used in the synthesis.It can be seen that carbon is the dominant element in the sample.The presence of O, Сl, and Zr atoms in the pores of ZrBDC can contribute to the methane adsorption, enhancing the СН 4 -ZrBDC interactions.Although, as shown by Kumar et al. [78] and Vandenbrande et al. [79], the narrowest constricted micropores are the preferential adsorption sites of methane in UiO-type MOFs.In contrast, Zr ions do not contribute any significant energy heterogeneity within the pores for methane adsorption.Moreover, the inorganic building units can hinder the access of methane to the Zr clusters [78].

The Density of the Shaped ZrBDC-Based Adsorbents.
Figure 5 shows the densities of the shaped ZrBDC samples prepared without any binder and mixed with CMC and PVA versus the compaction pressure.For comparison, the star symbol in the plot marks the packing density of assynthesized powder ZrBDC (4).
As follows from Figure 5, the packing density of the shaped samples increases with the compaction pressure.The most significant increase in the density from 0.3 to 1.6 g/cm 3 was achieved upon the shaping process without any binder (see curve 1 for ZrBDC-F) at the compaction pressure of 240 MPa.Under the same conditions, the compaction of ZrBDC mixed with PVA or CMC increased the density of the final materials by about 3.7 times (from 0.3 to 1.1 g/cm 3 ).It should be noted that when the applied pressure varied from 30 to 50 MPa, the densities of the shaped materials could be categorized into ascending order: ZrBDC − F < ZrBDC − CMC < ZrBDC − PVA, but when the pressure increased from 70 to 240 MPa, it transformed into descending order as follows: ZrBDC − F > ZrBDC − PVA > ZrBDC − CMC.As follows from the data in Figure 6(a), the compaction of the ZrBDC powder without a binder produces a noticeable degradation of its porous structure: when the compaction pressure was 240 MPa, the specific BET surface decreases almost threefold: from 1100 to 340 m 2 /g.The recent XRD experiments have illustrated the degradation of the ZrBDC structure upon its compaction under different pressures [58].Thus, one can expect that the high-pressure shaping of ZrBDC results in a significant loss in adsorption capacity for methane.However, when the compaction pressures are relatively low, 30-50 MPa, the resulting materials ZrBDC-F and ZrBDC-PVA retain the highly developed porous structure.This fact is illustrated by Table 2, which lists the porous structure parameters of the ZrBDC-based materials before and after the compaction at 30 MPa in the binder-free form and mixed with PVA and CMC.As is evident from Table 2, the compaction of ZrBDC with PVA at 30 MPa led to the increase in the density by more than two times, which was accompanied by a 13.6% reduction in the micropore volume and a 6.5% reduction in the specific surface area.Thus, we obtained an adsorbent, ZrBDC-PVA, with the lowest degradation of the porous structure compared to CMC.Moreover, note that the superiority in the density observed for ZrBDC-CMC is negligible compared to ZrBDC-PVA.
Here, E 0 is the characteristic energy of standard benzene vapors calculated from the ratio Е 0 = Е/β, where β = П/П 0 = 0:34 is the affinity coefficient determined as the ratio of parachors of the studied vapor (nitrogen) П to standard benzene vapor П 0 at the boiling points [60].
The comparative analysis of the textural properties of the as-prepared ZrBDC powder and shaped materials in Table    In contrast, CMC molecules with a larger size than PVA [80] were able to enter only mesopores, and hence the lowest mesopore volume was found for ZrBDC-CMC, whereas the micropore width in this material coincided with that observed for the as-prepared ZrBDC.As follows from the data for ZrBDC-PVA, the PVA molecules were not able to block mesopores due to the smaller molecular sizes [81] (iii) Among the three shaped MOFs, two materials can be considered promising adsorbents for ANG systems [2].The relatively developed porous structure of ZrBDC-F is expected to provide its significant adsorption capacity for methane, whereas the high packing density of ZrBDC-PVA ensures the high performance for methane storage 3.4.XRD Analysis of ZrBDC-Based Materials.The main reflections at 2θ~7, 8, 12, and 25 °on the XRD pattern for the as-prepared ZrBDC (Figure 7, curve 1) coincide with that calculated for UiO-66 [49].In the previous work [57], the XRD data for ZrBDC powder compacted without any binders under different pressures revealed that the main reflections became less pronounced as the applied compaction pressure increased, indicating a progressive degradation of the ZrBDC-PVA crystalline structure.The comparison of XRD patterns recorded before and after the shaping of the as-prepared ZrBDC powder with PVA and CMC (see Figure 7) enabled us to observe the effects of the binders on the MOF structure upon their compaction at 160 Pa.
The insignificant changes in the XRD patterns of the ZrBDC shaped with PVA and CMC lead us to conclude that the binder helps minimize the destruction of the MOF structure during the shaping process.It should be noted that both binders produced almost the same "protective" effect for the ZrBDC crystalline structure.Nevertheless, PVA was chosen   Figure 8 shows the isotherms of methane adsorption on the as-prepared ZrBDC sample measured up to the pressure of 10 MPa within a range of temperatures from 253 to 333 K.All isotherms show no adsorption-desorption hysteresis loop, indicating a reversible process of physical adsorption.The experimental isotherms were fitted by a formula derived by Bakaev using statistical thermodynamics for a model of adsorption in cavities considered as quasi-independent subsystems of the grand canonical ensemble [82]: where k 0 characterizes an adsorption system, k 1 , k 2 , and k 3 are the temperature-dependent and numerically adjusted coefficients, and P is the equilibrium pressure expressed in Pa.This formula has been successfully applied to describe the experimental adsorption isotherms of different gases on zeolites [82], polymer adsorbents [83], and activated carbons [84,85].In our experiments, the regression error was not more than 2%.
To assess the effectiveness of an ANG storage system, the total specific volumetric capacity of the system, V F , is employed.This parameter involves the amount of adsorbed gas and the amount of gas in the intergranular space outside the adsorbent: where ρ g is the density of a free gaseous phase in an ANG tank (kg/m 3 ) at the specified pressure P and temperature T ; d is the packing density of the adsorbent (kg/m 3 ), NTP is the standard thermodynamic conditions for gases: T = 293:15 K, P = 101:325 kPa, ρ NTP is the density of the gas at the NTP conditions (kg/m 3 ), and ε is the adsorbent bed void fraction in the ANG tank, which was calculated as follows: Here, d c ≈ 1650 kg/m 3 is the skeleton density of ZrBDC [26].
The deliverable volumetric capacity of the ANG system, V A , is more important for evaluating its performance since it determines the driving range of NG vehicles.The value of V A is evaluated as the amount of gas at standard thermodynamic conditions, which can be delivered from the tank with adsorbent upon isothermal gas discharge from the current level of pressure P in the tank to the minimal value P x (in our case, P x = 0:1 MPa) normalized to the geometric volume of the tank: 9 Adsorption Science & Technology compared with the specific volumetric capacity of a CNG system calculated for the same thermodynamic conditions (dashed lines in Figure 9).
It was found that the total specific volumetric capacity of the ANG system with ZrBDC achieved 175 m 3 (NTP)/m 3 at 293 K and 10 MPa (Figure 9(a)).This value occurred to be less than that for the ANG system with ZrBDC-F, which exceeded 190 m 3 (NTP)/m 3 (Figure 9(b)).The reason is an almost twofold increase in the packing density of the adsorbent due to the compaction (see Table 2).Moreover, the regime of the maximum absolute efficiency of the ANG system with ZrBDC-F, i.e., when V ANG F − V CNG F > 0, is observed within the pressure range from 1 to 6 MPa and is almost double this parameter for the system with ZrBDC.
A similar tendency is also observed for the relationship between the specific volumetric deliverable capacities of the ANG systems with as-synthesized powder ZrBDC and compacted binder-free ZrBDC-F (compare Figures 10(a) and 10(b)).It should be noted that at 253 K and when the pressures are above 10 MPa, the values of V A ANG for the ANG system with ZrBDC are less than V A CNG .Figure 10(c) shows the values of V A ANG obtained on the experimental ANG bench with the adsorber loaded with ZrBDC-PVA pellets.The data on V A ANG are represented as a function of pressure at four temperatures: 253, 273, 293, and 303 K.
It was found that at 10 MPa and 253 K V A achieved the value of 172 m 3 (NTP)/m 3 for the ANG system with ZrBDC-PVA and 155 m 3 (NTP)/m 3 for ZrBDC.Note that at the same P,T-conditions, the ANG system with the ZrBDC-F adsorbent shows the deliverable capacity of 160 m 3 (NTP)/m 3 , whereas this value exceeds that for ZrBDC-PVA by 3-5% within the range of pressures from 2 to 6 MPa.However, the relatively poor mechanical properties of ZrBDC-F will limit its practical application.Indeed, the Shore hardness measurements [86] carried out for the shaped ZrBDC-PVA, ZrBDC-CMC, and ZrBDC-F materials gave the values of 95, 96.5, and 50 Shore A units, respectively.
As shown above, the compaction of the ZrBDC powder, on the one hand, increased its density, but, on the other hand, the parameters of the porous structure deteriorated.These factors had opposite effects on methane adsorption.As follows from the obtained results, the ZrBDC-PVA pellets prepared under the pressure of 30 MPa are characterized by increased density and partially degraded porosity, but they demonstrated higher performance than the MOF powder.The compaction at 30 MPa with PVA improved the deliverable methane capacity of ZrBDC at 10 MPa and 253 K by ~11%.

Thermodynamic Properties of the СН4/ZrBDC
Adsorption System 3.6.1.Thermodynamic State Functions.As known [63][64][65][66]69], the effectiveness of the ANG system is determined not only by the adsorption properties of the porous material but also by the heat effects arising during the methane adsorption/desorption processes.The adsorbent bed in the ANG tank releases heat during the charging process.In contrast, it gets cool during discharge.In addition, if the temperature in the adsorber exceeds the maximum endurance temperature, then the adsorber would not be able to withstand the pressure, and there may be a release of the fuel.Thermal management of an ANG system (heat exchangers, material construction, and geometry of the adsorber) requires the knowledge of thermodynamic characteristics of the adsorption system, namely, the heat of adsorption.The latter provided a rationale to evaluate the thermodynamic state functions of the investigated adsorption systems.
The differential molar isosteric heat of adsorption, q st , is calculated as the difference between the molar enthalpy of   5), 323 (6), and 333 (7).Dashed lines show the dependences for a CNG system at the same P, T-conditions.11 Adsorption Science & Technology an equilibrium gas phase, h g , and the differential molar enthalpy of an adsorption system, h 1 [87]: Here, H 1 is the enthalpy of the adsorption system.
We employed an approach developed by Bakaev [88] and reported in details in [66,89] for calculating the differential molar isosteric heat of adsorption from the experimental isosteres of absolute methane adsorption with consideration of the nonideality of the gaseous phase and noninertness of the adsorbent: Here Z = Pv g /RT is the compressibility of the equilibrium gaseous phase at the specified P, T-conditions (Pa, K); v g is the specific volume of the gaseous phase, m 3 /kg; R is the universal gas constant, J/(mol•K); V a = V 0 ðР, ТÞ/m 0 is the specific reduced volume of the adsorption system, cm 3 /g; V 0 ðР, ТÞ is the volume of the adsorbent including micropores; m 0 is the mass of regenerated adsorbent.Thus, formula (7) includes the term associated with the isothermal adsorption-induced deformation ð∂V a /∂aÞ T and isosteric deformation ð∂V a /∂TÞ a of the adsorbent, the slope of the adsorption isotherms ð∂P/∂aÞ T and isosteres ½∂lnP/∂ð1/TÞ a , and the nonideality of the gaseous phase Z ðР, ТÞ.
When the effects associated with thermal and adsorption noninertness (deformation) of an adsorption system are insignificant, Eq. ( 7) simplifies to It should be noted that ZrBDC-based MOF is considered exceptionally rigid and thermally stable due to the high degree of coordination of the Zr-O metal centers to the organic linkers [90].Therefore, the adsorption-and temperature-stimulated deformations of the ZrBDC material, (∂V a /∂a) T and (∂V a /∂T) a could be neglected when calculating the differential isosteric heat of adsorption.
The heat of adsorption is determined by the adsorbateadsorbent and adsorbate-adsorbate interactions in the adsorption system.Therefore, we plotted the heats of methane adsorption in the as-prepared ZrBDC powder (see Figure 11(a)) and the ZrBDC-PVA pellets (Figure 11(b)) as a function of methane uptake over a wide range of temperatures in order to compare the mechanism of methane adsorption in these Zr-based materials.
As follows from Figures 11(a) and 11(b), at all measured temperatures, the initial differential molar isosteric heats of adsorption are about 19.3-18.9kJ/mol.These values are close or exceeded that determined for methane adsorption in most MOF materials [67,91], including Zr-based MOFs as Nu-1000 (16-18 kJ/mol [92] and UiO-type MOFs [78,93]).The comparison of Figures 11(a) and 11(b) revealed the discrepancy between the functions q st = f ðaÞ for the ZrBDC powder and ZrBDC-PVA pellets.The variation in the q st with an increase in the amount of adsorbed methane molecules reflects the changes in the mechanism of methane adsorption determined by the ratio between СH 4 -MOF and СH 4 -СH 4 interactions.The decrease in the value of q st from the initial value at the early stage of adsorption (up to a ~1 mmol/g) to about 12-16 kJ/mol (at a = 4 − 5 mmol/ g) observed in the pristine ZrBDC powder is caused by the gradual occupation of the adsorption sites with different interaction energies by methane molecules, starting the narrowest pores (so-called adsorption pockets with one methane molecule [94]) and Zr ions [78].
In contrast, the methane adsorption in the ZrBDC-PVA pellets is characterized by the extreme function q st = f ðaÞ with a wide maximum, the amplitude, width and position of which depend on temperature (see Figure 11(b)).Then, the heat of adsorption declines rapidly at higher values of methane uptake.In our opinion, these behaviors can be ascribable to the variations in the porous structure of ZrBDC caused by the compaction with PVA, namely the distribution and energy level of the adsorption sites.We suppose that the interactions between the binder macromolecules and methane molecules contribute to the heat of adsorption.It should be noted that Stadie et al. observed similar anomalous behaviors for methane adsorption in zeolite-templated carbons and attributed it to the relatively high homogeneity of binding site energies caused by a narrow pore size distribution, periodic spacing, and high content of sp 2 -hybridized carbon [95].We agree with the authors of this work that correct estimations of the impacts of intermolecular interactions to the isosteric heat of methane adsorption in ZrBDC-PVA require additional investigations of the adsorption energy of the different adsorption sites.
As follows from Figures 11(a) and 11(b)(curves 2-8) and (curve 7), the drop of q st slows down within the range of medium values of 3-5 or 3 mmol/g, which indicates an increase in the contribution from the СH 4 -СH 4 interactions.The effect of lateral СH 4 -СH 4 interactions during adsorption in different microporous adsorbents, including NaX zeolite [96], activated carbons [97,98], and MOFs [99,100], was shown experimentally and confirmed by molecular dynamics.The decrease in the heat of adsorption observed at the adsorption values above 6-8 mmol/g and 253-333 K can be attributed to a decrease in the intensity of the adsorption process (∂P/∂a increases), accompanied by the enhancement of repulsion due to the reduced distances between methane molecules and pore walls [101].As the temperature rises, these effects are observed at lower pore loadings due to the increasing contribution from kinetic energy associated with molecular mobility.
At the initial stage of micropore loading: a < 1:2 mmol/g for as-prepared ZrBDC and a < 0:7 mmol/g for ZrBDC-PVA, the isosteric heat of adsorption is independent of temperature (see Figures 11(a) and 11(b)).Adsorption Science & Technology However, with a subsequent increase in methane adsorption, the temperature dependence of the isosteric heat of adsorption is observed for both adsorbents, and it becomes increasingly noticeable as the pore space is filled with methane.The observed divergence in the curves q st = f ða Þ can be interpreted in terms of Eq. ( 8) by the significant changes in the contributions from the slope of isosteres and the compressibility of the gaseous phase Z with temperature.
These conclusions are confirmed by the temperature dependences of the differential molar isosteric enthalpy of the ZrBDC-CH 4 adsorption system, h 1 , which can be calculated from Eq. ( 6) using the data of differential isosteric heat of adsorption.Figure 12 shows the temperature dependences of the differential molar enthalpy of the ZrBDC-CH 4 system for different values of methane adsorption.At the initial stage of adsorption, the enthalpy of the adsorption system increases almost linearly with temperature, which is associated with the temperature invariance of the heat of adsorption, q st ≠ f ðTÞ, and the linear character of the temperature dependence of the molar enthalpy of an equilibrium gas phase, h g .However, the h 1 ðTÞ dependences become nonlinear during the methane adsorption: the slope ∂h 1 ðTÞ/∂T increases with temperature.
The observed variations of h 1 during the process of methane adsorption in ZrBDC, which were evaluated from the experimental isotherms of adsorption, are consistent with the computational studies of the enthalpy variations for methane in nanoporous solids.In particular, using the continuous fractional component Monte Carlo algorithm for calculating thermodynamics for the methanezeolite adsorption system, Torres-Knoop et al. showed that the repulsive interactions between an adsorbate molecule and the previously adsorbed species and the adsorbent were responsible for the rearrangements of methane molecules, which lead to the increase in the enthalpy at high micropore loadings close to saturation [101].
The differential molar entropy of the adsorbed methane molecules calculated as and represented as a function of methane uptake in Figure 13 reflects the variations in the state of methane molecules during adsorption compared to a gaseous phase.
As can be seen in Figure 13, at the adsorption value of about 1 mmol/g, there is a sharp decrease in entropy, indicating the transfer of methane molecules with three translational, three rotational, and six vibrational degrees of freedom to the bound state on the high-energy adsorption sites.The state of methane molecules adsorbed onto MOF is characterized by reduced mobility, which can be described as probably three rotational and three vibrational degrees of freedom or six vibrational degrees of freedom [99].In the range of adsorption values of 2-6 mmol/g, the drop in entropy slows down during the gradual methane loading of the pore space and the impact of the CH 4 -CH 4 interactions.An increase in the s 1 ðaÞ values is observed when the pore loadings are above ~6 mmol/g, which we attribute to the rearrangement of adsorbed methane molecules close to saturation.
The heat capacity of the CH 4 /ZrBDC adsorption system is essential for the characterization of the system since it determines the thermal stability due to thermal fluctuations during adsorption/desorption cycles.Unlike the heat capacity of any bulk phase, this parameter of the adsorption system depends not only on the pressure and temperature but also on the amount of adsorbed substance.Therefore, a complete thermodynamic description of the adsorption system involves the analysis of the differential molar isosteric heat capacity, C a , as a function of all these parameters.The value of C a was calculated using the Kirchhoff equation obtained by differentiating Eq. ( 11) with respect to  6), and 333 (7), calculated using Eq. ( 7).The dashed line 8 depicts the value of q st methane adsorption onto as-prepared powder ZrBDC when Z = 1 and ð∂P/∂aÞ ≈ 0.
13 Adsorption Science & Technology temperature [87]: Figure 14 demonstrates the temperature dependence of the differential molar isosteric heat capacity of the ZrBDC-CH 4 for various amounts of adsorbed methane.
As follows from Figure 14, the differential molar isosteric heat capacity of the CH 4 /ZrBDC adsorption system depends on temperature and increases with adsorption.At low values of methane adsorption, the increase in the isosteric heat capacity of the adsorption system is determined by the temperature dependence of the enthalpy of methane in a gaseous state since, as was shown above: q st ≠ f ðTÞ.However, as follows from Figure 11(a), during the process of methane adsorption, the temperature invariance of q st is replaced by the explicit temperature dependence of q st , which begins to  Adsorption Science & Technology contribute to the heat capacity of the adsorption system C a = f ðTÞ.The heat capacity of the CH 4 /ZrBDC adsorption system significantly exceeds (by ~2-3 times) the isochoric heat capacity of the gaseous methane phase, C V .This fact indicates a specific highly dispersed state of adsorbed methane in a strong adsorption field generated by a microporous adsorbent [87].
3.7.Thermal Effects in the Model ANG Systems with the ZrBDC-Based Adsorbents.In addition to the thermodynamic state functions, the data on differential isosteric heat of methane adsorption onto the ZrBDC-based materials were used to access thermal effects in the model ANG vessels arising during the adsorption (charging) process.For this purpose, we calculate the integral heat of adsorption Q released during adsorption [85]: Then, the values of Q were used to calculate the variation in temperature (ΔТ) of a model ANG system filled with the ZrBDC powder and the ZrBDC-PVA pellets during adsorption under adiabatic conditions.The heat capacities of the ZrBDC-based materials were evaluated using the empirical Neumann-Koop rule equation [102] and the data in Table 1.The results of the calculations are listed in Table 3.
According to Table 3, an increase in the temperature of the adsorption process is accompanied by a decrease in the integral heat released in both adsorption systems.A difference between the values of Q observed in the ANG systems with the powdered ZrBDC adsorbent and ZrBDC-PVA pellets reached 1.1 kJ at 273 K, and it decreased to 0.4 kJ at 333 K.A slight excess in the heating-up of the ANG system with the pellets can be ascribed to the lower heat capacity of ZrBDC-PVA.Also, the data in Table 3 testified that the heating-up of the ANG system is maximal at low temperatures.Thus, the ZrBDC-PDA pellets demonstrated superiority in the deliverable capacity for methane compared to  ZrBDC powder (see Figures 10(a) and 10(c)), which was not accompanied by a considerable increase in the temperature of the adsorbent.Indeed, the heating-up of the ANG system loaded with the MOF-based materials occurred much less than that for activated carbons: ΔT = 53 K, obtained for the values of methane adsorption of 8 mmol/g at 5 MPa and 273 K [85].Thus, the thermal effects arising during the adsorption (charging) process do not threaten the integrity of the MOF-based adsorbent during cyclic operation.

Conclusions
To conclude, we explored two factors affecting the performance of the ANG system with a MOF-based adsorbent: textural characteristics of the material, including its porous structure and packing density, and thermodynamic behaviors of the entire adsorption system as a whole.These factors affect the deliverable storage capacity and thermal management of the ANG system.The exploration and comparison of the structural and adsorption properties of a set of Zr-based MOFs: asprepared ZrBDC and three shaped ZrBDC materials, which were prepared by compacting without a binder and mixed with polymer binders such as PVA and CMC, over a range of compaction pressure from 30 to 240 MPa and at 298 K, combined with the data of the tests of the lab-scale ANG tank loaded with the shaped ZrBDC-based pellets brought us to the following conclusions: (i) The increase in the ZrBDC material density depends on the compaction pressure and the use of a binder.The highest density of the ZrBDCbased materials was achieved by compacting the as-prepared MOF powder under the pressure of 240 MPa without any binders: 1.6 g/cm 3 .The compaction of the ZrBDC powders mixed with the binders PVA or CMC under 240 MPa also led to a more than three-fold increase in their density: up to 1.1 g/cm 3   (ii) Although the binders retain the crystal structure of ZrBDC, they affect the degree of degradation of the porous structure of the shaped ZrBDC material.The proportion of micro-and mesopores in the shaped ZrBDC-based material occurred to be dependent on the ratio between the molecular size of a binder and pore diameter.PVA was found to be the best binder since it provided a possibility to achieve a minimal loss of micropore volume and specific surface area of the final ZrBDC-based upon compaction (iii) The compaction of the ZrBDC powder at a minimum pressure of 30 MPa and the use of PVA as a binder, leading to a more than two-fold increase in the packing density and lowest degradation of the porous structure, was sufficient to improve the methane adsorption capacity of the final adsorbent by ~11%: up to 172 m 3 (NTP)/m 3 between 0.1 and 10 MPa at 253 K.This value exceeded that observed for the as-prepared ZrBDC powder and compacted ZrBDC-F material.Therefore, the ZrBDC-PVA material can be considered as a promising adsorbent for small-scale ANG systems Thus, the results obtained in this study demonstrated the possibility of such techniques of fabrication of MOF-based adsorbents as compaction for improving their performance for ANG systems.In this regard, it is expected that further progress in fabricating the optimal ZrBDC-based materials which meet the DOE target can be ensured by combining two options: (i) functionalization aimed at an increase in the energy of adsorption and (ii) shaping of these functionalized materials.
The results of measurements of the methane adsorption onto the as-prepared ZrBDC powder and ZrBDC-PVA pellets over the temperature range from 253 to 333 K served to evaluate the thermodynamic state functions of the adsorption systems, primarily the differential molar isosteric heat of adsorption.The analysis of the thermodynamic functions versus the value of methane uptake within the range of temperatures yielded the following results: (i) The compaction with PVA altered the behaviors of the differential molar isosteric heat of methane adsorption onto ZrBDC as a function of methane uptake, which can be ascribed to the changes in the distribution and energy level of adsorption sites in the adsorbent, including narrow pores and access to metal ions Evaluating the key thermodynamic quantities is essential for selecting an efficient adsorbent for the ANG systems, novel gas-power drones and vehicles, and corresponding heat-exchange facilities to achieve the maximum storage density and deliverable capacity.

Figure 1 :
Figure 1: Photographs of the ZrBDC-PVA pellets (a) prepared by compacting the ZrBDC powder under the shaping pressure of 30 MPa at the temperature of 298 K for loading into the ANG tank (b).

Figure 4 :
Figure 4: SEM images of the as-prepared ZrBDC sample.The scale bars are 20 (a) and 5 μm (b).

Figure 5 :
Figure 5: Density of the ZrBDC-based materials prepared without any binders, ZrBDC-F (1), and mixed with the binders: ZrBDC-PVA (2) and ZrBDC-CMC (3), versus the compaction pressure.Orange star denotes the packing density of as-prepared powder ZrBDC.The temperature of measurements was 298 K.

Figure 6 :
Figure 6: Specific surface area S BET (a) and micropore volume W 0 of the ZrBDC-based materials: ZrBDC-F (1), ZrBDC-PVA (2) and ZrBDC-CMC (3), versus the compaction pressure.The star symbol (4) shows the values of S BET and W 0 of the as-prepared ZrBDC powder.

Figures 8 (Figure 7 :
Figures 8(a) and 8(b) show the total specific volumetric capacities of the adsorption systems with the ZrBDC (a) and ZrBDC-F (b) materials used as adsorbents represented as a function of pressure within the temperature range of 253-333 K.The data were calculated from the experimental isotherms of methane adsorption according to Eq. (3) and

(
ii) The variations in the thermodynamic state functions of the CH 4 /ZrBDC adsorption system during adsorption, including the differential molar isosteric heat of adsorption, enthalpy, entropy, and heat capacity, can be attributed to the changes in the ratio of the CH 4 -adsorbent and CH 4 -CH 4 interactions, which altered the state of adsorbed molecules from the partially-localized state near the highenergy adsorption sites with restricted mobility to the adsorption molecular clusters, which are rearranged close to saturation (iii) The temperature changes of ANG systems loaded with the ZrBDC powder and ZrBDC-PVA pellets upon the adsorption (charging) process differ insignificantly.The relatively low values of temperature changes in the ANG tank loaded with the ZrBDC-PVA pellets indicate the high operating properties during cyclic work

Table 1 :
Elemental chemical composition of ZrBDC evaluated from the SEM-EDX data and represented in weight (wt.%) and atomic (at.%) percentage.

Table 2 :
Parameters of the porous structure and the packing density of the as-prepared ZrBDC powder and shaped materials: ZrBDC-F, ZrBDC-PVA, and ZrBDC-CMC prepared under the compaction pressure of 30 MPa and temperature of 298 K.as the most optimal binder for preparing monolithic ZrBDC-based materials on the ground of the cumulative contribution of a binder to the density, porous, and crystalline structures of the resulting MOF upon the shaping process (Table2, Figures5, 6

Table 3 :
Thermodynamic characteristics of the model ANG vessel with a mass of 10 kg loaded with the ZrBDC powder and ZrBDC-PVA pellets at the pressure of 10 MPa and the temperatures of 273, 293, and 333 K.