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 235 records · Page 13

Fast Simulations of Gas Sloshing and Cold Front Formation

We present a simplified and fast method for simulating minor mergers between galaxy clusters. Instead of following the evolution of the dark matter halos directly by the N-body method, we employ a rigid potential approximation for both clusters. The simulations are run in the rest frame of the more massive cluster and account for the resulting inertial accelerations in an optimised way. We test the reliability of this method for studies of minor merger induced gas sloshing by performing a one-to-one comparison between our simulations and hydro+N-body ones. We find that the rigid potential approximation reproduces the sloshing-related features well except for two artefacts: the temperature just outside the cold fronts is slightly over-predicted, and the outward motion of the cold fronts is delayed by typically 200 Myr. We discuss reasons for both artefacts.

Roediger, E.↗

Fast Simulations of Gas Sloshing and Cold Front Formation

We present a simplified and fast method for simulating minor mergers between galaxy clusters. Instead of following the evolution of the dark matter halos directly by the N-body method, we employ a rigid potential approximation for both clusters. The simulations are run in the rest frame of the more massive cluster and account for the resulting inertial accelerations in an optimised way. We test the reliability of this method for studies of minor merger induced gas sloshing by performing a one-to-one comparison between our simulations and hydro+N-body ones. We find that the rigid potential approximation reproduces the sloshing-related features well except for two artifacts: the temperature just outside the cold fronts is slightly over-predicted, and the outward motion of the cold fronts is delayed by typically 200 Myr. We discuss reasons for both artifacts.

Roediger, E.↗

Galaxy formation with wave/fuzzy dark matter: The core-halo structure and the solitonic imprint

Dark matter-dominated cores have long been claimed for the well-studied local group dwarf galaxies. More recently, extended stellar halos have been uncovered around several of these dwarfs through deeper imaging and spectroscopy. Such core-halo structures (inner flat core and a characteristic r −3 asymptotic outer halo profile) are not a feature of conventional cold dark matter (CDM). In contrast, smooth and prominent dark matter cores are predicted for wave/fuzzy dark matter (ψDM). The question arises as to what extent the visible stellar profiles should reflect this dark matter core structure. Here we compare cosmological hydrodynamical simulations of CDM, “WDM” (model used as a proxy for ψDM) & ψDM, aiming to predict the stellar profiles for these three DM scenarios. We show that cores surrounded by extended halos are distinguishable for ψDM, where the stellar density is enhanced in the core due to the presence of the relatively dense soliton. Our analysis demonstrates that, in our simulations, a distinctive core-halo structure does not appear in the case of CDM in the DM, gas, or stars. Whereas we do find a core-halo transition for DM, gas, and stars for ψDM, and the scale of this transition is in line with the predicted core radius set by the soliton scale anticipated for the adopted boson mass of 2.5×10 −22 eV. The presence of a core-halo structure in the stellar profile for Galaxy 1 for ψDM is visible for the most massive and the first galaxy to form in the simulation. Clearly, further simulations are needed to establish how strict this possible relationship is between the DM and stellar core-halo profile as a potential observational discriminator. Furthermore, we observe the anticipated asymmetry for ψDM due to the soliton's motion (jumping and random walk), a distinctive characteristic not found in the symmetric distributions of stars in the warm and CDM models.

dark matter↗

Constraining dark matter microphysics with the annihilation signal from subhalos

In the cold dark matter scenario, galactic dark matter halos are populated with a large number of smaller subhalos. Previous work has shown that dark matter annihilations in subhalos can generate a distinctive, non-Poisson signal in the gamma-ray photon counts probability distribution function (PDF). Here we show that the gamma-ray PDF also carries information about the velocity dependence of the dark matter annihilation cross section. After calculating the PDF assuming s-wave and Sommerfeld-enhanced annihilation, we perform a mock data analysis to illustrate how current and future observations can constrain the microphysics of the dark matter annihilation. We find that, with current Fermi data, and assuming a dark matter annihilation cross section roughly at the limit of current bounds from annihilation in dwarf spheroidal galaxies, one can potentially distinguish the non-Poissonian fluctuations expected from dark matter annihilation in subhalos from Poisson sources, as well as from dark matter models with an incorrect velocity-dependence. Here we explore how robust these results are to assumptions about the modeling of the galactic gamma-ray background, but further work is needed to determine the impact of realistic astrophysical source populations on our results. We also point out a four-parameter degeneracy between the velocity dependence of the dark matter annihilation, the minimum subhalo mass, the power law index of the subhalo mass function, and the normalization of the dark matter signal. This degeneracy can be broken with priors from N-body simulations or from observational constraints on the subhalo mass function.

79 ASTRONOMY AND ASTROPHYSICS↗

Illuminating dark matter halo density profiles without subhaloes

Cold dark matter haloes consist of a relatively smooth dark matter component as well as a system of bound subhaloes. It is the prevailing practice to include all mass, including mass in subhaloes, in studies of halo density profiles in simulations. However, often in observational studies satellites are treated as having their own distinct dark matter density profiles in addition to the profile of the host. This difference can make comparisons between theoretical and observed results difficult. In this work, we investigate density profiles of the smooth components of host haloes by excluding mass contained within subhaloes. We find that the density profiles of the smooth halo component (without subhaloes) differ substantially from the conventional halo density profile, declining more rapidly at large radii. We also find that concentrations derived from smooth density profiles exhibit less scatter at fixed mass and a weaker mass dependence than standard concentrations. Both smooth and standard halo profiles can be described by a generalized Einasto profile, an Einasto profile with a modified central slope, with smaller residuals than either a Navarro–Frenk–White or Einasto profile. Furthermore, these results hold for both Milky Way-mass and cluster-mass haloes. This new characterization of smooth halo profiles can be useful for many analyses, such as lensing and dark matter annihilation, in which the smooth and clumpy components of a halo should be accounted for separately.

79 ASTRONOMY AND ASTROPHYSICS↗

From Images to Dark Matter: End-to-end Inference of Substructure from Hundreds of Strong Gravitational Lenses

Abstract Constraining the distribution of small-scale structure in our universe allows us to probe alternatives to the cold dark matter paradigm. Strong gravitational lensing offers a unique window into small dark matter halos (<10 10 M ⊙ ) because these halos impart a gravitational lensing signal even if they do not host luminous galaxies. We create large data sets of strong lensing images with realistic low-mass halos, Hubble Space Telescope (HST) observational effects, and galaxy light from HST’s COSMOS field. Using a simulation-based inference pipeline, we train a neural posterior estimator of the subhalo mass function (SHMF) and place constraints on populations of lenses generated using a separate set of galaxy sources. We find that by combining our network with a hierarchical inference framework, we can both reliably infer the SHMF across a variety of configurations and scale efficiently to populations with hundreds of lenses. By conducting precise inference on large and complex simulated data sets, our method lays a foundation for extracting dark matter constraints from the next generation of wide-field optical imaging surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Maximum entropy distributions of dark matter in ΛCDM cosmology

Context. Small-scale challenges to ΛCDM cosmology require a deeper understanding of dark matter physics. Aims. This paper aims to develop the maximum entropy distributions for dark matter particle velocity (denoted by X ), speed (denoted by Z ), and energy (denoted by E ) that are especially relevant on small scales where system approaches full virialization. Methods. For systems involving long-range interactions, a spectrum of halos of different sizes is required to form to maximize system entropy. While the velocity in halos can be Gaussian, the velocity distribution throughout the entire system, involving all halos of different sizes, is non-Gaussian. With the virial theorem for mechanical equilibrium, we applied the maximum entropy principle to the statistical equilibrium of entire system, such that the maximum entropy distribution of velocity (the X distribution) could be analytically derived. The halo mass function was not required in this formulation, but it did indeed result from the maximum entropy. Results. The predicted X distribution involves a shape parameter α and a velocity scale, v 0 . The shape parameter α reflects the nature of force ( α → 0 for long-range force or α → ∞ for short-range force). Therefore, the distribution approaches Laplacian with α → 0 and Gaussian with α → ∞. For an intermediate value of α , the distribution naturally exhibits a Gaussian core for v ≪ v 0 and exponential wings for v ≫ v 0 , as confirmed by N -body simulations. From this distribution, the mean particle energy of all dark matter particles with a given speed, v , follows a parabolic scaling for low speeds (∝ v 2 for v ≪ v 0 in halo core region, i.e., “Newtonian”) and a linear scaling for high speeds (∝ v for v ≫ v 0 in halo outskirt, i.e., exhibiting “non-Newtonian” behavior due to long-range gravity). We compared our results against N -body simulations and found a good agreement.

79 ASTRONOMY AND ASTROPHYSICS↗

Insights into cosmological structure formation with machine learning.

Our modern understanding of cosmological structure formation posits that small matter density fluctuations present in the early Universe, as traced by the cosmic microwave background, grow via gravitational instability to form extended haloes of dark matter. A theoretical understanding of the structure, evolution and formation of dark matter haloes is an essential step towards unravelling the intricate connection between halo and galaxy formation, needed to test our cosmological model against data from upcoming galaxy surveys. Physical understanding of the process of dark matter halo formation is made difficult by the highly non-linear nature of the haloes' evolution. I describe a new approach to gain physical insight into cosmological structure formation based on machine learning. This approach combines the ability of machine learning algorithms to learn non-linear relationships, with techniques that enable us to physically interpret the learnt mapping. I describe applicatio ns of th e method, with the aim of investigating which aspects of the early universe density field impact the later formation of dark matter haloes. First I present a case where the process of halo formation is turned into a binary classification problem; the algorithm predicts whether or not dark matter `particles' in the initial conditions of a simulation will collapse into haloes of a given mass range. Second, I present its generalization to regression, where the algorithm infers the final mass of the halo to which each particle will later belong. I show that the initial tidal shear does not play a significant role compared to the initial density field in establishing final halo masses. Finally, I demonstrate that extending the framework to deep learning algorithms such as convolutional neural networks allows us to explore connections between the early universe and late time haloes beyond those studied by existing analytic approximations of halo collapse.

79 ASTRONOMY AND ASTROPHYSICS↗

The cosmic Mach number - Direct comparisons of observations and models

A robust method is used to derive the cosmic Mach number, M, from existing peculiar velocity data sets. For the IR Tully-Fisher distances of spirals in the Local Supercluster of Aaronson et al. (1979, 1982), M of 1.03 and a characteristic distance from the Local Group of 1464 km/s are found. For the elliptical galaxy sample of Faber et al. (1987, 1988), M of 0.57, with a characteristic distance of 2572 km/s is found. These results are compared with Monte Carlo simulations of the observational realizations drawn from numerical simulations of the universe based on various scenarios, including the standard cold dark matter (CDM) scenario. The effect of velocity bias on the derived Mach number is found to be small. Only 5 percent of the CDM simulations have M as large as or larger than that derived for the Aaronson et al. sample.

Strauss, Michael A.↗

Hot gas in superclusters and microwave background distortions

Variations in the electron pressure can produce significant distortions of the cosmic microwave background radiation (CMBR) spectrum through the Sunyaev-Zel'dovich effect. In this paper, we compute the amplitude of these distortions in various cosmological scenarios using the results of hydrodynamical simulations combined with analytic theory. For reionized cold dark matter (CDM) with h = 0.5 probed with the OVRO experiment (Theta approximately = 0.1 deg), we find that the amplitude of these fluctuations is (Delta T/T)(sub rms) = 1 x 10(exp -5) (Omega(sub i)/0.06) if the experiment looks at a random part of the sky and 6 x 10(exp -6)(Omega(sub i)/(0.06) if bright areas are avoided. Here Omega is the density of ionized baryons in units of the critical density, rho(sub c) = 2.78 x 10(exp 11)h(exp 2)/cu Mpc. Primary fluctuations would generate (Delta T/T)(sub rms) approximately = 2 x 10(exp -6) for this model. These distortions could be a dominant source of CMBR anisotropics on small and intermediate angular scales. It is important to note that these signals are generated relatively recently, when nonlinear effects are improtant, and do not require early reionization. We also examine the contribution to anisotropies from the Doppler effect, in both the linear and nonlinear regimes. In this case, the bulk of the signal is generated before a redshift of 10, provided that the intercluster medium is ionized back to the surface of last scattering. On the smallest scales thermal effects are bigger, while, on larger scales the Doppler-induced fluctuations are greater with the exact boundary determined by the scenario. However, over a wide range of scales, the two signals are equal to within a factor of a few.

Persi, Fred M.↗

Precision Kinematic Sunyaev--Zel'dovich Measurements Across Halo Mass and Redshift with DESI DR2 and ACT DR6: Part I. Luminous Red Galaxies

We present the most precise measurements of the kinetic Sunyaev-Zel'dovich (kSZ) effect around luminous red galaxies to date, detecting the signal at $18σ$ significance in both harmonic and configuration space. Our analysis cross-correlates 2.4 million spectroscopic LRGs from the Dark Energy Spectroscopic Instrument (DESI) DR2 sample with Data Release 6 (DR6) of the Atacama Cosmology Telescope (ACT). We develop a novel harmonic-space cross-correlation approach using momentum-weighted kSZ templates, yielding nearly uncorrelated bandpowers within a framework consistent with other large-scale structure analyses. By incorporating the LRG halo occupation distribution (HOD) and its uncertainty, we convert measured galaxy gas profiles into halo gas profiles and provide generalized Navarro-Frenk-White (GNFW) fitting profiles, providing empirical targets for tuning feedback efficiency in hydrodynamical simulations and for baryonic modeling in large-scale structure analyses. We find strong evidence that gas profiles do not trace dark matter, providing direct evidence for gas redistribution beyond gravitational collapse. Comparing to hydrodynamical simulations, our measurements favor feedback efficiencies exceeding those in the Battaglia profile, suggesting more efficient gas ejection in group-scale halos than previously predicted. Splitting by redshift, we detect the kSZ signal at SNR $\approx 5$--$10$ in each of four bins and find amplitude evolution consistent with the expected decline in mean halo mass at fixed comoving number density. Splitting by stellar mass, we study the scaling of kSZ amplitude with galaxy properties. Together with BGS and ELG measurements in Paper II, these results span $0.1 \lesssim z \lesssim 1.6$ across three galaxy populations, demonstrating the potential of spectroscopic kSZ to map circumgalactic gas and constrain baryonic feedback.

Qu, F. J. [KIPAC, Menlo Park; Stanford U., HEPL; S↗

Strong gravitational lensing statistics as a test of cosmogonic scenarios

Gravitational lensing statistics can provide a direct and powerful test of cosmic structure formation theories. Since lensing tests, directly, the magnitude of the nonlinear mass density fluctuations on lines of sight to distant objects, no issues of 'bias' (of mass fluctuations with respect to galaxy density fluctuations) exist here, although lensing observations provide their own ambiguities of interpretation. We develop numerical techniques for generating model density distributions with the very large spatial dynamic range required by lensing considerations and for identifying regions of the simulations capable of multiple image lensing in a conservative and computationally efficient way that should be accurate for splittings significantly larger than 3 seconds. Applying these techniques to existing standard Cold dark matter (CDM) (Omega = 1) and Primeval Baryon Isocurvature (PBI) (Omega = 0.2) simulations (normalized to the Cosmic Background Explorer Satellite (COBE) amplitude), we find that the CDM model predicts large splitting (greater than 8 seconds) lensing events roughly an order-of-magnitude more frequently than the PBI model. Under the reasonable but idealized assumption that lensing structrues can be modeled as singular isothermal spheres (SIS), the predictions can be directly compared to observations of lensing events in quasar samples. Several large splitting (Delta Theta is greater than 8 seconds) cases are predicted in the standard CDM model (the exact number being dependent on the treatment of amplification bias), whereas none is observed. In a formal sense, the comparison excludes the CDM model at high confidence (essentially for the same reason that CDM predicts excessive small-scale cosmic velocity dispersions.) A very rough assessment of low-density but flat CDM model (Omega = 0.3, Lambda/3H(sup 2 sub 0) = 0.7) indicates a far lower and probably acceptable level of lensing. The PBI model is consistent with, but not strongly tested by, the available lensing data, and other open models would presumably do as well as PBI. These preliminary conclusions and the assumptions on which they are based can be tested and the analysis can be applied to other cosmogonic models by straightforward extension of the work presented here.

Cen, Renyue↗

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↗

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↗

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↗

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↗

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↗