Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter space exploration”

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

Artificial correlation heating in PIC simulations

The Particle-in-Cell (PIC) method, a cornerstone in plasma modeling, is widely employed for its ability to simulate kinetic phenomena in device-scale domains. Part of what makes this possible is that computational macroparticles represent many physical particles. It converges under certain constraints, including a grid spacing that resolves the Debye length and a time step small enough to respect the Courant–Friedrichs–Lewy condition and plasma frequency stability limit. Here, we introduce a new constraint necessary to avoid Artificial Correlation Heating (ACH). This requires that the macroparticle coupling strength be smaller than one, Γ ω < 1, where Γ ω ≡ Γω 2/3 ⁠, Γ = Z 2 e 2 /(4πε 0 ak B T) is the physical coupling strength, and w is the macroparticle weight. This is particularly relevant to 3D simulations of dense plasmas, which are becoming common with modern computing power. If this condition is violated, the finite macroparticle weight artificially enhances the coupling strength and causes the plasma to heat until the macroparticle coupling strength is near unity, depending on the grid resolution. A comprehensive model of ACH is developed that incorporates electron density, temperature, macroparticle weight, and grid resolution. It is then tested using PIC simulations, delineating the boundaries of the method's applicability and offering a predictive framework for ACH. Moreover, the research explores a runaway heating process induced by ACH in the presence of ionization, which can lead to numerical instability. A conclusion of this study is that the onset of ACH can impose a more stringent constraint on the macroparticle weight and average number of macroparticles per cell than what is typically expected, particularly in 3D simulations of dense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Early time solution as an alternative to the late time evolving dark energy with DESI DR2 BAO

Recently the Dark Energy Spectroscopic Instrument (DESI) provided constraints on the expansion history from their Data Release 2. The DESI baryon acoustic oscillation measurements are well described by a flat Λ CDM model, but the preferred parameters are in mild ( 2.3 σ ) tension with those determined from the cosmic microwave background. The DESI Collaboration has already explored a variety of solutions to this tension relying on variations in the late-time evolution of dark energy. Here we test an alternative—the introduction of an “early dark energy” (EDE) component. We find that EDE models can alleviate the tension, though they lead to differences in other cosmological parameters that have observational implications. Particularly the EDE models that fit the acoustic datasets prefer lower Ω m , higher H 0 , n s and σ 8 in contrast to the late-time solutions. Finally, we discuss the current status and near-future prospects for distinguishing amongst these solutions.

Baryon acoustic oscillations↗
DOE OSTI2022

Cosmology with the Roman Space Telescope – Synergies with CMB lensing

ABSTRACT We explore synergies between the Nancy Grace Roman Space Telescope and CMB lensing data to constrain dark energy and modified gravity scenarios. A simulated likelihood analysis of the galaxy clustering and weak lensing data from the Roman Space Telescope High Latitude Survey combined with CMB lensing data from the Simons Observatory is undertaken, marginalizing over important astrophysical effects and calibration uncertainties. Included in the modelling are the effects of baryons on small-scale clustering, scale-dependent growth suppression by neutrinos, as well as uncertainties in the galaxy clustering biases, in the intrinsic alignment contributions to the lensing signal, in the redshift distributions, and in the galaxy shape calibration. The addition of CMB lensing roughly doubles the dark energy figure-of-merit from Roman photometric survey data alone, varying from a factor of 1.7 to 2.4 improvement depending on the particular Roman survey configuration. Alternatively, the inclusion of CMB lensing information can compensate for uncertainties in the Roman galaxy shape calibration if it falls below the design goals. Furthermore, we report the first forecast of Roman constraints on a model-independent structure growth, parametrized by σ8(z), and on the Hu-Sawicki f(R) gravity as well as an improved forecast of the phenomenological (Σ0, μ0) model. We find that CMB lensing plays a crucial role in constraining σ8(z) at z > 2, with percent-level constraints forecasted out to z = 4. CMB lensing information does not improve constraints on the f(R) model substantially. It does, however, increase the (Σ0, μ0) figure-of-merit by a factor of about 1.5.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Plasma Dynamics and Nonthermal Particle Acceleration in 3D Nonrelativistic Magnetic Reconnection

Understanding plasma dynamics and nonthermal particle acceleration in 3D magnetic reconnection has been a long-standing challenge. In this paper, we explore these problems by performing large-scale fully kinetic simulations of multi-X-line plasmoid reconnection with various parameters in both the weak- and strong-guide-field regimes. In each regime, we have identified its unique 3D dynamics that lead to field-line chaos and efficient acceleration, and we have achieved nonthermal acceleration of both electrons and protons into power-law spectra. The spectral indices agree well with a simple Fermi acceleration theory that includes guide-field dependence. In the low-guide-field regime, the flux rope kink instability governs the 3D dynamics for efficient acceleration. The weak dependence of the spectra on the ion-to-electron mass ratio and β ($\ll$1) implies that the particles are sufficiently magnetized for Fermi acceleration in our simulations. While both electrons and protons are injected at reconnection exhausts, protons are primarily injected by perpendicular electric fields through Fermi reflections and electrons are injected by a combination of perpendicular and parallel electric fields. The magnetic power spectra agree with in situ magnetotail observations, and the spectral index may reflect a reconnection-driven size distribution of plasmoids instead of the Goldreich–Sridhar vortex cascade. As the guide field becomes stronger, the oblique flux ropes of large sizes capture the main 3D dynamics for efficient acceleration. Intriguingly, the oblique flux ropes can also experience flux rope kink instability, to drive extra 3D dynamics. This work has broad implications for 3D reconnection dynamics and particle acceleration in heliophysics and astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Estimating List-Mode Data Sensitivities to Nuclear Data with MCNP6

Nuclear data are a vital component of predictive simulations used in applications like experiment design, stockpile stewardship, nuclear nonproliferation/safeguards, health physics, and criticality safety. A singular simulation requires the coalescence of different areas of nuclear data such as cross sections, angular distributions, and energy distributions of emitted neutrons for different materials and energy ranges. Improving nuclear data and thus reducing the uncertainty in simulated parameters could enable smaller, better-informed safety factors and ultimately reduce operational and procedural costs. There is a constant effort to garner a better understanding of the physical quantities represented by nuclear data through experiments. Integral experiment benchmarks use simulated and measured results to validate current nuclear data values. In the past, benchmarks primarily focused on the effective multiplication factor (k eff ); however, this limited scope has caused compensating errors and areas of nuclear data that lack validation. Compensating errors are inaccuracies in nuclear data that are obfuscated by cancellation when observing integrated values such as k eff . Diverse integral benchmark experiments that look for quantities of interest other than k eff and include multiple responses minimize the possibility of compensating errors and provides validation to areas of nuclear data previously lacking experimental validation. Benchmark experiments can be optimized during the design process to be highly dependent on specific areas of nuclear data. The dependence of a response in an experiment to a specific area/type of nuclear data is defined as sensitivity. A larger sensitivity means that nuclear data uncertainties will play a larger role in the response(s) resulting in larger bias. Currently, the sensitivity capabilities of the Monte Carlo N-Particle (MCNP ®1 ) transport code are limited to responses of k eff and tallied values (e.g., flux, surface current). As a part of the EUCLID project, this work explores estimating list-mode nuclear data sensitivities that can be used to design experiments aimed to constrain and reduce compensating errors in nuclear data by focusing on responses other than k eff . Tallied values are ideal quantities that are estimated with detectors during experiments. List-mode data (a list of neutron collection times) are the direct output of detector systems in subcritical neutron noise experiments. Expanding MCNP sensitivity capabilities to include the sensitivity of responses estimated from list-mode data, such as the prompt neutron decay constant (α) and multiplicity estimates (S and D), enables more direct comparison of simulated and measured experimental quantities. Additionally, deterministic tools such as SENSMG are capable of obtaining sensitivities to a wide variety of responses; however, these tools cannot handle complex geometries due to the assumptions made in discretizing the phase-space variables of the Boltzman transport equation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

On the Existence of the Plateau Emission in High-energy Gamma-Ray Burst Light Curves Observed by Fermi-LAT

The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) shows long-lasting high-energy emission in many gamma-ray bursts (GRBs), similar to X-ray afterglows observed by the Neil Gehrels Swift Observatory (Swift). Some LAT light curves (LCs) show a late-time flattening reminiscent of X-ray plateaus. We explore the presence of plateaus in LAT temporally extended emission analyzing GRBs from the second Fermi-LAT GRB Catalog from 2008 to 2016 May with known redshifts, and check whether they follow closure relations corresponding to four distinct astrophysical environments predicted by the external forward shock model. We find that three LCs can be fit by the same phenomenological model used to fit X-ray plateaus and show tentative evidence for the existence of plateaus in their high-energy extended emission. The most favorable scenario is a slow-cooling regime, whereas the preferred density profile for each GRBs varies from a constant-density interstellar medium to an r –2 wind environment. We also compare the end time of the plateaus in γ-rays and X-rays using a statistical comparison with 222 Swift GRBs with plateaus and known redshifts from 2005 January to 2019 August. Within this comparison, the case of GRB 090510 shows an indication of chromaticity at the end time of the plateau. Finally, we update the 3D fundamental plane relation among the rest-frame end time of the plateau, its correspondent luminosity, and the peak prompt luminosity for 222 GRBs observed by Swift. We find that these three LAT GRBs follow this relation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Evaluating Variable-Impedance Magnetically-Insulated Transmission Lines as a Risk-Mitigation Measure for Next-Generation Pulsed Power

This project has produced the first detailed characterizations of power flow resulting from applying the “variable-impedance MITL” concept to real-life systems in Sandia’s pulsed power program (Z and next-generation pulsed power (NGPP)). We present simulation results and analyses for constant-impedance versions of both Z and NGPP and survey the operational viability of several variable-impedance re-designs in the parameter space of linear tapers. Circuit modeling (SCREAMER/Bertha) was used to pinpoint promising candidate designs, and EM-PIC (Empire) simulations were used to evaluate these candidates more rigorously. This approach was particularly successful in the Z regime which resulted in the identification of several viable variable-impedance MITL designs for each level. The approach was more challenged in the operating space NGPP occupies, producing data points that speak to a more restrictive design space due to anode plasma turn-on. In the end, we were able to converge on one viable variable-impedance design for the highest inductance line (level “F”) and one for the highest current line (level “A”). Altogether, the body of simulation evidence presented in this report suggest there does exist flexibility in operating space for magnetically-insulated transmission lines (MITLs) having variable geometric impedance to be a potential enabling technology for safely increasing current delivery (and potentially lowering stack voltage) in pulsed-power drivers by manipulating electron losses; however, operating points for a particular design must be carefully screened. Circuit and EM-PIC modeling provided consistent verdicts in safe operating regimes for operational viability, but additional physics such as anode plasma turn-on which is included in Empire but not in SCREAMER/Bertha was found to be a critical factor affecting power flow that lead to different assessments between the codes. It is not always the case that the occurrence of anode plasma caused a design to fail (some designs turned on anode plasma yet still delivered load currents meeting design targets); the details matter such as how early in the pulse anode surfaces break down (and how large a region). However, in every case that it did fail it was found that the feedback from anode plasma was the cause (i.e., turning off the anode plasma model in Empire restored agreement with the circuit model prediction). As circuit simulations represent an efficient and practical means of surveying design space compared to more computationally-expensive approaches such as EM-PIC, it could be prudent to invest in the research and development of models to include the effects of anode plasma such as ion emission in circuit codes. The variable-impedance MITL design is a new concept that enables controlled manipulation of the initial electron losses in the outer MITL and can be tested on Z today. We encourage follow-on work to explore further optimization (including alternative variable-impedance profiles, e.g., having constant dZ/dR), and to confirm the major findings presented in this report by fielding test hardware on actual Z shots.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Advancement of Commercial Intumescent Expanding Foams for Deactivation and Decommissioning in the Nuclear Sector - 20198

Florida International University (FIU), in collaboration with The Department of Energy's Office of Environmental Management (DOE-EM), Savannah River National Laboratory (SRNL), and sites across the Savannah River complex, have identified an operational requirement for a fixative technology that is intended to immobilize and/or isolate residual contamination within a 3-dimensional space. Fixation of radiological contamination can reduce worker risk and mitigate potentially hazardous conditions, however nearly every marketed contamination fixative has been found to be flammable; a significant concern in radiological facilities. Coupled with this, industry fixatives are normally used as a thin coating which can present problems when attempting to stabilize irregular geometry or areas that are difficult to access whilst ensuring full coverage. The technical evaluation and advancement of commercial-off-the-shelf (COTS) polyurethane foams has yielded a down-selected candidate that shows potential in meeting the requirements to support deactivation and decommissioning activities. Several performance criteria have been established and tested to progress the technology readiness level towards an active field demonstration (TRL-7). Such criteria include: mechanical failure limits, adhesive and cohesive properties, thermal/fire resilience, determining thermal behavior, ability to immobilize contamination, and a means of non-destructive evaluation of applications. The test scenario examined was targeted towards an application for decommissioning nuclear pipework, in which the down-selected polyurethane foam would act as a barrier to segregate pipework and mitigate the potential for release during cutting, packaging, and storage operations. Testing carried out at SRNL included: mechanical evaluation of tensile, compressive, and adhesion strength by dynamic mechanical analysis (DMA), as well as thermogravimetric analysis (TGA). FIU examined the foam's fixative properties by utilizing phosphorescent europium-dysprosium doped strontium aluminate powder to investigating the extent to which contamination can be immobilized. FIU has also exploited previous successes in the field of intumescent technologies to assess the down-selected foam's tolerance to an extreme fire scenario, while maintaining the ability to effectively mitigate a contamination release. Parallel to this, extensive thermal investigations were carried out to determine the upper boundary of anticipated heat generation during the curing process as heat generation has the potential to compromise rubber parts of contaminated enclosures. These investigations subsequently yielded a promising method for a non-destructive application evaluation by means of infrared thermography. Utilizing the high sensitivity of modern IR cameras, coupled with the heat generated during the curing process of the polyurethane foam, FIU has been exploring the concept of monitoring the external pipe surface for indications of an irregular or abnormal application, thus informing operational decision making. The testing carried out utilized several current 'best fit' ASTM standards, which serve as helpful guidelines for testing, however, a precise definition of the operational parameters and requirements is still necessary. With continued collaboration with SRNL, FIU aims refine said definitions and develop new standards by which this, and other decommissioning technologies, can be accredited by relevant standards based testing. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Using the Carbon Capture Simulation Initiative (CCSI) Tool to Design the Experiments in the Parametric Campaign of a Novel Compact Absorber for Carbon Capture

Gas absorption towers with structured packing and solvent have been used for Carbon Dioxide (CO 2 ) Capture for about many decades. To overcome process limitations and practical disadvantages for CO 2 capture from the stationary emitter (e.g. NG and coal power plant), many new designs have been proposed and explored at the various scales in the last decade with aim of either low energy penalty or low capital cost. To reduce the size of the absorption tower and hence the total cost of CO 2 capture, the University of Kentucky Center for Applied Energy Research Center (UK CAER) has designed and built a novel CO 2 capture absorption tower or Compact Absorber, integrated into an existing large-bench scale CO 2 capture unit. The Compact Absorber has three sections. The top of the column is a fogging section where the solvent is sprayed through a nozzle producing droplets flowing downward in a co-current fashion with the flue gas. The center of the column is a frothing section where the solvent and flue gas flow through regenerative frothing screens designed by Industrial Climate Solutions, Inc. The bottom of the column is a typical structured packing section were the flue gas and solvent flow in a counter-current fashion. The parametric campaign will be conducted in order to optimize the operating parameters for CO 2 capture including liquid/gas ratio, lean loading, and temperature, liquid residence time. A simulated flue gas with 14% CO 2 will be used along with a UK CAER developed proprietary solvent. The 100-hour parametric campaign is designed using a statistical approach of the Sequential Design of Experiments (sDOE). sDOE is one of the CCSI tools that provides an adaptive statistical approach for designing future experiments based on the results of previous experiments. Application of a typical DOE provides the user with the minimum number of experiments required to get the same data, but sDOE allows the user to make an informed choice of experiments based on the results of previous experiments. The complete absorption column has been constructed and has been partially commissioned. Initial data has been collected by operating using the fogging section and the frothing section. The fogging section produces solvent droplets of about 100 μm sauter mean diameter and as small as 25 μm using a hydraulic nozzle by BETE. The frothing section produces bubbles of about 5mm with high mixing of solvent promoting the higher mass transfer from gas to liquid. The absorber reaches the capture efficiency of about 50% with only two sections in operation. Based on the current results, it can be deduced that increasing the solvent feed temperature and including the packed section for absorption the capture efficiency will increase further. Initial data will be collected using all three sections of the absorber and will be used for sDOE. Non-Uniform Space Filling model of sDOE will be used to prioritize the input conditions resulting into maximum capture efficiency. sDOE is performed using the platform called Framework Optimization, Quantification of Uncertainty, and Surrogates (FOQUS). The method and results demonstrating the progress of the parametric campaign from the initial set of experiments to the final stage of obtaining optimized parameters using sDOE tool will be presented in detail.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

A Numerical Model for the Dynamics of Pickup Ions Outside the Heliopause and IBEX “Ribbon” Observation

The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

XRISM spectroscopy on orbital modulation of Fe Ly α lines in Cygnus X-3

To understand physical processes such as mass transfer and binary evolution in X-ray binaries, the orbital parameters of the system are fundamental and crucial information. Cygnus X-3 is a high-mass X-ray binary composed of a compact object of unknown nature and a Wolf–Rayet star, which is of great interest in the context of wind-fed mass accretion and binary evolution. Here we present XRISM/Resolve high-resolution spectroscopy, focusing on the Fe Ly$\alpha$ lines when the source was in its hypersoft state. We perform an orbital phase-resolved spectral analysis of the lines to study the orbital modulation of the emission and absorption lines. It is found that the emission lines reflect the orbital motion of the compact object whose estimated velocity amplitude is $430^{+150}_{-140}$ km s$^{-1}$, while the absorption lines show a variation that can be interpreted as originating from the stellar wind. We discuss possible mass ranges for the binary components using the mass function with the estimated value of the velocity amplitude in this work, combined with the relation between the mass-loss rate and the orbital period derivative and the empirical mass and mass-loss rate relation for Galactic Wolf–Rayet stars. They are constrained to be (1.3–5.1) and (9.3–12) M$_\odot$, respectively, for the assumed inclination angle of $i =25^\circ$, and become more relaxed, (1.3–24) and (9.3–16) M$_\odot$, for $i = 35^\circ$. Thus, it remains unclear whether the system harbors a black hole or a neutron star.

X-rays: binaries↗
DOE OSTI2026

ZTF SN Ia DR2 follow-up: Exploring the origin of the Type Ia supernova host galaxy step through Si II velocities

The relation between Type Ia supernovae (SNe Ia) and the stellar masses of their host galaxy is well documented. In particular, Hubble residuals display a distinct luminosity shift based on host mass. This is known as the mass step. This effect is widely used as an additional correction factor in the standardisation of SN Ia luminosities. We investigate the Hubble residuals and the mass step of normal SNe Ia in the context of Si IIλ6355 velocities based on 277 normal SNe Ia that are near their peak in the second data release (DR2) of the Zwicky Transient Facility (ZTF). We divided the sample into high-velocity (HV) and normal-velocity (NV) SNe Ia, separated at 12,000 km s −1 . This produced a sample of 70 HV and 207 NV objects. We then explored potential environment- and/or progenitor-related effects by investigating the Si IIλ6355 velocities with parameters such as the light-curve stretch x 1 , the colour c, and the host galaxy properties. Although we only find a marginal difference between the Hubble residuals of HV and NV SNe Ia, the NV mass step is 0.149 ± 0.024 mag (6.3σ). The HV mass step is smaller, 0.046 ± 0.041 mag (1.1σ), and is consistent with zero. The difference between the NV and HV mass steps is modest, at ∼2.2σ. Moreover, the clearest subtype difference appears for SNe in central regions (d DLR < 1), where NV SNe Ia show a large mass step, whereas HV SNe Ia are consistent with no step, yielding a difference of 3.1–3.6σ between NV and HV SNe Ia. We observe a host-colour step for both subtypes. NV SNe Ia show a step of 0.142 ± 0.024 mag (5.9σ), while HV SNe Ia show a step of 0.158 ± 0.042 mag (3.8σ), where the HV SNe Ia step appears to be larger, but the significance is lower because the sample size is smaller. Overall, the NV and HV colour steps are statistically consistent. HV SNe Ia also show modest (∼2.5–3σ) steps in certain subsets, such as those in outer regions (d DLR > 1), whereas NV SNe display stronger environmental trends. Our results indicate that NV SNe Ia appear to be more environmentally sensitive, particularly in central likely metal-rich and older regions, while HV SNe Ia show weaker and subset-dependent trends. This suggests that applying a universal mass-step correction might introduce biases, and that incorporating refined classifications and/or environment-dependent factors, such as the location within the host, might improve future cosmological analyses beyond the standard x 1 and c cuts.

supernovae: general↗
DOE OSTI2022

Dark energy survey year 3 results: cosmological constraints from the analysis of cosmic shear in harmonic space

ABSTRACT We present cosmological constraints from the analysis of angular power spectra of cosmic shear maps based on data from the first three years of observations by the Dark Energy Survey (DES Y3). Our measurements are based on the pseudo-Cℓ method and complement the analysis of the two-point correlation functions in real space, as the two estimators are known to compress and select Gaussian information in different ways, due to scale cuts. They may also be differently affected by systematic effects and theoretical uncertainties, making this analysis an important cross-check. Using the same fiducial Lambda cold dark matter model as in the DES Y3 real-space analysis, we find ${S_8 \equiv \sigma _8 \sqrt{\Omega _{\rm m}/0.3} = 0.793^{+0.038}_{-0.025}}$, which further improves to S8 = 0.784 ± 0.026 when including shear ratios. This result is within expected statistical fluctuations from the real-space constraint, and in agreement with DES Y3 analyses of non-Gaussian statistics, but favours a slightly higher value of S8, which reduces the tension with the Planck 2018 constraints from 2.3σ in the real space analysis to 1.5σ here. We explore less conservative intrinsic alignments models than the one adopted in our fiducial analysis, finding no clear preference for a more complex model. We also include small scales, using an increased Fourier mode cut-off up to $k_{\rm max}={5}\, {h}\, {\rm Mpc}^{-1}$, which allows to constrain baryonic feedback while leaving cosmological constraints essentially unchanged. Finally, we present an approximate reconstruction of the linear matter power spectrum at present time, found to be about 20 per cent lower than predicted by Planck 2018, as reflected by the lower S8 value.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

XRISM high-resolution X-ray spectroscopy of Cygnus X-1: Highly ionized iron absorption structures

We present the first high-resolution X-ray spectral analysis of Cygnus X-1 using XRISM. The observation wa3s carried out from 2024 April 7 to 10, covering the orbital phase range 0.65–0.17 during its low/hard state. Taking advantage of the exceptional energy resolution of the Resolve instrument, we examined highly ionized iron absorption lines and characterized the ionization states, column densities, and line-of-sight velocities of the absorbing plasma. Spectral analysis revealed an ionization parameter of $\xi \sim 3$, column densities of a few $\times 10^{21}$ cm$^{-2}$, and a blueshifted velocity of $\sim$100 km s$^{-1}$. The observation was divided into two phases: before and after orbital phase $\phi _{\rm {orb}} = 0.9$, corresponding to non-dipping and dipping intervals. While only weak absorption features were present before $\phi _{\rm {orb}} = 0.9$, strong absorption by He-like and H-like Fe appeared during the dipping phase. We measured equivalent widths of 2.3, 0.4, and 1.2 eV for He-like Fe K$\alpha$ and H-like Ly$\alpha _1$ and Ly$\alpha _2$, respectively—demonstrating the capability of XRISM Resolve to securely detect narrow absorption features of only a few eV. These measurements trace the motion of the absorbing material and offer insight into the kinematics and spatial distribution of the wind in the vicinity of the black hole. These findings enhance our understanding of wind-fed accretion in Cygnus X-1 and highlight the importance of continued high-resolution X-ray observations to further constrain the physical properties of winds and accretion flows in high-mass X-ray binaries.

Astronomy and AstroPhysics↗
DOE OSTI2025

How initial conditions-, structural-, and parameter-based model uncertainty interact and influence predictions in permafrost ecosystems: Modeling Archive

This dataset contains model output and input data, as well as source code examples for the Terrestrial Ecosystem Model with the Dynamic Vegetation Model and Dynamic Organic Soil (DVM-DOS-TEM) for the field sites Imnavait creek and the Bonanza creek Long Term Ecological Research Network (LTER). The data covers simulations from the last glacial maximum (LGM) until 2100 for a selection of paleo scenarios, setting the mean temperature of the LGM up to 10°C lower than pre-industrial conditions. The model structure was modulated to represent various model versions, and this dataset contains the relevant changes in the source code. The raw output data, the processed statistical data, the setup and processing scripts as well as parameter value distribution files from a parameter sensitivity analysis are included as well. Model outputs include active layer depth, organic soil carbon, soil layer depths, gross primary productivity (GPP) with and without nitrogen limitation, net primary productivity (NPP), soil liquid water content, heterotrophic, maintenance, and growth respiration, soil temperature, and vegetation carbon (*.nc files). The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research.Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Thermal and kinematic properties of ejecta in SN1987A revealed by XRISM

We present an analysis of high-resolution spectra from the shock-heated plasmas in SN 1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM). The 1.7–10 keV Resolve spectra are accurately represented by a single-component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58}$ × $10^{11}\,\,{\rm s\,\, cm}^{-3}$. The Resolve spectra are also well reproduced by the 3D magneto-hydrodynamic simulation presented by Orlando et al. (2020, A&A, 636, A22) suggesting substantial contribution from the ejecta. The metal abundances obtained with Resolve align with the Large Magellanic Cloud value, indicating that the X-rays in 2024 originate from “non-metal-rich” shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of ejecta. Doppler widths of the atomic lines from Si, S, and Fe correspond to velocities of 1500–1700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a $90\%$ upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV range, aligning with NuSTAR findings by Greco et al. (2022, ApJ, 931, 132). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4}\, M_{\odot }$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015, Science, 348, 670; Grebenev et al. 2012, Nature, 490, 373; Leising 2006, ApJ, 651, 1019).

ISM: supernova remnants↗
DOE OSTI2023

Polyhedral Distortions and Unusual Magnetic Order in Spinel FeMn 2 O 4

Spinel compounds AB 2 X 4 consist of both tetrahedral (AX 4 ) and octahedral (BX 6) environments with the former forming a diamond lattice and the latter a geometrically frustrated pyrochlore lattice. Exploring the fascinating physical properties and their correlations with structural features is critical in understanding these materials. FeMn 2 O 4 has been reported to exhibit one structural transition and two successive magnetic transitions. In this work, we report the polyhedral distortions and their correlations to the structural and two magnetic transitions in FeMn 2 O 4 by employing the high-resolution neutron powder diffraction. The cation distribution is found to be ($Mn_{0.9}^{2+}Fe_{0.1}^{3+}$)$_{A}$($Mn^{3+}Fe_{0.9}^{3+}Mn_{0.1}^{2+}$)$_{B}$O 4 . While large trigonal distortion is found even in the high-temperature cubic phase, the first-order cubic-tetragonal structural transition associated with the elongation of both tetrahedra and octahedra with shared oxygen atoms along the c axis occurs at T S ≈ 750 K, driven by the Jahn–Teller effect of the orbital active B-site Mn 3+ cation. Strong magnetoelastic coupling is unveiled at T N1 ≈ 400 K as manifested by the appearance of Néel-type collinear ferrimagnetic order, an anomaly in both tetrahedral and octahedral distortions, as well as an anomalous decrease of the lattice constants c and a weak anomaly of a. Upon cooling to T N2 ≈ 65 K, it evolves to a noncollinear ferrimagnetic order accompanied by the different moments at the split magnetic sites B1 and B2. Only one-half of the B-site Mn 3+ /Fe 3+ spins, i.e., the B2-site spins in the pyrochlore lattice, are canted, which is a unique magnetic order among spinels. The canting angle between A-site and B2-site moments is ~25°, but the B1-site moment stays antiparallel to the A-site moment even at 10 K. This noncollinear order is accompanied by a modification of the O–B–O bond angles in the octahedra without significant change in lattice constants or tetrahedral/octahedral distortion parameters, indicating a distinct magnetoelastic coupling. We demonstrate distinct roles of the A-site and B-site magnetic cations in the structural and magnetic properties of FeMn 2 O 4 . Our study indicates that FeMn 2 O 4 is a wonderful platform to unveil interesting magnetic order and to investigate their correlations with polyhedral distortions and lattice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗