Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrogen binding energy”

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

The Surface Chemistry of Methanol on Cu 3 Pd(111): Effects of Metal Alloying and Reaction with Hydrogen

Synchrotron-based ambient-pressure X-ray photoelectron spectroscopy (AP-XPS) was used to study the adsorption and surface chemistry of methanol on a Cu 3 Pd(111) model surface. The composition and morphological properties of the pristine Cu 3 Pd(111) substrate were analyzed using a combination of low-energy electron diffraction (LEED), scanning tunneling microscopy (STM), and low-energy ion scattering (LEIS). The results of ion scattering showed segregation of Pd toward the surface, with a Pd/Cu ratio close to 0.5, a larger value than the ratio of 0.33 expected for a Cu 3 Pd bimetallic structure. The surface of the alloy exhibited a good LEED pattern with three-fold long-range periodicities. In STM, clusters of palladium with hexagonal arrays and Pd- Pd distances of 2.7-2.8 Å were detected. Bonding to copper perturbed the electronic properties of the atoms in the Pd clusters, shifting their 4d states toward higher binding energy with respect to the Fermi level. Further, the valence band spectrum of Cu 3 Pd(111) exhibited a line shape that was very different from those displayed by Cu(111) or Pd(111). At low pressures, the adsorption of methanol on Cu 3 Pd(111) at 300 K mainly produced CH 3 O, CO and CH x species. AP-XPS showed that most Pd atoms in the surface of the Cu 3 Pd(111) alloy interacted with the decomposition products of methanol. No significant changes were observed in the core levels of copper upon the adsorption and dissociation of methanol, suggesting that the molecule mainly interacted with Pd sites of the alloy. Reaction with hydrogen led to fast removal of CH x , C and PdC x species from Cu 3 Pd(111) at moderate (< 450 K) temperatures and prevented a CH x → C transformation. If the stability of adsorbed CH 3 O is used as a descriptor for the hydrogenation of CO 2 to methanol, Cu 3 Pd should be a much better catalyst than monometallic palladium. This may be a consequence of electronic and ensemble effects in the alloy that moderate the reactivity of Pd sites.

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

Critical role of solvent-modulated hydrogen-binding strength in the catalytic hydrogenation of benzaldehyde on palladium

Protic and non-protic solvents influence the rate of benzaldehyde hydrogenation on Pd by one order of magnitude in the following sequence of increasing reactivity, dioxane < tetrahydrofuran < water < methanol. Despite these differences, the reaction pathway does not change; the majority of turnovers occurs via a Langmuir-Hinshelwood mechanism by stepwise addition of sorbed H to sorbed benzaldehyde, first to the carbonyl O and then to the C atom of the formyl group, forming benzyl alcohol. An analysis of the solvation energies shows that both ground and transition states are destabilized by the solvents compared to the same reaction at the gas-solid interface. The destabilization extent of the reacting organic substrates in both states are similar and, therefore, compensate each other, making the net kinetic effects inconsequential. In consequence, the marked reactivity differences arise solely from the differences in the surface stabilization of adsorbed hydrogen adatoms moderated by the four solvents-more weakly bound hydrogen adatoms, because of their solvation by solvent, exhibit higher rates of their hydrogenation reaction in kinetically relevant attack on the C atom.

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

The crystal structure of Vibrio cholerae (6-4) photolyase reveals interactions with cofactors and a DNA-binding region

Photolyases (PLs) reverse UV-induced DNA damage using blue light as an energy source. Of these PLs, (6-4) PLs repair (6-4)-lesioned photoproducts. We recently identified a gene from Vibrio cholerae (Vc) encoding a (6-4) PL, but structural characterization is needed to elucidate specific interactions with the chromophore cofactors. Here, we determined the crystal structure of Vc (6-4) PL at 2.5 Å resolution. Our high-resolution structure revealed that the two well-known cofactors, flavin adenine dinucleotide and the photoantenna 6,7-dimethyl 8-ribityl-lumazin (DMRL), stably interact with an α-helical and an α/β domain, respectively. Additionally, the structure has a third cofactor with distinct electron clouds corresponding to a [4Fe-4S] cluster. Moreover, we identified that Asp106 makes a hydrogen bond with water and DMRL, which indicates further stabilization of the photoantenna DMRL within Vc (6-4) PL. Further analysis of the Vc (6-4) PL structure revealed a possible region responsible for DNA binding. The region located between residues 478 to 484 may bind the lesioned DNA, with Arg483 potentially forming a salt bridge with DNA to stabilize further the interaction of Vc (6-4) PL with its substrate. Our comparative analysis revealed that the DNA lesion could not bind to the Vc (6-4) PL in a similar fashion to the Drosophila melanogaster (Dm, (6-4)) PL without a significant conformational change of the protein. The 23rd helix of the bacterial (6-4) PLs seems to have remarkable plasticity, and conformational changes facilitate DNA binding. In conclusion, our structure provides further insight into DNA repair by a (6-4) PL containing three cofactors.

(6-4) lesioned photoproduct↗
DOE OSTI2022

Dynamics of deuterium retention and desorption from plasma-facing materials in fusion reactor-relevant conditions

Hydrogen isotopes retention and desorption during and after discharges in fusion devices are still not well understood due to the complex device conditions and limitations of in-situ diagnostics and measurements. Here, we simulated well-diagnosed recent experiments at the DIII-D facility to benchmark our ITMC-DYN integrated package of modeling deuterium diffusion, retention, and desorption during and after D discharge irradiation. Modeling results were compared with detail experimental data of D desorption fluxes for various irradiation conditions. We predicted the temporal evolution of free and trapped D distribution in tungsten (W) plasma-facing material (PFM). Effects of key parameters namely diffusion coefficient, recombination rate, trapping energies against different defect types, were examined in these simulations. Existing experimental data of these parameters in literature varies significantly which makes it harder to identify key mechanisms and physics responsible for hydrogen isotope retention and desorption. The purpose of this work is to accurately simulate recent well-diagnosed reactor experiments given the uncertainties in such parameters and identify mechanisms responsible for the retention and desorption. We implemented the best identified diffusion, recombination, and trapping parameters in ITMC-DYN package that integrate both various collisional and thermal processes. We predicted, for example, that sample cooling between discharges in DIII-D operations can significantly affect the spatial distribution of trapped D in W under reactor irradiation conditions. Correct prediction of desorption spectra from samples irradiated during 10 DIII-D discharges showed that up to 35% of D can be retained in high binding energy defects such as vacancy clusters or voids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Dual Enhancement of Thermostability and Activity of Xylanase through Computer-Aided Rational Design

In the realm of enzyme engineering, the dual enhancement of thermostability and activity remains a challenge. Herein, we employed a computer-aided approach integrating folding free energy calculations and evolutionary analysis to engineer Paecilomyces thermophila xylanase into a hyperthermophilic enzyme for application in the paper and pulp industry. Through the computational rational design, XynM9 with superior thermostability and enhanced activity was designed. Its optimal reaction temperature increases by 10 °C to 85 °C, its T m increases by 10 °C to 93 °C, and its half-life increases 11-fold to 5.8 h. Additionally, its catalytic efficiency improves by 57% to 3926 s –1 mM –1 . Molecular dynamics simulations revealed that XynM9 is stabilized by more hydrogen bonds and salt bridges than wild-type xylanase. The mutant’s narrower catalytic cleft enhances the substrate-binding affinity, thus improving the catalytic efficiency. In harsh conditions at 80 °C and pH 10, using XynM9 significantly reduced both hemicellulose and lignin, which makes it a good candidate for use in the paper and pulp process. Finally, our study presents an accurate and efficient strategy for the dual enhancement of enzyme properties, guiding further improvement of computational tools for protein stabilization.

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

Cobalt(II) Phthalocyanine Substituents Tune the Electrocatalytic CO 2 Conversion to Methanol

Cobalt phthalocyanine (Co(II)Pc) and its derivatives are promising molecular electrocatalysts for the electrochemical reduction of CO 2 to CO and methanol (CH 3 OH). Despite increasing interest, a detailed mechanistic understanding of how ligand substituents influence catalytic activity, selectivity, and efficiency remains limited. In this study, we employ density functional theory (DFT) to systematically investigate the influence of electron-donating groups (EDGs) and electron-withdrawing groups (EWGs) on the electronic structure and redox properties of the Co(II)Pc electrocatalyst, and to elucidate CO 2 RR mechanistic pathways. Our results reveal that EWGs cause a positive shift in the reduction potentials, favor CO 2 binding over protonation of the Co metal center and promote downstream methanol formation at mild potentials. EDGs show opposite trends including favorable protonation steps, promoting a negative shift in the reduction potential, and facilitating the hydrogen evolution reaction (HER), which competes with the desired CO 2 RR pathway. Notably, CO dissociation is thermodynamically and kinetically unfavorable across all systems, positioning the redox potential versus CO dissociation energy as a key factor for methanol selectivity. Furthermore, these insights provide a predictive framework for rational catalyst design and underscore the critical role of electronic tuning in advancing molecular electrocatalysts for sustainable CO 2 conversion.

Alcohols↗
DOE OSTI2026

Modulating Operational Conditions to Mitigate Deactivation in Formate Dehydrogenation on Pd Phases

In heterogeneous catalysis, poisoning by surface-bound intermediates poses a major barrier to sustained catalyst performance in (de)hydrogenation reactions. Formate/bicarbonate systems, as liquid organic hydrogen carriers (LOHCs), offer a CO 2 -integrated, low-temperature pathway for hydrogen storage and release, making them attractive for circular energy applications. However, their lower hydrogen density and susceptibility to catalyst deactivation limit their competitiveness compared to conventional LOHCs like methylcyclohexane. Here, this study investigates the mechanistic origins of formate (HCOO – ) dehydrogenation and associated deactivation on Pd interfaces. Using density functional theory (DFT) simulations, we show that under thermocatalytic conditions, strongly bound formate accumulates on the catalyst surface (Pd(111)), blocking active sites, raising activation barriers, and leading to progressive performance loss. Because formate adsorption involves charge transfer, we exploit its sensitivity to electronic structure by modulating the electrochemical potential of the catalyst. Our results reveal that hydrogen transfer from water and formate exhibits opposing potential dependencies, providing insights into the opposing driving forces behind both catalytic activity and poisoning. To further probe this phenomenon, we examine electrochemically induced phase transitions in Pd, focusing on PdO(100) and PdH(110), which are stable under oxidizing and reducing potentials, respectively, and demonstrate enhanced dehydrogenation activity between −0.4 and 0.2 V vs standard hydrogen electrode (SHE). Complementary thermal treatments help decouple kinetic and thermodynamic contributions to intermediate binding. These findings underscore the critical role of the catalyst phase and external stimuli in dictating poison-active site interactions and highlight phase engineering as a promising strategy to mitigate deactivation. This work offers mechanistic insights and design principles for developing more resilient and efficient catalysts for LOHC applications under realistic operating conditions.

Pd phase↗
DOE OSTI2020

Manipulating the Geometric and Electronic Structures of Manganese Nitrides for Ammonia Synthesis

Abstract Manganese (Mn) nitrides are important nitrogen (N) carriers for small‐scale intermittent ammonia (NH 3 ) synthesis. However, only 3∼8 % of lattice N are converted into NH 3 . In this study, the geometric and electronic structures of well‐defined Mn 4 N and Mn 2 N lattices were altered using transition metal heteroatoms (Cr, Fe, Co, Ni, Mo) to understand the driving force behind lattice N diffusion and extraction. Density Functional Theory (DFT) revealed that the binding of early hydrogenation product (NH) follows a linear relationship with lattice N over a wide range of close‐packed surfaces. But, the binding of NH 2 and NH 3 are more sensitive to the geometric and electronic structures. Further, the chemical bonding can be quantitatively characterized with the covalency derived from Crystal Orbital Hamiltonian Population (COHP). In the Eley Rideal‐Mars van Krevelan pathway, the overall NH 3 formation free energy ( ) and lattice N diffusion barrier ( ) are the respective thermodynamic and kinetic determining factors. Aided by a rate‐determining step (RDS) model, Mn 4 N modified with Fe, Co, Ni single‐atom dopants all show enhanced rates for NH 3 formation.

Shan, Nannan↗
DOE OSTI2022

Hydrogen trapping energetics at BCC iron-helium interfaces

Density functional theory (DFT) calculations have been performed to assess the trapping and segregation strength of hydrogen (H) to helium (He) interfaces in BCC (Fe) as a function of He density and surface orientation. The He density ranges from 0 to 2 He/V which corresponds to equilibrium bubbles that are expected to form in the structural components in a fusion reactor environment. The DFT calculations consist of a slab of BCC metal and an initial lattice of FCC He. Here, the report on the binding energy of H to these interfaces as well as the migration barriers into and away from the surface, which together provide information on the de-trapping energy for H at He bubbles. It was found that the binding energy of H to He-Fe interfaces decreases with increasing He density. Migration barriers are also modified due to the He density, which are sensitive to the surface orientation as well as the diffusion pathway. These results provide valuable first-principles energetics that are necessary for higher-order, mesoscale models, and provide insight into the extent to which He bubbles may trap tritium in fusion structural materials.

42 ENGINEERING↗
DOE OSTI2026

Ensemble Effects on Hydroxide Bond Dissociation Free Energies in Polyoxovanadate Clusters

Understanding structure-property relationships is foundational to numerous modern chemistries, such as proton-coupled electron transfer (PCET). However, an experimentally measured property is the result of the behavior from an ensemble of molecules. Neglecting ensemble effects, especially under complex chemical environments, may obfuscate these relationships and lead to discrepancies between theory and experiment. In this work, we demonstrate the impact of configurational entropy and local chemical environments on hydroxide bond dissociation free energies [BDFE- (O−H)] for a set of polyoxovanadate nanoclusters, at ambient conditions. The O−H bond strengths are investigated via density functional theory (DFT) coupled with statistical thermodynamic analysis and bilinear modeling, and compared with previous experimental results on the same systems, namely electrochemical solutions of: [V 6 O 13−x (OH) x (TRIOL R ) 2 ] −2 (x = 2, 4, 6; R = NO 2 , Me) and [V 6 O 11−x (OMe) 2 (OH) x (TRIOL NO 2 ) 2 ] −2 (x = 2, 4). Interestingly, we find that ensemble effects, even at room temperature, can account for a significant portion of the BDFE(O−H) trend with the degree of reduction via H atom binding, which cannot be fully captured by single-structure, static DFT calculations. Moreover, we find that the ensemble effects may be replicated statistically, requiring only enumeration of energetically accessible H-binding sites. With the ensemble effects resolved, we present a simple bilinear model to reconcile remaining biases between experiment and ensemble-informed theory, which corelate with clusterspecific electronic environment differences. The bilinear model achieves outstanding accuracy vs experiments with a root-mean squared error of 0.4 kcal/mol. Finally, based on the physicochemical characteristics of hydrogen interaction with polyoxometalates, we present a simple methodology that captures the BDFE(O−H) trend while dramatically reducing required DFT calculations by 98% and achieving accuracy within 1 kcal/mol. Overall, this work elucidates the roles and structural origins of configurational entropy and chemical effects on polyoxometalate hydroxide bond energies, with potential applicability to various atomically precise metal oxide systems. Importantly, it introduces models for rapid and highly accurate property calculations in connection with experiments.

Cluster chemistry↗
DOE OSTI2020

A self-adjusting platinum surface for acetone hydrogenation

We show that platinum displays a self-adjusting surface that is active for the hydrogenation of acetone over a wide range of reaction conditions. Reaction kinetics measurements under steady-state and transient conditions at temperatures near 350 K, electronic structure calculations employing density-functional theory, and microkinetic modeling were employed to study this behavior over supported platinum catalysts. The importance of surface coverage effects was highlighted by evaluating the transient response of isopropanol formation following either removal of the reactant ketone from the feed, or its substitution with a similarly structured species. The extent to which adsorbed intermediates that lead to the formation of isopropanol were removed from the catalytic surface was observed to be higher following ketone substitution in comparison to its removal, indicating that surface species leading to isopropanol become more strongly adsorbed on the surface as the coverage decreases during the desorption experiment. This phenomenon occurs as a result of adsorbate–adsorbate repulsive interactions on the catalyst surface which adjust with respect to the reaction conditions. Reaction kinetics parameters obtained experimentally were in agreement with those predicted by microkinetic modeling when the binding energies, activation energies, and entropies of adsorbed species and transition states were expressed as a function of surface coverage of the most abundant surface intermediate (MASI, C 3 H 6 OH*). It is important that these effects of surface coverage be incorporated dynamically in the microkinetic model (e.g., using the Bragg–Williams approximation) to describe the experimental data over a wide range of acetone partial pressures.

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

Molecular dynamics simulations of reflection and sputtering behavior of boron under deuterium ion irradiation

Boronization is a commonly used method of wall conditioning in fusion reactors. The application of boron films to the plasma-facing materials results in enhanced plasma performance due to the reduction of intrinsic impurities. This is primarily driven by a reduction in oxygen content that is chemically trapped in the boron film. The reactive nature of these boron films also raises questions concerning interactions with hydrogen isotopes. In this work, boron-deuterium interactions were studied using molecular dynamics (MD). Reactive force field potentials were used to model the chemical interactions between B and D. An amorphous boron substrate was irradiated by D atoms at varying incident energies, 10 eV < E i < 150 eV and angles, 0° < α < 85°. The reflection probability was calculated and compared to results from the commonly-used binary collision approximation (BCA) method. This comparison found that the BCA underestimated the reflection probability at E i < 35 eV and α > 45 ∘ . The source of this discrepancy was found to be the surface binding energy model. The BCA calculation with an isotropic surface binding energy model was more closely aligned to the MD result. This, in combination with a correction function based on the MD results allows for corrections to the reflection probability of deuterium impinging on boron surfaces. The sputtering of the substrate material was also studied. While this study did not contain sufficient events to quantitatively describe the sputtering behavior, some qualitative results emerged: namely, chemical sputtering of B and D-containing molecules (BD, BD 2 , BD 3 ) at low ( < 20 eV) incident deuterium energies. This result suggests that chemical sputtering could be a significant factor in limiting boron coating lifetime when exposed to lower ion energies, such as those in detached plasmas. The results show that chemical interactions should be taken into account when modeling the interactions between D ions and B surfaces.

Boron↗
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↗
DOE OSTI2023

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the “gold standard” coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller−Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

Chemistry↗
DOE OSTI2021

Structural and Catalytic Properties of Isolated Pt 2+ Sites in Platinum Phosphide (PtP 2 )

This article describes the synthesis and catalytic properties of supported, 2-3 nm platinum phosphide (PtP 2 ) nanoparticles (NPs). Depending on the P loading, two PtP 2 structures are formed, that is, a PtP 2 surface on a (metallic) Pt core (Pt@PtP 2 ) and single-phase PtP 2 NPs. The structures were determined using extended X-ray absorption fine structure , in situ synchrotron X-ray diffraction, and scanning transmission electron microscopy. In PtP 2 NPs, Pt 2+ ions are geometrically isolated by P 2 2- ions, at a Pt-Pt distance of 4.02 Å, which is much longer than 2.78 Å in (metallic) Pt NPs. The oxidation state of Pt in PtP 2 NPs was determined by in situ X-ray absorption near-edge structure and in situ X-ray photoelectron spectroscopy and was found to be consistent with Pt 2+ ions even after treatment in H 2 at 550 °C. Unlike Pt NPs, which are highly active for propylene hydrogenation at room temperature, PtP 2 NPs are not active below about 150 °C, suggesting the absence of metallic surface Pt. In contrast to metallic Pt, which is poorly selective for acetylene hydrogenation, PtP 2 NPs display high selectivity toward ethylene. PtP 2 , also has high olefin selectivity for propane dehydrogenation, although the rate per g Pt is about 7 times lower than that of metallic Pt NPs of the same size. In situ resonant inelastic X-ray scattering spectroscopy shows that the energy of the filled Pt 5d valence orbitals is 1.5 eV lower than that of metallic Pt, which leads to weaker adsorbate binding consistent with its catalytic properties. Here, a H 2 -stable Pt 2+ site suggests different catalytic applications for these catalysts as compared to Pt NPs.

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

Capturing CO 2 in Quadrupolar Binding Pockets: Broadband Microwave Spectroscopy of Pyrimidine-(CO 2 ) n , n = 1,2

Pyrimidine has two in-plane CH(δ+)/N̈(δ–)/CH(δ+) binding sites that are complementary to the (δ–/2δ+/δ–) quadrupole moment of CO 2 . For this study, we recorded broadband microwave spectra over the 7.5–17.5 GHz range for pyrimidine-(CO 2 ) n with n = 1 and 2 formed in a supersonic expansion. Based on fits of the rotational transitions, including nuclear hyperfine splitting due to the two 14 N nuclei, we have assigned 313 hyperfine components across 105 rotational transitions for the n = 1 complex and 208 hyperfine components across 105 rotational transitions for the n = 2 complex. The pyrimidine-CO 2 complex is planar, with CO 2 occupying one of the quadrupolar binding sites, forming a structure in which the CO 2 is stabilized in the plane by interactions with the C–H hydrogens adjacent to the nitrogen atom. This structure is closely analogous to that of the pyridine-CO 2 complex studied previously by (Doran, J. L. J. Mol. Struct. 2012, 1019, 191–195). The fit to the n = 2 cluster gives rotational constants consistent with a planar cluster of C 2v symmetry in which the second CO 2 molecule binds in the second quadrupolar binding pocket on the opposite side of the ring. The calculated total binding energy in pyrimidine-CO 2 is –13.7 kJ mol –1 , including corrections for basis set superposition error and zero-point energy, at the CCSD(T)/ 6-311++G(3df,2p) level, while that in pyrimidine-(CO 2 ) 2 is almost exactly double that size, indicating little interaction between the two CO 2 molecules in the two binding sites. The enthalpy, entropy, and free energy of binding are also calculated at 300 K within the harmonic oscillator/rigid-rotor model. This model is shown to lack quantitative accuracy when it is applied to the formation of weakly bound complexes.

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

Multi-Scale Modeling for Plasma-Enhanced Ammonia Decomposition over Carbides and Nitrides

Ammonia is a carbon-free hydrogen carrier, but its decomposition typically requires high temperatures over costly Ru-based catalysts due to the large barrier for N≡N bond formation. We develop a multiscale framework combining density functional theory, zero-dimensional plasma kinetics, and microkinetic modeling to elucidate how non-thermal plasma (NTP) enables low-temperature NH 3 decomposition over Co-based carbides and nitrides, benchmarked against Ru and Co. Under thermal conditions, all catalysts are limited by N≡N bond formation, with Co 3 C(001) most active owing to its negatively charged surface, strong N* binding, and low activation barriers of N≡N bond formation. Plasma-induced vibrational excitation of NH 3 and its reactive radicals promotes a radical-driven •NH 2 –N* coupling pathway that dominates on Co 3 C(001) and Co 3 N(001), shifting the rate-limiting step to NH 3 (v1) dissociation, increasing turnover frequencies by up to 6 orders of magnitude, and reducing the temperature needed to reach a turnover frequency of 5 s –1 from >680 °C (Ru and Co under thermal condition) to 267 °C (Co 3 C) and 415 °C (Co 3 N). These results identify Co-based carbides and nitrides as promising plasma-active catalysts for energy-efficient hydrogen production from ammonia.

ammonia decomposition↗