Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flux Reconstruction”

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 289 records · Page 16

A high-order finite difference method for moving immersed domain boundaries and material interfaces

Here, we present a high-order sharp treatment of immersed moving domain boundaries and material interfaces, and apply it to the advection-diffusion equation in two and three dimensions. The spatial discretization combines dimension-split finite difference schemes with an immersed boundary treatment based on a weighted least-squares reconstruction of the solution, providing stable discretizations with up to sixth order accuracy for diffusion terms and third order accuracy for advection terms. The temporal discretization relies on a novel strategy for maintaining high-order temporal accuracy in problems with moving boundaries that minimizes implementation complexity and allows arbitrary explicit or diagonally-implicit Runge-Kutta schemes. The approach is broadly compatible with popular PDE-specialized Runge-Kutta time integrators, including low-storage, strong stability preserving, and diagonally implicit schemes. Through numerical experiments we demonstrate that the full discretization maintains high-order spatial and temporal accuracy in the presence of complex 3D geometries and for a range of boundary conditions, including Dirichlet, Neumann, and flux conditions with large jumps in coefficients.

97 MATHEMATICS AND COMPUTING↗

Time-resolved 3D imaging of two-phase fluid flow inside a steel fuel injector using synchrotron X-ray tomography

The multiphase flow inside a diesel injection nozzle is imaged using synchrotron X-rays from the Advanced Photon Source at Argonne National Laboratory. Through acquisitions performed at several viewing angles and subsequent tomographic reconstruction, in-situ 3D visualization is achieved for the first time inside a steel injector at engine-like operating conditions. The morphology of the internal flow reveals strong flow separation and vapor-filled cavities (cavitation), the degree of which correlates with the nozzle's asymmetric inlet corner profile. Micron-scale surface features, which are artifacts of manufacturing, are shown to influence the morphology of the resulting liquid-gas interface. The data obtained at 0.1 ms time resolution exposes transient flow features and the flow development timescales are shown to be correlated with in-situ imaging of the fuel injector's hydraulically-actuated valve (needle). As more than 98.5% of the X-ray photon flux is attenuated within the steel injector body itself, we are posed with a unique challenge for imaging the flow within. Time-resolved imaging under these low-light conditions is achieved by exploiting both the refractive and absorptive properties of X-ray photons. The data-processing strategy converted these images with a signal-to-noise ratio of similar to 10 into a meaningful dataset for understanding internal flow and cavitation in a nozzle of diameter 200 mu m enclosed within 1-2 millimeters of steel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A convergence metric for counting statistics in time-resolved small angle neutron scattering

Here, this work introduces a model-independent, dimensionless metric for predicting optimal measurement duration in time-resolved small-angle neutron scattering using early-time data. Built on a Gaussian process regression framework, the method reconstructs scattering profiles with quantified uncertainty, even from sparse or noisy measurements. Demonstrated on the EQ-SANS instrument at the Spallation Neutron Source, the approach generalizes to general SANS instruments with a two-dimensional detector. A key result is the discovery of a dimensionless convergence metric revealing a universal power-law scaling in profile evolution across soft matter systems. When time is normalized by a system-specific characteristic time t*, the variation in inferred profiles collapses onto a single curve with an exponent between −2 and −1. This trend emerges within the first ten time steps, enabling early prediction of measurement sufficiency. The method supports real-time experimental optimization and is especially valuable for maximizing efficiency in low-flux environments such as compact accelerator-based neutron sources.

Tung, Chi-Huan [Oak Ridge National Laboratory (ORN↗

Constraining the Milky Way’s ultraviolet-to-infrared SED with Gaussian process regression

ABSTRACT Improving our knowledge of global Milky Way (MW) properties is critical for connecting the detailed measurements only possible from within our Galaxy to our understanding of the broader galaxy population. We here train Gaussian process regression (GPR) models on SDSS (Sloan Digital Sky Survey) galaxies to map from galaxy properties (stellar mass, apparent axial ratio, star formation rate, bulge-to-total ratio, disc scale length, and bar vote fraction) to ultraviolet (UV; GALEX FUV/NUV), optical (SDSS ugriz), and infrared (IR; 2MASS JHKs and WISE W1/W2/W3/W4) fluxes and uncertainties. With these models, we estimate the photometric properties of the MW, resulting in a full UV-to-IR spectral energy distribution (SED) as it would be measured externally, viewed face-on. We confirm that the MW lies in the green valley in optical diagnostic diagrams, but show for the first time that the MW is in the star-forming region in standard UV and IR diagnostics – characteristic of the population of red spiral galaxies. Although our GPR method predicts one band at a time, the resulting MW UV–IR SED is consistent with SEDs of local spirals with characteristics broadly similar to the MW, suggesting that these independent predictions can be combined reliably. Our UV–IR SED will be invaluable for reconstructing the MW’s star formation history using the same tools employed for external galaxies, allowing comparisons of results from in situ measurements to those from the methods used for extragalactic objects.

79 ASTRONOMY AND ASTROPHYSICS↗

A Satellite-Derived Upper-Ocean Stratification Data Set for the Tropical North Atlantic With Potential Applications for Hurricane Intensity Prediction

Upper-ocean stratification strongly impacts vertical mixing and the heat flux between the ocean and atmosphere, especially under extreme conditions of tropical cyclones (TCs). Knowledge of prestorm stratification is important for accurate TC intensity prediction. In situ observations of the tropical ocean have significantly increased in the past decade. However, they are still too sparse to resolve ocean stratification variability in near-real time and on small spatial scales. Here, in this study, based on long-term observations and an ocean reanalysis data set from 2004–2017, we investigate the possibility of retrieving upper-ocean stratification from sea surface temperature (SST), sea surface salinity (SSS), and sea surface height (SSH) using a simple regression method. It is found that more than 90% of the mean seasonal cycle and about 30% to 80% of temperature and salinity stratification anomalies can be reconstructed using surface data from either observations or an ocean reanalysis. Simple regression can be used with satellite observations to create a high-resolution, near-real-time-gridded ocean stratification data set that successfully reproduces both the large and mesoscale variability of ocean stratification. When used in a simple expression for TC-induced SST cooling, the satellite-derived stratification shows improvements over an ocean analysis in terms of variance explained of SST cooling, offering promise as a near-real-time indicator of the ocean's impact on TC intensification.

54 ENVIRONMENTAL SCIENCES↗

Analysis of the inter-species power balance in JET plasmas

Most auxiliary heating methods provide heating to more than one particle species (electrons, ions, impurities) in a fusion plasma. This can lead to substantial temperature differences between species, depending on conditions such as heating power to the different species and collisionality, with temperature differences between species limited by inter-species thermal equipartition and transport. The analysis of the steady-state electron-ion and impurity-ion power balances presented in this paper are used for consistency-checking experimental ion and electron temperature measurements and for inferring the main ion temperature from measured impurity temperatures. As ion temperature measurements by charge exchange spectroscopy (CXS) based on impurity ions have become more difficult and time-consuming since the installation of the ITER-like wall (ILW) with Be and W PFC's, knowing the maximum sustainable temperature difference between ions and electrons, |T i -T e | allows rejecting erroneous measurements. It also obviates the need for an ion temperature measurement, if an electron temperature measurement is available and |T i - T e | cannot be larger than the combined errors of the underlying measurements. A power balance analysis is also required for estimating the errors of the ion and electron heat fluxes prior to any species-resolved transport analysis. The ion-impurity temperature differences are usually found to be small due to strong thermal equipartition between ion species. However, they can approach 10% in JET-ILW low density, high power discharges, such the ones under development for a future JET deuterium–tritium campaign (Joffrin et al 2019 Nucl. Fusion). This has a generally small, but not always negligible effect on the calculation of fusion reaction rates, which depend on main ion temperatures. An important outcome of this analysis is that temperature differences between impurity species are always much smaller than between the impurities and hydrogenic species and can usually be neglected. The article presents two methods for calculating the impurity-to-main ion temperature ratio. Finally, this analysis leads to a method for the reconstruction ion temperature profiles from ion temperature data available at only one or a small number of spatial locations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of the COMPASS results *

COMPASS addressed several physical processes that may explain the behaviour of important phenomena. This paper presents results related to the main fields of COMPASS research obtained in the recent two years, including studies of turbulence, L–H transition, plasma material interaction, runaway electron, and disruption physics: Tomographic reconstruction of the edge/SOL turbulence observed by a fast visible camera allowed to visualize turbulent structures without perturbing the plasma. Dependence of the power threshold on the X-point height was studied and related role of radial electric field in the edge/SOL plasma was identified. The effect of high-field-side error fields on the L–H transition was investigated in order to assess the influence of the central solenoid misalignment and the possibility to compensate these error fields by low-field-side coils. Results of fast measurements of electron temperature during ELMs show the ELM peak values at the divertor are around 80% of the initial temperature at the pedestal. Liquid metals were used for the first time as plasma facing material in ELMy H-mode in the tokamak divertor. Good power handling capability was observed for heat fluxes up to 12 MW m –2 and no direct droplet ejection was observed. Partial detachment regime was achieved by impurity seeding in the divertor. The evolution of the heat flux footprint at the outer target was studied. Runaway electrons were studied using new unique systems—impact calorimetry, carbon pellet injection technique, wide variety of magnetic perturbations. Radial feedback control was imposed on the beam. Forces during plasma disruptions were monitored by a number of new diagnostics for vacuum vessel (VV) motion in order to contribute to the scaling laws of sideways disruption forces for ITER. Current flows towards the divertor tiles, incl. possible short-circuiting through PFCs, were investigated during the VDE experiments. The results support ATEC model and improve understanding of disruption loads.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Lyman-$α$ forest holography: 3D predictions from 1D measurements

Cosmological analyses of Lyman-$α$ forest clustering rely on either one-dimensional correlations along individual sightlines or three-dimensional correlations between different sightlines. Because these observables probe the matter distribution on very different scales, they have traditionally been analyzed independently. In this work, we bridge this gap using ForestFlow, an emulator trained on a suite of cosmological hydrodynamical simulations that provides a unified description of Lyman-$α$ forest clustering from linear to nonlinear scales. This framework enables us to determine the range of three-dimensional clustering models compatible with the DESI one-dimensional flux power spectrum ($P_{\rm 1D}$). The resulting predictions successfully reproduce the large-scale clustering measured by the DESI BAO analysis and provide physically motivated priors on nonlinear clustering that are used in a companion paper presenting the full-shape analysis of the DESI DR2 Lyman-$α$ forest. We validate our methodology using the large-volume, high-resolution hydrodynamical simulation ACCEL-2, demonstrating excellent agreement across the full range of scales considered. Finally, we combine constraints from the $P_{\rm 1D}$ and BAO analyses on the parameter combinations $b_δσ_8$ and $b_ηf σ_8$, finding that the two probes provide comparable constraining power while exhibiting complementary parameter degeneracies. Our results establish a direct connection between one- and three-dimensional Lyman-$α$ forest measurements through ForestFlow, an approach we term Lyman-$α$ holography by analogy with the reconstruction of higher-dimensional structure from lower-dimensional information.

Chaves-Montero, J. [Barcelona, IFAE] (ORCID:000000↗

Performance evaluation of cosmic ray muon trajectory estimation algorithms

Muons, being elementary particles with minimal interaction with nuclear materials and abundant at sea level, have sparked interest in utilizing them for imaging various applications, such as mining [Borselli et al., Sci. Rep. 12, 22329 (2022)], volcano imaging [Nagamine et al., Nucl. Instrum. Meth. A, 356, 585(1995)], and underground tunnel detection [Guardincerri et al., Pure Appl. Geophys. 174, 2133 (2017)]. Recently, their use in nuclear nonproliferation and safeguard verification has gained attention, particularly in cargo screening for nuclear waste smuggling [Baesso et al., J. Instrum. 9, C10041 (2014)], source localization [L. J. Schultz et al., Nucl. Instrum. Meth. A 519, 687 (2004)], and locating nuclear fuel debris in reactors [Borozdin et al., Phys. Rev. Let. 109, 152501 (2012)]. However, the resolution of muon image reconstruction techniques is limited due to multiple Coulomb scattering (MCS) within the target object. To achieve robust muon tomography, it is crucial to develop efficient and flexible physics-based algorithms that can model the MCS process accurately and estimate the most probable trajectory of muons as they pass through the target object. To address this limitation, in this study, a novel algorithmic approach utilizing the Bayesian probability theory and Gaussian approximation of MCS is chosen. Different energy levels, materials, and target sizes were considered in the evaluations. The results demonstrate that the Generalized Muon Trajectory Estimation (GMTE) algorithm offers significant improvements over currently used algorithms. Across all test scenarios, the GMTE algorithm demonstrated ~50% and 38% increase in precision compared to Straight Line Path (SLP) and Point of Closest Approach (PoCA) algorithms, respectively. Furthermore, it exhibited 10%–35% and 10%–15% increases in muon flux utilization for high and medium Z materials, respectively, compared to the PoCA algorithm. In conclusion, the extensive simulations confirm the enhanced performance and efficiency of the GMTE algorithm, offering improved resolution and reduced measurement time for cosmic ray muon imaging compared to the current SLP and PoCA algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗

Planck intermediate results

In this work, we describe an extension of the most recent version of the Planck Catalogue of Compact Sources (PCCS2), produced using a new multi-band Bayesian Extraction and Estimation Package (BeeP). BeeP assumes that the compact sources present in PCCS2 at 857 GHz have a dust-like spectral energy distribution (SED), which leads to emission at both lower and higher frequencies, and adjusts the parameters of the source and its SED to fit the emission observed in Planck’s three highest frequency channels at 353, 545, and 857 GHz, as well as the IRIS map at 3000 GHz. In order to reduce confusion regarding diffuse cirrus emission, BeeP’s data model includes a description of the background emission surrounding each source, and it adjusts the confidence in the source parameter extraction based on the statistical properties of the spatial distribution of the background emission. BeeP produces the following three new sets of parameters for each source: (a) fits to a modified blackbody (MBB) thermal emission model of the source; (b) SED-independent source flux densities at each frequency considered; and (c) fits to an MBB model of the background in which the source is embedded. BeeP also calculates, for each source, a reliability parameter, which takes into account confusion due to the surrounding cirrus. This parameter can be used to extract sub-samples of high-frequency sources with statistically well-understood properties. We define a high-reliability subset (BeeP/base), containing 26 083 sources (54.1% of the total PCCS2 catalogue), the majority of which have no information on reliability in the PCCS2. We describe the characteristics of this specific high-quality subset of PCCS2 and its validation against other data sets, specifically for: the sub-sample of PCCS2 located in low-cirrus areas; the Planck Catalogue of Galactic Cold Clumps; the Herschel GAMA15-field catalogue; and the temperature- and spectral-index-reconstructed dust maps obtained with Planck’s Generalized Needlet Internal Linear Combination method. The results of the BeeP extension of PCCS2, which are made publicly available via the Planck Legacy Archive, will enable the study of the thermal properties of well-defined samples of compact Galactic and extragalactic dusty sources.

79 ASTRONOMY AND ASTROPHYSICS↗

Spectral characterization of flash and high flux x-ray radiographic sources with a magnetic Compton spectrometer

In this article, we present a new analysis method applied to revitalize permanent magnet Compton spectrometers used to measure photon energy spectra in the MeV range. The inversion of the measured electron distribution to determine the original photon distribution is achieved via a method of consistent coupled radiation transport and magnetic field mapping of the input photon spectra to the measured electron distribution. The method of linear least squares was used to perform the unfolding of the electron distribution to the initial photon spectra, without any assumptions made regarding the electron distribution. We present an application of this method to data from a nominal 19.4 MeV flash radiographic source (the first axis of the Dual Axis Radiographic Hydro-Test Facility) capable of generating 500 R @ 1 m in ~60 ns and a medical therapy source (a Scanditronix M22, Microtron) capable of variable energies with nominal endpoints of 6, 10, 15, and 20 MeV and an output of ~1000–2000 R/min @ 1 m. The results provide agreement between the modeled and unfolded experimentally measured photon spectra as quantified by statistical tests, from 1.5 to 20 MeV. Experimental results are presented as well as a discussion of the novel MCNP6-based simulations and methods for reconstruction of the spectra.

47 OTHER INSTRUMENTATION↗

Application of recoil-imaging time projection chambers to directional neutron background measurements in the SuperKEKB accelerator tunnel

Gaseous time projection chambers (TPCs) with high readout segmentation are capable of reconstructing detailed 3D ionization distributions of nuclear recoils resulting from elastic neutron scattering. Here, using a system of six compact TPCs with pixel ASIC readout, filled with a 70:30 mixture of He:CO 2 gas, we analyze the first directional measurements of beam-induced neutron backgrounds in the tunnel regions surrounding the Belle II detector at the SuperKEKB e + e – collider. With the use of 3D recoil tracking, we show that these TPCs are capable of maintaining nearly 100% nuclear recoil purity to reconstructed ionization energies (E reco ) as low as 5 keV ee . Using a large sample of Monte-Carlo (MC)-simulated 4 He, 12 C, and 16 O recoil tracks, we find consistency between predicted and measured recoil energy spectra in five of the six TPCs, providing useful validation of the neutron production mechanisms modeled in simulation. Restricting this sample to 4 He recoil tracks with E reco > 40 keV ee , we further demonstrate axial angular resolutions within 8° and we introduce a procedure that under suitable conditions, correctly assigns the vector direction to 91% of these simulated 4 He recoils. Applying this procedure to assign vector directions to measured 4 He recoil tracks, we observe consistency between the angular distributions of observed and simulated recoils, providing first experimental evidence of localized neutron “hotspots” in the accelerator tunnel. Observed rates of nuclear recoils in these TPCs suggest that simulation overestimates the neutron flux from these hotspots. Despite this, we estimate these hotspots to produce the majority of neutron backgrounds in the accelerator tunnel at SuperKEKB’s target luminosity of 6.3 x 10 35 cm –2 s –1 , making them important regions to continue to monitor.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Λ Baryon Production in ν¯µ Interactions in the MicroBooNE Detector

The Cabibbo suppressed production of $\Lambda$ baryons in anti-neutrino interactions with nuclei is a rare process that is yet to be measured with a modern neutrino detector with automated reconstruction. The cross section for this process is sensitive to a number of unique nuclear effects, most notably the secondary interactions of the produced hyperon while attempting to escape from the nucleus. Other interactions within the nuclear remnant can impact the estimation of neutrino energy in oscillation measurements, and thus an accurate description of the nuclear environment is required. The strangeness violating hyperon production process is only available to anti-neutrinos. The model of this interaction is implemented into the NuWro neutrino interaction Monte Carlo simulation, and some predictions are presented, focusing on the role of nuclear effects. This model introduces a hyperon-nucleus potential, which calculations from hypernuclear theory permit to be strongly repulsive in th e case of $\Sigma$ baryons. The presence of this potential is found to sculpt the shape of the differential cross section in some variables. The MicroBooNE detector will be described, followed by a description of a measurement of the flux averaged, restricted phase space cross section of Cabibbo suppressed $\Lambda$ baryon production. A sophisticated event selection is employed, as a very large quantity of background neutrino interactions must be removed to perform the measurement with any sensitivity. This selection introduces some novel techniques such as the island finding method, and achieves a background reduction of $\sim 10^6$, with an efficiency of around 7\%. The calculation of the systematic uncertainties will be explained, including two procedures explored to handle sources of background with extremely poor simulation statistics: an in-situ constraint using data from sidebands, and a visual inspection of the data and simulation to remove the troublesome background events. The sensitivity to the $\Lambda$ baryon production cross section is calculated in the form of Bayesian posterior probability distributions, combining the systematic uncertainties with data and simulation statistical uncertainties. As a rare process, the statistical uncertainties are highly non-Gaussian, and the Bayesian approach is applied to include the full shapes of these uncertainties. Data corresponding to $2.2 \times 10^{20}$ protons on target of neutrino mode running and $4.9 \times 10^{20}$ protons on target of anti-neutrino running is analysed. When the data was unblinded, five $\Lambda$ production candidates were selected from the data, consistent with the MC simulation prediction of $5.3 \pm 1.1$ events. The final estimated cross section is $1.8^{+2.0}_{-1.6} \times 10^{-40}$cm$^2/$Ar when employing the sideband constraint procedure. A similar result of $2.0^{+2.2}_{-1.8} \times 10^{-40}$cm$^2/$Ar is obtained when performing the visual scan instead. The methods used in t his analysis are intended to be easily exported to other LArTPC detectors such as the Short Baseline Near Detector.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A scalable compressible volume of fluid solver using a stratified flow model

Here, in this paper, a Volume of Fluid method which uses a Stratified Flow model for flux calculation is proposed to solve a compressible multiphase flow in a way that has high parallel efficiency. The solver is part of an open-source computing tool to study fuel entrainment and combustion in the application of paraffin-based hybrid rocket motors. This paper focuses on the compressible Volume of Fluid solver for the liquid and gaseous phases in this environment. This solver uses perfect gas and stiffened gas models for equation of state, with the ability to easily add other models as desired. Additionally a comparison of a Riemann solver versus an AUSM+up scheme is conducted, showing improvement with the later in the solution of multiphase flows. This solver successfully reproduces shock tube cases with the expected accuracy, as well as simulating two-dimensional shear and gravity-driven flows. Results also show that the code is able to utilize a Stratified Flow model to evaluate shear flow without the need for interface reconstruction or gradient calculations. Finally, it is shown that the solver has near-ideal scaling under strong scaling tests as well as good performance in static scaling, giving improved performance over current options and implying future advancement of high-performance computing efficiency for multiphase flow solvers.

97 MATHEMATICS AND COMPUTING↗

Nonlocal, Pattern-Aware Response and Feedback Framework for Regional Climate Response

We devise a pattern-aware feedback framework for representing the forced climate response using a suite of Green’s function experiments with solar radiation perturbations. By considering the column energy balance, a comprehensive linear response function (CLRF) for important climate variables and feedback quantities such as moist static energy, sea surface temperature, albedo, cloud optical depth, and lapse rate is learned from Green’s function data. The learned CLRF delineates the effects of the energy diffusion in both the ocean and atmosphere and the pattern-aware feedbacks from the aforementioned radiatively active processes. The CLRF can then be decomposed into forcing–response mode pairs, which are in turn used to construct a reduced-order model describing the dominant dynamics of climate responses. These mode pairs capture nonlocal effects and teleconnections in the climate and thus make the reduced-order model apt for capturing regional features of climate response. A key observation is that the CLRF captures the polar-amplified response as the most excitable mode of the climate system, and this mode is explainable in the data-learned pattern-aware feedback framework. The reduced-order model can be used for predicting the response for a given forcing and for reconstructing the forcing from a given response; we demonstrate these capabilities for multiple independent forcing scenarios.

Feedback↗

An automated fast neutron computed tomography instrument with on-line focusing for non-destructive evaluation

A fast neutron tomography imaging instrument has been designed, built, and tested at The Ohio State University 500 kW Research Reactor on a fast neutron beamline with a peak neutron flux ≈5.4 × 107 n·cm−2·s−1 at 1.6 MeV median neutron energy. The instrument and beamline are also configurable for thermal neutron imaging. The imaging apparatus is composed of a lens coupled, water-cooled Electron Multiplying Charge Coupled Device camera, a front-surface mirror, and a high light yield plastic Polyvinyl toluene scintillator. The instrument sits on a mobile cart. A total of 5 motion-control stages are built into the system for XYZ and rotational degrees of freedom for sample positioning; the fifth stage fine tunes the focal distance between the camera and the scintillator to achieve on-line focusing. A Python code with a user-friendly graphical user interface controls the fully automated image acquisition, not requiring user interaction, yet facilitating tracking of the image acquisition. A complete fast neutron computed tomography dataset with 360 projections requires less than 3 h, with 30 s per projection. On-line focusing is accomplished with a commercial, off-the-shelf, dielectrically actuated liquid lens. Finally, tomographic reconstructions are visualized using the Livermore Tomography Tools software package. The effective pixel size (width and height) is ≈0.1058 mm, yielding a minimum voxel size of 0.1058 × 0.1058 × 0.1058 mm3, and produces a spatial resolution of 231 μm when calculated from knife-edge measurements.

Bisbee, M. G. (ORCID:0000000313466697)↗

How to Build an Empirical Speed Distribution for Dark Matter in the Solar Neighborhood

The dark matter flux in a direct detection experiment depends on its local speed distribution. This distribution has been inferred from simulations of Milky Way–like galaxies, but such models serve only as proxies, given that no simulation directly captures the detailed evolution of our own Galaxy. This motivates alternative approaches that obtain this distribution directly from observations. In this work, we utilize 98 Milky Way analogues from the TNG50 simulation to develop and validate a procedure for inferring the dark matter speed distribution using the kinematics of nearby stars. We find that the dark matter that originated from old mergers, plus that from recent nonluminous accretions, is well described by a Maxwell–Boltzmann speed distribution centered at the local standard-of-rest velocity. Meanwhile, recently accreted dark matter from massive mergers has speeds that can be traced from the associated stellar debris of these events. The stellar populations systematically underestimate the velocity dispersion of their dark matter counterparts, but a simple kinematic boost brings the two into good alignment. Using the TNG50 host galaxies, we demonstrate that combining these two contributions provides an accurate reconstruction of the local dark matter speeds. As an application of the procedure to our own Galaxy, we utilize stellar kinematic data from Gaia to quantify how the dark matter remnants from the Milky Way’s last major merger impact its speed distribution in the solar neighborhood.

Shpigel, Tal [Princeton Univ., NJ (United States)]↗

Double-differential Measurements of Mesonless Charged-current Muon Neutrino Interactions on Argon with Final-State Protons using the MicroBooNE Detector

This note reports the first double-differential measurements of charged-current $ν_µ$ scattering on argon leading to final states containing zero mesons and one or more protons. This event topology (hereafter abbreviated as CC0$πNp$) is the most common at the neutrino energies typically produced by the Fermilab Booster Neutrino Beam. A detailed understanding of neutrino-argon scattering in the CC0$πNp$ channel is therefore crucial for the success of the precision neutrino oscillation analyses planned for the Short-Baseline Neutrino (SBN) program. This remains true for the upcoming Deep Underground Neutrino Experiment (DUNE), but the higher mean neutrino energy used there will ensure that more inelastic reaction modes, such as single pion production, will also play a major role. The analysis described in this note builds on a previous MicroBooNE study of the CC0$πNp$ channel which obtained the first single-differential cross-section measurements on an argon target. Since that foundational work, significant improvements have been made to MicroBooNE’s simulation software, event reconstruction algorithms, and procedure for calculating systematic uncertainties. When combined with a larger dataset (corresponding to a beam exposure of 6.79 × 10 20 protons-on-target versus 1.60 × 10 20 in Ref. [3]), these enhancements allow the important CC0$πNp$ channel to be studied in more detail. This note begins with a description of the data and simulation samples used as input to the analysis. Section 3 then defines the CC0$πNp$ signal event topology, and Section 4 describes a set of selection criteria designed to identify these events in MicroBooNE data. Binning schemes are then defined in Section 5 for two double-differential measurements of event rates. The first of these considers the momentum and scattering cosine of the outgoing muon, while the second reports the same observables for the leading proton, i.e., the final-state proton with the largest momentum. After a discussion of systematic uncertainties in Section 6, the note concludes by comparing the predictions of MicroBooNE Monte Carlo (MC) simulations to the measured double-differential distributions. These results will form the basis for a future extraction of flux-averaged double-differential CC0$πNp$ cross sections that will be immediately comparable to the theoretical predictions of multiple neutrino event generators. The selection described herein may also be used to study various other observables in CC0$πNp$ events, including those which are sensitive to correlations between leptonic and hadronic kinematics in the final state.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗