Cold dark matter heats up

Explaining Dark Matter Halo Density Profiles with Neural Networks

We use explainable neural networks to connect the evolutionary history of dark matter halos with their density profiles. The network captures independent factors of variation in the density profiles within a low-dimensional representation, which we physically interpret using mutual information. Without any prior knowledge of the halos' evolution, the network recovers the known relation between the early time assembly and the inner profile and discovers that the profile beyond the virial radius is described by a single parameter capturing the most recent mass accretion rate. The results illustrate the potential for machine-assisted scientific discovery in complicated astrophysical datasets.

Search for Gamma-Ray Spectral Lines from Dark Matter Annihilation up to 100 TeV toward the Galactic Center with MAGIC

Linelike features in TeV γ rays constitute a "smoking gun" for TeV-scale particle dark matter and new physics. Probing the Galactic Center region with ground-based Cherenkov telescopes enables the search for TeV spectral features in immediate association with a dense dark matter reservoir at a sensitivity out of reach for satellite γ-ray detectors, and direct detection and collider experiments. We report on 223 hours of observations of the Galactic Center region with the MAGIC stereoscopic telescope system reaching γ-ray energies up to 100 TeV. We improved the sensitivity to spectral lines at high energies using large-zenith-angle observations and a novel background modeling method within a maximum-likelihood analysis in the energy domain. No linelike spectral feature is found in our analysis. Therefore, we constrain the cross section for dark matter annihilation into two photons to ⟨σv⟩≲5×10^{-28}  cm^{3} s^{-1} at 1 TeV and ⟨σv⟩≲1×10^{-25}  cm^{3} s^{-1} at 100 TeV, achieving the best limits to date for a dark matter mass above 20 TeV and a cuspy dark matter profile at the Galactic Center. Finally, we use the derived limits for both cuspy and cored dark matter profiles to constrain supersymmetric wino models.

The Cusp–Core Problem in Gas-Poor Dwarf Spheroidal Galaxies

This review deals with the inconsistency of inner dark matter density profiles in dwarf galaxies, known as the cusp–core problem. In particular, we aim to focus on gas-poor dwarf galaxies. One of the most promising solutions to this cold dark matter small-scale issue is the stellar feedback, but it seems to be only designed for gas-rich dwarfs. However, in the regime of classical dwarfs, this core mechanism becomes negligible. Therefore, it is required to find solutions without invoking these baryonic processes as dark matter cores tend to persist even for these dwarfs, which are rather dark-matter-dominated. Here, we have presented two categories of solutions. One consists of creating dark matter cores from cusps within cold dark matter by altering the dark matter potential via perturbers. The second category gathers solutions that depict the natural emergence of dark matter cores in alternative theories. Given the wide variety of solutions, it becomes necessary to identify which mechanism dominates in the central region of galaxies by finding observational signatures left by them in order to highlight the true nature of dark matter.

Review of Solutions to the Cusp-Core Problem of the ΛCDM Model

This review aims at proposing to the field an overview of the Cusp-core problem, including a discussion of its advocated solutions, assessing how each can satisfactorily provide a description of central densities. Whether the Cusp-core problem reflects our insufficient grasp on the nature of dark matter, of gravity, on the impact of baryonic interactions with dark matter at those scales, as included in semi-analytical models or fully numerical codes, the solutions to it can point either to the need for a paradigm change in cosmology, or to to our lack of success in ironing out the finer details of the ΛCDM paradigm.

Fundamental Properties of the Dark and the Luminous Matter from the Low Surface Brightness Discs

Dark matter (DM) is one of the biggest mystery in the Universe. In this review, we start reporting the evidences for this elusive component and discussing about the proposed particle candidates and scenarios for such phenomenon. Then, we focus on recent results obtained for rotating disc galaxies, in particular for low surface brightness (LSB) galaxies. The main observational properties related to the baryonic matter in LSBs, investigated over the last decades, are briefly recalled. Next, these galaxies are analyzed by means of the mass modelling of their rotation curves both individual and stacked. The latter analysis, via the universal rotation curve (URC) method, results really powerful in giving a global or universal description of the properties of these objects. We report the presence in LSBs of scaling relations among their structural properties that result comparable with those found in galaxies of different morphologies. All this confirms, in disc systems, the existence of a strong entanglement between the luminous matter (LM) and the dark matter (DM). Moreover, we report how in LSBs the tight relationship between their radial gravitational accelerations g and their baryonic components gb results to depend also on the stellar disk length scale and the radius at which the two accelerations have been measured. LSB galaxies strongly challenge the ΛCDM scenario with the relative collisionless dark particle and, alongside with the non-detection of the latter, contribute to guide us towards a new scenario for the DM phenomenon.

Strong Bound on Canonical Ultralight Axion Dark Matter from the Lyman-Alpha Forest

We present a new bound on the ultralight axion (ULA) dark matter mass m_{a}, using the Lyman-alpha forest to look for suppressed cosmic structure growth: a 95% lower limit m_{a}>2×10^{-20}  eV. This strongly disfavors (>99.7% credibility) the canonical ULA with 10^{-22}  eV<m_{a}<10^{-21}  eV, motivated by the string axiverse and solutions to possible tensions in the cold dark matter model. We strengthen previous equivalent bounds by about an order of magnitude. We demonstrate the robustness of our results using an optimized emulator of improved hydrodynamical simulations.

First Star-Forming Structures in Fuzzy Cosmic Filaments

In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter, we investigate how the first galaxies are assembled in a "fuzzy" dark matter (FDM) cosmology where dark matter is an ultralight ∼10^{-22}  eV boson and the primordial stars are expected to form along dense dark matter filaments. Using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical solitonlike cores, which are unstable under gravity and collapse into kiloparsec-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter-a smoking gun signature of FDM.

A universal constant for dark matter-baryon interplay

Recent studies point out that there exists some rough scaling relations for dark matter and some tight connections between dark matter and baryons. However, most of the relations and tight connections can only be found in galaxies, but not in galaxy clusters. In this article, we consider a new expression that can characterize the properties of dark matter-baryon interplay for both galactic and galaxy cluster scales. By using the archival observational data of galaxies and galaxy clusters, we show that the value [Formula: see text] is almost a constant and scale independent within the optical radius ro, where [Formula: see text] is the average dark matter number density, [Formula: see text] is the average baryon number density, v is the characteristic velocity and V is the interacting volume. This would be the first universal relation between dark matter and baryons on both galactic and galaxy cluster scales. We anticipate this result to be a starting point to explain the small-scale problem and the scaling relations for dark matter in galaxies. The constant K discovered may reveal some underlying global interaction between dark matter and baryons.

Where are the most ancient stars in the Milky Way?

The oldest stars in the Milky Way (MW) bear imprints of the Galaxy's early assembly history. We use FIRE cosmological zoom-in simulations of three MW-mass disc galaxies to study the spatial distribution, chemistry, and kinematics of the oldest surviving stars (z form ≳ 5) in MW-like galaxies. We predict the oldest stars to be less centrally concentrated at z = 0 than stars formed at later times as a result of two processes. First, the majority of the oldest stars are not formed in situ but are accreted during hierarchical assembly. These ex situ stars are deposited on dispersion-supported, halo-like orbits but dominate over old stars formed in situ in the solar neighbourhood, and in some simulations, even in the galactic centre. Secondly, old stars formed in situ are driven outwards by bursty star formation and energetic feedback processes that create a time-varying gravitational potential at z ≳ 2, similar to the process that creates dark matter cores and expands stellar orbits in bursty dwarf galaxies. The total fraction of stars that are ancient is more than an order of magnitude higher for sight lines away from the bulge and inner halo than for inward-looking sight lines. Although the task of identifying specific stars as ancient remains challenging, we anticipate that million-star spectral surveys and photometric surveys targeting metal-poor stars already include hundreds of stars formed before z = 5. We predict most of these targets to have higher metallicity (-3 < [Fe/H] < -2) than the most extreme metal-poor stars.

The origin of ultra diffuse galaxies: stellar feedback and quenching

We test if the cosmological zoom-in simulations of isolated galaxies from the FIRE project reproduce the properties of ultra diffuse galaxies (UDGs). We show that outflows that dynamically heat galactic stars, together with a passively aging stellar population after imposed quenching, naturally reproduce the observed population of red UDGs, without the need for high spin haloes, or dynamical influence from their host cluster. We reproduce the range of surface brightness, radius, and absolute magnitude of the observed red UDGs by quenching simulated galaxies at a range of different times. They represent a mostly uniform population of dark matter-dominated dwarf galaxies with M * ~ 108 M☉, low metallicity, and a broad range of ages; the more massive the UDGs, the older they are. The most massive red UDG in our sample(M * ~ 3 × 108 M☉) requires quenching at z ~ 3 when its halo reached M h ~ 1011M☉. Our simulated UDGs form with normal stellar-to-halo ratios and match the central enclosed masses and the velocity dispersions of the observed UDGs. Enclosed masses remain largely fixed across a broad range of quenching times because the central regions of their dark matter haloes complete their growth early. If our simulated dwarfs are not quenched, they evolve into bluer low surface brightness galaxies with M/L similar to observed field dwarfs. While our simulation sample covers a limited range of formation histories and halo masses, we predict that UDG is a common, and perhaps even dominant, galaxy type around M * ~ 108 M☉, both in the field and in clusters.

Gas kinematics in FIRE simulated galaxies compared to spatially unresolved HI observations

The shape of a galaxy's spatially unresolved, globally integrated 21-cm emission line depends on its internal gas kinematics: galaxies with rotationally supported gas discs produce double-horned profiles with steep wings, while galaxies with dispersion-supported gas produce Gaussian-like profiles with sloped wings. Using mock observations of simulated galaxies from the FIRE project, we show that one can therefore constrain a galaxy's gas kinematics from its unresolved 21-cm line profile. In particular, we find that the kurtosis of the 21-cm line increases with decreasing V/σ and that this trend is robust across a wide range of masses, signal-to-noise ratios, and inclinations. We then quantify the shapes of 21-cm line profiles from a morphologically unbiased sample of ~2000 low-redshift, HI-detected galaxies with Mstar = 107-11 M☉ and compare to the simulated galaxies. At Mstar ≳ 1010 M☉, both the observed and simulated galaxies produce double-horned profiles with low kurtosis and steep wings, consistent with rotationally supported discs. Both the observed and simulated line profiles become more Gaussian like (higher kurtosis and less-steep wings) at lower masses, indicating increased dispersion support. However, the simulated galaxies transition from rotational to dispersion support more strongly: at Mstar 108-10 M, most of the simulations produce more Gaussian-like profiles than typical observed galaxies with similar mass, indicating that gas in the low-mass simulated galaxies is, on average, overly dispersion supported. Most of the lower-mass-simulated galaxies also have somewhat lower gas fractions than the median of the observed population. The simulations nevertheless reproduce the observed line-width baryonic Tully-Fisher relation, which is insensitive to rotational versus dispersion support.

The different baryonic Tully-Fisher relations at low masses

We compare the Baryonic Tully-Fisher relation (BTFR) of simulations and observations of galaxies ranging from dwarfs to spirals, using various measures of rotational velocity V_rot. We explore the BTFR when measuring V_rot at the flat part of the rotation curve, V_flat, at the extent of H i gas, V_last, and using 20 per cent (W₂₀) and 50 per cent (W₅₀) of the width of H i line profiles. We also compare with the maximum circular velocity of the parent halo, [Formula: see text], within dark matter only simulations. The different BTFRs increasingly diverge as galaxy mass decreases. Using V_last one obtains a power law over four orders of magnitude in baryonic mass, with slope similar to the observed BTFR. Measuring V_flat gives similar results as V_last when galaxies with rising rotation curves are excluded. However, higher rotation velocities would be found for low-mass galaxies if the cold gas extended far enough for V_rot to reach a maximum. W₂₀ gives a similar slope as V_last but with slightly lower values of V_rot for low-mass galaxies, although this may depend on the extent of the gas in your galaxy sample. W₅₀ bends away from these other relations towards low velocities at low masses. By contrast, [Formula: see text] bends towards high velocities for low-mass galaxies, as cold gas does not extend out to the radius at which haloes reach [Formula: see text]. Our study highlights the need for careful comparisons between observations and models: one needs to be consistent about the particular method of measuring V_rot, and precise about the radius at which velocities are measured.

Adaptive techniques for clustered N-body cosmological simulations

ChaNGa is an N-body cosmology simulation application implemented using Charm++. In this paper, we present the parallel design of ChaNGa and address many challenges arising due to the high dynamic ranges of clustered datasets. We propose optimizations based on adaptive techniques. We evaluate the performance of ChaNGa on highly clustered datasets: a $z \sim0$ z ∼ 0 snapshot of a 2 billion particle realization of a 25 Mpc volume, and a 52 million particle multi-resolution realization of a dwarf galaxy. For the 25 Mpc volume, we show strong scaling on up to 128K cores of Blue Waters. We also demonstrate scaling up to 128K cores of a multi-stepping run of the 2 billion particle simulation. While the scaling of the multi-stepping run is not as good as single stepping, the throughput at 128K cores is greater by a factor of 2. We also demonstrate strong scaling on up to 512K cores of Blue Waters for two large, uniform datasets with 12 and 24 billion particles.

The nongravitational interactions of dark matter in colliding galaxy clusters

Collisions between galaxy clusters provide a test of the nongravitational forces acting on dark matter. Dark matter's lack of deceleration in the "bullet cluster" collision constrained its self-interaction cross section σ(DM)/m < 1.25 square centimeters per gram (cm(2)/g) [68% confidence limit (CL)] (σ(DM), self-interaction cross section; m, unit mass of dark matter) for long-ranged forces. Using the Chandra and Hubble Space Telescopes, we have now observed 72 collisions, including both major and minor mergers. Combining these measurements statistically, we detect the existence of dark mass at 7.6σ significance. The position of the dark mass has remained closely aligned within 5.8 ± 8.2 kiloparsecs of associated stars, implying a self-interaction cross section σ(DM)/m < 0.47 cm(2)/g (95% CL) and disfavoring some proposed extensions to the standard model.

Constraints on an annihilation signal from a core of constant dark matter density around the milky way center with H.E.S.S

An annihilation signal of dark matter is searched for from the central region of the Milky Way. Data acquired in dedicated on-off observations of the Galactic center region with H.E.S.S. are analyzed for this purpose. No significant signal is found in a total of ∼9  h of on-off observations. Upper limits on the velocity averaged cross section, ⟨σv⟩, for the annihilation of dark matter particles with masses in the range of ∼300  GeV to ∼10  TeV are derived. In contrast to previous constraints derived from observations of the Galactic center region, the constraints that are derived here apply also under the assumption of a central core of constant dark matter density around the center of the Galaxy. Values of ⟨σv⟩ that are larger than 3×10^{-24}  cm^{3}/s are excluded for dark matter particles with masses between ∼1 and ∼4  TeV at 95% C.L. if the radius of the central dark matter density core does not exceed 500 pc. This is the strongest constraint that is derived on ⟨σv⟩ for annihilating TeV mass dark matter without the assumption of a centrally cusped dark matter density distribution in the search region.

Understanding the core-halo relation of quantum wave dark matter from 3D simulations

We examine the nonlinear structure of gravitationally collapsed objects that form in our simulations of wavelike cold dark matter, described by the Schrödinger-Poisson (SP) equation with a particle mass ∼10(-22)  eV. A distinct gravitationally self-bound solitonic core is found at the center of every halo, with a profile quite different from cores modeled in the warm or self-interacting dark matter scenarios. Furthermore, we show that each solitonic core is surrounded by an extended halo composed of large fluctuating dark matter granules which modulate the halo density on a scale comparable to the diameter of the solitonic core. The scaling symmetry of the SP equation and the uncertainty principle tightly relate the core mass to the halo specific energy, which, in the context of cosmological structure formation, leads to a simple scaling between core mass (Mc) and halo mass (Mh), Mc∝a(-1/2)Mh(1/3), where a is the cosmic scale factor. We verify this scaling relation by (i) examining the internal structure of a statistical sample of virialized halos that form in our 3D cosmological simulations and by (ii) merging multiple solitons to create individual virialized objects. Sufficient simulation resolution is achieved by adaptive mesh refinement and graphic processing units acceleration. From this scaling relation, present dwarf satellite galaxies are predicted to have kiloparsec-sized cores and a minimum mass of ∼10(8)M⊙, capable of solving the small-scale controversies in the cold dark matter model. Moreover, galaxies of 2×10(12)M⊙ at z=8 should have massive solitonic cores of ∼2×10(9)M⊙ within ∼60  pc. Such cores can provide a favorable local environment for funneling the gas that leads to the prompt formation of early stellar spheroids and quasars.

Galactic center gamma-ray excess from dark matter annihilation: is there a black hole spike?

If the supermassive black hole Sgr A* at the center of the Milky Way grew adiabatically from an initial seed embedded in a Navarro-Frenk-White dark matter (DM) halo, then the DM profile near the hole has steepened into a spike. We calculate the dramatic enhancement to the gamma-ray flux from the Galactic center (GC) from such a spike if the 1-3 GeV excess observed in Fermi data is due to DM annihilations. We find that for the parameter values favored in recent fits, the point-source-like flux from the spike is 35 times greater than the flux from the inner 1° of the halo, far exceeding all Fermi point source detections near the GC. We consider the dependence of the spike signal on astrophysical and particle parameters and conclude that if the GC excess is due to DM, then a canonical adiabatic spike is disfavored by the data. We discuss alternative Galactic histories that predict different spike signals, including (i) the nonadiabatic growth of the black hole, possibly associated with halo and/or black hole mergers, (ii) gravitational interaction of DM with baryons in the dense core, such as heating by stars, or (iii) DM self-interactions. We emphasize that the spike signal is sensitive to a different combination of particle parameters than the halo signal and that the inclusion of a spike component to any DM signal in future analyses would provide novel information about both the history of the GC and the particle physics of DM annihilations.