Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “adsorption isotherm”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Discrepancy quantification between experimental and simulated data of CO 2 adsorption isotherm using hierarchical Bayesian estimation

Here, to quantitatively analyze the inconsistencies commonly observed between experimental and simulated adsorption isotherms, parameter estimation of adsorption isotherm models was conducted by hierarchical Bayesian estimation with parameter uncertainties being quantified as probability distributions. The estimation method was implemented using Markov Chain Monte Carlo (MCMC) to analyze multiple data sets obtained from different sources, including a publicly available database. To describe the discrepancies of experimental and simulated adsorption data, the simulation data was set as the reference to which experimental measurements were compared. We applied the proposed approach to analyze CO 2 adsorption isotherms that are measured and simulated on zeolite 13X and MIL-101(Cr). In these case studies, the discrepancy of CO 2 adsorption isotherm was successfully quantified between experimental measurements and predictions given by molecular simulations using Grand Canonical Monte Carlo (GCMC), where uncertainties were quantified as probability distributions. Furthermore, experimental data sets that agree well with the GCMC simulation have been identified, providing insights into experimental and measurement methods as well as choosing the right assumptions in the molecular simulation.

42 ENGINEERING↗
DOE OSTI2022

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

High-Temperature Water Adsorption Isotherms and Ambient Temperature Water Diffusion Rates on Water Harvesting Metal–Organic Frameworks

Water adsorption isotherms from 25 to 125 °C were measured for three metal−organic frameworks (MOFs), MOF-303, MOF-LA2-1, and MIL-100(Fe), which are frequently studied for water harvesting applications. The results show how the step in the water adsorption isotherm varies as a function of temperature and detail the combination of pressure and temperature necessary to remove adsorbed water. Furthermore, isobaric−isothermal Gibbs ensemble Monte Carlo simulations performed for MOF-303 shed light on the change in occupation numbers of the different known water adsorption sites with increasing temperature. Additionally, the diffusion rates of water through these materials were measured using concentration swing frequency response, and micropore diffusion was identified as the controlling mechanism. The Darken relation was used to show the dependence of the diffusion rate on the concentration and the impact of the adsorption isotherm slope. The adsorption of water on MOF-LA2-1 is faster than that on MIL-100(Fe). These data show that MOF-LA2-1 with its high-water adsorption capacity, quick adsorption rate, and favorable desorption energetics is a leading candidate for atmospheric water harvesting.

Adsorption↗
DOE OSTI2025

Understanding Differences in Water Adsorption Isotherms: Structural Variations, Force Fields, and Monte Carlo Simulation Approaches

Accurate prediction of water adsorption in micro- and mesoporous materials with hydrophobic pores is essential for the design and characterization of advanced adsorbent materials for separation and energy applications. Here, we assess the reproducibility and consistency of water adsorption isotherms in two microporous all-silica MFI zeolite structures (MFI-K and MFI-O) using two different zeolite force fields and three simulation approaches: grand canonical Monte Carlo (GCMC), Gibbs ensemble Monte Carlo (GEMC), and transition matrix Monte Carlo (TMMC). We demonstrate that consistent treatment of the bulk fluid phase in GCMC and TMMC simulations is critical for reconciling isotherms across methods, and we construct simulation-based equations of state for the TIP4P water model to enable rigorous fugacity-to-pressure conversions. Large shifts in the isotherms are observed for two zeolite force fields developed using different parametrization strategies, with the GCS force field representing implicitly a defect-containing all-silica zeolite, whereas the TraPPE-zeo force field accurately represents an essentially defect-free all-silica zeolite. While water in the van Koningsveld structure of MFI exhibits a first-order phase transition and condensation-like step for adsorption near room temperature, water in the Olson structure of MFI displays continuous adsorption, attributed to differences in the adsorption free energy landscapes. Structural analysis reveals that small geometric variations, particularly Si–O–Si bond angles near the strongest adsorption sites, lead to these substantial differences in adsorption behavior. Furthermore, our results highlight the sensitivity of simulated water adsorption isotherms in hydrophobic frameworks to seemingly small differences in the framework structures, force field parametrization, and simulation approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The overlooked role of adsorption isotherms in electrocatalysis

Electrocatalysts enable the efficient interconversion of electrical and chemical energy for the sustainable production of fuels and chemicals. Here, in this Comment, we highlight the importance of developing electrochemical adsorption isotherms to demystify complex reaction mechanisms and rationalize catalytic activity.

electrocatalysis↗
DOE OSTI2022

Experimentally Verified Alkane Adsorption Isotherms in Nanoporous Materials from Literature Meta-Analysis

Identification of independent experimental replicate measurements of well-defined physical properties provides a useful way to assess the reproducibility of literature reports of these properties. Prior meta-analysis of this kind for single-component adsorption of CO 2 and alcohols in nanoporous materials showed that ~20% of reported isotherms for these molecules were inconsistent with available replicate measurements. Here, we extend this analysis to the single-component adsorption of alkanes in nanoporous materials using a collection of 2998 experimental isotherms. Among these isotherms, 846 replicates were identified and analyzed. Many more replicates are available for CH 4 adsorption than for any other alkane. Our analysis finds that 15% of the replicate alkane isotherms are inconsistent with other replicate measurements. Furthermore, after discarding these outlier isotherms we report graphical and numerical consensus isotherms for 93 independent combinations of adsorbate, adsorbent, and temperature as well as an assessment of another 124 independent systems for which two replicate isotherms are available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Elucidating the Influence of Metal Surface Composition on Organic Adsorbate Binding Using Active Particle Dynamics

The adsorption strengths of organic compounds on metal surfaces are sensitive to the metal composition, and they play a central role in many catalytic reactions, helping to control the coverage of the reactant and altering the overall reaction rate. While adsorption energies are straightforward to measure and calculate in vacuum and gas-phase environments, adsorption energetics can be dramatically altered by the presence of solvent in liquid-phase reactions. However, the effects of metal composition on binding strengths in a liquid environment are less well understood, primarily due to the difficulty of accurate in situ measurements of organic binding on metal surfaces in the liquid phase. Here, we utilize the motion of active particles in water to probe the adsorption energies of an organic adsorbate (furfural) on a range of metal surfaces (pure Pd, pure Pt, and four PdAu alloy compositions) to elucidate the effect of metal composition. Janus particles with catalytic caps of particular metal compositions all exhibited active motion resulting from consumption of H 2 O 2 ; adsorbate binding was inferred through the decrease in the velocity of active motion and was modeled by a Langmuir adsorption isotherm. The measured adsorption affinities were used to extract the adsorption enthalpy of furfural on the different metals. The Pd surface was found to bind furfural more strongly than the Pt surface by some 10 kJ/mol. Furthermore, the adsorption of furfural on the alloys was found to increase monotonically in magnitude with Pd content. As a result, the data reported herein aid the development of accurate understanding of organic adsorption in the presence of solvent and the role of the metal surface in tuning adsorption strengths to optimize catalytic processes in the liquid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thiol-Functionalized Adsorbents through Atomic Layer Deposition and Vapor-Phase Silanization for Heavy Metal Ion Removal

The removal of toxic heavy metal ions from water resources is crucial for environmental protection and public health. In this study, we address this challenge by developing a surface functionalization technique for the selective adsorption of these contaminants. Our approach involves atomic layer deposition (ALD) followed by vapor-phase silanization of porous substrates. We utilized porous silica gel powder (similar to 100 mu m particles, 89 m 2 /g surface area, similar to 30 nm pores) as an initial substrate. This powder was first coated with similar to 0.5 nm ALD Al 2 O 3 , followed by vapor-phase grafting of a thiol-functional silane. The modified powder, particularly in acidic conditions (pH = 4), showed high selectivity in adsorbing Cd(II), As(V), Pb(II), Hg(II), and Cu(II) heavy metal ions in mixed ion solutions over common benign ions (e.g., Na, K, Ca, and Mg). Langmuir adsorption isotherms and breakthrough adsorption studies were conducted to assess heavy metal binding affinity and revealed the order of Cd(II) < Pb(II) < Cu(II) < As(V) < Hg(II), with a significantly higher affinity for As(V) and Hg(II) ions. Time-dependent uptake studies demonstrated rapid removal of heavy metal ions from aqueous environments, with Hg(II) exhibiting the fastest adsorption kinetics on thiol-modified surfaces. Finally, these findings highlight the potential of ALD and vapor-phase silanization to create effective adsorbents for the targeted removal of hazardous contaminants from water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Liquid-Phase Effects on Adsorption Processes in Heterogeneous Catalysis

Aqueous solvation free energies of adsorption have recently been measured for phenol adsorption on Pt(111). Endergonic solvent effects of ~1 eV suggest solvents dramatically influence a metal catalyst’s activity with significant implications for the catalyst design. However, measurements are indirect and involve adsorption isotherm models, which potentially reduces the reliability of the extracted energy values. Computational, implicit solvation models predict exergonic solvation effects for phenol adsorption, failing to agree with measurements even qualitatively. In this study, an explicit, hybrid quantum mechanical/molecular mechanical approach for computing solvation free energies of adsorption is developed, solvation free energies of phenol adsorption are computed, and experimental data for solvation free energies of phenol adsorption are reanalyzed using multiple adsorption isotherm models. Explicit solvation calculations predict an endergonic solvation free energy for phenol adsorption that agrees well with measurements to within the experimental and force field uncertainties. Computed adsorption free energies of solvation of carbon monoxide, ethylene glycol, benzene, and phenol over the (111) facet of Pt and Cu suggest that liquid water destabilizes all adsorbed species, with the largest impact on the largest adsorbates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Investigation of the Effect of Framework Flexibility on Adsorption in SIFSIX-3-Cu using a Machine-Learned Force Field

Metal-organic frameworks (MOFs) are a promising class of adsorbents. The performance of MOF sorbents relies on high selectivity and low regeneration energy. This work focuses on the use of machine learned force fields (MLFFs) to model adsorption in a flexible MOF, SIFSIX-3-Cu. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies, forces, and stresses, using an iterative sampling scheme combining sampling based on molecular dynamics, Monte Carlo, and geometry optimization to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included in this model. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations using the MLFF predicted adsorption isotherms in good agreement with experimental data for a range of pressures (40 Pa – 104 Pa) in contrast to rigid models, which overpredict CO2 adsorption at low pressures. The improvement was the result of a description of the variability of fluorine-fluorine diagonal distances at adsorption sites. This detailed description of flexibility afforded by the MLFF resulted in more accurate predictions adsorption isotherms when compared to the experimentally measured values. These results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low pressure applications and provide a robust procedure for training MLFF models for MOFs.

Atomistic Simulation↗
DOE OSTI2023

Analysis of gas storage and transport in Eagle Ford shale using pressure pulse decay measurements with He, Kr and CO 2

Detailed characterization of gas transport and storage in low-permeability and organic-rich shales is associated with an array of challenges due to the complex morphology of the pore space, representing a broad range of pore sizes, combined with the heterogeneous fabric of the shale matrix. These factors, and their interplay during gas transport and sorption, complicate a) the analysis of shale samples at the laboratory scale (~ up to a ft. In length), and b) the prediction of natural gas production and carbon sequestration potential at larger scales. Here, in this work, tandem experiments with inert (helium – He) and adsorbing (krypton - Kr and carbon dioxide - CO 2 ) gases were performed and analyzed to develop an efficient workflow for characterizing and modeling transport and sorption at the laboratory scale. In particular, pressure pulse-decay (PPD) measurements were conducted on an Eagle Ford shale core sample at room temperature using He, Kr, and CO 2 . PPD measurements with He (a non-sorbing gas) were employed to probe the overall porosity, including natural fractures, microcracks, mesopores, and micropores. A triple-porosity model (TPM) was adopted to interpret the gas transport in the shale sample: The pore space is represented by three interacting continua, including macropores (larger fractures), mesopores (including microcracks), and micropores. To facilitate the application of the TPM, a modified analytical approach is introduced to extract effective transport parameters (in terms of the characteristic time for transport at relevant porosity levels) directly from the PPD measurements with He. To validate the modeling approach, model parameters extracted from two He PPD measurements are demonstrated to provide excellent agreement with a 3rd He PPD measurement performed at a higher pressure. The effective transport parameters, extracted from the He experiments, were subsequently converted for application to the Kr and CO 2 PPD experiments, by accounting for relevant transport modes and differences in fluid properties. Excess adsorption isotherms were extracted from the equilibrium pressures of the Kr and CO 2 PPD experiments using the He pore volume as a baseline. These adsorption isotherms were then integrated into the TPM to predict combined gas transport and sorption for Kr and CO 2 , with transport coefficients translated from the He measurements. The predictions for Kr and CO 2 are demonstrated to be in excellent agreement with the experimental observations. This, in turn, demonstrates that the proposed analytical approach provides for an effective characterization of mass transfer rates in shales, that can be applied directly in a TPM representation of mass transfer and sorption.

58 GEOSCIENCES↗
DOE OSTI2023

Molecular Dynamics Simulation of Pore-Size Effects on Gas Adsorption Kinetics in Zeolites

Strong gas-mineral interactions or slow adsorption kinetics require a molecular-level understanding of both adsorption and diffusion for these interactions to be properly described in transport models. In this combined molecular simulation and experimental study, noble gas adsorption and mobility is investigated in two naturally abundant zeolites whose pores are similar in size (clinoptilolite) and greater than (mordenite) the gas diameters. Simulated adsorption isotherms obtained from grand canonical Monte Carlo simulations indicate that both zeolites can accommodate even the largest gas (Rn). However, gas mobility in clinoptilolite is significantly hindered at pore-limiting window sites, as seen from molecular dynamics simulations in both bulk and slab zeolite models. Further, experimental gas adsorption isotherms for clinoptilolite confirm the presence of a kinetic barrier to Xe uptake, resulting in the unusual property of reverse Kr/Xe selectivity. Finally, a kinetic model is used to fit the simulated gas loading profiles, allowing a comparison of trends in gas diffusivity in the zeolite pores.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Insights into Controls of Mineralogy and Pore Structure on the Density of Methane Adsorption Phase in Shales under Supercritical Conditions

With gradually increasing enhancing population and industrialization, exploitation of shale gas is enhancing in the world to satisfy the growing demand for energy worldwide. To exploit shale gas efficiently, accurate assessment of shale gas in place (GIP) is necessary for determining production strategies. Unlike free and dissolved gas, adsorbed gas contributes to the shale GIP up to 85% due to the well-developed surface area and micropores. Measured excess adsorption of methane (up to over 95% of subsurface shale gas) is often corrected to absolute adsorption to obtain the actual amount of adsorption under supercritical (or geological) conditions. During the correction, adsorption phase density (APD) is critical. However, the APD of subsurface shale gas and the effects of shale properties (e.g., mineralogy and pore structure) on APD remain poorly understood. A series of high-pressure methane adsorption isotherms (HPMAI) on Caney Shales were collected and analyzed in conjunction with other United States and Chinese shales from the literature at 35–125 °C and up to 15 MPa. Here, a three-layer Ono–Kondo (OK3) model is utilized to derive the temperature- and pressure-dependent APD coupling low-pressure nitrogen adsorption isotherms (LPNAI) and HPMAI. X-ray diffraction and organic geochemistry are combined to reveal the mineralogy. Brunauer–Emmett–Teller, Barrett–Joyner–Halenda, and Horvath–Kawazoe analyses via LPNAI are used to investigate the pore structures. Results show that APD increases with organic matter (OM) proxied by total organic carbon and decreases with clay minerals and the sum of quartz and feldspar. OM dramatically contributes to the APD as multiple-layer adsorption exists, and the APD for OM could be 1.4–8.5 times that for clay minerals. Other inorganic minerals contribute less to APD. The properties and constitution of the surface area instead of the volume fraction contribute to the APD in shales. APD does not show an obvious correlation with micropore volume, likely related to the ratio of micropore volume to the total pore volume. Here, we provide a significant and comprehensive study of petrological factors that impact the APD of subsurface shale gas, which will improve the estimation of supercritical adsorption and shale GIP under reservoir conditions. Also, the findings in this work can provide applications for subsurface carbon dioxide adsorption and storage.

58 GEOSCIENCES↗
DOE OSTI2024

Exploitation of Pore Structure for Increased CO 2 Selectivity in Type 3 Porous Liquids

CO 2 capture requires materials with high adsorption selectivity and an industrial ease of implementation. To address these needs, a new class of porous materials was recently developed that combines the fluidity of solvents with the porosity of solids. Type 3 porous liquids (PLs) composed of solvents and metal–organic frameworks (MOFs) offer a promising alternative to current liquid carbon capture methods due to the inherent tunability of the nanoporous MOFs. However, the effects of MOF structural features and solvent properties on CO 2 –MOF interactions within PLs are not well understood. Herein experimental and computational data of CO 2 gas adsorption isotherms were used to elucidate both solvent and pore structure influences on ZIF-based PLs. The roles of the pore structure including solvent size exclusion, structural environment, and MOF porosity on PL CO 2 uptake were examined. A comparison of the pore structure and pore aperture was performed using ZIF-8, ZIF-L, and amorphous-ZIF-8. Adsorption experiments here have verified our previously proposed solvent size design principle for ZIF-based PLs (1.8× ZIF pore aperture). Furthermore, the CO 2 adsorption isotherms of the ZIF-based PLs indicated that judicious selection of the pore environment allows for an increase in CO 2 selectivity greater than expected from the individual PL components or their combination. This nonlinear increase in the CO 2 selectivity is an emergent behavior resulting from the complex mixture of components specific to the ZIF-L + 2'-hydroxyacetophenone-based PL.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Effect of layer bending on montmorillonite hydration and structure from molecular simulation

Conceptual models of smectite hydration include planar (flat) clay layers that undergo stepwise expansion as successive monolayers of water molecules fill the interlayer regions. However, X-ray diffraction (XRD) studies indicate the presence of interstratified hydration states, suggesting non-uniform interlayer hydration in smectites. Additionally, recent theoretical studies have shown that clay layers can adopt bent configurations over nanometer-scale lateral dimensions with minimal effect on mechanical properties. Therefore, in this study we used molecular simulations to evaluate structural properties and water adsorption isotherms for montmorillonite models composed of bent clay layers in mixed hydration states. Results are compared with models consisting of planar clay layers with interstratified hydration states (e.g. 1W–2W). The small degree of bending in these models (up to 1.5 Å of vertical displacement over a 1.3 nm lateral dimension) had little or no effect on bond lengths and angle distributions within the clay layers. Except for models that included dry states, porosities and simulated water adsorption isotherms were nearly identical for bent or flat clay layers with the same averaged layer spacing. Similar agreement was seen with Na- and Ca-exchanged clays. In conclusion, while the small bent models did not retain their configurations during unconstrained molecular dynamics simulation with flexible clay layers, we show that bent structures are stable at much larger length scales by simulating a 41.6×7.1 nm 2 system that included dehydrated and hydrated regions in the same interlayer.

58 GEOSCIENCES↗
DOE OSTI2021

Manipulating Pore Topology and Functionality to Promote Fluorocarbon-Based Adsorption Cooling

With the worldwide demand for refrigeration and cooling expected to triple, it is increasingly important to search for alternative energy resources to drive the refrigeration cycles with reduced electricity consumption. Recently, adsorption cooling has gained increased attention since energy reallocation in such systems is based on gas adsorption/desorption, which can be driven by waste/natural heat sources. Eco-friendly sorption-based cooling relies on the cyclic transfer of refrigerant gas from a high to low energy state by the pseudo-compression effect resulting from adsorption and desorption. The driving force for energy transfer relies on heat rather than electricity. The performance of a sorption chiller is primarily influenced by this cyclic sorption behavior, which is characterized as the working capacity of the porous sorbent. Thus increases in this working capacity directly translate to a more compact and efficient cooling system. How-ever, a lack of highly effective sorbent/refrigerant pairs lowers cooling performance, and therefore has limited applicability. To this end, synthetic metal-organic frameworks (MOFs) and covalent organic polymers (COPs) possess higher porosity and greater tunability leading to more substantial potential benefits for adsorption, compared to traditional sorbent materials. Similarly, hydrofluorocarbon refrigerants have more favorable applicability given ease of operation above atmospheric pressures due to suitable saturated vapor pressures and boiling points. For these reasons, our work focuses on an ongoing strategy to promote sorption cooling via improvements in the sorbent/refrigerant pair. Specifically, we target the interaction of hydrofluorocarbon refrigerants with MOF/COP materials at a molecular level by interpreting the host-guest chemistry and the role of framework pore topology. These molecular level differences translate to cooling performance, which is described herein. These strategies include engineering framework porosity (i.e., pore size, pore volume) by using elongated organic linkers and stereochemistry control during synthesis; manipulating the sorb-ate/sorbent interaction by introducing functional moieties or unsaturated metal centers to enhance working capacities in narrow pressure ranges; varying pore topology/morphology to impact adsorption isotherm behavior; and leveraging defective sites within the frameworks to further enhance adsorption capability. Here this atomic level understanding of sorb-ate-sorbent interactions is conducted using various in situ experimental techniques such as synchrotron-based X-ray diffraction, X-ray absorption spectroscopy, in situ Fourier transform infrared spectroscopy, and direct sorption energies determinization with calorimetry. Moreover, the experimentally studied interactions and the corresponding adsorption mechanism are corroborated by computational studies using density functional theory (DFT) and grand canonical Monte Carlo (GCMC) simulations. Using this approach, we envision the capability to engineer designed frameworks with precise molecular control to target refrigerant molecules and thereby enhance the performance of desired working pairs for sorption-based cooling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗