Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dark matter simulations”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Core formation in high-z massive haloes: heating by post-compaction satellites and response to AGN outflows

ABSTRACT Observed rotation curves in star-forming galaxies indicate a puzzling dearth of dark matter in extended flat cores within haloes of mass ≥1012M⊙ at z ∼ 2. This is not reproduced by current cosmological simulations, and supernova-driven outflows are not effective in such massive haloes. We address a hybrid scenario where post-compaction merging satellites heat up the dark-matter cusps by dynamical friction, allowing active galactic nucleus (AGN)-driven outflows to generate cores. Using analytic and semi-analytic models (SatGen), we estimate the dynamical friction heating as a function of satellite compactness for a cosmological sequence of mergers. Cosmological simulations (VELA) demonstrate that satellites of initial virial masses >1011.3M⊙, which undergo wet compactions, become sufficiently compact for significant heating. Constituting a major fraction of the accretion on to haloes ≥1012M⊙, these satellites heat up the cusps in half a virial time at z ∼ 2. Using a model for outflow-driven core formation (CuspCore), we demonstrate that the heated dark-matter cusps develop extended cores in response to removal of half the gas mass, while the more compact stellar systems remain intact. The mergers keep the dark matter hot, while the gas supply, fresh and recycled, is sufficient for the AGN outflows. AGNs indeed become effective in haloes ≥1012M⊙, where the black hole growth is no longer suppressed by supernovae and its compaction-driven rapid growth is maintained by a hot circumgalactic medium. For simulations to reproduce the dynamical friction effects, they should resolve the compaction of the massive satellites and avoid artificial tidal disruption. AGN feedback could be boosted by clumpy black hole accretion and clumpy response to AGN.

Dekel, Avishai↗
DOE OSTI2021

Wave Dark Matter

Here we review the physics and phenomenology of wave dark matter: a bosonic dark matter candidate lighter than about 30 eV. Such particles have a de Broglie wavelength exceeding the average interparticle separation in a galaxy like the Milky Way and are, thus, well described as a set of classical waves. We outline the particle physics motivations for such particles, including the quantum chromodynamics axion as well as ultralight axion-like particles such as fuzzy dark matter. The wave nature of the dark matter implies a rich phenomenology: Wave interference gives rise to order unity density fluctuations on de Broglie scale in halos. One manifestation is vortices where the density vanishes and around which the velocity circulates. There is one vortex ring per de Broglie volume on average. For sufficiently low masses, soliton condensation occurs at centers of halos. The soliton oscillates and undergoes random walks, which is another manifestation of wave interference. The halo and subhalo abundance is expected to be suppressed at small masses, but the precise prediction from numerical wave simulations remains to be determined. For ultralight ∼10 −22 eV dark matter, the wave interference substructures can be probed by tidal streams or gravitational lensing. The signal can be distinguished from that due to subhalos by the dependence on stream orbital radius or image separation. Axion detection experiments are sensitive to interference substructures for wave dark matter that is moderately light. The stochastic nature of the waves affects the interpretation of experimental constraints and motivates the measurement of correlation functions. Current constraints and open questions, covering detection experiments and cosmological, galactic, and black hole observations, are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Mixed warm dark matter constraints using Milky Way satellite galaxy counts

Warm dark matter has been strongly constrained in recent years as the sole component of dark matter. However, a less-explored alternative is that dark matter consists of a mixture of warm and cold dark matter (MWDM). In this work, we use observations of Milky Way satellite galaxies to constrain MWDM scenarios where the formation of small-scale structure is suppressed either by generic thermal relic warm dark matter or a sterile neutrino produced through the Shi-Fuller mechanism. To achieve this, we model satellite galaxies by combining numerical simulations with semianalytical models for the subhalo population, and use a galaxy-halo connection model to match galaxies onto dark matter subhalos. By comparing the number of satellites predicted by MWDM models to the observed satellite population from the Dark Energy Survey and Pan-STARRS1, we constrain the fraction of warm dark matter, f WDM , as a function of its mass, m WDM . We exclude dark matter being composed entirely of thermal relic warm dark matter with m WDM ≤ 6.6 keV at a posterior ratio of 10 ∶ 1 , consistent with previous works. However, we find that warm dark matter with smaller mass is allowed when mixed with cold dark matter, and that the f WDM constraints strengthen with decreasing m WDM until they plateau at f WDM ≲ 0.45 for m WDM ≲ 1.5 keV . Likewise, in the case of a sterile neutrino with mass of 7 keV produced through the Shi-Fuller mechanism, we exclude a fraction of f ν s ≲ 0.45 , independent of mixing angle. Our results extend constraints on MWDM to a region of parameter space that has been relatively unconstrained with previous analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Earth-catalyzed detection of magnetic inelastic dark matter with photons in large underground detectors

Inelastic dark matter with moderate splittings, $\mathcal{O}$ (few to 150) keV, can upscatter to an excited state in the Earth, with the excited state subsequently decaying, leaving a distinctive monoenergetic photon signal in large underground detectors. The photon signal can exhibit sidereal-daily modulation, providing excellent separation from backgrounds. Using a detailed numerical simulation, we examine this process as a search strategy for magnetic inelastic dark matter with the dark matter mass near the weak scale, where the upscatter to the excited state and decay proceed through the same magnetic dipole transition operator. At lower inelastic splittings, the scattering is dominated by moderate mass elements in the Earth with high spin, especially 27 Al, while at larger splittings, 56 Fe becomes the dominant target. We show that the proposed large volume gaseous detector CYGNUS will have excellent sensitivity to this signal. Xenon detectors also provide excellent sensitivity through the inelastic nuclear recoil signal, and if a future signal is seen, we show that the synergy among both types of detection can provide strong evidence for magnetic inelastic dark matter. In the course we have calculated nuclear response functions for elements relevant for scattering in the Earth, which are publicly available on GitHub.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Comparative Analysis of Cosmological Simulations via Angular Power Spectra

Cosmology is the study of the universe's origin, large-scale structure, and dynamics. One such area of study is understanding dark energy and dark matter. Dark matter has not been observed to interact with luminous matter other than through gravitational attractions, forming large scale structures observed throughout the universe. In this area, the primary model is the Lambda Cold Dark Matter (ΛCDM) model. This standard model consists of cosmological dark energy (Λ), cold dark matter (CDM), and normal luminous matter. We use angular power spectra (APS) to compare simulated ΛCDM models with observed large-scale galaxy distributions from the DECam Local Volume Exploration Survey (DELVE) [1]. After ensuring that the simulations can produce APS that match observational results, we vary cosmological parameters to constrain the allowed ranges of parameter values that are consistent with the observed data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Comparative analysis of cosmological simulations via angular power spectra

Cosmology is the study of the universe's origin, large-scale structure, and dynamics. One such area of study is understanding dark energy and dark matter. Dark matter has not been observed to interact with luminous matter other than through gravitational attractions, forming large scale structures observed throughout the universe. In this area, the primary model is the Lambda Cold Dark Matter (ΛCDM) model. This standard model consists of cosmological dark energy (Λ), cold dark matter (CDM), and normal luminous matter. We use Angular Power Spectra (APS) to compare simulated ΛCDM models with observed large-scale galaxy distributions from the DECam Local Volume Exploration Survey (DELVE). After ensuring that the simulations can produce APS which match observational results, we vary cosmological parameters to constrain the allowed ranges of parameter values that are consistent with the observed data.

Feithen, William↗
DOE OSTI2023

Convergence of small scale Ly α structure at high- z under different reionization scenarios

ABSTRACT The Ly α forest (LAF) at z > 5 probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-z, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting for a late end to reionization, small structures down to the sub-kpc level may survive to later times than conventionally thought due to the reduced time for pressure smoothing to impact the gas, further tightening simulation resolution requirements. We perform a suite of simulations using the Eulerian cosmological hydrodynamics code Nyx, spanning domain sizes of 1.25 − 10 h−1 Mpc and 5 − 80 h−1 kpc cells, and explore the interaction of these variables with the timing of reionization on the properties of the matter distribution and the simulated LAF at z = 5.5. In observable Ly α power, convergence within 10 per cent is achieved for k < 0.1 s km–1, but larger k shows deviation of up to 20 per cent. While a later reionization retains more small structure in the density field, because of the greater thermal broadening there is little difference in the convergence of LAF power between early (z = 9) and later (z = 6) reionizations. We conclude that at z ∼ 5.5, resolutions of 10 kpc are necessary for convergence of LAF power at k < 0.1 s km–1, while higher-k modes require higher resolution, and that the timing of reionization does not significantly impact convergence given realistic photoheating.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A pair of early- and late-forming galaxy cluster samples: A novel way of studying halo assembly bias assisted by a constrained simulation

The halo assembly bias, a phenomenon referring to dependencies of the large-scale bias of a dark matter halo other than its mass, is a fundamental property of the standard cosmological model. First discovered in 2005 from the Millennium Run simulation, it has been proven very difficult to be detected observationally, with only a few convincing claims of detection so far. The main obstacle lies in finding an accurate proxy of the halo formation time. In this study, by utilizing a constrained simulation that can faithfully reproduce the observed structures larger than 2 Mpc in the local universe, for a sample of 634 massive clusters at z ≤ 0.12, we found their counterpart halos in the simulation and used the mass growth history of the matched halos to estimate the formation time of the observed clusters. This allowed us to construct a pair of early- and late-forming clusters, with a similar mass as measured via weak gravitational lensing, and large-scale biases differing at the ≈3σ level, suggestive of the signature of assembly bias, which is further corroborated by the properties of cluster galaxies, including the brightest cluster galaxy and the spatial distribution and number of member galaxies. Our study paves a way to further detect assembly bias based on cluster samples constructed purely on observed quantities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Analytical estimation of the signal to noise ratio efficiency in axion dark matter searches using a Savitzky-Golay filter

The signal to noise ratio efficiency ϵ SNR in axion dark matter searches has been estimated using large-statistic simulation data reflecting the background information and the expected axion signal power obtained from a real experiment. This usually requires a lot of computing time even with the assistance of powerful computing resources. Employing a Savitzky-Golay filter for background subtraction, in this work, we estimated a fully analytical ϵ SNR without relying on large-statistic simulation data, but only with an arbitrary axion mass and the relevant signal shape information. Hence, our work can provide ϵ SNR using minimal computing time and resources prior to the acquisition of experimental data, without the detailed information that has to be obtained from real experiments. Axion haloscope searches have been observing the coincidence that the frequency independent scale factor ξ is approximately consistent with the ϵ SNR . This was confirmed analytically in this work, when the window length of the Savitzky-Golay filter is reasonably wide enough, i.e., at least 5 times the signal window.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Stellar splashback: the edge of the intracluster light

We examine the outskirts of galaxy clusters in the C-EAGLE simulations to quantify the ‘edges’ of the stellar and dark matter distribution. The radius of the steepest slope in the dark matter, commonly used as a proxy for the splashback radius, is located at $\sim \, r_{200 \rm m}$; the strength and location of this feature depends on the recent mass accretion rate, in good agreement with previous work. Interestingly, the stellar distribution (or intracluster light, ICL) also has a well-defined edge, which is directly related to the splashback radius of the halo. Thus, detecting the edge of the ICL can provide an independent measure of the physical boundary of the halo, and the recent mass accretion rate. We show that these caustics can also be seen in the projected density profiles, but care must be taken to account for the influence of substructures and other non-diffuse material, which can bias and/or weaken the signal of the steepest slope. This is particularly important for the stellar material, which has a higher fraction bound in subhaloes than the dark matter. Finally, we show that the ‘stellar splashback’ feature is located beyond current observational constraints on the ICL, but these large projected distances (>> 1 Mpc) and low surface brightnesses (μ>>32 mag arcsec -2 ) can be reached with upcoming observational facilities such as the Vera C. Rubin Observatory, the Nancy Grace Roman Space Telescope, and Euclid.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Axions: From magnetars and neutron star mergers to beam dumps and BECs

We review topics in searches for axion-like-particles (ALPs), covering material that is complementary to other recent reviews. The first half of our review covers ALPs in the extreme environments of neutron star cores, the magnetospheres of highly magnetized neutron stars (magnetars), and in neutron star mergers. The focus is on possible signals of ALPs in the photon spectrum of neutron stars and gravitational wave/electromagnetic signals from neutron star mergers. We then review recent developments in laboratory-produced ALP searches, focusing mainly on accelerator-based facilities including beam-dump type experiments and collider experiments. We provide a general-purpose discussion of the ALP search pipeline from production to detection, in steps, and our discussion is straightforwardly applicable to most beam-dump type and reactor experiments. We end with a selective look at the rapidly developing field of ultralight dark matter, specifically the formation of Bose–Einstein Condensates (BECs). We review the properties of BECs of ultralight dark matter and bridge these properties with developments in numerical simulations, and ultimately with their impact on experimental searches.

Astronomy & Astrophysics↗
DOE OSTI2020

New Measurements on the Neutrino Mass: KATRIN and Project 8 (Final Technical Report for DE-SC0014130)

The overarching goal of the funded research is to build experiments and develop techniques towards the major goal of measuring the mass of the neutrino, particularly by the method of high-resolution spectroscopy of the endpoint of tritium. The goals of this project were to make hardware and intellectual contributions towards the KATRIN experiment, an international project under construction in Germany, and towards the "Project 8" R&D program to develop new electron spectroscopy methods and apply them to tritium. A minor goal of the project was to support simulation work and R&D on future neutrino and dark matter detectors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Measuring the Density Fields around Bright Quasars at z ~ 6 with XQR-30 Spectra

Measuring the density of the intergalactic medium using quasar sight lines in the epoch of reionization is challenging due to the saturation of Ly α absorption. Near a luminous quasar, however, the enhanced radiation creates a proximity zone observable in the quasar spectra where the Ly α absorption is not saturated. In this study, we use 10 high-resolution ( R ≳ 10,000) z ~ 6 quasar spectra from the extended XQR-30 sample to measure the density field in the quasar proximity zones. We find a variety of environments within 3 pMpc distance from the quasars. We compare the observed density cumulative distribution function (CDF) with models from the Cosmic Reionization on Computers simulation and find a good agreement between 1.5 and 3 pMpc from the quasar. This region is far away from the quasar hosts and hence approaching the mean density of the universe, which allows us to use the CDF to set constraints on the cosmological parameter σ 8 = 0.6 ± 0.3. The uncertainty is mainly due to the limited number of high-quality quasar sight lines currently available. Utilizing the more than 200 known quasars at z ≳ 6, this method will allow us to tighten the constraint on σ 8 to the percent level in the future. In the region closer to the quasar within 1.5 pMpc, we find that the density is higher than predicted in the simulation by 1.23 ± 0.17, suggesting that the typical host dark matter halo mass of a bright quasar ( M 1450 < -26.5) at z ~ 6 is .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Measurement of charge and light yields for 127 Xe L -shell electron captures in liquid xenon

Dark matter searches using dual-phase xenon time-projection chambers (LXe-TPCs) rely on their ability to reject background electron recoils (ERs) while searching for signal-like nuclear recoils (NRs). ER response is typically calibrated using β -decay sources, such as tritium, but these calibrations do not characterize events accompanied by an atomic vacancy, as in solar neutrino scatters off inner-shell electrons. Such events lead to emission of x rays and Auger electrons, resulting in higher electron-ion recombination and thus a more NR-like response than inferred from β -decay calibration. We present a cross-calibration of tritium β -decays and Xe 127 electron-capture decays (which produce inner-shell vacancies) in a small-scale LXe-TPC and give the most precise measurements to date of light and charge yields for the Xe 127 L -shell electron-capture in liquid xenon. We observe a 6.9 σ ( 9.2 σ ) discrepancy in the L -shell capture response relative to tritium β decays, measured at a drift field of 363 ± 14 V / cm ( 258 ± 13 V / cm ), when compared to simulations tuned to reproduce the correct β -decay response. In dark matter searches, use of a background model that neglects this effect leads to overcoverage (higher limits) for background-only multi-kiloton-year exposures, but at a level much less than the 1 - σ experiment-to-experiment variation of the 90% C.L. upper limit on the interaction rate of a 50 GeV / c 2 dark matter particle.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Enhancing the sensitivity of the LUX-ZEPLIN (LZ) dark matter experiment to low energy signals

Two-phase xenon detectors, such as that at the core of the forthcoming LZ dark matter experiment, use photomultiplier tubes to sense the primary (S1) and secondary (S2) scintillation signals resulting from particle interactions in their liquid xenon target. This paper describes a simulation study exploring two techniques to lower the energy threshold of LZ to gain sensitivity to low-mass dark matter and astrophysical neutrinos, which will be applicable to other liquid xenon detectors. The energy threshold is determined by the number of detected S1 photons; typically, these must be recorded in three or more photomultiplier channels to avoid dark count coincidences that mimic real signals. To lower this threshold: a) we take advantage of the double photoelectron emission effect, whereby a single vacuum ultraviolet photon has a $\sim20\%$ probability of ejecting two photoelectrons from a photomultiplier tube photocathode; and b) we drop the requirement of an S1 signal altogether, and use only the ionization signal, which can be detected more efficiently. For both techniques we develop signal and background models for the nominal exposure, and explore accompanying systematic effects, including the dependence on the free electron lifetime in the liquid xenon. When incorporating double photoelectron signals, we predict a factor of $\sim 4$ sensitivity improvement to the dark matter-nucleon scattering cross-section at $2.5$ GeV/c$^2$, and a factor of $\sim1.6$ increase in the solar $^8$B neutrino detection rate. Dropping the S1 requirement may allow sensitivity gains of two orders of magnitude in both cases. Finally, we apply these techniques to even lower masses by taking into account the atomic Migdal effect; this could lower the dark matter particle mass threshold to $80$ MeV/c$^2$.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Enhancing the sensitivity of the LUX-ZEPLIN (LZ) dark matter experiment to low energy signals

Two-phase xenon detectors, such as that at the core of the forthcoming LZ dark matter experiment, use photomultiplier tubes to sense the primary (S1) and secondary (S2) scintillation signals resulting from particle interactions in their liquid xenon target. This paper describes a simulation study exploring two techniques to lower the energy threshold of LZ to gain sensitivity to low-mass dark matter and astrophysical neutrinos, which will be applicable to other liquid xenon detectors. The energy threshold is determined by the number of detected S1 photons; typically, these must be recorded in three or more photomultiplier channels to avoid dark count coincidences that mimic real signals. To lower this threshold: a) we take advantage of the double photoelectron emission effect, whereby a single vacuum ultraviolet photon has a $\sim20\%$ probability of ejecting two photoelectrons from a photomultiplier tube photocathode; and b) we drop the requirement of an S1 signal altogether, and use only the ionization signal, which can be detected more efficiently. For both techniques we develop signal and background models for the nominal exposure, and explore accompanying systematic effects, including the dependence on the free electron lifetime in the liquid xenon. When incorporating double photoelectron signals, we predict a factor of $\sim 4$ sensitivity improvement to the dark matter-nucleon scattering cross-section at $2.5$ GeV/c$^2$, and a factor of $\sim1.6$ increase in the solar $^8$B neutrino detection rate. Dropping the S1 requirement may allow sensitivity gains of two orders of magnitude in both cases. Finally, we apply these techniques to even lower masses by taking into account the atomic Migdal effect; this could lower the dark matter particle mass threshold to $80$ MeV/c$^2$.

Akerib, DS↗
DOE OSTI2026

Probabilistic Inference of Low-Surface-Brightness Galaxy Morphological Parameters Using Simulation-Based Inference

Low-surface-brightness galaxies (LSBGs) are diffuse, often dark-matter-dominated systems whose faintness makes their structural parameters difficult to measure reliably in wide-field imaging surveys. Robust parameter inference, including uncertainty quantification, is important for population studies and for comparisons with models of galaxy formation, as future surveys are expected to produce increasingly large samples of diffuse galaxies. In practice, LSBG profile modeling is sensitive to sky- background errors, masking choices, contaminating background sources, and the computational cost of obtaining posterior-level uncertainties for large samples. Motivated by these questions, we develop a simulation-based inference (SBI) framework for estimating posterior distributions of LSBG morphological parameters from simulated galaxy images. Using PyImfit, we generate DES-like single-Sersic profile LSBG images with known position angle, ellipticity, Sersic index, effective surface brightness, and effective radius. We then train a normalizing-flow-based neural posterior estimator using the sbi package to infer these parameters from the simulated images. For isolated simulated galaxies, the SBI posterior recovers the true input parameters, produces posterior predictive residuals consistent with the assumed noise model, and shows good empirical calibration in a DES-motivated test regime. We also compare SBI with PyImfit-based MCMC inference and find broadly comparable posterior constraints, while SBI enables substantially faster posterior sampling after training. Finally, we test robustness to compact background contaminants. A model trained only on isolated galaxies produces undercovered posteriors on contaminated images, whereas training on simulations with variable contaminant positions and fluxes improves calibration across contaminated test sets. These results demonstrate the promise of SBI for scalable, uncertainty-aware LSBG morphology inference, while emphasizing that posterior reliability strongly depends on whether training simulations include relevant observational complications.

Batbayar, Bilguun [U. Chicago (main)]↗
DOE OSTI2023

Cosmological structure formation and soliton phase transition in fuzzy dark matter with axion self-interactions

Herein we investigate cosmological structure formation in fuzzy dark matter (FDM) with the attractive self-interaction (SI) with numerical simulations. Such a SI would arise if the FDM boson were an ultra-light axion, which has a strong CP symmetry-breaking scale (decay constant). Although weak, the attractive SI may be strong enough to counteract the quantum ‘pressure’ and alter structure formation. We find in our simulations that the SI can enhance small-scale structure formation, and soliton cores above a critical mass undergo a phase transition, transforming from dilute to dense solitons.

79 ASTRONOMY AND ASTROPHYSICS↗