Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STRENGTH CALCULATIONS”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

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

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

25 ENERGY STORAGE↗
DOE OSTI2024

Melting and solidification dynamics during laser melting of reaction-based metal matrix composites uncovered by in-situ synchrotron X-ray diffraction

Laser additive manufacturing (AM) of reaction-based metal matrix composites (MMCs) involves highly complex and non-equilibrium material transformation behavior, including melting, dissolution, precipitation, and solidification. Yet, the dynamics and interplay of these phase transformation processes remain poorly understood, posing substantial challenges in identifying the microstructure formation mechanism, and predicting and controlling the microstructure in the printed parts. Here, in this study, we performed the in-situ X-ray diffraction experiment to characterize the phase evolution dynamics of the 316L + 10 vol.%TiC system during laser melting, which provides direct and quantitative insights of the complex phase reaction and evolution dynamics under rapid heating and cooling conditions relevant to additive manufacturing of reaction-based MMCs. Further in-depth thermodynamic and kinetic calculations revealed that most of the phase evolution behavior observed in the in-situ X-ray diffraction experiment cannot be solely explained by widely used equilibrium thermodynamic models, and diffusion-controlled nonequilibrium dissolution and precipitation kinetics must be considered to elucidate the complex phase evolution behavior, including incomplete TiC dissolution, and three-step TiC precipitation. The three distinct types of precipitates generate unique hierarchical TiC micro- and nanostructures, which enhances the yield strength from 513 MPa to 877 MPa by 71 %, tensile strength from 628 MPa to 1054 MPa by 68 %, and Young's modulus from 193 GPa to 221 GPa by 14 %. The findings of our research provide the knowledge foundation for the design of unique microstructures and advanced MMC materials through laser AM.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multi-stage load partitioning in additively manufactured Al6061+TiC nanocomposite characterized by in-situ neutron diffraction

Metal matrix composites (MMCs), combining metal matrix and ceramic particles, exhibit high mechanical strength and stiffness compared to conventional metallic materials. During fusion-based additive manufacturing (AM), MMCs undergo a rapid heating and cooling thermal history, which affects the interfacial bonding, thermal misfit stress and mechanical load transfer behavior between the metal matrix and particles, thus impacting the bulk mechanical performance. Here, we investigate the deformation dynamics and phase-specific load transfer in fusion-based AMed Al6061 +TiC MMCs through in-situ neutron diffraction. By calculating the phase-specific lattice strains, a multi-stage load transfer and deformation behavior during uniaxial compression is revealed: (1) elastic deformation in Stage I for both phases; (2) sudden stress rebalance between two phases in Stage II, evidenced by a sudden decrease in Al lattice strain and increase in TiC; (3) active load carrying by both phases in Stage III, with both phases experiencing an increase in lattice strains. Further, the deformation mechanism of each stage is deducted by correlating the evolutions of the interphase stresses and peaks’ broadening and intensity of the Al matrix. The local plastic deformation of the Al matrix near the phase interface is triggered and leads to stress rebalancing in Stage II. The global plastic deformation subsequently propagates throughout the Al matrix in Stage III, and the composite is further hardened through both dislocation multiplication and load transfer. The findings offer valuable insights into deformation and load-sharing behavior in AMed MMCs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Investigation into Na and Cs activity coefficients in high salt solutions to support Cs removal in Hanford tank waste

The treatment of Hanford tank waste is one of the most technically challenging environmental cleanup activities for the U.S. Department of Energy to date. To expedite the processing of liquid waste stored in underground tanks in southeastern Washington state, it is necessary to remove the significant dose contributor, 137 Cs. Toward this effort, ion exchange with crystalline silicotitanate (CST) has been employed as part of the Tank Side Cesium Removal system. The model used to predict Cs exchange onto CST was developed using activity coefficients calculated from the Bromley equation. A series of batch contact tests that varied in [Na] were conducted to look at the impact of Na concentration on Cs distribution. Experimental distribution ratios ( K d ) were compared to the distribution ratios predicted using three different activity coefficient models: (1) commercially available HSC software, (2) the Bromley equation, and (3) a simplified approach adapted from Marcos-Arroyo et al. Ultimately, the Bromley method underpredicted the effect of ionic strength on the Na activity coefficient ( γ Na+ ), HSC overestimated the impact of ionic strength on the expected performance due to the Cs activity coefficient ( γ Cs+ ), but the simplified approach predicted the experimental K d values quite well in a binary matrix. In conclusion, expansion of this approach in complex matrices is necessary for application to Hanford tank waste.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Tensile Property and Lifetime Prediction for Low-Temperature Aged Uranium-Niobium Alloys

Thermal aging models and lifetime predictions for uranium-niobium (U-Nb) alloys were created using an approach similar to those previously employed. Lifetime estimates for generic U-6Nb components were thus updated; the reported value being 800 years. This update represents a small change in lifetime vs. that of the 2012 assessment (540 years). This lifetime estimate emerged from consideration of several model fits specific to the aging datasets and properties chosen. Aging was quantified using quasi-static tensile properties measured on specimens artificially aged for up to 10 years. The major change relative to the most recent 2012 LANL assessment was that a more comprehensive body of U-Nb literature data was mined, in addition to being augmented by the latest LANL and UK AWE data. The tensile data compilation was published separately (LANL report LA-14493, December 2016). Recognizing the chemical banding of industrially produced U-6Nb, models were developed for the mid-range (6 wt.%) and extrema (4 and 8 wt.%) compositions. Lifetime estimates were calculated for all three alloy classes (4, 6, 8 wt.% nominal) and two measures of total tensile elongation (TE) to failure, namely TE-ext. — extensometer method, and TE-NCD — normalized crosshead displacement method. The conservative assumption was made that whichever composition (4, 6, or 8 wt.% Nb) and property (TE-ext or TE-NCD) was the first to cross the ductility failure threshold would limit the lifetime of the entire component. Tensile strength properties did not figure into the lifetime predictions, but could be useful as age-sensitive diagnostics and were also modeled. Of these, only first yield strength is expected to show a change at 40°C aging vs. time = 0 over the ~100-year timespan of engineering interest. Second yield strength evolves more slowly, and ultimate tensile strength slower still. Among all the models, the apparent activation energies for aging were mostly in the narrow 29– 37 kcal/mol range, which is close to that for diffusion of Nb in gamma-uranium. This agreement may be coincidental. The data from recent long-term aging studies substantially improved the model fit quality and robustness of the lifetimes. Appendices document sensitivity studies of the model fits and lifetimes with respect to using more limited datasets. These results highlight the limitations of relying solely on data from scattered literature studies and smaller datasets more generally.

36 MATERIALS SCIENCE↗
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 OSTI2021

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗
DOE OSTI2025

High-Speed and High-Quality Field Welding Repair Based on Advanced Non-Destructive Evaluation and Numerical Modeling

Creep strength-enhanced ferritic (CSEF) steels such as Grade 91 (9Cr-1Mo-V) and Grade 92 (Fe-9Cr-2W-0.5Mo) steels are widely used in the fossil-fuel-fired and nuclear power plants. The weld integrity of these steels is crucial for power plants' safe and reliable operations. Due to harsh service conditions, the steel weld can become susceptible to environmental degradation. Field welding repair is used to restore the degraded weld’s performance where a controlled temper-bead welding technique is commonly used to temper the freshly formed martensite during welding. However, knowledge of weld repairability is limited and experimental trial and error optimization to achieve desired microstructure and joint properties is expensive and time-consuming. Many existing computational models, e.g., finite element models, are limited to solving heat conduction equation and ignoring convective heat transfer due to molten metal flow. These models can result in over-prediction of peak temperatures of weld pool and heat-affected zone (HAZ), which in turn can affect the accuracy of tempering prediction. Moreover, these finite element models require an input of the deposit profiles in advance and thus limits the usability of these models. Here, a molten pool-based, multi-pass multi-layer model has been developed based on computational fluid dynamics (CFD) approach with the Volume of Fluid (VOF) method. The model calculates the bead formation, thereby eliminating the need for pre-determined bead profiles required by finite element models. For computational efficiency, a coordinate system attached to the moving heat source is utilized. A subroutine is developed to convert the temperature profiles in the reference frame stationary to the heat source to that stationary to the workpiece. The converted thermal cycles are then imported into a microstructure model to compute the tempering kinetics and resultant hardness using a Johnson-Mehl-Avrami-Kolmogorov (JMAK), and modified Grange-Baughman parameter. The modeling approach is first developed and validated on single- and multi-pass deposition of stainless steel filler metal onto a SA-533 high strength steel substrate. The models are then applied to a multi-pass V-groove repair weld of Grade 91 steel plate as well as directed energy deposition of Grade 92 steel. Non-destructive characterization of microstructures was performed on Grade 91 and 92 steel welds. Two welding processes, cold metal transfer (CMT) and flux-cored arc welding (FCAW), were investigated for the Grade 91 steel weld samples. For the Grade 92 weld samples, three different heat inputs (low, medium, and high) of gas tungsten arc welding (GTAW) were utilized to replicate traditional field welding processes. The non-destructive evaluation (NDE) method used for this research was immersion ultrasonic testing (UT) using a micro-resolution ultrasonic imaging methodology specifically designed to operate in the through-transmission configuration operating at 20 MHz of frequency. The system used a focused ultrasonic beam spot size diameter between 250-300 μm, and a 6 μm laser vibrometer spot size for detection, to produce highly defined images with longitudinal and mode-converted shear waves. From the micro-resolution ultrasonic C-scan images, three microstructural regions, i.e., weld metal (WM), HAZ, and base metal (BM), were clearly identifiable. Various levels of ultrasonic amplitudes distributed over the three regions were correlated with electron beam backscattered diffraction (EBSD) images using grain size, grain boundaries, and dislocation densities. The results showed that areas with relatively higher ultrasonic amplitude levels were associated with smaller grains and higher dislocation densities, while areas with lower amplitude levels were associated with larger grains and lower dislocation densities. In addition, ultrasonic velocity data obtained across the three different weld microstructural regions of Grade 91 test samples were correlated with optical metallographic images and hardness measurements. The results showed distinctive decreases in ultrasonic velocity and hardness over the HAZ region, where weld failures often occur during service.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tuning Spin-Orbit Torques Across the Phase Transition in VO 2 /NiFe Heterostructure

The emergence of spin-orbit torques as a promising approach to energy-efficient magnetic switching has generated large interest in material systems with easily and fully tunable spin-orbit torques. We report current-induced spin-orbit torques in VO 2 /NiFe heterostructures are investigated using spin-torque ferromagnetic resonance, where the VO 2 layer undergoes a prominent insulator-metal transition. A roughly twofold increase in the Gilbert damping parameter, α, with temperature is attributed to the change in the VO 2 /NiFe interface spin absorption across the VO 2 phase transition. More remarkably, a large modulation (±100%) and a sign change of the current-induced spin-orbit torque across the VO 2 phase transition suggest two competing spin-orbit torque generating mechanisms. The bulk spin Hall effect in metallic VO 2 , corroborated by the first-principles calculation of the spin Hall conductivity σ SH ≈ -10 4 ($\frac{\hbar}{e}$) Ω -1 m -1 , is verified as the main source of the spin-orbit torque in the metallic phase. The self-induced/anomalous torque in NiFe, with opposite sign and a similar magnitude to the bulk spin Hall effect in metallic VO 2 , can be the other competing mechanism that dominates as temperature decreases. For applications, the strong tunability of the torque strength and direction opens a new route to tailor spin-orbit torques of materials that undergo phase transitions for new device functionalities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Exploring the Effects of Node Topology, Connectivity, and Metal Identity on the Binding of Nerve Agents and Their Hydrolysis Products in Metal–Organic Frameworks

Recent studies have shown that metal–organic frameworks (MOFs) built from hexanuclear M(IV) oxide cluster nodes are effective catalysts for nerve agent hydrolysis, where the properties of the active sites on the nodes can strongly influence the reaction energetics. The connectivity and metal identity of these M6 nodes can be easily tuned, offering extensive opportunities for computational screening to predict promising new materials. Thus, we used density functional theory (DFT) to examine the effects of node topology, connectivity, and metal identity on the binding energies of multiple nerve agents and their corresponding hydrolysis products. By computing an optimization metric based on the relative binding strengths of key hydrolysis reaction species (water, agent, and bidentate-bound products), we predicted optimal M6 nodes for hydrolyzing specific nerve agent and simulant molecules, where our results are in qualitative agreement with observed experimental trends. This analysis highlighted the notion that no single metal or node topology is optimal for all possible organophosphates, suggesting that MOFs should be selected based on the agent of interest. Using the large amount of data generated from our DFT calculations, we then derived quantitative structure–activity relationship (QSAR) models to help explain the complex trends observed in the binding energies. Through linear regression, we identified the most important descriptors for describing the binding of nerve agents and their hydrolysis products to M6 nodes. These results suggested that both molecular and node properties, including both structural and chemical features, collectively contribute to the binding energetics. By performing a thorough statistical analysis, we showed that our QSAR models are capable of making quantitatively accurate binding energy predictions for nerve agents and their hydrolysis products in a wide variety of M(IV)-MOFs. In conclusion, the insights gained herein can be used to guide future experiments for the synthesis of MOFs with enhanced catalytic activity for organophosphate hydrolysis.

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

Single-photon superradiance in individual caesium lead halide quantum dots

The brightness of an emitter is ultimately described by Fermi’s golden rule, with a radiative rate proportional to its oscillator strength times the local density of photonic states. As the oscillator strength is an intrinsic material property, the quest for ever brighter emission has relied on the local density of photonic states engineering, using dielectric or plasmonic resonators. By contrast, a much less explored avenue is to boost the oscillator strength, and hence the emission rate, using a collective behaviour termed superradiance. Recently, it was proposed that the latter can be realized using the giant oscillator-strength transitions of a weakly confined exciton in a quantum well when its coherent motion extends over many unit cells. Here we demonstrate single-photon superradiance in perovskite quantum dots with a sub-100 picosecond radiative decay time, almost as short as the reported exciton coherence time. The characteristic dependence of radiative rates on the size, composition and temperature of the quantum dot suggests the formation of giant transition dipoles, as confirmed by effective-mass calculations. The results aid in the development of ultrabright, coherent quantum light sources and attest that quantum effects, for example, single-photon emission, persist in nanoparticles ten times larger than the exciton Bohr radius.

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

Temperature evolution of the phonon dynamics in the Kitaev spin liquid

In this work, we present a study of the phonon dynamics in the honeycomb Kitaev spin model at finite temperatures. We show that the fractionalized spin excitations of the Kitaev spin liquid, the itinerant Majorana fermions, and static Z 2 fluxes have distinct effects on the phonon dynamics, which makes the phonon dynamics a promising tool for exploring the Kitaev spin liquid candidate materials. In particular, we focus on the signature of the fractionalized excitations in the thermodynamic behavior of the sound attenuation and the phonon Hall viscosity: the former describes the phonon decay into the fractionalized excitations, and the latter is the leading order time-reversal symmetry breaking effect on the acoustic phonon. We find that the angular dependence of the attenuation coefficient and its magnitude is modified by the thermal excitation of the Z 2 fluxes. The strength of this effect strongly depends on the relative magnitude of the sound velocity and the Fermi velocity characterizing the low-energy Majorana fermions. We also show that the Hall viscosity is strongly suppressed by the increase of the density of the Z 2 fluxes at finite temperatures. All these observations reflect the effects of the emergent disorder on the Majorana fermions introduced by the Z 2 fluxes. Our analysis is based on the complementary analytical calculations in the low-temperature zero-flux sector, and numerical calculations in the inhomogeneous flux sectors at intermediate and high temperatures with stratified Monte Carlo (strMC) method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Dynamic mode decomposition for extrapolating nonequilibrium Green's-function dynamics

The Hartree-Fock generalized Kadanoff-Baym ansatz (HF-GKBA) offers an approximate numerical procedure for propagating the two-time nonequilibrium Green's function (NEGF). Here, using the GW self-energy, we compare the HF-GKBA to exact results for a variety of systems with long- and short-range interactions, different two-body interaction strengths, and various nonequilibrium preparations. We find excellent agreement between the HF-GKBA and exact time evolution in models when more realistic long-range exponentially decaying interactions are considered. Furthermore, this agreement persists for long times and for intermediate to strong interaction strengths. In large systems, HF-GKBA becomes prohibitively expensive for long-time evolutions. For this reason, we look at the use of dynamical mode decomposition (DMD) to reconstruct long-time NEGF trajectories from a sample of the initial trajectory. Using no more than 16% of the total time evolution, we reconstruct the total trajectory with high fidelity. Our results show the potential for DMD to be used in conjunction with HF-GKBA to calculate long-time trajectories in large-scale systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

9 Li($\textit{d,p}$) reaction as a specific probe of 10 Li, the paradigm of parity-inverted nuclei around the $\textit{N}$ = 6 closed shell

In this work, we show, within the framework of renormalized nuclear field theory and of the induced reaction surrogate formalism, that the highly debated 10 Li structure, studied in a recent high statistics 9 Li ($\textit{d,p}$) 10 Li one-neutron transfer experiment, is consistent with, or better, requires, the presence of a virtual 1/2 + state of similar single-particle strength than that of the 1/2 – resonance at 0.45 ± 0.03 MeV. Based on continuum spectroscopy self-energy techniques, we find that the physical mechanism responsible for parity inversion in $^{10}_3$Li is the same as that at the basis of the similar phenomenon observed in $^{11}_4$Be and as that needed in 11 Li to have an important $\textit{s}$-wave ground-state component. In particular the strong dynamical coupling between the $s_{1/2}$ and the $d_{5/2}$ states, mediated by the quadrupole vibration of the core 9 Li. A phenomenon which also affects the strength distribution of the $d_{5/2}$ state, in particular, in the energy range of 3–4.5 MeV. Furthermore, this mechanism is also consistent with the (normal) sequence of the $1_{p_{1/2}}$ and $2_{s_{1/2}}$ levels in the $\textit{N}$ = 7 isotones $^{12}_5$B and $^{13}_6$C. The main aim of the present Rapid Communication is that of treating structure and reactions on equal footing and in a common language. In other words, the calculation of the 9 Li ($\textit{d,p}$) 10 Li reaction as a single conceptual step from individual single-particle motion and collective vibrations to absolute double differential cross sections of renormalized virtual and resonant final states, which can be directly compared with experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Collapse of neutrino wave functions under Penrose gravitational reduction

Models of spontaneous wave function collapse have been postulated to address the measurement problem in quantum mechanics. Their primary function is to convert coherent quantum superpositions into incoherent ones, with the result that macroscopic objects cannot be placed into widely separated superpositions for observably prolonged times. Many of these processes will also lead to loss of coherence in neutrino oscillations, producing observable signatures in the flavor profile of neutrinos at long travel distances. The majority of studies of neutrino oscillation coherence to date have focused on variants of the continuous state localization model, whereby an effective decoherence strength parameter is used to model the rate of coherence loss with an assumed energy dependence. Another class of collapse models that have been proposed posit connections to the configuration of gravitational field accompanying the mass distribution associated with each wave function that is in the superposition. A particularly interesting and prescriptive model is Penrose’s description of gravitational collapse which proposes a decoherence time τ determined through E g τ ∼ ℏ , where E g is a calculable function of the Newtonian gravitational potential. Here we explore application of the Penrose collapse model to neutrino oscillations, reinterpreting previous experimental limits on neutrino decoherence in terms of this model. We identify effects associated with both spatial collapse and momentum diffusion, finding that the latter is ruled out in data from the IceCube South Pole Neutrino Observatory so long as the neutrino wave packet width at production is σ ν , x ≤ 2 × 10 − 12 m . Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2026

Impact of nitrogen on lithium-induced corrosion and cracking of RAFM steel

The emergence of private fusion enterprises, combined with advancements in U.S. and EU ITER-TBM as well as EU DEMO liquid breeder blanket concepts, has renewed focus on the materials challenges. This study examines the compatibility of reduced activation ferritic/martensitic (RAFM) steel in liquid Lithium (Li) containing various Nitrogen (N) concentrations, tested at 600 and 700 °C for 500 and 1000 h under static conditions. Microstructural characterization revealed the corrosion-induced dissolution and precipitation in Li containing 0.082 and 0.17 wt.% N at 600 °C, while limited dissolution and precipitation were observed in N-gettered Li using zirconium (Zr) and titanium (Ti) foils at 600 °C. Tensile testing at room temperature after Li exposure revealed a decrease in strength and ductility of the RAFM steel as a function of N content in Li. Severe embrittlement was observed after 500 h in Li containing 0.17 wt.% N at 600 °C and after 500 h in hot-gettered Li at 700 °C. Chromium (Cr) depletion from the steel matrix was observed, leading to the formation of M6C-type carbides on the surface identified by electron microscopy (SEM), X-ray diffraction (XRD), in agreement with thermodynamic calculations. These findings highlight the critical materials challenges posed by liquid Li breeder blanket concepts, emphasizing the need for improved strategies to address Li-induced dissolution, precipitation, and steel embrittlement due to N contamination of the coolant.

Romedenne, Marie [ORNL] (ORCID:0000000317936561)↗
DOE OSTI2025

High-Temperature Active Magnetic Bearing Development for Supercritical CO 2 Machinery Applications

Hermetic machinery utilizing gas bearings for MW-scale supercritical CO 2 (sCO 2 ) machinery applications can have significantly lower power loss and enable improved cycle efficiency compared to conventional machinery with oil-lubricated bearings. Active magnetic bearings (AMBs) are another option anticipated to have similar power loss and load capacity to gas bearings as well as offering larger mechanical tolerances, the ability to tune properties, and high reliability due to lack of mechanical wear. AMBs also have proven commercial experience at MW-scale, though environments for high-temperature sCO 2 power cycle machinery conditions are novel. Besides the potential impact of AMBs for sCO 2 turbomachinery, the technology also offers promising benefits for steam and gas turbines for power generation, compressors and expanders for industrial heat and power, and in other oil and gas and space applications. The goals of this project were to conceptual design an AMB and perform material testing. Conceptual designs for radial and thrust AMBs were produced based on a hermetically-sealed sCO 2 machinery waste-heat recovery (WHR) application for sizing and loads, and choosing a target design temperature of 540°C useful for high-temperature sCO 2 turbines for concentrating solar power (CSP) applications. Conceptual designs were initially developed for multiple radial and thrust AMBs with spreadsheet-based calculations before selecting one of each to develop further using higher-fidelity design methods for magnetic and structural performance. It was found that the radial AMB at 540°C was feasible, but the thrust AMB needed to be limited to 315°C for high-speed operation. The decision for a reduced-temperature thrust AMB was the result of several significant conclusions: 1) Hiperco 50A, originally chosen for good magnetic performance at high temperature, had insufficient strength for high-speed operation, so it was replaced with 17-4 PH. 2) The reduced magnetic performance from 17-4 PH yielded a larger bearing size, reducing the strength margin. 3) This ultimately led to a creative design implementing a more-compact E-core topology, compared to the original (conventional) C-core, and an integral shaft-disk with Hirth joint connection. These conceptual designs are unique for the size and temperature in CO 2 , relevant for MW-scale CSP applications. Long term, high temperature test data was generated for several materials, filling a void in the current body of literature. Corrosion and magnetic performance measurements were produced for PM materials (Alnico 5-7C, Alnico 9C, and SmCo) with and without nickel-coating for environments of high-temperature CO 2 up to 550°C at atmospheric pressure and 450°C at 103 bar for up to 6,000 hours. Corrosion measurements were also produced for Hiperco 50, a SM material relevant for AMB laminations, with and without C5 coating. Comparisons were also made for 450°C and 550°C, atmospheric pressure air exposures up to 5,000 hours. Results generally show that coatings can be effective at improving oxidation resistance of the bare materials, and Alnicos generally outperformed SmCo after high-temperature exposure. In addition to technical feasibility demonstrated by the design, economic feasibility was demonstrated by updating the TEA from the reference machine, re-evaluating it with CAPEX and OPEX to reflect estimated changes from process-lubricated bearings to AMBs. AMBs were shown to be comparable in performance to process-lubricated bearings, still showing a notable improvement over conventional machinery architecture with oil-lubricated bearings.

42 ENGINEERING↗
DOE OSTI2024

Multicode benchmark on simulated Ti K-edge x-ray absorption spectra of Ti-O compounds

X-ray absorption spectroscopy (XAS) is an element-specific materials characterization technique that is sensitive to structural and electronic properties. First-principles simulated XAS has been widely used as a powerful tool to interpret experimental spectra and draw physical insights. Recently, there has also been growing interest in building computational XAS databases to enable data analytics and machine learning applications. However, there are non-trivial differences among commonly used XAS simulation codes, both in underlying theoretical formalism and in technical implementation. Reliable and reproducible computational XAS databases require systematic benchmark studies. In this work, we benchmarked Ti K-edge XAS simulations of ten representative Ti-O binary compounds, which we refer to as the Ti-O-10 dataset, using three state-of-the-art codes: XSPECTRA, OCEAN and exciting. We systematically studied the convergence behavior with respect to the input parameters and developed a workflow to automate and standardize the calculations to ensure converged spectra. Our benchmark comparison considers a 35 eV spectral range starting from the K-edge onset, representative of widely used near-edge spectra. Quantitative comparison over this range is based on Spearman’s rank correlation score (r sp ). Our results show: (1) the two Bethe-Salpeter equation (BSE) codes (OCEAN and exciting) have excellent agreement with an average $r$ $sp$ of 0.998; (2) good agreement is obtained between the core-hole potential code (XSPECTRA) and BSE codes (OCEAN and exciting) with an average $r$ $sp$ of 0.990, and this smaller $r$ $sp$ reflects the noticeable differences in the main edge spectral shape that can be primarily attributed to the difference in the strength of the screened core-hole potential; (3) simulations from both methods overall reproduce well the main experimental spectral features of rutile and anatase, and the different treatments of the screened core-hole potential have visible impact on pre-edge intensities and the peak ratio of the main edge; (4) there exist moderate differences in the relative edge alignment of the three codes with a standard deviation of about 0.2 eV, which arise from multiple contributions including the frozen core approximation, final state effects, and different approximations used for the self-energy correction. In conclusion, our benchmark study provides important standards for first-principles XAS simulations with broad impact in data-driven XAS analysis.

36 MATERIALS SCIENCE↗