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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Modeling of Ultrasonically Generated Liquid-Liquid Dispersions During Controlled Directional Solidification

There are innumerable two-component systems in which two very different liquid phases co-exist in equilibrium over a range of temperature and composition, e.g., oil and water, salt fluxes and solders, aluminum and lead. Often it is of practical concern to fabricate a solid component consisting of a uniform dispersion of one phase in the other. Unfortunately, uniform microstructural development during solidification of two immiscible liquids is hampered by inherent, often large, density differences between the phases that lead to severe segregation. Uniformity is also compromised by preferential wetting and coalescence phenomena. It is, however, well known that ultrasonic energy can initiate and maintain a fine liquid-liquid dispersion. The work presented here extends that observation by application of ultrasonic energy to promote uniform phase incorporation during controlled directional solidification. To this end experiments with the transparent organic, immiscible, succinonitrile-glycerol system were conducted and the numerous processing parameters associated with this technique were evaluated in view of optimizing dispersion uniformity. In view of the initial experimental results a model that predicts the dispersed liquid droplet size as a function of material properties, sample geometry, and applied energy has been developed. In the mathematical model we consider the ultrasonic field in an experimental ampoule of length L and diameter D induced by a probe having a vibration frequency of f=2OKhz (circular frequency omega = 2 pi f). The amplitude is adjustable from A=65 to 13Omicrons. The probe tip diameter is d, the liquid has a density of p, in which the speed of sound and surface tension are, respectively, c and sigma. The mathematical model and numerical investigation for the experiments [1] is done using the following assumptions: (i) The droplet size is small in comparison to the sound wave length; (ii) The forces between droplets are neglected (relative concentration is small); (iii) The droplet is stable if the kinetic energy, E(sub K), of the liquid motion due to ultrasonic field influence is less then the binding energy, E(sub S), due to the surface tension (it is easy to show that the surface energy of two droplets resulting from one is larger by about a factor of two.); (iv) The stability limit is characterized by E(sub S) to approx. E(sub K).

Grugel, R. N.↗
DOE OSTI2021

Tensor force role in β decays analyzed within the Gogny-interaction shell model

The half-life of the famous C 14 β decay is anomalously long, with different mechanisms: the tensor force, cross-shell mixing, and three-body forces, proposed to explain the cancellations that lead to a small transition matrix element. In this study, we revisit and analyze the role of the tensor force for the β decay of C 14 as well as of neighboring isotopes. We add a tensor force to the Gogny interaction, and derive an effective Hamiltonian for shell-model calculations. The calculations were carried out in a p – s d model space to investigate cross-shell effects. Furthermore, we decompose the wave functions according to the total orbital angular momentum L in order to analyze the effects of the tensor force and cross-shell mixing. The inclusion of the tensor force significantly improves the shell-model calculations of the β -decay properties of carbon isotopes. In particular, the anomalously slow β decay of C 14 can be explained by the isospin T = 0 part of the tensor force, which changes the components of N 14 with the orbital angular momentum L = 0 , 1 , and results in a dramatic suppression of the Gamow-Teller transition strength. At the same time, the description of other nearby β decays are improved. Decomposition of wave function into L components illuminates how the tensor force modifies nuclear wave functions, in particular suppression of β -decay matrix elements. Cross-shell mixing also has a visible impact on the β -decay strength. Inclusion of the tensor force does not seem to significantly change, however, binding energies of the nuclei within the phenomenological interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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↗
NASA NTRS2019

X-Ray Photoelectron Spectroscopy

The interaction of photon and the electron goes back to the early part of 19th century emanating from the photo-electric effect depicted by none other than Albert Einstein (Ref 1) described in 1905, and the redistribution of kinetic energy resulting from the interaction of x-ray and solids reported during early part of the century (Ref.2). The spectrum resolutions obtained at that time was not sufficient to observe distinct peaks in spectra for materials. Thus, these phenomena hardly attracted any attention for many years following these discoveries. The modern X-ray Photoelectron Spectroscopy (XPS) has been possible by the extensive and significant contribution from Kai Siegbahn and others (Ref.3, 4) of Uppsala University. Siegbahn developed and employed a high-resolution electron spectrometer that revealed electron peaks in a spectrum emerging from the interaction of x-rays and solids. Eventually, Kai Siegbahn received Nobel Prize in 1981 for his contributions to XPS. Around 1958, shifts in elemental peaks were realized in compounds when the same elements are bound to other but different elements. This discovery resulted in the chemical state identification in various chemicals as well as the oxidation states of atoms in compounds. Because of these useful physical effects, the Uppsala group named XPS with a synonymous name of ESCA (Electron Spectroscopy for Chemical Analysis) used widely today and will be used here alternatively. Therefore, XPS or ESCA not only identifies the element, but also the compound these elements form, from their chemical shifts. Compared to other micro-analytical techniques such as Energy Dispersive (EDS) or Wavelength Dispersive (WDS) techniques, XPS analyzes only few atomic layers present on the surface. This was discovered early in 1966 (Ref. 5). While this has awarded a merit to the analytical technique to analyze very thin layers such as films and coatings, it often analyzes the adsorbed superficial gases and contaminations on a sample introduced to its analytical chamber. This necessitates the surface is cleaned and the underlying material, material of interest, is exposed in a clean environment such that the material of interest is analyzed. The cleaning is accomplished by a scanning ion gun within the analytical chamber of the instrument. Ion gun uses an argon gas and is commonly attached in most modern machines. Reliable and efficient vacuum systems employed in modern machines does not allow adsorbed layers to rebuild after the surface is cleaned. Development of efficient and reliable vacuum pumps over these developmental years is yet another important step in the commercialization of XPS machines. Vacuum levels of better than 10-7 torr are essential to increase the mean free path of electrons released from the sample surface. Thus, modern machines are equipped with high capacity ion, turbo or cryogenic pumps in their analytical chambers. Today, XPS has advanced from an applied physics laboratory to industry for use in quality control as well as analysis of contaminants and has taken a dominant role in microanalysis. Its uniqueness arises from the fact that it is considered non-destructive compared to other common micro-analytical techniques using the electron and ion excitation sources. Polymers and plastics could be analyzed since the binding energies of saturated and unsaturated bonds in atoms could be separated. Extremely thin layers could be analyzed including materials with layered structures. The technique, though did not advance for many years, has now opened a new window for research as well as applications in industry due to its ability to separate and measure the chemical shifts in bound elements. Principles

Panda, Binayak↗
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↗