Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mix diagnostics”

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

Advanced Simulation of ITER Core X-ray Crystal Spectroscopy

X-Ray Simulation Analysis (XRSA) is an analytical ray-tracing mixed code developed specifically for the ITER Core X-Ray Crystal Spectroscopy (XRCS-Core) diagnostic, which employs a dual-reflection configuration incorporating multiple pre-reflectors made of Highly Oriented Pyrolytic Graphite (HOPG) and spherically curved analyzing crystals. The ITER XRCS-Core is designed for high spectral resolution measurement in specific wavelength ranges, including narrow bands around 1.354 Å for W 64+ , 2.19 Å for Xe 51+ , and 2.555 Å for Xe 44+ and Xe 47+ , enabling diagnostic capability across a broad electron temperature range in the ITER plasma. XRSA facilitates efficient simulation of the spectral performance of this complex X-ray spectroscopic system. Recent updates to the XRSA code have incorporated two critical effects: auto-focusing, which specifically applies to HOPG, and polarization. These two effects are particularly important in the dual-reflection configuration used in the ITER XRCS-Core system to provide more accurate modeling results. Here, simulations conducted with the updated code demonstrate that polarization has a substantial impact on the performance of the dual-reflection system. Additionally, the combined influence of polarization and system layout introduces performance variations across channels through the same crystal.

Computer graphics↗
DOE OSTI2025

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: I. Stability

Ion cyclotron emission (ICE) driven by fast ions has the potential to be a reactor relevant diagnostic of confined and escaping fast ions; moreover, controlled experiments in magnetic fusion devices can shed light on similar magnetospheric instabilities below and above the ion cyclotron frequency f ci . Dedicated DIII-D experiments studied instabilities below f ci and at cyclotron harmonics l f ci by injecting thirteen different neutral beam populations into plasmas with different thermal compositions of H, D, and 3He at five different values of toroidal field. High-frequency magnetic loops diagnose mode properties. Each beam type is individually injected for an approximate slowing-down time (~100 ms), resulting in a database of 2529 conditions that are analyzed in frequency bands corresponding to ICE harmonics and < f ci global Alfvén eigenmodes. The sub-cyclotron modes are more unstable at low field and in plasmas with large H concentration. For D injection, first harmonic ICE at f cD is more unstable in hydrogen plasmas but, for H injection, the converse is true (ICE at f cH is more unstable in D), suggesting that thermal same-species ions damp fundamental ICE. ICE stability also depends upon the pitch-angle anisotropy gradient and upon the ‘bump-on-tail’ gradient ∂f/∂v. Although the frequency of unstable ICE in these plasmas always occurs near the central cyclotron frequency, the precise value also depends on the Alfvén speed and beam direction, consistent with the idea that instability occurs for modes that simultaneously satisfy ω ≃ 2πlf ci , the magnetosonic dispersion relation, and the Doppler-shifted resonance condition.

electromagnetic ion cyclotron instabilities↗
DOE OSTI2024

Diagnosing up-scattered deuterium–tritium fusion neutrons produced in burning plasmas at the National Ignition Facility (invited)

In the push to higher performance fusion plasmas, two critical quantities to diagnose are α-heat deposition that can improve and impurities mixed into the plasma that can limit performance. In high-density, highly collisional inertial confinement fusion burning plasmas, there is a significant probability that deuterium–tritium (DT) fusion products, 14.1 MeV neutrons and 3.5 MeV α-particles, will collide with and deposit energy onto (“up-scatter”) surrounding deuterium and tritium fuel ions. These up-scattered D and T ions can then undergo fusion while in-flight and produce an up-scattered neutron (15–30 MeV). These reaction-in-flight (RIF) neutrons can then be uniquely identified in the measured neutron energy spectrum. Further, the magnitude, shape, and relative size of this spectral feature can inform models of stopping-power in the DT plasma and hence is directly proportional to α-heat deposition. In addition, the RIF spectrum can be related to mix into the burning fuel, particularly relevant for high-Z shell and other emerging National Ignition Facility platforms. The neutron time-of-flight diagnostic upgrades needed to obtain this small signal, ~10 –5 times the primary DT neutron peak, will be discussed. Results from several gain > 1 implosions will be shown and compared to previous RIF spectra. Finally, comparisons of experimental data to a simplified computational model will be made.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Thermonuclear performance variability near ignition at the National Ignition Facility

We describe our current understanding of the variability and degradation mechanisms observed through a series of five indirectly driven inertial fusion implosions fielded at the National Ignition Facility in the fall of 2021, four of which attempted to reproduce the first experiment to achieve Lawson's criterion for ignition with a thermonuclear yield of 1.35 MJ on August 8, 2021. A large number of absolutely calibrated (imaging, time-resolved, and spectrally resolved) x-ray and neutron diagnostics are fielded on the NIF along multiple lines of sight for each experiment. This allows for a reconstruction of the DT fuel and ablator mix injected into the hotspot around peak burn. We show that nuclear yield variations are well reproduced by numerical modeling when the measured low mode asymmetries and mix mass are included. Furthermore, these observed perturbations during burn are linked to small variations in laser delivery and capsule defects. Stringent specifications are then set to achieve robust ignition with the implosion design studied in this paper.

Divol, L. (ORCID:0000000269699898)↗
DOE OSTI2026

First-principles simulation of a shocked H-He mixture along the principal Hugoniot

Recent laser-shock experiments on an H–He mixture containing 11% helium (atomic fraction) have suggested the presence of an immiscibility region inside Jupiter. Reflectivity measurements were used as the primary diagnostic of H–He demixing, with discontinuities in the optical reflectivity proposed as a signature of phase separation under conditions relevant to Jupiter's interior. Here, we investigate shock-compressed H–He using ab initio molecular dynamics simulations with optical properties evaluated within the Kubo–Greenwood formalism. The equation of state and ionic configurations were obtained using the thermal Tr 2 SCANL meta-GGA exchange–correlation (XC) functional, while optical properties were computed using the recently developed RS-KDT0 range-separated thermal hybrid XC, which provides state-of-the-art accuracy for band-gap predictions in the warm dense matter regime. The calculated reflectivity shows overall good agreement with experimental measurements; however, no discontinuity is observed at elevated temperatures. Moreover, the reflectivity predictions for the mixed system are consistent with the experimental measurements in the temperature range where the mixture is inferred to be demixed. Furthermore, these results suggest that reflectivity alone may not provide a unique or sensitive diagnostic of H–He demixing at low helium concentrations under these conditions.

Density functional calculations↗
DOE OSTI2025

Free-spray characteristics and spray-wall interactions of methanol on a gasoline direct injector under flash-boiling and non-flash-boiling conditions

Methanol is considered a promising alternative fuel for internal combustion engines (ICEs) due to its high-octane number, fast laminar flame speed, and elevated latent heat of vaporization, all of which support higher compression ratios and improved thermal efficiency. However, its substantial latent heat of vaporization also poses cold-start challenges, such as misfire and fuel film deposition. This study aims to investigate methanol spray morphology and spray-wall interaction using the Spray M injector from the Engine Combustion Network within a constant-pressure flow vessel. A recently developed unified numerical framework capable of modeling both flash and non-flash boiling sprays is validated against experimental liquid volume fraction data acquired via 3-D computed tomography. Here, the results reveal that flash boiling significantly alters the spray morphology, leading to smaller droplets and spray collapse due to enhanced air-entrainment-induced turbulence. Quantitative agreement between experiments and simulations confirms this behavior. Coupled 0-D equilibrium and 3-D computational fluid dynamics analyses show that flash boiling accelerates evaporation and reduces fuel residence time, while non-flash conditions maintain a persistent liquid core more susceptible to wall wetting. Wall temperature diagnostics reveal that spray collapse alters heat transfer patterns by shifting cooling effects. Mixture fraction analysis indicates that evaporation is primarily governed by shear-layer turbulence, though deviations from adiabatic equilibrium mixing emerge under low-turbulence conditions. Finally, increasing fuel, ambient, and wall temperatures reduces wall wetting and film thickness, mitigating cold-start risks. These findings enhance the understanding of methanol sprays’ behavior and support its adoption as a viable, alternative fuel for ICEs.

Engine Combustion Network↗
DOE OSTI2025

First observations of edge instabilities in strongly shaped negative triangularity plasmas on DIII-D

Throughout the 2023 DIII-D negative triangularity (NT) campaign, unexpected magnetic fluctuations indicative of an edge localized instability were observed in the majority of strongly shaped (–0.55 ≤ δ ≤ –0.45) discharges. Moreover, these same fluctuations were found to be absent in weaker NT plasma shapes. Based on a database study, three separate mode behaviors are observed based primarily on the edge safety factor and normalized plasma beta β N . The main regime is the bursty state, in which short (~5 ms) magnetic bursts occur every 5–10 ms. Each burst drives a small reduction in the edge electron temperature (T e ) as measured by electron cyclotron emissions (ECE), which then rapidly recovers between the bursts. The quiescent regime appears to have governing mechanics similar to the bursty regime, but with multiple unique mode signatures observed in the magnetics. These additional modes greatly reduce the bursty nature of the plasma. Finally, a saturated state is observed in a number of discharges, particularly at early time slices. The saturated regime is often short-lived (maximum 500 ms), but can have a large impact on the dynamics of the plasma edge region, causing significant drops in electron temperature as observed in the ECE diagnostic. 2D ECE imaging (ECE-I) is utilized, showing large kink-like perturbations localized near the last closed flux surface that extend into the SOL indicative of an external kink/peeling behavior. Additionally, evidence of phase mixing and observed tearing-parity perturbations localized near the edge rational surfaces is observed in the quiescent regime plasmas. As a result, the various modes do not appear to have a large impact on global confinement but may play a role in regulating the transport of the NT edge and the observed edge gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Characterization of Flashback and Flame-Holding in a Jet-in-Crossflow Mixing Configuration With Methane–Hydrogen Fuel Blends

An approach that combines experimental and numerical analyses has been implemented to characterize the fundamentals of flashback events and flame-holding phenomena during high-hydrogen combustion in a jet-in-crossflow (JICF) configuration. Such flame dynamics are visualized experimentally using nanosecond (ns)-based hydroxyl planar laser-induced fluorescence (OH-PLIF) and chemiluminescence diagnostics techniques. The JICF burner has an optically accessible premixing tube allowing the optical diagnostics. The testing was conducted for varied premixer velocities (V) and equivalence ratio (ϕ) for 90–100% (H 2 , by mole) H 2 /CH 4 reactant mixtures at atmospheric temperature and pressure conditions. Two distinct flashback events were identified – conventional rich flashback and lean flashback, recorded while increasing ϕ and decreasing ϕ, respectively. The cause of lean flashback was attributed to lower momentum flux ratio which bends the jet sharply, closer toward the injection plane. The mean OH-PLIF images characterized the flame-holding behavior, where the flame was found to be stabilized on the leeward side only or on both windward and leeward sides as a lifted flame near the fuel port. A large eddy simulation (LES) with detailed chemistry combustion modeling approach was implemented along with an OH* submechanism, and it showed qualitative agreement with the integrated line-of-sight chemiluminescence results as well as the planar OH-PLIF measurements.

chemiluminescence diagnostic↗
DOE OSTI2025

Water Mass Transformation Budgets in Finite‐Volume Generalized Vertical Coordinate Ocean Models

Water Mass Transformation (WMT) theory provides conceptual tools that in principle enable innovative analyses of numerical ocean models; in practice, however, these methods can be challenging to implement and interpret, and therefore remain under-utilized. Our aim is to demonstrate the feasibility of diagnosing all terms in the water mass budget and to exemplify their usefulness for scientific inquiry and model development by quantitatively relating water mass changes, overturning circulations, boundary fluxes, and interior mixing. We begin with a pedagogical derivation of key results of classical WMT theory. We then describe best practices for diagnosing each of the water mass budget terms from the output of Finite-Volume Generalized Vertical Coordinate (FV-GVC) ocean models, including the identification of a non-negligible remainder term as the spurious numerical mixing due to advection scheme discretization errors. We illustrate key aspects of the methodology through the analysis of a polygonal region of the Greater Baltic Sea in a regional demonstration simulation using the Modular Ocean Model v6 (MOM6). We verify the convergence of our WMT diagnostics by brute-force, comparing time-averaged (“offline”) diagnostics on various vertical grids to timestep-averaged (“online”) diagnostics on the native model grid. Finally, we briefly describe a stack of xarray-enabled Python packages for evaluating WMT budgets in FV-GVC models (culminating in the new xwmb package), which is intended to be model-agnostic and available for community use and development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Predictive links between microbial communities and biological oxygen utilization in the Arctic Ocean

Microbial metabolism influences rates of net community production (NCP), exerting a direct biological control on marine oxygen and carbon fluxes. In the Arctic, it is increasingly important to understand and quantify this process, as ecological and oceanographic conditions shift due to changing climate. Here, we describe potential ecological links between pelagic microbial diversity and an NCP precursor, biological oxygen utilization, using machine learning and paired observations of community structure and metabolic activity from a seasonally and spatially variable transect of the Arctic Ocean (2019–2020 MOSAiC Expedition). Community structure was determined using 16S (prokaryotic) and 18S (eukaryotic) rRNA gene amplicon sequencing, and metabolic activity was derived from ΔO 2 /Ar. Using self-organizing maps, we identified clear successional patterns in observed microbial community structure that were seasonally driven in the upper ocean and vertically stratified with depth. Metabolic activity was also stratified, with a primarily net heterotrophic water column (median −1.5% biological oxygen saturation), excepting periodic oxygen supersaturation (maximum: 13.6%) within the mixed layer. Using DNA sequences as predictor variables, we then constructed a random forest regression model that reliably reconstructed biological oxygen concentrations (root mean squared error = 4.14 μmol kg −1 ). Top predictors from this model were from heterotrophic (bacteria) or potentially mixotrophic (dinoflagellate) taxa. These analyses highlight biologically driven diagnostic tools that can be used to expand biogeochemical datasets and improve the microbial perspectives and metabolisms represented in ecological models of net productivity and carbon flux in a changing Arctic Ocean.

Chamberlain, Emelia J. [Univ. of San Diego, San Di↗
DOE OSTI2025

Detection of trace contaminants released from rapidly heated pulsed power electrodes by laser diagnostics

Development of a capability to measure trace gas contaminants released from pulsed power electrodes would immediately impact Sandia pulsed power research. These releases occur during ultra-fast heating of metal electrodes during pulsed power discharges and can lead to substantial power losses through plasma formation. Detection of contaminants is a formidable challenge due to the need for in-situ spatially and temporally resolved measurements of trace gases in the extreme environment of ultra-fast heated metal surfaces. We investigate the feasibility of laser diagnostics for detecting contaminants, including H-atom, OH, and H 2 O. Laser-induced fluorescence and photofragmentation fluorescence showed significant plasma emission interferences and did not yield any detectable H-atom, OH, or H 2 O. Our newly developed H-atom detection using femtosecond degenerate four-wave mixing suppressed interferences and enabled detection of H-atoms. A few shots showed large signals in the near-surface region of metal foils, suggesting the formation of a wave of H-atoms from the metal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Detection and Attribution Model Intercomparison Project (DAMIP v2.0) contribution to CMIP7

The first version of the Detection and Attribution Model Intercomparison Project (DAMIP v1.0) coordinated key simulations exploring the role of individual forcings in past, current and future climate as part of the Coupled Model Intercomparison Project, Phase 6 (CMIP6). The simulations have been used extensively in the literature for detection and attribution of long-term changes, constraining projections of climate change, attributing extreme events and understanding drivers of past and future simulated climate changes. Attribution studies using DAMIP v1.0 simulations underpinned prominent assessments of human-induced warming in the Intergovernmental Panel on Climate Change (IPCC) Sixth Assessment Report. Here, we describe the set of DAMIP v2.0 simulations, proposed for the next phase of CMIP, CMIP7. Detection and attribution studies rely on pre-industrial control simulations and historical simulations, which will be part of the Diagnostic, Evaluation and Characterization of Klima (DECK) set of simulations for CMIP7. In addition, we identify the three highest-priority single-forcing experiments for CMIP7 to be run as “Assessment Fast Track” simulations in support of the Seventh Assessment Report of the IPCC: simulations with natural forcings only, anthropogenic well-mixed greenhouse gases only and anthropogenic aerosols only. Beyond this, the DAMIP v2.0 experimental design includes full-column ozone-only simulations and land-use-only simulations, such that the set of individual forcing experiments, when these are considered together, represents the full set of historical forcings. While concentration-driven simulations are prioritised for attribution, emissions-driven versions of the DAMIP experiments are also proposed to support understanding of the influence of carbon-cycle feedbacks on the simulated responses to individual forcings.

Gillett, Nathan P. [Environment and Climate Change↗
DOE OSTI2025

Application of deep learning to single-shot gas-phase laser-induced breakdown spectroscopy

Single-shot fs laser-induced breakdown spectroscopy (LIBS) has the potential to capture ns-scale electrode desorption phenomena in pulsed power fusion drivers. However, the successful implementation of the diagnostic for this purpose is challenging, as it requires interpreting single-shot measurements collected from low-density gas mixtures. In this work, we demonstrate the efficacy of a Bayesian-optimized convolutional neural network (CNN) to interpret these measurements. We generated 256 distinct measurement conditions at relevant gas pressures ranging from 80–530 mTorr by mixing 100–250 sccm H 2 and 50–200 sccm CH 4 in increments of 10 sccm. Despite the considerable overlap between signals separated by 20 sccm, the CNN is able to predict the H 2 flow rate with a root-mean-square error (RMSE) of 15.9 sccm and the CH 4 flow rate with an RMSE of 12.0 sccm. The average relative prediction error is <9% for each gas and largely remains below or near 10%.

Brown, Nathan Parnell [Sandia National Lab. (SNL-N↗
DOE OSTI2024

Quantifying Quantum Chaos through Microcanonical Distributions of Entanglement

A characteristic feature of “quantum chaotic” systems is that their eigenspectra and eigenstates display universal statistical properties described by random matrix theory (RMT). However, eigenstates of local systems also encode structure beyond RMT. To capture this feature, we introduce a framework that allows us to compare the properties of eigenstates in local systems with those of pure random states. In particular, our framework defines a notion of distance between quantum state ensembles that utilizes the Kullback-Leibler divergence to compare the microcanonical distribution of entanglement entropy (EE) of eigenstates with a reference RMT distribution generated by pure random states (with appropriate constraints). This notion gives rise to a quantitative metric for quantum chaos that not only accounts for averages of the distributions but also higher moments. The differences in moments are compared on a highly resolved scale set by the standard deviation of the RMT distribution, which is exponentially small in system size. As a result, the metric can distinguish between chaotic and integrable behaviors and, in addition, quantify and compare the of chaos (in terms of proximity to RMT behavior) between two systems that are assumed to be chaotic. We implement our framework in local, minimally structured, Floquet random circuits, as well as a canonical family of many-body Hamiltonians, the mixed-field Ising model (MFIM). Importantly, for Hamiltonian systems, we find that the reference random distribution must be appropriately constrained to incorporate the effect of energy conservation in order to describe the ensemble properties of midspectrum eigenstates. The metric captures deviations from RMT across all models and parameters, including those that have been previously identified as strongly chaotic, and for which other diagnostics of chaos such as level spacing statistics look strongly thermal. In Floquet circuits, the dominant source of deviations is the second moment of the distribution, and this persists for all system sizes. For the MFIM, we find significant variation of the KL divergence in parameter space. Notably, we find a small region where deviations from RMT are minimized, suggesting that “maximally chaotic” Hamiltonians may exist in fine-tuned pockets of parameter space. Published by the American Physical Society 2024

Physics↗
DOE OSTI2026

MRV Challenge 3: velocity and passive scalar comparison in a 3D turbulent flow

The third iteration of a magnetic resonance velocimetry (MRV) challenge comparison activity consisting of 3D flow measurements for a turbulent flow in a water channel moving past centrally positioned cubic obstacles is reported. In this challenge iteration, MRV measurements are coupled with either temperature or concentration measurements to extend the diagnostic utility. Five research teams from around the world conducted the measurements on a single shared apparatus. The water channel included partial elements along the channel sidewalls that varied in height and precluded easy optical accessibility. In addition to a turbulent mainstream flow, a secondary flow entered the channel between the second and third cubic elements from a square injector hole at the channel bottom wall. The injector jet interacts with the mainstream flow and mixes turbulently as it advects downstream. For the selected flow regime and water solutions used as working fluids, the mixing of the higher temperature or concentration through the secondary flow with the mainstream flow satisfies the same dynamics so that both temperature and concentration fields behave as passive scalars and can be directly compared. The measurements are explained in detail by each participating team, and the results are interpolated onto a common grid and compared using line profiles, contour plots, and isosurfaces.

Benson, Mike [ORNL] (ORCID:000000023210116X)↗
DOE OSTI2026

Xanthos-Lake Model Source Code

This repository contains the source code for Xanthos-Lake, a lake-modeling extension of the Xanthos framework that introduces a coupled lake component comprising the Xanthos-Lake Snow and Ice Model (xLSIM) and the Xanthos-Lake Water Balance Model (xLWBM). xLSIM is a basin-aware machine-learning model for lake snow, ice, and thermal conditions. It predicts monthly lake ice thickness, snow depth, snow-cover fraction, mixing-layer temperature, and lake ice fraction from meteorological forcing and lake surface-area information. It uses sequence-based deep-learning architectures, including Transformer and hybrid Long Short-Term Memory–Transformer (LSTM–Transformer) models, together with seasonal encoding, multi-lake learning, physical masking, and basin-level cryospheric and non-cryospheric classification. The training workflow uses Ray for scalable execution and includes optional Ray Tune hyperparameter optimization. Model predictions, observations, diagnostics, and feature-importance outputs are written in NetCDF. xLWBM is the water-balance component of the new lake framework. It simulates monthly lake storage, surface area, evaporation, inflow, outflow, and lake–groundwater exchange. It combines physical water-balance equations with calibrated bathymetric relationships, weir-based outlet flow, modified Penman open-water evaporation, groundwater head relaxation, Penman–Monteith snow and ice sublimation, and snow, ice, and thermal conditions supplied by xLSIM. The model calibrates lake parameters against satellite-derived surface-area data, using evaporation-based calibration where surface-area data are unavailable, and supports small, medium, and large lake classes. For large lakes, xLWBM is integrated with the managed-routing workflow so that lake storage and outflow interact directly with downstream river routing and reservoir operations. Together, xLSIM and xLWBM provide Xanthos with a coupled lake-modeling capability. xLSIM supplies the snow, ice, and thermal conditions that affect lake evaporation and snow- and ice-related water exchanges, while xLWBM translates those conditions into dynamic lake storage, surface area, evaporation, and discharge. In return, xLWBM supplies evolving lake surface area to xLSIM. This coupling enables Xanthos to represent lakes as active hydrologic components within basin-scale water-availability and routing simulations.

Machine Learning↗
DOE OSTI2026

Fingerprinting Uranium Oxides with Electron Energy Loss Spectroscopy Supported by Theoretical Computations

Uranium oxides occur in a variety of phases that differ in crystal structure and uranium oxidation states. Electron energy loss spectroscopy (EELS) is one of the few techniques that has sufficient spatial resolution and sensitivity to electronic structure to distinguish amongst phases at the nanoscale. However, beam-sensitive materials such as uranium oxides are subject to spectral modification due to interactions with the electron beam. Therefore, theory support is essential to reliably exclude the impact of beam damage and generate true reference datasets. Here we use a comparison of theoretical and experimental spectra to probe the impact of beam damage on O K-edge and U N-edge (N6,7 and N4,5) EELS spectra of various single-valent and mixed-valence uranium oxide bulk phases. Using a low-dose experimental set-up, we show that the O K-edge theoretical spectra are in excellent agreement with experiment for both peak positions and relative intensities of respective peaks. In contrast, U N-edge features are less distinguishing due to the partially localized nature of the U 5f orbitals and overlapping multiplet and spin–orbit coupling effects. This work demonstrates that O K-edge EELS is sufficiently diagnostic to distinguish a wide range of uranium oxides and that the experimental approach used here minimizes beam damage and allows valence state discrimination across the U(IV), U(V) and U(VI) series. When combined with imaging modes available in electron mi-croscopy, the work enables detailed investigation and characterization of uranium redox transformations at the nanoscale.

Carbone, Jacopo↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗