Search NASA⌕ Search

Engineering topics

Snurr, Randall Q.

Publications and source records attributed to Snurr, Randall Q..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

gRASPA

GPU Monte Carlo Simulation Code with a taste of RASPA We present enhancements in Monte Carlo simulation speed and functionality within an open-source code, gRASPA, which uses graphical processing units (GPUs) to achieve significant performance improvements compared to serial, CPU implementations of Monte Carlo. The code supports a wide range of Monte Carlo simulations, including canonical ensemble (NVT), grand canonical, NVT Gibbs, Widom test particle insertions, and continuous-fractional component Monte Carlo. Implementation of grand canonical transition matrix Monte Carlo (GC-TMMC) and a novel feature to allow different moves for the different components of metal-organic framework (MOF) structures exemplify the capabilities of gRASPA for precise free energy calculations and enhanced adsorption studies, respectively. The introduction of a High-Throughput Computing (HTC) mode permits many Monte Carlo simulations on a single GPU device for accelerated materials discovery. The code can incorporate machine learning (ML) potentials. The open-source nature of gRASPA promotes reproducibility and openness in science, and users may add features to the code and optimize it for their own purposes. The code is written in CUDA/C++ and SYCL/C++ to support different GPU vendors. The gRASPA code is publicly available at https://github.com/snurr-group/gRASPA.

Li, Zhao [Purdue/Northwestern/Notre Dame Universit↗
DOE OSTI2024

Photocatalytic Semi‐Hydrogenation of Acetylene to Polymer‐Grade Ethylene with Molecular and Metal–Organic Framework Cobaloximes

Abstract The semi‐hydrogenation of acetylene in ethylene‐rich gas streams is a high‐priority industrial chemical reaction for producing polymer‐grade ethylene. Traditional thermocatalytic routes for acetylene reduction to ethylene, despite progress, still require high temperatures and high H 2 consumption, possess relatively low selectivity, and use a noble metal catalyst. Light‐powered strategies are starting to emerge, given that they have the potential to use directly the abundant and sustainable solar irradiation, but are ineffective. Here an efficient, >99.9% selective, visible‐light powered, catalytic conversion of acetylene to ethylene is reported. The catalyst is a homogeneous molecular cobaloxime that operates in tandem with a photosensitizer at room temperature and bypasses the use of non‐environmentally friendly and flammable H 2 gas feed. The reaction proceeds through a cobalt‐hydride intermediate with ≈100% conversion of acetylene under competitive (ethylene co‐feed) conditions after only 50 min, and with no evolution of H 2 or over‐hydrogenation to ethane. The cobaloxime is further incorporated as a linker in a metal–organic framework; the result is a heterogeneous catalyst for the conversion of acetylene under competitive (ethylene co‐feed) conditions that can be recycled up to six times and remains catalytically active for 48 h, before significant loss of performance is observed.

Stone, Aaron E. B. S.↗
DOE OSTI2024

Predicting Partial Atomic Charges in Metal–Organic Frameworks: An Extension to Ionic MOFs

Molecular simulation is an invaluable tool to predict and understand the usage of metal–organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges, commonly obtained from density functional theory (DFT) calculations, are often required to model the electrostatic interactions between the MOF and adsorbates, especially when the adsorbates have dipole or quadrupole moments, such as water and CO 2 . Machine learning (ML) models have been previously employed to predict partial charges and avoid the computational cost associated with DFT calculations. However, previous ML models suffer from small training data sets, which limit their scope of application. In this work, we introduce two novel machine learning models, PACMOF2-neutral and PACMOF2-ionic, aimed at predicting the density-derived electrostatic and chemical (DDEC6) partial atomic charges for both neutral and ionic MOFs. These models not only yield DFT-level accuracy at a fraction of the computational cost but also demonstrate a remarkable improvement in prediction of adsorption, as validated with grand canonical Monte Carlo simulations. Furthermore, the robustness and fast computational time of the PACMOF2 models, along with their transferability to other porous materials such as covalent organic frameworks and zeolites, underscores their potential in high-throughput screening of MOFs for diverse applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Titanium-Based Metal–Organic Framework For Tandem Metallaphotocatalysis

Metal–organic frameworks (MOFs) have garnered substantial attention for their unique properties, such as high porosity and tunable structures, making them versatile for various applications. This paper constructs photoactive titanium–organic frameworks by combining Ti(IV) clusters and a bipyridine linker. The MOF is synthesized in situ through imine condensation, resulting in NU-2300. Subsequent ex situ nickel salt complexation results in NU-2300-Ni, which is then used for light-mediated carbon–heteroatom cross-couplings. The photophysical properties of the metallaphotocatalyst were investigated by UV-vis and EPR analysis, and both the Ti cluster and the bipyridine linker were found to contribute to successful catalysis, making it a tandem catalyst. The heterogeneous material retained its performance through five cycles of thioetherification. As a result, this work contributes not only to MOF synthetic strategies, but also to expanding MOF applications as recyclable, tandem metallaphotocatalysts.

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

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2024

An active, stable cubic molybdenum carbide catalyst for the high-temperature reverse water-gas shift reaction

Although technologically promising, the reduction of carbon dioxide (CO 2 ) to produce carbon monoxide (CO) remains economically challenging owing to the lack of an inexpensive, active, highly selective, and stable catalyst. We show that nanocrystalline cubic molybdenum carbide (α-Mo 2 C), prepared through a facile and scalable route, offers 100% selectivity for CO 2 reduction to CO while maintaining its initial equilibrium conversion at high space velocity after more than 500 hours of exposure to harsh reaction conditions at 600°C. The combination of operando and postreaction characterization of the catalyst revealed that its high activity, selectivity, and stability are attributable to crystallographic phase purity, weak CO-Mo 2 C interactions, and interstitial oxygen atoms, respectively. Mechanistic studies and density functional theory (DFT) calculations provided evidence that the reaction proceeds through an H 2 -aided redox mechanism.

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

Constraining Flexibility in the MIL-88 Topology through Integration of 3-Dimensional Linkers

Metal–organic frameworks (MOFs) are a class of crystalline, nanoporous materials that are recognized for their tunability. While some MOFs demonstrate flexibility, this characteristic can pose challenges in achieving precise pore control or establishing permanent porosity. Specifically, MIL-88B, is notable for its high flexibility, as it is constructed from metal trimer clusters and 2-dimensional linkers (2DLs) featuring planar, aromatic cores, allowing significant structural changes. In this study, we synthesized two new MOFs, NU-2010 and NU-2011, which are structurally analogous to MIL-88B but incorporate ditopic 3-dimensional linkers (3DLs) with sterically bulky cores and higher symmetry. Our aim was to investigate whether the introduction of 3DLs could mitigate the flexibility observed in MIL-88B. We employed a combination of single crystal and powder X-ray diffraction techniques to assess the flexibility of MIL-88B, NU-2010, and NU-2011 under various conditions including thermal activation, solvent exchange, and temperature changes. Furthermore, our findings indicate that incorporating 3DLs significantly reduces framework flexibility in NU-2010 and NU-2011 relative to MIL-88B.

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

Exploring the Structural, Dynamic, and Functional Properties of Metal–Organic Frameworks through Molecular Modeling

This review spotlights the role of atomic-level modeling in research on metal-organic frameworks (MOFs), especially the key methodologies of density functional theory (DFT), Monte Carlo (MC) simulations, and molecular dynamics (MD) simulations. The discussion focuses on how periodic and cluster-based DFT calculations can provide novel insights into MOF properties, with a focus on predicting structural transformations, understanding thermodynamic properties and catalysis, and providing information or properties that are fed into classical simulations such as force field parameters or partial charges. Classical simulation methods, highlighting force field selection, databases of MOFs for high-throughput screening, and the synergistic nature of MC and MD simulations, are described. By predicting equilibrium thermodynamic and dynamic properties, these methods offer a wide perspective on MOF behavior and mechanisms. Additionally, the incorporation of machine learning (ML) techniques into quantum and classical simulations is discussed. These methods can enhance accuracy, expedite simulation setup, reduce computational costs, as well as predict key parameters, optimize geometries, and estimate MOF stability. By charting the growth and promise of computational research in the MOF field, the aim is to provide insights and recommendations to facilitate the incorporation of computational modeling more broadly into MOF research.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular Simulations of the Chain Length Dependent Adsorption of C7‐C14 n ‐Alkanes in ZIF‐8

Recent experiments show that the diffusivities of C8 and C10 n-alkanes in ZIF-8 are higher than those of C7 and C9, respectively. We investigated this unusual ‘odd-even' effect by simulating the adsorption of C7-C14 n-alkanes in ZIF-8 using hybrid Monte Carlo molecular simulations. The resultant adsorption isotherms, guest-host energies, isosteric heats, and chain length distributions are analyzed for trends among the n-alkanes. ZIF-8 cages filled with n-alkanes are characterized using a combination of image processing and data science to quantify the differences in packing that occur with chain length. In conclusion, results indicate that packing changes drastically from C7 to C8 and from C9 to C10, which is consistent with the diffusion trends previously reported in experiments.

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

Single-Step Selective Oxidation of Methane by Iron-Oxo Species in the Metal–Organic Framework MFU-4l

The direct and selective conversion of methane to methanol can be considered a holy grail for catalysis research. In this work, we study a metal-organic framework known as MFU-4l, modified by design to include highly reactive iron-oxo species for the catalytic C-H bond activation of methane. We investigate the oxidation of methane and the further potential oxidation of the product methanol using N 2 O as an oxygen source and map the potential energy landscape of these reactions using density functional theory calculations. We show that the highest energy barrier encountered during the methane oxidation process is not the C-H bond breaking, but the activation of the iron center by N 2 O. Furthermore, the potential energy landscape for the C-H bond activation exhibits a large, high-energy plateau region instead of a sharp transition state, thus differing from the traditional radical rebound mechanism. Furthermore, this insight offers interesting potential routes to enhance the catalytic activity of the catalyst, to hinder unwanted deactivation pathways, and to reduce the activity towards the over-oxidation of the product.

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

Defect engineering of porous aromatic frameworks via end capping improves dioxane removal from water

Amorphous porous organic polymers show promise for energy-efficient adsorptive separations, but it is difficult to understand or improve their performance through intentional structural modification. Herein, we report the synthesis of porous aromatic frameworks (PAFs) with pore structures tailored by the incorporation of a monofunctionalized end-capping monomer that disrupts framework topology. Combining experimental characterization with molecular simulations, we show that this defect engineering strategy yields less densely crosslinked networks, which leads to pore collapse and the presence of unique adsorption sites. These defect-engineered PAFs exhibit enhanced removal of 1,4-dioxane, an important environmental pollutant, from water. Increasing the concentration of end-capping monomers produces PAFs with narrower pore size distributions and improved 1,4-dioxane uptake. These results illustrate that defect engineering can effectively modulate polymer connectivity and porosity for applications in selective adsorptive separations. In conclusion, this technique avoids post-synthetic treatments and presents another approach to tailor amorphous polymeric adsorbents.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evidence of a Uranium‐Paddlewheel Node in a Catecholate‐Based Metal–Organic Framework

Abstract The interactions between uranium and non‐innocent organic species are an essential component of fundamental uranium redox chemistry. However, they have seldom been explored in the context of multidimensional, porous materials. Uranium‐based metal–organic frameworks (MOFs) offer a new angle to study these interactions, as these self‐assembled species stabilize uranium species through immobilization by organic linkers within a crystalline framework, while potentially providing a method for adjusting metal oxidation state through coordination of non‐innocent linkers. We report the synthesis of the MOF NU‐1700 , assembled from U 4+ ‐paddlewheel nodes and catecholate‐based linkers. We propose this highly unusual structure, which contains two U 4+ ions in a paddlewheel built from four linkers—a first among uranium materials—as a result of extensive characterization via powder X‐ray diffraction (PXRD), sorption, transmission electron microscopy (TEM), and thermogravimetric analysis (TGA), in addition to density functional theory (DFT) calculations.

Knapp, Julia G.↗
DOE OSTI2023

Evidence of a Uranium‐Paddlewheel Node in a Catecholate‐Based Metal–Organic Framework

Abstract The interactions between uranium and non‐innocent organic species are an essential component of fundamental uranium redox chemistry. However, they have seldom been explored in the context of multidimensional, porous materials. Uranium‐based metal–organic frameworks (MOFs) offer a new angle to study these interactions, as these self‐assembled species stabilize uranium species through immobilization by organic linkers within a crystalline framework, while potentially providing a method for adjusting metal oxidation state through coordination of non‐innocent linkers. We report the synthesis of the MOF NU‐1700 , assembled from U 4+ ‐paddlewheel nodes and catecholate‐based linkers. We propose this highly unusual structure, which contains two U 4+ ions in a paddlewheel built from four linkers—a first among uranium materials—as a result of extensive characterization via powder X‐ray diffraction (PXRD), sorption, transmission electron microscopy (TEM), and thermogravimetric analysis (TGA), in addition to density functional theory (DFT) calculations.

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

Computational Screening of Supported Metal Oxide Nanoclusters for Methane Activation: Insights into Homolytic versus Heterolytic C–H Bond Dissociation

Since its discovery in zeolites, the [CuOCu] 2+ motif has played an important role in our understanding of selective methane activation over supported metal oxide nanoclusters. Although there are two known C-H bond dissoci-ation mechanisms, namely homolytic and heterolytic cleavage, most computational studies on optimizing metal oxide nanoclusters for improved methane activation reactivity have focused only on the homolytic mechanism. In this work, both mechanisms were examined for a set of 21 mixed metal oxide complexes of the form of [M 1 OM 2 ] 2+ (M 1 , M 2 = Mn, Fe, Co, Ni, Cu, Zn). Except for pure copper, heterolytic cleavage was found to be the dominant C-H bond activation pathway for all systems. Furthermore, mixed systems including [CuOMn] 2+ , [CuONi] 2+ , and [CuOZn] 2+ are predicted to possess similar methane activation activity as pure [CuOCu] 2+ . Furthermore, these results suggest that both homolytic and heterolytic mechanisms should be considered in computing methane activation energies on supported metal oxide nanoclusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗