Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lateral interactions”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Determining Catalytically Relevant Surfaces through Coverage-Dependent Lattice Gas Models: Carbon Adsorption on Fe(100)

Here, we have quantified the C–C lateral interactions on Fe(100) using a density functional theory (DFT)-parameterized lattice gas cluster expansion (LG CE) model trained using 265 unique configurations spanning a C coverage from 0 to 1 monolayer (ML). Our LG CE model shows high predictive accuracy with a leave-multiple-out cross-validation score of 10.2 and 16.6 meV/site for systems with and without the top two layers of Fe atoms fixed, respectively. Electronic ground-state structures identified from the lattice gas model (including the structures at 0 and 1 monolayers) were further used to generate ab initio phase diagrams under a range of temperatures and pressures. At low temperatures (<400 K), we found that the 1.0 monolayer structure is dominant, whereas at higher temperatures (>500 K), the 0.88 ML structure is most likely to form on the Fe surface. Interestingly, our model identified a c (2 × 2) ordered structure at 1/2 ML, which correlates well with previous DFT studies for carbon adsorption on iron surfaces and matches with the experimentally observed low-energy electron diffraction structure. Overall, the DFT-parameterized energies for the C/Fe system including effects of coverage and configurational space can further help in developing multiscale models for various heterogeneous reactions involving C–C and C–Fe interactions.

08 HYDROGEN↗
DOE OSTI2019

Formulation of Multicomponent Lattice Gas Model Cluster Expansions Parameterized on Ab Initio Data: An Introduction to the Ab Initio Mean-Field Augmented Lattice Gas Modeling Code

A formalism is presented for the construction of multicomponent lattice gas models parameterized with ab initio (typically density functional theory) data. The leave-multiple-out and the leave-one-out cross validation score are showcased, and practical algorithms are developed and implemented in a new code called the Ab initio Mean-field Augmented Lattice Gas Modeling (AMALGM) code. Here, the functionality of these algorithms is demonstrated with a fully worked out example using the O/Fe(100) system. AMALGM and the formalism on which it is based are envisioned as a surface-oriented lattice gas alternative to other cluster expansion codes that are typically geared toward bulk systems whose lateral interactions between components are parameterized using an Ising model. While the formalism is created within the context of surfaces, it is equally applicable to bulk materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Incorporating Coverage-Dependent Reaction Barriers into First-Principles-Based Microkinetic Models: Approaches and Challenges

Mean-field microkinetic models (MKMs) are appealing for their relatively facile construction, computational tractability, and high-throughput catalyst screening capabilities. As such, they will continue to be a valuable tool for materials design in heterogeneous catalysis even as the field aims to describe more complex systems. Numerous prior reports have provided the groundwork for constructing first-principles-based MKMs, including the analysis of strategies for incorporating lateral interactions into thermodynamic parameters (e.g., adsorption energies). Yet, there remains a need for concerted dialogue on methods for calculating and incorporating coverage-dependent kinetic parameters into MKMs. In this Perspective, we assess strategies for doing so, including the corresponding key physical implications and computational challenges. Here, we emphasize that decoupling thermodynamic and kinetic parameters within MKMs can violate thermodynamic consistency and risk unphysical solutions. For some reactions and catalyst materials, scaling relationships can predict coverage-dependent activation energies, but there are several exceptions evident in the literature, indicating that this approach is not universally applicable and that the field could benefit from research aimed at elucidating the limitations. Conducting high-coverage transition state searches is a rigorous but computationally costly strategy, and the effects of various methods for mitigating this cost on resulting energetics have yet to be broadly explored and validated. The goal of this Perspective is to generate discussion on and inspire focused research into the physical relevance of approaches for describing coverage-dependent reaction barriers in MKMs, including the development of computationally tractable methodologies, to advance the applicability of MKMs across diverse reaction chemistries and conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Electrons to Reactors Multiscale Modeling: Catalytic CO Oxidation over RuO 2

First-principles kinetic Monte Carlo (1p-kMC) simulations for CO oxidation on two RuO 2 facets, RuO 2 (110) and RuO 2 (111), were coupled to the computational fluid dynamics (CFD) simulations package MFIX, and reactor-scale simulations were then performed. 1p-kMC coupled with CFD has recently been shown as a feasible method for translating molecular scale mechanistic knowledge to the reactor scale, enabling comparisons to in situ and online experimental measurements. Only a few studies with such coupling have been published. This work incorporates multiple catalytic surface facets into the scale-coupled simulation, and three possibilities were investigated: the two possibilities of each facet individually being the dominant phase in the reactor, and also the possibility that both facets were present on the catalyst particles in the ratio predicted by an ab initio thermodynamics-based Wulff construction. When lateral interactions between adsorbates were included in the 1p-kMC simulations, the two surfaces, RuO 2 (110) and RuO 2 (111), were found to be of similar order-of-magnitude in activity for the pressure range of 1 × 10 –4 bar to 1 bar, with the RuO 2 (110) surface-termination showing more simulated activity than the RuO 2 (111) surface-termination. Coupling between the 1p-kMC and CFD was achieved with a lookup table generated by the error-based modified Shepard interpolation scheme. Isothermal reactor scale simulations were performed and compared to two separate experimental studies, conducted with reactant partial pressures of ≤0.1 bar. Simulations without an isothermality restriction were also conducted and showed that the simulated temperature gradient across the catalytic reactor bed is <0.5 K, which validated the use of the isothermality restriction for investigating the reactor-scale phenomenological temperature dependences. The approach with the Wulff construction based reactor simulations reproduced a trend similar to one experimental data set relatively well, with the (110) surface being more active at higher temperaures; in contrast, for the other experimental data set, our reactor simulations achieve surprisingly and perhaps fortuitously good agreement with the activity and phenomenological pressure dependence when it is assumed that the (111) facet is the only active facet present. Lastly, the active phase of catalytic CO oxidation over RuO 2 remains unsettled, but the present study presents proof of principle (and progress) toward more accurate multiscale modeling from electrons to reactors and new simulation results.

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

Structural Diversity in the Complexes of Trimeric Perfluoro- o -phenylene Mercury with Tetrathia- and Tetramethyltetraselenafulvalene

Five potential charge transfer complexes of trimeric perfluoro-o-phenylene mercury (I) with tetrathiafulvalene (TTF) and tetramethyltetraselenefulvalene (TMTSF) were grown from different solvent mixtures. The adducts (I)2·TTF (1) and I·TTF (2) were grown by slow evaporation from the 1:1 mixture of dichloromethane (CH 2 Cl 2 , DCM) and carbon disulfide (CS 2 ). Use of the different 1:1 solvent mixtures of dichloromethane (CH 2 Cl 2 , DCM) and dichloroethane (C 2 H 4 Cl 2 , DCE) has led to the crystalline adducts I·TTF (3) and I·TTF·DCE (4). Adduct I.TMTSF (5) was grown by the interface crystallization on the border of two immiscible layers, ethyl acetate, and carbon disulfide. The cocrystals differ by the donor–acceptor ratio, molecular packing, and the solvent inclusion. The components in 1–5 form mixed donor–acceptor stacks. The stacks are stabilized by Hg···S and Hg···C short contacts, while the lateral interactions between stacks include F···F, CH···F, and S/Se···F short contacts.

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

Coupling of kinetic Monte Carlo simulations of surface reactions to transport in a fluid for heterogeneous catalytic reactor modeling

We have developed a method to couple kinetic Monte Carlo simulations of surface reactions at a molecular scale to transport equations at a macroscopic scale. This method is applicable to steady state reactors. We use a finite difference upwinding scheme and a gap-tooth scheme to efficiently use a limited amount of kinetic Monte Carlo simulations. In general the stochastic kinetic Monte Carlo results do not obey mass conservation so that unphysical accumulation of mass could occur in the reactor. We have developed a method to perform mass balance corrections that is based on a stoichiometry matrix and a least-squares problem that is reduced to a non-singular set of linear equations that is applicable to any surface catalyzed reaction. The implementation of these methods is validated by comparing numerical results of a reactor simulation with a unimolecular reaction to an analytical solution. Furthermore, the method is applied to two reaction mechanisms. The first is the ZGB model for CO oxidation in which inevitable poisoning of the catalyst limits the performance of the reactor. The second is a model for the oxidation of NO on a Pt(111) surface, which becomes active due to lateral interaction at high coverages of oxygen. This reaction model is based on ab initio density functional theory calculations from literature.

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

Structure of chalcogen overlayers on Au(111): Density functional theory and lattice-gas modeling

Ordering of different chalcogens, S, Se, and Te, on Au(111) exhibit broad similarities but also some distinct features, which must reflect subtle differences in relative values of the long-range pair and many-body lateral interactions between adatoms. We develop lattice-gas (LG) models within a cluster expansion framework, which includes about 50 interaction parameters. These LG models are developed based on density functional theory (DFT) analysis of the energetics of key adlayer configurations in combination with the Monte Carlo (MC) simulation of the LG models to identify statistically relevant adlayer motifs, i.e., model development is based entirely on theoretical considerations. The MC simulation guides additional DFT analysis and iterative model refinement. Given their complexity, development of optimal models is also aided by strategies from supervised machine learning. The model for S successfully captures ordering motifs over a broader range of coverage than achieved by previous models, and models for Se and Te capture the features of ordering, which are distinct from those for S. More specifically, the modeling for all three chalcogens successfully explains the linear adatom rows (also subtle differences between them) observed at low coverages of ~0.1 monolayer. The model for S also leads to a new possible explanation for the experimentally observed phase with a (5 × 5)-type low energy electron diffraction (LEED) pattern at 0.28 ML and to predictions for LEED patterns that would be observed with Se and Te at this coverage.

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

Irreversible oxygen poisoning: Modeling the initial water dissociation kinetics on δ-Pu(111) and δ-Pu(100) through density functional theory

In this work, the initial kinetics of water dissociation on two facets of δ-plutonium, δ-Pu(111) and δ-Pu(100), are explored through density functional theory in order to understand how water dissociation occurs on these facets. We explored the dissociation of water via the formation of hydroxyls, atomic hydrogen, and atomic oxygen species on each facet. We calculate low energetic barriers for water to split to adsorbed hydrogen and hydroxyl species at 0.19 eV for δ-Pu(111) and 0.07 eV for δ-Pu(100). The hydroxyl has a barrier of 0.64 and 0.37 eV to cleave the hydrogen–oxygen bond on δ-Pu(111) and δ-Pu(100), respectively. Due to the highly exergonic adsorption free energy of atomic oxygen of −2.10 eV, the metallic surfaces are found to be fully covered in oxygen, even with the inclusion of oxygen lateral interactions. When combined with the reaction thermodynamics, this free energy of adsorption forms a molecular oxygen desorption barrier greater than 9.51 eV (918 kJ/mol). These results, combined with simulated temperature programed desorption spectra, indicate that oxygen formed via water dissociation on δ-Pu (111) and (100) facets induces an irreversible poison that prevents further reaction of water directly on metallic plutonium surfaces, most likely due to the strong hybridization of the Pu and O valence states. Therefore, these results imply that another mechanism is responsible for the continuously experimentally measured water dissociation, which is most likely due to the Pu oxide.

36 MATERIALS SCIENCE↗
DOE OSTI2023

KinCat v.1.0

SAND2024-02099O The software is designed to allow researchers to perform kinetic Monte Carlo (KMC) simulations of catalytic reactions on a 2D lattice. The code is written efficiently to run on a variety of shared memory computing architectures (e.g. GPU, multi-core) and to natively express the full complexity of lateral interactions on reaction rates. The software allows researchers to perform KMC simulations of catalytic reactions on a 2D lattice. It uses parallel shared-memory computing architectures to reduce run-times and allows for simultaneous simulation of multiple independent runs. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Najm, Habib↗
DOE OSTI2015

Sterile Neutrino Search with MINOS

MINOS, Main Injector Neutrino Oscillation Search, is a long-baseline neutrino oscillation experiment in the NuMI muon neutrino beam at the Fermi National Accelerator Laboratory in Batavia, IL. It consists of two detectors, a near detector positioned 1 km from the source of the beam and a far detector 734 km away in Minnesota. MINOS is primarily designed to observe muon neutrino disappearance resulting from three flavor oscillations. The Standard Model of Particle Physics predicts that neutrinos oscillate between three active flavors as they propagate through space. This means that a muon-type neutrino has a certain probability to later interact as a different type of neutrino. In the standard picture, the neutrino oscillation probabilities depend only on three neutrino flavors and two mass splittings, Δm 2 . An anomaly was observed by the LSND and MiniBooNE experiments that suggests the existence of a fourth, sterile neutrino flavor that does not interact through any of the known Standard Model interactions. Oscillations into a theoretical sterile flavor may be observed by a deficit in neutral current interactions in the MINOS detectors. A distortion in the charged current energy spectrum might also be visible if oscillations into the sterile flavor are driven by a large mass-squared difference, m s 2 ~ 1 eV 2 . The results of the 2013 sterile neutrino search are presented here.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1997

Optical and Atomic Force Microscopy Characterization of PbI2 Quantum Dots

Lead iodide (PbI2) clusters were synthesized from the chemical reaction of NaI (or KI) with Pb(NO3)2 in H2O, D2O, CH3OH, and C3H7OH media. The observation of the absorption features above 350 nm with the help of integrating sphere accessory strongly suggests the quantum dot formation of PbI2 in solution. Spectral comparison between the synthesized PbI2 clusters in solution and PbI2 nanophase by impregnation of PbI2 in four different pore-sized porous silica indicates that the PbI2 cluster size in solution is less than 2.5 nm in lateral dimension. Atomic force microscopy (AFM) measurements show that the PbL clusters deposited onto three different molecularly flat surfaces are single-layered. The measured height is 1.0 - 0.1 nm. The swollen layer thickness can be attributed to the intralayer contraction from the strong lateral interaction among PbI2 molecules, which is supported by ab initio calculation. Raman scattering measurement of LO and TO modes of PbI2 in bulk and in the confined state were also conducted in 50-150 cu cm region. The observed three bands at 74, %, 106 1/cm are assigned to TO2, LO2, and LO, mode, respectively. The relatively small red-shift in LO modes may be caused by the surface phonon polaritons of PbI2 nanophase in the porous silica.

Mu, R.↗
NASA NTRS2013

A Panchromatic View of the Restless SN 2009ip Reveals the Explosive Ejection of a Massive Star Envelope

The double explosion of SN 2009ip in 2012 raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN 2009ip during its remarkable rebrightenings. High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the Very Large Array, Swift, Fermi, Hubble Space Telescope, and XMM) constrain SN 2009ip to be a low energy (E (is) approximately 10(exp 50) erg for an ejecta mass approximately 0.5 solar mass) and asymmetric explosion in a complex medium shaped by multiple eruptions of the restless progenitor star. Most of the energy is radiated as a result of the shock breaking out through a dense shell of material located at approximately 5 × 10( exp 14) cm with M approximately 0.1 solar mass, ejected by the precursor outburst approximately 40 days before the major explosion. We interpret the NIR excess of emission as signature of material located further out, the origin of which has to be connected with documented mass-loss episodes in the previous years. Our modeling predicts bright neutrino emission associated with the shock break-out if the cosmic-ray energy is comparable to the radiated energy. We connect this phenomenology with the explosive ejection of the outer layers of the massive progenitor star, which later interacted with material deposited in the surroundings by previous eruptions. Future observations will reveal if the massive luminous progenitor star survived. Irrespective of whether the explosion was terminal, SN 2009ip brought to light the existence of new channels for sustained episodic mass loss, the physical origin of which has yet to be identified.

view↗
NASA NTRS2013

A Panchromatic View of the Restless SN 2009ip Reveals the Explosive Ejection of a Massive Star Envelope

The double explosion of SN 2009ip in 2012 raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN 2009ip during its remarkable rebrightenings. High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the Very Large Array, Swift, Fermi, Hubble Space Telescope, and XMM) constrain SN 2009ip to be a low energy (E approximating 10(exp 50) ergs for an ejecta mass approximating 0.5 M solar mass) and asymmetric explosion in a complex medium shaped by multiple eruptions of the restless progenitor star. Most of the energy is radiated as a result of the shock breaking out through a dense shell of material located at approximately 5 times 10 (exp 14) cm with M approximating 0.1 solar mass, ejected by the precursor outburst approximately 40 days before the major explosion. We interpret the NIR (Near Infrared) excess of emission as signature of material located further out, the origin of which has to be connected with documented mass-loss episodes in the previous years. Our modeling predicts bright neutrino emission associated with the shock break-out if the cosmic-ray energy is comparable to the radiated energy. We connect this phenomenology with the explosive ejection of the outer layers of the massive progenitor star, which later interacted with material deposited in the surroundings by previous eruptions. Future observations will reveal if the massive luminous progenitor star survived. Irrespective of whether the explosion was terminal, SN 2009ip brought to light the existence of new channels for sustained episodic mass loss, the physical origin of which has yet to be identified.

massive↗
NASA NTRS1985

Transverse momentum distribution of pion deg in the fragmentation region of super high energy interactions

The lateral distribution of gamma-families observed by emulsion chamber is sensitive to test transverse momentum of high energy interaction. But most part gamma-families are succesive interation's results which makes it necessary to analyze the propagation of gamma-ray in atmosphere. A gamma-ray with energy Er and transverse momentum Pt is produced at the altitude h. In the emulsion chamber experiment, Eob is estimated by decascade method. There are 30 gamma-families observed by Ganbala emulsion Chamber. A group of Monte-Carlo simulation gamma-families are used to compare with experimental data with the same treatment. It is found that both distributions are consistent.

Ren, J. R.↗
DOE OSTI2021

Catalytic site ensembles: A context to reexamine the Langmuir-Hinshelwood kinetic description

The Langmuir-Hinshelwood formalism describes catalytic reactions of Langmuirian surface species under the assumption that all adsorbates are randomly-distributed – enabling adjacency of surface-bound intermediates to be determined solely by coverages of single-site occupants. We demonstrate herein that this approximation is inappropriate even for simple catalytic reactions (e.g. A + A → A 2 ) and manifestly neglects islanding of slowly-consumed species and partitioning of highly-reactive species inherently engendered by ≥ two-site elementary steps (e.g. A*–A* → A 2(g) + *–*). Rigorous description of kinetically-consequential islanding/partitioning phenomena requires explicit description of the coverage and chemical dynamics of all multi-site ensembles. Higher-order, ensemble-specific rate terms identify the particular microscopic events relevant to each ensemble, and, in doing so, reveal that each elementary step (e.g. A (g) adsorption) describes not one event (e.g. A (g) + * → A*) , but a sum over all ensemble-specific paths (e.g. A (g) + *–* → A*–* and A (g) + A*–* → A*–A*). De-convoluting each elementary step into its constituent multi-site paths proffers kinetic detail otherwise inaccessible – enabling (i) identification of rate- and selectivity-determining site ensembles, (ii) calculation of rates and degrees of rate control of ensemble-specific elementary steps, (iii) incorporation of adsorbate surface diffusion, (iv) incorporation of lateral adsorbate interactions, and (v) quantitative description of catalysis of multi-site-occupying intermediates (e.g. *C n H m * species in hydrocarbon (de-)hydrogenation and C-C bond coupling/cleavage reactions) which we demonstrate here is inaccessible to the Langmuir-Hinshelwood formalism even if adsorbate surface diffusion is infinitely-fast.

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

Investigation of shock-shock interaction and Mach reflection in laterally colliding laser-blow-off plasmas

Interactions of two Li plasma plumes and shock waves are investigated at various pressures (∼10{sup −5} to 3 mbar) in the argon gas ambient. Fast imaging and optical emission spectroscopy are used to study the plume dynamics and characteristic emission of plasmas. The plasma plumes are created in laser-blow-off geometry. The expansion of plasma plumes in the ambient gas leads to the formation of an interaction zone. The formation of interaction zone is dependent on the ambient pressure and below a certain pressure, no significant change is observed in the shape and size of the interaction plasma. In the higher pressure, formation of interaction zone and its shape are dependent on ambient pressure. Dynamics of seed plasmas and interaction zone are also affected by the shock-shock interactions. The shock-shock interaction depends on the angle of incidence (α) between two shock waves at the initial time of interaction but as the plumes expand, the shock-shock interaction does not follow α dependence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2009

Interaction of Aircraft Wakes From Laterally Spaced Aircraft

Large Eddy Simulations are used to examine wake interactions from aircraft on closely spaced parallel paths. Two sets of experiments are conducted, with the first set examining wake interactions out of ground effect (OGE) and the second set for in ground effect (IGE). The initial wake field for each aircraft represents a rolled-up wake vortex pair generated by a B-747. Parametric sets include wake interactions from aircraft pairs with lateral separations of 400, 500, 600, and 750 ft. The simulation of a wake from a single aircraft is used as baseline. The study shows that wake vortices from either a pair or a formation of B-747 s that fly with very close lateral spacing, last longer than those from an isolated B-747. For OGE, the inner vortices between the pair of aircraft, ascend, link and quickly dissipate, leaving the outer vortices to decay and descend slowly. For the IGE scenario, the inner vortices ascend and last longer, while the outer vortices decay from ground interaction at a rate similar to that expected from an isolated aircraft. Both OGE and IGE scenarios produce longer-lasting wakes for aircraft with separations less than 600 ft. The results are significant because concepts to increase airport capacity have been proposed that assume either aircraft formations and/or aircraft pairs landing on very closely spaced runways.

Proctor, Fred H.↗