The development of new water adsorbents that are hydrothermally stable and can operate more efficiently than existing materials is essential for the advancement of water adsorption-driven chillers. Most of the existing benchmark materials and related systems in this field suffer from clear limitations that must be overcome to meet global requirements for sustainable and green energy production and utilization. Here, we report the energy-efficient water sorption properties of three isostructural metal-organic frameworks (MOFs) based on the simple ligand pyridine-2,4-dicarboxylate, named M-CUK-1 [M3(μ3-OH)2(2,4-pdc)2] (where M = Co2+, Ni2+, or Mg2+). The highly hydrothermally stable CUK-1 series feature step-like water adsorption isotherms, relatively high H2O sorption capacities between P/P0 = 0.10-0.25, stable cycling, facile regeneration, and, most importantly, benchmark coefficient of performance values for cooling and heating at a low driving temperature. Furthermore, these MOFs are prepared under green hydrothermal conditions in aqueous solutions. Our joint experimental-computational approach revealed that M-CUK-1 integrates several optimal features, resulting in promising materials as advanced water adsorbents for adsorption-driven cooling and heating applications.

@en

The metal–organic frameworks (MOFs) MIL-101(Cr) and NH₂-MIL-125 offer high adsorption capacities and have therefore been suggested for sustainable energy conversion in adsorption chillers. Herein, these MOFs are benchmarked to commercial Siogel. The evaluation method combines small-scale experiments with dynamic modeling of full-scale adsorption chillers. For the common temperature set 10/30/80 °C, it is found that MIL-101(Cr) has the highest adsorption capacity, but considerably lower efficiency (−19%) and power density (−66%) than Siogel. NH₂-MIL-125 increases efficiency by 18% compared with Siogel, but reduces the practically important power density by 28%. From the results, guidelines for MOF development are derived: High efficiencies are achieved by matching the shape of the isotherms to the specific operating temperatures. By only adapting shape, efficiencies are 1.5 times higher. Also, higher power density requires matching the shape of the isotherms to create high driving forces for heat and mass transfer. Second, if MOFs’ heat and mass transfer coefficients could reach the level of Siogel, their maximum power density would double. Thus, development of MOFs should go beyond adsorption capacity, and tune the structure to the application requirements. As a result, MOFs could to serve as optimal adsorbents for sustainable energy conversion.

@en

A mathematical model is used to predict adsorption isotherms from experimentally measured breakthrough curves. Using this approach, by performing only breakthrough experiments for a mixture of two (or more) components, one can obtain pure component adsorption isotherms up to the pressure of the experiment. As a case study, the adsorption of an equimolar mixture of CO₂ and CH₄ in zeolite ITQ-29 is investigated. Pure component linear adsorption isotherms for CO₂ and CH₄ are predicted by fitting the theoretical breakthrough curves to the experimental ones. Henry coefficients obtained from our approach are in excellent agreement with those measured experimentally. A similar procedure is applied to predict the complete Langmuir adsorption isotherm from breakthrough curves at high pressures. The resulting adsorption isotherms are in very good agreement with those measured experimentally. In our model for transient adsorption, mass transfer from the gas phase to the adsorbed phase is considered using the Linear Driving Force model and dispersion of the gas phase in the packed bed is taken into account. IAST is used to compute the equilibrium loadings for a mixture of gases. The influence of the dispersion coefficient and the effective mass transfer coefficient on the shape of breakthrough curves is investigated and discussed. Rough estimations of these values are sufficient to predict adsorption isotherms from breakthrough curves.

@en

A computational screening of 2930 experimentally synthesized metal-organic frameworks (MOFs) is carried out to find the best-performing structures for adsorption-driven cooling (AC) applications with methanol and ethanol as working fluids. The screening methodology consists of four subsequent screening steps for each adsorbate. At the end of each step, the most promising MOFs for AC application are selected for further investigation. In the first step, the structures are selected on the basis of physical properties (pore limiting diameter). In each following step, points of the adsorption isotherms of the selected structures are calculated from Monte Carlo simulations in the grand-canonical ensemble. The most promising MOFs are selected on the basis of the working capacity of the structures and the location of the adsorption step (if present), which can be related to the applicable operational conditions in AC. Because of the possibility of reversible pore condensation (first-order phase transition), the mid-density scheme is used to efficiently and accurately determine the location of the adsorption step. At the end of the screening procedure, six MOFs with high deliverable working capacities (∼0.6 mL working fluid in 1 mL structure) and diverse adsorption step locations are selected for both adsorbates from the original 2930 structures. Because the highest experimentally measured deliverable working capacity to date for MOFs with methanol is ca. 0.45 mL mL^-1, the selected six structures show the potential to improve the efficiency of ACs.

@en

The influence of product shape selectivity on the bifunctional conversion of n-C₇ by zeolite catalysts is investigated. Three different zeolite catalysts with different pore sizes (MFI-type, MEL-type, and BEA-type zeolites) have been investigated experimentally. For all three catalysts, n-C₇ is isomerized to monobranched isomers which are further isomerized into dibranched isomers, and these dibranched molecules are converted into cracking products. More dibranched isomers and less cracking products are produced by BEA-type zeolite compared to MFI-type and MEL-type zeolites and clear differences are observed in the distribution of dibranched isomers produced by different catalysts. Molecular simulation is used to compute the adsorption isotherms and free energy barriers for diffusion of dibranched isomers in MFI-type, MEL-type, and BEA-type zeolites. Combining simulation results and experimental observations, it is shown that product shape selectivity can explain the distribution of dibranched molecules while transition state shape selectivity fails to do so. For the medium-pore zeolites (MFI-type and MEL-type zeolites), free energy barriers for diffusion of dibranched molecules are significant. For MFI-type and MEL-type zeolites, the dibranched molecule that has to overcome lower diffusion barrier is produced with a higher yield and the distribution of dimethylpentane molecules is determined by their diffusion rate. It is shown that there is almost no free energy barrier for the diffusion of any of these molecules in BEA-type zeolite. As BEA-type zeolite imposes no free energy barrier for diffusion of any of dibranched isomers, the distribution of dibranched isomers is very close to the equilibrium distribution in the gas phase. Due to the limited mobility of dimethylpentanes within the pores of MFI-type and MEL-type zeolites, most of the dimethylpentane molecules are trapped inside the zeolite and undergo consecutive cracking. Dimethylpentane molecules diffuse sufficiently fast in the large pores of BEA-type zeolite and transfer to the gas phase, before consecutive reaction converts these molecules into cracking products. Moreover, the effect of the MFI-type crystal size on the production of dibranched isomers is investigated. The yield of dibranched isomers reduces by increasing the size of the crystal and larger part of dibranched isomers are cracked as the crystal size of MFI-type is increased.

@en

Chemical clocks are often used as exciting classroom experiments, where an induction time is followed by rapidly changing colours that expose oscillating concentration patterns. This type of reaction belongs to a class of nonlinear chemical kinetics also linked to chaos, wave propagation and Turing patterns. Despite its vastness in occurrence and applicability, the clock reaction is only well understood for liquid-state processes. Here we report a chemical clock reaction, in which a solidifying entity, metal-organic framework UiO-66, displays oscillations in crystal dimension and number, as shown by X-ray scattering. In rationalizing this result, we introduce a computational approach, the metal-organic molecular orbital methodology, to pinpoint interaction between the tectonic building blocks that construct the metal-organic framework material. In this way, we show that hydrochloric acid plays the role of autocatalyst, bridging separate processes of condensation and crystallization.

@en

In this modeling study, the uses of nitrogen (77.3 K), probe molecule of choice for decades, and argon, opted as alternative in the 2015 IUPAC report on adsorptive characterization, as probe molecules for geometric surface area determination are compared. Graphene sheets possessing slit-shaped pores with varying size (width) are chosen as model porous solids, and different methods for the determination of specific surface areas are investigated. The BET method, which is the most commonly applied analysis, is compared to the Langmuir and relatively recently proposed ESW (excess sorption work) method. We show that either using argon or nitrogen as adsorptive, the physical meaningfulness of adsorption-derived surface areas highly depends on the pore size. When less than two full layers of adsorbate molecules can be formed within slitlike pores of a graphitic material (Dpore < 5.8 Å for Ar/N2), adsorption-derived surface areas are about half that of the geometric surface area. Between two and four layers (6.8 < Dpore < 12.8 Å), adsorption surface areas can be significantly larger (up to 75%) than the geometric surface area because monolayer−multilayer formation and pore filling cannot be distinguished. For four or more layers of adsorbate molecules (Dpore > 12.8 Å), adsorption-derived surface areas are comparable to their geometrically accessible counterparts. Note that for the Langmuir method this only holds if pore-filling effects are excluded during determination. This occurs in activated carbon materials as well. In the literature, this indistinguishability issue has been largely overlooked, and erroneous claims of materials with extremely large surface areas have been made. Both the BET and Langmuir areas, for Dpore > 12.8 Å, correspond to geometric surface areas, whereas the ESW method yields significantly lower values. For the 6.8 Å < Dpore < 12.8 Å range, all methods erroneously overestimate the specific surface area. For the energetically homogeneous graphene sheets, differences between argon and nitrogen for the assessment of surface areas are minor.@en