Results for “phase mapping” — Page 60

DOE OSTI2023

Mapping the Morphology of DNA on Carbon Nanotubes in Solution Using X-ray Scattering Interferometry

Single-walled carbon nanotubes (SWCNTs) with adsorbed single-stranded DNA (ssDNA) are applied as sensors to investigate biological systems, with potential applications ranging from clinical diagnostics to agricultural biotechnology. Unique ssDNA sequences render SWCNTs selectively responsive to target analytes such as (GT)n-SWCNTs recognizing the neuromodulator, dopamine. It remains unclear how the ssDNA conformation on the SWCNT surface contributes to functionality, as observations have been limited to computational models or experiments under dehydrated conditions that differ substantially from the aqueous biological environments in which the nanosensors are applied. We demonstrate a direct mode of measuring in-solution ssDNA geometries on SWCNTs via X-ray scattering interferometry (XSI), which leverages the interference pattern produced by AuNP tags conjugated to ssDNA on the SWCNT surface. We employ XSI to quantify distinct surface-adsorbed morphologies for two (GT) n ssDNA oligomer lengths (n = 6, 15) that are used on SWCNTs in the context of dopamine sensing and measure the ssDNA conformational changes as a function of ionic strength and during dopamine interaction. Here, we show that the shorter oligomer, (GT) 6 , adopts a more periodically ordered ring structure along the SWCNT axis (inter-ssDNA distance of 8.6 ± 0.3 nm), compared to the longer (GT) 15 oligomer (most probable 5'-to-5' distance of 14.3 ± 1.1 nm). During molecular recognition, XSI reveals that dopamine elicits simultaneous axial elongation and radial constriction of adsorbed ssDNA on the SWCNT surface. Our approach using XSI to probe solution-phase morphologies of polymer-functionalized SWCNTs can be applied to yield insights into sensing mechanisms and inform future design strategies for nanoparticle-based sensors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Radically Tunable n-Type Organic Semiconductor via Polymorph Control

Polymorphism has emerged as an important design consideration in organic semiconductors (OSCs). Previously, in many OSCs, even small changes in molecular stacking can cause drastic changes to the optical and electronic properties. However, investigation into n-type semiconductors has significantly lagged behind their p-type counterparts. In this work, we present the prolific polymorphism of 2-dimensional quinoidal terthiophene (2DQTT-o-B) and systematically investigate each of 5 polymorphs, 3 of which have been previously unreported. Grazing incidence X-ray diffraction provided a key method to understanding the structure of each polymorph. Via the polymorphic transitions mapped, we tuned the electron mobility by 5 orders of magnitude, from 5.63 × 10 –5 to 0.22 cm 2 V –1 s –1 . These were accompanied by modifications to the optical properties, namely we observed substantial differences in the refractive index noted by intensity differences under polarized optical microscopy and a large shift in optical band gap from 1.18 eV up to 1.40 eV. Lastly, we suggest that changes to these properties may be related to the unique quinoidal to aromatic transition observed in quinoidal molecules.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Artificial Intelligence-Enhanced CMIP6 Climate Projections Across the Conterminous United States

This dataset comprises high-resolution climate projections at 1/24 degree grid (~4km) over the conterminous United States (CONUS) based on ten Global Climate Models (GCMs) that are part of the Coupled Models Intercomparison Project phase 6 (CMIP6). The CMIP6 GCMs are downscaled using two artificial intelligence (AI) techniques, primarily based on the computer vision approach called super-resolution. We train two separate networks: super-resolution convolutional neural network (SRCNN) and super-resolution generative adversarial network (SRGAN). The networks are trained using Daymet observations, originally available at a 1 km resolution. For training purposes, the Daymet data is interpolated to 1/24 degree (~4km), 0.25 degree and 1 degree, which serve as high, intermediate and low-resolution inputs respectively. For each of the SRCNN and SRGAN network, we use a two-step resolution enhancement, the first step generates 4x refinement from 1 degree to 0.25 degree and the second step generates 6x refinement from 0.25 degree to 1/24 degree (~4km). We downscale daily scale precipitation, maximum temperature and minimum temperature for the six CMIP6 GCMs for 1980 to 2019 in the historical period and 2020 to 2059 in the near-term future under the shared socioeconomic pathway 585 and 245 (SSP585 and SSP245) emission scenarios. We also perform double bias-correction with Daymet observations using a quantile mapping approach, first for GCMs prior to making predictions at 1 degree grid and second after making final predictions at ~4km.

13 HYDRO ENERGY↗

DOE OSTI2024

Exciting DeePMD: Learning excited-state energies, forces, and non-adiabatic couplings

We extend the DeePMD neural network architecture to predict electronic structure properties necessary to perform non-adiabatic dynamics simulations. While learning the excited state energies and forces follows a straightforward extension of the DeePMD approach for ground-state energies and forces, how to learn the map between the non-adiabatic coupling vectors (NACV) and the local chemical environment descriptors of DeePMD is less trivial. Most implementations of machine-learning-based non-adiabatic dynamics inherently approximate the NACVs, with an underlying assumption that the energy-difference-scaled NACVs are conservative fields. We overcome this approximation, implementing the method recently introduced by Richardson [J. Chem. Phys. 158, 011102 (2023)], which learns the symmetric dyad of the energy-difference-scaled NACV. Furthermore, the efficiency and accuracy of our neural network architecture are demonstrated through the example of the methaniminium cation CH 2 NH 2 + .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Simplex‐based model for nanoparticle grain identification in four‐dimensional scanning transmission electron microscopy data

Grain identification in polycrystalline nanoparticles, for example, determining which crystal phases are present at each spatial location, is fundamental to materials characterisation. This is particularly challenging when grains overlap extensively, as commonly occurs in four-dimensional scanning transmission electron microscopy (4D-STEM) datasets. We propose a simplex-based model (SBM) in which each simplex vertex represents the diffraction pattern (DP) of a pure grain, and the simplex edges and interior represent overlapping grains. Our SBM grain identification algorithm operates on the Bragg disk (BD) data matrix distilled from the 4D-STEM data to identify the grain membership at each scan position, together with a BD feature matrix whose columns represent the DPs for each constituent grain, which is important for identifying the crystal structure of each grain. We solve the model using a two-stage algorithm. In Stage 1, we adapt a linear mixing algorithm to estimate an initial BD feature matrix whose columns represent DPs of potentially overlapping grains. Our Stage 2 algorithm incorporates sparsity considerations to transform the initial BD feature matrix so that its columns represent DPs of pure grains. Using simulated datasets with various grain configurations, we demonstrate that SBM recovers both the BD feature matrix and membership maps more accurately than existing methods, even when a grain lacks any pure region and completely overlaps with other grains.

4D-STEM segmentation↗

DOE OSTI2025

Dynamics of McMillan mappings II. axially symmetric map

Here, in this article, we investigate the transverse dynamics of a single particle in a model integrable accelerator lattice, based on a McMillan axially-symmetric electron lens. Although the McMillan e-lens has been considered as a device potentially capable of mitigating collective space charge forces, some of its fundamental properties have not been described yet. The main goal of our work is to close this gap and understand the limitations and potentials of this device. It is worth mentioning that the McMillan axially symmetric map provides the first-order approximations of dynamics for a general linear lattice plus an arbitrary thin lens with motion separable in polar coordinates. Therefore, advancements in its understanding should give us a better picture of more generic and not necessarily integrable round beams. In the first part of the article, we classify all possible regimes with stable trajectories and find the canonical action-angle variables. This provides an evaluation of the dynamical aperture, Poincaré rotation numbers as functions of amplitudes, and thus determines the spread in nonlinear tunes. Also, we provide a parameterization of invariant curves, allowing for the immediate determination of the map image forward and backward in time. The second part investigates the particle dynamics as a function of system parameters. We show that there are three fundamentally different configurations of the accelerator optics causing different regimes of nonlinear oscillations. Each regime is considered in great detail, including the limiting cases of large and small amplitudes. In addition, we analyze the dynamics in Cartesian coordinates and provide a description of observable variables and corresponding spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Metaplectic geometrical optics for ray-based modeling of caustics: Theory and algorithms

The optimization of radio frequency-wave (RF) systems for fusion experiments is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at caustics such as cutoffs and focal points, erroneously predicting the wave intensity to be infinite. This is a critical shortcoming of GO, since the caustic wave intensity is often the quantity of interest, e.g., RF heating. Full-wave modeling can be used instead, but the computational cost limits the speed at which such optimizations can be performed. Here, we have developed a less expensive alternative called metaplectic geometrical optics (MGO). Instead of evolving waves in the usual x (coordinate) or k (spectral) representation, MGO uses a mixed X$\equiv$Ax+Bk representation. By continuously adjusting the matrix coefficients A and B along the rays, one can ensure that GO remains valid in the X coordinates without caustic singularities. The caustic-free result is then mapped back onto the original x space using metaplectic transforms. Here, we overview the MGO theory and review algorithms that will aid the development of an MGO-based ray-tracing code. We show how using orthosymplectic transformations leads to considerable simplifications compared to previously published MGO formulas. We also prove explicitly that MGO exactly reproduces standard GO when evaluated far from caustics (an important property that until now has only been inferred from numerical simulations), and we relate MGO to other semiclassical caustic-removal schemes published in the literature. Finally this discussion is then augmented by an explicit comparison of the computed spectrum for a wave bounded between two cutoffs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Second Annual Progress Report on Correlation Between Microstructure and Mechanical Properties of Neutron-Irradiated Ferritic-Martensitic and Austenitic Steels

Ferritic-martensitic steels G92-2b (an optimized Grade 92 heat), NF616 and T91, and austenitic stainless steel 800H and its Grain Boundary Engineering (GBE)-treated version 800H-TMP (ThermoMechanical Processing) were irradiated in the High Flux Isotope Reactor (HFIR) of Oak Ridge National Laboratory (ORNL) and the Advanced Test Reactor (ATR) of Idaho National Laboratory (INL). Selected G92-2b samples were irradiated up to 14.66 dpa in the HFIR at two temperature ranges: 400–496.7°C and 683.3– 720°C. NF616 and T91 were irradiated in the ATR up to 8.16 dpa with the irradiation temperatures ranged from 241°C to 447.5°C. Alloy 800H and 800H-TMP samples were irradiated in both the HFIR and the ATR. Selected 800H and 800H-TMP samples had HFIR irradiation to 1.28 dpa at 580°C and ATR irradiation up to 9.12 dpa at 359°C to 431°C. Vickers hardness measurements, fractography, and microstructural characterization were performed on the selected samples in the Low Activation Materials Design and Analysis (LAMDA) laboratory. Radiation-hardening of G92-2b was observed at the lower doses and lower irradiation temperatures (400- 496.7°C), with GB03 (0.52 dpa at 400°C) and GB04 (7.44 dpa at ~490°C) showing ~12% and ~8% hardening, respectively. Softening by ~14% was observed for GB05 (14.66 dpa at 496.7°C). Radiationsoftening of G92-2b was more prevalent at the higher irradiation temperatures (683.3-~720°C), with GB10 (0.46 dpa at 683.3°C), GB11 (7.44 dpa at ~720°C), and GB12 (14.63 dpa at ~720°C) showing ~8%, ~8%, and ~40% softening, respectively. Radiation-hardening of NF616 and T91 was observed with the hardness increased by ~37% to ~65% depending on the irradiation doses and irradiation temperatures. Within the studied irradiation conditions of NF616 and T91, samples with a higher dose had a larger hardness after irradiation. All the tested alloy 800H and 800H-TMP samples in this work showed radiation-hardening by ~96±7% to ~152±10%. Alloy 800H-TMP tended to have slightly smaller radiation-hardening than alloy 800H. The fractography results of G92-2b sample GB03, GB10, and GB11, together with the previously characterized fractography of GB04, GB05, and GB12, indicated that the ductility of G92-2b was maintained up to 14.66 dpa at the lower irradiation temperatures of 400-496.7°C, while some loss of ductility (less necking) was observed for higher doses at the higher irradiation temperatures of 683.3-720°C. This agrees with the previously reported tensile test results of G92-2b, where the elongation of G92-2b was reduced at higher doses at the higher irradiation temperatures. Dimple sizes increased at higher doses, which are more evident at the higher irradiation temperatures of 683.3-~720°C. The fractography of NF616 sample D2 (2.96 dpa at 291.5°C), D4 (5.91 dpa at 359°C), and D6 (8.16 dpa at 431°C) indicated loss of ductility with negligible necking for sample D2, while ductile failure for samples D4 and D6. Fractography of alloy 800H and 800H-TMP samples in various irradiation conditions showed ductile failure with obvious necking. Dimples were observed, with some of them containing large Ti-rich particles, in all the 800H/800H-TMP samples. Electron backscatter diffraction characterization of GB12 indicated the recovery of the lath structure, which was generally replaced by an equiaxed grain structure. Transmission electron microscopy (TEM) characterization showed the presence of frequent M 23 C 6 (M = primarily Cr), MX (M = primarily V), spherical Nb(C,N) precipitates, and occasional Laves phase precipitates in the G92-2b samples. MX precipitates with sizes of 20-30 nm were observed at boundaries of smaller grains, indicating the pinning effect of the V-rich precipitates. The lath structure recovery was more evident at the higher irradiation temperatures (683.3-~720°C), with decreased densities of line dislocations and M23C6 precipitates. The irradiated T91 (TA04) showed the growth of M23C6 precipitates to 101 ± 40 nm from the initial 68 ± 22 nm in the unirradiated condition. Dislocation loops of both {100} and {111} types were present in TA04. TEM characterization was also performed on the irradiated 800H (N4, N5, N6, and AR2) and 800H-TMP (P4, P5, P6, and HG1). Accumulation of M 23 C 6 precipitates at grain boundaries was observed in all the xii 800H/800H-TMP samples, and the presence of Ti(C,N) precipitates at grain boundaries and in the matrix was observed in the irradiated 800H-TMP. Some Ti(C,N) precipitates are embedded in the M 23 C 6 precipitates, maintaining specific orientation relationships between the precipitates and between the precipitate and matrix. In addition, nanoscale Si-rich clusters were observed in the matrix of all the 800H/800H-TMP samples, with EDS Si maps tending to have a lower contrast in 800H-TMP samples. Atom probe tomography was conducted on the same samples, supported by a Rapid Turnaround Examination project under Nuclear Science User Facilities. The results are being analyzed to be integrated with the TEM results for a confident description of the γ’ precipitates. Dislocation loops also formed in all the 800H/800H-TMP samples. The density and the average size of dislocation loops were quantified to be in the order of 10 22 – 10 23 m -3 and 11.7 – 15.9 nm, respectively, in the ATR-irradiated 800H/800H-TMP samples. The loop densities in irradiated 800H were higher than that in irradiated 800H-TMP under the same irradiation conditions. Further systematic data analyses, together with some complementary experiments, will be pursued for these samples to foster peer-reviewed journal article publications.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Solving the auroral-arc-generator question by using an electron beam to unambiguously connect critical magnetospheric measurements to auroral images

Here, an active mapping mission is described that unambiguously connects measurements in the Earth's magnetosphere to visible aurora in the atmosphere. The core of the mission is an electron-beam source operated on a spacecraft in the equatorial magnetosphere, with the electron beam traveling along the Earth's magnetic-field lines to the atmosphere, depositing its energy to create an optical beam-spot in the atmosphere at the footpoint of the spacecraft's magnetic-field line. This optical spot can be imaged by ground-based cameras, putting the location of the spacecraft's magnetic footpoint into the context of the optical aurora. Scientific instruments carried on the spacecraft make critical measurements of the properties of the magnetosphere at the locations where the magnetosphere powers the aurora, allowing the determination of the plasma-physics mechanisms by which the magnetosphere drives the aurora, in particular answering the outstanding question of how the magnetosphere drives low-latitude auroral arcs. Long-standing questions in magnetosphere-ionosphere coupling that have not been answered because we could not unambiguously connect locations in the magnetosphere with their image in the ionosphere will finally be addressed. In this paper the properties of a “standard” growth-phase auroral arc are collected, theories of the magnetospheric generation of auroral arcs are reviewed, and critical magnetospheric measurements to discern the mechanisms that drive auroral arcs are determined. Further, the plasma physics of the experiment is investigated, including spacecraft-charging mitigation, beam stability, beam scattering, and electron orbit theory. Tradeoffs (keV versus MeV) concerning the energy of the electron beam are enumerated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗

DOE OSTI2024

Stress Relaxation Cracking in 347H Stainless Steel Arc Welds: Susceptibility Evaluation of Heat-Affected Zone

Stress relaxation cracking (SRC) is considered one of the major failure mechanisms for 347H stainless steel welds at elevated service temperatures or during post weld heat treatment (PWHT), especially within the heat-affected zone (HAZ). This work focuses on the characterization of SRC susceptibility within 347H physically simulated arc welded HAZ at elevated temperatures. A four-step SRC thermomechanical test in combination with finite element modeling (FEM) of the welding and testing processes is developed to establish a susceptibility map for HAZ. The test first runs a thermal cycle with three different peak temperatures (1335, 1275, and 1150 °C) to duplicate representative HAZ subzone microstructures, followed by time-to-failure examination under a variety of pre-stress (260–600 MPa) and pre-strain conditions (0.03–0.19) as a function of reheat temperatures between 750 and 1050 °C. With the aid of FEM, SRC susceptibility maps are generated to identify the threshold stress, plastic strain, and creep strain as a function of test temperature. It was found out that HAZ subzone with a lower peak temperature (1150 °C) appears to be slightly less susceptible to SRC than the other two subzones that experienced higher peak temperatures. Generally, time-to-fracture reduces with increasing initially applied stress and strain for all test temperatures. The pre-stress thresholds decrease from about 500 to 330 MPa as the testing temperature increases from 800 to 1050 °C, while the corresponding initial plastic strain thresholds reduces from 0.15 to 0.06. The SRC susceptibility was also evaluated through the Larson–Miller Parameter (LMP) analysis as a function of plastic strain, initial stress and starting stress upon reaching the testing temperature, respectively. The 1050 °C test with a high pre-applied strain (0.1) exhibits an extremely short time to failure (t = 3 s) that lies outside the general trend in LMP analysis. Additionally, it was identified that a plastic strain above 0.07 is identified to significantly reduce the bulk creep strain tolerance to fracture and therefore increases SRC susceptibility. Hardness measurement and fractography analysis indicated that the strain aging of niobium carbonitrides and other potential phases in conjunction with intergranular precipitates contributes to an increase in microhardness and increased intergranular cracking susceptibility.

14 SOLAR ENERGY↗

DOE OSTI2021

Energy-Dispersive X-ray Spectroscopy and Atom-probe Tomography Data Quantifying Component-Ratios of Multicomponent Nano-Precipitates in Ion-Irradiated Ceria

Samples of ~1 µm films of CeO2 doped with 2 wt% Mo, 1.5 wt% Ru, 0.75 wt% Pd, 0.5 wt% Re and 0.25 wt% Rh grown with pulsed laser deposition were irradiated with I2+ ions (610°C and 730°C, 1016 and 5×1016 I2+/cm2). For selected samples post-irradiation heat treatment was conducted (900°C, 1100°C). The specimens were sectioned with focussed ion beam milling and characterized in a transmission electron microscope with energy dispesive x-ray spectroscopy, and with atom-probe tomography. Energy dispersive x-ray spectroscopy was used to obtain elemental maps showing the distribution of dopants in the specimen after exposure. Some of these maps are discussed in detail in our companion article “Formation of multicomponent alloy particles in doped ceria under I2+ ion irradiation and thermal annealing” in the Journal of Nuclear Materials [REF]. Advanced computational analysis could be used to more accurately quantify local compositions. Data is provided for additional regions of interest and one additional irradiation condition. The doped Ceria film that was heat treated at 1100 °C delaminated from the substrate in most places. Samples were extracted from the underside of a delaminated piece and analyzed with atom-probe tomography. The resulting data show Ceria and a Mo-rich particle and demonstrate that this approach is feasable in principle to study local compositions in a sample exposed to such extreme conditions.

multicomponent alloys, epsilon phase, EDS, atom pr↗

DOE OSTI2025

Resolving the dynamic correlated disorder in KTa 1- x Nb x O 3

Understanding the complex temporal and spatial correlations of ions in disordered perovskite oxides is critical to rationalize their functional properties. Here, in this study, we provide new insights into the longstanding controversy regarding the off-centering of transition metal ions in the archetypal ferroelectric alloy KTa 1-x Nb x O 3 (KTN). By mapping the full energy (E) and wavevector (Q) dependence of the dynamical structure factor S(Q, E) using neutron scattering, and rationalizing our observations with atomistic simulations leveraging machine learning, we fully resolve the static vs dynamic nature of diffuse scattering sheets, as well as their composition (x) and temperature dependence. Our first-principles simulations, extended with machine-learning molecular dynamics, reproduce both inelastic neutron spectra and diffuse features, and establish how dynamically-correlated transition metal off-centerings couple to phonons, unifying local and collective viewpoints. This study sheds new light into an exemplary ferroelectric systems and shows the importance of mapping the full S(Q, E) to reveal critical spatio-temporal correlations of atomic disorder from which functional properties emerge.

42 ENGINEERING↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2026

In-Situ Magnetic Field Reconstruction in the MAGIS-100 Experiment

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.; Ferm↗

DOE OSTI2021

Spatially resolved properties of supernova host galaxies in SDSS-IV MaNGA

We crossmatch galaxies from Mapping Nearby Galaxies at Apache Point Observatory with the Open Supernova Catalog, obtaining a total of 132 SNe within MaNGA bundle. These 132 SNe can be classified into 67 Type Ia and 65 Type CC. We study the global and local properties of supernova host galaxies statistically. Type Ia SNe are distributed in both star-forming galaxies and quiescent galaxies, while Type CC SNe are all distributed along the star-forming main sequence. As the stellar mass increases, the Type Ia/CC number ratio increases. We find: (1) there is no obvious difference in the interaction possibilities and environments between Type Ia SN hosts and a control sample of galaxies with similar stellar mass and SFR distributions, except that Type Ia SNe tend to appear in galaxies which are more bulge-dominated than their controls. For Type CC SNe, there is no difference between their hosts and the control galaxies in galaxy morphology, interaction possibilities as well as environments; (2) compared to galaxy centers, the SN locations have smaller velocity dispersion, lower metallicity, and younger stellar population. This is a natural result of radius gradients for all these parameters. In conclusion, the SN location and its symmetrical position relative to the galaxy center, as well as regions with similar effective radii have very similar [Mg/Fe], gas-phase metallicity, gas velocity dispersion and stellar population age.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

The Cost of Floating Offshore Wind Energy in California Between 2019 and 2032

California’s energy planning is centered around meeting the emissions reduction and renewable energy requirements of Senate Bill 350 by 2030. However, state power system planning is expected to eventually address California’s requirement to achieve 100% of total retail electricity sales from renewable energy and zero-carbon resources by 2045, as mandated by Senate Bill 100. To comply with these directives, California needs to investigate the further development of energy efficiency, storage, and a diverse range of renewable energy, zero-carbon emission, and transmission resources, including offshore wind. Wind resources off the coast of California have the potential to generate a significant portion of the state’s electric energy as it moves toward a zero-carbon economy and can help diversify its energy mix. Floating offshore wind technology, which is suitable for the deep waters along the California coast, is currently in a precommercial phase, with approximately 84 megawatts (MW) installed worldwide at the end of 2019. Globally there are over 7,000 MW in planning and permitting phases of development, with the first commercial-scale projects expected to be operational in 2024. This study provides site-specific cost and performance data for floating offshore wind to inform California’s long-term energy planning. The identification of new resources to meet California’s policy goals at least cost is part of the Integrated Resource Planning (IRP) process, which is coordinated by the California Public Utilities Commission (CPUC). In 2019–2020 IRP modeling, offshore wind was included for the first time as a candidate resource in some sensitivity cases (CPUC 2019). The data and information presented in this report can be used to update offshore wind inputs in future IRP cycles. The authors conducted a geospatial cost analysis over portions of the offshore wind resource area of California. The analyzed spatial domain includes sites with a mean wind speed of at least 7 meters per second and water depths between 40 meters (m) and 1,300 m. Costs and energy production vary across this analysis domain. We calculated these parameters on a grid layout with over 750 sites, with each site representing a 1,000-MW commercial offshore wind power plant. Levelized cost of energy (LCOE) was calculated at each site over the analysis domain. The resulting variation in LCOE across the analysis domain is illustrated through heat maps in this report. Five study areas were selected within the analysis domain where more detailed cost analysis was conducted and cost parameters, such as annual energy production, capital cost expenditures (CapEx), operational cost expenditures (OpEx), and net capacity factors are reported. These five study areas include Morro Bay, Diablo Canyon, Humboldt, Cape Mendocino and Del Norte.

17 WIND ENERGY↗