Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface energetics”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Unraveling Exciton Trap Dynamics and Nonradiative Loss Pathways in Quantum Dots via Atomistic Simulations

Surface defects in colloidal quantum dots are a major source of nonradiative losses, yet the microscopic mechanisms underlying exciton trapping and recombination remain elusive. Here, we develop a model Hamiltonian based on atomistic electronic calculations to investigate exciton dynamics in CdSe/CdS core/shell QDs containing a single hole trap introduced by an unpassivated sulfur atom. By systematically varying the defect depth and reorganization energy, we uncover how defect-induced excitonic states mediate energy relaxation pathways. Our simulations reveal that a single localized defect can induce a rich spectrum of excitonic states, leading to multiple dynamical regimes, from slow, energetically off-resonant trapping to fast, cascaded relaxation through in-gap defect states. Crucially, we quantify how defect-induced polaron shifts and exciton-phonon couplings govern the balance between efficient radiative emission and rapid nonradiative decay. These insights clarify the microscopic origin of defect-assisted loss channels and suggest pathways for tailoring QD optoelectronic properties via surface and defect engineering.

Defects↗
DOE OSTI2021

Non-adiabatic quantum interference in the ultracold Li + LiNa → Li 2 + Na reaction

Electronically non-adiabatic effects play an important role in many chemical reactions. However, how these effects manifest in cold and ultracold chemistry remains largely unexplored. Here for the first time we present from first principles the non-adiabatic quantum dynamics of the reactive scattering between ultracold alkali-metal LiNa molecules and Li atoms. We show that non-adiabatic dynamics induces quantum interference effects that dramatically alter the ultracold rotationally resolved reaction rate coefficients. The interference effect arises from the conical intersection between the ground and an excited electronic state that is energetically accessible even for ultracold collisions. These unique interference effects might be exploited for quantum control applications such as a quantum molecular switch. The non-adiabatic dynamics are based on full-dimensional ab initio potential energy surfaces for the two electronic states that includes the non-adiabatic couplings and an accurate treatment of the long-range interactions. A statistical analysis of rotational populations of the Li2 product reveals a Poisson distribution implying the underlying classical dynamics are chaotic. The Poisson distribution is robust and amenable to experimental verification and appears to be a universal property of ultracold reactions involving alkali metal dimers.

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

Reshaping of Truncated Pd Nanocubes: Energetic and Kinetic Analysis Integrating Transmission Electron Microscopy with Atomistic-Level and Coarse-Grained Modeling

Stability against reshaping of metallic fcc nanocrystals synthesized with tailored far-from-equilibrium shapes is key to maintaining optimal properties for applications such as catalysis. Yet Arrhenius analysis of experimental reshaping kinetics, and appropriate theory and simulation, is lacking. As such, we use TEM to monitor the reshaping of Pd nanocubes of ~25 nm side length between 410 °C (over ~4.5 h) and 440 °C (over ~0.25 h), extracting a high effective energy barrier of E eff ≈ 4.6 eV. We also provide an analytic determination of the energy variation along the optimal pathway for reshaping that involves transfer of atoms across the nanocube surface from edges or corners to form new layers on side {100} facets. The effective barrier from this analysis is shown to increase strongly with the degree of truncation of edges and corners in the synthesized nanocube. Theory matches experiment for the appropriate degree of truncation. In addition, we perform simulations of a stochastic atomistic-level model incorporating a realistic description of diffusive hopping for undercoordinated surface atoms, thereby providing a visualization of the initial reshaping process.

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

First-Principles Study of n -Butane Monomolecular Cracking and Dehydrogenation on Two-Dimensional-Zeolite Model Systems: Reaction Mechanisms and Effects of Spatial Confinement

Two-dimensional (2D) ultrathin (~0.5 nm) aluminosilicate bilayer films, consisting of hexagonal prisms (a.k.a. double 6-membered rings D6R) with acidic bridging hydroxyl groups exposed on the surface, have been previously synthesized on a Ru(0001) surface as a zeolite model system. These structures are helpful for mimicking zeolite catalysts with D6R building blocks, such as chabazite. We performed density functional theory calculations to investigate the monomolecular cracking and dehydrogenation of n-butane molecules over the acidic hydroxyl groups of the 2D model system and compared the reaction energetics with that in bulk chabazite. The intrinsic activation energy barrier is the highest for dehydrogenation and lowest for central C–C bond cracking in bulk chabazite. The trend of intrinsic energy barriers for dehydrogenation and terminal and central C–C bond cracking is reproduced on the 2D aluminosilicate film. Overall, the activation barriers are higher on the 2D film than in bulk chabazite due to the lack of confinement in the former. We further explored the effects of the zeolite channel size on the n-butane adsorption and monomolecular cracking using different bulk nanoporous zeolite frameworks (TON, MEL, MEI, and VFI). We found that as the confinement of channels decreases, n-butane adsorption becomes weaker, and the intrinsic energy barrier of terminal C–C cracking increases. The activation energy barriers (dehydrogenation and terminal and central C–C cracking) on the 2D bilayer film surface, which may be considered as zeolite cages at the infinite cage size limit, are close to that in VFI with a relatively large channel size. Comparing the reaction pathway of n-butane terminal C–C cracking in 3D nanocages and on the surface of the 2D aluminosilicate film revealed that stabilizing the transition states in the 3D nanocages is responsible for the decrease in the intrinsic energy barriers for bulk zeolites.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Assessment of the literature about Be-W mixed material layer formation in the fusion reactor environment

Abstract All plasma facing surfaces in a fusion reactor, whether initially pure or an alloy, will rapidly evolve into a mixed material due to plasma-induced erosion, migration and redeposition. Beryllium (Be) erosion from the main chamber, and its transport and deposition on to a tungsten (W) divertor results in the growth of mixed Be-W layers, which can evolve to form beryllides. These Be-W mixed materials exhibit generally less desirable properties than pure tungsten or pure beryllium, such as lower melting points. In order to better understand the parameter space for growth of these alloys, this paper reviews the literature on Be-W mixed material formation experiments—in magnetically confined fusion reactors, in linear plasma test stands, and during thin-film deposition—and on computational modeling of Be-W interactions, as well as briefly assesses the Be-W growth kinetics. We conclude that the following kinetic steps drive the material mixing: adsorption of the implanted/deposited ion on the metal surface; diffusion of the implanted/deposited ion from surface into the bulk, which is accelerated by defects; and loss of deposited material through erosion. Adsorption dominates (or prevents) material mixing in thin-film deposition experiments, whereas diffusion drives material mixing in plasma exposures due to the energetic ion implantation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Assessment of the literature about Be-W mixed material layer formation in the fusion reactor environment

All plasma facing surfaces in a fusion reactor, whether initially pure or an alloy, will rapidly evolve into a mixed material due to plasma-induced erosion, migration and redeposition. Beryllium (Be) erosion from the main chamber, and its transport and deposition on to a tungsten (W) divertor results in the growth of mixed Be-W layers, which can evolve to form beryllides. These Be-W mixed materials exhibit generally less desirable properties than pure tungsten or pure beryllium, such as lower melting points. In order to better understand the parameter space for growth of these alloys, this paper reviews the literature on Be-W mixed material formation experiments—in magnetically confined fusion reactors, in linear plasma test stands, and during thin-film deposition—and on computational modeling of Be-W interactions, as well as briefly assesses the Be-W growth kinetics. We conclude that the following kinetic steps drive the material mixing: adsorption of the implanted/deposited ion on the metal surface; diffusion of the implanted/deposited ion from surface into the bulk, which is accelerated by defects; and loss of deposited material through erosion. Adsorption dominates (or prevents) material mixing in thin-film deposition experiments, whereas diffusion drives material mixing in plasma exposures due to the energetic ion implantation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Metadynamics investigation of lanthanide solvation free energy landscapes and insights into separations energetics

Lanthanide ion solvation chemistry in nonaqueous phases is key to understanding and developing effective separation processes for these critical materials. Due to the complexity and inherent disorder of the solution phase, a comprehensive picture of the solvated metal ion is often difficult to generate solely from conventional spectroscopic approaches and electronic structure calculations, particularly in the extractant phase. In this work, we use classical molecular dynamics (MD) simulation with an advanced sampling technique, metadynamics, supplemented by experimental spectroscopy and speciation analysis, to measure lanthanide solvation free energy landscapes. We define coordination-based collective variables to probe the entire range of solvation configurations in the organic phase of lanthanum (La), europium (Eu), and lutetium (Lu) nitrate salts bound with a commonly used extractant, N,N′-dimethyl, N,N′-dioctylhexylethoxymalonamide (DMDOHEMA). The known lanthanide extraction trend of La ≈ Eu > Lu is readily explained by the measured free energy surfaces, which show consistent DMDOHEMA coordination from La to Eu, followed by loss of DMDOHEMA coordination from Eu to Lu. These simulations suggest how ligand crowding at the metal center can control selectivity, in this case resulting in the opposite extraction trend as observed with other conventional extractants, where the enthalpic contribution from increasing lanthanide charge density across the series dominates the extraction energetics. We also find that the presence of inner-sphere water, verified by time-resolved fluorescence, diversifies the accessible solvation structures. As a result, understanding solvation requires consideration of an entire thermodynamic ensemble, rather than the single dominant lowest-energy structure, as is often considered out of necessity in interpretation of spectroscopic data or in electronic structure-based ligand design approaches. In general, we demonstrate how metadynamics uniquely enables investigation of complex, multidimensional solvation energetic landscapes, and how it can explain selectivity trends where extraction is controlled by more complex mechanisms than simple charge density-based selectivity.

Wang, Xiaoyu↗
DOE OSTI2024

Precision Calcination Mechanism of CaCO 3 to High-Porosity Nanoscale CaO CO 2 Sorbent Revealed by Direct In Situ Observations

Deploying energy storage and carbon capture at scale is hindered by the substantial endothermic penalty of decomposing CaCO 3 to CaO and CO 2 , and the rapid loss of CO 2 absorption capacity by CaO sorbent particles due to sintering at the high requisite decomposition temperatures. The decomposition reaction mechanism underlying sorbent deactivation remains unclear at the atomic level and nanoscale due to past reliance on postmortem characterization methods with insufficient spatial and temporal resolution. Thus, elucidating the important CaCO 3 decomposition reaction pathway requires direct observation by time-resolved (sub-)nanoscale methods. Here, chemical and structural dynamics during the decomposition of CaCO 3 nanoparticles to nanoporous CaO particles comprising high-surface-area CaO nanocrystallites are examined. Comparing in situ transmission electron microscopy (TEM) and synchrotron X-ray diffraction experiments gives key insights into the dynamics of nanoparticle calcination, involving anisotropic CaCO 3 thermal distortion before conversion to thermally dilated energetically stable CaO crystallites. Time-resolved TEM uncovered a novel CaO formation mechanism involving heterogeneous nucleation at extended CaCO 3 defects followed by sweeping reaction front motion across the initial CaCO 3 particle. These observations clarify longstanding, yet incomplete, reaction mechanisms and kinetic models lacking accurate information about (sub-)nanoscale dynamics, while also demonstrating calcination of CaCO 3 without sintering through rapid heating and precise temperature control.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Interaction–Deletion: A Composite Energy Method for the Optimization of Molecular Systems Selectively Removing Specific Nonbonded Interactions

The complex interactions between different portions of a large molecule can be challenging to analyze through traditional electronic structure calculations. Moreover, standard methods cannot easily quantify the physical consequences of individual pairwise interactions inside a molecule. By creating a set of molecular fragments, we propose a composite energy method to explore changes in a molecule caused by removing selected nonbonded interactions between different molecular portions. Energies and forces are easily obtained with this composite approach, allowing geometry optimizations that lead to chemically meaningful structures that describe how the omitted interactions contribute to changes in the local geometrical minima. We illustrate the application of our new hybrid scheme by computing the influence of intramolecular hydrogen-bonding interactions in two small molecules: 1,6-(tG + G + TG + G + g – )-hexanediol and a cyclic analogue, cis-1,4-cyclohexanediol. The resulting structural and energetic changes are interpreted to yield key physical insights and quantify concepts such as “preparation energy” or “reorganization energy”. We demonstrate that the composite method can be extended to larger molecular systems by showing its application on a Si(100) surface model containing interactions between dissociated ammonia molecules on adjacent surface dimers. The scheme’s efficacy is also tested by applying it to systems having multiple intramolecular interactions, viz., 310-polyglycine and H+GPGG. Furthermore, the cooperative nature of intramolecular hydrogen bonds is explored by using interaction–deletion in 2-nitrobenzene-1,3-diol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2023

Assessing Implicit and Explicit Polarizable Solvation Models for Nuclear–Electronic Orbital Systems: Quantum Proton Polarization and Solvation Energetics

Accurate simulations of many chemical processes require the inclusion of both nuclear quantum effects and a solvent environment. The nuclear–electronic orbital (NEO) approach, which treats electrons and select nuclei quantum mechanically on the same level, combined with a polarizable continuum model (PCM) for the solvent environment, addresses this challenge in a computationally practical manner. In this work, the NEO-PCM approach is extended beyond the IEF-PCM (integral equation formalism PCM) and C-PCM (conductor PCM) approaches to the SS(V)PE (surface and simulation of volume polarization for electrostatics) and ddCOSMO (domain decomposed conductor-like screening model) approaches. IEF-PCM, SS(V)PE, C-PCM, and ddCOSMO all exhibit similar solvation energies as well as comparable nuclear polarization within the NEO framework. The calculations show that the nuclear density does not leak out of the molecular cavity because it is much more localized than the electronic density. Finally, the polarization of quantized protons is analyzed in both continuum solvent and explicit solvent environments described by the polarizable MB-pol model, illustrating the impact of specific hydrogen-bonding interactions captured only by explicit solvation. Furthermore, these calculations highlight the relationship among solvation formalism, nuclear polarization, and energetics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Unraveling the transformation pathway of the 𝛽 to 𝛾 phase transition in Ga 2 ⁢O 3 from atomistic simulations

Defect spinel 𝛾−Ga 2 ⁢O 3 is the least stable polymorph of Ga 2 ⁢O 3 , so its frequent appearance as a structural defect within or on the surface of monoclinic 𝛽−Ga 2 ⁢O 3 remains a mystery. Through first-principles calculations, we explore potential pathways for the phase transition from 𝛽−Ga 2⁢ O 3 to 𝛾−Ga 2 ⁢O 3 , and examine two key driving forces: tensile strain and Ga deficiency. When configurational entropy contributions to phase energies are included, the 𝛾 phase becomes energetically competitive with the 𝛽 phase, with the free energy difference between these phases diminishing even further under Ga-deficient conditions. Notably, a stability crossover occurs at room temperature at high vacancy concentrations ([V$^{3−}_{Ga}$]>3%) . A simple model 𝛽 → 𝛾 transformation pathway is identified, comprising two primary reactions, that enables the formation of the 𝛾 phase via simultaneous migration of Ga atoms from tetrahedral lattice sites to octahedral interstitial positions. The transformation barriers are prohibitively large in pristine Ga 2 ⁢O 3 , but can be substantially reduced by: (1) the presence of Ga vacancies, (2) elongational strains along the crystallographic 𝑎-axis, and (3) when volumetric relaxations are possible during transformation. These results elucidate prior experimental observations, where 𝛾−Ga 2⁢ O 3 is seen on damaged surfaces or in highly 𝑛-type 𝛽−Ga 2⁢ O 3 environments, which support Ga deficiency and mechanical strain. The insights into the driving forces and mechanisms of 𝛾−Ga 2⁢ O 3 formation enhance understanding of how localized strain and nonequilibrium defect concentrations may facilitate its formation from the 𝛽 phase.

Defects↗
DOE OSTI2022

Ligand removal energetics control CO 2 electroreduction selectivity on atomically precise, ligated alloy nanoclusters

Atomically precise, thiolate-protected gold nanoclusters (TPNCs) exhibit remarkable catalytic performance for the electrochemical reduction of carbon dioxide (CO 2 R) to CO. The origin of their high CO 2 R activity and selectivity has been attributed to partial ligand removal from the thiolate-covered surfaces of TPNCs to expose catalytically active sulfur atoms. Recently, heterometal doped (alloy) TPNCs have been shown to exhibit enhanced CO 2 R activity and selectivity compared to their monometallic counterparts. However, systematic studies on the effect of doping (metal type and location on TPNC) on active site exposure and CO 2 R activity are missing in literature. Herein, we apply Density Functional Theory calculations to investigate the effect of heterometal (Pt, Pd, Hg and Cd) doping of Au 25 (SR) 18 TPNC on the active site exposure and CO 2 R activity and selectivity. In this study, we reveal that doping significantly modifies relevant TPNC electronic properties, such as electron affinity, while also altering partial ligand removal and carboxyl (*COOH) intermediate formation energies. Furthermore, we demonstrate that changing the dopant (e.g. Hg) position can change the selectivity of the TPNC towards CO (g) or H 2(g) formation, highlighting the importance of dopant locations in TPNC-based CO 2 R. Most notably, we report a universal (i.e. capturing different dopant types and positions) linear trend between the ligand removal energy and i) the *COOH formation energy, as well as, ii) the hydrogen (*H) formation energy on the different alloy TPNCs. Thus, utilizing the ligand removal energy as a descriptor for CO 2 RR activity and selectivity, our work opens new avenues for accelerated computational screening of different alloy TPNCs for electrocatalytic CO 2 R applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Boosting CO2R Performance of Ag Electrocatalysts by Sulfur-Doped Carbon Support

We find that S-doped carbon support can boost the CO2 reduction (CO2R) performance of Ag electrocatalysts. Firstly, surface science enabled electrocatalysis showed that Ag supported on S-doped highly oriented pyrolytic graphite (HOPG), a model electrocatalyst, demonstrated 100% higher CO turnover frequency (TOFCO = 3.6 ± 0.2 CO/atomAg/s) than that supported on S-free HOPG (TOFCO = 1.8 ± 0.2 CO/atomAg/s). Computational modeling based on density functional theory (DFT) revealed a more stabilized *COOH intermediate on Ag supported on S-doped carbon and thus a more favorable energetic pathway of CO2-to-CO, consistent with experimental results from the model electrocatalysts studies. Finally, this proof of concept was translated to the synthesis of powder electrocatalyst with 2 wt% Ag supported on S-doped carbon black, demonstrating > 40-fold high CO mass activity than a commercial Ag cathode with steady FECO ~ 96% at 100 mA/cm2 for 50 hours of continuous operation in a gas diffusion electrode (GDE) electrolyzer. For comparison, 2 wt% Ag supported on carbon black without S- doping showed a maximum FECO ~ 70% at 100 mA/cm2. This work demonstrates a successful bottom-up design of CO2R electrocatalysts guided by surface science enabled electrocatalysis.

CO2 conversion↗
DOE OSTI2023

Mechanism of Stoichiometrically Governed Titanium Oxide Brownian Tree Formation on Stepped Au(111)

Previously observed formation of substoichiometric titanium oxide dendritic structures across terraces of Au(111) is computationally studied and shown to follow the classical fractal formation mechanism of diffusion-limited aggregation (DLA). Globally optimized gas-phase oxide cluster structures are sampled in a variety of landing formations on gold surfaces and shown to favor isomers driving polymerization to Brownian tree fractal structures. Mobility of Ti 3 O 5 monomers is shown to be extremely high, with diffusion barriers of 0.21 eV or less. Through bonding stabilization, polymerization of these monomers is energetically favorable and irreversible on the 111 terrace but geometrically impossible to propagate along the step edge. Simulated scanning tunneling microscopy (STM) images show strong similarity to experiment. By contrast, observation of Ti 3 O 6 aggregating as wires along step edges is explained by the affinity of oxygen to step edges and statistical arguments for aggregation entropy at the step, in addition to low barriers for monomer diffusion and polymerization.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Toward Understanding and Controlling Organic Reactions on Metal Oxide Catalysts

Metal oxides have structurally complex surfaces on which a variety of adsorption site types can occur, including cation sites, anion sites, oxygen vacancy sites, and Brønsted acid sites. These sites can catalyze the catalytic transformation of organic molecules via diverse routes, thus enabling H abstraction, O abstraction, C–C bond formation, and other reactions. This Perspective provides an update on recent advances and future directions for various organic reactions on metal oxide catalyst surfaces, particularly for C–H activation of alkanes and for C–C bond formation with organic oxygenate reactants. Here, we put emphasis on the molecular scale details, on the active site structures required to enable the formation of kinetically relevant transition states, energetic descriptors, as well as contemporary ideas to enable low activation energies. This progress has been enabled by specialized experiments and the increased capabilities of modern electronic structure calculations.

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

Observing cosmic-ray extensive air showers with a silicon imaging detector

Extensive air showers induced from high-energy cosmic rays provide a window into understanding the most energetic phenomena in the universe. We present a new method for observing these showers using the silicon imaging detector Subaru Hyper Suprime-Cam (HSC). This method has the advantage of being able to measure individual secondary particles. When paired with a surface detector array, silicon imaging detectors like Subaru HSC will be useful for studying the properties of extensive air showers in detail. The following report outlines the first results of observing extensive air showers with Subaru HSC. The potential for reconstructing the incident direction of primary cosmic rays is demonstrated and possible interdisciplinary applications are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Multiple redox mechanisms for water-gas shift reaction on Fe 3 O 4 (1 1 1) surface: A density functional theory and mean-field microkinetic modeling study

In this work, we investigated redox mechanisms for the high temperature (HT) water–gas shift reaction (WGSR) on the Fe 3 O 4 (1 1 1) surface using density functional theory (DFT) calculations and mean-field microkinetic modeling (MKM). The redox pathways branched into three Langmuir-Hinshelwood processes (LH1, LH2 and LH3) and one Mars-van-Krevelen (MVK) process (in the presence of oxygen vacancy) from multidentate binding of CO and CO 2 over four available reactive sites (Fe oct2 , Fe tet1 , Fe bridge, and O1). We found that the LH1 and LH2 processes have CO adsorption at a single iron site (Fe oct2 or Fe tet1 ), while the LH3 and MVK processes have stronger chemisorption of CO or CO 2 by both the Fe oct2 site and the O 1 site. From the mechanistic study of these reaction path, we recognized that availability of O 1 sites was key to proceed to either CO oxidation by the single site (LH1, LH2) or by the dual site (LH3, MVK). We observed that the reaction energetics were significantly different in the CO oxidation steps, where the single or dual site results in markedly different apparent activation energy (73–281 kJ/mol) and reaction rates (10 -3 –10 -10 mol∙m -2 ∙s -1 ) among the four reaction mechanisms. The utilization of mean field MKM with DFT reaction energetics helps explain the experimental debates for the catalytic reaction details as well as providing a possible direction to engineer the catalyst with higher activity.

42 ENGINEERING↗