In Silico Screening of Metal-Organic Frameworks for Adsorption-Driven Heat Pumps and Chillers

Accelerated convergence via adiabatic sampling for adsorption and desorption processes

Under isothermal conditions, phase transitions occur through a nucleation event when conditions are sufficiently close to coexistence. The formation of a nucleus of the new phase requires the system to overcome a free energy barrier of formation, whose height rapidly rises as supersaturation decreases. This phenomenon occurs both in the bulk and under confinement and leads to a very slow kinetics for the transition, ultimately resulting in hysteresis, where the system can remain in a metastable state for a long time. This has broad implications, for instance, when using simulations to predict phase diagrams or screen porous materials for gas storage applications. Here, we leverage simulations in an adiabatic statistical ensemble, known as adiabatic grand-isochoric ensemble (μ, V, L) ensemble, to reach equilibrium states with a greater efficiency than its isothermal counterpart, i.e., simulations in the grand-canonical ensemble. For the bulk, we show that at low supersaturation, isothermal simulations converge slowly, while adiabatic simulations exhibit a fast convergence over a wide range of supersaturation. We then focus on adsorption and desorption processes in nanoporous materials, assess the reliability of (μ, V, L) simulations on the adsorption of argon in IRMOF-1, and demonstrate the efficiency of adiabatic simulations to predict efficiently the equilibrium loading during the adsorption and desorption of argon in MCM-41, a system that exhibits significant hysteresis. We provide quantitative measures of the increased rate of convergence when using adiabatic simulations. Adiabatic simulations explore a wide temperature range, leading to a more efficient exploration of the configuration space.

Engineered Nanoporous Frameworks for Adsorption Cooling Applications

The energy demand for traditional vapor-compressed technology for space cooling continues to soar year after year due to global warming and the increasing human population's need to improve living and working conditions. Thus, there is a growing demand for eco-friendly technologies that use sustainable or waste energy resources. This review discusses the properties of various refrigerants used for adsorption cooling applications followed by a brief discussion on the thermodynamic cycle. Next, sorbents traditionally used for cooling are reviewed to emphasize the need for advanced capture materials with superior properties to improve refrigerant sorption. The remainder of the review focus on studies using engineered nanoporous frameworks (ENFs) with various refrigerants for adsorption cooling applications. The effects of the various factors that play a role in ENF-refrigerant pair selection, including pore structure/dimension/shape, morphology, open-metal sites, pore chemistry and possible presence of defects, are reviewed. Next, in-depth insights into the sorbent-refrigerant interaction, and pore filling mechanism gained through a combination of characterization techniques and computational modeling are discussed. Finally, we outline the challenges and opportunities related to using ENFs for adsorption cooling applications and provide our views on the future of this technology.

Avoiding pitfalls in molecular simulation of vapor sorption: Example of propane and isobutane in metal-organic frameworks for adsorption cooling applications

This study introduces recommendations for conducting molecular simulations of vapor adsorption, with an emphasis on enhancing the accuracy, reproducibility, and comparability of results. The first aspect we address is consistency in the implementation of some details of typical molecular models, including tail corrections and cutoff distances, due to their significant influence on generated data. We highlight the importance of explicitly calculating the saturation pressures at relevant temperatures using methods such as Gibbs ensemble Monte Carlo simulations and illustrate some pitfalls in extrapolating saturation pressures using this method. For grand canonical Monte Carlo (GCMC) simulations, the input fugacity is usually calculated using an equation of state, which often requires the critical parameters of the fluid. We show the importance of using critical parameters derived from the simulation with the same model to ensure internal consistency between the simulated explicit adsorbate phase and the implicit bulk phase in GCMC. We show the advantages of presenting isotherms on a relative pressure scale to facilitate easier comparison among models and with experiment. Extending these guidelines to a practical case study, we evaluate the performance of various isoreticular metal-organic frameworks (MOFs) in adsorption cooling applications. This includes examining the advantages of using propane and isobutane as working fluids and identifying MOFs with a superior performance.

The Role of Interdigitated Electrodes in Printed and Flexible Electronics

Flexible electronics, also referred to as printable electronics, represent an interesting technology for implementing electronic circuits via depositing electronic devices onto flexible substrates, boosting their possible applications. Among all flexible electronics, interdigitated electrodes (IDEs) are currently being used for different sensor applications since they offer significant benefits beyond their functionality as capacitors, like the generation of high output voltage, fewer fabrication steps, convenience of application of sensitive coatings, material imaging capability and a potential of spectroscopy measurements via electrical excitation frequency variation. This review examines the role of IDEs in printed and flexible electronics since they are progressively being incorporated into a myriad of applications, envisaging that the growth pattern will continue in the next generations of flexible circuits to come.

Exploring the Potential of Metal-Organic Frameworks for Cryogenic Helium-Based Gas Gap Heat Switches via High-Throughput Computational Screening

With the advantages of a long lifetime and high reliability, gas gap heat switches (GGHSs) are attractive in many thermal management applications, especially in space-borne cryogenic systems. The performance of a GGHS is significantly affected by the adsorption characteristics of the adsorbent in the sorption pump. Compared with the commonly used adsorbent in the GGHSs (activated carbon), metal-organic frameworks (MOFs) have larger surface areas, higher pore volumes, and exceptional tunability, which motivates this study to explore their potential for application in cryogenic GGHSs. To this end, two performance metrics, the required volume of adsorbent (vsor) and total input heat (qtot), were computed for about 6000 MOFs via molecular simulations and compared with those of activated carbon. It is found that over 2300 MOFs possess a smaller vsor than activated carbon, and the smallest vsor of MOFs is about 12.7% of that of activated carbon. vsor and qtot generally change in the same direction, which implies it is possible to reduce both parameters simultaneously by choosing a suitable MOF. Structure-performance analysis reveals that 1/vsor consistently increases first and then decreases with pore limiting diameter, largest cavity diameter, available pore volume, accessible surface area, helium void fraction, and bulk density. Descriptor ranges corresponding to high-performing MOFs were identified based on Precision-Recall analysis. Notably, Zr-containing MOFs are particularly likely to have smaller vsor values than activated carbon. It is anticipated that the promising MOFs identified by this study will motivate further experimental investigations, and the insights into structure-performance relationships can serve to guide the rational design of novel MOF candidates for GGHSs.

Ammonia Adsorption Performance of Zeolitic Imidazolate Frameworks for Cooling

Promoting the cooling performance of adsorption chillers (ACs) greatly relies on the exploration of high-performance adsorbent/refrigerant working pairs. Ammonia is not only an environmentally friendly refrigerant but also favorable for heat and mass transfer in ACs owing to its large vapor pressure and enthalpy of evaporation. Zeolite imidazolate frameworks (ZIFs) with excellent ammonia stability are identified as a class of potential adsorbents for practical ammonia-based ACs. However, high-performing ZIF/ammonia working pairs with excellent AC performance are still to be developed. In this work, the cooling performance including the coefficient of performance for cooling (COPC) and the specific cooling effects (SCEs) of 26 ZIFs with the same composites but different topologies was evaluated by combining molecular simulation and mathematical modeling. Five high-performing ZIFs with COPC > 0.45 and SCE > 250 kJ/kg were identified, among which gis-ZIF with the highest COPC of 0.51 and lta-ZIF with the highest SCE of 354 kJ/kg both are promising to be synthesized and applied further. Besides, the quantitative structure-performance relationship (QSPR) was extracted that can help quickly identify and design high-performing ZIFs according to their ammonia adsorption isotherms and structural characteristics. Moreover, "S"-shaped adsorption isotherms with high saturation adsorption capacity (>0.2 g/g), suitable step position (0.2-0.4), and relatively low Henry's constant (<1 × 10-5 mol/(kg·Pa)) are more favorable for excellent COPC and SCE. From the perspective of structure characteristics, ZIFs possessing low crystal density (<0.9 g/cm3), high accessible surface area (>2000 m2/g), balanced largest cavity diameter (∼15 Å), and accessible pore volume (∼0.65 cm3/g) are beneficial for high-efficient cooling performance.

In Silico Screening of Zeolites for High-Pressure Hydrogen Drying

According to the ISO 14687-2:2019 standard, the water content of H₂ fuel for transportation and stationary applications should not exceed 5 ppm (molar). To achieve this water content, zeolites can be used as a selective adsorbent for water. In this work, a computational screening study is carried out for the first time to identify potential zeolite frameworks for the drying of high-pressure H₂ gas using Monte Carlo (MC) simulations. We show that the Si/Al ratio and adsorption selectivity have a negative correlation. 218 zeolites available in the database of the International Zeolite Association are considered in the screening. We computed the adsorption selectivity of each zeolite for water from the high-pressure H₂ gas having water content relevant to vehicular applications and near saturation. It is shown that due to the formation of water clusters, the water content in the H₂ gas has a significant effect on the selectivity of zeolites with a helium void fraction larger than 0.1. Under each operating condition, five most promising zeolites are identified based on the adsorption selectivity, the pore limiting diameter, and the volume of H₂ gas that can be dried by 1 dm³ of zeolite. It is shown that at 12.3 ppm (molar) water content, structures with helium void fractions smaller than 0.07 are preferred. The structures identified for 478 ppm (molar) water content have helium void fractions larger than 0.26. The proposed zeolites can be used to dry 400-8000 times their own volume of H₂ gas depending on the operating conditions. Our findings strongly indicate that zeolites are potential candidates for the drying of high-pressure H₂ gas.

Screening metal-organic frameworks for adsorption-driven osmotic heat engines via grand canonical Monte Carlo simulations and machine learning

Adsorption-driven osmotic heat engines offer an alternative way for harvesting low-grade waste heat below 80°C. In this study, we performed a high-throughput computational screening based on grand canonical Monte Carlo simulations to identify the high-performance metal-organic frameworks (MOFs) from 1322 computationally ready experimental MOF structures for adsorption-driven osmotic heat engines with LiCl-methanol as the working fluid. Structure-property relationship analysis reveals that MOFs exhibiting high energy efficiency possess large working capacity, pore size and surface area, and moderate adsorption enthalpy comparable to the evaporation enthalpy. Furthermore, machine learning is employed to accelerate the computational screening for satisfied MOFs via the structure properties. The optimal structure properties of the MOFs are further identified via the ensemble-based regression model by optimizing the energy efficiency via the genetic algorithm, which shed light on rationally designing and fabricating MOFs for desired heat-to-electricity conversion.

Too Many Materials and Too Many Applications: An Experimental Problem Waiting for a Computational Solution

Finding the best material for a specific application is the ultimate goal of materials discovery. However, there is also the reverse problem: when experimental groups discover a new material, they would like to know all the possible applications this material would be promising for. Computational modeling can aim to fulfill this expectation, thanks to the sustained growth of computing power and the collective engagement of the scientific community in developing more efficient and accurate workflows for predicting materials' performances. We discuss the impact that reproducibility and automation of the modeling protocols have on the field of gas adsorption in nanoporous crystals. We envision a platform that combines these tools and enables effective matching between promising materials and industrial applications.

Adsorber heat exchanger using Al-fumarate beads for heat-pump applications - a transport study

Metal-Organic Frameworks (MOFs), thanks to their type V water adsorption isotherms ("S-Shape") and large water capacities, are considered as potential breakthrough adsorbents for heat-pump applications. In particular, Al(OH)-fumarate could enable efficient regeneration at a lower temperature than silica-gel which would allow us to address the conversion of waste heat at low temperature such as found in data centers. Despite its greater adsorption capacity features, heat and mass transport limitations could jeopardize the potential performance of Al(OH)-fumarate. Heat and mass transport depend on the size of the bodies (mm range), their packing and on the pore structures, i.e. macro-mesopore volumes and sizes. This paper describes the cost-efficient and scalable synthesis and shaping processes of Al(OH)-fumarate beads of various sizes appropriate for use in water Adsorption Heat-Pumps (AHPs). The objective was to study transport limitations (i.e. mass and heat) in practical e beads which meet mechanical stability requirements. Dynamic data at the grain scale was obtained by the Large Temperature Jump method while dynamic data at the adsorber scale was obtained on a heat exchanger filled with more than 1 kg of Al(OH)-fumarate beads. Whereas the binder content had little impact on mass and heat transfer in this study, we found that Knudsen diffusion in mesopores of the grain may be the main limiting factor at the grain scale. At the adsorber scale, heat-transfer within the bed packing as well as to the heat exchanger is likely responsible for the slow adsorption and desorption kinetics which have been observed for very low desorption temperature. Finally, the dynamic aspects of the observed water adsorption isotherm shift with temperature are discussed in light of reported reversible structure modification upon temperature triggered water adsorption-desorption.

Screening of metal-organic frameworks for water adsorption heat transformation using structure-property relationships

It is of great importance to correlate the water adsorption performance of MOFs to their physicochemical features in order to design and prepare MOFs for applications in adsorption heat transformation. In this work, both data analysis from existing studies and Grand Canonical Monte Carlo molecular simulation investigations were carried out. The results indicated that the highest water adsorption capacity was determined by the pore volume of MOF adsorbents, while there was a linear correlation interrelationship between isosteric heats of adsorption and the water adsorption performance at a low relative pressure. More detailed analysis showed that the charge distribution framework and pore size of MOFs contributed together to the hydrophilicity. Electrostatic interaction between water molecules and the framework atoms played a key role at low relative water pressure. A quantitative structure-property relationship model that can correlate the hydrophilicity of MOFs to their pore size and atomic partial charge was established. Along with some qualitative considerations, the screening methodology is proposed and is used to screen proper MOFs in the CoRE database. Seven MOFs were detected, and four of them were synthesized to validate the screening principle. The results indicated that these four MOFs possessed outstanding water adsorption performance and could be considered as promising candidates in applications for adsorption heating and cooling.

Evaluation of Metal–Organic Frameworks as Potential Adsorbents for Solar Cooling Applications

The reduction of carbon dioxide emissions has become a need of the day to overcome different environmental issues and challenges. The use of alternative and renewable-based technologies is one of the options to achieve the target of sustainable development through the reduction of these harmful emissions. Among different technologies thermally activated cooling systems are one which can reduce the harmful emissions caused by conventional heating, ventilation, and air conditioning technology. Thermal cooling systems utilize different porous materials and work on a reversible adsorption/desorption cycle. Different advancements have been made for this technology but still a lot of work should be done to replace conventional systems with this newly developed technology. High adsorption capacity and lower input heat are two major requirements for efficient thermally driven cooling technologies. In this regard, it is a need of the day to develop novel adsorbents with high sorption capacity and low regeneration temperature. Due to tunable topologies and a highly porous nature, the hybrid porous crystalline materials known as metal–organic frameworks (MOFs) are a great inspiration for thermally driven adsorption-based cooling applications. Keeping all the above-mentioned aspects in mind, this paper presents a comprehensive overview of the potential use of MOFs as adsorbent material for adsorption and desiccant cooling technologies. A detailed overview of MOFs, their structure, and their stability are presented. This review will be helpful for the research community to have updated research progress in MOFs and their potential use for adsorption-based cooling systems.

Gibbs Ensemble Monte Carlo Simulation of Fluids in Confinement: Relation between the Differential and Integral Pressures

The accurate description of the behavior of fluids in nanoporous materials is of great importance for numerous industrial applications. Recently, a new approach was reported to calculate the pressure of nanoconfined fluids. In this approach, two different pressures are defined to take into account the smallness of the system: the so-called differential and the integral pressures. Here, the effect of several factors contributing to the confinement of fluids in nanopores are investigated using the definitions of the differential and integral pressures. Monte Carlo (MC) simulations are performed in a variation of the Gibbs ensemble to study the effect of the pore geometry, fluid-wall interactions, and differential pressure of the bulk fluid phase. It is shown that the differential and integral pressure are different for small pores and become equal as the pore size increases. The ratio of the driving forces for mass transport in the bulk and in the confined fluid is also studied. It is found that, for small pore sizes (i.e., < 5 σ fluid ), the ratio of the two driving forces considerably deviates from 1.

Screening of Covalent-Organic Frameworks for Adsorption Heat Pumps

Exploring high-performing adsorption-driven heat pumps (AHPs) remains a challenging task owing to the low working capacity, high regeneration temperature, and low energy efficiency of conventional adsorbents. Quick discovery of the novel promising adsorbents could help to improve the coefficient of performance of AHPs for heating (COP_H) and cooling (COP_C). Herein, we reported an approach to identify the high-performing covalent-organic frameworks (COFs) for heating, cooling, and ice making by high-throughput computational screening based on grand canonical Monte Carlo simulations and, for the first time, machine learning. It was demonstrated that compared with metal-organic frameworks (MOFs), COFs were more suitable adsorbents of AHPs for cooling because of their weak interaction toward ethanol that favors stepwise adsorption. Structure-property relationship analysis revealed that the average enthalpy of adsorption commensurate with the enthalpy of evaporation will benefit the performance of AHPs besides the high working capacity and low step positions of adsorption isotherms. In order to reduce the computational cost of screening, a random forest model was developed to successfully predict the COP_C of both COFs and MOFs.

Silica-Supported Ionic Liquids for Heat-Powered Sorption Desalination

This work investigates the application of novel sorption materials to heat-powered desalination systems. Two ionic liquids 1-ethyl-3-methylimidazolium acetate (Emim-Ac) and 1-ethyl-3-methylimidazolium methanesulfonate (Emim-Oms) were impregnated in two silica supports, namely, Syloid AL-1FP and Syloid 72FP. Emim-Ac and Emim-Oms composite sorbents have been compared on morphology, water vapor sorption equilibrium, and heat of sorption. Fourier-transform infrared spectroscopy shows that the ionic liquid partly self-organizes on the silica surface. When used in a sorption desalination process powered by low grade heat at 60 °C, these composites have exceptionally high theoretical working capacities ranging from 1 to 1.7 g_water g_sorbent^-1. Experimental tests on a lab scale desalinator show that Emim-Ac/Syloid 72FP in real operating conditions can produce 25 kg_water kg_sorbent^-1 day^-1. To date, this yield is 2.5 times higher than the best achieved with silica gel.

On the Efficient Separation of Gas Mixtures with the Mixed-Linker Zeolitic-Imidazolate Framework-7-8

A recently reported modification of the zeolitic-imidazolate framework-8 (ZIF-8) with partial replacement of the 2-methylimidazolate (mIm) linker with benzimidazolate (bIm), namely ZIF-7-8, is investigated with molecular simulations using a first-time reported force field. The size of the ZIF-7-8 aperture, which governs the gas-separation efficiency of this material and which has not been estimated before for this modification, is smaller than that of the original ZIF-8. The diffusivities of CO₂, N₂, and CH₄ estimated through transition state theory calculations result in remarkably high diffusion selectivities for CO₂/CH₄ and CO₂/N₂ mixtures. This performance enhancement is investigated in terms of structural flexibility in the form of the aperture motion through extensive estimation of the effective diameter, the total effective area, and the motion of the aperture linkers, of both ZIF-8 and ZIF-7-8. Both apertures exhibit an oscillation through the rotation of the linkers, which are adjusted according to the size of the penetrant molecules the moment they pass through it. Finally, a subsequent analysis reveals that there is strong dependency of the separation performance on the bIm-to-mIm ratio: below 33% bIm incorporation, the appearance of ZIF-8-alike wide apertures decreases dramatically the size-based selectivity of the mixtures in ZIF-7-8.

Polarizable Force Field for CO2 in M-MOF-74 Derived from Quantum Mechanics

On the short term, carbon capture is a viable solution to reduce human-induced CO₂ emissions, which requires an energy efficient separation of CO₂. Metal-organic frameworks (MOFs) may offer opportunities for carbon capture and other industrially relevant separations. Especially, MOFs with embedded open metal sites have been shown to be promising. Molecular simulation is a useful tool to predict the performance of MOFs even before the synthesis of the material. This reduces the experimental effort, and the selection process of the most suitable MOF for a particular application can be accelerated. To describe the interactions between open metal sites and guest molecules in molecular simulation is challenging. Polarizable force fields have potential to improve the description of such specific interactions. Previously, we tested the applicability of polarizable force fields for CO₂ in M-MOF-74 by verifying the ability to reproduce experimental measurements. Here, we develop a predictive polarizable force field for CO₂ in M-MOF-74 (M = Co, Fe, Mg, Mn, Ni, Zn) without the requirement of experimental data. The force field is derived from energies predicted from quantum mechanics. The procedure is easily transferable to other MOFs. To incorporate explicit polarization, the induced dipole method is applied between the framework and the guest molecule. Atomic polarizabilities are assigned according to the literature. Only the Lennard-Jones parameters of the open metal sites are parameterized to reproduce energies from quantum mechanics. The created polarizable force field for CO₂ in M-MOF-74 can describe the adsorption well and even better than that in our previous work.