Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “binding energy difference”

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

Probing surface-adsorbate interactions through active particle dynamics

Adsorbate molecules present in a reaction mixture may bind to and block catalytic sites. Measurement of the surface coverage of these molecules via adsorption isotherms is critical for modeling and design of catalytic reactions on surfaces. However, it is challenging to measure isotherms in solution in a way that is directly relevant to catalytic activity under reaction conditions, particularly since adsorbates may bind with an enormous range of surface affinity parameters. Here we used the motion of self-propelled catalytic Janus particles, which employ the decomposition of hydrogen peroxide fuel as a propulsion mechanism, to determine the effective surface coverage of thioglycerol, furfural, and ethanol on a platinum surface as a function of concentration in aqueous solution by measuring the decrease in active motion due to the blocking of active sites. For strongly adsorbing thioglycerol, this effective coverage was compared and contrasted to the total adsorbed amount measured using inductively-coupled plasma analysis. Demonstrating the broad applicability of this approach, the surface affinity of the three adsorbates spanned more than four orders of magnitude. For each species, the adsorbate-mediated attenuation of active motion occurred over a wide concentration range and was well-described by a Langmuir isotherm. The strongly interacting thioglycerol had the highest affinity towards the surface (K a = 15.5 ± 4.3 mM –1 ) and fully deactivated the active particle motion at surface saturation. Furfural had an intermediate affinity (K a = 0.42 ± 0.07 mM –1 ) but did not fully block H 2 O 2 access to the surface at apparent saturation, consistent with a maximum fractional surface coverage of θ max = 0.67. Ethanol exhibited even lower affinity (K a = 0.0025 ± 2x10 -4 mM –1 ) and its coverage saturated at only θ max = 0.38. Analysis of isotherms at elevated temperatures enabled direct extraction of the enthalpies of adsorption. The degree of surface coverage at adsorbate saturation appeared to correlate with the relative energies of adsorption for the different adsorbate species and was consistent with adsorbate saturation of one of multiple active site populations towards H 2 O 2 decomposition. Furthermore, computational investigations into solvent effects on furfural adsorption showed good quantitative agreement with the experimental results. This work leverages unique properties of active particles to explore fundamental catalysis questions and demonstrates a novel paradigm for significant and experimentally accessible multidisciplinary research.

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

The carboxylation status of osteocalcin has important consequences for its structure and dynamics

The carboxylation status of Osteocalcin (Ocn) not only influences formation and structure in bones but also has important endocrine functions affecting energy metabolism and expenditure. In this study, the role of γ-carboxylation of the glutamate residues in the structure-dynamics-function relationship in Ocn is investigated. Three forms of Ocn, differentially carboxylated at the Glu-17, 21 and 24 residues, along with a mutated form of Ocn carrying Glu/Ala mutations, are modeled and simulated using molecular dynamics (MD) simulation in the presence of calcium ions. Characterization of the global conformational dynamics of Ocn, described in terms of the orientational variations within its 3-helical domain, highlights large structural variations in the non-carboxylated osteocalcin (nOcn). The bi-carboxylated Ocn (bOcn) and tri-carboxylated (tOcn) species, in contrast, display relatively rigid tertiary structures, with the dynamics of most regions strongly correlated. Radial distribution functions calculated for both bOcn and tOcn show long-range ordering of the calcium ion distribution around the carboxylated glutamate (γGlu) residues, likely playing an important role in promoting stability of these Ocns. Additionally, the same calcium ions are observed to coordinate with neighboring γGlu, better shielding their negative charges and in turn stabilizing these systems more than do the singly coordinating calcium ions observed in the case of nOcn. bOcn is also found to exhibit a more helical C-terminal structure, that has been shown to activate its cellular receptor GPRC6A, highlighting the allosteric role of Ocn carboxylation in modulating the stability and binding potential of the active C-terminal. Furthermore, the carboxylation status of Ocn as well and its calcium coordination appear to have a direct influence on Ocn structure and dynamics, possibly leading to the known differences in Ocn biological function.

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

Formation of Gas-Phase Allyl Radicals from Glycerol on Rutile TiO2(110)

The reaction pathways of glycerol on partially reduced TiO2(110) have been studied by a combination of molecular beam dosing and temperature-programmed desorption (TPD). Glycerol is found to bind more strongly to TiO2(110) than glycols due to an additional hydroxyl group in the molecule. When the Ti surface sites are saturated by glycerol, the majority of glycerol molecules (~ 95%) undergo further reactions to yield carbon-free (water and hydrogen) and carbon-containing products that are desorbed during the TPD ramp. In contrast to glycols for which both water and hydrogen arise exclusively from the cleavage of glycol O-H bonds, we identify two distinct reaction channels resulting from the cleavage of O-H and C-H bonds, respectively, for both water and hydrogen evolution from glycerol on TiO2(110). Quantitative determination of the desorption yields of water and hydrogen reveals that the C-H bond scission channel which occurs at higher temperature than the O-H bond scission channel, is a minor reaction channel for water formation, while it is a major reaction channel responsible for hydrogen evolution. For the carbon-containing products, we exclude both carbon deposition and carbon oxides (CO and CO2) formation, and identify the allyl radical as the major carbon-containing product. The formation of the allyl radical is accompanied by the formation of various minor carbon-containing products including propylene, acrolein, allyl alcohol, acetone, formaldehyde and ethylene. The different surface chemistries observed between glycerol and simpler alcohols and glycols on TiO2(110), suggest that smaller oxygenates cannot serve as models for larger oxygenates in probing their reaction pathways on oxide surfaces. The new reaction channels observed for glycerol on TiO2(110) may open up a new way to convert bio-renewable glycerol directly to pure hydrogen and hydrocarbons on TiO2. This work was supported by the US Department of Energy, Office of Sciences, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences, and performed in EMSL, a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory (PNNL). PNNL is a multiprogram national laboratory operated for the DOE by Battelle.

Chen, Long↗
DOE OSTI2020

Solubility-Dependent Protective Effects of Binary Alloys for Lithium Anode

Recently, alloys have been widely utilized as protection layers for homogenized Li growth. However, the protective mechanisms for different alloys are not clearly understood. In this work, guided by the binary phase diagrams, Mg, Au, Zn, Al, Fe, and Cu protective layers have been selected for study of their protective abilities and mechanisms. The selected metals can be classified into three categories according to the nature of their associated alloy protective layers: (1) metals with high solubility in Li while not forming any intermediate line compounds, e.g., Mg, which is good for fast Li stripping; (2) metals with limited solid solubility but forming stable line compounds with Li upon the lithiation, e.g., Au, Zn, and Al, where the formed stable line compounds Li x Au, Li x Zn, and Li x Al will hamper Li kinetics; (3) metals with negligible mutual solubility and no line compound formation with Li, e.g., Fe and Cu, which block the Li stripping path and therefore result in the highest polarization voltage.

25 ENERGY STORAGE↗
DOE OSTI2022

Acid anion electrolyte effects on platinum for oxygen and hydrogen electrocatalysis

Platinum is an important material with applications in oxygen and hydrogen electrocatalysis. To better understand how its activity can be modulated through electrolyte effects in the double layer microenvironment, herein we investigate the effects of different acid anions on platinum for the oxygen reduction/evolution reaction (ORR/OER) and hydrogen evolution/oxidation reaction (HER/HOR) in pH 1 electrolytes. Experimentally, we see the ORR activity trend of HClO 4 > HNO 3 > H 2 SO 4 , and the OER activity trend of HClO 4 > HNO 3 ~ H 2 SO 4 . HER/HOR performance is similar across all three electrolytes. Notably, we demonstrate that ORR performance can be improved 4-fold in nitric acid compared to in sulfuric acid. Assessing the potential-dependent role of relative anion competitive adsorption with density functional theory, we calculate unfavorable adsorption on Pt(111) for all the anions at HER/HOR conditions while under ORR/OER conditions ClO$^{-}_{4}$ binds the weakest followed by NO$^{-}_{3}$ and SO$^{2-}_{4}$. Our combined experimental-theoretical work highlights the importance of understanding the role of anions across a large potential range and reveals nitrate-like electrolyte microenvironments as interesting possible sulfonate alternatives to mitigate the catalyst poisoning effects of polymer membranes/ionomers in electrochemical systems. These findings help inform rational design approaches to further enhance catalyst activity via microenvironment engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Satellites, core hole excitations, and spin-resolved electronic structure in the spectroscopy of half-metallic CrO 2

Photoelectron satellites—the structures appearing on the low kinetic or high binding-energy side of the “main” or “elastic” photopeak—betray the complex many-body interactions set in motion by the sudden creation of the core hole. In this work, we demonstrate, using the technologically important ferromagnetic half-metal CrO 2 , how such satellites can manifest themselves in other core-level spectroscopies of the material and how they can reveal important details pertinent to its electronic structure. Specifically, we identify a fluorescence satellite in the Cr 𝐿 3 resonant x-ray-emission spectra that radiates at a constant emission energy across the Cr 𝐿 3 x-ray edge with energy ≈1.3 eV above the ordinary valence fluorescence. Here, we provide evidence that this feature arises from the valence recombination of the Cr 2⁢𝑝 core hole “dressed” by the same shakeup charge-transfer process present in both the Cr x-ray photoelectron and the Cr x-ray absorption spectra with its energy uniquely measuring the exchange splitting of the Cr 3⁢𝑑 level. Further analysis of the x-ray emission data reveals three additional features that radiate at constant loss energy that are attributed to combinations of Cr 3⁢𝑑⁢(𝑡 2⁢𝑔 ) → Cr 3⁢𝑑⁢(𝑡 2⁢𝑔 ), charge-transfer O 2⁢𝑝→Cr 3⁢𝑑, and crystal-field Cr⁢ 3⁢𝑑⁡(𝑡 2⁢𝑔 )→Cr⁢ 3⁢𝑑⁡(𝑒 𝑔 ) excitations. These assignments and their energies are supported by density-functional theory calculations, the accuracy of which we demonstrate by hard x-ray valence-photoemission measurements. Atomic multiplet calculations, which include crystal-field effects, help interpret x-ray photoelectron and x-ray absorption spectra of the covalently mixed Cr ion. Resonant Cr K-𝐿 2,3 ⁢𝐿 2,3 Auger-electron emission spectra support a ligand-to-metal nature of the charge-transfer process while highlighting the charge sensitivity differences between photon-in/electron-out and photon-in/photon-out spectroscopies.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab initio benchmarks of neutrinoless double- β decay in light nuclei with a chiral Hamiltonian

Here, we report ab initio benchmark calculations of nuclear matrix elements (NMEs) for neutrinoless double-β (0νββ) decays in light nuclei with mass numbers ranging from A = 6 to A = 22. We use the transition operator derived from light-Majorana neutrino exchange and evaluate the NME with three different methods: two variants of in-medium similarity renormalization group (IMSRG) and importance-truncated no-core shell model (IT-NCSM). The same two plus three nucleon interaction from chiral effective field theory is employed, and both isospin-conserving (ΔT = 0) and isospin-changing (ΔT = 2) transitions are studied. We compare our resulting ground-state energies and NMEs to those of recent ab initio no-core shell model and coupled-cluster calculations, also with the same inputs. We show that the NMEs of ΔT = 0 transitions are in good agreement among all calculations, at the level of 10%. For ΔT = 2, relative deviations are more significant in some nuclei. The comparison with the exact IT-NCSM result allows us to analyze these cases in detail, and indicates the next steps toward improving the IMSRG-based approaches. The present study clearly demonstrates the power of consistent cross checks that are made possible by ab initio methodology. This capability is crucial for providing meaningful many-body uncertainties in the NMEs for the 0νββ decays in heavier candidate nuclei, where quasiexact benchmarks are not available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Star Block Copolymers at Homopolymer Interfaces: Conformation and Compatibilization

Star block copolymers (s-BCPs), comprised of multiple linear diblock copolymers joined at a central point, are shown to segregate to the interface between two immiscible homopolymers that are identical to the blocks of the s-BCPs. The s-BCPs undergo a configurational transition at the interface, with different blocks of copolymers being embedded in their respective homopolymers, thereby bridging the interface and promoting adhesion. A series of 4-arm s-BCPs were synthesized with hydrogenated or deuterated polystyrene (PS/dPS) as the core block and poly(2-vinylpyridine) (P2VP) as the corona block, which was directly placed at the interface between the two homopolymers. Neutron reflectivity (NR) was used to determine the concentration profiles of the PS homopolymer, s-BCP core blocks, and P2VP total segments under equilibrium. The investigation varies the molecular weight (MW) and the total number of s-BCPs at the interface. Self-consistent-field theory (SCFT) was also employed to calculate the concentration profiles of the components at the interface, which were in excellent agreement with experimental results. The NR showed that the interfacial width between the homopolymers increased with the increasing number of s-BCPs at the interface up to a saturation limit. Beyond this limit, additional s-BCPs were released into the corona-miscible phase as unimolecular micelles. For a comparable interlayer thickness of s-BCPs at the interface, lower MW s-BCPs generated a broader interface. SCFT analysis suggested that, at the same packing density, the arms of the low MW s-BCPs align more parallel to the interface, while the arms of high MW s-BCPs adopt a more normal orientation, like their linear BCP counterparts. Furthermore, it was also observed that the core blocks, constrained by the junction points, were oriented more parallel and closer to the interface than the corona blocks. The phase behavior of the polymer blends revealed that s-BCP additives can efficiently reduce the domain size, with the low MW yielding smaller domain sizes due to the greater reduction in the interfacial energy and the high MW arresting phase separation due to their higher binding energy and a jamming of the interfacial assemblies. Asymmetric double cantilever beam (ADCB) tests demonstrated that s-BCPs promoted adhesion more efficiently than their linear BCP counterparts due to stronger binding energy per molecule, suggesting a more efficient compatibilizer for polymer upcycling. In conclusion, the results from these studies provide fundamental insights into the assembly of s-BCPs at homopolymer interfaces, the reduction of domain size, and promotion of adhesion, providing a strategy for the use of s-BCPs as stealth surfactants and universal compatibilizers.

interface↗
DOE OSTI2020

Using Radioiodine Speciation to Address Environmental Remediation and Waste Stream Sequestration Problems at the Fukushima Daiichi Nuclear Power Plant and a DOE Site

Iodine-129 ( 129 I), with a half-life of half-life of 16 million years, is commonly considered the single greatest risk driver in high-level and low-level nuclear repositories. This risk stems from several basic properties of 129 I, and under many geochemical conditions, it can move as an anion at nearly the rate of water through the subsurface environment. 129 I is also extremely radiologically toxic because over 90% of body burden accumulates in the thyroid, which weighs only about 14g in an adult. There is also a large worldwide inventory of radioiodine as a result of its high fission yield and this inventory is rapidly increasing as a result of nuclear energy production. Radioiodine is produced at a rate of 40 GBq (1 Ci) per gigawatt of electricity produced by nuclear power. To illustrate how the properties of 129 I magnify its risk, 129 I accounts for only 0.00002% of the radiation released from the Savannah River Site in Aiken, South Carolina, but contributes 13% of the population dose, a six orders of magnitude magnification of risk with respect to its radioactivity. The currently favored solid phase for LLW immobilization is cement, while HLW immobilization is the incorporation of waste into glass (vitrification). However, so far, the incorporation into cement and subsequent leaching of only iodide has been seriously investigated. The major problem with this approach is that it ignores the complex speciation of iodine, i.e., it ignores iodate and organo-iodine which have different chemistries. Most of the past research was devoted to the mechanism of iodide uptake in cement hydrate phases that is sorption and/or incorporation. Very few data exist on iodate and organo-iodine incorporation in cement, even though large quantities of liquid waste containing also radioiodine have already been solidified in cement Iodine-129 from low-level waste is commonly disposed of in cementitious materials. Grout, a dense cementitious fluid, mixed with a reducing slag, is often used to immobilize radionuclides. However, the reducing environment might not be conducive to immobilize iodine. For example, the silver based immobilization technologies (e.g., AgCl, Ag-impregnated granular activated carbon, Ag-mordenite) remove iodine from the aqueous phase by promoting the formation of Ag-iodide precipitates. The solubility of AgI is eight orders of magnitude lower than it is for AgIO 3 . Similarly, coprecipitation of iodine into calcium carbonate phases occurs only with IO 3 - and not with I - and org-I. If one would want to immobilize iodine more effectively, different engineering approaches would need to be used to promote binding of I - , IO 3 - , or organo-I. Using laboratory experiments with grout, slag, and silver-based adsorbents, and GC-MS and I K-edge XANES and EXAFS and C K-edge XANES spectroscopy for identifying iodine speciation, the major problems with these methods have been identified as focused too much on just one of the iodine species for immobilization, while others, especially organo-I, remained mobile. Finally, we established that most of the adsorbents that are used contain sufficient amounts of organic matter to create organo-I . It is anticipated that increased attention directed at understanding and quantifying the speciation of radioiodine, as opposed to simply total radioiodine, will lead to improved remediation results to be used for long-term radioiodine disposal in cementitious waste forms.

129-Iodine↗
DOE OSTI2023

Modeling the Effect of Surface Platinum–Tin Alloys on Propane Dehydrogenation on Platinum–Tin Catalysts

Uncertainty analysis, reported experimental literature data, and density functional theory were synthesized to model the effect of surface tin coverage on platinum-based catalysts for nonoxidative propane dehydrogenation to propylene. Here, this study tests four different platinum–tin skin surface models as potential catalytic sites, Pt 3 Sn/Pt(100), PtSn/Pt(100), Pt 3 Sn/Pt(111), and Pt 2 Sn/Pt(211), and compares them to the corresponding pure Pt surface sites using an uncertainty analysis methodology that uses BEEF-vdW with its ensembles (BMwE) to generate the uncertainty for the energies of the intermediates and transition states. One experimental data set with two experimental observations, selectivity to propylene and turnover frequency of propylene, was used as a calibration data set to evaluate the impact of the experimental data on informing the models. This study finds that the prior model for Pt 3 Sn/Pt(100) is the most active and Pt 2 Sn/Pt(211) is the most selective toward propylene. Active sites on the (100) facet have the highest probability of being responsible for C 1 and C 2 product formations (C–C bond cleavage). Increasing the Sn coverage on the (100) surface facet to a PtSn/Pt(100) active site leads to a significantly reduced rate and might explain the experimentally observed higher selectivity of Sn-doped catalysts relative to pure Pt catalysts. Next, this study finds that for all surfaces, except PtSn/Pt(100), the rate-controlling steps are the initial dehydrogenation steps alongside some partially rate-controlling second dehydrogenation steps. For PtSn/Pt(100), only the initial terminal dehydrogenation step to CH 3 CH 2 CH 2 * and second dehydrogenation steps are rate-controlling. Next, the calibrated models for all surfaces were found to be selective toward propylene production and model the reported turnover frequency successfully. Nevertheless, Pt 2 Sn/Pt(211) emerges as the active site with some (minor) evidence as the main active site based on Jeffreys’ scale interpretation of Bayes factors. This observation agrees with prior studies that also found step sites to be most likely the most relevant active sites for pure Pt catalysts. Overall, the results indicate that tin, in addition to affecting the binding strength of the adsorbed species, prevents deeper dehydrogenation (reducing coking) and cracking reactions through increasing activation barriers for unwanted side reactions.

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

Cold neutron-deuteron capture and Wigner-SU(4) symmetry

We calculate the cold neutron-deuteron (nd) capture cross section,σ nd to next-to-next-to leading order (NNLO) using the model-independent approach of pionless effective-field theory [EFT(π)]. At leading order we find σ nd = 0.314 ± 0.217 mb, while the experimental result is 0.508(15) mb for a laboratory neutron velocity of 2200 m/s. At next-to-leading-order (NLO), we show that σnd is sensitive to the low-energy constant (LEC) $L$$^{(0)}_{1}$ of the two-nucleon isovector current appearing at NLO. A fit of $L$$^{(0)}_{1}$ at NLO to the triton magnetic moment yields a NLO prediction of σ nd = 0.393 ± 0.164 mb, where the error comes from propagating the error from the $L$$^{(0)}_{1}$ fit. At NNLO, we find that a new three-nucleon magnetic moment counterterm is required for renormalization-group invariance of both σnd and the triton magnetic moment. Fitting the NNLO correction to $L$$^{(0)}_{1}$ (denoted $L$$^{(1)}_{1}$) to cold neutron-proton capture (σnp) yields a NNLO prediction of σ nd = 0.447 ± 0.130 mb, where the error comes from propagating the error from the $L$$^{(1)}_{1}$ fit. We also study different fittings of $L$$^{(0)}_{1}$ and $L$$^{(1)}_{1}$ to σ np , σ nd , and/or the triton magnetic moment. For example, fitting $L$$^{(0)}_{1}$ simultaneously to σ np , σ nd , and the triton magnetic moment at NLO, and fitting $L$$^{(1)}_{1}$ simultaneously to σ np and σnd at NNLO, yields σ nd = 0.480 ± 0.114 mb and 0.511 ± 0.042 mb, respectively, where errors are naively estimated from EFT(π) power counting. Additionally, we discuss how Wigner SU(4) symmetry may alter the naive EFT(π) expansion of σ nd .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

H-cluster Intermediates and Catalytic Properties of Clostridium pasteurianum [FeFe]-Hydrogenase III

[FeFe]-Hydrogenases are structurally diverse enzymes that catalyze reversible H2 activation at a catalytic cofactor or H-cluster. The H-cluster is a [4Fe-4S] cubane linked by a cysteine thiolate to a diiron subsite containing unique CO, CN-, and dithiomethylamine ligands. The established H-cluster resting state of [4Fe-4S]2+-[FeII-FeI], or Hox, functions in H2 binding and oxidation, or by proton-coupled reduction initiates H2 evolution. In contrast, in Clostridium pasteurianum [FeFe]-hydrogenase III (CpIII) the resting state of the H-cluster is fully oxidized, [4Fe-4S]2+-[FeII-FeII], or Hox+1. To begin to understand if Hox+1 has a role in the mechanism of CpIII, we determined the spectroscopic and redox properties of CpIII H-cluster states under catalytic conditions. CpIII poised in Hox+1 and either equilibrated under 1 atm of H2 or reduced with sodium dithionite, resulted in a mixture of reduced states including Hox (Em8 = -407 mV), Htrans-like [4Fe-4S]+-[FeII-FeII] (Em8 = -418 mV), Hred [4Fe-4S]+-[FeII-FeI], and HredH+ [4Fe-4S]2+-[FeI-FeI] (Em8 = -455-480 mV). Under H2 the population of the Htrans-like state was >20-fold higher than Hox, implicating a role in CpIII catalysis. Unlike other enzymes, there was no spectral evidence of fully reduced states, such as HsredH+ ([4Fe-4S]+-[FeI-FeI]) or Hhyd ([4Fe-4S]+-[FeII-FeII]-H-). Thus, while the H-cluster states of CpIII encompass most of the catalytic intermediates, it is either unable to form HsredH+ and Hhyd, or these states are highly destabilized in CpIII. Thus, these results demonstrate that catalytic intermediates of reduced CpIII differ from the typical intermediates of other catalytic [FeFe]-hydrogenases and may explain the catalytic preference for H2 production.

08 HYDROGEN↗
DOE OSTI2021

Semi-Automated Creation of Density Functional Tight Binding Models through Leveraging Chebyshev Polynomial-Based Force Fields

Density functional tight binding (DFTB) is an attractive method for accelerated quantum simulations of condensed matter due to its enhanced computational efficiency over standard density functional theory (DFT) approaches. However, DFTB models can be challenging to determine for individual systems of interest, especially for metallic and interfacial systems where different bonding arrangements can lead to significant changes in electronic states. In this regard, we have created a rapid-screening approach for determining systematically improvable DFTB interaction potentials that can yield transferable models for a variety of conditions. Our method leverages a recent reactive molecular dynamics force field where many-body interactions are represented by linear combinations of Chebyshev polynomials. This allows for the efficient creation of multi-center representations with relative ease, requiring only a small investment in initial DFT calculations. Here, we have focused our workflow on TiH 2 as a model system and show that a relatively small training set based on unit-cell-sized calculations yields a model accurate for both bulk and surface properties. Our approach is easy to implement and can yield reliable DFTB models over a broad range of thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

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

The Uranium-Containing and Thorium-Containing Anions Studied by Photoelectron Spectroscopy

An in-depth knowledge of actinide chemistry is fundamental to many aspects of nuclear science and technology, including the synthesis and processing of materials and the remediation of waste disposal sites. Among the actinides, the chemical bonding behaviors of actinium and thorium resemble those of the transition metals; the 5f-electrons of protactinium, uranium, neptunium, and plutonium often play important roles in their bonding; and among the still heavier elements, their bonding tends to mimic the lanthanide elements in terms of electron shielding and their f-electron contributions. Bonding that involves 5f-electrons, however, is especially important, in part because of the significance of uranium and plutonium, but also because these elements are among the few where f-electron participation in bonding is relatively common. This work focused on studying uranium-containing and thorium-containing anions in the gas phase using negative ion photoelectron spectroscopy. Since this technique directly probed valence electrons, it was uniquely positioned to address open questions regarding molecular bonding and electron configurations. A particularly important issue concerned how bonding in actinide-containing molecules was affected by modifications to their actinide atoms’ environment, i.e., due to their interaction with ligands. A closely related question was how actinide atoms’ suborbitals were qualitatively reordered and their energies quantitatively shifted as a result of their ligated environments. These were especially relevant issues in regard to uranium due to it having multiple possible oxidation states (OS) and the potential for 5f electron participation in bonding. The effects of ligands on oxidation states and 5f-orbital energies in uranium bonding was expected to be pronounced. Both ligands and excess electrons were seen as probes of actinide atoms within actinide-containing molecules. Our strategy for advancing knowledge of chemical bonding in the actinide-containing species utilized the synergy between experiments and theory, where in some cases experimental results validated theory and where in others computational results assisted in interpreting experiments. Calculations on actinide systems are terrifically challenging due to large spin-orbit interactions, relativistic effects, and just the sheer number of electrons involved. Even in the simplest species, e.g., U and U2, the most sophisticated, modern calculations carried out by the most experienced theorists often only approximate experimentally-measured values, such as electron affinities. For theory to provide confident predictions that can be used to solve real problems it needed an iterative and ultimately corrective mechanism by which its methods can develop further. Experiments can be used to identify when theory has failed; whereupon the subsequent process of using the experiment-theory interplay can be used to find the cause of the failure. Upon fixing it in one case, different test species can be proposed and studied by the experiment-theory combination to determine whether the problem has been corrected. Thus, experiments not only measure the values of molecular properties, they also provide navigational 3 beacons that keep computations off the reefs in an otherwise dark sea with few reference points. Experimental measurements in the actinide field are not only important, they are in actuality essential to computational progress. While it was not always possible to compare the theoreticallydetermined quantity of interest directly with the same experimentally-measured observable, it was usually possible to compare consequential properties that are both calculable and measurable. In the work completed here electron affinities and electronic state spacings were often sensitive consequential parameters. Reasonable agreement between measured and computational values signaled that a calculation that was very likely to be on-track. We had established collaborative relationships with five computational groups, all of which have expertise in computational actinide chemistry. Their PI’s are L. Cheng, D. Dixon, L. Gagliardi, K. Peterson, and B. Vlaisavljevich. Our close interaction with our theory partners led to us suggesting systems to them and them to us. This reciprocal interaction between our experimental and their computational results was among the most important strengths of this work and was a thread woven throughout. Even though anion photoelectron spectroscopic studies are conducted on anions, much of the information that they provide, pertains to the electronic structure of the neutral counterparts of those anions; among these are electron affinities and electronically excited state spacings. Our experimental tools included several specialized ion sources for forming the anionic species of interest, a mass spectrometer for identifying and mass-selecting them, and an anion photoelectron spectrometer for determining their electron affinities (EA) and characterizing the electronic states of the selected anions’ neutral counterparts. Anion photoelectron spectroscopy is conducted by crossing a mass-selected beam of anions with a fixed-frequency laser beam and energy-analyzing the resultant photodetached electrons. The photodetachment process is governed by the energyconserving relationship: hν = EBE + EKE, where hν is the photon’s energy, EBE is the electron binding (photodetachment transition) energy, and EKE is the electron’s kinetic energy. In our apparatus mass-selection is accomplished via time-of-flight mass spectrometry (TOF-MS), electron energy analysis is achieved with either a magnetic bottle or by velocity mapped imaging. Photodetachment of electrons from anions is implemented via either Nd:YAG or excimer lasers. The photodetachment transition energy, i.e., the EBE, between the ground vibrational and electronic state of an anion and the ground vibrational and electronic state of that anion’s neutral counterpart is the adiabatic electron affinity (EA) of that neutral molecule. Likewise, photodetachment transitions between the ground vibrational and electronic state of an anion and the various electronically-excited states of that anion’s corresponding neutral map the electronic spectrum of that neutral species, i.e., the spectral spacings in the photoelectron spectrum are a mirror image of the neutral’s electronic spectrum. It was, of course, crucial to be able to form the anionic species of interest. There, we had a particularly broad field of anion sources from which to choose. These included several variants of pulsed laser vaporization (LV), laser photoemission, infrared desorption plus photoemission, pulsed arc discharge (PACIS), electrospray ionization (ESI), and Rydberg electron transfer (RET). Each of these anion sources were readily combined with, i.e., connected to, the anion photoelectron spectroscopic portion of our apparatus as described above. Among the sources that utilize lasers, visible light for LV sources as well as IR for desorption sources are provided by Nd:YAG lasers. Ultraviolet photons are provided by both Nd:YAG and excimer lasers, whereas the excitation wavelengths for RET experiments come from two Nd:YAG-pumped dye lasers.

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

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

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

PFSA-Ionomer Dispersions to Thin-Films: Interplay Between Sidechain Chemistry and Dispersion Solvent

In fuel-cell catalyst layers (CLs), ionomers exists as nanoscale thin-films binding catalyst particles, wherein ionic and gaseous species transport to catalyst sites. To improve film function, ionomers have been designed based on the prototypical perfluorosulfonic acid (PFSA). Since CLs are fabricated through solution-based processing of inks, understanding how ink/dispersion solvent impacts ionomer interactions and function is important to guide future ionomer design. Herein, PFSA and its chemical derivatives are characterized from dispersion (in solvent) to thin film (cast on support) to investigate the impact of solvent water content on ionomer dispersion (structure, acidity) and resulting film properties. Dispersion characterization reveals that secondary aggregate structure depends on the sidechain chemistry, but all PFSA-based ionomers evolve similarly as dispersion becomes water-rich. The differences in aggregation alter ionomer adsorption and thin film morphology. Further, hydrophilic domain spacing and orientation are affected primarily by sidechain chemistry, and modulated by solvent composition. Thin-film conductivity correlates with hydration and nanodomain alignment, signifying the role of morphology in properties. Overall, this work provides insights into how chemistry can be used to coarse-tune structure–property relationships, while processing solvent is for fine-tuning of dispersion-film interplay for controlling CL-ionomer function, which can be translated to other chemistries and ink formulations.

25 ENERGY STORAGE↗
DOE OSTI2023

2023 American Society for Mass Spectrometry (ASMS) 71st Annual Conference on Mass Spectrometry and Allied Topics

Introduction (120 words max) Understanding metal-cluster chemistry occuring at solvent boundaries in the aqueous and organic phases has applications in environments from cellular processes to nuclear fuel reprocessing. Transport of metal ions at the boundaries from aqueous to organic phases involves forming a metal-ligand complex, and wherever the initial metal coordination environment is significantly different from the final one, the metal transitions through a series of transient species in passing from one phase to another. Here an investigation of the role of coordination in the chemistry of the transient species using gas-phase measurements that are free of solvent effects to better understand the binding of complexes of metals with triphenylphosphine chalcogenide ligands, examining metal-ligand homo- and hetero-dimers to better understand transient species. Methods (120 word max) Mass spectrometry and collision induced dissociation (CID) experiments were performed with a Bruker (Billerica, MA, USA) micrOTOF-Q II quadrupole time-of-flight mass spectrometer (QTOF) and Bruker amaZon speed ETD (ion trap). High resolution/high mass accuracy spectra were generated using the QTOF. External calibration was performed with Agilent (Santa Clara, CA, USA) ESI-L Low Concentration tuning mix. Both mass spectrometers were equipped with either the electrospray ionization source or nanospray sources. Metal samples were prepared between 40 – 60 uM of the metal-ligand complex in 25% water and 75% acetonitrile. Metal-ion clusters were isolated and subjected to collision induced dissociation. Density functional theory calculations were performed. Preliminary Data or Plenary Speakers Abstract (300 words max) Metal ion clusters with triphenylphosphine chalcogenide ligands were observed for group I metals with triphenylphosphine chalcogenide samples in the mass spectrum upon electrospray ionization. For each metal ligand complex of interest, the parent ion was isolated and collision induced dissociation fragmentation spectra were acquired. We observed clusters of group I metal with triphenylphosphine oxide, triphenyl phosphine sulfide, and triphenylphosphine selenide, with homodimers and heterodimer formation. In samples where the ligands were mixed, we observed mixed sodium ligand clusters at varying amounts. These mixed ligand clusters were fragmented. Metal clusters of mixed ligand dimers containing triphenylphosphine oxide showed preferential loss of the other ligand, either triphenylphosphine sulfide or triphenyl selenide. In samples with mixed triphenylphosphine sulfide and triphenylphosphine selenium ligands, sodium bound similarly between the ligands, and losses were more evenly split, showing loss ratio upon CID with losses of triphenylphosphine sulfide:triphenylphosphine selenide 43:57 ratio observed on CID. These results suggest that the oxide binds significantly more strongly than either the selenium or sulfur triphenylphosphine ligand, and the sulfur and selenium ligands are more evenly bound. Calculations were performed using density functional theory to calculate likely structures and bond energies between the group I metal and the ligands. Novel Aspect Novel analysis of sodium bound dimers with chalcogenide triphenylphosphine ligands were investigated using mass spectrometry and theoretical calculations.

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

Advanced Pt-Based Core–Shell Electrocatalysts for Fuel Cell Cathodes

Proton-exchange membrane fuel cells (PEMFCs) are highly efficient energy storage and conversion devices. Thus, the platinum group metal (PGM)-based catalysts which are the dominant choice for the PEMFCs have received extensive interest during the past couple of decades. However, the drawbacks in the existing PGM-based catalysts (i.e., high cost, slow kinetics, poor stability, etc.) still limit their applications in fuel cells. The Pt-based core–shell catalysts potentially alleviate these issues through the low Pt loading with the associated low cost and the high corrosion resistance and further improve the oxygen reduction reaction’s (ORR’s) activity and stability. This Account focuses on the synthetic strategies, catalytic mechanisms, factors influencing enhanced ORR performance, and applications in PEMFCs for the Pt-based core–shell catalysts. We first highlight the synthetic strategies for Pt-based core–shell catalysts including the galvanic displacement of an underpotentially deposited non-noble metal monolayer, thermal annealing, and dealloying methods, which can be scaled-up to meet the requirements of fuel cell operations. Subsequently, catalytic mechanisms such as the self-healing mechanism in the Pt monolayer on Pd core catalysts, the pinning effect of nitrogen (N) dopants in N-doped PtNi core–shell catalysts, and the ligand effect of the ordered intermetallic structure in L10-Pt/CoPt core–shell catalysts and their synergistic effects in N-doped L10-PtNi catalysts are described in detail. Additionally, the core–shell structure in the Pt-based catalysts have two main effects for enhanced ORR performance: (i) the interaction between Pt shells and core substrates can tune the electronic state of the surface Pt, thus boosting the ORR activity and stability, and (ii) the outer Pt shell with modest thickness can enhance the oxidation and dissolution resistance of the core, resulting in improved durability. We then review the recent attempts to optimize the ORR performance of the Pt-based core–shell catalysts by considering the shape, composition, surface orientation, and shell thickness. The factors influencing the ORR performance can be grouped into two categories: the effect of the core and the effect of the shell. In the former, PtM core–shell catalysts which use different non-PGM element cores (M) are summarized, and in the latter, Pt-based core–shell catalysts with different shell structures and compositions are described. The modifications of the core and/or shell structure can not only optimize the intermediate-binding energetics on the Pt surface through tuning the strain of the surface Pt, which increases the intrinsic activity and stability, but also offer a significantly decreased catalyst cost. Finally, we discuss the membrane electrode assembly performance of Pt-based core–shell catalysts in fuel cell cathodes and evaluate their potential in real PEMFCs for light-duty and heavy-duty vehicle applications. Even though some challenges to the activity and lifetime in the fuel cells remain, the Pt-based core–shell catalysts are expected to be promising for many practical PEMFC applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗