Search NASASearch

SEARCH · Search NASA

Results for “adsorption”

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

Modeling Pb(II) Adsorption on Mineral Surfaces: Bridging Density Functional Theory and Experiment with Thermodynamic Insights

Despite decades of work on aqueous lead (Pb) adsorption on a-Fe2O3 (hematite) and a-Al2O3 (alumina), gaps between measurements and modeling obscure molecular-level understanding. Achieving well-matched geometries between theory and experimental for mineral-water interfaces is a hurdle, as surface functional group type and distribution must be accounted for in determining mechanisms. Additionally, computational methods that can describe the substrate are often not appropriate to capture aqueous effects. Progress requires focusing on well-studied and relevant systems, such as key facets (001),(012), and (110) of hematite and alumina, and ubiquitous contaminants such as aqueous Pb. In the past, bulk-parametrized bond-valence principles were used to rationalize Pb(II) adsorption trends. These approaches can break down at surfaces, where flexible bonding environments and adsorption-induced surface relaxations play a critical role. Here, we adapt and apply a density functional theory (DFT) and thermodynamics framework, integrating DFT-calculated energies with experimental data and electrochemical principles, to predict Pb(II) adsorption. Our model results capture trends across the full set of surfaces and predict that inner-sphere Pb(II) sorption on (001) alumina varies from unfavorable to weakly favorable across a range of pH conditions. This aligns with experiment insights that Pb(II) interacts at that surface through outer-sphere interactions. Extending to Fe(II) adsorption, we demonstrate a coverage-dependent site preference, potentially explaining disorder in overlayers grown by the oxidative adsorption of Fe(II) on hematite (001).

lead contamination

Preferential adsorption of noble gases in zeolitic tuff with variable saturation: A modeling study of counter-intuitive diffusive-adsorptive behavior

Noble gas transport through geologic media has important applications in the prediction and characterization of measured gas signatures related to underground nuclear explosions (UNEs). Retarding processes such as adsorption can cause significant species fractionation of radionuclide gases, which has implications for measured and predicted signatures used to distinguish radioxenon originating from civilian nuclear facilities or from UNEs. Accounting for the effects of variable water saturation in geologic media on tracer transport is one of the most challenging aspects of modeling gas transport because there is no unifying relationship for the associated tortuosity changes between different rock types, and reactive transport processes such as adsorption that are affected by the presence of water likewise behave differently between gas species. In this study, we perform numerical diffusive-adsorptive transport simulations to estimate gas transport parameters associated with bench-scale laboratory diffusion cell experiments measuring breakthrough in zeolitic and non-zeolitic rocks for a gaseous mixture of xenon, krypton, and SF 6 at varying degrees of water saturation (S w ). Counter-intuitive transport behavior was observed in the zeolitic rock experiments whereby breakthrough concentrations were significantly higher when the core was partially saturated (S w = 17 %) than under dry (S w = 0 %) conditions. Breakthrough of xenon was especially retarded in the dry core – likely due to comparatively high affinity of xenon for zeolitic adsorption sites – and estimated effective diffusion coefficients for all gases were approximately an order of magnitude lower than what is predicted by porosity-tortuosity models. We propose the counter-intuitive behavior observed is because water infiltration of zeolite nanopores reduces both the adsorptive capacity of the rock and the tortuosity of connected flow paths. We developed a two-site competitive kinetic Langmuir adsorption reaction for the porous media transport simulator in order to constrain transport parameters within zeolitic tuff, where differential adsorption to zeolite and non-zeolite pores was observed. We determined that liquid saturation-dependent diffusive-adsorptive transport is affected by subtle and at times competing processes that are specific to different gases, which have a significant overall influence on effective transport parameters.

58 GEOSCIENCES

Adsorptive denitrogenation of model aviation fuel using mesoporous silica in a packed bed adsorption system

This study aims to understand the effects of system process parameters such as flow rate, adsorbent particle size, and use of recycled adsorbent on denitrogenation performance of a model fuel using mesoporous silica gel. The goal is to reduce the nitrogen content of the model fuel from 1500 parts per million to single-digit ppm to meet ASTM specifications for drop-in fuels. This work was done with the intent of applying adsorptive denitrogenation to sustainable aviation fuel (SAF) product fractions produced via hydrothermal liquefaction (HTL). The adsorption performance of the silica is evaluated via packed column breakthrough data, with data generated from collecting from the column outlet and quantifying nitrogen content via gas chromatography. Select experiments use a significantly larger (2.5x column diameter and length) column to demonstrate linear scalability of the process. Thermogravimetric analysis data is collected to evaluate the effects of thermal calcination as a sorbent regeneration method. Effects of a more complex feed are also investigated using a known reference fuel with additional added nitrogen containing compounds. The results presented in this work successfully demonstrate up to 99.8 % removal of NCCs from a model fuel fraction at an original NCC concentration of approximately 1500 ppm to single-digit parts per million after treatment. We also examine calcination of sorbent materials to remove the adsorbed species to enable sorbent reuse and minimize waste generation and show that the calcined material can be reused up to 5 cycles with reduced adsorption capacity. Overall, this work indicates that adsorptive denitrogenation using silica gel is a viable solution to enable the integration of HTL-derived aviation fuels into existing fuel infrastructure.

Adsorption techniques

Facet-Dependent Adsorption of Pb(II) on Hematite (001), (116), and (104) Surfaces

Hematite’s common (001) and (012) facets are frequently used in model studies of lead (Pb) adsorption behavior, but there is a lack of research on the high-energy facets, e.g. (104), present in nature. Also, few studies have attempted to connect the molecular details of facet-specific Pb adsorption to macroscopic uptake behavior. Here, to address these knowledge gaps, we investigated Pb(II) adsorption behaviors on facet-engineered hematite nanoparticles dominated by (001), (104), and (116). Adsorption experiments revealed significant variations in Pb(II) uptake among the three samples, with (001) demonstrating the highest capacity and (116) showing the best adsorption efficiency when normalized to specific surface area. Adsorption kinetics followed the pseudo-second-order model, indicating the adsorption process is governed mostly by chemisorption. Adsorption isotherms were well fitted by the Langmuir model, indicating uptake proceeds until roughly monolayer adsorption. Detailed characterization revealed Pb(II) was adsorbed as single atoms with complex inner-sphere binding modes that varied across different facets, indicating adsorption is both structurally and energetically facet-dependent. Co-adsorption experiments further demonstrated Cu 2+ , Zn 2+ , and humic acid significantly promoted Pb(II) adsorption. This study advances the understanding of hematite surface reactivity in controlling macroscopic wet adsorption behaviors, providing valuable insights into the environmental fate of Pb(II).

Contaminant removal

Liquid Phase Modeling in Porous Media: Adsorption of Methanol and Ethanol in H-MFI in Condensed Water

Zeolites are used in the chemical and separation industries for their exceptional selectivity, adsorption capacity, regenerability, and stability in gas and liquid phase processing. Here, we developed an explicit solvation method for predicting solvent/condensed phase effects on adsorption free energies in microporous media such as zeolites based on the hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) technique. Our explicit solvation method for zeolite systems, called eSZS, aims to capture site-specific interactions during the adsorption process at the Brønsted acid sites of H-MFI zeolite while still considering the diverse configuration space of the solvent molecules. This strategy is ideal for chemical reactions or adsorbates that interact with the microporous medium in few distinct adsorbate/transition state configurations, i.e., the harmonic or similar approximations are acceptable for the adsorbate/transition state while such approximations break down for the solvent molecules that require extensive configuration space sampling. In this way, our approach effectively overcomes the limitations of implicit solvation models and classical force field methods for describing solvation effects on chemical reactions within porous materials such as zeolites. Specifically, in this study, we investigated various aspects of our hybrid QM/MM approach, including QM cluster size dependencies in a periodic electrostatically embedded cluster model (PEECM), rules for link atoms at the QM/MM boundary, and functional and basis set considerations for converged and reasonably accurate gas and aqueous phase methanol and ethanol adsorption free energy predictions in H-MFI. For gas phase adsorption of methanol and ethanol in H-MFI at a Brønsted acid site in T12 position, we compute adsorption free energies at 298 K of −0.61 and −0.75 eV, respectively, using a PEECM containing 50 Si and 1 Al atom with ωB97x-D/def2-TZVP level of theory. For solvent effect calculations, we sample the aqueous phase using grand canonical Monte Carlo (GCMC) simulations to (1) obtain a mean field of electrostatic interactions in the reaction system and (2) perform a rigorous free energy perturbation calculation. Similar to the experimentally and computationally observed endergonic solvation effects observed for hydrocarbon adsorption on metal surfaces, we also observe that a condensed aqueous environment destabilizes methanol and ethanol at these acid sites in H-MFI at 298 K. Specifically, the computed solvation free energies of adsorption (ΔΔG solv ) for methanol and ethanol are +0.44 and +0.54 eV, respectively. From this study, it is evident that adsorbates (methanol and ethanol) are competing with water for adsorption space inside the H-MFI zeolite, leading to an endergonic solvation effect. Here, we expect that the endergonic, aqueous solvent effect during adsorption in microporous zeolites is highly tunable by changing the pore size and hydrophobicity of the microporous material as this will affect the water density inside the pore structure.

Adsorption

First-principles DFT modeling of nitrobenzene adsorption on the Ag(111) surface at varying monolayer coverages

Here, the adsorption behavior of nitrobenzene on the Ag(111) surface as a function of coverage is investigated using density functional theory. Adsorption energies and optimized geometries are analyzed together with isolated intermolecular interaction calculations and electronic structure analysis, including Bader charge partitioning and projected density of states, to disentangle the roles of adsorbate–surface bonding and through-space adsorbate–adsorbate interactions. At low and intermediate coverages (θ = 1/9 and θ = 2/9), bidentate adsorption at top sites is favored due to strong adsorbate–surface interactions, with additional stabilization at θ = 2/9 arising from favorable intermolecular separations. At higher coverage (θ = 1/3), bidentate adsorption is destabilized by strong intermolecular repulsion, and monodentate adsorption becomes energetically preferred as rotational freedom allows more favorable intermolecular spacing. Charge density differences, Bader charge, and density of states analyses show that charge transfer from the Ag surface is localized primarily on the nitro group and increases with coverage and adsorption denticity, although this increase does not directly correlate with adsorption strength at high coverage due to competing intermolecular interactions. These results demonstrate that surface coverage can induce a transition in preferred adsorption denticity driven by intermolecular interactions, highlighting the importance of adsorbate packing in organic molecule adsorption at metal interfaces.

36 MATERIALS SCIENCE

PFSA-Ionomer Adsorption to C and Pt/C Particles in Fuel-Cell Inks

Catalyst inks used to make fuel-cell electrodes consist of Pt/C catalyst particles and a perfluorosulfonic acid (PFSA) ionomer dispersed in water/alcohol solvent mixtures. PFSA ionomer in the ink adsorbs to the surface of the catalyst particles, dictating the dispersion colloid properties. Following adsorption, the subsequent distribution of excess nonadsorbed ionomer in the ink then governs the final structure of the electrode. Here, we characterize the adsorption of the PFSA ionomer onto Pt/C catalyst particles. PFSA adsorption is largely irreversible. Adsorbed sulfonic-acid moieties impart a negative charge on the catalyst surface, causing electrostatic repulsion between the free ionomer in solution and the ionomer-covered Pt/C particle surface. The amount of adsorption is limited by the resulting electrostatic charge that grows as more ionomer adsorbs, and the catalyst surface becomes more negatively charged. Attenuating electrostatic repulsion by increasing the ink ionic strength promotes ionomer adsorption. Electrostatically limited adsorption is observed, irrespective of the solvent water/n-propanol ratio or the catalyst particle porosity and Pt loading. Experimentally measured ionomer adsorption isotherms are well predicted by a Smoluchowski-based kinetic adsorption model, in which the electrostatic energy barrier for adsorption is predicted from DLVO theory. These findings help to unravel the complex phenomena within these colloidal dispersions, allowing for subsequent tailoring of inks to optimize fuel-cell electrode structure and performance.

Adsorption

Adsorption of REEs to Kaolinite via Ion Exchange and Surface Complexation as a Function of Water Chemistry

Rare earth elements (REEs) are critical components of modern technology behind renewable energy, transportation, and electronics but have a limited current supply. A substantial portion of global REE production relies on ion adsorption deposits. A high abundance of kaolinite in REE enrichment zones within these deposits suggests that kaolinite controls the subsurface migration of REEs. This study aimed to improve the current understanding of REE binding to kaolinite under varying water chemistry conditions. We conducted batch experiments with kaolinite (KGa-2) and three REEs (Nd, Dy, and Yb) at varying pH, electrolyte concentration, dissolved inorganic carbon (DIC), low molecular weight organic acids (citric and oxalic acids), and total REE concentration conditions. Increasing electrolyte concentration inhibits REE adsorption at pH < 7, suggesting that ion exchange contributes to adsorption at these pH values. DIC affects adsorption above pH 7–8 by forming strong aqueous complexes with heavy REEs. Citric acid decreases REE adsorption via aqueous complexation of REEs at pH > 5 but does not affect adsorption at pH < 5. The surface complexation model captures the main adsorption trends with two mechanisms: ion exchange on basal planes at pH < ∼6 and inner-sphere surface complexation to edge sites at pH > ∼6. Equilibrium constants for surface complexation increase in the order of Yb > Dy > Nd, indicating a higher strength of adsorption for heavy REEs. This study demonstrates how water chemistry conditions control the adsorption mechanisms that may determine the mobility of REEs in subsurface environments rich in kaolinite.

58 GEOSCIENCES

Glycerol Adsorption on TiO 2 Surfaces: A Systematic Periodic DFT Study

Abstract Conversion of glycerol to added‐value products is desirable due to its surplus during biodiesel synthesis. TiO 2 has been the most explored catalyst. We performed a systematic study of glycerol adsorption on anatase (101), anatase (001), and rutile (110) TiO 2 at the Density Functional Theory level. We found several adsorption modes on these surfaces, with anatase (101) being the less reactive one, leading to adsorption energies between −0.8 and −0.4 eV, with all adsorptions molecular in nature. On the contrary, anatase (001) is the most reactive surface, leading to both molecular and dissociative adsorption modes, with energies ranging from −4 to −1 eV and undergoing severe surface reconstructions in some cases. Rutile (110) also shows both molecular and dissociative adsorptions, but it is less reactive than anatase (001). Surfaces with oxygen vacancies affects the adsorbed states and energies. The electronic structure analysis reveals that glycerol adsorption mainly affects the band gap of the material and not the individual contributions to the valence and conduction band. Bader charge analysis shows that strong adsorption modes on anatase (001) and rutile (110) are associated with large charge transfer from glycerol to the surface, while weak and molecular adsorption modes involve low charge transfer.

Chemistry

Control of solvent adsorption in porous liquids through solvent-solvent interactions

Type 3 Porous Liquids (PLs) are a new class of materials that offer transformative potential for gas capture and utilization. These PLs are comprised of a bulky solvent with a suspended empty nanoporous material that enables selective and high-capacity capture from dilute and complex gas steams. The relationship between the nanoporous material, the structure of the solvent molecules, and the time dependence of solvent adsorption into the nanoporous material governs long term adsorption properties. Herein, molecular dynamics (MD) calculations evaluated the solvent dynamics of nine neat solvents with a broad range of chemical identities and experimental density measurements probed the time-dependent adsorption of these solvents into the ZIF-8 metal-organic frameworks (MOFs) to form PL dispersions over three weeks. Two of the nine dispersion, using glyceryl triacetate and 2’-hydroxyacetopheneone as solvents, presented sufficiently low solvent infiltration—characterized by solvent adsorption less than 40% of the ZIF-8 pore volume—to be viable as PLs. The experimental solvent adsorption data for ZIF-8 dispersed in the four aromatic solvents (acetophenone, methylbenzoate, 2-isopropylphenol, and 2’-hydroxyacetophenone) was fit to a kinetic model to quantify rate of solvent adsorption into ZIF-8. The rates of adsorption of these four similar-sized solvents demonstrated that solvent adsorption is slowed and limited by solvent clustering dynamics that are, in turn, driven by intermolecular hydrogen bonding. Ultimately, experimental solvent sorption measurements indicated that glyceryl triacetate and 2’-hydroxyacetophenone with ZIF-8 formed PLs with the most stable porosity and the greatest potential for gas capture. In conclusion, the combined experimental-computational materials exploration framework used here reveals a novel relationship between molecular scale solvent dynamics and macroscale solvent adsorption kinetics critical for the discovery and rapid evaluation of Type 3 PL compositions.

Hydrogen bonding

Adsorption of Rare Earth Elements to Kaolinite [dataset]

Rare earth elements (REEs) are critical components of modern technology behind renewable energy, transportation, and electronics but have a limited current supply. A substantial portion of global REE production relies on ion adsorption deposits. A high abundance of kaolinite in REE enrichment zones within these deposits suggests that kaolinite controls the subsurface migration of REEs. This study aimed to improve the current understanding of REE binding to kaolinite under varying water chemistry conditions. We conducted batch experiments with kaolinite (KGa-2) and three REEs (Nd, Dy, and Yb) at varying pH, electrolyte concentration, dissolved inorganic carbon (DIC), low molecular weight organic acids (citric and oxalic acids), and total REE concentration conditions. Increasing electrolyte concentration inhibits REE adsorption at pH 7, suggesting that ion exchange contributes to adsorption at these pH values. DIC affects adsorption above pH 7−8 by forming strong aqueous complexes with heavy REEs. Citric acid decreases REE adsorption via aqueous complexation of REEs at pH 5 but does not affect adsorption at pH 5. The surface complexation model captures the main adsorption trends with two mechanisms: ion exchange on basal planes at pH ∼6 and inner-sphere surface complexation to edge sites at pH ∼6. Equilibrium constants for surface complexation increase in the order of Yb Dy Nd, indicating a higher strength of adsorption for heavy REEs. This study demonstrates how water chemistry conditions control the adsorption mechanisms that may determine the mobility of REEs in subsurface environments rich in kaolinite.

58 GEOSCIENCES

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

Adsorption of hydroxamic acid ligands for improved extraction of rare earth elements from monazite ores

Efficient separation of rare earth element (REE) ores via froth flotation requires the development of novel ligands with enhanced adsorption capacity and selectivity. To realize these advances, understanding the mechanisms underlying interactions between the ligand and mineral surfaces is essential. This study systematically evaluates the adsorption behavior of alkyl- and aromatic alkyl-substituted hydroxamic acid ligands on monazite surfaces using complementary spectroscopic techniques, including UV–visible (UV–vis) spectroscopy, Raman spectroscopy, infrared spectroscopy, and vibrational sum frequency generation (SFG) spectroscopy, together with the ab initio molecular dynamics (AIMD) simulations. Among the studied ligands, octanohydroxamic acid (OHA) and 4-ethoxy-N,2-dihydroxybenzamide (EDHBA) exhibit high adsorption capacity under basic pH (8–10) by forming multilayers on the surface. OHA has a higher equilibrium adsorption capacity compared to EDHBA, but it forms a less stable multilayer susceptible to disruption in the presence of interfering ions. AIMD results show that OHA adopts a single stable chelating geometry, while EDHBA exhibits multiple binding modes involving distinct interactions with La surface atoms and phosphate-bound oxygens, resulting in more complex adsorption kinetics. The variations in surface binding and intermolecular interactions observed between alkyl and aromatic molecules influence the differences in adsorption kinetics, equilibrium adsorption capacities on the mineral surface, and their flotation performance. This work provides valuable insight into the adsorption mechanism of ligands at mineral interfaces, which is crucial for guiding the design of new ligands with enhanced separation performance.

Zhou, Muchu [ORNL] (ORCID:0000000182650215)

The Temperature Dependence of the Langmuir Adsorption Model for a Single-Site Metal–Organic Framework

The single-site Langmuir adsorption model, also known as the Langmuir isotherm equation, is one of the simplest possible descriptions of adsorption phenomena, and yet finds widespread applicability across a range of disciplines. In its simplest form, it is deployed to treat adsorption equilibria at constant temperature (i.e., along isotherms); however, at the heart of its derivation is a more general class of models that each incorporate an explicit temperature dependence, subject to assumptions about the spatial/translational degrees of freedom of the adsorbed species. In this work, measurements of the temperature dependence of supercritical adsorption of H 2 on a single-site metal-organic framework (MOF) are presented, and fitted using a range of Langmuir models with distinct treatments of degrees of freedom in the adsorbed phase. Surprisingly, all of the models can be used to adequately represent the measured data, despite yielding significantly different values for binding energy and the temperature dependence of the isosteric enthalpy of adsorption (i.e., the isosteric heat, qst). Furthermore, a critical finding is that the mean-temperature isosteric enthalpy of adsorption remains consistent across all models within experimental error, highlighting its reliability for evaluating adsorption thermodynamics.

08 HYDROGEN

Complex Adsorption Behavior of Neodymium and Ytterbium on Structurally-Distinct Alumina Surfaces

New sources of rare earth elements (REEs) are needed to support a green energy transition. REEs adsorbed to aluminum-rich clays in weathering deposits represent important resources but the mechanisms responsible for their retention and ease of extraction are unresolved. Disordered coordination and co-occurrence of multiple species pose challenges to investigating REE adsorption processes via established spectroscopic methods. In this study, we applied element-specific surface crystallography methods to obtain a new perspective on the complexity of REE adsorption mechanisms and affinities. Alumina (001) and (012) crystal surfaces were utilized to evaluate surface-specific controls on Nd(III) and Yb(III) adsorption behavior. The REEs displayed similar total adsorption to alumina (001) as a mixture of inner- and outer-sphere complexes, but Nd displayed a greater proportion of inner-sphere binding. Adsorption of ordered inner- and outer-sphere REE species was substantially lower on alumina (012). These distinct behaviors reflect differences in the surface functional group charging and topography of the two surfaces. However, alumina (012) also hosted a substantial population of disordered adsorbed species, especially for Nd, potentially associated with Al vacancy surface defects. Here, the accumulation of light versus heavy REEs via adsorption in weathering deposits likely results from multiple, competing reactions affected by clay particle morphology. Leaching procedures for resource recovery should account for differential rates of desorption by coexisting inner- and outer-sphere REE surface complexes.

58 GEOSCIENCES

Adsorption Studies of Dilute Krypton and Xenon from Nitrogen on SBMOF-1 and Activated Charcoal for Applications in Isotope Harvesting

Adsorptive partitioning of dilute krypton (Kr) and xenon (Xe) onto Stony Brook Metal–Organic Framework (SBMOF-1) and activated charcoal (AC) from carrier nitrogen was experimentally measured at temperatures ranging from 195 to 293 K. For this purpose, a closed-loop system for gas adsorption experiments was developed. From the Kr adsorption measurements, the adsorption equilibrium constant for Kr on SBMOF-1 was calculated, yielding a value for the enthalpy of adsorption of −19 ± 1 kJ·mol –1 . The partition coefficients were utilized to estimate the extraction rates of 76 Kr, 77 Kr, and 122 Xe isotopes during isotope harvesting at the Facility for Rare Isotope Beams (FRIB). We conclude that both materials showed promising results for the extraction of noble gases from FRIB effluents using temperature swing adsorption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

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

Metal–organic frameworks (MOFs) offer promise as selective CO 2 sorbents, but successful MOF sorbent materials need high CO 2 binding affinity and selectivity for CO 2 over water. This work focuses on the use of machine-learned force fields (MLFFs) to model CO 2 adsorption in flexible MOFs, with a focus on SIFSIX-3-Cu, an anion-pillared MOF known for its high CO 2 affinity. A preliminary high-throughput screening of over 900 anion-pillared MOFs was performed using rigid UFF+DDEC6 force fields to predict zero-loading heats of adsorption for CO 2 and H 2 O. SIFSIX-3-Cu was selected for further computational study due to its predicted CO 2 heat of adsorption and experimental relevance. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies and forces, with an iterative sampling scheme combining molecular dynamics, geometry optimization, random geometric insertion, and NVT Monte Carlo-based configuration generation to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included, contrasting with previous models that approximated the MOF as rigid. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations with the MLFF produced CO 2 adsorption isotherms in good agreement with experimental data at direct air capture (DAC) pressures (e.g., 40 Pa), in contrast to previous overestimations of CO 2 sorption by models with rigid structures. Bond and angle histogram analysis showed that MOF flexibility increased the variance of fluorine–fluorine diagonal distances at adsorption sites, resulting in a lower predicted sorption for flexible, asymmetric SIFSIX-3-Cu pore geometries compared to the rigid, symmetric DFT-optimized SIFSIX-3-Cu pore geometry. A detailed description of flexibility afforded by the MLFF resulted in an accurately predicted CO 2 uptake (0.88 mmol/g) at low pressure (40 Pa) compared to the experimentally measured value (1.24 mmol/g). In conclusion, these results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low-pressure applications.

adsorption