-
Constraints on the properties of $ν$MSM dark matter using the satellite galaxies of the Milky Way
Authors:
Oliver Newton,
Mark R. Lovell,
Carlos S. Frenk,
Adrian Jenkins,
John C. Helly,
Shaun Cole,
Andrew J. Benson
Abstract:
Low-mass galaxies provide a powerful tool with which to investigate departures from the standard cosmological paradigm in models that suppress the abundance of small dark matter structures. One of the simplest metrics that can be used to compare different models is the abundance of satellite galaxies in the Milky Way. Viable dark matter models must produce enough substructure to host the observed…
▽ More
Low-mass galaxies provide a powerful tool with which to investigate departures from the standard cosmological paradigm in models that suppress the abundance of small dark matter structures. One of the simplest metrics that can be used to compare different models is the abundance of satellite galaxies in the Milky Way. Viable dark matter models must produce enough substructure to host the observed number of Galactic satellites. Here, we scrutinize the predictions of the neutrino Minimal Standard Model ($ν{\rm MSM}$), a well-motivated extension of the Standard Model of particle physics in which the production of sterile neutrino dark matter is resonantly enhanced by a lepton asymmetry in the primordial plasma. This process enables the model to evade current constraints associated with non-resonantly produced dark matter. Independently of assumptions about galaxy formation physics we rule out, with at least 95 per cent confidence, all parameterizations of the $ν{\rm MSM}$ with sterile neutrino rest mass, $M_{\rm s} \leq 1.4\, {\rm keV}$. Incorporating physically motivated prescriptions of baryonic processes and modelling the effects of reionization strengthen our constraints, and we exclude all $ν{\rm MSM}$ parameterizations with $M_{\rm s} \leq 4\, {\rm keV}$. Unlike other literature, our fiducial constraints do not rule out the putative 3.55 keV X-ray line, if it is indeed produced by the decay of a sterile neutrino; however, some of the most favoured parameter space is excluded. If the Milky Way satellite count is higher than we assume, or if the Milky Way halo is less massive than $M^{\rm MW}_{200} = 8 \times 10^{11}\, {\rm M_\odot}$, we rule out the $ν{\rm MSM}$ as the origin of the 3.55 keV excess. In contrast with other work, we find that the constraints from satellite counts are substantially weaker than those reported from X-ray non-detections.
△ Less
Submitted 28 August, 2024;
originally announced August 2024.
-
A warm dark matter cosmogony may yield more low-mass galaxy detections in 21-cm surveys than a cold dark matter one
Authors:
Kyle A. Oman,
Carlos S. Frenk,
Robert A. Crain,
Mark R. Lovell,
Joel Pfeffer
Abstract:
The 21-cm spectral line widths, $w_{50}$, of galaxies are an approximate tracer of their dynamical masses, such that the dark matter halo mass function is imprinted in the number density of galaxies as a function of $w_{50}$. Correcting observed number counts for survey incompleteness at the level of accuracy needed to place competitive constraints on warm dark matter (WDM) cosmological models is…
▽ More
The 21-cm spectral line widths, $w_{50}$, of galaxies are an approximate tracer of their dynamical masses, such that the dark matter halo mass function is imprinted in the number density of galaxies as a function of $w_{50}$. Correcting observed number counts for survey incompleteness at the level of accuracy needed to place competitive constraints on warm dark matter (WDM) cosmological models is very challenging, but forward-modelling the results of cosmological hydrodynamical galaxy formation simulations into observational data space is more straightforward. We take this approach to make predictions for an ALFALFA-like survey from simulations using the EAGLE galaxy formation model in both cold (CDM) and WDM cosmogonies. We find that for WDM cosmogonies more galaxies are detected at the low-$w_{50}$ end of the 21-cm velocity width function than in the CDM cosmogony, contrary to what might naïvely be expected from the suppression of power on small scales in such models. This is because low-mass galaxies form later and retain more gas in WDM cosmogonies (with EAGLE). While some shortcomings in the treatment of cold gas in the EAGLE model preclude placing definitive constraints on WDM scenarios, our analysis illustrates that near-future simulations with more accurate modelling of cold gas will likely make strong constraints possible, especially in conjunction with new 21-cm surveys such as WALLABY.
△ Less
Submitted 19 July, 2024; v1 submitted 22 January, 2024;
originally announced January 2024.
-
Dark matter decay in the Milky Way halo
Authors:
Mark R. Lovell
Abstract:
Dark matter may be detected in X-ray decay, including from the decay of the dark matter particles that make up the Milky Way (MW) halo. We use a range of density profiles to compute X-ray line intensity profiles, with a focus on the resonantly produced sterile neutrino dark matter candidate. Compared to the Navarro--Frenk--White density profile, we show that using an adiabatically contracted halo…
▽ More
Dark matter may be detected in X-ray decay, including from the decay of the dark matter particles that make up the Milky Way (MW) halo. We use a range of density profiles to compute X-ray line intensity profiles, with a focus on the resonantly produced sterile neutrino dark matter candidate. Compared to the Navarro--Frenk--White density profile, we show that using an adiabatically contracted halo profile suppresses the line intensity in the halo outskirts and enhances it in the Galactic Centre (GC), although this enhancement is eliminated by the likely presence of a core within 3~kpc. Comparing our results to MW halo observations, other X-ray observations, and structure formation constraints implies a sterile neutrino mixing angle parameter $s_{11}\equiv\sin^{2}(2θ)\times10^{11}\sim[3,4]$ (particle lifetime $τ_{28}\equivτ/(10^{28}\mathrm{sec})\sim[1.0,1.3]$), which is nevertheless is strong tension with some reported non-detections. We make predictions for the likely decay flux that the XRISM satellite would measure in the GC, plus the Virgo and Perseus clusters, and outline further steps to determine whether the dark matter is indeed resonantly produced sterile neutrinos as detected in X-ray decay.
△ Less
Submitted 10 January, 2024;
originally announced January 2024.
-
The impact of free-streaming on dwarf galaxy counts in low-density regions
Authors:
Tamar Meshveliani,
Mark R. Lovell,
Robert A. Crain,
Joel Pfeffer
Abstract:
We study the statistics of dwarf galaxy populations as a function of environment in the cold dark matter (CDM) and warm dark matter (WDM) cosmogonies, using hydrodynamical simulations starting from initial conditions with matched phases but differing power spectra, and evolved with the EAGLE galaxy formation model. We measure the abundance of dwarf galaxies within 3~Mpc of DM haloes with a present…
▽ More
We study the statistics of dwarf galaxy populations as a function of environment in the cold dark matter (CDM) and warm dark matter (WDM) cosmogonies, using hydrodynamical simulations starting from initial conditions with matched phases but differing power spectra, and evolved with the EAGLE galaxy formation model. We measure the abundance of dwarf galaxies within 3~Mpc of DM haloes with a present-day halo mass similar to that of the Milky Way (MW), and find that the radial distribution of galaxies $M_{*}>10^7$\Msun is nearly identical for WDM and CDM. However, the cumulative mass function becomes shallower for WDM at lower masses, yielding 50~per~cent fewer dwarf galaxies of $M_{*}\gtrsim10^{5}$~\Msun than CDM. The suppression of low-mass halo counts in WDM relative to CDM increases significantly from high-density regions to low-density regions for haloes in the region of the half-mode mass, $M_\rm{hm}$. The luminous fraction in the two models also diverges from the overdense to the underdense regions for $M>2M_\rm{hm}$, as the increased collapse delay at small densities pushes the collapse to after the reionization threshold. However, the stellar mass--halo mass relation of WDM haloes relative to CDM increases towards lower-density regions. Finally, we conclude that the suppression of galaxies with $M_{*}\gtrsim10^5$\Msun between WDM and CDM is independent of density: the suppression of halo counts and the luminous fraction is balanced by an enhancement in stellar mass--halo mass relation.
△ Less
Submitted 30 November, 2023;
originally announced December 2023.
-
Halo assembly in cold and warm dark matter during the JWST frontier epoch
Authors:
Mark R. Lovell
Abstract:
The JWST mission is in the process of probing the galaxy mass function at $z>10$, when conceivably any delay in halo assembly due to the presence of a dwarf galaxy-scale power spectrum cutoff may drastically suppress the number of galaxies relative to the cold dark matter (CDM) expectation. We employ N-body simulations of CDM and warm dark matter (WDM) to explore how the difference in halo collaps…
▽ More
The JWST mission is in the process of probing the galaxy mass function at $z>10$, when conceivably any delay in halo assembly due to the presence of a dwarf galaxy-scale power spectrum cutoff may drastically suppress the number of galaxies relative to the cold dark matter (CDM) expectation. We employ N-body simulations of CDM and warm dark matter (WDM) to explore how the difference in halo collapse time between these models scales with $z=0$ descendant halo mass. We demonstrate that collapse begins first for the most massive haloes, and the delay in collapse time between CDM and WDM haloes correlates inversely with descendant mass. We thus infer that only present-day dwarf galaxies exhibit any difference in their assembly history between CDM and WDM at $z=10$, and therefore support previous studies that have found JWST is unlikely to determine whether our Universe is better described by the CDM cosmology or the WDM cosmology without favourable lensing studies.
△ Less
Submitted 1 June, 2023;
originally announced June 2023.
-
Anticipating the XRISM search for the decay of resonantly produced sterile neutrino dark matter
Authors:
Mark R. Lovell
Abstract:
The sterile neutrino ($N_1$) features in multiple extensions of the standard model and is a compelling dark matter candidate, especially as the decay of $N_1$ with mass $m_\mathrm{s}=7.1\mathrm{keV}$ is a possible source for the unexplained $3.55\mathrm{keV}$ X-ray line reported in galaxy clusters. This particle will be accessible to the XRISM X-ray mission over the next 12 months. We revisit the…
▽ More
The sterile neutrino ($N_1$) features in multiple extensions of the standard model and is a compelling dark matter candidate, especially as the decay of $N_1$ with mass $m_\mathrm{s}=7.1\mathrm{keV}$ is a possible source for the unexplained $3.55\mathrm{keV}$ X-ray line reported in galaxy clusters. This particle will be accessible to the XRISM X-ray mission over the next 12 months. We revisit the physics behind $N_1$ and the uncertainty in its parameters. We compare predictions for the $m_{s}=7.1\mathrm{keV}$ $N_1$ mixing angle, $\sin^2(2θ)$, and half-mode mass, $M_\mathrm{hm}$, as described in the $ν$MSM standard model extension to existing X-ray observations and structure formation constraints. The strongest available constraints rule out $N_1$ as a dark matter candidate, and a more optimistic reading of the data prefers $\sin^2(2θ)=5\times10^{-11}$ and $M_\mathrm{hm}=3.5\times10^{8}M_{\odot}$. We highlight that the most promising upcoming opportunity for a detection is to find a line of velocity dispersion $\sim500\mathrm{kms}^{-1}$ in the Virgo cluster with XRISM, and then draw up a list of future objects of study to determine: (i) whether the line is from dark matter generally, and (ii) if from dark matter, whether that candidate is indeed $N_1$.
△ Less
Submitted 3 August, 2023; v1 submitted 27 March, 2023;
originally announced March 2023.
-
Sensitivity of strong lensing observations to dark matter substructure: a case study with Euclid
Authors:
Conor M. O'Riordan,
Giulia Despali,
Simona Vegetti,
Mark R. Lovell,
Ángeles Moliné
Abstract:
We introduce a machine learning method for estimating the sensitivity of strong lens observations to dark matter subhaloes in the lens. Our training data include elliptical power-law lenses, Hubble Deep Field sources, external shear, and noise and PSF for the Euclid VIS instrument. We set the concentration of the subhaloes using a $v_\mathrm{max}$-$r_\mathrm{max}$ relation. We then estimate the da…
▽ More
We introduce a machine learning method for estimating the sensitivity of strong lens observations to dark matter subhaloes in the lens. Our training data include elliptical power-law lenses, Hubble Deep Field sources, external shear, and noise and PSF for the Euclid VIS instrument. We set the concentration of the subhaloes using a $v_\mathrm{max}$-$r_\mathrm{max}$ relation. We then estimate the dark matter subhalo sensitivity in $16{,}000$ simulated strong lens observations with depth and resolution resembling Euclid VIS images. We find that, with a $3σ$ detection threshold, $2.35$ per cent of pixels inside twice the Einstein radius are sensitive to subhaloes with a mass $M_\mathrm{max}\leq 10^{10}M_\odot$, $0.03$ per cent are sensitive to $M_\mathrm{max}\leq 10^{9}M_\odot$, and, the limit of sensitivity is found to be $M_\mathrm{max}=10^{8.8\pm0.2}M_\odot$. Using our sensitivity maps and assuming CDM, we estimate that Euclid-like lenses will yield $1.43^{+0.14}_{-0.11}[f_\mathrm{sub}^{-1}]$ detectable subhaloes per lens in the entire sample, but this increases to $35.6^{+0.9}_{-0.9}[f_\mathrm{sub}^{-1}]$ per lens in the most sensitive lenses. Estimates are given in units of the inverse of the substructure mass fraction $f_\mathrm{sub}^{-1}$. Assuming $f_\mathrm{sub}=0.01$, one in every $70$ lenses in general should yield a detection, or one in every $\sim$ three lenses in the most sensitive sample. From $170,000$ new strong lenses detected by Euclid, we expect $\sim 2500$ new subhalo detections. We find that the expected number of detectable subhaloes in warm dark matter models only changes relative to cold dark matter for models which have already been ruled out, i.e., those with half-mode masses $M_\mathrm{hm}>10^8M_\odot$.
△ Less
Submitted 27 February, 2023; v1 submitted 28 November, 2022;
originally announced November 2022.
-
Gravothermal collapse of Self-Interacting Dark Matter halos as the Origin of Intermediate Mass Black Holes in Milky Way satellites
Authors:
Tamar Meshveliani,
Jesús Zavala,
Mark R. Lovell
Abstract:
Milky Way (MW) satellites exhibit a diverse range of internal kinematics, reflecting in turn a diverse set of subhalo density profiles. These profiles include large cores and dense cusps, which any successful dark matter model must explain simultaneously. A plausible driver of such diversity is self-interactions between dark matter particles (SIDM) if the cross section passes the threshold for the…
▽ More
Milky Way (MW) satellites exhibit a diverse range of internal kinematics, reflecting in turn a diverse set of subhalo density profiles. These profiles include large cores and dense cusps, which any successful dark matter model must explain simultaneously. A plausible driver of such diversity is self-interactions between dark matter particles (SIDM) if the cross section passes the threshold for the gravothermal collapse phase at the characteristic velocities of the MW satellites. In this case, some of the satellites are expected to be hosted by subhalos that are still in the classical SIDM core phase, while those in the collapse phase would have cuspy inner profiles, with a SIDM-driven intermediate mass black hole (IMBH) in the centre as a consequence of the runaway collapse. We develop an analytical framework that takes into account the cosmological assembly of halos and is calibrated to previous simulations; we then predict the timescales and mass scales ($M_{\rm BH}$) for the formation of IMBHs in velocity-dependent SIDM (vdSIDM) models as a function of the present-day halo mass, $M_0$. Finally, we estimate the region in the parameter space of the effective cross section and $M_0$ for a subclass of vdSIDM models that result in a diverse MW satellite population, as well as their corresponding fraction of SIDM-collapsed halos and those halos' inferred IMBH masses. We predict the latter to be in the range $0.1-1000~ {\rm M_\odot}$ with a $M_{\rm BH}-M_0$ relation that has a similar slope, but lower normalization, than the extrapolated empirical relation of super-massive black holes found in massive galaxies.
△ Less
Submitted 24 May, 2023; v1 submitted 4 October, 2022;
originally announced October 2022.
-
Matching the mass function of Milky Way satellites in competing dark matter models
Authors:
Mark R. Lovell,
Jesús Zavala
Abstract:
Any successful model of dark matter must explain the diversity of observed Milky Way (MW) satellite density profiles, from very dense ultrafaints to large, low density satellites such as Crater~II that appear to be larger their anticipated host dark matter haloes. We consider cold dark matter (CDM), warm dark matter (WDM, 3.3keV thermal relic power spectrum), and a self-interacting dark matter mod…
▽ More
Any successful model of dark matter must explain the diversity of observed Milky Way (MW) satellite density profiles, from very dense ultrafaints to large, low density satellites such as Crater~II that appear to be larger their anticipated host dark matter haloes. We consider cold dark matter (CDM), warm dark matter (WDM, 3.3keV thermal relic power spectrum), and a self-interacting dark matter model (SIDM) that induces gravothermal collapse in low mass subhaloes. Predictions for these density profiles are complicated by the limitations of simulation resolution in the stripping of subhaloes by the MW system, therefore we make predictions for satellite properties in these three models using $N$-body simulations combined with a semi-analytic halo stripping algorithm. We find that most CDM and WDM subhaloes of mass $>10^{8}$$M_{\odot}$ are large enough after stripping to fit most satellites; however, the required amount of stripping often requires a stronger tidal field than is available on the subhalo's orbit. The lower concentrations of WDM subhaloes enable this model to explain the required satellite masses with less stripping than is necessary for CDM, and is thus consistent with orbits of larger pericentres. SIDM cores offer the best fits to massive, low density satellites at the expense of predicting many $>10^{9}$$M_{\odot}$ subhaloes to host low density satellites with no observed analogue. We conclude that an SIDM model must have a very high velocity-dependent cross-section in order to match all satellites, and that WDM offers a marginally better fit than CDM to the MW satellite mass function.
△ Less
Submitted 8 February, 2023; v1 submitted 14 September, 2022;
originally announced September 2022.
-
Differentiating small-scale subhalo distributions in CDM and WDM models using persistent homology
Authors:
Jessi Cisewski-Kehe,
Brittany Terese Fasy,
Wojciech Hellwing,
Mark R. Lovell,
Pawel Drozda,
Mike Wu
Abstract:
The spatial distribution of galaxies at sufficiently small scales will encode information about the identity of the dark matter. We develop a novel description of the halo distribution using persistent homology summaries, in which collections of points are decomposed into clusters, loops and voids. We apply these methods, together with a set of hypothesis tests, to dark matter haloes in MW-analog…
▽ More
The spatial distribution of galaxies at sufficiently small scales will encode information about the identity of the dark matter. We develop a novel description of the halo distribution using persistent homology summaries, in which collections of points are decomposed into clusters, loops and voids. We apply these methods, together with a set of hypothesis tests, to dark matter haloes in MW-analog environment regions of the cold dark matter (CDM) and warm dark matter (WDM) Copernicus Complexio $N$-body cosmological simulations. The results of the hypothesis tests find statistically significant differences (p-values $\leq$ 0.001) between the CDM and WDM structures, and the functional summaries of persistence diagrams detect differences at scales that are distinct from the comparison spatial point process functional summaries considered (including the two-point correlation function). The differences between the models are driven most strongly at filtration scales $\sim100$~kpc, where CDM generates larger numbers of unconnected halo clusters while WDM instead generates loops. This study was conducted on dark matter haloes generally; future work will involve applying the same methods to realistic galaxy catalogues.
△ Less
Submitted 1 April, 2022;
originally announced April 2022.
-
Could fresh lava be (warm) dark matter?
Authors:
Mark R. Lovell
Abstract:
Dark matter models can be classified according to their impact on the properties of galaxies, including cold dark matter (CDM), warm dark matter (WDM), self-interacting dark matter (SIDM) and fuzzy dark matter (FDM). In celebration of April Fool's Day, and also of the 1-year anniversary of the start of the 2022 volcanic eruption at Fagradalsfjall here in Iceland, we explore fresh lava as a candida…
▽ More
Dark matter models can be classified according to their impact on the properties of galaxies, including cold dark matter (CDM), warm dark matter (WDM), self-interacting dark matter (SIDM) and fuzzy dark matter (FDM). In celebration of April Fool's Day, and also of the 1-year anniversary of the start of the 2022 volcanic eruption at Fagradalsfjall here in Iceland, we explore fresh lava as a candidate for WDM specifically. We verify first hand that lava is indeed warm (exhibits free-streaming and retains temperature for several months after the eruption ends, is 1000K, sets fire to grass, can feel one's eyebrows singe at a distance of 4m) and dark once sufficiently decoupled from its source of production.
△ Less
Submitted 30 March, 2022;
originally announced March 2022.
-
Snowmass2021 Cosmic Frontier White Paper: Dark Matter Physics from Halo Measurements
Authors:
Keith Bechtol,
Simon Birrer,
Francis-Yan Cyr-Racine,
Katelin Schutz,
Susmita Adhikari,
Mustafa Amin,
Arka Banerjee,
Simeon Bird,
Nikita Blinov,
Kimberly K. Boddy,
Celine Boehm,
Kevin Bundy,
Malte Buschmann,
Sukanya Chakrabarti,
David Curtin,
Liang Dai,
Alex Drlica-Wagner,
Cora Dvorkin,
Adrienne L. Erickcek,
Daniel Gilman,
Saniya Heeba,
Stacy Kim,
Vid Iršič,
Alexie Leauthaud,
Mark Lovell
, et al. (19 additional authors not shown)
Abstract:
The non-linear process of cosmic structure formation produces gravitationally bound overdensities of dark matter known as halos. The abundances, density profiles, ellipticities, and spins of these halos can be tied to the underlying fundamental particle physics that governs dark matter at microscopic scales. Thus, macroscopic measurements of dark matter halos offer a unique opportunity to determin…
▽ More
The non-linear process of cosmic structure formation produces gravitationally bound overdensities of dark matter known as halos. The abundances, density profiles, ellipticities, and spins of these halos can be tied to the underlying fundamental particle physics that governs dark matter at microscopic scales. Thus, macroscopic measurements of dark matter halos offer a unique opportunity to determine the underlying properties of dark matter across the vast landscape of dark matter theories. This white paper summarizes the ongoing rapid development of theoretical and experimental methods, as well as new opportunities, to use dark matter halo measurements as a pillar of dark matter physics.
△ Less
Submitted 24 April, 2023; v1 submitted 14 March, 2022;
originally announced March 2022.
-
The feasibility of constraining DM interactions with high-redshift observations by JWST
Authors:
Ali Kurmus,
Sownak Bose,
Mark Lovell,
Francis-Yan Cyr-Racine,
Mark Vogelsberger,
Christoph Pfrommer,
Jesús Zavala
Abstract:
Observations of the high redshift universe provide a promising avenue for constraining the nature of the dark matter (DM). This will be even more true with the advent of the James Webb Space Telescope (JWST). We run cosmological simulations of galaxy formation as part of the Effective Theory of Structure Formation (ETHOS) project to compare high redshift galaxies in Cold (CDM) and alternative DM m…
▽ More
Observations of the high redshift universe provide a promising avenue for constraining the nature of the dark matter (DM). This will be even more true with the advent of the James Webb Space Telescope (JWST). We run cosmological simulations of galaxy formation as part of the Effective Theory of Structure Formation (ETHOS) project to compare high redshift galaxies in Cold (CDM) and alternative DM models which have varying relativistic coupling and self-interaction strengths. The interacting DM scenarios produce a cutoff in the linear power spectrum on small-scales, followed by a series of "dark acoustic oscillations". We find that DM interactions suppress the abundance of galaxies below $M_\star \sim 10^8\,M_\odot$ for the models considered. The cutoff in the power spectrum delays structure formation relative to CDM. Objects in ETHOS that end up at the same final masses as their CDM counterparts are characterised by a more vigorous phase of early star formation. While galaxies with $M_\star \lesssim 10^6\,M_\odot$ make up more than 60 per cent of star formation in CDM at $z\approx 10$, they contribute only about half the star formation density in ETHOS. These differences diminish with decreasing redshift. We find that the effects of DM self-interactions are negligible compared to effects of relativistic coupling (i.e. the effective initial conditions for galaxy formation) in all properties of the galaxy population we examine. Finally, we show that the clustering strength of galaxies at high redshifts depends sensitively on DM physics, although these differences are manifest on scales that may be too small to be measurable by JWST.
△ Less
Submitted 29 July, 2022; v1 submitted 9 March, 2022;
originally announced March 2022.
-
High and low Sérsic index bulges in Milky Way- and M31-like galaxies: origin and connection to the bar with TNG50
Authors:
Ignacio D. Gargiulo,
Antonela Monachesi,
Facundo A. Gómez,
Dylan Nelson,
Annalisa Pillepich,
Rüdiger Pakmor,
R. J. J. Grand,
Francesca Fragkoudi,
Lars Hernquist,
Mark Lovell,
Federico Marinacci
Abstract:
We study bulge formation in MW/M31-like galaxies in a $Λ$-cold dark matter scenario, focusing on the origin of high- and low-Sersic index bulges. For this purpose we use TNG50, a simulation of the IllustrisTNG project that combines a resolution of $\sim 8 \times 10^4 M_{\odot}$ in stellar particles with a cosmological volume 52 cMpc in extent. We parametrize bulge surface brightness profiles by th…
▽ More
We study bulge formation in MW/M31-like galaxies in a $Λ$-cold dark matter scenario, focusing on the origin of high- and low-Sersic index bulges. For this purpose we use TNG50, a simulation of the IllustrisTNG project that combines a resolution of $\sim 8 \times 10^4 M_{\odot}$ in stellar particles with a cosmological volume 52 cMpc in extent. We parametrize bulge surface brightness profiles by the Sérsic index and the bulge-to-total (B/T) ratio obtained from two-component photometric decompositions. In our sample of 287 MW/M31-like simulated galaxies, $17.1\%$ of photometric bulges exhibit high-Sérsic indices and $82.9\%$ show low-Sérsic indices. We study the impact that the environment, mergers and bars have in shaping the surface brightness profiles. We explore two different definitions for local environment and find no correlation between bulge properties and the environment where they reside. Simulated galaxies with higher Sérsic indices show, on average, a higher fraction of ex-situ stars in their kinematically selected bulges. For this bulge population the last significant merger (total mass ratio $m_{\rm sat}/m_{\rm host} > 0.1$) occurs, on average, at later times. However, a substantial fraction of low-Sérsic index bulges also experience a late significant merger. We find that bars play an important role in the development of the different types of photometric bulges. We show that the fraction of simulated galaxies with strong bars is smaller for the high- than for the low-Sérsic index population, reaching differences of $20\%$ at $z > 1$. Simulated galaxies with high fractions of ex-situ stars in the bulge do not develop strong bars. Conversely, simulated galaxies with long-lived strong bars have bulges with ex-situ fractions, $f_{\rm ex-situ} < 0.2$.
△ Less
Submitted 26 November, 2021;
originally announced November 2021.
-
The spatial distribution of Milky Way satellites, gaps in streams and the nature of dark matter
Authors:
Mark R. Lovell,
Marius Cautun,
Carlos S. Frenk,
Wojciech A. Hellwing,
Oliver Newton
Abstract:
The spatial distribution of Milky Way (MW) subhaloes provides an important set of observables for testing cosmological models. These include the radial distribution of luminous satellites, planar configurations, and the abundance of dark subhaloes whose existence or absence is key to distinguishing amongst dark matter models. We use the COCO $N$-body simulations of cold dark matter (CDM) and 3.3ke…
▽ More
The spatial distribution of Milky Way (MW) subhaloes provides an important set of observables for testing cosmological models. These include the radial distribution of luminous satellites, planar configurations, and the abundance of dark subhaloes whose existence or absence is key to distinguishing amongst dark matter models. We use the COCO $N$-body simulations of cold dark matter (CDM) and 3.3keV thermal relic warm dark matter (WDM) to predict the satellite spatial distribution in the limit that the impact of baryonic physics is minimal. We demonstrate that the radial distributions of CDM and 3.3keV-WDM luminous satellites are identical if the minimum pre-infall halo mass to form a galaxy is $>10^{8.5}\mathrm{M}_{\odot}$. The distribution of dark subhaloes is significantly more concentrated in WDM due to the absence of low mass, recently accreted substructures that typically inhabit the outer parts of a MW halo in CDM. We show that subhaloes of mass $[10^{7},10^{8}]\mathrm{M}_{\odot}$ and within 30kpc of the centre are the stripped remnants of larger haloes in both models. Therefore their abundance in WDM is $3\times$ higher than one would anticipate from the overall WDM subhalo population. We estimate that differences between CDM and WDM concentration-mass relations can be probed for subhalo-stream impact parameters $<2$kpc. Finally, we find that the impact of WDM on planes of satellites is likely negligible. Comprehensive comparisons with observations will require further work with high resolution, self-consistent hydrodynamical simulations.
△ Less
Submitted 8 February, 2023; v1 submitted 7 April, 2021;
originally announced April 2021.
-
Constraints on the properties of warm dark matter using the satellite galaxies of the Milky Way
Authors:
Oliver Newton,
Matteo Leo,
Marius Cautun,
Adrian Jenkins,
Carlos S. Frenk,
Mark R. Lovell,
John C. Helly,
Andrew J. Benson
Abstract:
The satellite galaxies of the Milky Way (MW) are effective probes of the underlying dark matter (DM) substructure, which is sensitive to the nature of the DM particle. In particular, a class of DM models have a power spectrum cut-off on the mass scale of dwarf galaxies and thus predict only small numbers of substructures below the cut-off mass. This makes the MW satellite system appealing to const…
▽ More
The satellite galaxies of the Milky Way (MW) are effective probes of the underlying dark matter (DM) substructure, which is sensitive to the nature of the DM particle. In particular, a class of DM models have a power spectrum cut-off on the mass scale of dwarf galaxies and thus predict only small numbers of substructures below the cut-off mass. This makes the MW satellite system appealing to constrain the DM properties: feasible models must produce enough substructure to host the number of observed Galactic satellites. Here, we compare theoretical predictions of the abundance of DM substructure in thermal relic warm DM (WDM) models with estimates of the total satellite population of the MW. This produces conservative robust lower limits on the allowed mass, $m_\mathrm{th}$, of the thermal relic WDM particle. As the abundance of satellite galaxies depends on the MW halo mass, we marginalize over the corresponding uncertainties and rule out $m_\mathrm{th} \leq 2.02\, \mathrm{keV}$ at 95 per cent confidence independently of assumptions about galaxy formation processes. Modelling some of these - in particular, the effect of reionization, which suppresses the formation of dwarf galaxies - strengthens our constraints on the DM properties and excludes models with $m_\mathrm{th} \leq 3.99\, \mathrm{keV}$ in our fiducial model. We also find that thermal relic models cannot produce enough satellites if the MW halo mass is $M_{200}\leq 0.6\times 10^{12}\, \mathrm{M_\odot}$, which imposes a lower limit on the MW halo mass in CDM. We address several observational and theoretical uncertainties and discuss how improvements in these will strengthen the DM mass constraints.
△ Less
Submitted 1 September, 2021; v1 submitted 17 November, 2020;
originally announced November 2020.
-
Joint constraints on thermal relic dark matter from strong gravitational lensing, the Lyman-$α$ forest, and Milky Way satellites
Authors:
Wolfgang Enzi,
Riccardo Murgia,
Oliver Newton,
Simona Vegetti,
Carlos Frenk,
Matteo Viel,
Marius Cautun,
Christopher D. Fassnacht,
Matt Auger,
Giulia Despali,
John McKean,
Léon V. E. Koopmans,
Mark Lovell
Abstract:
We derive joint constraints on the warm dark matter (WDM) half-mode scale by combining the analyses of a selection of astrophysical probes: strong gravitational lensing with extended sources, the Lyman-$α$ forest, and the number of luminous satellites in the Milky Way. We derive an upper limit of $λ_{\rm hm}=0.089{\rm~Mpc~h^{-1} }$ at the 95 per cent confidence level, which we show to be stable fo…
▽ More
We derive joint constraints on the warm dark matter (WDM) half-mode scale by combining the analyses of a selection of astrophysical probes: strong gravitational lensing with extended sources, the Lyman-$α$ forest, and the number of luminous satellites in the Milky Way. We derive an upper limit of $λ_{\rm hm}=0.089{\rm~Mpc~h^{-1} }$ at the 95 per cent confidence level, which we show to be stable for a broad range of prior choices. Assuming a Planck cosmology and that WDM particles are thermal relics, this corresponds to an upper limit on the half-mode mass of $M_{\rm hm }< 3 \times 10^{7} {\rm~M_{\odot}~h^{-1}}$, and a lower limit on the particle mass of $m_{\rm th }> 6.048 {\rm~keV}$, both at the 95 per cent confidence level. We find that models with $λ_{\rm hm}> 0.223 {\rm~Mpc~h^{-1} }$ (corresponding to $m_{\rm th }> 2.552 {\rm~keV}$ and $M_{\rm hm }< 4.8 \times 10^{8} {\rm~M_{\odot}~h^{-1}}$) are ruled out with respect to the maximum likelihood model by a factor $\leq 1/20$. For lepton asymmetries $L_6>10$, we rule out the $7.1 {\rm~keV}$ sterile neutrino dark matter model, which presents a possible explanation to the unidentified $3.55 {\rm~keV}$ line in the Milky Way and clusters of galaxies. The inferred 95 percentiles suggest that we further rule out the ETHOS-4 model of self-interacting DM. Our results highlight the importance of extending the current constraints to lower half-mode scales. We address important sources of systematic errors and provide prospects for how the constraints of these probes can be improved upon in the future.
△ Less
Submitted 22 June, 2021; v1 submitted 26 October, 2020;
originally announced October 2020.
-
The Onset of Gravothermal Core Collapse in Velocity Dependent Self-Interacting Dark Matter Subhaloes
Authors:
Hannah C. Turner,
Mark R. Lovell,
Jesús Zavala,
Mark Vogelsberger
Abstract:
It has been proposed that gravothermal collapse due to dark matter self-interactions (i.e. self-interacting dark matter, SIDM) can explain the observed diversity of the Milky Way (MW) satellites' central dynamical masses. We investigate the process behind this hypothesis using an $N$-body simulation of a MW-analogue halo with velocity dependent self-interacting dark matter (vdSIDM) in which the lo…
▽ More
It has been proposed that gravothermal collapse due to dark matter self-interactions (i.e. self-interacting dark matter, SIDM) can explain the observed diversity of the Milky Way (MW) satellites' central dynamical masses. We investigate the process behind this hypothesis using an $N$-body simulation of a MW-analogue halo with velocity dependent self-interacting dark matter (vdSIDM) in which the low velocity self-scattering cross-section, $σ_{T}/m_{x}$, reaches 100 cm$^{2}$g$^{-1}$; we dub this model the vd100 model. We compare the results of this simulation to simulations of the same halo that employ different dark models, including cold dark matter (CDM) and other, less extreme SIDM models. The masses of the vd100 haloes are very similar to their CDM counterparts, but the values of their maximum circular velocities, $V_{max}$, are significantly higher. We determine that these high $V_{max}$ subhaloes were objects in the mass range [$5\times10^{6}$, $1\times10^{8}$] $M_\odot$ at $z=1$ that undergo gravothermal core collapse. These collapsed haloes have density profiles that are described by single power laws down to the resolution limit of the simulation, and the inner slope of this density profile is approximately $-3$. Resolving the ever decreasing collapsed region is challenging, and tailored simulations will be required to model the runaway instability accurately at scales $<1$ kpc.
△ Less
Submitted 4 August, 2021; v1 submitted 6 October, 2020;
originally announced October 2020.
-
Cosmological Dark Matter: a Review (the April Fool Edition)
Authors:
M. R. Lovell
Abstract:
Evidence has continued to accumulate over the last few decades as to the existence and nature of dark matter. Depending on the particle candidate, the dark matter can exhibit one of several cosmologically defined models: hot dark matter, cold dark matter, warm dark matter, self-interacting dark matter, and fuzzy dark matter. In this paper I review the relevance and status of these models, whether…
▽ More
Evidence has continued to accumulate over the last few decades as to the existence and nature of dark matter. Depending on the particle candidate, the dark matter can exhibit one of several cosmologically defined models: hot dark matter, cold dark matter, warm dark matter, self-interacting dark matter, and fuzzy dark matter. In this paper I review the relevance and status of these models, whether it is possible for more than one of these models to each constitute some fraction of the dark matter, and discuss the prospects for determining if any of these models can successfully describe the properties and evolution of our own Universe.
△ Less
Submitted 9 April, 2020; v1 submitted 30 March, 2020;
originally announced March 2020.
-
Towards a general parametrization of the warm dark matter halo mass function
Authors:
Mark R. Lovell
Abstract:
Studies of flux anomalies statistics and perturbations in stellar streams have the potential to constrain models of warm dark matter (WDM), including sterile neutrinos. Producing these constraints requires a parametrization of the WDM mass function relative to that of the cold dark matter (CDM) equivalent. We use five WDM models with half-mode masses, $M_\mathrm{hm}=[1.3,35]\times10^{8}$~…
▽ More
Studies of flux anomalies statistics and perturbations in stellar streams have the potential to constrain models of warm dark matter (WDM), including sterile neutrinos. Producing these constraints requires a parametrization of the WDM mass function relative to that of the cold dark matter (CDM) equivalent. We use five WDM models with half-mode masses, $M_\mathrm{hm}=[1.3,35]\times10^{8}$~$M_{\odot}$, spread across simulations of the Local Group, lensing ellipticals and the $z=2$ universe, to generate such a parametrization: we fit parameters to a functional form for the WDM-to-CDM halo mass function ratio, $n_\mathrm{WDM}(M_{X})/n_\mathrm{CDM}(M_{X})$, of ($1+(αM_\mathrm{hm}/M_{X})^β)^γ$. For $M_{X}\equiv$ virial mass of central halos we obtain $α=2.3$, $β=0.8$, and $γ=-1.0$, and this fit is steeper than the extended Press-Schechter formalism predicts. For $M_{X}\equiv$ mass of subhalos we instead obtain $α=4.2$, $β=2.5$ and $γ=-0.2$; in both mass definitions the scatter is $\sim20$~per~cent. The second fit typically underestimates the relative abundance of $z=2$ WDM subhaloes at the tens of per cent level. We caution that robust constraints will require bespoke simulations and a careful definition of halo mass, particularly for subhalos of mass $<10^{8}M_{\odot}$.
△ Less
Submitted 1 June, 2020; v1 submitted 2 March, 2020;
originally announced March 2020.
-
Local Group star formation in warm and self-interacting dark matter cosmologies
Authors:
Mark R. Lovell,
Wojciech Hellwing,
Aaron Ludlow,
Jesús Zavala,
Andrew Robertson,
Azadeh Fattahi,
Carlos S. Frenk,
Jennifer Hardwick
Abstract:
The nature of the dark matter can affect the collapse time of dark matter haloes, and can therefore be imprinted in observables such as the stellar population ages and star formation histories of dwarf galaxies. In this paper we use high resolution hydrodynamical simulations of Local Group-analogue (LG) volumes in cold dark matter (CDM), sterile neutrino warm dark matter (WDM) and self-interacting…
▽ More
The nature of the dark matter can affect the collapse time of dark matter haloes, and can therefore be imprinted in observables such as the stellar population ages and star formation histories of dwarf galaxies. In this paper we use high resolution hydrodynamical simulations of Local Group-analogue (LG) volumes in cold dark matter (CDM), sterile neutrino warm dark matter (WDM) and self-interacting dark matter (SIDM) models with the EAGLE galaxy formation code to study how galaxy formation times change with dark matter model. We are able to identify the same haloes in different simulations, since they share the same initial density field phases. We find that the stellar mass of galaxies depends systematically on resolution, and can differ by as much as a factor of two in haloes of a given dark matter mass. The evolution of the stellar populations in SIDM is largely identical to that of CDM, but in WDM early star formation is instead suppressed. The time at which LG haloes can begin to form stars through atomic cooling is delayed by $\sim$200~Myr in WDM models compared to CDM. It will be necessary to measure stellar ages of old populations to a precision of better than 100~Myr, and to address degeneracies with the redshift of reionization -- and potentially other baryonic processes -- in order to use these observables to distinguish between dark matter models.
△ Less
Submitted 25 September, 2020; v1 submitted 25 February, 2020;
originally announced February 2020.
-
Exploiting flux ratio anomalies to probe warm dark matter in future large scale surveys
Authors:
David Harvey,
Wessel Valkenburg,
Amelie Tamone,
Alexey Boyarsky,
Frederic Courbin,
Mark Lovell
Abstract:
Flux ratio anomalies in strong gravitationally lensed quasars constitute a unique way to probe the abundance of non-luminous dark matter haloes, and hence the nature of dark matter. In this paper we identify double imaged quasars as a statistically efficient probe of dark matter, since they are 20 times more abundant than quadruply imaged quasars. Using N-body simulations that include realistic ba…
▽ More
Flux ratio anomalies in strong gravitationally lensed quasars constitute a unique way to probe the abundance of non-luminous dark matter haloes, and hence the nature of dark matter. In this paper we identify double imaged quasars as a statistically efficient probe of dark matter, since they are 20 times more abundant than quadruply imaged quasars. Using N-body simulations that include realistic baryonic feedback, we measure the full distribution of flux ratios in doubly imaged quasars for cold (CDM) and warm dark matter (WDM) cosmologies. Through this method, we fold in two key systematics - quasar variability and line-of-sight structures. We find that WDM cosmologies predict a ~6 per cent difference in the cumulative distribution functions of flux ratios relative to CDM, with CDM predicting many more small ratios. Finally, we estimate that ~600 doubly imaged quasars will need to be observed in order to be able to unambiguously discern between CDM and the two WDM models studied here. Such sample sizes will be easily within reach of future large scale surveys such as Euclid. In preparation for this survey data we require discerning the scale of the uncertainties in modelling lens galaxies and their substructure in simulations, plus a strong understanding of the selection function of observed lensed quasars.
△ Less
Submitted 4 December, 2019;
originally announced December 2019.
-
The halo mass function in alternative dark matter models
Authors:
M. R. Lovell
Abstract:
The claimed detection of large amounts of substructure in lensing flux anomalies, and in Milky Way stellar stream gaps statistics, has lead to a step change in constraints on simple warm dark matter models. In this study we compute predictions for the halo mass function both for these simple models and also for comprehensive particle physics models of sterile neutrinos and dark acoustic oscillatio…
▽ More
The claimed detection of large amounts of substructure in lensing flux anomalies, and in Milky Way stellar stream gaps statistics, has lead to a step change in constraints on simple warm dark matter models. In this study we compute predictions for the halo mass function both for these simple models and also for comprehensive particle physics models of sterile neutrinos and dark acoustic oscillations. We show that the mass function fit of Lovell et al. underestimates the number of haloes less massive than the half-mode mass, $M_\mathrm{hm}$ by a factor of 2, relative to the extended Press-Schechter (EPS) method. The alternative approach of applying EPS to the Viel et al. matter power spectrum fit instead suggests good agreement at $M_\mathrm{hm}$ relative to the comprehensive model matter power spectra results, although the number of haloes with mass $<M_\mathrm{hm}$ is still suppressed due to the absence of small scale power in the fitting function. Overall, we find that the number of dark matter haloes with masses $<10^{8}M_{\odot}$ predicted by competitive particle physics models is underestimated by a factor of $\sim2$ when applying popular fitting functions, although careful studies that follow the stripping and destruction of subhaloes will be required in order to draw robust conclusions.
△ Less
Submitted 19 August, 2020; v1 submitted 26 November, 2019;
originally announced November 2019.
-
Baryonic clues to the puzzling diversity of dwarf galaxy rotation curves
Authors:
Isabel M. E. Santos-Santos,
Julio F. Navarro,
Andrew Robertson,
Alejandro Benítez-Llambay,
Kyle A. Oman,
Mark R. Lovell,
Carlos S. Frenk,
Aaron D. Ludlow,
Azadeh Fattahi,
Adam Ritz
Abstract:
We use a compilation of disc galaxy rotation curves to assess the role of the luminous component ("baryons") in the rotation curve diversity problem. As in earlier work, we find that rotation curve shape correlates with baryonic surface density: high surface density galaxies have rapidly-rising rotation curves consistent with cuspy cold dark matter halos; slowly-rising rotation curves (characteris…
▽ More
We use a compilation of disc galaxy rotation curves to assess the role of the luminous component ("baryons") in the rotation curve diversity problem. As in earlier work, we find that rotation curve shape correlates with baryonic surface density: high surface density galaxies have rapidly-rising rotation curves consistent with cuspy cold dark matter halos; slowly-rising rotation curves (characteristic of galaxies with inner mass deficits or "cores") occur only in low surface density galaxies. The correlation, however, seems too weak to be the main driver of the diversity. In addition, dwarf galaxies exhibit a clear trend, from "cuspy" systems where baryons are unimportant in the inner mass budget to "cored" galaxies where baryons actually dominate. This trend constrains the various scenarios proposed to explain the diversity, such as (i) baryonic inflows and outflows during galaxy formation; (ii) dark matter self-interactions; (iii) variations in the baryonic mass structure coupled to rotation velocities through the "mass discrepancy-acceleration relation" (MDAR); or (iv) non-circular motions in gaseous discs. Together with analytical modeling and cosmological hydrodynamical simulations, our analysis shows that each of these scenarios has promising features, but none seems to fully account for the observed diversity. The MDAR, in particular, is inconsistent with the observed trend between rotation curve shape and baryonic importance; either the trend is caused by systematic errors in the data or the MDAR does not apply. The origin of the dwarf galaxy rotation curve diversity and its relation to the structure of cold dark matter halos remains an open issue.
△ Less
Submitted 28 April, 2020; v1 submitted 20 November, 2019;
originally announced November 2019.
-
The lensing properties of subhaloes in massive elliptical galaxies in sterile neutrino cosmologies
Authors:
Giulia Despali,
Mark Lovell,
Simona Vegetti,
Robert A. Crain,
Benjamin D. Oppenheimer
Abstract:
We use high-resolution hydrodynamical simulations run with the EAGLE model of galaxy formation to study the differences between the properties of - and subsequently the lensing signal from - subhaloes of massive elliptical galaxies at redshift 0.2, in Cold and Sterile Neutrino (SN) Dark matter models. We focus on the two 7 keV SN models that bracket the range of matter power spectra compatible wit…
▽ More
We use high-resolution hydrodynamical simulations run with the EAGLE model of galaxy formation to study the differences between the properties of - and subsequently the lensing signal from - subhaloes of massive elliptical galaxies at redshift 0.2, in Cold and Sterile Neutrino (SN) Dark matter models. We focus on the two 7 keV SN models that bracket the range of matter power spectra compatible with resonantly-produced SN as the source of the observed 3.5 keV line. We derive an accurate parametrisation for the subhalo mass function in these two SN models relative to CDM, as well as the subhalo spatial distribution, density profile, and projected number density and the dark matter fraction in subhaloes. We create mock lensing maps from the simulated haloes to study the differences in the lensing signal in the framework of subhalo detection. We find that subhalo convergence is well described by a log-normal distribution and that signal of subhaloes in the power spectrum is lower in SN models with respect to CDM, at a level of 10 to 80 per cent, depending on the scale. However, the scatter between different projections is large and might make the use of power-spectrum studies on the typical scales of current lensing images very difficult. Moreover, in the framework of individual detections through gravitational imaging a sample of ~30 lenses with an average sensitivity of M_sub=5x10^7M_sun would be required to discriminate between CDM and the considered sterile neutrino models.
△ Less
Submitted 29 October, 2019; v1 submitted 15 July, 2019;
originally announced July 2019.
-
The diverse dark matter density at sub-kiloparsec scales in Milky Way satellites:implications for the nature of dark matter
Authors:
Jesús Zavala,
Mark R. Lovell,
Mark Vogelsberger,
Jan D. Burger
Abstract:
Milky Way (MW) satellites reside within dark matter (DM) subhalos with a broad distribution of circular velocity profiles. This diversity is enhanced with the inclusion of ultra-faint satellites, which seemingly have very high DM densities, albeit with large systematic uncertainties. We argue that if confirmed, this large diversity in the MW satellite population poses a serious test for the struct…
▽ More
Milky Way (MW) satellites reside within dark matter (DM) subhalos with a broad distribution of circular velocity profiles. This diversity is enhanced with the inclusion of ultra-faint satellites, which seemingly have very high DM densities, albeit with large systematic uncertainties. We argue that if confirmed, this large diversity in the MW satellite population poses a serious test for the structure formation theory with possible implications for the DM nature. For the Cold Dark Matter model, the diversity might be a signature of the combined effects of subhalo tidal disruption by the MW disk and strong supernova feedback. For models with a dwarf-scale cutoff in the power spectrum, the diversity is a consequence of the lower abundance of dwarf-scale halos. This diversity is most challenging for Self-Interacting Dark Matter (SIDM) models with cross sections $σ/m_χ\gtrsim1~$cm$^2$g$^{-1}$ where subhalos have too low densities to explain the ultra-faint galaxies. We propose a novel solution to explain the diversity of MW satellites based on the gravothermal collapse of SIDM haloes. This solution requires a velocity-dependent cross section that predicts a bimodal distribution of cuspy dense (collapsed) subhaloes consistent with the ultra-faint satellites, and cored lower density subhaloes consistent with the brighter satellites.
△ Less
Submitted 30 August, 2019; v1 submitted 22 April, 2019;
originally announced April 2019.
-
Simulating the dark matter decay signal from the Perseus galaxy cluster
Authors:
Mark R. Lovell,
Dmytro Iakubovskyi,
David Barnes,
Sownak Bose,
Carlos S. Frenk,
Tom Theuns,
Wojciech A. Hellwing
Abstract:
The nearby Perseus galaxy cluster is a key target for indirect detection searches for decaying dark matter. We use the C-EAGLE simulations of galaxy clusters to predict the flux, width and shape of a dark matter decay line, paying particular attention to the unexplained 3.55keV line detected in the spectra of some galaxies and clusters, and the upcoming XRISM X-ray observatory mission. We show tha…
▽ More
The nearby Perseus galaxy cluster is a key target for indirect detection searches for decaying dark matter. We use the C-EAGLE simulations of galaxy clusters to predict the flux, width and shape of a dark matter decay line, paying particular attention to the unexplained 3.55keV line detected in the spectra of some galaxies and clusters, and the upcoming XRISM X-ray observatory mission. We show that the line width in C-EAGLE clusters similar to Perseus is typically [600-800]$\mathrm{kms^{-1}}$, and therefore narrower than the amplitude of the velocity dispersion of galaxies in the cluster. Haloes that are significantly disturbed can, however, exhibit galaxy velocity dispersions higher than $1000\mathrm{kms^{-1}}$, and in this case will show a large difference between the line profiles of on- and off-center observations. We show that the line profile is likely to be slightly asymmetric, but still well approximated by a Gaussian at the 10% level, and that the halo asymmetry can lead to fluxes that vary by a factor of two. In summary, we predict that, if the previously reported 3.55keV line detections do originate from dark matter decay, the XRISM mission will detect a line with a roughly Gaussian profile at a rest frame energy of 3.55keV, with a width $>600\mathrm{kms^{-1}}$ and flux approximately in the range $[4-9]\times10^{-8}\mathrm{counts/sec/cm^{2}}$.
△ Less
Submitted 27 March, 2019;
originally announced March 2019.
-
The IllustrisTNG Simulations: Public Data Release
Authors:
Dylan Nelson,
Volker Springel,
Annalisa Pillepich,
Vicente Rodriguez-Gomez,
Paul Torrey,
Shy Genel,
Mark Vogelsberger,
Ruediger Pakmor,
Federico Marinacci,
Rainer Weinberger,
Luke Kelley,
Mark Lovell,
Benedikt Diemer,
Lars Hernquist
Abstract:
We present the full public release of all data from the TNG50, TNG100 and TNG300 simulations of the IllustrisTNG project. IllustrisTNG is a suite of large volume, cosmological, gravo-magnetohydrodynamical simulations run with the moving-mesh code Arepo. TNG includes a comprehensive model for galaxy formation physics, and each TNG simulation self-consistently solves for the coupled evolution of dar…
▽ More
We present the full public release of all data from the TNG50, TNG100 and TNG300 simulations of the IllustrisTNG project. IllustrisTNG is a suite of large volume, cosmological, gravo-magnetohydrodynamical simulations run with the moving-mesh code Arepo. TNG includes a comprehensive model for galaxy formation physics, and each TNG simulation self-consistently solves for the coupled evolution of dark matter, cosmic gas, luminous stars, and supermassive blackholes from early time to the present day, z=0. Each of the flagship runs -- TNG50, TNG100, and TNG300 -- are accompanied by lower-resolution and dark-matter only counterparts, and we discuss scientific and numerical cautions and caveats relevant when using TNG. Full volume snapshots are available at 100 redshifts; halo and subhalo catalogs at each snapshot and merger trees are also released. The data volume now directly accessible online is ~1.1 PB, including 2,000 full volume snapshots and ~110,000 high time-resolution subbox snapshots. Data access and analysis examples are available in IDL, Python, and Matlab. We describe improvements and new functionality in the web-based API, including on-demand visualization and analysis of galaxies and halos, exploratory plotting of scaling relations and other relationships between galactic and halo properties, and a new JupyterLab interface. This provides an online, browser-based, near-native data analysis platform which supports user computation with fully local access to TNG data, alleviating the need to download large simulated datasets.
△ Less
Submitted 29 January, 2021; v1 submitted 13 December, 2018;
originally announced December 2018.
-
ETHOS - an effective theory of structure formation: formation of the first haloes and their stars
Authors:
Mark R. Lovell,
Jesús Zavala,
Mark Vogelsberger
Abstract:
A cutoff in the linear matter power spectrum at dwarf galaxy scales has been shown to affect the abundance, formation mechanism and age of dwarf haloes and their galaxies at high and low redshift. We use hydrodynamical simulations of galaxy formation within the ETHOS framework in a benchmark model that has such a cutoff, and that has been shown to be an alternative to the cold dark matter (CDM) mo…
▽ More
A cutoff in the linear matter power spectrum at dwarf galaxy scales has been shown to affect the abundance, formation mechanism and age of dwarf haloes and their galaxies at high and low redshift. We use hydrodynamical simulations of galaxy formation within the ETHOS framework in a benchmark model that has such a cutoff, and that has been shown to be an alternative to the cold dark matter (CDM) model that alleviates its dwarf-scale challenges. We show how galaxies in this model form differently to CDM on a halo-by-halo basis, at redshifts $z\ge6$. We show that ETHOS haloes at the half-mode mass scale form with 50~per~cent less mass than their CDM counterparts due to their later formation times, yet they retain more of their gas reservoir due to the different behaviour of gas and dark matter during the monolithic collapse of the first haloes in models with a galactic-scale cutoff. As a result, galaxies in ETHOS haloes near the cutoff scale grow rapidly between $z=10-6$ and by $z=6$ end up having very similar stellar masses, higher gas fractions and higher star formation rates relative to their CDM counterparts. We highlight these differences by making predictions for how the number of galaxies with old stellar populations is suppressed in ETHOS for both $z=6$ galaxies and for gas-poor Local Group fossil galaxies. Interestingly, we find an age gradient in ETHOS between galaxies that form in high and low density environments.
△ Less
Submitted 11 December, 2018;
originally announced December 2018.
-
The signal of decaying dark matter with hydrodynamical simulations
Authors:
Mark R. Lovell,
David Barnes,
Yannick Bahé,
Joop Schaye,
Matthieu Schaller,
Tom Theuns,
Sownak Bose,
Robert A. Crain,
Claudio dalla Vecchia,
Carlos S. Frenk,
Wojciech Hellwing,
Scott T. Kay,
Aaron D. Ludlow,
Richard G. Bower
Abstract:
Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate tha…
▽ More
Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the the dark matter decay line for Perseus-analogue galaxy clusters lies in the range 1300-1700~\kms. Therefore, the non-detection of the 3.55~keV line in the centre of the Perseus cluster by the {\it Hitomi} collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with {\it XMM-Newton} and with the planned X-ray telescopes {\it XRISM} and {\it ATHENA}. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter.
△ Less
Submitted 11 October, 2018;
originally announced October 2018.
-
The fraction of dark matter within galaxies from the IllustrisTNG simulations
Authors:
Mark R. Lovell,
Annalisa Pillepich,
Shy Genel,
Dylan Nelson,
Volker Springel,
Rüdiger Pakmor,
Federico Marinacci,
Rainer Weinberger,
Paul Torrey,
Mark Vogelsberger,
Adebusola Alabi,
Lars Hernquist
Abstract:
We use the IllustrisTNG (TNG) cosmological simulations to provide theoretical expectations for the dark matter mass fractions (DMFs) and circular velocity profiles of galaxies. TNG predicts flat circular velocity curves for $z = 0$ Milky Way (MW)-like galaxies beyond a few kpc from the galaxy centre, in better agreement with observational constraints than its predecessor, Illustris. TNG also predi…
▽ More
We use the IllustrisTNG (TNG) cosmological simulations to provide theoretical expectations for the dark matter mass fractions (DMFs) and circular velocity profiles of galaxies. TNG predicts flat circular velocity curves for $z = 0$ Milky Way (MW)-like galaxies beyond a few kpc from the galaxy centre, in better agreement with observational constraints than its predecessor, Illustris. TNG also predicts an enhancement of the dark matter mass within the 3D stellar half-mass radius ($r_\mathrm{half}$; $M_\mathrm{200c} = 10^{10}-10^{13}\mathrm{M}_{\odot}$, $z \le2$) compared to its dark matter only and Illustris counterparts. This enhancement leads TNG present-day galaxies to be dominated by dark matter within their inner regions, with $f_\mathrm{DM}(<r_\mathrm{half})\gtrsim0.5$ at all masses and with a minimum for MW-mass galaxies. The 1$σ$ scatter is $\lesssim$ 10~per~cent at all apertures, which is smaller than that inferred by some observational datasets, e.g. 40 per cent from the SLUGGS survey. TNG agrees with the majority of the observationally inferred values for elliptical galaxies once a consistent IMF is adopted (Chabrier) and the DMFs are measured within the same apertures. The DMFs measured within $r_\mathrm{half}$ increase towards lower redshifts: this evolution is dominated by the increase in galaxy size with time. At $z\sim2$, the DMF in disc-like TNG galaxies decreases with increasing galaxy mass, with $f_\mathrm{DM}(<r_\mathrm{half}) \sim 0.10-0.65$ for $10^{10} \lesssim M_{\rm stars}/\mathrm{M}_{\odot} \lesssim 10^{12}$, and are two times higher than if TNG galaxies resided in Navarro-Frenk-White dark matter haloes unaffected by baryonic physics. It remains to be properly assessed whether recent observational estimates of the DMFs at $z\sim2$ rule out the contraction of the dark matter haloes predicted by the TNG model.
△ Less
Submitted 15 January, 2019; v1 submitted 30 January, 2018;
originally announced January 2018.
-
Constraining sterile neutrino cosmologies with strong gravitational lensing observations at redshift z~0.2
Authors:
S. Vegetti,
G. Despali,
M. R. Lovell,
W. Enzi
Abstract:
We use the observed amount of subhaloes and line-of-sight dark matter haloes in a sample of 11 gravitational lens systems from the Sloan Lens ACS Survey to constrain the free-streaming properties of the dark matter particles. In particular, we combine the detection of a small-mass dark matter halo by Vegetti et al. 2010 with the non-detections by Vegetti et al. 2014 and compare the derived subhalo…
▽ More
We use the observed amount of subhaloes and line-of-sight dark matter haloes in a sample of 11 gravitational lens systems from the Sloan Lens ACS Survey to constrain the free-streaming properties of the dark matter particles. In particular, we combine the detection of a small-mass dark matter halo by Vegetti et al. 2010 with the non-detections by Vegetti et al. 2014 and compare the derived subhalo and halo mass functions with expectations from cold dark matter (CDM) and resonantly produced sterile neutrino models. We constrain the half-mode mass, i.e. the mass scale at which the linear matter power spectrum is reduced by 50 per cent relatively to the CDM model, to be $\log M_{\rm{hm}} \left[M_\odot\right] < 12.0$ (equivalent thermal relic mass $m_{\rm th} > 0.3$ keV) at the 2$σ$ level. This excludes sterile neutrino models with neutrino masses $m_{\rm s} < 0.8$ keV at any value of $L_{\rm 6}$. Our constraints are weaker than currently provided by the number of Milky Way satellites, observations of the 3.5 keV X-ray line, and the Lyman $α$ forest. However, they are more robust than the former as they are less affected by baryonic processes. Moreover, unlike the latter, they are not affected by assumptions on the thermal histories for the intergalactic medium. Gravitational lens systems with higher data quality and higher source and lens redshift are required to obtain tighter constraints.
△ Less
Submitted 9 October, 2018; v1 submitted 4 January, 2018;
originally announced January 2018.
-
ETHOS -- an effective theory of structure formation: predictions for the high-redshift Universe -- abundance of galaxies and reionization
Authors:
Mark R. Lovell,
Jesús Zavala,
Mark Vogelsberger,
Xuejian Shen,
Francis-Yan Cyr-Racine,
Christoph Pfrommer,
Kris Sigurdson,
Michael Boylan-Kolchin,
Annalisa Pillepich
Abstract:
We contrast predictions for the high-redshift galaxy population and reionization history between cold dark matter (CDM) and an alternative self-interacting dark matter model based on the recently developed ETHOS framework that alleviates the small-scale CDM challenges within the Local Group. We perform the highest resolution hydrodynamical cosmological simulations (a 36~Mpc$^3$ volume with gas cel…
▽ More
We contrast predictions for the high-redshift galaxy population and reionization history between cold dark matter (CDM) and an alternative self-interacting dark matter model based on the recently developed ETHOS framework that alleviates the small-scale CDM challenges within the Local Group. We perform the highest resolution hydrodynamical cosmological simulations (a 36~Mpc$^3$ volume with gas cell mass of $\sim10^5\mathrm{M}_{\odot}$ and minimum gas softening of $\sim180$~pc) within ETHOS to date -- plus a CDM counterpart -- to quantify the abundance of galaxies at high redshift and their impact on reionization. We find that ETHOS predicts galaxies with higher ultraviolet (UV) luminosities than their CDM counterparts and a faster build-up of the faint end of the UV luminosity function. These effects, however, make the optical depth to reionization less sensitive to the power spectrum cut-off: the ETHOS model differs from the CDM $τ$ value by only 10 per cent and is consistent with Planck limits if the effective escape fraction of UV photons is 0.1-0.5. We conclude that current observations of high-redshift luminosity functions cannot differentiate between ETHOS and CDM models, but deep JWST surveys of strongly-lensed, inherently faint galaxies have the potential to test non-CDM models that offer attractive solutions to CDM's Local Group problems.
△ Less
Submitted 9 January, 2019; v1 submitted 28 November, 2017;
originally announced November 2017.
-
The galaxy population in cold and warm dark matter cosmologies
Authors:
Lan Wang,
Violeta Gonzalez-Perez,
Lizhi Xie,
Andrew P. Cooper,
Carlos S. Frenk,
Liang Gao,
Wojciech A. Hellwing,
John Helly,
Mark R. Lovell,
Lilian Jiang
Abstract:
We use a pair of high resolution N-body simulations implementing two dark matter models, namely the standard cold dark matter (CDM) cosmogony and a warm dark matter (WDM) alternative where the dark matter particle is a 1.5keV thermal relic. We combine these simulations with the GALFORM semi-analytical galaxy formation model in order to explore differences between the resulting galaxy populations.…
▽ More
We use a pair of high resolution N-body simulations implementing two dark matter models, namely the standard cold dark matter (CDM) cosmogony and a warm dark matter (WDM) alternative where the dark matter particle is a 1.5keV thermal relic. We combine these simulations with the GALFORM semi-analytical galaxy formation model in order to explore differences between the resulting galaxy populations. We use GALFORM model variants for CDM and WDM that result in the same z=0 galaxy stellar mass function by construction. We find that most of the studied galaxy properties have the same values in these two models, indicating that both dark matter scenarios match current observational data equally well. Even in under-dense regions, where discrepancies in structure formation between CDM and WDM are expected to be most pronounced, the galaxy properties are only slightly different. The only significant difference in the local universe we find is in the galaxy populations of "Local Volumes", regions of radius 1 to 8Mpc around simulated Milky Way analogues. In such regions our WDM model provides a better match to observed local galaxy number counts and is five times more likely than the CDM model to predict sub-regions within them that are as empty as the observed Local Void. Thus, a highly complete census of the Local Volume and future surveys of void regions could provide constraints on the nature of dark matter.
△ Less
Submitted 4 April, 2017; v1 submitted 14 December, 2016;
originally announced December 2016.
-
Searching for Primordial Black Holes in the radio and X-ray sky
Authors:
Daniele Gaggero,
Gianfranco Bertone,
Francesca Calore,
Riley M. T. Connors,
Mark Lovell,
Sera Markoff,
Emma Storm
Abstract:
We model the accretion of gas onto a population of massive primordial black holes in the Milky Way, and compare the predicted radio and X-ray emission with observational data. We show that under conservative assumptions on the accretion process, the possibility that ${\cal O}(10) \, M_\odot$ primordial black holes can account for all of the dark matter in the Milky Way is excluded at $5σ$ by a com…
▽ More
We model the accretion of gas onto a population of massive primordial black holes in the Milky Way, and compare the predicted radio and X-ray emission with observational data. We show that under conservative assumptions on the accretion process, the possibility that ${\cal O}(10) \, M_\odot$ primordial black holes can account for all of the dark matter in the Milky Way is excluded at $5σ$ by a comparison with a VLA radio catalog at $1.4$ GHz, and at $\simeq 40σ$ by a comparison with a Chandra X-ray catalog ($0.5 - 8$ keV). We argue that this method can be used to identify such a population of primordial black holes with more sensitive future radio and X-ray surveys.
△ Less
Submitted 28 June, 2017; v1 submitted 1 December, 2016;
originally announced December 2016.
-
Properties of Local Group galaxies in hydrodynamical simulations of sterile neutrino dark matter cosmologies
Authors:
Mark R. Lovell,
Sownak Bose,
Alexey Boyarsky,
Robert A. Crain,
Carlos S. Frenk,
Wojciech A. Hellwing,
Aaron D. Ludlow,
Julio F. Navarro,
Oleg Ruchayskiy,
Till Sawala,
Matthieu Schaller,
Joop Schaye,
Tom Theuns
Abstract:
We study galaxy formation in sterile neutrino dark matter models that differ signifi- cantly from both cold and from `warm thermal relic' models. We use the EAGLE code to carry out hydrodynamic simulations of the evolution of pairs of galaxies chosen to resemble the Local Group, as part of the APOSTLE simulations project. We compare cold dark matter (CDM) with two sterile neutrino models with 7 ke…
▽ More
We study galaxy formation in sterile neutrino dark matter models that differ signifi- cantly from both cold and from `warm thermal relic' models. We use the EAGLE code to carry out hydrodynamic simulations of the evolution of pairs of galaxies chosen to resemble the Local Group, as part of the APOSTLE simulations project. We compare cold dark matter (CDM) with two sterile neutrino models with 7 keV mass: one, the warmest among all models of this mass (LA120) and the other, a relatively cold case (LA10). We show that the lower concentration of sterile neutrino subhalos compared to their CDM counterparts makes the inferred inner dark matter content of galaxies like Fornax (or Magellanic Clouds) less of an outlier in the sterile neutrino cosmologies. In terms of the galaxy number counts the LA10 simulations are indistinguishable from CDM when one takes into account halo-to-halo (or `simulation-to-simulation') scatter. In order for the LA120 model to match the number of Local Group dwarf galaxies, a higher fraction of low mass haloes is required to form galaxies than is predicted by the EAGLE simulations. As the census of the Local Group galaxies nears completion, this population may provide a strong discriminant between cold and warm dark matter models.
△ Less
Submitted 22 August, 2017; v1 submitted 31 October, 2016;
originally announced November 2016.
-
Addressing the too big to fail problem with baryon physics and sterile neutrino dark matter
Authors:
Mark R. Lovell,
Violeta Gonzalez-Perez,
Sownak Bose,
Alexey Boyarsky,
Shaun Cole,
Carlos S. Frenk,
Oleg Ruchayskiy
Abstract:
N-body dark matter simulations of structure formation in the $Λ$CDM model predict a population of subhalos within Galactic halos that have higher central densities than inferred for satellites of the Milky Way, a tension known as the `too big to fail' problem. Proposed solutions include baryonic effects, a smaller mass for the Milky Way halo, and warm dark matter. We test these three possibilities…
▽ More
N-body dark matter simulations of structure formation in the $Λ$CDM model predict a population of subhalos within Galactic halos that have higher central densities than inferred for satellites of the Milky Way, a tension known as the `too big to fail' problem. Proposed solutions include baryonic effects, a smaller mass for the Milky Way halo, and warm dark matter. We test these three possibilities using a semi-analytic model of galaxy formation to generate luminosity functions for Milky Way halo-analogue satellite populations, the results of which are then coupled to the Jiang & van den Bosch model of subhalo stripping to predict the subhalo $V_\mathrm{max}$ functions for the 10 brightest satellites. We find that selecting the brightest satellites (as opposed to the most massive) and modelling the expulsion of gas by supernovae at early times increases the likelihood of generating the observed Milky Way satellite $V_\mathrm{max}$ function. The preferred halo mass is $6\times10^{11}M_{\odot}$, which has a 14 percent probability to host a $V_\mathrm{max}$ function like that of the Milky Way satellites. This probability is reduced to 8 percent for a $1.0\times10^{12}M_{\odot}$ halo and to 3 percent for a $1.4\times10^{12}M_{\odot}$ halo. We conclude that the Milky Way satellite $V_\mathrm{max}$ function is compatible with a CDM cosmology, as previously found by Sawala et al. using hydrodynamic simulations. Sterile neutrino-warm dark matter models achieve a higher degree of agreement with the observations, with a maximum 35 percent chance of generating the observed Milky Way satellite $V_\mathrm{max}$ function. However, more work is required to check that the semi-analytic stripping model is calibrated correctly in the sterile neutrino cosmology, and to check if our sterile neutrino models produce sufficient numbers of faint satellites.
△ Less
Submitted 22 August, 2017; v1 submitted 31 October, 2016;
originally announced November 2016.
-
Reionisation in sterile neutrino cosmologies
Authors:
Sownak Bose,
Carlos S. Frenk,
Hou Jun,
Cedric G. Lacey,
Mark R. Lovell
Abstract:
We investigate the process of reionisation in a model in which the dark matter is a warm elementary particle such as a sterile neutrino. We focus on models that are consistent with the dark matter decay interpretation of the recently detected line at 3.5 keV in the X-ray spectra of galaxies and clusters. In warm dark matter models the primordial spectrum of density perturbations has a cut-off on t…
▽ More
We investigate the process of reionisation in a model in which the dark matter is a warm elementary particle such as a sterile neutrino. We focus on models that are consistent with the dark matter decay interpretation of the recently detected line at 3.5 keV in the X-ray spectra of galaxies and clusters. In warm dark matter models the primordial spectrum of density perturbations has a cut-off on the scale of dwarf galaxies. Structure formation therefore begins later than in the standard cold dark matter (CDM) model and very few objects form below the cut-off mass scale. To calculate the number of ionising photons, we use the Durham semi-analytic model of galaxy formation, GALFORM. We find that even the most extreme 7 keV sterile neutrino we consider is able to reionise the Universe early enough to be compatible with the bounds on the epoch of reionisation from Planck. This, perhaps surprising, result arises from the rapid build-up of high redshift galaxies in the sterile neutrino models which is also reflected in a faster evolution of their far-UV luminosity function between $10>z>7$ than in CDM. The dominant sources of ionising photons are systematically more massive in the sterile neutrino models than in CDM. As a consistency check on the models, we calculate the present-day luminosity function of satellites of Milky Way-like galaxies. When the satellites recently discovered in the DES survey are taken into account, strong constraints are placed on viable sterile neutrino models.
△ Less
Submitted 10 May, 2016;
originally announced May 2016.
-
Substructure and galaxy formation in the Copernicus Complexio warm dark matter simulations
Authors:
Sownak Bose,
Wojciech A. Hellwing,
Carlos S. Frenk,
Adrian Jenkins,
Mark R. Lovell,
John C. Helly,
Baojiu Li,
Violeta Gonzalez-Perez,
Liang Gao
Abstract:
We use the Copernicus Complexio (COCO) high resolution $N$-body simulations to investigate differences in the properties of small-scale structures in the standard cold dark matter (CDM) model and in a model with a cutoff in the initial power spectrum of density fluctuations consistent with both a thermally produced warm dark matter (WDM) particle or a sterile neutrino with mass 7 keV and leptogene…
▽ More
We use the Copernicus Complexio (COCO) high resolution $N$-body simulations to investigate differences in the properties of small-scale structures in the standard cold dark matter (CDM) model and in a model with a cutoff in the initial power spectrum of density fluctuations consistent with both a thermally produced warm dark matter (WDM) particle or a sterile neutrino with mass 7 keV and leptogenesis parameter $L_6=8.7$. The latter corresponds to the "coldest" model with this sterile neutrino mass compatible with the identification of the recently detected 3.5 keV X-ray line as resulting from particle decay. CDM and WDM predict very different number densities of subhaloes with mass $\leq 10^9\,h^{-1}\,M_\odot$ although they predict similar, nearly universal, normalised subhalo radial density distributions. Haloes and subhaloes in both models have cuspy NFW profiles, but WDM subhaloes below the cutoff scale in the power spectrum (corresponding to maximum circular velocities $V_{\mathrm{max}}^{z=0} \leq50~\mathrm{kms}^{-1}$) are less concentrated than their CDM counterparts. We make predictions for observable properties using the GALFORM semi-analytic model of galaxy formation. Both models predict Milky Way satellite luminosity functions consistent with observations, although the WDM model predicts fewer very faint satellites. This model, however, predicts slightly more UV bright galaxies at redshift $z>7$ than CDM, but both are consistent with observations. Gravitational lensing offers the best prospect of distinguishing between the models.
△ Less
Submitted 26 October, 2016; v1 submitted 25 April, 2016;
originally announced April 2016.
-
Simulated Milky Way analogues: implications for dark matter direct searches
Authors:
Nassim Bozorgnia,
Francesca Calore,
Matthieu Schaller,
Mark Lovell,
Gianfranco Bertone,
Carlos S. Frenk,
Robert A. Crain,
Julio F. Navarro,
Joop Schaye,
Tom Theuns
Abstract:
We study the implications of galaxy formation on dark matter direct detection using high resolution hydrodynamic simulations of Milky Way-like galaxies simulated within the EAGLE and APOSTLE projects. We identify Milky Way analogues that satisfy observational constraints on the Milky Way rotation curve and total stellar mass. We then extract the dark matter density and velocity distribution in the…
▽ More
We study the implications of galaxy formation on dark matter direct detection using high resolution hydrodynamic simulations of Milky Way-like galaxies simulated within the EAGLE and APOSTLE projects. We identify Milky Way analogues that satisfy observational constraints on the Milky Way rotation curve and total stellar mass. We then extract the dark matter density and velocity distribution in the Solar neighbourhood for this set of Milky Way analogues, and use them to analyse the results of current direct detection experiments. For most Milky Way analogues, the event rates in direct detection experiments obtained from the best fit Maxwellian distribution (with peak speed of 223 - 289 km/s) are similar to those obtained directly from the simulations. As a consequence, the allowed regions and exclusion limits set by direct detection experiments in the dark matter mass and spin-independent cross section plane shift by a few GeV compared to the Standard Halo Model, at low dark matter masses. For each dark matter mass, the halo-to-halo variation of the local dark matter density results in an overall shift of the allowed regions and exclusion limits for the cross section. However, the compatibility of the possible hints for a dark matter signal from DAMA and CDMS-Si and null results from LUX and SuperCDMS is not improved.
△ Less
Submitted 10 May, 2016; v1 submitted 18 January, 2016;
originally announced January 2016.
-
Satellite galaxies in semi-analytic models of galaxy formation with sterile neutrino dark matter
Authors:
Mark R. Lovell,
Sownak Bose,
Alexey Boyarsky,
Shaun Cole,
Carlos S. Frenk,
Violeta Gonzalez-Perez,
Rachel Kennedy,
Oleg Ruchayskiy,
Alex Smith
Abstract:
The sterile neutrino is a viable dark matter candidate that can be produced in the early Universe via non-equilibrium processes, and would therefore possess a highly non-thermal spectrum of primordial velocities. In this paper we analyse the process of structure formation with this class of dark matter particles. To this end we construct primordial dark matter power spectra as a function of the le…
▽ More
The sterile neutrino is a viable dark matter candidate that can be produced in the early Universe via non-equilibrium processes, and would therefore possess a highly non-thermal spectrum of primordial velocities. In this paper we analyse the process of structure formation with this class of dark matter particles. To this end we construct primordial dark matter power spectra as a function of the lepton asymmetry, $L_6$, that is present in the primordial plasma and leads to resonant sterile neutrino production. We compare these power spectra with those of thermally produced dark matter particles and show that resonantly produced sterile neutrinos are much colder than their thermal relic counterparts. We also demonstrate that the shape of these power spectra is not determined by the free-streaming scale alone. We then use the power spectra as an input for semi-analytic models of galaxy formation in order to predict the number of luminous satellite galaxies in a Milky Way-like halo. By assuming that the mass of the Milky Way halo must be no more than $2\times10^{12}M_{\odot}$ (the adopted upper bound based on current astronomical observations) we are able to constrain the value of $L_6$ for $M_s\le 8$~keV. We also show that the range of $L_6$ that is in best agreement with the 3.5~keV line (if produced by decays of 7~keV sterile neutrino) requires that the Milky Way halo has a mass no smaller than $1.5\times10^{12}M_{\odot}$. Finally, we compare the power spectra obtained by direct integration of the Boltzmann equations for a non-resonantly produced sterile neutrino with the fitting formula of Viel~et~al. and find that the latter significantly underestimates the power amplitude on scales relevant to satellite galaxies.
△ Less
Submitted 7 November, 2016; v1 submitted 12 November, 2015;
originally announced November 2015.
-
Mass assembly history and infall time of the Fornax dwarf spheroidal galaxy
Authors:
Mei-Yu Wang,
Louis E. Strigari,
Mark R. Lovell,
Carlos S. Frenk,
Andrew R. Zentner
Abstract:
We use cosmological simulations to identify dark matter subhalo host candidates of the Fornax dwarf spheroidal galaxy using the stellar kinematic properties of Fornax. We consider cold dark matter (CDM), warm dark matter (WDM), and decaying dark matter (DDM) simulations for our models of structure formation. The subhalo candidates in CDM typically have smaller mass and higher concentrations at z =…
▽ More
We use cosmological simulations to identify dark matter subhalo host candidates of the Fornax dwarf spheroidal galaxy using the stellar kinematic properties of Fornax. We consider cold dark matter (CDM), warm dark matter (WDM), and decaying dark matter (DDM) simulations for our models of structure formation. The subhalo candidates in CDM typically have smaller mass and higher concentrations at z = 0 than the corresponding candidates in WDM and DDM. We examine the formation histories of the ~ 100 Fornax candidate subhalos identified in CDM simulations and, using approximate luminosity-mass relationships for subhalos, we find two of these subhalos that are consistent with both the Fornax luminosity and kinematics. These two subhalos have a peak mass over ten times larger than their z = 0 mass. We suggest that in CDM the dark matter halo hosting Fornax must have been severely stripped of mass and that it had an infall time into the Milky Way of ~ 9 Gyr ago. In WDM, we find that candidate subhalos consistent with the properties of Fornax have a similar infall time and a similar degree of mass loss, while in DDM we find a later infall time of ~ 3 - 4 Gyr ago and significantly less mass loss. We discuss these results in the context of the Fornax star formation history, and show that these predicted subhalo infall times can be linked to different star formation quenching mechanisms. This emphasizes the links between the properties of the dark matter and the mechanisms that drive galaxy evolution.
△ Less
Submitted 22 September, 2015; v1 submitted 14 September, 2015;
originally announced September 2015.
-
Simulated Milky Way analogues: implications for dark matter indirect searches
Authors:
F. Calore,
N. Bozorgnia,
M. Lovell,
G. Bertone,
M. Schaller,
C. S. Frenk,
R. A. Crain,
J. Schaye,
T. Theuns,
J. W. Trayford
Abstract:
We study high-resolution hydrodynamic simulations of Milky Way type galaxies obtained within the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project, and identify the those that best satisfy observational constraints on the Milky Way total stellar mass, rotation curve, and galaxy shape. Contrary to mock galaxies selected on the basis of their total virial mass, the Milky Wa…
▽ More
We study high-resolution hydrodynamic simulations of Milky Way type galaxies obtained within the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project, and identify the those that best satisfy observational constraints on the Milky Way total stellar mass, rotation curve, and galaxy shape. Contrary to mock galaxies selected on the basis of their total virial mass, the Milky Way analogues so identified consistently exhibit very similar dark matter profiles inside the solar circle, therefore enabling more accurate predictions for indirect dark matter searches. We find in particular that high resolution simulated haloes satisfying observational constraints exhibit, within the inner few kiloparsecs, dark matter profiles shallower than those required to explain the so-called Fermi GeV excess via dark matter annihilation.
△ Less
Submitted 4 January, 2016; v1 submitted 7 September, 2015;
originally announced September 2015.
-
The COpernicus COmplexio: Statistical Properties of Warm Dark Matter Haloes
Authors:
Sownak Bose,
Wojciech A. Hellwing,
Carlos S. Frenk,
Adrian Jenkins,
Mark R. Lovell,
John C. Helly,
Baojiu Li
Abstract:
The recent detection of a 3.5 keV X-ray line from the centres of galaxies and clusters by Bulbul et al. (2014a) and Boyarsky et al. (2014a) has been interpreted as emission from the decay of 7 keV sterile neutrinos which could make up the (warm) dark matter (WDM). As part of the COpernicus COmplexio (COCO) programme, we investigate the properties of dark matter haloes formed in a high-resolution c…
▽ More
The recent detection of a 3.5 keV X-ray line from the centres of galaxies and clusters by Bulbul et al. (2014a) and Boyarsky et al. (2014a) has been interpreted as emission from the decay of 7 keV sterile neutrinos which could make up the (warm) dark matter (WDM). As part of the COpernicus COmplexio (COCO) programme, we investigate the properties of dark matter haloes formed in a high-resolution cosmological $N$-body simulation from initial conditions similar to those expected in a universe in which the dark matter consists of 7 keV sterile neutrinos. This simulation and its cold dark matter (CDM) counterpart have $\sim13.4$bn particles, each of mass $\sim 10^5\, h^{-1} M_\odot$, providing detailed information about halo structure and evolution down to dwarf galaxy mass scales. Non-linear structure formation on small scales ($M_{200}\, \leq\, 2 \times 10^9\,h^{-1}\,M_\odot$) begins slightly later in COCO-Warm than in COCO-Cold. The halo mass function at the present day in the WDM model begins to drop below its CDM counterpart at a mass $\sim 2 \times 10^{9}\,h^{-1}\,M_\odot$ and declines very rapidly towards lower masses so that there are five times fewer haloes of mass $M_{200}= 10^{8}\,h^{-1}\,M_\odot$ in COCO-Warm than in COCO-Cold. Halo concentrations on dwarf galaxy scales are correspondingly smaller in COCO-Warm, and we provide a simple functional form that describes its evolution with redshift. The shapes of haloes are similar in the two cases, but the smallest haloes in COCO-Warm rotate slightly more slowly than their CDM counterparts.
△ Less
Submitted 9 July, 2015; v1 submitted 7 July, 2015;
originally announced July 2015.
-
Decaying dark matter: the case for a deep X-ray observation of Draco
Authors:
Mark R. Lovell,
Gianfranco Bertone,
Alexey Boyarsky,
Adrian Jenkins,
Oleg Ruchayskiy
Abstract:
Recent studies of M31, the Galactic centre, and galaxy clusters have made tentative detections of an X-ray line at ~3.5 keV that could be produced by decaying dark matter. We use high resolution simulations of the Aquarius project to predict the likely amplitude of the X-ray decay flux observed in the GC relative to that observed in M31, and also of the GC relative to other parts of the Milky Way…
▽ More
Recent studies of M31, the Galactic centre, and galaxy clusters have made tentative detections of an X-ray line at ~3.5 keV that could be produced by decaying dark matter. We use high resolution simulations of the Aquarius project to predict the likely amplitude of the X-ray decay flux observed in the GC relative to that observed in M31, and also of the GC relative to other parts of the Milky Way halo and to dwarf spheroidal galaxies. We show that the reported detections from M31 and Andromeda are compatible with each other, and with upper limits arising from high galactic latitude observations, and imply a decay time τ ~10^28 seconds. We argue that this interpretation can be tested with deep observations of dwarf spheroidal galaxies: in 95 per cent of our mock observations, a 1.3 Msec pointed observation of Draco with XMM-Newton will enable us to discover or rule out at the 3σ level an X-ray feature from dark matter decay at 3.5 keV, for decay times τ < 0.8x10^28 sec.
△ Less
Submitted 20 November, 2015; v1 submitted 2 November, 2014;
originally announced November 2014.
-
The properties of warm dark matter haloes
Authors:
Mark R. Lovell,
Carlos S. Frenk,
Vincent R. Eke,
Adrian Jenkins,
Liang Gao,
Tom Theuns
Abstract:
Well-motivated elementary particle candidates for the dark matter, such as the sterile neutrino, behave as warm dark matter (WDM).For particle masses of order a keV, free streaming produces a cutoff in the linear fluctuation power spectrum at a scale corresponding to dwarf galaxies. We investigate the abundance and structure of WDM haloes and subhaloes on these scales using high resolution cosmolo…
▽ More
Well-motivated elementary particle candidates for the dark matter, such as the sterile neutrino, behave as warm dark matter (WDM).For particle masses of order a keV, free streaming produces a cutoff in the linear fluctuation power spectrum at a scale corresponding to dwarf galaxies. We investigate the abundance and structure of WDM haloes and subhaloes on these scales using high resolution cosmological N-body simulations of galactic haloes of mass similar to the Milky Way's. On scales larger than the free-streaming cutoff, the initial conditions have the same power spectrum and phases as one of the cold dark matter (CDM) haloes previously simulated by Springel et al as part of the Virgo consortium Aquarius project. We have simulated four haloes with WDM particle masses in the range 1.4-2.3keV and, for one case, we have carried out further simulations at varying resolution. N-body simulations in which the power spectrum cutoff is resolved are known to undergo artificial fragmentation in filaments producing spurious clumps which, for small masses (<10^7Msun in our case) outnumber genuine haloes. We have developed a robust algorithm to identify these spurious objects and remove them from our halo catalogues. We find that the WDM subhalo mass function is suppressed by well over an order magnitude relative to the CDM case for masses <10^9Msun. Requiring that there should be at least as many subhaloes as there are observed satellites in the Milky Way leads to a conservative lower limit to the (thermal equivalent) WDM particle mass of ~1.5\rmn{keV}. WDM haloes and subhaloes have cuspy density distributions that are well described by NFW or Einasto profiles. Their central densities are lower for lower WDM particle masses and none of the models we have considered suffer from the "too big to fail" problem recently highlighted by Boylan-Kolchin et al.
△ Less
Submitted 19 May, 2014; v1 submitted 6 August, 2013;
originally announced August 2013.
-
Dark Matter Halo Merger Histories Beyond Cold Dark Matter: I - Methods and Application to Warm Dark Matter
Authors:
Andrew J. Benson,
Arya Farahi,
Shaun Cole,
Leonidas A. Moustakas,
Adrian Jenkins,
Mark Lovell,
Rachel Kennedy,
John Helly,
Carlos Frenk
Abstract:
We describe a methodology to accurately compute halo mass functions, progenitor mass functions, merger rates and merger trees in non-cold dark matter universes using a self-consistent treatment of the generalized extended Press-Schechter formalism. Our approach permits rapid exploration of the subhalo population of galactic halos in dark matter models with a variety of different particle propertie…
▽ More
We describe a methodology to accurately compute halo mass functions, progenitor mass functions, merger rates and merger trees in non-cold dark matter universes using a self-consistent treatment of the generalized extended Press-Schechter formalism. Our approach permits rapid exploration of the subhalo population of galactic halos in dark matter models with a variety of different particle properties or universes with rolling, truncated, or more complicated power spectra. We make detailed comparisons of analytically derived mass functions and merger histories with recent warm dark matter cosmological N-body simulations, and find excellent agreement. We show that, once the accretion of smoothly distributed matter is accounted for, coarse-grained statistics such as the mass accretion history of halos can be almost indistinguishable between cold and warm dark matter cases. However, the halo mass function and progenitor mass functions differ significantly, with the warm dark matter cases being strongly suppressed below the free-streaming scale of the dark matter. We demonstrate the importance of using the correct solution for the excursion set barrier first-crossing distribution in warm dark matter - if the solution for a flat barrier is used instead the truncation of the halo mass function is much slower, leading to an overestimate of the number of low mass halos.
△ Less
Submitted 7 October, 2012; v1 submitted 13 September, 2012;
originally announced September 2012.
-
The haloes of bright satellite galaxies in a warm dark matter universe
Authors:
Mark Lovell,
Vincent Eke,
Carlos Frenk,
Liang Gao,
Adrian Jenkins,
Tom Theuns,
Jie Wang,
Simon White,
Alexey Boyarsky,
Oleg Ruchayskiy
Abstract:
High resolution N-body simulations of galactic cold dark matter haloes indicate that we should expect to find a few satellite galaxies around the Milky Way whose haloes have a maximum circular velocity in excess of 40 kms. Yet, with the exception of the Magellanic Clouds and the Sagittarius dwarf, which likely reside in subhaloes with significantly larger velocities than this, the bright satellite…
▽ More
High resolution N-body simulations of galactic cold dark matter haloes indicate that we should expect to find a few satellite galaxies around the Milky Way whose haloes have a maximum circular velocity in excess of 40 kms. Yet, with the exception of the Magellanic Clouds and the Sagittarius dwarf, which likely reside in subhaloes with significantly larger velocities than this, the bright satellites of the Milky Way all appear to reside in subhaloes with maximum circular velocities below 40 kms. As recently highlighted by Boylan-Kolchin et al., this discrepancy implies that the majority of the most massive subhaloes within a cold dark matter galactic halo are much too concentrated to be consistent with the kinematic data for the bright Milky Way satellites. Here we show that no such discrepancy exists if haloes are made of warm, rather than cold dark matter because these haloes are less concentrated on account of their typically later formation epochs. Warm dark matter is one of several possible explanations for the observed kinematics of the satellites.
△ Less
Submitted 16 January, 2012; v1 submitted 14 April, 2011;
originally announced April 2011.
-
The Link between Galactic Satellite Orbits and Subhalo Accretion
Authors:
Mark Lovell,
Vincent Eke,
Carlos Frenk,
Adrian Jenkins
Abstract:
We calculate the orbital angular momentum of dark matter subhaloes in the Aquarius simulations of cold dark matter galactic haloes. We calculate the orientation of their angular momentum relative to that of the spin vector of their host halo and find a variety of different configurations. All six Aquarius haloes contain statistically significant populations of subhalo orbits that are aligned with…
▽ More
We calculate the orbital angular momentum of dark matter subhaloes in the Aquarius simulations of cold dark matter galactic haloes. We calculate the orientation of their angular momentum relative to that of the spin vector of their host halo and find a variety of different configurations. All six Aquarius haloes contain statistically significant populations of subhalo orbits that are aligned with the main halo spin. All haloes posses a population of subhaloes that rotates in the same direction as the main halo and three of them possess, in addition, a population that rotates in the opposite direction. These configurations arise from the filamentary accretion of subhaloes. Quasi-planar distributions of coherently rotating satellites, such as those inferred in the Milky Way and other galaxies, arise naturally in simulations of a $Λ$CDM universe.
△ Less
Submitted 18 April, 2011; v1 submitted 3 August, 2010;
originally announced August 2010.