Search NASA⌕ Search

SEARCH · Search NASA

Results for “embedding 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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Anisotropic interactions for continuum modeling of protein–membrane systems

In this work, a model for anisotropic interactions between proteins and cellular membranes is proposed for large-scale continuum simulations. The framework of the model is based on dynamic density functional theory, which provides a formalism to describe the lipid densities within the membrane as continuum fields while still maintaining the fidelity of the underlying molecular interactions. Within this framework, we extend recent results to include the anisotropic effects of protein–lipid interactions. As applications, we consider two membrane proteins of biological interest: a RAS–RAF complex tethered to the membrane and a membrane embedded G protein-coupled receptor. A strong qualitative and quantitative agreement is found between the numerical results and the corresponding molecular dynamics simulations. Combining the scope of continuum level simulations with the details from molecular level particle simulations enables research into protein–membrane behaviors at a more biologically relevant scale, which crucially can also be accessed via experiment.

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

Efficient 3D Molecular Design with an E(3) Invariant Transformer VAE

This work introduces a three-dimensional (3D) invariant graph-to-string transformer variational autoencoders (VAE) (Vagrant) for generating molecules with accurate density functional theory (DFT)-level properties. Vagrant learns to model the joint probability distribution of a 3D molecular structure and its properties by encoding molecular structures into a 3D-aware latent space. Directed navigation through this latent space implicitly optimizes the 3D structure of a molecule, and the latent embedding can be used to condition a generative transformer to predict the candidate structure as a one-dimensional (1D) sequence. Additionally, we introduce two novel sampling methods that exploit the latent characteristics of a VAE to improve performance. Here, we show that our method outperforms comparable 3D autoregressive and diffusion methods for predicting quantum chemical property values of novel molecules in terms of both sample quality and computational efficiency.

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

Multifidelity Neural Network Formulations for Prediction of Reactive Molecular Potential Energy Surfaces

Here, this paper focuses on the development of multifidelity modeling approaches using neural network surrogates, where training data arising from multiple model forms and resolutions are integrated to predict high-fidelity response quantities of interest at lower cost. We focus on the context of quantum chemistry and the integration of information from multiple levels of theory. Important foundations include the use of symmetry function-based atomic energy vector constructions as feature vectors for representing structures across families of molecules and single-fidelity neural network training capabilities that learn the relationships needed to map feature vectors to potential energy predictions. These foundations are embedded within several multifidelity topologies that decompose the high-fidelity mapping into model-based components, including sequential formulations that admit a general nonlinear mapping across fidelities and discrepancy-based formulations that presume an additive decomposition. Methodologies are first explored and demonstrated on a pair of simple analytical test problems and then deployed for potential energy prediction for C 5 H 5 using B2PLYP-D3/6-311++G(d,p) for high-fidelity simulation data and Hartree–Fock 6-31G for low-fidelity data. For the common case of limited access to high-fidelity data, our computational results demonstrate that multifidelity neural network potential energy surface constructions achieve roughly an order of magnitude improvement, either in terms of test error reduction for equivalent total simulation cost or reduction in total cost for equivalent error.

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

High pressure melt line of nickel using a generalized embedded atomic method potential

As the second most abundant metal in the Earth's core, nickel plays an important role in determining the structure and temperature of the Earth's core. Yet, the melt line of Ni at pressures corresponding to the Earth's core has not been explored in the literature. Many previous experimental and simulation efforts have reported the melting point of Ni at pressures below 100 GPa, but there exist large discrepancies, most of which have persisted due to various experimental and simulation bottlenecks in handling extreme pressure and temperature conditions. We adopted the generalized embedded atom method, which overcomes the limitations of existing interatomic potentials, to probe phase stability and phase boundaries of Ni at pressures between 50 and 500 GPa. Further, the potential was validated by comparing the cold curves, phonon dispersion curves, and enthalpies of fusion with ab initio density functional theory calculations. Our analysis shows that face centered cubic (FCC) is stable, and the hexagonal close packed (HCP) and body centered cubic (BCC) phases are metastable close to the melt line. Melting temperatures at different pressures were obtained from two-phase co-existence simulations and take the following functional form: $T$ m = $1969.23+19.15P-0.012P$ 2 . In contrast to iron, differences between the melting points of the stable and metastable phases of Ni are less than 250 K at 300 GPa, and the difference in melting points of the metastable BCC and HCP phases changes sign at 500 GPa, which implies that the phase transition mechanisms during solidification can be very complex.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Distributed strain sensing using Bi-metallic coated fiber optic sensors embedded in stainless steel

Silica fiber optic sensors are resistant to corrosive environments and high temperatures, making them attractive candidates for harsh conditions found in nuclear and aerospace industries. Moreover, fibers can be deployed remotely for continuous measuring of spatially distributed temperatures and strains. Here, this study investigated embedding a Ni/Cu bi-metallic coated fiber in a stainless-steel 316 (SS316) matrix using laser powder bed fusion towards functionalizing metal components for site-specific health monitoring. The embedded fiber was continuously interrogated during controlled heating to 1000°C. The measured fiber strains were similar to the expected differential thermal strains between the fiber and the SS316 matrix, until divergent behavior was observed at temperatures >500°C. No debonding at the matrix–coating–fiber interfaces was observed during microscopy, but significant interactions between the coatings and matrix resulted in diffusion-driven chemistry variations and Kirkendall void formation. Applying the strain-lag theory revealed plastic behavior in the Ni coating at temperatures >500°C, limiting the strain transfer to the fiber at higher temperatures. It was estimated that the elastic modulus in the Ni coating had decreased from ~200 GPa at room temperature to below 40 GPa, starting at 600°C. The low elastic modulus above 600°C is within the margin of what the tangent modulus would be in the case of bilinear isotropic hardening. Regardless of the divergent strain transfer at higher temperatures, the fiber was exposed to the equivalent of 1.9 % engineering strain at 1000°C, but measured only a 0.7 % engineering strain due to the poor strain transfer. Although compensating for the plastic behavior of Ni proved challenging, the bonding of a brittle silica fiber to a metal matrix surviving to 1000°C invites potential iterations on coating material for future application. For example, the embedded fiber is sufficient for acoustic energy transfer, realizing high temperature distributed acoustic sensing.

36 MATERIALS SCIENCE↗
DOE OSTI2021

In situ preparation of visible-light-driven carbon quantum dots/NaBiO 3 hybrid materials for the photoreduction of Cr(VI)

Here, different carbon quantum dots (CQDs)/NaBiO 3 hybrid materials were synthesized as photocatalysts to effectively utilize visible light for the photocatalytic degradation of contaminants effectively. These hybrid materials exhibit an enhanced photocatalytic reduction of hexavalent chromium (Cr(VI)) in the aqueous medium. Zero-dimensional nanoparticles of CQDs were embedded within the two-dimensional NaBiO 3 nanosheets by the hydrothermal process. Compared with that of the pure NaBiO 3 nanosheets, the photocatalytic performance of the hybrid catalysts was significantly high and 6 wt.% CQDs/NaBiO 3 catalyst exhibited better photocatalytic performance. We performed the first-principles density functional theory calculations to study the interfacial properties of pure NaBiO 3 nanosheets and hybrid photocatalysts, and confirmed the CQDs played an important role in the CQDs/NaBiO 3 composites. The experimental results indicated that the enhanced reduction of Cr(VI) was probably due to the high loading of CQDs (electron acceptor) on NaBiO 3 , which made NaBiO 3 nanomaterials to respond in visible light and significantly improved their electron-hole separation efficiency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The Triple Catalytic Action of Tertiary Nitrogen Catalysts in Recyclable Epoxy-Anhydride Thermosets

The thermosetting polymer matrix in fiber reinforced composites is an important component for energy related applications, such as the lightweighting of vehicles or their use in wind and waterpower turbine blades, due to their ability to provide superior adhesion, stiffness, and applicability to a wide range of manufacturing processes. Despite these benefits, today's thermosets are widely considered to be unrecyclable; thus, there is a large interest in redesigning these materials to be inherently recyclable so that energy intensive production of fibers and monomers can be circumvented, bolstering composite manufacture supply chains. Polyester covalent adaptable networks (PECANs) are one such promising alternative to the incumbent, nonrecyclable epoxy-amine thermosets. PECANs can be formed from the ring-opening co-polymerization (ROCOP) of epoxy-anhydride monomer mixtures and subsequent curing at mild temperatures to exhibit similar performance to conventional epoxies while also possessing unique dynamic chemistries along the ester-hydroxyl backbone that are capable of transesterification and thus reprocessability. While significant advancements have been made in formulating these materials for improved mechanical properties or optimizing solvolysis and reprocessing strategies, less attention has been placed on the impact of the residing amine catalyst used to generate the polyester network. In this work, we evaluated the triple-catalytic efficacy of 12 tertiary amines that act as a curing (bulk ROCOP), a transesterification (internal bond exchange), and a deconstruction (methanolysis) catalyst for PECAN thermosets. Specifically, we first distinguish between chain-growth and step-growth polymerization mechanisms for epoxy-amine and epoxy-anhydride mechanisms. We also utilized density functional theory (DFT) to estimate the basicity (pKb) of each catalyst. Of the tested catalysts, the ROCOP of the studied PECAN network can be completed between 95 and 247 min (at 80 degrees C), with variable gelation phenomena. Additionally, the stress relaxation (transesterification metric) efficiency of the tested PECAN networks with alternative embedded catalysts ranged from 95% to 15% reduction in stress after 5 h at 200 degrees C, and the depolymerization efficacy ranged from 2.5% to 9.8% deconstruction after 36 h at 130 degrees C. Overall, the nitrogen-based moieties were demonstrated to influence polymerization kinetics, catalyze the dynamic transesterification exchange mechanism, and aid in the solvolysis of the thermosets at end-of-life.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Symbolic diagnostics to interpret and analyze neural network models

Embedded machine-learned models (EMLMs) have the promise to improve the predictive accuracy of engineering simulators in environments of national interest. EMLMs often comprise complex input-output maps (e.g., neural networks), which make them unamenable to rigorous analysis and generally difficult to interpret. In the face of decades of theory, this lack of interpretability is a significant barrier to building confidence in these models. This work outlines an approach to interpret EMLMs using sparse polynomial regression for comparison with theoretical understanding. To do so, we build on the concept of Locally Interpretable Model-agnostic Explanations (LIME) using physics-informed clustering, prototype selection, and library construction. While general, we demonstrate our method on tensor-basis neural networks used in Reynolds-Averaged Navier-Stokes simulations of hypersonic fluid flows. Results are presented for a simulated toy model and for direct numerical simulations (DNS) of turbulent flows over a flat plate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Enhancing perovskite electrocatalysis through synergistic functionalization of B-site cation for efficient water splitting

The family of perovskite oxides is a promising class of catalysts for diverse energy conversion processes including water splitting. In this work, a facile two-step manipulation (in-situ exsolution and post-sulfurization) strategy was proposed and applied to LaCo 0.2 Fe 0.8 O 3 (LCF) perovskite parent, through which, the electronic state, spatial immersion and intrinsic activity of B-site cobalt (Co) were stepwise tuned at nanoscale proximity accordingly (i.e., lattice Co ions segregated Co 0 → embedded CoS 2 ). Impressively, the as-prepared catalyst (S-LCF) obtains an emergent oxygen deficient microstructure seamlessly pinned with uniformly distributed CoS 2 nanoparticles (NPs), which demonstrates enhanced performance toward both oxygen evolution reaction (OER) and hydrogen evolution reaction (HER), and shows good stability in overall water splitting. The density functional theory (DFT) calculations illustrate the optimized metal-oxygen covalency and hydrogen adsorption Gibbs free energy (Δ G H *) on S-LCF, which further buttresses the prominence of our B-site cation engineering tactics.

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

Design Rules for Open-Shell Molecular Wires: Insights from Correlated Many-Body Transport

Open-shell π-conjugated systems have emerged as promising molecular wires that can sustain unusually high low-bias conductance over tens of nanometers. However, predictive design rules remain limited because near-degeneracy and spin polarization give rise to correlated electronic states that are not reliably captured by standard density functional theory and Landauer transport descriptions. Here, we establish physically transparent design rules for long-range, weakly length-dependent charge transport based on the interplay between bond-length alternation and electron–electron interactions. Using a fully correlated transport framework that combines density matrix renormalization group calculations with nonequilibrium Green’s function embedding, we show that weak dimerization maintains spatially extended edge states, while moderate on-site interaction stabilizes open-shell character without excessively separating transport-relevant resonances. This balance yields zero-bias transmission that is largely insensitive to molecular length, providing a route to connect chemically tunable structure to correlated transport beyond mean-field descriptions.

Charge transport↗
DOE OSTI2026

Residual symmetries and scalar multiplet vacuum alignment in non-Abelian flavour models

We demonstrate that, upon minimizing a renormalizable, single-scalar potential invariant under a non-Abelian symmetry, special orientations in the associated vacuum alignment of the scalar multiplet correspond to the preservation of a discrete residual flavour symmetry in the broken phase of the theory. Conversely, we show that these special scalar alignments are perturbed when additional Lagrangian operators (e.g. renormalizable, multi-flavon operators and/or effective, higher-dimensional operators) are present that break said residual symmetry, leading to a vacuum reorientation and phenomenological consequences. We therefore construct a one-to-one correspondence principle between broken residual symmetries and vacuum alignment corrections, providing a mechanism to identify (and correct) a subtle but persistent form of phenomenologically relevant fine-tuning embedded in — but often ignored by — most successful non-Abelian flavour models. We first establish this correspondence in a set of toy models based on the S4 permutation symmetry, and then apply the lessons learned to the more realistic A4 Altarelli-Feruglio and ∆(27) Universal Texture Zero models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

On the Sensitivity to Density-Functional Approximations for CO Binding Energies of Single-Atom Catalysts in Nitrogen-Doped Graphene

Density functional theory (DFT) methods are the working horse in screening new catalytic materials. They are widely used to predict trends in binding energies, which are then used to compare the activity of different materials. Here, the binding strength of CO is an important descriptor to the CO 2 reduction catalytic activity of the single transition metal atoms embedded on nitrogen-doped graphene (TM/NG). In this work, however, we show that CO binding strengths in different TM/NG has very different sensitivity to DFT methods. Specifically, Fe/NG CO binding energy changes dramatically with the percentage of exact exchange in the functional; Co/NG does less so, while Ni/NG nearly has no change. Such varying behaviors is a direct result of different local spin configurations, similar to the performance of DFT methods for metal porphyrin complexes. Therefore, caution should be exercised when using DFT binding energies for quantitative predictions in TM/NG single atom catalysis.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Hierarchical Model-Free Transactive Control of Building Loads to Support Grid Services

Residential buildings consume 4.4 quads of electricity annually, approximately 37% of the total electricity consumption in the United States. This represents a vast resource that can be used for demand management and other ancillary services. This project aims to develop a robust, scalable hierarchical transactional control mechanism incorporating elements of model-free control (MFC) and game theory to harness buildings to provide ancillary services to the grid. This approach is being taken to address the challenges of incorporating traditional transactional control schemes into existing buildings. The challenges include small individual building sizes requiring aggregation of many buildings, unpredictable energy usage that makes model identification difficult, and satisfying the sensitive occupant comfort constraints. In the proposed approach, by separating the control mechanism into two layers above and below the load aggregator, MFC can be used below the aggregator to modulate flexible building loads in response to pricing signals with guaranteed performance. This allows the burden of identifying an accurate model of the system to be shifted to the above-aggregator layer, where fluctuations in individual building usage have less impact on predicted building system behavior. Game theory concepts can then be used to determine pricing curves and control signals among regional aggregators. Managing this control in a game-theoretic approach will allow us to build in financial incentives that increase customer engagement. Additionally, the usage of MFC necessitates less burdensome computational and communication requirements, thus, it is easily deployable on small, embedded devices. In a broader sense, developing a strategy capable of effectively incorporating residential and small commercial buildings will allow greater throughput of existing and emerging grid services in addition to future transactive energy grid management methods. Using MFC within a hierarchical control architecture will allow the shifting of existing forecasting challenges to an aggregate level, where dynamics are slower and more predictable. This will enable a smooth interface between the grid services requests of utilities and the reliable control required by participating buildings. MFC, which supports distributed control architecture, permits a scalable solution that can be deployed to neighborhood-size systems as well as individual buildings. This project focuses on three objectives: (1) developing the mathematical framework, algorithm toolkit, and software toolset of the two-layer transactive control testbed; (2) developing a scalable solution for application over many residential and small-size commercial buildings with sparse distributed communication; and (3) field testing and implementation on hardware of the control strategies developed in the previous two objectives. The research and development activities are focused and designed to be impactful within the relevant 2025 targets timeframe. An open-source control framework for exploiting variability and dispatchability of building loads will be delivered as the outcome of the project. This capability enables greater participation of loads in electricity markets and ancillary services that are both useful for the utility and financially beneficial for building owners.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Actinide Molten Salts: A Machine-Learning Potential Molecular Dynamics Study

We know that actinide molten salts represent a class of important materials in nuclear energy. Understanding them at a molecular level is critical to proper and optimal design of relevant technological applications. Yet, owing to the complexity of electronic structure due to the 5f orbitals, computational studies of heavy elements in condensed phases using ab initio potentials to study the structure and dynamics of these elements embedded in molten salts are difficult. This lack of efficient computational protocols makes it difficult to obtain information on properties that require extensive statistical sampling like transport. To tackle this problem, we adopted a machine-learning approach to study ThCl 4 -NaCl and UCl 3 -NaCl binary systems. The machine-learning potential, with the density functional theory accuracy, allows us to obtain long molecular dynamics trajectories (ns) for large systems (10 3 atoms) at a considerably low computing cost, thereby efficiently gaining information about their bonding structures, thermodynamics, and dynamics at a range of temperature. We observed a considerable change in the coordination environments of actinide elements and their characteristic coordination-sphere lifetime. Our study also suggests that actinides in molten salts may not follow well known entropy-scaling laws.

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

Dynamical moment fragmentation in the all-in all-out pyrochlore Nd 2 ⁢Sn 2⁢ O 7

We report single crystal neutron spectroscopy and bulk characterization on hydrothermally grown Nd 2 Sn 2 O 7 , revealing magnetic moment fragmentation embedded within the all-in–all-out ordered state. The spectra reveal a nearly flat band with pinch-point momentum dependence accompanied by dispersive branches that produce half-moon features across multiple Brillouin zones. These defining signatures are captured quantitatively by a minimal dipolar-octupolar spin Hamiltonian, demonstrating excellent agreement between experiment and theory. The higher flat-mode energy helps account for the absence of dynamical interference in prior μSR studies, while the lack of any photon-like excitation imposes strict constraints on the proposed Coulombic antiferromagnet scenario. Furthermore, our results extend moment fragmentation to Nd 2 Sn 2 O 7 and identify it as a clean, tractable platform for quantitative exploration of emergent gauge field physics in frustrated magnets.

DC susceptibility measurements↗
DOE OSTI2026

Substrate-Directed Underlayer Growth of Bilayer MoS 2 Revealed by Mo Isotope Labeling

Direct control over the vertical formation sequence and stacking registry in van der Waals (vdW) bilayers is essential for device performance and moiré engineering yet difficult to resolve unambiguously with conventional probes. Here, we use Mo isotope labeling in a two-step chemical vapor deposition process to synthesize bilayer MoS 2 and trace its vertical formation on common substrates. By combining site-selective laser thinning, Raman spectroscopy, time-of-flight secondary ion mass spectrometry, and atomic-resolution scanning transimission electron microscopy (STEM), we find a clear substrate dependence: on SiO 2 /Si, the second layer nucleates and grows beneath the first (underlayer), whereas on sapphire, it forms on top (overlayer). Density functional theory indicates that a larger equilibrium interfacial separation and weaker MoS 2 –substrate interactions on amorphous SiO 2 permit confined interfacial diffusion and underlayer nucleation, whereas stronger interactions and smaller separations on sapphire favor overlayer growth. On SiO 2 , confined epitaxy templates commensurate 2H, 3R, and mixed bilayers, as confirmed by second harmonic generation spectroscopy and STEM. During underlayer coalescence, embedded mirror-twin grain boundaries stitch atomically sharp 2H|3R junctions via alternating 4|8 ring motifs. Molecular-dynamics simulations reveal that these alternating 4|8 motifs accommodate interlayer vdW coupling and locally modulate the stacking registry. These results provide mechanistic insight into confined epitaxial growth and establish isotope labeling as a powerful probe of two-dimensional materials synthesis.

MoS2↗
DOE OSTI2023

Methane Partial Oxidation over Multifunctional 2-D Materials

The objective of this research is to design, synthesize, and evaluate highly selective, active, and stable multifunctional catalysts for the low temperature (< 500 Kelvin (K)) partial oxidation of methane to methanol (MTM) with molecular oxygen: CH 4 (g) + $\frac{1}{2}$O 2 (g) → CH 3 OH(g). Methane, the primary component of natural gas, is a source of energy and economic growth as well as an environmental concern. Recent developments in horizontal drilling and enhanced extraction methods have resulted in production of an estimated 62.4 trillion m 3 of ‘stranded’, or uneconomic, natural gas. Uneconomical natural gas is often flared or vented at remote oil production sites. Leaked, flared, and/or vented gas represents a "lost opportunity”, and this research project aims to maximize the value of the resource. Conventional catalysts for MTM suffer from low methanol selectivity since they exhibit ~0.55 eV higher barrier for C-H bond activation of methane compared to methanol. Without breaking these scaling relations, methanol oxidation is orders of magnitude faster than methane oxidation and it is very challenging to envision a process with economically viable single-pass yield. Here, we chose to investigate single-atom catalysts embedded and stabilized in two-dimensional materials such as graphene (GR) and "supported" on Group VIII and IB transition metals such as nickel. The electronic atomic monolayer-metal support interaction (EAMSI) present in these systems could promote methanol selectivity by breaking the scaling relations of the C-H bond activation of methane and methanol. A density functional theory (DFT) based computational study focused on predicting families of GR-based catalysts that could be active and selective for MTM. The catalyst systems predicted by the computational study were synthesized and evaluated for the gas phase MTM under relevant conditions. Unfortunately, the experimental activity and selectivity was lower than computationally predicted. The origin for the discrepancy is likely related to difficulties in synthesizing single atom catalysts in a threecomponent catalyst system at high density and with high selectivity. Future work in our groups is thus focused on reducing the system complexity to a two-component catalyst system. Finally, a techno-economic analysis (TEA) was also conducted to identify critical bottlenecks that inhibit future commercialization.

03 NATURAL GAS↗
DOE OSTI2026

The presence of Ru metal modifies the behavior of deuterium bound to pyridinic nitrogen in doped graphene-like materials

We used a combination of x-ray photoelectron spectroscopy (XPS) and density functional theory (DFT) to investigate how the presence of a ruthenium (Ru) metal substrate modifies the N-D bond strength of the pyridinic nitrogen species in nitrogen (N)-doped graphene (Gr) on Ru(0001) compared to metal-free highly-oriented pyrolytic graphite (HOPG) substrate. The N-dopants were introduced through low-energy N+/N2+ irradiation. N is embedded in the carbon materials in two predominant configurations: as graphitic N (GN: N substituted in the hexagonal C lattice) and pyridinic N (PN: substitutional N adjacent to a C vacancy). XPS data showed that upon atomic deuterium (D) exposure at 220 K, the PN peak shifted to a higher binding energy by +1.2 eV for HOPG and +1.0 eV for N-doped graphene on Ru(0001), while the GN peak remained unchanged, indicating that the D atoms bound solely to pyridinic N. 85% of PN sites on HOPG can be saturated with D atoms, whereas only ~30% of pyridinic N sites are able to bind D atoms in N-doped Gr/Ru(0001). Our DFT analysis shows that this difference is due to the coordination of PN to Ru atoms, which necessitates bond cleavage of the N-Ru interaction prior to D atom adsorption. D begins desorbing from N-HOPG at ~573 K and is fully desorbed by ~973 K, whereas desorption from N-doped graphene on Ru(0001) begins at ~290 K, with complete desorption observed at ~700 K, indicating that the Ru metal weakens the N-D bond strength. We also studied a high-surface-area, layered, and porous N-doped carbon material, internally labeled NC900, which was synthesized by pyrolysis of glucose and graphitic carbon nitride (g-C3N4) at 1173 K. D exposure caused a +1.1 eV PN peak shift in NC900, with no change in the GN peak. Moreover, the NC900 exhibited the same desorption behavior as HOPG, demonstrating that well-defined model systems can effectively capture the behavior of more complex N-doped carbon materials.

Alupothe Gedara, Buddhika S.↗