Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles density functional theory”

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

First-Principles Study of Ferroelectric and Optical Properties in Derivatives of Thiourea

Organic ferroelectric materials in the form of molecular crystals are promising alternatives to well-known ferroelectric perovskite oxides due to their structural flexibility, tunability, and ease of processing. First-principles density functional theory calculations combined with group theory analyses are employed in this work to discuss polar nature of a class of derivatives of thiourea: CH 4 N 2 A where A = O, S, Se, and Te. Further, optical properties in the UV range along with thermodynamic and dynamic stability of such materials are also evaluated in the study. We find that CH 4 N 2 Te with bandgap energy of ~3.33 eV shows large spontaneous polarization of P = 10.92 μC/cm 2 . The electronic and optical properties become tunable under strain.

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

Defect Properties, Anion Ordering, and Photochromic Mechanism in Yttrium Oxyhydride

Yttrium oxyhydride (YHO) undergoes a reversible photochromic transition when exposed to ultraviolet light. However, the mechanism for this transformation is not fully understood, and the structure and precise chemical composition of YHO remain under debate. Here, we use first-principles density functional theory calculations with a hybrid functional to study the structure, chemical stability, and point defect properties of YHO. As experiments have shown, we find that YHO prefers a cubic structure, with H and O anions present in equal concentrations and located on tetrahedral sites. Stoichiometric and ordered YHO is chemically stable, but it has a wide band gap of 5.01 eV, considerably larger than that measured in experiments (2.4–3.8 eV). On the other hand, Y4H10O has a smaller band gap of 2.97 eV and also has a region of chemical stability; thus, the actual material may include some fraction of this H-rich structure. The defect chemistry of YHO is dominated by anionic antisite species (H O and O H ), with hydrogen interstitials (H i ) and vacancies (V H ) also present in reasonably high concentrations. We show that antisite disorder lowers the band gap relative to the perfectly ordered structure, bringing the magnitude of the gap into closer agreement with experiment. Based on our calculations of defect migration and the positions of defect states relative to the band edges, we link the onset of photochromic behavior to the reaction H O – → V O 0 + H i – , which follows photoexcitation of a H O + defect. H i – can subsequently migrate away and be trapped by additional H O + defects, contributing to the persistence of the reaction, while the resultant oxygen vacancy, V O 0 , introduces an occupied defect state that leads to optical absorption at visible wavelengths. Our results can explain reported discrepancies between experimental and computational results for YHO, and they allow us to propose specific atomic-scale processes that can lead to photochromism. In conclusion, understanding these mechanisms is key for unlocking YHO’s application in devices ranging from smart windows and optoelectronics to electrochemical synapses for neural networks.

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

Atomic-Scale Insights into Carbon Dissolution in α-, γ-, and θ-Al 2 O 3 : Phase-dependent Transport Dynamics from First-Principles Calculations

α-Al 2 O 3 exhibits superior carburizing corrosion resistance compared to metastable γ-Al 2 O 3 and θ-Al 2 O 3 phases in high-temperature CO 2 environments, yet its atomic-scale origins remain unclear. Using first-principles density functional theory, we systematically investigate carbon dissolution and diffusion in α-Al 2 O 3 , γ-Al 2 O 3 , and θ-Al 2 O 3 , including the effects of oxygen (O) and aluminum (Al) vacancies. Our results show that α-Al 2 O 3 consistently exhibits higher carbon solution enthalpies than γ-Al 2 O 3 and θ-Al 2 O 3 in both pristine and defective structures, indicating lower intrinsic carbon solubility in α-Al 2 O 3 . Vacancies significantly enhance carbon incorporation: O vacancies reduce solution enthalpy, while Al vacancies further amplify this effect, with a strong preference for carbon at Al vacancy sites. Carbon diffusion barriers are also highest in α-Al 2 O 3 , reflecting slower carbon mobility. Al vacancies increase diffusion barriers across all phases, while O vacancies raise barriers in α- and γ-Al 2 O 3 but slightly lower them in θ-Al 2 O 3 . These results reveal a dual mechanism behind the carburizing resistance of α-Al 2 O 3 : reduced carbon solubility and elevated diffusion barriers. Furthermore, this work provides atomic-scale insights to guide the design of alumina-based materials with improved carburizing resistance through phase selection and defect engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An Atomistic Study of Reactivity in Solid-State Electrolyte Interphase Formation for Li/Li7P3S11

Lithium metal batteries offer superior volumetric and gravimetric specific capacities compared to those based on traditional graphite anodes. Although advancements in solid-state electrolytes address safety concerns, challenges remain, particularly regarding interphase formation in lithium metal anodes. This work presents a computational framework based on high-throughput first-principles density functional theory and machine-learning interatomic potentials (MLIPs) including automated iterative, active learning to enable robust computational exploration of interphase formation between lithium metal anodes and an inorganic solid-state electrolyte. As a demonstration, we apply the framework to a Li/Li7P3S11 interface and find that it accurately identifies the experimentally observed, thermodynamically stable interphase products as well as their overall spatial arrangement within a heterogeneous, amorphous layered structure, with Li2S domains of nanocrystallinity. Our simulations show two stages, a fast and slow diffusion reaction regime, that corroborate the relative phase formation rate of Li x P, Li2S, and Li3P. Using the Onsager transport theory, we capture time-dependent ionic diffusion within the reacting interface, including cross-correlation effects. We found that cross-correlation effects between Li-P and P-S ionic motion significantly influence P-ion diffusion, making it highly sensitive to the local environment and potentially leading to "kinetic trapping" of Li-P phases. The passivation of the interface is shown as the ionic fluxes all approach zero, effectively halting interphase growth.

Diffusion↗
DOE OSTI2020

Thermodynamic Preference for Atom Adsorption on versus Intercalation into Multilayer Graphene

The thermodynamic preference of a foreign atom for adsorption on versus intercalation into a graphitic surface is of fundamental and widespread interest. From an exhaustive first-principles density functional theory investigation for 38 typical elements over the periodic table, we reveal a quasilinear correlation between the Shannon effective ionic radius and the chemical-potential difference for a single atom from adsorption to intercalation at multilayer graphene surfaces. A critical Shannon radius is found to be around 0.10 nm, below (above) which intercalation (adsorption) is more favorable for elements with ionic-like bonding after intercalation. Single atoms with van der Waals-biased bonding show some deviation from the linear relationship, while single atoms for the elements with covalent-like bonding do not favor intercalation relative to adsorption. An energy decomposition analysis indicates that the chemical-potential difference determining the thermodynamic preference of a foreign atom for adsorption versus intercalation results from the competition between the electronic and elastic strain effects.

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

Thermodynamically Driven Formation of Intercalated Cu Carpets from Supported Cu Pyramids on MoS 2

Thermodynamic and kinetic analyses based on our first-principles density functional theory calculations are used to interpret the experimentally observed formation of Cu carpets intercalated under the top layer of a 2H-MoS 2 substrate. Spontaneous Cu transport from Cu pyramids on top of the MoS 2 substrate through surface point defects to the growing Cu carpet is shown to be driven by a slightly lower chemical potential for the Cu carpet. We demonstrate that the competition between a preference for a thicker Cu carpet and the cost of elastic stretching of the top MoS 2 layer results in a selected Cu carpet thickness. Furthermore, we also propose that Cu transport occurs primarily via vacancy-mediated diffusion through constricting point defect portals.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Degeneracy in Intercalated Pb Phases under Buffer-Layer Graphene on SiC(0001) and Diffuse Moiré Spots in Surface Diffraction

First-principles density functional theory (DFT) is used to analyze the stability of Pb intercalated phases under buffer layer graphene on SiC(0001) as a function of the supercell size, Pb coverage, and degree of Pb ordering. By comparing the chemical potentials of such two-dimensional Pb structures, here we find that there is a family of structurally distinct thermodynamically preferred Pb subsurface configurations with minute stability differences. These differences are comparable to the thermal energies at about 450 °C, where the Pb intercalated phases are grown. High-resolution surface-diffraction experiments using Spot Profile Analysis Low-Energy Electron Diffraction (SPA-LEED) confirm this high degree of degeneracy of the Pb intercalated phases from broad, low-intensity moiré spots observed exclusively from intercalated Pb. The low intensity of the moiré spots implies the coexistence of structurally different subsurface Pb phases.

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

Engineering Small HOMO–LUMO Gaps in Polycyclic Aromatic Hydrocarbons with Topologically Protected States

Topological phases in laterally confined low-dimensional nanographenes have emerged as versatile design tools that can imbue otherwise unremarkable materials with exotic band structures ranging from topological semiconductors and quantum dots to intrinsically metallic bands. The periodic boundary conditions that define the topology of a given lattice have thus far prevented the translation of this technology to the quasi-zero-dimensional (0D) domain of small molecular structures. Here, we describe the synthesis of a polycyclic aromatic hydrocarbon (PAH) featuring two localized zero modes (ZMs) formed by the topological junction interface between a trivial and nontrivial phase within a single molecule. First-principles density functional theory calculations predict a strong hybridization between adjacent ZMs that gives rise to an exceptionally small HOMO–LUMO gap. Scanning tunneling microscopy and spectroscopy corroborate the molecular structure of 9/7/9-double quantum dots and reveal an experimental quasiparticle gap of 0.16 eV, corresponding to a carbon-based small molecule long-wavelength infrared (LWIR) absorber.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Kinetic Isotope Effect in the Radical Step-Growth Termination of On-Surface Synthesized Graphene Nanoribbons

Bottom-up on-surface synthesis has emerged as a versatile tool to access and finely tune the electronic structure of nanographenes. The controlled generation of reactive intermediates catalyzed and stabilized by a supporting substrate has enabled the design, assembly, and characterization of a wide range of exotic tailor-made quantum materials. Even under these tightly controlled conditions, the growth of extended structures remains limited by termination processes and undesired side reactions. Here, we identify an H atom transfer as one principal contributor to the radical step growth termination that limits the on-surface growth of N = 7 armchair graphene nanoribbons (7-AGNRs) on Au(111) surfaces. Analysis of 7-AGNR lengths grown from protiated and deuterated molecular precursors reveals a primary kinetic isotope effect of KIE ∼ 1.4. First-principles density functional theory calculations suggest that a concerted H atom transfer mechanism that involves the breaking of a C–H/D bond in the transition state is associated with radical chain termination.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

ALD-Derived WO 3– x Leads to Nearly Wake-Up-Free Ferroelectric Hf 0.5 Zr 0.5 O 2 at Elevated Temperatures

Breaking the memory wall in advanced computing architectures will require complex 3D integration of emerging memory materials such as ferroelectrics─either within the back-end-of-line (BEOL) of CMOS front-end processes or through advanced 3D packaging technologies. Achieving this integration demands that memory materials exhibit high thermal resilience, with the capability to operate reliably at elevated temperatures, such as 125°C, due to the substantial heat generated by front-end transistors. However, silicon-compatible HfO 2 -based ferroelectrics tend to exhibit antiferroelectric-like behavior in this temperature range, accompanied by a more pronounced wake-up effect, posing significant challenges to their thermal reliability. Here, we report that by introducing a thin tungsten oxide (WO 3–x ) layer─known as an oxygen reservoir─and carefully tuning its oxygen content, ultrathin Hf 0.5 Zr 0.5 O 2 (5 nm) films can be made robust against the ferroelectric-to-antiferroelectric transition at elevated temperatures. This approach not only minimizes polarization loss in the pristine state but also effectively suppresses the wake-up effect, reducing the required wake-up cycles from 10 5 to only 10 at 125°C, a qualifying temperature for back-end memory integrated with front-end logic, as defined by the JEDEC standard. First-principles density functional theory (DFT) calculations reveal that WO 3 enhances the stability of the ferroelectric orthorhombic phase (o-phase) at elevated temperatures by increasing the tetragonal-to-orthorhombic phase energy gap and promoting favorable phonon mode evolution, thereby supporting o-phase formation under both thermodynamic and kinetic constraints.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Design of Graphene/Ionic Liquid Composites for Carbon Capture

Pore size is a crucial factor impacting gas separation in porous separation materials, but how to control the pore size to optimize the separation performance remains a challenge. Here, we propose a design of graphene/ionic liquid composites with tunable slit pore sizes, where cations and anions of ionic liquids are intercalated between graphene layers. By varying the sizes of the ions, we show from first-principles density functional theory calculations that the accessible pore size can be tuned from 3.4 to 6.0 Å. Grand canonical Monte Carlo simulations of gas sorption find that the composite materials possess high CO 2 uptake at room temperature and 1 bar (up to ~8.5 mmol/g). Further simulations of the sorption of gas mixtures reveal that high CO 2 /N 2 and CO 2 /CH 4 adsorption selectivities can be obtained when the accessible pore size is <5 Å. This work suggests a new strategy to achieve tunable pore sizes via the graphene/IL composites for highly selective CO 2 /N 2 and CO 2 /CH 4 adsorption.

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

Raman Digital Twin of Monolayer Janus Transition Metal Dichalcogenides

Monolayer transition metal dichalcogenides (TMDs) are a key class of two-dimensional (2D) materials with broad technological potential. Their Janus counterparts exhibit unique properties due to broken out-of-plane symmetry and further enrich the functionalities of TMDs. However, experimental synthesis and identification of Janus TMDs remain challenging. It is thus highly desirable to have a rapid, simple, and in situ characterization technique to monitor, in real time, the conversion process from the parent to Janus structure. Raman spectroscopy stands out for such a task as it is a powerful, nondestructive, and very commonly used tool to characterize 2D materials both in situ and ex situ. To realize the full potential of Raman spectroscopy on rapid characterization of Janus TMDs, we present a computational “Raman digital twin” library for various monolayer Janus TMDs in both 2H and Td phases. We focus on group-6 TMDs: MoS 2 , WS 2 , MoSe 2 , WSe 2 , MoTe 2 , WTe 2 and their Janus variants: MoSSe, MoSTe, MoSeTe, WSSe, WSTe, and WSeTe. Using first-principles density functional theory (DFT), we calculate their vibrational properties and predict distinct Raman fingerprints. These phonon and Raman signatures reflect each material’s structural symmetry and atomic composition, enabling clear identification via Raman spectroscopy. Our theoretical work supports experimental efforts by providing benchmarks for material identification, structural analysis, and quality control. In conclusion, the computational library expedites the discovery and development of Janus 2D materials, facilitating tighter integration between theoretical predictions and experimental validation.

Chalcogenides↗
DOE OSTI2026

2D in-Plane Ordered MXene Nanosheets Derived from (Mo 2/3 Er 1/3 ) 2 AlC Rare-Earth i-MAX for Energy Storage Applications

MXenes have become one of the most versatile families of two-dimensional (2D) materials due to their high conductivity, hydrophilicity, and remarkable electrochemical performance. This has stimulated intense efforts to design and synthesize MXenes, including structurally unique in-plane ordered 2D MXenes called i-MXenes. Here, we have synthesized the quaternary rare earth (RE)-based i-MAX phase (Mo 2/3 Er 1/3 ) 2 AlC using an arc melting method, and the corresponding 2D i-MXene was then obtained through a LiF/HCl soft etching process. Literature studies have shown that Al and the RE element are etched out during the etching process, leading to the formation of pure vacancy-ordered Mo1.33C 2D i-MXene. However, our investigation reveals that upon exposure to a fluorine solution, the i-MAX phase forms RE fluoride impurities, which are challenging to remove through HCl−DI water washing and persist in the final product, resulting in impure Mo 1.33 C@Er i-MXene. These results were confirmed by various characterizations such as X-ray diffraction, Raman spectroscopy, X-ray photoelectron spectroscopy, and scanning transmission electron microscopy. Although the Mo 1.33 C@Er electrode showed a 24-fold increase in specific capacitance compared to its parent i-MAX phase, it still exhibited a high charge-transfer resistance arising from the insulating nature of RE fluoride byproducts, which adversely influence the overall capacitance behavior of the synthesized 2D Mo 1.33 C@Er i-MXenes. This study contributes to identifying pathways for the preparation of pure 2D i-MXenes from RE-based i-MAX phases and developing improved synthesis methods. With additional process optimization, the 2D i-MXene holds a strong potential for electrochemical energy storage applications. Additionally, the electronic structures of Mo 1.33 C were theoretically studied using first-principles density functional theory calculations, which revealed that pristine Mo 1.33 C is metallic, and this metallic nature is preserved even with −O, −F, and mixed functionalization.

chemical structure↗
DOE OSTI2020

Formic Acid: A Hydrogen-Bonding Cocatalyst for Formate Decomposition

Hydrogen bonding accelerates many catalytic reactions by orienting intermediates, stabilizing transition states, and even opening reaction pathways. However, most mechanistic studies regarding the decomposition of formic acid (FA), a promising hydrogen storage material, neglect hydrogen-bonding interactions even though FA is a strong hydrogen-bond donor and acceptor. Here, we probe the formation of bimolecular hydrogen-bonded complexes between FA and formate (FA–HCOO complexes) adsorbed on metal surfaces and how these complexes affect HCOO* decomposition. Using first-principles density functional theory (DFT) calculations on 12 close-packed (111)/(0001) and 8 open (100) surfaces of 12 transition metals—Ag, Au, Co, Cu, Ir, Ni, Os, Re, Pd, Pt, Rh, and Ru, we—show that FA–HCOO complexes are generally thermodynamically stable, even at elevated temperatures and pressures. We then illustrate that these complexes produce infrared spectroscopic signatures consistent with as yet unassigned experimental peaks. We last demonstrate that by stabilizing the dangling bond of monodentate HCOO*, these complexes significantly lower the barriers for rotation of HCOO* from a bidentate to a monodentate configuration, the rate-limiting step for HCOO* decomposition on many surfaces. FA thus acts as a cocatalyst for HCOO* decomposition. Our results may guide the community toward improved catalysts for reactions involving HCOO* such as FA decomposition, methanol steam reforming, and the water gas shift reaction. More broadly, our work highlights the ability of hydrogen bonding to modify the adsorbed structures of intermediates and lower the barriers for their reaction on heterogeneous catalysts. Lastly, this phenomenon can be relevant for other reactions involving ammonia, alcohols, and carboxylic acids.

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

Insights into the Oxygen Evolution Reaction on Graphene-Based Single-Atom Catalysts from First-Principles-Informed Microkinetic Modeling

Single-atom transition metals embedded in nitrogen-doped graphene have emerged as promising electrocatalysts due to their high activity and low material cost. These materials have been shown to catalyze a variety of electrochemical reactions, but their active sites under reaction conditions remain poorly understood. Using first-principles density functional theory calculations, we develop a pH-dependent microkinetic model to evaluate the relative performance of transition metal catalysts embedded in fourfold N-substituted double carbon vacancies in graphene for the oxygen evolution reaction. We find that reaction pathways involving intermediates co-adsorbed on the metal site are preferred on all transition metals. These pathways lead to enhancements in catalytic activity and broaden the activity peak when compared with purely thermodynamics-based predictions. Furthermore, these findings demonstrate the importance of investigating reaction pathways on graphene-based catalysts and other two-dimensional (2D) materials that involve metal active centers decorated by spectator intermediate species.

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

Acetylene Semi-Hydrogenation on a Perovskite Oxyhydride Surface: Insights from First Principles and Microkinetic Modeling

Although perovskite oxyhydrides (POHs) have shown promise as a support for selective hydrogenation, it is unknown if they themselves can catalyze alkyne semi-hydrogenation. Here, we use first-principles density functional theory, coupled with microkinetic modeling, to investigate acetylene semi-hydrogenation on a prototypical POH, cubic BaTiO 2.5 H 0.5 (BTOH). Two different mechanisms are examined on a representative surface of BTOH under the reaction conditions: both are based on the Horiuti–Polanyi mechanism, but the way of H 2 dissociation is different. In mechanism 1, a lattice hydride H atom and then a surface adsorbed H atom sequentially hydrogenate the adsorbed acetylene. In mechanism 2, two lattice hydride H atoms from the BTOH sequentially hydrogenate the adsorbed acetylene. In both mechanisms, the H atoms are replenished from gas phase H 2 dissociation. Using density functional theory (DFT), we have calculated pathways for these two mechanisms in hydrogenation of acetylene to ethylene, as well as further to ethane. Microkinetic modeling based on the DFT energetics indicates that at 523 K, hydrogenation to ethylene and ethane via mechanism 1 occurs at a rate up to two orders of magnitude faster than via mechanism 2. A selectivity analysis for the temperature range of 373–673 K shows that the product observed is essentially only ethylene within the more active mechanism 1. This is because at the steady state, there is a significant amount of surface anion vacancies that facilitate heterolytic H 2 dissociation and help stabilize a vinyl intermediate. In conclusion, the present findings suggest that POHs are capable of selective hydrogenation, thanks to their rich hydrogen surface-redox chemistry.

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

Pt Nanoparticles on Atomic-Metal-Rich Carbon for Heavy-Duty Fuel Cell Catalysts: Durability Enhancement and Degradation Behavior in Membrane Electrode Assemblies

Proton exchange membrane fuel cells (PEMFCs) are a promising zero-emission power source for heavy-duty vehicles (HDVs). However, long-term durability of up to 25,000 h is challenging because current carbon support, catalyst, membrane, and ionomer developed for traditional light-duty vehicles cannot meet the stringent requirement. Therefore, understanding catalyst degradation mechanisms under the HDV condition is crucial for rationally designing highly active and durable platinum group metal (PGM) catalysts for high-performance membrane electrode assemblies (MEAs). Herein, we report a PGM catalyst consisting of platinum nanoparticles with a high content (40 wt %) on atomic-metal-site (e.g., MnN 4 )-rich carbon support. MEAs with the Pt (40 wt %)/Mn–N–C cathode catalyst achieved significantly enhanced performance and durability, generating 1.41 A cm –2 at 0.7 V under HDV conditions (0.25 mgPt cm –2 and 250 kPa abs pressure) and retaining 1.20 A cm –2 after an extended and accelerated stress test up to 150,000 voltage cycles. Electron microscopy studies indicate that most fine Pt nanoparticles are retained on or/and in the carbon support covered with the ionomer throughout the catalyst layer at the end of life. During the long-term stability test, the observed electrochemical active surface area reduction and performance loss primarily result from Pt depletion in the catalyst layer due to Pt dissolution and redeposition at the interface of the cathode and membrane. Importantly, the first-principle density functional theory calculations further reveal a support entrapment effect of the Mn–N–C, in which the MnN 4 site can specifically adsorb the Pt atom and further retard the Pt dissolution and migration, therefore enhancing long-term MEA durability.

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

Unveiling the Electrocatalytic Hydrogen Evolution Reaction Pathway on RuP 2 through Ab Initio Grand Canonical Monte Carlo

In this study, the high catalytic reactivity of ruthenium phosphide (RuP 2 ) has been identified by first-principles density functional theory (DFT) calculations for the electrocatalytic hydrogen evolution reaction (HER). Complex surface reconstructions are considered by applying the ab initio grand canonical Monte Carlo (ai-GCMC) algorithm, efficiently providing a sufficient phase-space exploration of possible surfaces. Combined with surface-phase Pourbaix diagrams, we are able to identify the actual surfaces that obtained under specific experimental environments, thus leading to a more accurate understanding of the nature of the active sites and the binding strength of adsorbates. Specifically, through hundreds of surface reconstructions and hydrogenation states generated with ai-GCMC, we identify the most favorable surface phases of RuP 2 under aqueous acidic conditions. We discover that the HER activity is determined by multiple surfaces with different stoichiometries within a narrow electrode potential window. Low HER overpotential (η) has been found for each of the identified surfaces, as low as 0.04 V. High H-coverage reconstructed surfaces have been discovered under acidic conditions, and the surface Ru sites introduced by additional Ru adatoms or exposed by P-vacancies serve as the active sites for HER based on their nearly reversible H binding. Furthermore, this work provides atomistic insights into the origin of high HER activity on RuP 2 by exploring the dynamic surface phases of electrocatalysts and features a generalizable method to explore the reconstructed/hydrogenated surface space as a function of experimental conditions.

25 ENERGY STORAGE↗