Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Water Uptake by Silica Nanopores: Impacts of Surface Hydrophilicity and Pore Size

Water uptake at 298 K by two mesoporous silicas with different pore sizes was studied using volumetric vapor sorption. Here, through variation of sample pretreatment temperature and time, the number of surface hydroxyl groups was varied, leading to marked changes in the water sorption behavior. With increasing surface hydroxyl density, the adsorption of water to the pore surfaces increases, leading to formation of a thicker and denser water sorption layer. The BET model was used to measure surface hydroxyl density from low pressure parts of water adsorption isotherms. Adsorption of water molecules scales with the number density of surface hydroxyl groups, and monolayer coverage is reached at reduced pressure of ca. 0.3. The onset of adsorption pore condensation of water shifts to lower reduced pressure with increasing surface hydroxyl density, indicating growing thickness of adsorption layers. However, the water pore condensation step of the desorption branch shifts to smaller reduced pressure with increasing surface hydroxylation, reducing the adsorption-desorption hysteresis width. An adsorbed phase model is utilized to calculate the distribution of pore water molecules between adsorbed and pore condensation phases, and characterize the density and thickness of the water sorption phase. In the smaller pores, the adsorbed phase forms a sparse monolayer, while approximately a bilayer is formed in the wider pores. The previously reported silica nanopore underfilling by water is confirmed and rationalized by a significantly reduced water density in the adsorbed phase with respect to bulk, while the pore core is filled with water of approximately bulk density.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting the Rate of Homogeneous Intermetallic Nucleation within Steep Composition Gradients

Simulations of isothermal homogeneous nucleation from deeply undercooled amorphous melts exhibit systematic variations in nucleation behavior depending upon the strength of an imposed composition gradient. Data from molecular dynamics (MD) simulations in a model Ni/Al system permit quantification of the nucleation rate of the NiAl-B2 intermetallic phase and indicate that nucleation proceeds in a polymorphous mode. The nucleation rate decreases with increasing gradient, and nucleation is completely suppressed above a critical gradient. Based on an argument of crystal nucleus stability, a simple estimate provides good prediction of the critical gradient. Here, a modified classical nucleation model parameterized with thermodynamic and kinetic quantities calculated independently predicts the nucleation rates and their variation with the imposed gradient, matching the MD results very well.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Adsorbate Partition Functions via Phase Space Integration: Quantifying the Effect of Translational Anharmonicity on Thermodynamic Properties

A new method for computing anharmonic thermophysical properties for adsorbates on metal surfaces is presented. Classical Monte Carlo phase space integration is performed to calculate the partition function for the motion of a hydrogen atom on Cu(111). Here, a minima-preserving neural network potential energy surface is used within the integration routine. Two different sampling schema for generating the training data are presented, and two different density functionals are used. The results are benchmarked against direct state counting results by using discrete variable representation. The phase space integration results are in excellent quantitative agreement with the benchmark results. Additionally, both the discrete variable representation and the phase space integration results confirm that the motion of H on Cu(111) is highly anharmonic. The results were applied to calculate the free energy of dissociative adsorption of H 2 and the resulting Langmuir isotherms at 400, 800, and 1200 K in a partial pressure range of 0–1 bar. It shows that the anharmonic effects lead to significantly higher predicted surface site fractions of hydrogen.

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

Adsorption and Catalytic Activity of Gold Nanoparticles in Mesoporous Silica: Effect of Pore Size and Dispersion Salinity

The assembled state of nanoparticles (NPs) within porous matrices plays a governing role in directing their biological, electronic, and catalytic properties. However, the effects of the spatial confinement and environmental factors, such as salinity, on the NP assemblies within the pores are poorly understood. In this study, we use adsorption isotherms, spectrophotometry, and small-angle neutron scattering to develop a better understanding of the effect of spatial confinement on the assembled state and catalytic performance of gold (Au) NPs in propylamine-functionalized SBA-15 and MCM-41 mesoporous silica materials (mSiO 2 ). We carry out a detailed investigation of the effect of pore diameter and ionic strength on the packing and spatial distribution of AuNPs within mSiO 2 to get a comprehensive insight into the structure, functioning, and activity of these NPs. We demonstrate the ability of the adsorbed AuNPs to withstand aggregation under high salinity conditions. We attribute the observed preservation of the adsorbed state of AuNPs to the strong electrostatic attraction between oppositely charged pore walls and AuNPs. The preservation of the structure allows the AuNPs to retain their catalytic activity for a model reaction in high salinity aqueous solution, here, the reduction of p-nitrophenol to p-aminophenol, which otherwise is significantly diminished due to bulk aggregation of the AuNPs. This fundamental study demonstrates the critical role of confinement and dispersion salinity on the adsorption and catalytic performance of NPs.

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

Molecular Structure of Adsorbed Water Phases in Silica Nanopores

The adsorption of water vapor in silica nanopores with different pore morphologies and surface hydrophilicities was studied to quantify the densities and thicknesses of the water sorption layers and deduce their molecular structures. Water adsorption to surface hydroxyls is described by a multilayer sorption model. At low pressure, the water adsorption isotherms are largely independent of pore size and the adsorbed amounts scale with the surface hydroxyl density. Adsorbed-phase densities corresponding to the adsorption of two water molecules per surface hydroxyl group are found in the first adsorbed water layer for a wide range of surface hydroxyl densities. The densities and layer thickness values found in narrow pores indicate that patchy adsorbed layers form if not enough water molecules exist for a full layer, which coexist with dry pore surface regions. This behavior indicates cooperative adsorption effects, i.e., a preference for the formation of hydrogen bonds between water molecules bound to surface hydroxyls. In narrow pores, pore condensation limits further growth of the sorption layer, and in larger pores and at the planar quartz surface, a second adsorbed water layer is formed, which can hold up to approximately 4 additional water molecules per surface hydroxyl group. The water molecules in these thicker adsorption layers arrange such that the sorption layer density is similar to the bulk water density. Pore confinement limitations on the sorption layers are observed in pores with radii of as large as 8 nm. Molecular dynamics modeling reveals two preferential orientations for water molecules adsorbing to surface hydroxyl groups and suggests an intralayer structuring in the adsorbed monolayer. Furthermore, adsorbed water molecules in the sorption layer are bonded to the surface hydroxyl group via the donation and acceptance of hydrogen bonds.

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

Hydrogen Uptake Kinetics of 1,4-Bis(phenylethynyl)benzene Rubberized O-Rings: Measurements, Modeling, and Comparison with Additional Forms of Organic Getters

A class of molecules called “hydrogen getters” can react with, or scavenge, H 2 in applications where the hydrogen presence and/or buildup are not desirable. One such “getter”, 1,4-bis(phenylethynyl)benzene or DEB, can be incorporated into a silicone matrix in the form of an O-ring for added flexibility, environmental resiliency, and convenient use as a gasket in sealed applications. However, the performance and kinetics of this DEB-loaded rubberized O-ring have not yet been characterized. Here in this work, the hydrogen uptake kinetics of the rubberized DEB O-ring were extracted by the isoconversional analysis from isothermal isobaric data under conditions of 13,332 Pa H 2 and 305–325 K. The isoconversional and cylindrical diffusion approximations were then used to predict and model the hydrogen uptake of rubberized DEB O-rings under any arbitrary condition, as illustrated in this work for the simple case of a constant rate of hydrogen generation/input. In comparison with other Pd/carbon-based/organic getter systems, these rubberized DEB O-rings provide a type of hydrogen getter material with enhanced flexibility and catalyst protection in applications where an O-ring seal is needed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modifying Interface Solvation and Oxygen Reduction Electrocatalysis with Hydrophobic Species

Ionic liquid interlayers improve the oxygen reduction reaction (ORR) kinetics on bulk and nanostructured catalysts for both Pt and alloyed-Pt materials. Despite the demonstrated performance enhancement at the half-cell and membrane electrode assembly level, the mechanism of the improvement is not fully understood. In this work, we combine single-crystal experiments with microkinetic modeling to uncover the origin of the ORR kinetic improvement on Pt(111) in the presence of ionic liquids and hydrophobic cations. With the incorporation of a modified Frumkin isotherm, our model accurately simulates the disorder–order transition observed in hydroxyl and bisulfate adsorption on Pt(111) under acidic conditions. Voltametric analysis shows that ionic liquids impact solvation to break so-called scaling relations between the adsorption strength of OH ad and O ad , but these effects have little impact on ORR activity. Instead, destabilized OH ad reduces the overall hydroxyl (spectator) coverage, resulting in higher availability of active sites.

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 OSTI2022

Theoretical Pathway toward Improved Reverse Osmosis Membrane Selectivity for Neutral Solutes: Inspiration from Gas Separations

We report highly rigid membrane materials with tailored structures have exhibited permeabilities and selectivities that exceed the polymer upper bound in gaseous and organic solvent separations for challenging mixtures containing species that are similar in size and shape. One potential question is whether such membrane materials can maintain meaningful guest diffusivities in situations where the microporous spaces are essentially full of guest molecules. Here, we use a simplified transition state theory approach to estimate the diffusivity of water and small organics within a microporous membrane. The transition state theory model is parameterized using experimental values from zeolites and carbon molecular sieve materials found in the literature. We demonstrate the differences in transport and Maxwell–Stefan diffusivities based on guest species loading with different isotherm behaviors. These calculations theorize a path forward for highly selective reverse osmosis membranes for aqueous phase separations.

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

Multiscale Reactive Model for 1,3,5-Triamino-2,4,6-trinitrobenzene Inferred by Reactive MD Simulations and Unsupervised Learning

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

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

Synthesis and Properties of the Helium Clathrate and Defect Perovskite [He 2– x $^\square$ x ][CaNb]F 6

The defect double perovskite [He 2– x $^\square$ x ][CaNb]F 6 , with helium on its A-site, can be prepared by the insertion of helium into ReO 3 -type CaNbF 6 at high pressure. Upon cooling from 300 to 100 K under 0.4 GPa helium, ~60% of the A-sites become occupied. Helium uptake was quantified by both neutron powder diffraction and gas insertion and release measurements. After the conversion of gauge pressure to fugacity, the uptake of helium by CaNbF 6 can be described by a Langmuir isotherm. The enthalpy of absorption for helium in [He 2– x $^\square$ x ][CaNb]F 6 is estimated to be ~+3(1) kJ mol –1 , implying that its formation is entropically favored. Helium is able to diffuse through the material on a time scale of minutes at temperatures down to ~150 K but is trapped at 100 K and below. The insertion of helium into CaNbF 6 reduces the magnitude of its negative thermal expansion, increases the bulk modulus, and modifies its phase behavior. On compressing pristine CaNbF 6 , at 50 and 100 K, a cubic (Fm$\overline{3}$m) to rhombohedral (R$\overline{3}$) phase transition was observed at <0.20 GPa. However, a helium-containing sample remained cubic at 0.4 GPa and 50 K. CaNbF 6 , compressed in helium at room temperature, remained cubic to >3.7 GPa, the limit of our X-ray diffraction measurements, in contrast to prior reports that upon compression in a nonpenetrating medium, a phase transition is detected at ~0.4 GPa.

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

Revealing the Influence of Binding Motifs on Electron Transfer and Recombination Kinetics for CdSe Quantum Dots Functionalized with a Modified Viologen

Anchoring of molecules to the surfaces of semiconductor nanocrystals (NCs) presents an opportunity to leverage the precise synthetic tunability of molecular function and the remarkable light harvesting properties of NCs to drive photochemical reactions. However, charge transfer between the two species depends not only on the energy level alignments but also on the details of their binding interactions, which are difficult to probe. Here, in this work, we characterize the binding between CdSe quantum dots (QDs) and a new phosphonated derivative of the electron acceptor methyl viologen, designed to attach to the QD surface via the phosphonate group. We use isothermal titration calorimetry to probe the thermodynamics of the QD–molecule interaction and use the parameters determined therein to analyze transient absorption spectroscopy measurements of forward and back electron transfer from QDs to the viologen. We find that the ligand-like phosphonate binding leads to an electron-transfer rate constant that is 3 orders of magnitude smaller than that for the face-on binding of the bipyridine ring of methyl viologen. Back electron transfer is also significantly slower in the derivative. Interestingly, a minor fraction of the phosphonated derivative also binds in the face-on configuration, with similar forward and back electron transfer kinetics as methyl viologen. Numerical simulations show that the ligand-like binding will lead to significantly improved quantum yields of photocatalysis over a wide range of reaction rates. By independently characterizing binding thermodynamics and charge transfer kinetics, this work reveals how the complexities underlying electron transfer at the NC–molecule interface determine photocatalytic outcomes. This work also represents a step toward controlling forward and back electron transfer kinetics via rational molecular design.

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

Competitive Sorption of H 2 O and CO 2 in Clay Mineral Interlayers

The first-reported competitive CO 2 –H 2 O sorption isotherms on a naturally occurring expandable clay mineral (saponite) exchanged with Na, K, and Cs at near ambient CO 2 concentrations (1% CO 2 ) are reported by combining ATR-FTIR spectroscopy with gravimetric sorption methods. Little CO 2 sorption was observed at high relative humidity (RH) but increased significantly at low RH. Unlike prior work where CO 2 sorption was reported on samples that were previously subjected to heating and/or evacuation treatment to remove sorbed H 2 O, this study manipulated CO 2 sorption by controlling RH. CO 2 sorption increased in the order Cs- > K- > Na-exchanged saponite and was anticorrelated with the Gibbs energy of hydration of the exchangeable cation. As expected, H 2 O sorption followed the opposite trend, with increased H 2 O sorption occurring in the order Na > K > Cs. The amount of CO 2 sorbed using 1 atm CO 2 (p/p o = 1.0) ranged from 50 to 30 mg CO 2 /g clay and from 1.25 to 0.8 mg CO2 /g clay for 1% CO 2 (p/p o = 0.01) in N 2 . The sensitivity of the ATR-FTIR measurements is demonstrated by detecting the CO 2 stretch of adsorbed CO 2 from air containing 500 ppmv CO 2 . The spectral features of sorbed CO 2 were minimally affected by the nature of the exchangeable cation or variations in H 2 O content. However, the position of the HOH bending band shifted significantly from 1638 to 1610 cm –1 as H 2 O content decreased, reflecting changes in intermolecular hydrogen bonding between H 2 O molecules in the saponite interlayer. The shift in position of the HOH bending band to lower energy was coincident with increased CO 2 sorption. CO 2 sorption in clay mineral interlayers is consistent with sorption on weakly hydrated partially hydrophobic sites found on low charge density smectites and is enhanced by lower RH and by the presence of weakly hydrated exchangeable cations like Cs + or K + . In conclusion, the greatest CO 2 sorption occurs on Cs-saponite because this is the least hydrated cation of the three used in the study and indicates that CO2 sorption occurs on the neutral portion of the siloxane surface.

Johnston, Cliff T. [Purdue Univ., West Lafayette, ↗
DOE OSTI2025

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↗
DOE OSTI2026

Identification of Carbonyl Species on Palladium Supported on Ceria in Complex Microenvironments

Herein, we present a systematic comparison between Pd carbonyl (Pd-CO) species, specifically over Pd/CeO 2 based catalysts, observed during isothermal adsorption and in several prototypical catalytic reactions to identify and understand CO adsorption on palladium-ceria based catalysts. Pd-CO is observed via DRIFTS to probe the gas-solid conditions, while ATR-IR is used to probe the affinity of Pd-CO under more complex solvated gas-solid-liquid conditions to discern the influence of the microenvironments for carbonyl adsorption. Here, we explore the presence of Pd-CO under several reactive environments, including CO adsorption, CO 2 + H 2 , CO + H 2 , CH 4 + CO 2 and CO under gas-solid-liquid media, highlighting reactions with notable Pd-CO formation. The differences between palladium carbonyls and carbonate species show that carbonyl species are much more affected via a shifting of the peak position than carbonates, which remain static irrespective of the immediate chemical environment. By following the rate of CO accumulation via K-M mode DRIFTS, we observe migration from linear, 2095 cm -1 , to bridge site, 1978 cm -1 , as a function of time under a static CO atmosphere. With the use of DFT, we discerned changes in Pd-carbonyl stretches due to both coverage effects of CO under simulated reaction conditions and temperature effects. Regardless of whether CO is formed as an intermediate or a reactant, the competitive adsorption of *H and *CO affects the binding strength of *CO at all temperatures, with low temperature favoring atop binding and high temperature favoring the more stable FCC Pd-CO site.

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

Transient Pulse-Response Time-of-Flight Mass Spectrometry for Complex, Deactivating Heterogeneous Catalytic Systems: Application to Ethane Dehydroaromatization

The study of complex, multistep bond-forming and -breaking reactions in heterogeneous catalytic systems often encounters challenges associated with the involvement of large numbers of intermediates among branching pathways. Kinetic information obtained from traditional steady-state measurements can be complemented with that from time-resolved methods to uncover details of the underlying chemistry. Herein, we describe an approach for tracking the complete time-resolved chemical composition (ca. 4–200 u) of a reactor effluent in response to a reactant pulse. We use a six-port rotary valve with a metered sampling loop to pulse reactants at ambient pressure into a flow reactor packed with a catalyst bed within the isothermal region of a heated furnace. The temporal evolution of effluent species is tracked using time-resolved molecular-beam time-of-flight mass spectrometry. We highlight the possibilities that this method has to offer by studying the complex bifunctional mechanism of ethane dehydroaromatization over an HZSM-5-supported platinum catalyst. We demonstrate that energy-tunable ionization sources, which facilitate isomer resolution, enable the measurement of the full mass spectral time-dependent system response. This includes the evolution of major products and mechanistically relevant reactive intermediates such as 1,3-butadiene and cyclopentadiene; these species have not previously been observed from this reaction. In additional studies, we also assess the role of platinum in the catalyst by examining temporal responses to ethane and ethylene feeds. Results show that two temporally distinct formation pathways exist for methane and benzene, and that their importance depends on both catalyst composition and reactant identity. Additionally, characteristics of catalyst deactivation are uniquely observable in the time-resolved mass spectral response, including the selective deactivation of a benzene formation pathway. The combination of time-of-flight mass spectrometry with tunable ionization enables the simultaneous observation of all effluents in a complex mixture of intermediates with isomer/isobar differentiation capabilities that can be applied to any complex reaction system.

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

Machine Learning Interatomic Potentials for Modeling Framework Flexibility and Water Uptake in NbOFFIVE-1-Ni Metal–Organic Framework

Metal–organic frameworks (MOFs), with their distinctive porous structures and tunable chemical properties, have shown immense promise in the separation and storage of gases. Currently, the accurate simulation of their adsorptive properties remains challenging, especially for systems where the molecules fit very tightly into the pores. Traditional simulation methods often approximate the frameworks as rigid and do not account for the framework flexibility seen in materials such as NbOFFIVE-1-Ni. First-principles molecular dynamics (FPMD) simulations offer the desired accuracy in modeling this flexibility but are limited by their extensive computational demands, rendering them impractical for long simulations. Conversely, classical force field-based simulations offer computational efficiency but lack the necessary accuracy. Here, to break this accuracy-efficiency trade-off, we have developed machine learning interatomic potentials trained on energies and forces from FPMD to model the framework flexibility of NbOFFIVE-1-Ni in the presence of water over nanosecond time scales. Furthermore, by integrating MLIP-driven molecular dynamics (MLIP-MD) with grand canonical Monte Carlo (GCMC) simulations, we further incorporated framework flexibility into adsorption predictions, yielding water adsorption isotherms that better align with experimental data compared to those of conventional GCMC simulations. These advances offer new opportunities for the design and optimization of MOFs in gas storage and separation applications.

adsorption↗
DOE OSTI2022

Anomalies and Local Structure of Liquid Water from Boiling to the Supercooled Regime as Predicted by the Many-Body MB-pol Model

For the past 50 years, researchers have sought molecular models that can accurately reproduce water’s microscopic structure and thermophysical properties across broad ranges of its complex phase diagram. Herein, molecular dynamics simulations with the many-body MB-pol model are performed to monitor the thermodynamic response functions and local structure of liquid water from the boiling point down to deeply supercooled temperatures at ambient pressure. The isothermal compressibility and isobaric heat capacity show maxima near 223 K, in excellent agreement with recent experiments, and the liquid density exhibits a minimum at ~208 K. A local tetrahedral arrangement, where each water molecule accepts and donates two hydrogen bonds, is found to be the most probable hydrogen-bonding topology at all temperatures. Further, this work suggests that MB-pol may provide predictive capability for studies of liquid water’s physical properties across broad ranges of thermodynamic states, including the so-called water’s “no man’s land” which is difficult to probe experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗