A critical reassessment of particle Dark Matter limits from dwarf satellites

Dwarf satellite galaxies are ideal laboratories for identifying particle Dark Matter signals. When setting limits on particle Dark Matter properties from null searches, it becomes however crucial the level at which the Dark Matter density profile within these systems is constrained by observations. In the limit in which the spherical Jeans equation is assumed to be valid for a given tracer stellar population, we study the solution of this equation having the Dark Matter mass profile as an output rather than as a trial parametric input. Within our new formulation, we address to what level dwarf spheroidal galaxies feature a reliable mass estimator. We assess then possible extrapolation of the density profiles in the inner regions and -- keeping explicit the dependence on the orbital anisotropy profile of the tracer population -- we derive general trends on the line-of-sight integral of the density profile squared, a quantity commonly dubbed J-factor and crucial to estimate fluxes from prompt Dark Matter pair annihilations. Taking Ursa Minor as a study case among Milky Way satellites, we perform Bayesian inference using the available kinematical data for this galaxy. Contrary to all previous studies, we avoid marginalization over quantities poorly constrained by observations or by theoretical arguments. We find minimal J-factors to be about 2 to 4 times smaller than commonly quoted estimates, approximately relaxing by the same amount the limit on Dark Matter pair annihilation cross section from gamma-ray surveys of Ursa Minor. At the same time, if one goes back to a fixed trial parametric form for the density, e.g. using a NFW or Burkert profile, we show that the minimal J can hardly be reduced by more than a factor of 1.5. \ua9 2016 IOP Publishing Ltd and Sissa Medialab srl

Vertical disc heating in Milky Way-sized galaxies in a cosmological context

Vertically extended, high velocity dispersion stellar distributions appear to be a ubiquitous feature of disc galaxies, and both internal and external mechanisms have been proposed to be the major driver of their formation. However, it is unclear to what extent each mechanism can generate such a distribution, which is likely to depend on the assembly history of the galaxy. To this end, we perform 16 high-resolution cosmological-zoom simulations of Milky Way-sized galaxies using the state-of-the-art cosmological magnetohydrodynamical code AREPO, and analyse the evolution of the vertical kinematics of the stellar disc in connection with various heating mechanisms. We find that the bar is the dominant heating mechanism in most cases, whereas spiral arms, radial migration and adiabatic heating from mid-plane density growth are all subdominant. The strongest source, though less prevalent than bars, originates from external perturbations from satellites/subhaloes of masses log10(M/M⊙) ≳ 10. However, in many simulations the orbits of newborn star particles become cooler with time, such that they dominate the shape of the age–velocity dispersion relation and overall vertical disc structure unless a strong external perturbation takes place

How sensitive are predicted galaxy luminosities to the choice of stellar population synthesis model?

We present a new release of the GALFORM semi-analytical model of galaxy formation and evolution, which exploits a Millennium Simulation-class N-body run performed with the Wilkinson Microwave Anisotropy Probe 7 cosmology. We use this new model to study the impact of the choice of stellar population synthesis (SPS) model on the predicted evolution of the galaxy luminosity function. The semi-analytical model is run using seven different SPS models. In each case, we obtain the rest-frame luminosity function in the far-ultraviolet, optical and near-infrared (NIR) wavelength ranges. We find that both the predicted rest-frame ultraviolet and optical luminosity function are insensitive to the choice of SPS model. However, we find that the predicted evolution of the rest-frame NIR luminosity function depends strongly on the treatment of the thermally pulsating asymptotic giant branch (TP-AGB) stellar phase in the SPS models, with differences larger than a factor of 2 for model galaxies brighter than MAB(K) − 5 log h < −22 (∼L* for 0 ≤ z ≤ 1.5). We have also explored the predicted number counts of galaxies, finding remarkable agreement between the results with different choices of SPS model, except when selecting galaxies with very red optical–NIR colours. The predicted number counts of these extremely red galaxies appear to be more affected by the treatment of star formation in discs than by the treatment of TP-AGB stars in the SPS models

Quenching and ram pressure stripping of simulated Milky Way satellite galaxies

We present predictions for the quenching of star formation in satellite galaxies of the Local Group from a suite of 30 cosmological zoom simulations of Milky Way-like host galaxies. The Auriga simulations resolve satellites down to the luminosity of the classical dwarf spheroidal galaxies of the Milky Way. We find strong mass-dependent and distance-dependent quenching signals, where dwarf systems beyond 600 kpc are only strongly quenched below a stellar mass of 107 M⊙. Ram pressure stripping appears to be the dominant quenching mechanism and 50% of quenched systems cease star formation within 1 Gyr of first infall. We demonstrate that systems within a host galaxy’s R200 radius are comprised of two populations: (i) a first infall population that has entered the host halo within the past few Gyrs and (ii) a population of returning ‘backsplash’ systems that have had a much more extended interaction with the host. Backsplash galaxies that do not return to the host galaxy by redshift zero exhibit quenching properties similar to galaxies within R200 and are distinct from other external systems. The simulated quenching trend with stellar mass has some tension with observations, but our simulations are able reproduce the range of quenching times measured from resolved stellar populations of Local Group dwarf galaxies

The Auriga Project: the properties and formation mechanisms of disc galaxies across cosmic time

We introduce a suite of 30 cosmological magneto-hydrodynamical zoom simulations of the formation of galaxies in isolated Milky Way mass dark haloes. These were carried out with the moving mesh code arepo, together with a comprehensive model for galaxy formation physics, including active galactic nuclei (AGN) feedback and magnetic fields, which produces realistic galaxy populations in large cosmological simulations. We demonstrate that our simulations reproduce a wide range of present-day observables, in particular, two-component disc-dominated galaxies with appropriate stellar masses, sizes, rotation curves, star formation rates and metallicities. We investigate the driving mechanisms that set present-day disc sizes/scalelengths, and find that they are related to the angular momentum of halo material. We show that the largest discs are produced by quiescent mergers that inspiral into the galaxy and deposit high-angular momentum material into the pre-existing disc, simultaneously increasing the spin of dark matter and gas in the halo. More violent mergers and strong AGN feedback play roles in limiting disc size by destroying pre-existing discs and by suppressing gas accretion on to the outer disc, respectively. The most important factor that leads to compact discs, however, is simply a low angular momentum for the halo. In these cases, AGN feedback plays an important role in limiting central star formation and the formation of a massive bulge

Magnetic field formation in the Milky Way like disc galaxies of the Auriga project

The magnetic fields observed in the Milky Way and nearby galaxies appear to be in equipartition with the turbulent, thermal and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogether. Here, we analyse the magnetic fields predicted by the simulations of the Auriga Project, a set of 30 high-resolution cosmological zoom simulations of Milky Way like galaxies, carried out with a moving-mesh magnetohydrodynamics code and a detailed galaxy formation physics model. We find that the magnetic fields grow exponentially at early times owing to a small-scale dynamo with an e-folding time of roughly 100 Myr in the centre of haloes until saturation occurs around z = 2–3, when the magnetic energy density reaches about 10 per cent of the turbulent energy density with a typical strength of 10--50μG 10--50μG . In the galactic centres, the ratio between magnetic and turbulent energies remains nearly constant until z = 0. At larger radii, differential rotation in the discs leads to linear amplification that typically saturates around z = 0.5–0. The final radial and vertical variations of the magnetic field strength can be well described by two joint exponential profiles, and are in good agreement with observational constraints. Overall, the magnetic fields have only little effect on the global evolution of the galaxies as it takes too long to reach equipartition. We also demonstrate that our results are well converged with numerical resolution

A new methodology to test galaxy formation models using the dependence of clustering on stellar mass

We present predictions for the two-point correlation function of galaxy clustering as a function of stellar mass, computed using two new versions of the GALFORM semi-analytic galaxy formation model. These models make use of a high resolution, large volume N-body simulation, set in the 7-year Wilkinson Microwave Anisotropy Probe cosmology. One model uses a universal stellar initial mass function (IMF), while the other assumes different IMFs for quiescent star formation and bursts. Particular consideration is given to how the assumptions required to estimate the stellar masses of observed galaxies (such as the choice of IMF, stellar population synthesis model, and dust extinction) influence the perceived dependence of galaxy clustering on stellar mass. Broad-band spectral energy distribution fitting is carried out to estimate stellar masses for the model galaxies in the same manner as in observational studies. We show clear differences between the clustering signals computed using the true and estimated model stellar masses. As such, we highlight the importance of applying our methodology to compare theoretical models to observations. We introduce an alternative scheme for the calculation of the merger time-scales for satellite galaxies in GALFORM, which takes into account the dark matter subhalo information from the simulation. This reduces the amplitude of small-scale clustering. The new merger scheme offers improved or similar agreement with observational clustering measurements, over the redshift range 0 < z < 0.7. We find reasonable agreement with clustering measurements from the Galaxy and Mass Assembly Survey, but find larger discrepancies for some stellar mass ranges and separation scales with respect to measurements from the Sloan Digital Sky Survey and the VIMOS Public Extragalactic Redshift Survey, depending on the GALFORM model used