SEARCH · Search NASA

Results for “Relaxation method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Early-time γ -ray constraints on cosmic-ray acceleration in the core-collapse SN 2023ixf with the Fermi Large Area Telescope

Context. While supernova remnants (SNRs) have been considered the most relevant Galactic cosmic ray (CR) accelerators for decades, core-collapse supernovae (CCSNe) could accelerate particles during the earliest stages of their evolution and hence contribute to the CR energy budget in the Galaxy. Some SNRs have indeed been associated with TeV γ -rays, yet proton acceleration efficiency during the early stages of an SN expansion remains mostly unconstrained. Aims. The multi-wavelength observation of SN 2023ixf, a Type II supernova (SN) in the nearby galaxy M 101 (at a distance of 6.85 Mpc), opens the possibility to constrain CR acceleration within a few days after the collapse of the red super-giant stellar progenitor. With this work, we intend to provide a phenomenological, quasi-model-independent constraint on the CR acceleration efficiency during this event at photon energies above 100 MeV. Methods. We performed a maximum-likelihood analysis of γ -ray data from the Fermi Large Area Telescope up to one month after the SN explosion. We searched for high-energy, non-thermal emission from its expanding shock, and estimated the underlying hadronic CR energy reservoir assuming a power-law proton distribution consistent with standard diffusive shock acceleration. Results. We do not find significant γ -ray emission from SN 2023ixf. Nonetheless, our non-detection provides the first limit on the energy transferred to the population of hadronic CRs during the very early expansion of a CCSN. Conclusions. Under reasonable assumptions, our limits would imply a maximum efficiency on the CR acceleration of as low as 1%, which is inconsistent with the common estimate of 10% in generic SNe. However, this result is highly dependent on the assumed geometry of the circumstellar medium, and could be relaxed back to 10% by challenging spherical symmetry. Consequently, a more sophisticated, inhomogeneous characterisation of the shock and the progenitor’s environment is required before establishing whether or not Type II SNe are indeed efficient CR accelerators at early times.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

The properties of the nitrogen-vacancy center in milled chemical vapor deposition nanodiamonds

Fluorescent nanodiamonds (FNDs) containing negatively charged nitrogen-vacancy (NV − ) centers are vital for many emerging quantum sensing applications from magnetometry to intracellular sensing in biology. However, developing a scalable fabrication method for FNDs hosting color centers with consistent bulk-like photoluminescence (PL) and spin coherence properties remains a highly desired but unrealized goal. Here, we investigate optimized ball milling of single-crystal diamonds produced via chemical vapor deposition (CVD) and containing 2 ppm of substitutional nitrogen and 0.3 ppm of NV − to achieve this goal. The NV charge state, PL lifetime, and spin properties of bulk CVD diamond samples are directly compared to milled CVD FNDs and commercial high-pressure high-temperature (HPHT) FNDs. We find that on average, the relative contribution of the NV − charge state to the total NV PL is lower and the NV PL lifetime is longer in CVD FNDs compared to HPHT FNDs, both likely due to the lower N s 0 concentration in CVD FNDs. The CVD bulk and CVD FNDs on average show similar average T 1 spin relaxation times of 3.2 ± 0.7 ms and 4.7 ± 1.6 ms, respectively, compared to 0.17 ± 0.01 ms for commercial HPHT FNDs. Our results demonstrate that ball milling of CVD diamonds enables the large-scale fabrication of NV ensembles in FNDs with bulk-like T 1 spin relaxation properties.

chemical vapor deposition↗

DOE OSTI2025

Constraints on neutrino physics from DESI DR2 BAO and DR1 full shape

The Dark Energy Spectroscopic Instrument (DESI) Collaboration has obtained robust measurements of baryon acoustic oscillations in the redshift range 0.1 < 𝑧 < 4.2, based on the Lyman-𝛼 forest and galaxies from data release 2. We combine these measurements with cosmic microwave background (CMB) data from Planck and the Atacama Cosmology Telescope to place our tightest constraints yet on the sum of neutrino masses. Assuming the cosmological Λ⁢ CDM model and three degenerate neutrino states, we find ∑𝑚 𝜈 < 0.0642 eV (95%) with a marginalized error of 𝜎⁡(∑𝑚 𝜈 ) = 0.020 eV. We also constrain the effective number of neutrino species, finding 𝑁 eff = 3.2⁢3$^{+0.35}_{−0.34}$ (95%), in line with the Standard Model prediction. When accounting for neutrino oscillation constraints, we find a preference for the normal mass ordering and an upper limit on the lightest neutrino mass of 𝑚 𝑙 < 0.023 eV (95%). However, we determine using frequentist and Bayesian methods that our constraints are in tension with the lower limits derived from neutrino oscillations. Correcting for the physical boundary at zero mass, we report a 95% Feldman-Cousins upper limit of ∑𝑚 𝜈 < 0.053 eV, breaching the lower limit from neutrino oscillations. Considering a more general Bayesian analysis with an effective cosmological neutrino mass parameter, ∑𝑚 𝜈,eff , that allows for negative energy densities and removes unsatisfactory prior weight effects, we derive constraints that are in 3⁢𝜎 tension with the same oscillation limit, while the error rises to 𝜎⁡(∑𝑚 𝜈,eff ) = 0.053 eV. In the absence of unknown systematics, this finding could be interpreted as a hint of new physics not necessarily related to neutrinos. The preference of DESI and CMB data for an evolving dark energy model offers one possible solution. In the 𝑤 0 ⁢𝑤 𝑎 ⁢CDM model, we find ∑𝑚 𝜈 < 0.163 eV (95%), relaxing the neutrino tension. These constraints all rely on the effects of neutrinos on the cosmic expansion history. Using full-shape power spectrum measurements of data release 1 galaxies, we place complementary constraints that rely on neutrino free streaming. Our strongest such limit in Λ ⁢CDM, using selected CMB priors, is ∑𝑚 𝜈 < 0.193 eV (95%).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Thermal Activation of Zirconium(IV) Acetylacetonate Catalysts to Enhance Polyurethane Synthesis and Reprocessing

Carbamate formation and exchange catalysts enable efficient polyurethane (PU) manufacturing, as well as emerging recycling and reprocessing methods for PU thermosets. Zirconium β-diketonate complexes, such as Zr acetylacetonate [Zr(acac) 4 ], are effective alternatives to toxic organotin catalysts that have been used for PU reprocessing. Here, we report that Zr(acac) 4 undergoes a thermally activated process in the PU network during reprocessing that transforms it into a more active carbamate exchange catalyst. This process is associated with the irreversible loss of acetylacetonate ligands and is not observed for the more sterically hindered Zr 2,2,6,6-tetramethyl-3,5-heptanedione [Zr(tmhd) 4 ] complex. Crossover experiments between PU thermoplastics indicated enhanced carbamate exchange after the thermal activation of Zr(acac) 4 in the presence of one of the PUs, whereas a sample of Zr(acac) 4 activated in the absence of the PU had no catalytic activity. Thermal gravimetric analysis suggested that this process is associated with the loss of one protonated acac ligand. Stress relaxation analysis of PU thermosets indicated a distinct change in the characteristic relaxation time associated with the thermal activation of Zr(acac) 4 at temperatures above 140 °C; no such change was observed for samples reprocessed using Zr(tmhd) 4 . Density functional theory and molecular experiments suggest that irreversible ligand exchange of acac with alkoxide or carbamate reduces the activation energy for urethane formation and reversion. Furthermore, the Zr(acac) 4 catalyst activated in the presence of a PU’s polyol precursor provided more porous and less dense PU foams compared to those made using the unactivated Zr(acac) 4 catalyst. Furthermore, these findings are important for developing improved PU synthesis and recycling processes. Thermally activating a catalyst during reprocessing may provide more nuanced control of the in-use and reprocessing characteristics of PU thermosets.

Alcohols↗

DOE OSTI2025

High-order limiting methods using maximum principle bounds derived from the Boltzmann equation I: Euler equations

The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution bounds and limiting for the Euler equations through the kinetic representation provided by the Boltzmann equation, which allows for extending limiters designed for linear advection directly to the Euler equations. Given an arbitrary set of solution values to compute bounds over (e.g., numerical stencil) and a desired linear advection limiter, the proposed approach yields an analytic expression for the admissible region of particle distribution function values, which may be numerically integrated to yield a set of bounds for the density, momentum, and total energy. Further, these solution bounds are shown to preserve positivity of density/pressure/internal energy and, when paired with a limiting technique, can robustly resolve strong discontinuities while recovering high-order accuracy in smooth regions without any ad hoc corrections (e.g., relaxing the bounds). This approach is demonstrated in the context of an explicit unstructured high-order discontinuous Galerkin/flux reconstruction scheme for a variety of difficult problems in gas dynamics, including cases with extreme shocks and shock-vortex interactions. Furthermore, this work presents a foundation for limiting techniques for more complex macroscopic governing equations that can be derived from an underlying kinetic representation for which admissible solution bounds are not well-understood.

42 ENGINEERING↗

DOE OSTI2026

Lattice vacancy migration barriers in Fe-Ni alloys, and an indication as to why Ni atoms diffuse slowly: A first-principles study

Lattice vacancy migration barriers in ferromagnetic Fe 𝑥 ⁢Ni 1−𝑥 alloys (0.4 ≤ 𝑥 ≤ 0.6) are accurately quantified within the framework of ab initio electronic structure calculations using the nudged elastic band (NEB) method. Both the atomically disordered (A1) fcc phase, as well as the atomically ordered, tetragonal L⁢1 0 phase—which is under consideration as a material for a rare-earth-free gap magnet for advanced engineering applications—are investigated. Across an ensemble of NEB calculations performed on supercell configurations spanning a range of compositions and containing disordered, partially ordered, and fully ordered structures, we find that Ni-vacancy interchanges encounter significantly higher energetic barriers than do Fe-vacancy interchanges. We contend that this aspect is a key factor in determining the differences in mobility between Fe and Ni atoms in this ferromagnetic alloy. Moreover, we are able to interpret these findings in terms of the ferromagnetic alloy's underlying spin-polarized electronic structure. Specifically, we report a coupling between the size of local lattice distortions and the magnitude of the local electronic spin polarization around vacancies. This causes Fe atoms to relax into lattice vacancies, while Ni atoms remain rigidly fixed to their original lattice positions. These results give atomic-scale insight into the longstanding experimental observation that Ni exhibits remarkably slow atomic diffusion in Fe-Ni alloys.

density functional theory↗

DOE OSTI2026

Enhanced coercivity in Fe5C2/SiO2 core/shell nanocrystals

Rod-shaped Fe5C2 and core/shell Fe5C2/SiO2 nanocrystals were synthesized via a solution-based chemical method. Structural analysis confirmed the monoclinic phase of Fe5C2 with space group C2/c. Zero-field-cooling (ZFC) and field-cooling (FC) magnetization curves revealed distinct magnetic behaviors: uncoated Fe5C2 exhibited a low-temperature FC plateau indicative of strong dipolar interactions, while Fe5C2/SiO2 showed a monotonic increase in FC magnetization, suggesting reduced dipolar interactions due to SiO2 surface passivation. Isothermal remanent magnetization (IRM) and DC demagnetization (DCD) measurements supported this trend, with δM plots confirming weaker dipolar interactions in the coated sample. Bloch’s law fitting of temperature-dependent saturation magnetization showed a smaller Bloch’s constant for pure Fe5C2 and a larger value for Fe5C2/SiO2, reflecting enhanced surface disorder and reduced exchange coupling in the latter. Notably, Fe5C2/SiO2 demonstrated increased coercivity, attributed to decreased dipolar interaction and elevated surface anisotropy. Kneller’s law fitting yielded higher blocking temperatures for Fe5C2 (476 K) than Fe5C2/SiO2 (456 K), highlighting the impact of dipolar interactions on magnetic relaxation. These findings illustrate how SiO2 coatings effectively modulate dipolar interactions and enhance coercivity in Fe5C2 nanocrystals.

Joshi, Pramanand [Department of Physics, Universit↗

DOE OSTI2026

Systematic Construction of Time-Dependent Hamiltonians for Microwave-Driven Josephson Circuits

Time-dependent electromagnetic drives are fundamental for controlling complex quantum systems, including superconducting Josephson circuits. In these devices, accurate time-dependent Hamiltonian models are imperative for predicting their dynamics and designing high-fidelity quantum operations. Existing numerical methods, such as black-box quantization (BBQ) and energy-participation ratio (EPR), excel at modeling the static Hamiltonians of Josephson circuits. However, these techniques do not fully capture the behavior of driven circuits stimulated by external microwave drives, nor do they include a generalized approach to account for the inevitable noise and dissipation that enter through microwave ports. Here, we introduce numerical techniques that leverage classical microwave simulations, efficiently executable in finite-element solvers, to obtain the time-dependent Hamiltonian of microwave-driven superconducting circuits with arbitrary geometries under charge, flux, or mixed electromagnetic modulation. Importantly, our techniques do not rely on a lumped-element description of the superconducting circuit, in contrast to previous approaches to tackling this problem. We demonstrate the versatility of our approach by characterizing the driven properties of realistic circuit devices in complex electromagnetic environments, including coherent dynamics due to charge and flux modulation, as well as drive-induced relaxation and dephasing. Our techniques offer a powerful toolbox for optimizing circuit designs and advancing practical applications in superconducting quantum computing.

Lu, Yao [Yale U.; Yale U. (main); Fermilab] (ORCID↗

DOE OSTI2025

Modeling cross-beam energy transfer with sector ray tracing

Ray-based cross-beam energy transfer (CBET) models are an essential feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion implosions, but full 3D ray-based CBET calculations can have a prohibitively high computational cost. Sector ray tracing can be used to reduce the cost by orders of magnitude in cases where the coronal plasma and laser drive can be approximated as spherically symmetric. An extension of sector ray tracing (section ray tracing) can be used to relax the assumption of a spherically symmetric drive while still retaining most of the computational savings of a sector ray trace. We discuss the foundations of sector and section ray tracing and compare them to full ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Advances in Molecular Beam Epitaxy Growth of Ultra-Wide Bandgap Ga2O3 Based Alloys

Gallium oxide (Ga2O3) is an emerging ultra-wide bandgap semiconductor material that has attracted attention for its potential to outperform existing SiC and GaN based devices operating at high breakdown voltages and high temperature. Isovalent alloying of In and Al in Ga2O3 provides the ability to engineer bandgap energy and strain of the material. Alloying with Al increases the bandgap energy and the theoretically achievable Baliga's figure of merit, a key measure of a material's ultimate performance limits for high power switching devices. Alloying with In introduces compressive strain and can be used to counteract the tensile strain of Al incorporation. The resulting (AlxGa1-x-yIny)2O3 alloy can be lattice-matched to commercially available Ga2O3 wafers and has a tunable bandgap energy greater than that of Ga2O3, 4.76 eV. Such lattice-matched material can be grown arbitrarily thick without the detrimental effects of elastic strain and relaxation, making it suitable for high voltage diodes and transistors. However, efforts to synthesize isovalent alloys are complicated by their tendency to phase separate into corundum Al2O3 or bixbyite In2O3. Literature reports of the quaternary (AlxGa1-x-yIny)2O3 are limited to <1% unintentional indium incorporation in In-catalyzed (AlxGa1-x)2O3. The primary limitation to quaternary growth is the limited incorporation of indium at elevated growth temperatures. This limited incorporation is due to both the volatility of indium oxide and Al and Ga cation exchange reactions which replace indium in In2O3. We report on the development of a novel high-throughput molecular beam epitaxy (MBE) technique to screen the growth conditions for the ternary alloy (InyGa1-y)2O3, and the application of these findings to the first successful synthesis of phase pure monoclinic (AlxGa1-x-yIny)2O3 by MBE. By leveraging the unique sub-oxide chemistry of Ga2O3 and in-situ monitoring of crystal properties by reflection high-energy electron diffraction (RHEED), a cyclical growth and etch-back method is developed and applied to rapidly characterize the (InyGa1-y)2O3 growth space. This cyclical method provides approximately 10x increase in experimental throughput and up to 46x improvement in Ga2O3 substrate utilization. Appropriate growth conditions for monoclinic (InyGa1- y)2O3 are identified by machine learning analysis of RHEED patterns and targeted growths are characterized ex-situ to confirm improved In incorporation. These growth conditions are then combined with established (AlxGa1-x)2O3 growth conditions to grow quaternary (AlxGa1-x-yIny)2O3 with Al mole fractions ranging from 1.4% - 24.4% and In mole fractions ranging from 3.1% to 15.5%. The chemical and optical properties of the alloys are investigated by XRD, XPS, and spectroscopic ellipsometry. A lattice-matched (AlxGa1-x-yIny)2O3 alloy is examined by 4D-STEM and the chemical and physical uniformity of Al and In incorporation are discussed.

alloy↗

DOE OSTI2025

Optimization of La 2 NiO 4+δ Electrolysis Cell Oxygen Electrode through Surfactant-Enabled LaCoO 3±δ Nanocatalyst Deposition

Lanthanum nickelate (LNO) has shown promise as a Cr-resistant air electrode material for SOECs but has suboptimal surface oxygen exchange properties. Nanocoating of the LNO surface with lanthanum cobaltite (LCO) was chosen to improve cell performance as a surface oxygen conductor. The work focused on the implementation of a two-step nano-LCO film deposition utilizing catechol molecules in a porous LNO electrode. The subgoals of the work were to maintain nanosized LCO particles/ grains to increase active surface area and to control the regularity/ homogeneity of the coating across the microstructure. To achieve these goals, a novel surfactant-enhanced liquid infiltration method was utilized, where nucleation sites were spread across the electrode structure to control the location and size of LCO particles. Various catechol surfactant compositions were evaluated for their ability to control the kinetics of nanoparticle deposition and the homogeneity of the coating. Chelated LCO was characterized by X-ray diffraction (XRD), which found a substantial improvement in LCO formation with surfactant addition and determined polymerized norepinephrine to be the best-performing surfactant, with 88.4% pure LCO formed at low temperature. X-ray photoelectron spectroscopy (XPS) confirmed LCO nanostructures formed by the two-step infiltration process, showing no impurities and a stable perovskite structure. Deposition kinetics were analyzed using atomic force microscopy (AFM), correlating infiltration times and solution molarity to nanoparticle size and distribution, the results of which were confirmed in symmetrical cell samples by scanning electron microscopy (SEM). Electrochemical impedance spectroscopy (EIS) testing demonstrated substantial improvements in polarization resistance, where the nanocoating reduced the resistance by ∼55% to 0.152 Ω·cm 2 at 700 °C and 0.039 Ω·cm 2 at 800 °C. Electrical conductivity relaxation (ECR) at this temperature confirmed an improved surface oxygen exchange coefficient of the LCO + LNO heterostructure predicted by the Bode data from EIS, alongside a reduction in activation energy by about 30%.

Deposition↗

DOE OSTI2025

A Performance Portable, Fully Implicit Landau Collision Operator with Batched Linear Solvers

Modern accelerators use hierarchical parallel programming models that enable massive multithreading within a processing element (PE), with multiple PEs per device driven by traditional processes. Batching is a technique for exposing PE-level parallelism in algorithms that have traditionally run on MPI processes or multiple threads within a single process. Opportunities for batching arise in, for example, kinetic discretizations of magnetized plasmas where collisions are advanced in velocity space at each spatial point independently. This paper builds on previous work on a high-performance, fully nonlinear, Landau collision operator by batching the linear solver, as well as batching the spatial point problems and adding new support for multiple grids for multiscale, multispecies problems. An anisotropic relaxation verification test that agrees well with previously published results and analytical models is presented. The performance results from NVIDIA A100 and AMD MI250X nodes are presented with hardware utilization analysis for each architecture. Finally, the entire implicit Landau operator time advance is implemented in Kokkos for performance portability, running entirely on the device and is available in the PETSc numerical library.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Establishing criteria for the transition from kinetic to fluid modeling in hollow cathode analysis

Hollow cathodes for plasma switch applications are investigated via 2D3V particle-in-cell simulations of the channel and plume region. The kinetic nature of the plasma within the channel is dependent on the thermalization rate of electrons, emitted from the insert. When Coulomb collisions occur at a much greater rate than ionization or excitation collisions, the electron energy distribution function rapidly relaxes to a Maxwellian and the plasma within the channel can be described accurately via a fluid model. In contrast, if inelastic processes are much faster than Coulomb collisions, then the electron energy distribution function in the channel exhibits a notable high-energy tail, and a kinetic treatment is required. This criterion is applied to hollow cathodes from the literature, revealing that a fluid approach is suitable for most electric propulsion applications, whereas a kinetic treatment can be more critical to accurate modeling of plasma switches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

A Computational Framework for Simulations of Dissipative Nonadiabatic Dynamics on Hybrid Oscillator-Qubit Quantum Devices

Here, we introduce a computational framework for simulating nonadiabatic vibronic dynamics on circuit quantum electrodynamics (cQED) platforms. Our approach leverages hybrid oscillator-qubit quantum hardware with midcircuit measurements and resets, enabling the incorporation of environmental effects such as dissipation and dephasing. To demonstrate its capabilities, we simulate energy transfer dynamics in a triad model of photosynthetic chromophores inspired by natural antenna systems. We specifically investigate the role of dissipation during the relaxation dynamics following photoexcitation, where electronic transitions are coupled to the evolution of quantum vibrational modes. Our results indicate that hybrid oscillator-qubit devices, operating with noise levels below the intrinsic dissipation rates of typical molecular antenna systems, can achieve the simulation fidelity required for practical computations on near-term and early fault-tolerant quantum computing platforms.

Hamiltonians↗

DOE OSTI2024

Initial Testing of an In Situ Load Retention Aging Vessel

A thermal aging vessel instrumented with load cells was fabricated. The primary function of the vessel is to continuously monitor the in situ load retention of up to three compressed polymer coupons undergoing thermally accelerated aging under nitrogen. A secondary function is to enable gas sampling of the vessel headspace during thermal aging. Heating of the vessel is achieved using a custom heater jacket. To improve upon our conventional aging study methods which require periodic interruption of aging to perform load testing in an Instron machine at room temperature, this technology aims to automate/facilitate data acquisition/analysis, improve data quality, and enable uninterrupted compression of the polymer which represents the service condition. As an example case to assess functionality of the in situ vessel, the load retention of a siloxane elastomer material additively manufactured by direct-ink-writing (DIW) was measured at three different isothermal aging temperatures for ~1 month. Initial compression of the coupons while near the aging temperature was achieved by temporarily opening the heated vessel to access the interior chamber and manually tightening four nuts to drive the heated compression plate down onto the heated coupons. Initial testing demonstrated achievement of the primary load retention monitoring function. Unfortunately, the vessel leaked which prevented gas sampling; an active purge was used to maintain a nitrogen atmosphere. Welded or otherwise sealed joints, which could be implemented in a future design, would likely eliminate leak paths. To apply time-temperature superposition (TTS), a technique used to provide long-term prediction of the load retention from short-term isothermal data, the load retention needed to be calculated relative to the load at an estimated “equilibrium” time, after most of the transient viscoelastic physical relaxation occurred. The peak load immediately after compression could not be used as the load retention basis for two reasons: (1) age-related changes must be isolated from non-age-related physical relaxation before applying TTS and (2) the manual mechanism used to compress the specimens at the aging temperature was neither smooth nor repeatable which affected the peak load value. To better understand the effect of the mode of initial compression on the measured load, and possibly better estimate “equilibrium” physical relaxation times, systematic stress relaxation experiments were performed using an Instron machine with a thermal chamber. At a given temperature, the DIW polymer was compressed to a fixed strain in either a stepped or continuous manner at two different rates, then held at that strain for 24 hrs. The results indicated that, at a given temperature, the different stress relaxation curves appeared to converge to the same curve at some “equilibrium” time when the non-age-related physical relaxation was mostly complete. Though this observation suggests that the discontinuous manual compression employed by the vessel is feasible, a compression mechanism that is rapid, smooth, and repeatable would enhance its use.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Intermediate time sub-diffusion and stress relaxation in ring polymer melts

The slow dynamics of non-concatenated ring melts remains a frontier problem in polymer science with implications for many soft material environments including cellular biophysics. Here, in this work, we report large-scale simulations of model ring melts that analyze the monomer and center-of-mass (CM) mean square displacements (MSD) and stress relaxation function on intermediate time and length scales. The degree of dynamical slowing down is characterized by the maximally sub-diffusive fractional time scaling exponents. The data span an exceptionally wide range of ring degrees of polymerization and stiffnesses and are not successfully organized based on the classic measure linear chain entanglement, N/N e . Rather, we find that the crossover degree of polymerization, N D , based on ring macromolecular caging that successfully allows master curves to be constructed for the long-time CM self-diffusion constant also collapses these temporal dynamic scaling exponents. Different properties display different exponents and exhibit one or two regimes of linear variation with the logarithm of N D / N . A distinct crossover of the CM-MSD and stress relaxation exponents emerges at sufficiently large N or stiffness that is not found for the monomer MSD, indicating a novel form of dynamic decoupling. This crossover aligns with the predicted critical degree of polymerization for transitioning from a weak to strong caging regime, indicative of activated transport. The latter may reflect the emergence of an intermolecular collective contribution to stress in analogy with dense soft colloidal matter. Suggestions are made for future theoretical work to address the rich patterns of behavior discovered.

Anomalous diffusion↗

DOE OSTI2025

Stretched exponential magnetic relaxation dynamics in artificial square ice revealed through x-ray photon correlation spectroscopy

X-ray photon correlation scattering measurements are undertaken on a thermally active artificial spin ice based on the square lattice, referred to as artificial square ice, to probe the fluctuation timescales as a function of temperature as the system passes through the paramagnetic-antiferromagnetic phase transition, which belongs to the two-dimensional Ising universality class. In the paramagnetic regime, a single exponential timescale is seen, whereas at and below the critical temperature, a stretched exponential decorrelation is observed, with the stretching exponent decreasing from unity down to below one-half as the temperature reduces. This trend is confirmed by kinetic Monte Carlo simulations of a simplified point-dipolar square ice system, and is in agreement with past theoretical work on the kinetic Ising model where stretched exponential relaxation due to equilibrium domain wall dynamics below the critical temperature was found.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Advancing the Performance of Anion Exchange Membrane Electrolysis by Employing a Powder-Based Ionomer during Anode Catalyst Layer Fabrication

The performance of anion exchange membrane water electrolysis (AEMWE) can be significantly improved by utilizing powdered ionomers during the fabrication of the anode catalyst layer (CL) to modify the CL properties. When comparing powdered ionomers to dispersed ionomers across various catalysts including cobalt oxide (Co 3 O 4 ), nickel−iron oxide (NiFe 2 O 4 ), and iridium oxide (IrO 2 ) the anode fabricated with powdered ionomers demonstrates improved performance in polarization curves, enhanced charge transfer kinetics, and reduced ohmic and transport losses, as evidenced by voltage breakdown and electrochemical impedance spectroscopy analyses. Optimal performance is achieved using a Co 3 O 4 catalyst with a 10 wt % powdered ionomer via the catalystcoated substrate method. Microscopy analyses reveal that electrodes formed with powdered ionomers during fabrication exhibit a more uniform catalyst and ionomer distribution, increased porosity with smaller pore areas, improved electronic conduction with less catalyst agglomeration isolated by a nonconductive ionomer, and enhanced interfacial contact with the membrane and transport layer. These findings highlight that ionomers in a powdered form can promote beneficial properties and are a promising approach to improving AEMWE efficiency.

08 HYDROGEN↗