Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic barrier”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Design and Optimization of Structured Multi-Functional Trapping Catalysts for Conversion of Hydrocarbons and NOx from Diesel and Advanced Combustion Engines

Oxides of nitrogen in the form of nitric oxide (NO) and nitrogen dioxide (NO 2 ) commonly referred to as NOx, is one of the two chemical precursors that lead to ground-level ozone, a ubiquitous air pollutant in urban areas. A major source of NOx is generated by equipment and vehicles powered by diesel engines, which have a combustion exhaust that contains NOx in the presence of excess O 2 . Vehicular emission control catalysts are ineffective in eliminating CO, hydrocarbons, and NO x during engine cold-start when exhaust temperatures are below 200°C. The objective of the project was to develop and demonstrate a multi-functional, catalyzed trap that enables vehicles with advanced combustion strategies to meet Tier 3 emissions standards while achieving the 150 °C challenge for sustained co-oxidation of HCs and CO and ≥90% NO trapping and release during warmup. Specifically, the multi-functional Lean HC+NOx (LHCNT) was developed for application in the exhaust aftertreatment of conventional diesel engines and engines having low temperature combustion (LTC) regimes. Activities included the design and synthesis of adsorbents and catalysts, screening and evaluation. Passive NOx absorbers (PNA), hydrocarbon (HC) traps, and oxidation catalysts (OC) were evaluated for use in series or as integrated devices. Predictive tools were developed utilizing the characterization and analysis of these materials, and an emission system was designed and optimized utilizing the catalyst systems. Microkinetic models were developed for the PNA for the simple NO-only feed and complex feed containing CO, H 2 , and model hydrocarbons (ethylene and dodecane). A first-principles, mechanistic-based model of the PNA was developed which utilized molecular-scale estimates (density functional theory) of energy barriers, mechanistic-based kinetics and realistic treatments of the flow and transport processes. Two new oxidation catalysts were developed (PdCu alloy, mixed copper-ceria-cobalt oxide), both of which significantly lessened the detrimental inhibition by CO on hydrocarbon and NO oxidation. A method for lessening the detrimental impact of CO on PNA activity was developed that involves use of an oxidation catalyst upstream of the PNA. The SwRI Ectolab TM burner system was applied to evaluate the baseline PNA material and confirmed performance comparable to the benchflow PNA studies using simulated exhaust. Spatially-resolved mass spectrometry (SpaciMS) was used to measure the transient spatial profiles of reacting species spanning the length of a three-function LHCNT containing PNA, HCT, and OC. The findings from this study provide diesel vehicle and catalyst companies valuable information to develop more cost effective emission control catalysts which helps to expand the use of more fuel efficient diesel power. The fundamental modeling and experimental tools and findings from this project can be applied to catalyst technologies used in the energy and chemical industries. The project led to 14 publications in the peer-reviewed literature with 2 additional currently under review. Finally, the project also led to training of several doctoral students who were placed in research jobs in industry and academia. Specifically, Mugdha Ambast (UH) has joined Cummins, Kevin Gu (UVa) has joined GM, and Abhay Gupta (UH) is to join Caterpillar.

42 ENGINEERING↗
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Recruiting Perovskites to Degrade Toxic Trinitrotoluene

Everybody knows TNT, the most widely used explosive material and a universal measure of the destructiveness of explosions. A long history of use and extensive manufacture of toxic TNT leads to the accumulation of these materials in soil and groundwater, which is a significant concern for environmental safety and sustainability. Reliable and cost-efficient technologies for removing or detoxifying TNT from the environment are lacking. Despite the extreme urgency, this remains an outstanding challenge that often goes unnoticed. We report here that highly controlled energy release from explosive molecules can be accomplished rather easily by preparing TNT–perovskite mixtures with a tailored perovskite surface morphology at ambient conditions. These results offer new insight into understanding the sensitivity of high explosives to detonation initiation and enable many novel applications, such as new concepts in harvesting and converting chemical energy, the design of new, improved energetics with tunable characteristics, the development of powerful fuels and miniaturized detonators, and new ways for eliminating toxins from land and water.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Spontaneous and Ion-Specific Formation of Inverted Bilayers at Air/Aqueous Interface

Developing better separation technologies for rare earth metals, an important aspect of a sustainable materials economy, is challenging due to their chemical similarities. Identifying molecular scale interactions that amplify the subtle differences between the rare earths can be useful in developing new separation technologies. Here, we describe ion-dependent monolayer to inverted bilayer transformation of extractant molecules at the air/aqueous interface. The inverted bilayers form with Lu 3+ ions but not with Nd 3+ . By introducing Lu 3+ ions to preformed monolayers, we extract kinetic parameters corresponding to the monolayer to inverted bilayer conversion. Further, temperature-dependent studies show Arrhenius behavior with an energy barrier of 40 kcal/mol. The kinetics of monolayer to inverted bilayer conversion is also affected by the character of the background anion, although anions are expected to be repelled from the interface. Our results show the outsized importance of ion-specific effects on interfacial structure and kinetics, pointing to their role in chemical separation methods.

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

Quantum Sieving for Isotopic Separations of Gases Using Porous Materials─30 Years of Progress

Here, this paper reviews theoretical and experimental efforts to establish microporous materials that can be used to separate isotopologues of molecular gases such as H 2 and D 2 . Emphasis is placed on use of simplified models to highlight the quantum phenomena that make these separations possible. In equilibrium adsorption, differences in zero point energy in the adsorbed states of molecules favor binding of heaver isotopologues (e.g. D 2 relative to H 2 ). Experimental data showing this effect was reported as early as 1933, but the theoretical work of Beenakkers et al. in 1995 [Chem. Phys. Lett. 232 (1995) 379-382] spurred modern efforts to develop materials with strong so-called quantum sieving. In porous materials where the transition state for molecular diffusion is more strongly confined than the energy minima associated with equilibrium adsorption, differences in zero point energies between these two sites can lead to isotopic differences in molecular diffusivities. This effect favors diffusion of heavier species (e.g. D 2 ) relative to lighter species (e.g. H 2 ). This so-called kinetic quantum sieving has been observed experimentally in porous carbons and in porous organic cage materials. We show that quantum tunneling, which favors hopping of lighter species across energy barriers, diminishes the strength of kinetic quantum sieving but that it appears to make only a small contribution to the net molecular diffusivities in many porous materials of practical interest.

Sholl, David S. [Oak Ridge National Laboratory (OR↗
DOE OSTI2023

Unveiling the catalytic mechanism of GTP hydrolysis in microtubules

Microtubules (MTs) are large cytoskeletal polymers, composed of αβ-tubulin heterodimers, capable of stochastically converting from polymerizing to depolymerizing states and vice versa. Depolymerization is coupled with hydrolysis of guanosine triphosphate (GTP) within β-tubulin. Hydrolysis is favored in the MT lattice compared to a free heterodimer with an experimentally observed rate increase of 500- to 700-fold, corresponding to an energetic barrier lowering of 3.8 to 4.0 kcal/mol. Mutagenesis studies have implicated α-tubulin residues, α:E254 and α:D251, as catalytic residues completing the β-tubulin active site of the lower heterodimer in the MT lattice. The mechanism for GTP hydrolysis in the free heterodimer, however, is not understood. Additionally, there has been debate concerning whether the GTP-state lattice is expanded or compacted relative to the GDP state and whether a “compacted” GDP-state lattice is required for hydrolysis. In this work, extensive quantum mechanics/molecular mechanics simulations with transition-tempered metadynamics free-energy sampling of compacted and expanded interdimer complexes, as well as a free heterodimer, have been carried out to provide clear insight into the GTP hydrolysis mechanism. α:E254 was found to be the catalytic residue in a compacted lattice, while in the expanded lattice, disruption of a key salt bridge interaction renders α:E254 less effective. The simulations reveal a barrier decrease of 3.8 ± 0.5 kcal/mol for the compacted lattice compared to a free heterodimer, in good agreement with experimental kinetic measurements. Additionally, the expanded lattice barrier was found to be 6.3 ± 0.5 kcal/mol higher than compacted, demonstrating that GTP hydrolysis is variable with lattice state and slower at the MT tip.

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

Energy barriers for Dy and H penetrating graphene on 6 H -SiC(0001) and freestanding bilayer graphene from first-principles calculations

Currently, intercalation of foreign guest atoms into two-dimensional (2D) layered van der Waals materials is an active research area motivated in part by the development of next-generation energy-storage technologies and optoelectronic devices. One such extensively studied 2D material is the graphene-on-SiC system. To realize and control the desired intercalated structures, it is fundamentally important to understand the kinetic process of intercalation. For the intercalation of a guest atom into graphene layers on SiC substrate, a critical kinetic parameter is the energy barrier of a guest atom penetrating the perfect graphene top layer into the gallery under it. However, accurate theoretical calculations for such penetration barriers are unavailable in literature. From our first-principles density functional theory calculations, we obtain the global energy barriers of 3.47 and 1.80 eV for single Dy and H atoms penetrating the graphene top layer on a graphene buffer layer supported by a Si-terminated 6H-SiC(0001) substrate, respectively. Furthermore, for comparison as well as for examining the lateral strain effects, we also obtain the global barriers of 5.05 and 1.50 eV for single Dy and H atoms penetrating freestanding bilayer graphene with a tensile strain of about 8.8% to match our model for supported graphene, as well as the global barriers of 7.21 and 4.18 eV for penetrating unstrained freestanding bilayer graphene, respectively. From corresponding minimum energy paths with multiple energy minima and saddle points, we can also obtain various local energy barriers and the global backward barrier from the graphene gallery back to the top surface.

2D materials↗
DOE OSTI2026

Closed‐Loop Recycling of Poly(pentylene Adipate‐ co ‐Terephthalate) via Amine‐Catalyzed Methanolysis

Advancing circularity through effective chemical recycling strategies is essential for enabling the sustainable use of emerging biodegradable polyesters in film applications. In this context, poly(pentylene adipate-co-terephthalate) with a 40/60 adipate-to-terephthalate ratio (PPeAT60), a significantly stiffer and partially biobased copolyester recently reported by our team, offers improved mechanical performance and processability compared to commercial poly(butylene adipate-co-terephthalate), while providing a platform to evaluate advanced recycling approaches for film applications. Here, this study describes depolymerization of PPeAT60 through amine-catalyzed methanolysis, followed by the repolymerization of the recovered monomers. The resulting repolymerized material exhibited a molar mass distribution, chemical structure, and thermal stability comparable to the virgin polymer, confirming the technical feasibility of the methanolysis-based recycling route. Both polymers displayed nearly identical melting, glass transition, and decomposition temperatures, crystallinity, and modulus, demonstrating full structural and mechanical reproducibility after recycling. Variations in crystallization kinetics, elongation at break, and barrier properties likely arise from trace impurities (<1%) introduced during processing acting as nucleating and stress-inducing agents. Overall, these results demonstrate that PPeAT60 retains its key properties after chemical recycling, establishing methanolysis as a promising platform for recycling of this partially biobased polymer that may be able to replace both conventional polyolefins and PBAT in flexible packaging applications.

09 BIOMASS FUELS↗
DOE OSTI2025

Water‐mediated Decomposition Pathways of ε‐Hexanitrohexaazaisowurtzitane

The dominant initial decomposition mechanism of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (HNIW) is understood from the literature to be the N─NO 2 bond scission unimolecular mechanism. Here we have found a water-mediated HONO release initial decomposition mechanism energetically competitive with the N─NO 2 bond scission mechanism modeled using ab initio nudged elastic band simulations. The activation energy of the water-mediated HONO mechanism was calculated to be 138.9 kJ/mol with the Perdew, Burke, and Ernzerhof (PBE) functional (177.8 with PBE0), while the double NO 2 release mechanism is 171.9 kJ/mol (203.3 with PBE0). Branching secondary, tertiary, and quaternary decomposition steps were also discovered including oxidation of HONO or H 2 O causing ring-opening that leads to C─N bond breaking and N 2 O or NO release. The reaction rate of HONO oxidation is much faster than HONO release, making the HONO concentration low. The Helmholtz free energy barrier of the double NO 2 release mechanism is lower than the water-mediated HONO release above 438 K (165°C) due to the vibrational contribution to the free energy. This helps explain why the literature reports more NO 2 release at higher temperatures. The present study reveals a physical mechanism for how water can catalyze the decomposition of HNIW at low enough temperatures, and reveals secondary, tertiary, and quaternary reaction mechanisms leading to NO x release. In conclusion, as water is realistically always present, the reported barriers are important for building kinetic models needed to ensure the safety and reliability of HNIW.

Steele, Brad A. [Lawrence Livermore National Labor↗
DOE OSTI2023

Plasma-assisted deflagration to detonation transition in a microchannel with fast-frame imaging and hybrid fs/ps coherent anti-Stokes Raman scattering measurements

Our study examines kinetic enhancement by nanosecond dielectric barrier discharge (ns-DBD) plasma on fuellean dimethyl ether (DME), oxygen (O 2 ), and argon (Ar) premixtures during deflagration to detonation transition (DDT) experiments in a microchannel. Non-equilibrium plasma produces active species and radicals and creates fast and slow heating of a mixture to promote ignition due to electronic and vibrational excitation. Experiments have been conducted to examine the influence of the plasma discharge on the premixture and on the resultant deflagration to detonation transition (DDT) onset time and distance through the use of high speed imaging and one-dimensional, two-beam, femtosecond/picosecond, coherent anti-Stokes Raman scattering (CARS). A highspeed camera is used to trace the time histories of flame front position and velocity and to identify the dynamics and onset of DDT. The results show that plasma discharge can nonlinearly affect the onset time and distance of DDT. It is shown that a small number of plasma discharge pulses prior to ignition result in reduced DDT onset time and distance by 60% and 40%, respectively, when compared to the results without pre-excitation by ns discharges. The results also show that an increase of plasma discharge pulses results in an extended DDT onset time and distance of 224% and 94%, respectively. Time history of the deflagration wave speed of DME and the analysis of ignition timescale under the choking condition of the deflagration front suggest low temperature ignition may play a role for DME near the isobaric choking condition of the burned gas and the DDT. Plasmainduced conversion of the reactive mixture was assessed via the O 2 to CO 2 ratio as measured through fs/ps CARS during the DBD discharges. The present experiments demonstrate the ability of non-equilibrium plasma to alter the chemistry of DME/O 2 /Ar premixtures in order to control DDT for applications in advanced propulsion engines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Harnessing nanoscale densification for controlling gas selectivity in flexible zeolitic imidazolate frameworks

Flexible metal-organic frameworks (MOFs) are promising for gas separation due to their molecular recognition capabilities. However, gate-opening behavior in certain flexible MOFs can reduce gas selectivity, enabling simultaneous adsorption of multiple gases. We introduce a binder-free, pressure-free densification process to enhance gas selectivity in flexible zeolitic imidazolate frameworks (ZIF-7 and ZIF-9) by optimizing nanoparticle packing. Unlike powders, the densified ZIF-7 and ZIF-9 retain rhombohedral crystal structures after solvent removal, confirmed via synchrotron powder X-ray diffraction. This process modifies adsorption behavior by introducing a diffusion barrier, enabling diffusivity-driven kinetic separation and enhanced CO2/hydrocarbon selectivity. At 0.5 bar, ZIF-7 mono shows 2×, 4×, and 9.5× higher CO 2 selectivity over C 2 H 2 , C 2 H 4 , and C 2 H 6 , respectively, versus ZIF-7 pwd . ZIF-9 mono similarly shows 6.5×, 3×, and 39× improvements over ZIF-9 pwd . Both monoliths outperform powders in ternary gas mixtures. This study provides a novel strategy to control gate opening and enhance gas selectivity in flexible MOFs for practical separation applications.

CO2/hydrocarbon separation↗
DOE OSTI2026

Insights into Nonelectroactive C–C Bond Formation on Cu(100) during Electrochemical CO 2 Reduction from Multiconfigurational Wavefunction Theory

Carbon–carbon (C–C) bond formation is necessary for hydrocarbon (and oxygenate) synthesis beyond methane (and formate/formic acid) during electrochemical CO and CO 2 reduction (ECOR and ECO 2 R). Cu has notable ability to form hydrocarbons compared to other pure metals. In particular, the (100) facet of face-centered cubic Cu forms ethylene competitively with H 2 and methane during both ECOR and ECO 2 R. Past simulations based on density functional theory (DFT) with standard exchange-correlation functional approximations predict fast nonelectroactive C–C bond formation channels involving adsorbed (*) CO together with another *CO, formyl (*CHO), or hydroxymethylidyne (*COH), forming OC*–*CO, OC*–CHO*, and OC*–*COH, respectively. Such simulations support the prevailing hypothesis that emergence of C 2 products is kinetically determined at the early stages of the reduction chemistry. Here we show, via simulations with more accurate many-body, i.e., “correlated”, wavefunction theory (enabled by an embedding scheme), that the coupling of *CO with a *CO or a *COH (previously predicted at the same level of theory to kinetically dominate over *CHO as the one-electron reduction product of *CO) is highly activated (kinetically impeded), with free energy barriers >1 eV, in contradiction to previous DFT-based simulations. Intriguingly, we find that the coupling of two adjacent *COHs incurs only a small barrier (<0.3 eV) and is exoergic (< –1 eV); however, given the predicted low surface mobility of *COH, the emergence of HOC*–*COH is also improbable, at least at low *COH coverages. We therefore conclude that it is highly unlikely for *CO to participate in nonelectroactive C–C bond formation on pristine Cu(100), contrary to conventional wisdom, and that the energetically favorable *COH dimerization may occur only after substantial buildup of *COH on the surface.

adsorption↗
DOE OSTI2020

Prediction of Highly Selective Electrocatalytic Nitrogen Reduction at Low Overpotential on a Mo-Doped g-GaN Monolayer

Identifying efficient electrocatalysts with low overpotential and high selectivity for producing ammonia from nitrogen gas is essential for any future electrocatalytic nitrogen reduction reaction (NRR)-based ammonia synthesis. Via density functional theory calculations and the computational hydrogen electrode model, we systematically examine the prospect of using a single-transition-metal (TM)-atom-doped graphene-like GaN (g-GaN) monolayer as an electrocatalyst for artificial nitrogen reduction. Among 15 TMs investigated, the Mo-doped g-GaN (Mo@g-GaN) monolayer is the only electrocatalyst predicted to be feasible for the NRR. The Mo@g-GaN monolayer satisfies all screening criteria considered for activating the inert N≡N triple bond effectively, including stabilization of the adsorbed (*) NRR intermediate *NNH and destabilization of the *NH 2 species. This monolayer also possesses sufficient overall stability. A complete analysis of the likely mechanisms involved in the NRR on this catalyst suggests that the Mo@g-GaN monolayer could exhibit promising NRR catalytic activity. It achieves this via one specific (distal) pathway, which has a very low onset potential of –0.33 V vs the reversible hydrogen electrode (RHE), corresponding to a low overpotential of 0.42 V vs the RHE, defined using the measured equilibrium potential for NRR of 0.09 V vs the RHE. The potential-determining step, conversion of *NH 2 to *NH 3 , also exhibits a surmountable barrier of 0.42 eV, suggesting kinetics will be facile. Finally, the Mo@g-GaN monolayer is predicted to exhibit substantial selectivity (~31%) toward ammonia synthesis over the competing hydrogen evolution reaction. The finding presented in this work may open a potential route for artificial ammonia synthesis using a single-atom catalyst under ambient conditions.

25 ENERGY STORAGE↗
DOE OSTI2024

Formation of the oxyl’s potential energy surface by the spectral kinetics of a vibrational mode

One of the most reactive intermediates for oxidative reactions is the oxyl radical, an electron-deficient oxygen atom. The discovery of a new vibration upon photoexcitation of the oxygen evolution catalysis detected the oxyl radical at the SrTiO3 surface. The vibration was assigned to a motion of the sub-surface oxygen underneath the titanium oxyl (Ti–O●−) created upon hole transfer to (or electron extraction from) a hydroxylated surface site. Evidence for such an interfacial mode is derived from its spectral shape, which exhibited a Fano resonance—a coupling of a sharp normal mode to continuum excitations. Here, this Fano resonance is utilized to derive precise formation kinetics of the oxyl radical and its associated potential energy surface (PES). From the Fano lineshape, the formation kinetics are obtained from the anti-resonance (the kinetics of the coupling factor), the resonance (the kinetics of the coupled continuum excitations), and the frequency integrated spectrum (the kinetics of the normal mode’s cross-section). All three perspectives yield logistic function growth with a half-rise of 2.3 ± 0.3 ps and a time constant of 0.48 ± 0.09 ps. A non-equilibrium transient associated with photoexcitation is separated from the rise of the equilibrated PES. The logistic function characterizes the oxyl coverage at the very initial stages (t ∼ 0) to have an exponential growth rate that quickly decreases toward zero as a limiting coverage is reached. Such time-dependent reaction kinetics identify a dynamic activation barrier associated with the formation of a PES and quantify it for oxyl radical coverage.

Chemistry↗
DOE OSTI2023

Baseline CGMF Calculations for Pu Isotopes

Updates were made to CGMF to run neutron-induced fission calculations on Pu isotopes of A =238, 239, 240, 241, and 242 (the release version of CGMF can calculate 239,241 Pu(n,f) and 238,240,242 Pu(sf) only). Here, we show the baseline calculations from CGMF for these neutron-induced fission reactions, including updates made for total kinetic energy systematics and fission barriers fit to be more consistent with the neutron-induced fission cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Capture efficiency and bias from the defect dynamics near grain boundaries in BCC Fe using mesoscale simulations

The capture efficiency describes the capability of a sink, such as a grain boundary (GB), dislocation, and void, to absorb point defects (PDs). The bias defines the difference in capture efficiency between the absorption of a vacancy and dumbbell at a sink. Complete kinetic information on PDs, including diffusion barriers and diffusion orientations, as well as accurate saddle points, are needed to determine the capture efficiency and bias at a sink accurately, which is computationally demanding. Here, in the present study, the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC) method was used to investigate the defect dynamics of PDs near different types of grain boundaries (GBs) (with both $\langle$100$\rangle$ and $\langle$110$\rangle$ families) accurately in body-centered cubic (BCC) iron (Fe). The capture efficiency, sink strength, and bias factor of different types of GBs were determined in Fe, which, different from traditional rate theory estimation, showed a distinct capture efficiency, sink strength, and bias in different GBs. The results demonstrate a strong positive correlation between the capture efficiency and the GB strain width, instead of the GB misorientation, GB energy, or GB-PD binding energy, which have been investigated previously. This work provides valuable insight into the radiation-induced microstructural evolution of GBs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stability of heterogeneous single-atom catalysts: a scaling law mapping thermodynamics to kinetics

Heterogeneous single-atom catalysts (SACs) hold the promise of combining high catalytic performance with maximum utilization of often precious metals. We extend the current thermodynamic view of SAC stability in terms of the binding energy (E bind ) of single-metal atoms on a support to a kinetic (transport) one by considering the activation barrier for metal atom diffusion. A rapid computational screening approach allows predicting diffusion barriers for metal–support pairs based on E bind of a metal atom to the support and the cohesive energy of the bulk metal (E c ). Metal–support combinations relevant to contemporary catalysis are explored by density functional theory. Assisted by machine-learning methods, we find that the diffusion activation barrier correlates with (E bind ) 2 /E c in the physical descriptor space. This diffusion scaling-law provides a simple model for screening thermodynamics to kinetics of metal adatom on a support.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mechanistic study of direct coupling of CO 2 and C 2 H 4 over atomically dispersed metal at graphene edges

Direct coupling of CO 2 and ethylene (hereinafter DCCE) to acrylic acid is valuable for valorizing CO 2 to manufacture acrylate-derived products. However, previous studies in DCCE have been limited on molecular catalysts with challenges in improving catalytic performance. In this work, we employed density functional theory calculations and ab initio molecular dynamics simulations to investigate the heterogeneous catalysis of DCCE over atomically dispersed metal centers at nitrogen-doped zigzag edge of graphene. Based on competitive adsorption and structural stability, Mo, Cr, V, Ru, and Ni active sites are chosen to explore the reaction kinetics. Here, we find that the activation barriers are determined by the charge redistribution at transition states, which explains the trend of activity for the C-C coupling and the hydrogen transfer, two key steps in DCCE. Furthermore, we show that the intramolecular hydrogen transfer (rate-limiting step) is hindered due to the lack of local coordinate at the active sites. We thus propose to use co-adsorbed water as a “proton-exchanger” following a water-assisted route, and show that the activation barriers are reduced over all metal centers. Particularly, water promotes the hydrogen transfer over metals with strong CO 2 -ethylene co-activation and facile C-C coupling kinetics, which could be considered promising for DCCE. In both mechanisms, the stability of metallactone intermediate can be used to predict the catalytic activity. It is anticipated that the insights from this work can provide guidelines for mimicking well-defined multifunctional active sites in molecular catalysts to design heterogeneous catalysts for such C-C coupling, which advances catalytic utilization of CO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗