Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “peak finding”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Responsiveness of miscanthus and switchgrass yields to stand age and nitrogen fertilization: A meta‐regression analysis

Abstract Optimal management of the perennial bioenergy crops, miscanthus and switchgrass, requires an understanding of their responsiveness to nitrogen (N) fertilizer at different maturity stages across locations and growing conditions. Earlier studies that have examined the yield response of these crops to N and stand age using field experiments or meta‐analysis techniques provide mixed evidence. We extend earlier studies by applying a multi‐level mixed‐effects (MLME) meta‐regression model to conduct a more extensive multivariate regression of yield response of these crops to N and stand age, while controlling for climate and location conditions and unobserved factors related to study design. Our findings are based on 1403 and 2811 yield observations for miscanthus and switchgrass, respectively, from experiments conducted between 2002 and 2019 across the rainfed region in the United States. We find statistically significant evidence that an additional year of maturity increases miscanthus and switchgrass yields but at a decreasing rate; yields peak at the 7th and 6th year respectively, for the observed range of applied N rates and stands. We also find that an increase in N application increases yield by a statistically significant level, but at a declining rate; the magnitude of the yield response to N is, however, small and varies with the age of the crop. The impact of N is larger on older compared to younger and middle‐aged stands of miscanthus. In contrast, the impact of N on switchgrass is larger on middle‐aged compared to younger and older stands of switchgrass. We do not find a statistically significant effect of soil productivity on yield for either crop. This analysis provides a basis for developing N application recommendations and optimal rotation age for miscanthus and switchgrass and shows that these energy crops can grow just as productively on low productivity land as on high productivity land.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Signature of a Massive Rotating Metal-poor Star Imprinted in the Phoenix Stellar Stream

The Phoenix stellar stream has a low intrinsic dispersion in velocity and metallicity that implies the progenitor was probably a low-mass globular cluster. In this work we use Magellan/Magellan Inamori Kyocera Echelle (MIKE) high-dispersion spectroscopy of eight Phoenix stream red giants to confirm this scenario. In particular, we find negligible intrinsic scatter in metallicity ($σ$([FE II/H]) = ${0.04}_{-0.03}^{+0.11}$) and a large peak-to-peak range in [Na/Fe] and [Al/Fe] abundance ratios, consistent with the light element abundance patterns seen in the most metal-poor globular clusters. However, unlike any other globular cluster, we also find an intrinsic spread in [Sr II/Fe] spanning ~1 dex, while [Ba II/Fe] shows nearly no intrinsic spread ($σ$([Ba II/H]) = ${0.03}_{-0.02}^{+0.10}$). This abundance signature is best interpreted as slow-neutron-capture element production from a massive fast-rotating metal-poor star (15–20 $M$ ⊙ , $v$ ini /$v$ crit = 0.4, [Fe/H] = -3.8). The low inferred cluster mass suggests the system would have been unable to retain supernovae ejecta, implying that any massive fast-rotating metal-poor star that enriched the interstellar medium must have formed and evolved before the globular cluster formed. Furthermore, neutron-capture element production from asymptotic giant branch stars or magneto-rotational instabilities in core-collapse supernovae provide poor fits to the observations. We also report one Phoenix stream star to be a lithium-rich giant ($A$(Li) = 3.1 ± 0.1). At [Fe/H ] = -2.93; it is among the most metal-poor lithium-rich giants known.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

MARS-Q modeling of kink-peeling instabilities in DIII-D QH-mode plasma

In quiescent H-mode (QH-mode) regime, edge harmonic oscillations (EHOs) are believed to provide necessary radial transport to prevent occurrence of large edge localized modes. A systematic modeling study is performed here on the low- n EHOs in a DIII-D QH-mode plasma (Chen et al 2016 Nucl. Fusion 56 076011), by utilizing the MARS-Q code (Liu et al 2013 Phys. Plasmas 20 042503). Both the n = 1 and n = 2 instabilities are found to be strongly localized near the plasma edge, exhibiting the edge-peeling characteristics. The DIII-D resistive wall is found to have minor effects on these instabilities. The plasma resistivity is found to strongly modify the mode growth rate. Assuming the Spitzer model for the plasma resistivity, the computed mode growth rate scales as S sup>-1/3 with S being the Lundquist number. Toroidal flow of the plasma slightly stabilizes these edge localized kink-peeling modes. Drift kinetic effects all have a destabilization effect on these modes. Non-perturbative magneto-hydrodynamic-kinetic hybrid computations find that the drift kinetic effects associated with thermal particle species push the peak location of the eigenmode radially inward but still in the pedestal region. The modeled plasma temperature and density fluctuations in the plasma edge region, as well as the poloidal magnetic field perturbations along both the low and high field sides of the plasma surface, are in good agreement with experimental measurements. Finally, the quasi-linear initial value simulations find a strong non-linear interplay between the kink-peeling instability and the toroidal flow near the plasma edge. The combined effect of the damping of the flow amplitude and change of the edge flow shear is found to be the stabilizing factor for the kink-peeling mode, leading to the mode saturation and thus EHOs.

Physics↗
DOE OSTI2024

Cold-Season Precipitation Sensitivity to Microphysical Parameterizations: Hydrologic Evaluations Leveraging Snow Lidar Datasets

Abstract Cloud microphysical processes are an important facet of atmospheric modeling, as they can control the initiation and rates of snowfall. Thus, parameterizations of these processes have important implications for modeling seasonal snow accumulation. We conduct experiments with the Weather Research and Forecasting (WRF V4.3.3) Model using three different microphysics parameterizations, including a sophisticated new scheme (ISHMAEL). Simulations are conducted for two cold seasons (2018 and 2019) centered on the Colorado Rockies’ ∼750-km 2 East River watershed. Precipitation efficiencies are quantified using a drying-ratio mass budget approach and point evaluations are performed against three NRCS SNOTEL stations. Precipitation and meteorological outputs from each are used to force a land surface model (Noah-MP) so that peak snow accumulation can be compared against airborne snow lidar products. We find that microphysical parameterization choice alone has a modest impact on total precipitation on the order of ±3% watershed-wide, and as high as 15% for certain regions, similar to other studies comparing the same parameterizations. Precipitation biases evaluated against SNOTEL are 15% ± 13%. WRF Noah-MP configurations produced snow water equivalents with good correlations with airborne lidar products at a 1-km spatial resolution: Pearson’s r values of 0.9, RMSEs between 8 and 17 cm, and percent biases of 3%–15%. Noah-MP with precipitation from the PRISM geostatistical precipitation product leads to a peak SWE underestimation of 32% in both years examined, and a weaker spatial correlation than the WRF configurations. We fall short of identifying a clearly superior microphysical parameterization but conclude that snow lidar is a valuable nontraditional indicator of model performance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

MUSEQuBES: calibrating the redshifts of Ly α emitters using stacked circumgalactic medium absorption profiles

ABSTRACT Ly α emission lines are typically found to be redshifted with respect to the systemic redshifts of galaxies, likely due to resonant scattering of Ly α photons. Here, we measure the average velocity offset for a sample of 96 z ≈ 3.3 Ly α emitters (LAEs) with a median Ly α flux (luminosity) of ${\approx}10^{-17}~\rm erg~cm^{-2}~s^{-1}$ (${\approx}10^{42}~\rm erg~s^{-1}$) and a median star formation rate (SFR) of ≈1.3 $\rm M_{\odot }~yr^{-1}$ (not corrected for possible dust extinction), detected by the Multi-Unit Spectroscopic Explorer as part of our MUSEQuBES circumgalactic medium (CGM) survey. By postulating that the stacked CGM absorption profiles of these LAEs, probed by eight background quasars, must be centred on the systemic redshift, we measure an average velocity offset, Voffset = 171\pm 8 km s−1, between the Ly α emission peak and the systemic redshift. The observed Voffset is lower by factors of ≈1.4 and ≈2.6 compared to the velocity offsets measured for narrow-band-selected LAEs and Lyman break galaxies, respectively, which probe galaxies with higher masses and SFRs. Consistent with earlier studies based on direct measurements for individual objects, we find that the Voffset is correlated with the full width at half-maximum of the red peak of the Ly α line, and anticorrelated with the rest-frame equivalent width. Moreover, we find that Voffset is correlated with SFR with a sub-linear scaling relation, $V_{\rm offset}\propto \rm SFR^{0.16\pm 0.03}$. Adopting the mass scaling for main-sequence galaxies, such a relation suggests that Voffset scales with the circular velocity of the dark matter haloes hosting the LAEs.

Muzahid, Sowgat↗
DOE OSTI2026

Characterization of Precipitation-Induced Radon Progeny Deposition Events Using a City-Scale Sensor Network

Networks of radiation detectors provide a platform for real-time radioactive source detection and identification in urban environments. Detection algorithms in these systems must adapt to naturally-occurring changes in background, which requires well-characterized relationships between precipitation events and their corresponding radiological signature. Here, we present a quantitative and qualitative description of rain-induced radon progeny deposition events occurring in Chicago from September 2023 to February 2024. We measure ambient gamma radiation levels, precipitation rate, temperature, pressure, and relative humidity in a network of sensor nodes. For each identified precipitation period, we decompose spectra into static- and radon-associated components as defined by a non-negative matrix factorization (NMF) algorithm. We find a consistent power-law relationship between a precipitation-dependent peak of the radon progeny proxy (RPP) and the peak strength of the radon-associated NMF component for most precipitation events. We conduct a case study of a rainfall period with abnormally high levels of implied radon progeny concentration and describe its temporal and spatial evolution. We hypothesize that this phenomenon is due to the air mass path that intersects a uranium-rich region of Wyoming. Finally, we cluster precipitation events into three distinct categories. One category roughly corresponds to events with deep low-pressure systems and high relative radon concentration, while another is characteristic of light stratiform rain with slightly higher temperatures and intermediate relative radon concentration. The third category appears to contain weak-gradient or lake breeze convection showers with intermittent precipitation and low relative radon concentration. These findings suggest that radiological anomaly detection could be improved by training unique background models corresponding to each category of meteorological event.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Understanding Anomalous Gas-Phase Peak Shifts in Dip-and-Pull Ambient Pressure XPS Experiments

Dip-and-pull ambient pressure X-ray photoelectron spectroscopy (AP-XPS) holds promise to uncover elementary processes of (photo)electrochemistry. We show, however, that the sample for such experiments should preferably be nonporous and the potential on the surface homogeneous. We carried out dip-and-pull AP-XPS experiments on a hematite thin film sample under the photoelectrochemical oxygen evolution reaction (OER) and find unexpected O 1s core level shifts. Upon electrochemical biasing under simulated solar light illumination, the gas-phase water peak shifted more than the electrolyte peak. To uncover the origin of the unexpected larger shift of the gas-phase peak, we performed electrostatic simulations using COMSOL, to map the potential field in the relevant volume between the sample and the first aperture of the XPS spectrometer. A number of geometric models were considered. We find that when the potential on the sample surface is inhomogeneous, e.g., with ionically isolated electrolyte patches, the gas-phase peak of the spectrum can shift more than the peak due to the electrolyte film. This suggests that at the measured sample position, the local potential was not as set by the potentiostat. Despite this, we find reversible consumption and replenishment of hydroxide in the spectra, which, due to OH– being the reactant of the OER in alkaline electrolyte, makes sense chemically. We propose that this is linked to OH– diffusion across the measured sample position, driven by the potential dependent consumption and replenishment of OH– at the nearby well-connected surface regions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Ultrasound Characterization of Plastic Bonded Explosives With Varied Binder Content and Density

Advanced explosive formulations and methods of production require rapid characterization for both processes and materials. Here, we present results for ultrasound characterization of a plastic bonded explosive (PBX) formulation as a function of binder weight percent and pressing density using through-transmission sound speed measurements. The explosive in this study was composed of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) combined with the fluoropolymer binder Kynar Ultraflex B. Binder content ranged from 2.6 wt.% to 11.5 wt.% with associated explosive formulation densities from 1.69 g/cc to 1.89 g/cc. We find linear and nonlinear changes in longitudinal and shear sound speeds, peak frequencies, and measured elastic properties for the different PBX formulations. These measured changes can differentiate the tested formulations, highlighting the potential for ultrasound as a nondestructive characterization method for newly-formulated PBXs. We also demonstrate through-transmission measurements using different coupling media and present results from an ultrasound technique called resonant ultrasound spectroscopy (RUS). Both transmission and RUS approaches may allow for rapid and even automated characterization of PBX components suitable for high-throughput manufacturing operations.

Engineering↗
DOE OSTI2026

Bipolar Membrane Electrodialyzers as Flexible Demand Response Resources: Co-Optimization of Cost Savings and Product Formation

Bipolar membrane electro dialyzers (BPMED) are widely used for chemical production and processing, including in the emerging ocean alkalinity enhancement (OAE) industry. In this paper, we explore the potential of BPMED devices as flexible electrochemical loads within power system operations. Using a multi-objective optimization framework, we evaluate BPMED operation across 24-hour and monthly horizons to examine how dispatch strategies respond to electricity price and grid conditions. Simulation results show that altering the relative weights of the choices in the objective function strongly shape the operating patterns, with cost-focused strategies that suppress the operation during peak prices. Furthermore, we propose alternative formulations that optimize operations to achieve both cost savings and alignment with periods of lower grid-side carbon intensity (CI), as low grid-side CI is key to maximize OAE efficiency. Additionally, a detailed sensitivity analysis highlights the importance of device properties, where low area-specific resistance (ASR) of membrane and high current efficiency (CE) are observed to jointly unlock cost-effective operation. However, even modest shunt efficiency losses are observed to erode performance and decrease system value. Importantly, the analysis demonstrates that BPMED can serve as a controllable and flexible demand response resource, shifting load to support multiple grid-side objectives, including (but not limited to) renewable integration, alleviate peak demand, and provide co-benefits for system reliability. These findings underscore BPMED’s dual role as a process technology and a grid-supporting asset, pointing to promising pathways for operational optimization of multiple objectives.

Bhattacharya, Saptarshi (ORCID:0000000308902060)↗
DOE OSTI2024

Jacobian-scaled K-means clustering for physics-informed segmentation of reacting flows

This work introduces Jacobian-scaled K-means (JSK-means) clustering, which is a physicsinformed clustering strategy centered on the K-means framework. The method allows for the injection of underlying physical knowledge into the clustering procedure through a distance function modification: instead of leveraging conventional Euclidean distance vectors, the JSKmeans procedure operates on distance vectors scaled by matrices obtained from dynamical system Jacobians evaluated at the cluster centroids. The goal of this work is to show how the JSKmeans algorithm - without modifying the input dataset - produces clusters that capture regions of dynamical similarity, in that the clusters are redistributed towards high-sensitivity regions in phase space and are described by similarity in the source terms of samples instead of the samples themselves. The algorithm is demonstrated on a complex reacting flow simulation dataset (a channel detonation configuration), where the dynamics in the thermochemical composition space are known through the highly nonlinear and stiff Arrhenius-based chemical source terms. Interpretations of cluster partitions in both physical space and composition space reveal how JSK-means shifts clusters produced by standard K-means towards regions of high chemical sensitivity (e.g., towards regions of peak heat release rate near the detonation reaction zone). Furthermore, the findings presented here illustrate the benefits of utilizing Jacobian-scaled distances in clustering techniques, and the JSK-means method in particular displays promising potential for improving former partition-based modeling strategies in reacting flow (and other multi-physics) applications.

Clustering↗
DOE OSTI2025

(Un)Explained EMIC Waves: Understanding Quiet Time EMIC Wave Drivers

Electromagnetic ion cyclotron (EMIC) waves are known to be generated through cyclotron resonance with the local ion particle population and grow when there is a large enough temperature anisotropy. In general, these temperature anisotropies necessary for wave growth are found to be associated with either solar wind pressure pulses or particle injections during geomagnetic storms or substorms. However, some EMIC events do not show any clear association with these known drivers and appear unexplained. In our analysis of high-amplitude (>1 nT) non-storm time EMIC waves, we find that 24 out of the 223 (~11%) EMIC events with peak amplitude greater than 1 nT were not found to be associated with any clear EMIC wave driver. This raises a compelling question: What magnetospheric or solar wind driver provides the free energy to grow these quiet-time EMIC waves? Here, we examine two EMIC events on 13 April 2017 and 13 February 2014, which were excited during extremely quiet solar wind and geomagnetic conditions. Furthermore, an in-depth analysis of field and particle measurements from multiple data sets, including ground and in situ data for these two events, indicates that extremely weak and otherwise insignificant pressure values and/or very weak substorm injections occurring multiple hours before the event play a significant role in quiet time wave generation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Sample glue layer investigation and mitigation for laser induced prompt impulse experiments

Understanding longer timescale material reactions under dynamic stress loading is critical for applications in materials engineering, shock physics, and planetary science. Prompt impulse experiments generate lower pressures since the ablator—the material first removed by the laser—is thicker and farther from the diagnostic plane, capturing aggregate material responses from the initial shock wave, rarefaction waves, and later time effects. This complexity demands thorough material characterization and simulation support. Since traditional sample construction is specific to supported shock experiments, designing prompt impulse experiments requires reconsideration around target design and sample engineering. Here, we present sample preparation techniques, experimental investigations, and theoretical simulations to investigate glue layer impacts, aiming to standardize samples for consistent data at lower laser fluences. We find that glue layers <30 μm have a minimal impact on peak velocity and pulse shape. The peak velocity scales linearly with glue layer thickness until a glue layer of 75 μm. For glue layers >75 μm, the peak velocity no longer scales with thickness; however, the pulse shape continues to degrade as described by simulations.

Lasers↗
DOE OSTI2020

Thickness dependence of electronic structure and optical properties of a correlated van der Waals antiferromagnetic NiPS 3 thin film

In this work, we study the thickness dependence of the electronic, magnetic, and optical properties of a NiPS 3 thin film, which is an antiferromagnetic charge-transfer insulator. Utilizing state-of-the-art advanced density functionals, we find the antiferromagnetic zigzag order, the band gap, and the main peaks in the dielectric tensor are all in good agreement with the corresponding experimental values. Upon thinning, the zigzag antiferromagnetic order becomes virtually degenerate with a competing Néel order, consistent with the suppression of long-range order observed by Raman spectroscopy due to strong magnetic fluctuations. Additionally, due to the robustness of the electronic band gap observed by spectroscopic ellipsometry measurements above T N , we suggest that the persistence of the band gap is driven by strong electronic correlations. Other systematic changes in electronic dispersion, effective mass, and Kerr angle with thickness are also discussed. Finally, an applied external electric field is found to suppress the band gap by up to 13%, until precipitating an insulator-metal transition at a critical field value of 0.7 eV/Å.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dirac spectral density in N f = 2 + 1 QCD at T = 230 MeV

We compute the renormalized Dirac spectral density in N f = 2 + 1 QCD at physical quark masses, temperature T = 230 MeV , and system size L s = 3.4 fm . To that end, we perform a pointwise continuum limit of the staggered density in lattice QCD with staggered quarks. We find, for the first time, that a clear infrared structure (IR peak) emerges in the density of the Dirac operator describing dynamical quarks. We also provide numerical evidence that a component of this peak, which becomes dominant in the thermodynamic limit, is due to a nontrivial accumulation of near-zero modes. Features of this structure are consistent with those previously attributed to the recently proposed IR phase of thermal QCD. Our results (i) provide the only complete first-principles evidence that these IR features exist and are physical, (ii) improve the upper bound for IR phase transition temperature T IR so that the new window is 200 < T IR < 230 MeV , and (iii) are consistent with the nonrestoration of anomalous U A ( 1 ) symmetry (chiral limit) below T = 230 MeV . Published by the American Physical Society 2024

Alexandru, Andrei (ORCID:0000000345471554)↗
DOE OSTI2021

Active matter commensuration and frustration effects on periodic substrates

In this report we show that self-driven particles coupled to a periodic obstacle array exhibit active matter commensuration effects that are absent in the Brownian limit. As the obstacle size is varied for sufficiently large activity, a series of commensuration effects appear in which the motility induced phase separation produces commensurate crystalline states, while for other obstacle sizes we find frustrated or amorphous states. The commensuration effects are associated with peaks in the amount of sixfold ordering and the maximum cluster size. When a drift force is added to the system, the mobility contains peaks and dips similar to those found in transport studies for commensuration effects in superconducting vortices and colloidal particles.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Optimal Demand Response Incorporating Distribution LMP with PV Generation Uncertainty

The utilization of aggregated demand-side flexibility via demand response (DR) has become a promising pathway for the integration of renewable energy resources in power systems. Nowadays, there are several management strategies for DR such as the price-based transactive control strategies. However, many of such existing price-based control strategies neglect the physics and operational constraints of the underlying distribution networks when computing the price, raising concerns regarding their theoretical and practical values. This paper studies this issue and investigates optimal DR (ODR) by incorporating the distribution locational marginal price (DLMP). In particular, we discuss DR in connection with DLMPs and propose a multiperiod bilevel optimization problem to find the ODR strategy. Here, the objective is to minimize the peak load, load fluctuation, and payments of load aggregators. In addition, a robust bilevel ODR model is formulated to provide a robust ODR strategy while minimizing operating costs under the worst-case realization of uncertainties; this mitigates the impact of forecasting errors on renewable energy resources. Then, we propose an efficient solution approach by employing the Karush-Kuhn-Tucker conditions and strong duality. Simulation results are presented to illustrate the mutual impacts of the interaction between DR and DLMP and the benefits of the robust ODR strategy.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Volcanic Imprints in Last-Millennium Land Summer Temperatures in the Circum-North Atlantic Area

Summer cooling is one of the most direct consequences of explosive volcanic eruptions that can affect ecosystems and human societies. Recent studies revealed a multiyear cooling impact on hemispheric and global summer temperatures after tropical eruptions, yet the volcanic responses appear to vary on regional scales. Here, in this work, we revisit volcano-induced summer cooling in eastern Canada and northern and central Europe by applying superposed epoch analysis on CMIP6-PMIP4 simulations and millennial temperature reconstructions based on tree-ring density. We then examine potential causes modulating region-specific volcanic impact. While confirming that, on average, tropical eruptions over the last millennium have induced a longer cooling (>4 yr) than eruptions from extratropical Northern Hemisphere in all three North Atlantic regions, we show that the peak magnitude of cooling is stronger in eastern Canada. We also find that the detected volcanic temperature anomalies can be strongly affected by the selection and number of volcanic events and nonvolcanic signals embedded in the climate time series. This study highlights the risks of using highly noisy proxy records to investigate volcanic impacts, especially in regions with strong unforced climate variability. The CMIP6-PMIP4 simulations generally agree with the three reconstructions on the average response to tropical eruptions, but their performance is poorer regarding the production of significant cooling after extratropical eruptions. Our results further suggest that the particular sensitivity to tropical eruptions in eastern Canada is likely related to increased sea ice surrounding Quebec–Labrador associated with the positive Arctic Oscillation and North Atlantic Oscillation formed during the first posteruption winter.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Structure Formation Limits on Axion-Like Dark Matter

We derive structure formation limits on dark matter (DM) composed of keV-scale axion-like particles (ALPs), produced via freeze-in through the interactions with photons and Standard Model (SM) fermions. We employ Lyman-alpha (Ly-α) forest data sets as well as the observed number of Milky Way (MW) subhalos. We compare results obtained using Maxwell-Boltzmann and quantum statistics for describing the SM bath. It should be emphasized that the presence of logarithmic divergences complicates the calculation of the production rate, which can not be parameterized with a simple power law behaviour. The obtained results, in combination with X-ray bounds, exclude the possibility for a photophilic “frozen-in” ALP DM with mass below ~ 19. For the photophobic ALP scenario, in which DM couples primarily to SM fermions, the ALP DM distribution function is peaked at somewhat lower momentum and hence for such realization we find weaker limits on DM mass. Future facilities, such as the upcoming Vera C. Rubin observatory, will provide measurements with which the current bounds can be significantly improved to ~ 80.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗