Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic kinetics”

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

Comparing ab initio and quantum-kinetic approaches to electron transport in warm dense matter

Accurate knowledge of the electronic transport properties of warm dense matter is one of the main concerns of research in high-energy-density physics. Three modern approaches with vastly different levels of fidelity are reviewed and compared: the Kubo–Greenwood (KG) approach based on density-functional-theory molecular dynamics simulations (QMD), quantum kinetic theory based on average-atom models, and time-dependent density functional theory. Throughout, emphasis is placed on the connection between static properties of the electrons (e.g., density of states) and transport properties. Overall, it is found that whenever the conduction electrons can be modeled as being nearly free, fair to excellent agreement is found between QMD and kinetic theory approaches. Such a circumstance is required for modeling warm dense matter as a plasma of ions and free electrons, which is assumed in most kinetic theory approaches. The sensitivity of transport properties to the electronic structure is further highlighted by comparing different exchange–correlation approximations in QMD and KG calculations. It is found that the inclusion of exact exchange via thermal hybrid functionals can make a pronounced impact on electrical and thermal conduction in warm dense matter. We also investigate dynamic screening physics via kinetic theory and time-dependent density functional theory calculations of the mean free path of an electron in a hot dense plasma. In sum, we identify three axes along which to make progress in predicting electron transport in warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cascade electrocatalysis via AgCu single-atom alloy and Ag nanoparticles in CO 2 electroreduction toward multicarbon products

Electrocatalytic CO 2 reduction into value-added multicarbon products offers a means to close the anthropogenic carbon cycle using renewable electricity. However, the unsatisfactory catalytic selectivity for multicarbon products severely hinders the practical application of this technology. In this paper, we report a cascade AgCu single-atom and nanoparticle electrocatalyst, in which Ag nanoparticles produce CO and AgCu single-atom alloys promote C-C coupling kinetics. As a result, a Faradaic efficiency (FE) of 94 ± 4% toward multicarbon products is achieved with the as-prepared AgCu single-atom and nanoparticle catalyst under ~720 mA cm -2 working current density at -0.65 V in a flow cell with alkaline electrolyte. Density functional theory calculations further demonstrate that the high multicarbon product selectivity results from cooperation between AgCu single-atom alloys and Ag nanoparticles, wherein the Ag single-atom doping of Cu nanoparticles increases the adsorption energy of *CO on Cu sites due to the asymmetric bonding of the Cu atom to the adjacent Ag atom with a compressive strain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamics and kinetics of interconversion between platinum nanoparticles and cations in zeolites

In metal-containing zeolites, sintering and redispersion processes exercise control over the identity of metal site structures, but the thermodynamic and kinetic factors that influence these molecular-level processes are not completely understood. Furthermore, we assess the ability of first principles informed free energy models (for supported and unsupported nanoparticles) and kinetic models integrating Ostwald ripening with atom trapping to describe the interconversion between Pt cations and nanoparticles encapsulated in chabazite (CHA) zeolites. Density functional theory-derived thermodynamic phase diagrams show that the interconversion between cations, favored in oxidizing environments, and particles, favored in reducing environments, is fully reversible within a wide range of their respective conditions (temperatures and pressures) for CHA and several other zeolite topologies. Kinetic Monte Carlo simulations of Pt redispersion are consistent with experimentally observed redispersion kinetics of encapsulated Pt nanoparticles in CHA zeolites, and model results suggest the zeolite host imparts additional stability for Pt nanoparticles. We envision our thermodynamic and kinetic models for Pt-CHA are also capable of describing nanoparticle and cation interconversions for other zeolite frameworks under similar conditions.

36 MATERIALS SCIENCE↗

Isotopic Fractionation and Kinetic Isotope Effects of a Purified Bacterial Nitric Oxide Reductase (NOR)

Nitrous oxide (N 2 O) is a serious concern due to its role in global warming and ozone destruction. Agricultural practices account for ~80% of all anthropogenic N 2 O produced in the US, due in large part to the stimulation of microbial denitrification. Stable isotopes are uniquely suited to examine both microbial N 2 O sources and the mechanism of N 2 O biosynthesis through the use of 2 Site Preference (δ 15 N SP ; the difference in δ 15 N between the central and outer N atoms in N 2 O) and kinetic isotope effects (KIEs), respectively. Using trace gas isotope ratio mass spectrometry (TG-IRMS), we determined the δ 15 N, δ 15 N α , δ 15 N β , and δ 18 O of N 2 O produced by a purified cytochrome c nitric oxide reductase (cNOR) from Paracoccus denitrificans. We also calculated δ 15 N SP , the KIEs, and associated isotopic enrichment factors (ε) for N bulk , N α , and N β . A normal isotope effect was observed for bulk 15 N, with a KIE value of 1.0086 ± 0.0009 (ε = -8.6 ± 0.9‰). The isotope effects for both 15 N α and 15 N β were also normal, with position-specific KIEs of 1.0072 ± 0.0010 (ε = -7.2 ± 1.0‰) and 1.0100 ± 0.0010 (ε = -9.9 ± 1.0‰), respectively, and δ 15 N SP values ranged from 0.5 to 8.7‰ with no significant trend as the reaction proceeded. Values of δ 18 O increased with N 2 O production (slope of δ 18 O against [-flnf/(1-f)] = -19.9 ± 1.9‰). We present implications for the mechanism of N 2 O production from cNOR based on our data.

Rivett, Elise D. [Michigan State Univ., East Lansi↗

Fe-single-atom catalyst nanocages linked by bacterial cellulose-derived carbon nanofiber aerogel for Li-S batteries

Li-S battery (LSB) is promising for achieving high capacity. Still, its development is hindered by the complex redox process with sluggish kinetics and particularly the resulting lithium polysulfides (LiPS) shuttle effects. Single-atom catalysts (SACs), with their maximized atom utilization, could effectively chemisorb soluble LiPSs and expedite the sulfide conversion reaction kinetics. Here we report incorporating Fe single metal atom catalyst (Fe-SAC) in the sulfur cathode design and its electrocatalytic effects. Fe-doped ZIF-8 nanocages were introduced into a cheap biomass bacteria cellulose. A pyrolysis process converted them into an aerogel structure with Fe-SAC-functionalized N-doped carbon nanocages linked by a carbon nanofiber network (FeSA-NC@CBC), which was applied as a scaffold to fabricate freestanding and binder-free sulfur cathodes. He we conducted electrochemical measurements to reveal Fe-SAC functions including lowering energy barriers for S 8 reduction to liquid-phase LiPSs and further to solid-phase Li 2 S 2 /Li 2 S and accelerating Li 2 S 2 /Li 2 S nucleation and deposition, as corroborated by our theoretical calculation results. Benefiting from the synergistic effects of highly active Fe-SAC and three-dimensional conductive network, the sulfide reaction kinetics is improved, which can diminish LiPS shuttle effects and therefore improve LBS rate performance and cycling stability. Accordingly, the fabricated FeSA-NC@CBC composite cathode delivers an excellent rate capability at 2C with a reversible capacity of 840 mAh/g and a long-term cyclic stability of 800 mAh/g at 1C after 500 cycles.

25 ENERGY STORAGE↗

Mechanistic insights into nitrogen activation on atomic Ru clusters in self-pillared pentasil using operando atomistic models and experimental kinetics

Alternative catalysts to the industrial Haber Bosch process have been of significant interest in the field of heterogeneous catalysis, yet realizing ammonia synthesis under mild conditions (e.g., 300 °C and 10 bar) is challenging due to the low per-pass conversion. One strategy is to promote the associative ammonia synthesis mechanism which eschews direct N-N bond cleavage. Here, in this work, we use self-pillared pentasil, a self-pillared hierarchical zeolite built by thin MFI zeolite nanosheets, as a support for subnanometric Ru clusters to synthesize ammonia. We show that Ru remains well-dispersed during reaction and further demonstrate that ammonia synthesis rates are higher than Cs-Ru/MgO. Reaction kinetics show a positive order in H 2 providing evidence for the associative mechanism, which then becomes negative in H 2 if Ru is allowed to aggregate into nanoparticles. Operando Density Functional Theory models for Ru speciation in SPP, free energy diagrams, and microkinetic modeling were then applied to develop a reaction mechanism that involves sequential hydrogenation of N 2 from metallic Ru clusters. For this hydrogenation to occur, there are site requirements for N 2 to adopt a bridge-bound configuration that facilitates sequential hydrogenation on single sites and metal clusters. These site requirements in turn inform the design of improved zeolite-supported ammonia synthesis catalysts.

36 MATERIALS SCIENCE↗

Strain release by 3D atomic misfit in fivefold twinned icosahedral nanoparticles with amorphization and dislocations

Multiple twinning to form fivefold twinned nanoparticles in crystal growth is common and has attracted broad attention ranging from crystallography research to physical chemistry and materials science. Lattice-misfit strain and defects in multiple twinned nanoparticles (MTP) are key to understand and tailor their electronic properties. However, the structural defects and related strain distributions in MTPs are poorly understood in three dimensions (3D). Here, we show the 3D atomic misfit and strain relief mechanism in fivefold twinned icosahedral nanoparticles with amorphization and dislocations by using atomic resolution electron tomography. We discover a two-sided heterogeneity in variety of structural characteristics. A nearly ideal crystallographic fivefold face is always found opposite to a less ordered face, forming Janus-like icosahedral nanoparticles with two distinct hemispheres. The disordered amorphous domains release a large amount of strain. Molecular dynamics simulations further reveal the Janus-like icosahedral nanoparticles are prevalent in the MTPs formed in liquid-solid phase transition. This work provides insights on the atomistic models for the modelling of formation mechanisms of fivefold twinned structures and computational simulations of lattice distortions and defects. We anticipate it will inspire future studies on fundamental problems such as twin boundary migration and kinetics of structures in 3D at atomic level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Data from Reactive Species and Reaction Pathways for the Oxidative Cleavage of 4-Octene and Oleic Acid with H2O2 over Tungsten Oxide Catalysts

Oxidative cleavage of carbon–carbon double bonds (C═C) in alkenes and fatty acids produces aldehydes and acids valued as chemical intermediates. Solid tungsten oxide catalysts are low cost, nontoxic, and selective for the oxidative cleavage of C═C bonds with hydrogen peroxide (H2O2) and are, therefore, a promising option for continuous processes. Despite the relevance of these materials, the elementary steps involved and their sensitivity to the form of W sites present on surfaces have not been described. Here, we combine in situ spectroscopy and rate measurements to identify significant steps in the reaction and the reactive species present on the catalysts and examine differences between the kinetics of this reaction on isolated W atoms grafted to alumina and on those exposed on crystalline WO3 nanoparticles. Raman spectroscopy shows that W–peroxo complexes (W–(η2-O2)) formed from H2O2 react with alkenes in a kinetically relevant step to produce epoxides, which undergo hydrolysis at protic surface sites. Subsequently, the CH3CN solvent deprotonates diols to form alpha-hydroxy ketones that react to form aldehydes and water following nucleophilic attack of H2O2. Turnover rates for oxidative cleavage, determined by in situ site titrations, on WOx–Al2O3 are 75% greater than those on WO3 at standard conditions. These differences reflect the activation enthalpies (ΔH‡) for the oxidative cleavage of 4-octene that are much lower than those for the isolated WOx sites (36 ± 3 and 60 ± 6 kJ·mol–1 for WOx–Al2O3 and WO3, respectively) and correlate strongly with the difference between the enthalpies of adsorption for epoxyoctane (ΔHads,epox), which resembles the transition state for epoxidation. The WOx–Al2O3 catalysts mediate oxidative cleavage of oleic acid with H2O2 following a mechanism comparable to that for the oxidative cleavage of 4-octene. The WO3 materials, however, form only the epoxide and do not cleave the C–C bond or produce aldehydes and acids. These differences reflect the distinct site requirements for these reaction pathways and indicate that acid sites required for diol formation are strongly inhibited by oleic acids and epoxides on WO3 whereas the Al2O3 support provides sites competent for this reaction and increase the yield of the oxidative cleavage products.

Catalysis↗

Deep reinforcement learning for predicting kinetic pathways to surface reconstruction in a ternary alloy

Abstract The majority of computational catalyst design focuses on the screening of material components and alloy composition to optimize selectivity and activity for a given reaction. However, predicting the metastability of the alloy catalyst surface at realistic operating conditions requires an extensive sampling of possible surface reconstructions and their associated kinetic pathways. We present CatGym, a deep reinforcement learning (DRL) environment for predicting the thermal surface reconstruction pathways and their associated kinetic barriers in crystalline solids under reaction conditions. The DRL agent iteratively changes the positions of atoms in the near-surface region to generate kinetic pathways to accessible local minima involving changes in the surface compositions. We showcase our agent by predicting the surface reconstruction pathways of a ternary Ni 3 Pd 3 Au 2 (111) alloy catalyst. Our results show that the DRL agent can not only explore more diverse surface compositions than the conventional minima hopping method, but also generate the kinetic surface reconstruction pathways. We further demonstrate that the kinetic pathway to a global minimum energy surface composition and its associated transition state predicted by our agent is in good agreement with the minimum energy path predicted by nudged elastic band calculations.

36 MATERIALS SCIENCE↗

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↗

Theoretical and kinetic modeling study of H 2 S pyrolysis

Hydrogen sulfide pyrolysis was investigated theoretically and through chemical kinetic modeling. Reactions on the SHH potential energy surface, primarily S + H 2 (+Ar) ⇌ H 2 S (+ Ar) (R1) and S + H 2 ⇌ SH + H (R6b) were characterized by ab initio calculations. Results for k 1 were in good agreement with experiment, but deviated strongly below 2000 K from values previously used in modeling. Collider efficiencies for H 2 S, S 2 , and N 2 compared to Ar were calculated for R1. Hydrogen sulfide decomposition experiments reported in literature were re-examined in terms of an updated detailed chemical kinetic model. Concentration profiles for the atomic S at high temperature in shock tubes supported the present value of k 1 and served to constrain the rate constants for reaction of S with SH and H 2 S. To explain results from batch and flow reactors, conducted at high H 2 S concentrations in the 900–1400 K range, a very fast rate constant was required for HSS + H ⇌ SH + SH. Under dilute conditions, the gas-phase chemistry was too slow to compete and the decomposition of H 2 S was controlled by loss on the reactor surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reactive Species and Reaction Pathways for the Oxidative Cleavage of 4-Octene and Oleic Acid with H 2 O 2 over Tungsten Oxide Catalysts

The oxidative cleavage of carbon-carbon double bonds (C=C) in alkenes and fatty acids produces aldehydes and acids valued as chemical intermediates. Solid tungsten oxide catalysts are low cost, non-toxic, and selective for the oxidative cleavage of C=C bonds with hydrogen peroxide (H 2 O 2 ) and are, therefore, a promising option for continuous processes. Despite the relevance of these materials, the elementary steps involved and their sensitivity to the form of W sites present on surfaces have not been described. In this work, we combine in situ spectroscopy and rate measurements to identify significant steps in the reaction, the reactive species present on catalysts, and examine differences between the kinetics of this reaction on isolated W atoms grafted to alumina and on those exposed on crystalline WO 3 nanoparticles. Raman spectroscopy shows that W-peroxo complexes (W-(η 2 -O 2 ) formed from H 2 O 2 ) react in a kinetically relevant step to produce epoxides, which undergo hydrolysis at protic surface sites. Subsequently, the CH 3 CN solvent deprotonates diols to form alpha hydroxy ketones that react to form aldehydes and water following nucleophilic attack of H 2 O 2 . Turnover rates for oxidative cleavage, determined by in situ site titrations, on WO x -Al 2 O 3 are 75% greater than those on WO 3 at standard conditions. These differences reflect activation enthalpies (ΔH ‡ ) for oxidative cleavage of 4-octene that are much lower for the isolated WO x sites (36±3 and 60±6 kJ∙mol -1 for WO x -Al 2 O 3 and WO 3 , respectively) and correlate strongly with the difference between enthalpies of adsorption for epoxyoctane (ΔH ads,epox ), which resembles the transition state for epoxidation. The WO x -Al 2 O 3 catalysts mediate oxidative cleavage of oleic acid with H 2 O 2 following a mechanism comparable to that for 4-octene. The WO 3 materials, however, form only the epoxide and do not cleave the C-C bond or produce aldehydes and acids. These differences reflect the distinct site requirements for these reaction pathways and indicate that acid sites required for diol formation are strongly inhibited by oleic acids and epoxides on WO 3 whereas the Al 2 O 3 support provides sites competent for this reaction and increase the yield of the oxidative cleavage products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomically sharp jagged edges of chemical vapor deposition-grown WS 2 for electrocatalysis

Chemical vapor deposition (CVD)-grown 2D transition metal dichalcogenides can adopt faceted edges. To investigate how sharp these sites can be, we utilize aberration-corrected annular dark-field scanning transmission electron microscopy (ADF-STEM) to resolve the atomic structure of two-dimensional (2D) WS 2 domains that show jagged edges. Nanoscale triangular edge structures with S zigzag terminations are observed. Both the peak and valley regions exhibit near-atomic sharpness. The peaks are as sharp as two atoms in width. Highly ordered valley sites display a minimum width of three atoms, and prospective single-atom valleys appear as two-atom-wide sites in the ADF-STEM contrast. Regarding the kinetics, density-functional theory (DFT) calculations indicate that the 2-W-atom peak would not evolve into a single-W-atom peak even though the latter configuration is also stable with a high enough binding energy under the growth conditions. These results help deepen our understanding of the possible structuring at the nanoscale and the atomic-scale limits of peaks and valleys formed via intersection of two zigzag edges. The enriched edge sites lead to higher catalytic activities for the hydrogen evolution reaction (HER).

36 MATERIALS SCIENCE↗

Neural network kinetics for exploring diffusion multiplicity and chemical ordering in compositionally complex materials

Diffusion involving atom transport from one location to another governs many important processes and behaviors such as precipitation and phase nucleation. The inherent chemical complexity in compositionally complex materials poses challenges for modeling atomic diffusion and the resulting formation of chemically ordered structures. Here, we introduce a neural network kinetics (NNK) scheme that predicts and simulates diffusion-induced chemical and structural evolution in complex concentrated chemical environments. The framework is grounded on efficient on-lattice structure and chemistry representation combined with artificial neural networks, enabling precise prediction of all path-dependent migration barriers and individual atom jumps. To demonstrate the method, we study the temperature-dependent local chemical ordering in a refractory NbMoTa alloy and reveal a critical temperature at which the B2 order reaches a maximum. The atomic jump randomness map exhibits the highest diffusion heterogeneity (multiplicity) in the vicinity of this characteristic temperature, which is closely related to chemical ordering and B2 structure formation. The scalable NNK framework provides a promising new avenue to exploring diffusion-related properties in the vast compositional space within which extraordinary properties are hidden.

36 MATERIALS SCIENCE↗

Adsorption and reaction of methanol on Fe 3 O 4 (001)

The interaction of methanol with iron oxide surfaces is of interest due to its potential in hydrogen storage, and from a fundamental perspective as a chemical probe of the reactivity. We present here a study examining the adsorption and reaction of methanol on magnetite Fe 3 O 4 (001) at cryogenic temperatures using a combination of temperature programmed desorption, x-ray photoelectron spectroscopy, and scanning tunneling microscopy. We report the methanol desorption profile from Fe 3 O 4 (001) is complex, exhibiting peaks at 145, 175, 238, and 273 K, corresponding to the desorption of intact methanol, as well as peaks at 350 and 500 K due to the reaction of methoxy intermediates. The saturation of a monolayer of methanol corresponds to ~5 molecules/unit cell (u.c.), which is slightly higher than the number of surface octahedral iron atoms of 4/u.c. We probe the kinetics and thermodynamics of the desorption of molecular methanol using inversion analysis. The deconvolution of the complex desorption profile into individual peaks allows for calculations of both the desorption energy and the prefactor of each feature. The initial 0.7 methanol/u.c. reacts to form a methoxy and hydroxy intermediates by 180 K, which remains on the surface above room temperature after intact methanol has desorbed. The methoxy species react via one of two channels, a recombination reaction with surface hydroxyls to form additional methanol at ~350 K, and a disproportionation reaction to form methanol and formaldehyde at ~500 K. Only 20% of the methoxy species undergo the disproportionation reaction, with most of them reacting via the 350 K pathway.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dy adsorption on and intercalation under graphene on 6 H -SiC(0001) surface from first-principles calculations

Previous experimental observations motivate clarification of configuration stabilities and kinetic processes for intercalation of guest atoms into a layered van der Waals material such as a graphene-SiC system. From our first-principles density functional theory (DFT) calculations, we analyze Dy adsorption and intercalation for graphene on a 6H-SiC(0001) surface, where the system includes two single-atom-thick graphene layers: the top-layer graphene (TLG) and the underling buffer-layer graphene (BLG) above the terminal Si layer. Our chemical potential analysis shows that intercalation of a single Dy atom into the gallery between TLG and BLG is more favorable than adsorption on TLG but that intercalation into the gallery underneath BLG is highly unfavorable. We obtain diffusion barriers of ~0.45 and 0.54 eV for a Dy atom diffusing on and under TLG, respectively. We find that the direct penetration of a Dy atom from the graphene top into the gallery under TLG is almost inhibited below a temperature of ~1400 K due to a large global barrier of at least ~3.5 eV. Instead, we find that a single Dy atom on TLG can easily intercalate by crossing a TLG step (e.g., a zigzag step presaturated by a Dy chain or a reconstructed zigzag step zz57). Additionally, we also perform DFT calculations for different Dy coverages to demonstrate how the favorability of Dy intercalation, as well as the corresponding interlayer spacings, depend on the coverage. Consequently, we can provide general insight and guidance for extensively studied systems involving intercalation of foreign atoms into graphene on a SiC substrate.

36 MATERIALS SCIENCE↗

Origin of Selective Production of Hydrogen Peroxide by Electrochemical Oxygen Reduction

Oxygen reduction reaction (ORR) is one of the most important electrochemical reactions. Starting from a common reaction intermediate *–O–OH, the ORR splits into two pathways, either producing hydrogen peroxide (H 2 O 2 ) by breaking the *–O bond or leading to water formation by breaking the O–OH bond. However, it is puzzling why many catalysts, despite the strong thermodynamic preference for the O–OH breaking, exhibit high selectivity for hydrogen peroxide. Moreover, the selectivity is dependent on the potential and pH, which remain not understood. Here, for this work, we develop an advanced first-principles model for effective calculation of the electrochemical reaction kinetics at the solid–water interface, which were not accessible by conventional models. Using this model to study representative catalysts for H 2 O 2 production, we find that breaking the O–OH bond can have a higher energy barrier than breaking *–O, due to the rigidity of the O–OH bond. Importantly, we reveal that the selectivity dependence on potential and pH is rooted into the proton affinity to the former/later O in *–O–OH. For single cobalt atom catalyst, decreasing potential promotes proton adsorption to the former O, thereby increasing the H 2 O 2 selectivity. In contrast, for the carbon catalyst, the proton prefers the latter O, resulting in a lower H 2 O 2 selectivity in acid condition. These findings explain the experiments and highlight the kinetic origins of the selectivity. Our work improves the understanding of ORR by uncovering the proton affinity as a new factor and provides a new model to effectively simulate the atomic-level kinetics of heterogeneous electrochemistry.

08 HYDROGEN↗

Discovering Ni/Cu Single-Atom Alloy as a Highly Active and Selective Catalyst for Direct Methane Conversion to Ethylene: A First-Principles Kinetic Study

Direct methane conversion to liquid fuels or value-added chemicals is a promising technology to utilize natural resources without resorting to further petroleum extraction. However, discovering efficient catalysts for this reaction is challenging due to either coke formation or unfavorable C-H bond activation. Herein, we design single-atom alloy (SAA) catalysts to simultaneously eliminate the above two bottlenecks based on mechanism-guided strategies: (1) the active single atom enables favorable C-H bond breaking and (2) the less reactive host metal facilitates C-C coupling and thus avoids strong binding of carbonaceous species. Employing electronic structure theory calculations, we screened the stability of multiple SAAs with 3d-5d transition metals atomically dispersed on a copper surface in terms of avoiding dopant aggregation and segregation. We then evaluated reactivities of the stable SAAs as catalysts for direct methane conversion to C2 products, including methane dehydrogenation and C-C coupling mechanisms. Combining selectivity analysis with kinetic modeling, we predicted that nickel dispersed on copper, i.e., Ni/Cu SAA, is a highly active and selective catalyst that can efficiently transform methane to ethylene. This work designs efficient SAA catalysts for direct methane activation and provides chemical insights into engineering compositions of SAAs to tune their catalytic performances.

Kothakonda, Manish↗