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 19 records

Predicting the Electric Field Effect on the Lateral Interactions Between Adsorbates: O/Fe(100) from First Principles

A density functional theory parameterized lattice gas cluster expansion of oxygen on an Fe(100) surface is developed in the presence and absence of an applied positive and negative external electric field, characterizing the heterogeneity in oxygen's surface distribution and the effect of an external electric field on the lateral interactions between the adsorbates. We show that the presence of a negative electric field tends to weaken both attractive and repulsive interactions while the positive electric field tends to strengthen these interactions, altering surface distributions and ground state configurations. Since lateral interactions have been shown to play a critical role in defining catalytic behavior, the application of an applied electric field has the potential to be a useful tool in adjusting critical chemical properties of Fe-based hydrodeoxygenation catalysts, by decreasing the repulsive interactions between the adspecies in the presence of a negative field and thereby mitigating the formation of an oxide.

Bray, Jacob E.↗

Root lateral interactions drive water uptake patterns under water limitation

The plasticity of root water uptake determines the maintenance of transpiration during periods of water limitation and drought. However, the mechanistic basis of plant water uptake, as well as the implications of water uptake strategies at the individual and ecosystem scale remain elusive. We model three-dimensional root water uptake under variably saturated conditions for a one-hectare temperate forest plot for a growing season with a pronounced mid-season dry period. Variations in root architecture, hydraulic properties, and degree of lateral interaction between root systems produce divergent local responses to water limitation and provide insights on individual and community response to meteorological conditions. Results demonstrate the plasticity of ecosystem-scale responses to surface drying, where interacting roots shift regions of active uptake to deeper soil layers with less abundant root biomass. These shifts, a product of both root system and soil hydraulic properties, illustrate intimate links between root and soil hydraulics in determining plant water sourcing. Furthermore, we demonstrate that root lateral interactions are beneficial at the ecosystem-scale, even when trees compete for water. Specifically, a more spatially extensive root system facilitates access to a larger soil water reservoir, often ameliorating water limitation and reducing sharp water potential gradients. While the reduction of water stress is a benefit, it can be offset by increased root construction and maintenance costs associated with the larger rooting system. A plausible “viability” region of root communal co-existence is therefore implied where competitive pressures and root production costs are balanced by a potential water benefit.

54 ENVIRONMENTAL SCIENCES↗

Evaluating the benefits of kinetic Monte Carlo and microkinetic modeling for catalyst design studies in the presence of lateral interactions

Popular computational catalyst design strategies rely on the identification of reactivity descriptors, which can be used along with Brønsted–Evans–Polanyi (BEP) and scaling relations as input to a microkinetic model (MKM) to make predictions for activity or selectivity trends. The main benefit of this approach is related to the inherent dimensionality reduction of the large material space to just a few catalyst descriptors. Conversely, it is well documented that a small set of descriptors is insufficient to capture the intricacies and complexities of a real catalytic system. The inclusion of coverage effects through lateral adsorbate-adsorbate interactions can narrow the gap between simplified descriptor predictions and real systems, but mean-field MKMs cannot properly account for local coverage effects. This shortcoming of the mean-field approximation can be rectified by switching to a lattice-based kinetic Monte Carlo (kMC) method using cluster expansion representation of adsorbate–adsorbate lateral interactions. Using the prototypical CO oxidation reaction as an example, we critically evaluate the benefits of kMC over MKM in terms of trend predictions and computational cost when using only a small set of input parameters. After confirming that in the absence of lateral interactions the kMC and MKM approaches yield identical trends and mechanistic information, we observed substantial differences between the two kinetic models when lateral interactions were introduced. Here, the mean-field implementation applies coverage corrections directly to the descriptors, causing an artificial overprediction of the activity of strongly binding metals. In contrast, the cluster expansion in kMC implementation can differentiate among the highly active metals but it is very sensitive to the set of included interaction parameters. Considering that computational screening relies on a minimal set of descriptors, for which MKM makes reasonable trend predictions at a ca. three orders of magnitude lower computational cost than kMC, the MKM approach does provide a better entry point for computational catalyst design.

36 MATERIALS SCIENCE↗

Parallel kinetic Monte Carlo simulation framework incorporating accurate models of adsorbate lateral interactions

Ab initio kinetic Monte Carlo (KMC) simulations have been successfully applied for over two decades to elucidate the underlying physico-chemical phenomena on the surfaces of heterogeneous catalysts. These simulations necessitate detailed knowledge of the kinetics of elementary reactions constituting the reaction mechanism, and the energetics of the species participating in the chemistry. The information about the energetics is encoded in the formation energies of gas and surface-bound species, and the lateral interactions between adsorbates on the catalytic surface, which can be modeled at different levels of detail. The majority of previous works accounted for only pairwise-additive first nearest-neighbor interactions. More recently, cluster-expansion Hamiltonians incorporating long-range interactions and many-body terms have been used for detailed estimations of catalytic rate [C. Wu, D. J. Schmidt, C. Wolverton, and W. F. Schneider, J. Catal. 286, 88 (2012)]. In view of the increasing interest in accurate predictions of catalytic performance, there is a need for general-purpose KMC approaches incorporating detailed cluster expansion models for the adlayer energetics. We have addressed this need by building on the previously introduced graph-theoretical KMC framework, and we have developed Zacros, a FORTRAN2003 KMC package for simulating catalytic chemistries. To tackle the high computational cost in the presence of long-range interactions we introduce parallelization with OpenMP. We further benchmark our framework by simulating a KMC analogue of the NO oxidation system established by Schneider and co-workers [J. Catal. 286, 88 (2012)]. We show that taking into account only first nearest-neighbor interactions may lead to large errors in the prediction of the catalytic rate, whereas for accurate estimates thereof, one needs to include long-range terms in the cluster expansion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Desorption of Benzene, 1,3,5-Trifluorobenzene, and Hexafluorobenzene from a Graphene Surface: The Effect of Lateral Interactions on the Desorption Kinetics

The desorption of benzene, 1,3,5-trifluorobenzene (TFB), and hexafluorobenzene (HFB) from a graphene covered Pt(111) substrate was investigated using temperature programmed desorption (TPD). All three species have well-resolved monolayer and second layer desorption peaks. The desorption spectra for submonolayer coverages of benzene and hexafluorobenzene are consistent with first-order desorption kinetics. In contrast, the submonolayer TPD spectra for 1,3,5- trifluorobenzene align on a common leading-edge which is indicative of zero-order desorption kinetics. The desorption behavior of the three molecules can be correlated with the strength of the quadrupole moments. Calculations (second-order Møller-Plesset perturbation and density functional theory) show that the potential minimum for coplanar TFB dimers is more than a factor of two greater than that for either benzene or HFB dimers. The calculations support the interpretation that benzene and HFB are less likely to form the two-dimensional islands that are needed for submonolayer zero-order desorption kinetics.

Smith, Ronald S.↗

Benchmarking the accuracy of coverage-dependent models: adsorption and desorption of benzene on Pt (111) and Pt 3 Sn (111) from first principles

Bimetallic catalysts have demonstrated properties favorable for upgrading biofuel through catalytic hydrodeoxygenation. However, the design and optimization of such bimetallic catalysts requires the ability to construct accurate, predictive models of these systems. To generate a model that predicts the kinetic behavior of benzene adsorbed on Pt (111) and a Pt 3 Sn (111) surface alloy (Pt 3 Sn (111)), the adsorption of benzene was studied for a wide range of benzene coverages on both surfaces using density functional theory (DFT) calculations. The adsorption energy of benzene was found to correlate linearly with benzene coverage on Pt (111) and Pt 3 Sn (111); both surfaces exhibited net repulsive lateral interactions. Through an analysis of the d-band properties of the metal surface, it was determined that the coverage dependence is a consequence of the electronic interactions between benzene and the surface. The linear coverage dependence of the adsorption energy allowed us to quantify the influence of the lateral interactions on the heat of adsorption and temperature programmed desorption (TPD) spectra using a mean-field model. A comparison of our simulated TPD to experiment showed that this mean-field model adequately reproduces the desorption behavior of benzene on Pt (111) and Pt 3 Sn (111). In particular, the TPD correctly exhibits a broadening desorption peak as the initial coverage of benzene increases on Pt (111) as well as and a low temperature desorption peak on Pt 3 Sn (111). Furthermore, due to the sensitivity of the TPD peak temperature to the desorption energy, precise alignment of experimental and theoretical TPD spectra demands an accurate calculation of the adsorption energy. Therefore, an analysis of the effect of the exchange-correlation functional on TPD modeling is presented. Through this work, we show the necessity of incorporating lateral interactions into theoretical models in order to correctly predict experimental behavior.

09 BIOMASS FUELS↗

On factors controlling activity of submonolayer bimetallic catalysts: Nitrogen desorption

We model N{sub 2} desorption on submonolayer bimetallic surfaces consisting of Co clusters on Pt(111) via first-principles density functional theory-based kinetic Monte Carlo simulations. We find that submonolayer structures are essential to rationalize the high activity of these bimetallics in ammonia decomposition. We show that the N{sub 2} desorption temperature on Co/Pt(111) is about 100 K higher than that on Ni/Pt(111), despite Co/Pt(111) binding N weaker at low N coverages. Co/Pt(111) has substantially different lateral interactions than single metals and Ni/Pt. The lateral interactions are rationalized with the d-band center theory. The activity of bimetallic catalysts is the result of heterogeneity of binding energies and reaction barriers among sites, and the most active site can differ on various bimetallics. Our results are in excellent agreement with experimental data and demonstrate for the first time that the zero-coverage descriptor, used until now, for catalyst activity is inadequate due not only to lacking lateral interactions but importantly to presence of multiple sites and a complex interplay of thermodynamics (binding energies, occupation) and kinetics (association barriers) on those sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unifying principles for catalytic hydrotreating processes (Final Technical Report)

This project builds on the hypothesis that the hydrotreating processes for the removal of oxygen and sulfur are fundamentally similar at the atomic-scale and existing knowledge from the treatment of petroleum derived feedstock can be leveraged for the design of novel catalysts for the upgrade of bio-oil. We tested this hypothesis by comparing computed potential energy diagrams for hydrodesulfurization (HDS) of thiophene over MoS 2 with hydrodeoxygenation (HDO) of furan over MoO 3 and concluded that certain aspects, such as catalyst promotion with transition metals, are valid strategies for both reactions. On the other hand, we also noticed significant differences in the mechanism for hydrogen (H 2 ) activation, which requires sites with metallic character. While MoS 2 is known to have metallic edge states that can catalyze H 2 dissociation, this elementary step is prohibitively slow on defect-free oxides. Only in the presence of vacancies or by creating metal/oxide interfaces can efficient H 2 activation sites during HDO be formed. The need for bifunctional catalyst when it comes to efficient and selective HDO or dehydrogenation reactions was further corroborated in joint experimental and theoretical studies of the Guerbet reaction for the coupling of biomass derived oxygenates over PdCu alloys, nitrate reduction over In-promoted Pd nanoparticles, and ethylene dehydroaromatization over Ga-exchanged ZSM-5 zeolites. All of these catalytic systems have in common that catalytic sites with distinct functional requirements are needed to create a working catalyst. Detailed computational studies were carried out for HDO of m -cresol and phenol on Ru-modified TiO 2 surfaces, which allowed us to attribute catalytic activity to the metal/oxide interface. A surprising finding was that heterolytic cleavage of the H-H bond across the Ru/TiO 2 interface was critically important, despite lower barriers for homolytic H 2 activation on Ru metal. The explanation lies in the high barriers for hydrogen spillover from Ru to TiO 2 , which becomes unnecessary in the heterolytic activation pathway. Moreover, we also reported that proton transfer steps between metal and oxide sites are mediated by weakly adsorbed surface water. During attempts to develop and validate a kinetic Monte Carlo (kMC) model for HDO reactions at the Ru/TiO 2 interface, it became clear that lateral interactions are paramount to describe realistic surface chemistry and without these interactions, the reduction and hydroxylation behavior from our simulations was inconsistent with reported experiments. To assess the importance of lateral interactions in popular computational catalyst design strategies relying on the identification of reactivity descriptors, which can be used along with Brønsted–Evans–Polanyi (BEP) and scaling relations as input to a microkinetic model (MKM) to make predictions for activity or selectivity trends, we compared predicted trends with those obtained from descriptor-based kMC models. We critically evaluated the benefits of kMC over MKM in terms of trend predictions and computational cost when using only a small set of input parameters. After confirming that in the absence of lateral interactions the kMC and MKM approaches yield identical trends and mechanistic information, we observed substantial differences between the two kinetic models when lateral interactions were introduced. The mean-field implementation applies coverage corrections directly to the descriptors, causing an artificial overprediction of the activity of strongly binding metals. In contrast, the cluster expansion in kMC implementation can differentiate among the highly active metals but it is very sensitive to the set of included interaction parameters. Considering that computational screening relies on a minimal set of descriptors, for which MKM makes reasonable trend predictions at a ca. three orders of magnitude lower computational cost than kMC, we concluded that the MKM approach does provide an overall better entry point for computational catalyst design. Overall, this project has led to 11 peer-reviewed publications, and their scientific is impact is well illustrated by their combined 702 citations.

08 HYDROGEN↗

Atomistic and phase field simulations of three dimensional interactions of {$10\bar12$} twins with grain boundaries in Mg: twin transmission and dislocation emission

In polycrystals, the interaction of dislocations and twins with grain boundaries (GBs) plays a role in hardening and formability during plastic deformation. While dislocation-GB interactions are relatively well-understood, twin-GB interactions remain mostly unknown. In this work, an approach using molecular dynamics and phase-field simulations is followed to study the forward and lateral interactions between {$10\bar12$} twins and tilt grain boundaries in Mg. Molecular dynamics results show that the resolved shear stress on slip/twinning modes of the neighboring grain, not the geometric alignment, is the dominant factor in determining the outcome of the twin-GB interactions. For some lateral interaction configurations, as the misorientation angle increases, the resolved shear stress on the same {$10\bar12$} twin variant of the neighboring grain reduces while it increases for slip or I 2 stacking fault emissions or other twin modes such as {$11\bar12$} and {$10\bar11$}, explaining why twin transmission is not seen at high misorientation angles. Furthermore, lateral and forward interactions of the twin with tilt grain boundaries whose misorientation axes are normal to the coherent twin boundary show significantly different outcomes. For the forward interaction, the twin is absorbed and stacking faults are emitted when interacting for low misorientation angles (up to 30°) while the lateral interaction results in twin transmission, nucleation of a {$11\bar12$} twin, and emission of I 2 stacking faults. Finally, comparisons between twin interactions with symmetric and asymmetric tilt GBs with different GB structures show similar outcomes.

36 MATERIALS SCIENCE↗

Coverage-Dependent Adsorption of Hydrogen on Fe(100): Determining Catalytically Relevant Surface Structures via Lattice Gas Models

Hydrogen adatoms are a critical surface species for several reactions catalyzed by Fe surfaces such as Fischer–Tropsch synthesis, ammonia synthesis, and the hydrodeoxygenation of biomass-derived oxygenates. Parameterizing the energetics for H/Fe in terms of both coverage and configuration space can significantly aid in the development of multiscale models as well as provide atomic level insight into the dominant surface structures present under realistic reaction conditions. Here, we construct a lattice gas model for H/Fe(100), where the lateral interactions are determined from first-principles using density functional theory. Using 950 symmetrically unique H/Fe(100) configurations, we generate a cluster expansion with a predictive accuracy in terms of surface energy of 3.8 meV/site over a coverage range from 0 to 3 monolayers. Ten electronic ground state structures are identified from this thorough scan (including the structures at 0 and 3 monolayers), which were subsequently used to generate ab initio phase diagrams under a range of temperatures and pressures. Under reaction conditions typical of Fischer–Tropsch synthesis, ammonia synthesis, and biomass oxygenate hydrodeoxygenation, we find that the 1.0 monolayer structure is dominant. Furthermore, examination of the total H–H lateral interactions for the H/Fe(100) electronic ground state structures shows that H/Fe(100) can be accurately modeled via a mean-field ideal lattice gas model for coverages less than 1.0 monolayers. Altogether, this work enables the incorporation of H–H lateral interactions on Fe(100) into multiscale models, via either mean-field or site-dependent techniques, and provides atomic insight into the catalytically relevant H/Fe(100) structures for a range of heterogeneous reactions.

09 BIOMASS FUELS↗

Capturing Surface Coverage Effects in Heterogeneous Catalysis

Adsorbate–adsorbate lateral interactions at relevant surface coverages have a significant effect on chemical kinetics, thereby influencing the activity of a heterogeneous catalyst. Coverage-dependent kinetic and thermodynamic parameters therefore must be included in studies of such complex systems to properly predict the turnover frequencies and kinetic trends. Thus, it becomes extremely important to accurately capture the strength of lateral interactions between neighboring species under realistic reaction conditions. In this Perspective, we discuss the various existing computational and experimental methods for determining adspecies coverage and configurational effects. The choice of the tools and methods employed in such studies depends on factors such as time, length scales, computational cost, the presence of solvents, and reaction conditions. The applications of each method and the respective challenges are also discussed here. As a result, we discuss the recent developments and future of the state-of-the-art for inclusion of surface coverage and configuration into a holistic picture for accurate predictions of catalytic behavior.

09 BIOMASS FUELS↗

Report on the Status of the SFWST Campaign International Activities in Disposal Research at SNL

The following summaries are provided as fulfillment of milestone M4SF-17SN080305022 and represent international coordination activities in disposal research funded by the US DOE Spent Fuel and Waste Storage and Technologies (SFWST) Campaign during Fiscal Year 2017: SFWST funded bi-lateral interactions with Taiwan, OECD-NEA Repository Metadata (RepMet) project, SFWST funded bi-lateral interactions with the Republic of Korea.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quantitative Effects of Disorder on Chemically Modified Amorphous Carbon Electrodes

Real materials are disordered. This disorder influences the properties of these materials and the chemical processes that occur at their interfaces. Gaining a molecular-level understanding of the underlying physical manifestations caused by disordered materials is crucial to unraveling and ultimately controlling the efficiency and performance of these materials in a range of energy-related devices. This understanding necessitates measurement techniques through which disorder can be detected, quantified, and monitored. However, such quantitative measurements are notoriously difficult, as effects often average out in ensemble measurements. In this work, we describe how a combination of electrochemical and spatially resolved surface spectroscopy measurements illuminate a molecular-level picture of disorder in materials. Using amorphous carbon as an intrinsically disordered material, we covalently attached a monolayer of ferrocene. Interfacial electron transfer across the amorphous carbon–ferrocene interface is highly sensitive to disruptions of order. Here, by systematically varying linker properties and surface loadings, the influence of lateral interactions between nonuniformly distributed ferrocene headgroups on ensemble electrochemical measurements is demonstrated. Electrochemical and imaging data collectively indicate that conformational flexibility of the ferrocene moieties provides a mechanism to elude repulsive and unbalanced lateral interactions, while rigid linkages provide direct information about the underlying disorder of the material. This study is the first of its kind to quantify and visualize molecular disorder and heterogeneity with an experimental model accessed through ensemble measurements.

36 MATERIALS SCIENCE↗

Coverage-Dependent Adsorption of Phenol on Pt(111) from First Principles

In this work, we quantify the coverage dependence of the adsorbate-adsorbate and metal-adsorbate interactions for phenol on Pt(111) using density functional theory. For the four most favorable adsorption sites, we find that the adsorption energy of phenol decreases linearly as a function of coverage. As such, the repulsive phenol-phenol lateral interactions are strongest near saturation when the inter-molecular distances are less than ~4.5 Å, manifesting in a decrease of the C-O dihedral bond angle of ~4°. The linear dependence of the adsorption energy on coverage allows for the construction of a mean-field model. We validate our mean-field model for phenol adsorption by comparing the theoretically predicted differential heat of adsorption to the experimental data, where a maximum deviation of 0.11 eV is found. This compares well to earlier reported results that used a similar mean-field approach for the adsorption of benzene on Pt(111), suggesting that lateral interactions between functionalized aromatics can be parameterized using a simpler approximation. Overall, this work demonstrates a simple method toward the improvement in current state-of-the-art hydrodeoxygenation modeling at surfaces by allowing for the incorporation of coverage effects in systems with large, aromatic compounds without resorting to overly complex models.

09 BIOMASS FUELS↗

The LINKS motif zippers trans-acyltransferase polyketide synthase assembly lines into a biosynthetic megacomplex

Polyketides such as the clinically-valuable antibacterial agent mupirocin are constructed by architecturally-sophisticated assembly lines known as trans-acyltransferase polyketide synthases. Organelle-sized megacomplexes composed of several copies of trans-acyltransferase polyketide synthase assembly lines have been observed by others through transmission electron microscopy to be located at the Bacillus subtilis plasma membrane, where the synthesis and export of the antibacterial polyketide bacillaene takes place. In this work we analyze ten crystal structures of trans-acyltransferase polyketide synthases ketosynthase domains, seven of which are reported here for the first time, to characterize a motif capable of zippering assembly lines into a megacomplex. While each of the three-helix LINKS (Laterally-INteracting Ketosynthase Sequence) motifs is observed to similarly dock with a spatially-reversed copy of itself through hydrophobic and ionic interactions, the amino acid sequences of this motif are not conserved. Such a code is appropriate for mediating homotypic contacts between assembly lines to ensure the ordered self-assembly of a noncovalent, yet tightly-knit, enzymatic network. Furthermore, LINKS-mediated lateral interactions would also have the effect of bolstering the vertical association of the polypeptides that comprise a polyketide synthase assembly line.

59 BASIC BIOLOGICAL SCIENCES↗

Structure sensitivity in oxide catalysis: First-principles kinetic Monte Carlo simulations for CO oxidation at RuO{sub 2}(111)

We present a density-functional theory based kinetic Monte Carlo study of CO oxidation at the (111) facet of RuO{sub 2}. We compare the detailed insight into elementary processes, steady-state surface coverages, and catalytic activity to equivalent published simulation data for the frequently studied RuO{sub 2}(110) facet. Qualitative differences are identified in virtually every aspect ranging from binding energetics over lateral interactions to the interplay of elementary processes at the different active sites. Nevertheless, particularly at technologically relevant elevated temperatures, near-ambient pressures and near-stoichiometric feeds both facets exhibit almost identical catalytic activity. These findings challenge the traditional definition of structure sensitivity based on macroscopically observable turnover frequencies and prompt scrutiny of the applicability of structure sensitivity classifications developed for metals to oxide catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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 ∼ 10{sup 50} erg for an ejecta mass ∼0.5 M {sub ☉}) 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 ∼5 × 10{sup 14} cm with M ∼ 0.1 M {sub ☉}, ejected by the precursor outburst ∼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.

79 ASTRONOMY AND ASTROPHYSICS↗