Enriching the hot circumgalactic medium

Simple models of galaxy formation in a cold darkmatter universe predict that massive galaxies are surrounded by a hot, quasi-hydrostatic circumgalactic corona of slowly cooling gas, predominantly accreted from the intergalactic medium (IGM). This prediction is borne out by the recent cosmological hydrodynamical simulations of Crain et al., which reproduce observed scaling relations between the X-ray and optical properties of nearby disc galaxies. Such coronae are metal poor, but observations of the X-ray emitting circumgalactic medium (CGM) of local galaxies typically indicate enrichment to near-solar iron abundance, potentially signalling a shortcoming in current models of galaxy formation.We show here that, while the hotCGMof galaxies formed in the simulations is typically metal poor in a mass-weighted sense, its X-ray luminosity-weighted metallicity is often close to solar. This bias arises because the soft X-ray emissivity of a typical ∼0.1 keV corona is dominated by collisionally excited metal ions that are synthesized in stars and recycled into the hot CGM. We find that these metals are ejected primarily by stars that form in situ to the main progenitor of the galaxy, rather than in satellites or external galaxies. The enrichment of the hot CGM therefore proceeds in an ‘inside–out’ fashion throughout the assembly of the galaxy: metals are transported from the central galaxy by supernova-driven winds and convection over several Gyr, establishing a strong negative radial metallicity gradient. Whilst metal ions synthesized by stars are necessary to produce the X-ray emissivity that enables the hot CGM of isolated galaxies to be detected with current instrumentation, the electrons that collisionally excite them are equally important. Since our simulations indicate that the electron density of hot coronae is dominated by the metal-poor gas accreted from the IGM, we infer that the hot CGM observed via X-ray emission is the outcome of both hierarchical accretion and stellar recycling

Alignments between galaxies, satellite systems and haloes

The spatial distribution of the satellite populations of the Milky Way and Andromeda are puzzling in that they are nearly perpendicular to the discs of their central galaxies. To understand the origin of such configurations we study the alignment of the central galaxy, satellite system and dark matter halo in the largest of the ‘Evolution and Assembly of GaLaxies and their Environments’ (EAGLE) simulation. We find that centrals and their satellite systems tend to be well aligned with their haloes, with a median misalignment angle of 33◦ in both cases. While the centrals are better aligned with the inner 10 kpc halo, the satellite systems are better aligned with the entire halo indicating that satellites preferentially trace the outer halo. The central–satellite alignment is weak (median misalignment angle of 52◦) and we find that around 20 per cent of systems have a misalignment angle larger than 78◦, which is the value for the Milky Way. The central–satellite alignment is a consequence of the tendency of both components to align with the dark matter halo. As a consequence, when the central is parallel to the satellite system, it also tends to be parallel to the halo. In contrast, if the central is perpendicular to the satellite system, as in the case of the Milky Way and Andromeda, then the central–halo alignment is much weaker. Dispersion-dominated (spheroidal) centrals have a stronger alignment with both their halo and their satellites than rotation-dominated (disc) centrals. We also found that the halo, the central galaxy and the satellite system tend to be aligned with the surrounding large-scale distribution of matter, with the halo being the better aligned of the three

The Milky Way's total satellite population and constraining the mass of the warm dark matter particle

The Milky Way’s (MW) satellite population is a powerful probe of warm dark matter (WDM) models as the abundance of small substructures is very sensitive to the properties of the WDM particle. However, only a partial census of the MW’s complement of satellite galaxies exists because surveys of the MW’s close environs are incomplete both in depth and in sky coverage. We present a new Bayesian analysis that combines the sample of satellites recently discovered by the Dark Energy Survey (DES) with those found in the Sloan Digital Sky Survey (SDSS) to estimate the total satellite galaxy luminosity function down to Mv = 0. We find that there should be at least 124+40−27124−27+40 (68% CL, statistical error) satellites as bright or brighter than Mv = 0 within 300 kpc of the Sun, with only a weak dependence on MW halo mass. When it comes online the Large Synoptic Survey Telescope should detect approximately half of this population. We also show that WDM models infer the same number of satellites as in ΛCDM, which will allow us to rule out those models that produce insufficient substructure to be viable

Colours and luminosities of z=0.1 simulated galaxies in the EAGLE simulations

We calculate the colours and luminosities of redshift z = 0.1 galaxies from the EAGLE simulation suite using the GALAXEV population synthesis models. We take into account obscuration by dust in birth clouds and diffuse ISM using a two-component screen model, following the prescription of Charlot and Fall. We compare models in which the dust optical depth is constant to models where it depends on gas metallicity, gas fraction and orientation. The colours of EAGLE galaxies for the more sophisticated models are in broad agreement with those of observed galaxies. In particular, EAGLE produces a red sequence of passive galaxies and a blue cloud of star forming galaxies, with approximately the correct fraction of galaxies in each population and with g-r colours within 0.1 magnitudes of those observed. Luminosity functions from UV to NIR wavelengths differ from observations at a level comparable to systematic shifts resulting from a choice between Petrosian and Kron photometric apertures. Despite the generally good agreement there are clear discrepancies with observations. The blue cloud of EAGLE galaxies extends to somewhat higher luminosities than in the data, consistent with the modest underestimate of the passive fraction in massive EAGLE galaxies. There is also a moderate excess of bright blue galaxies compared to observations. The overall level of agreement with the observed colour distribution suggests that EAGLE galaxies at z = 0.1 have ages, metallicities and levels of obscuration that are comparable to those of observed galaxies

The masses and density profiles of halos in a LCDM galaxy formation simulation

We investigate the internal structure and density profiles of halos of mass $10^{10}-10^{14}~M_\odot$ in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These follow the formation of galaxies in a $\Lambda$CDM Universe and include a treatment of the baryon physics thought to be relevant. The EAGLE simulations reproduce the observed present-day galaxy stellar mass function, as well as many other properties of the galaxy population as a function of time. We find significant differences between the masses of halos in the EAGLE simulations and in simulations that follow only the dark matter component. Nevertheless, halos are well described by the Navarro-Frenk-White (NFW) density profile at radii larger than ~5% of the virial radius but, closer to the centre, the presence of stars can produce cuspier profiles. Central enhancements in the total mass profile are most important in halos of mass $10^{12}-10^{13}M_\odot$, where the stellar fraction peaks. Over the radial range where they are well resolved, the resulting galaxy rotation curves are in very good agreement with observational data for galaxies with stellar mass $M_*<5\times10^{10}M_\odot$. We present an empirical fitting function that describes the total mass profiles and show that its parameters are strongly correlated with halo mass

The Fermi GeV excess: challenges for the dark matter interpretation

One of the most exciting recent results in the field of dark matter indirect searches has been the discovery of an excess emission in gamma rays from the Galactic centre above the standard astrophysical background. We show that current hydrodynamic simulations, namely simulated Milky Way-like galaxies within the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project, challenge the possibility to interpret the GeV excess as due to annihilation of dark matter particles in the halo if the Milky Way

Simulated Milky Way analogues: implications for dark matter direct searches

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 MilkyWay 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 _t 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

The impact of angular momentum on black hole accretion rates in simulations of galaxy formation

Feedback from energy liberated by gas accretion onto black holes (BHs) is an attractive mechanism to explain the exponential cut-off at the massive end of the galaxy stellar mass function (SMF). Semi-analytic models of galaxy formation in which this form of feedback is assumed to suppress cooling in haloes where the gas cooling time is large compared to the dynamical time do indeed achieve a good match to the observed SMF. Furthermore, hydrodynamic simulations of individual halos in which gas is assumed to accrete onto the central BH at the Bondi rate have shown that a self-regulating regime is established in which the BH grows just enough to liberate an amount of energy comparable to the thermal energy of the halo. However, this process is efficient at suppressing the growth not only of massive galaxies but also of galaxies like the Milky Way, leading to disagreement with the observed SMF. The Bondi accretion rate, however, is inappropriate when the accreting material has angular momentum. We present an improved accretion model that takes into account the circularisation and subsequent viscous transport of infalling material and include it as a "subgrid" model in hydrodynamic simulations of the evolution of halos with a wide range of masses. The resulting accretion rates are generally low in low mass (\lsim 10^{11.5} \msun) halos, but show outbursts of Eddington-limited accretion during galaxy mergers. During outbursts these objects strongly resemble quasars. In higher mass haloes, gas accretion occurs continuously, typically at $~10$ % of the Eddington rate, which is conducive to the formation of radio jets. The resulting dependence of the accretion behaviour on halo mass induces a break in the relation between galaxy stellar mass and halo mass in these simulations that matches observations

The properties of ”dark” ΛCDM halos in the Local Group

We examine the baryon content of low-mass Λ cold dark matter (ΛCDM) haloes (108 < M200/M⊙ < 5 × 109) using the APOSTLE cosmological hydrodynamical simulations. Most of these systems are free of stars and have a gaseous content set by the combined effects of cosmic reionization, which imposes a mass-dependent upper limit, and of ram-pressure stripping, which reduces it further in high-density regions. Haloes mainly affected by reionization (RELHICS; REionization-Limited H I Clouds) inhabit preferentially low-density regions and make up a population where the gas is in hydrostatic equilibrium with the dark matter potential and in thermal equilibrium with the ionizing UV background. Their thermodynamic properties are well specified, and their gas density and temperature profiles may be predicted in detail. Gas in RELHICS is nearly fully ionized but with neutral cores that span a large range of H I masses and column densities and have negligible non-thermal broadening. We present predictions for their characteristic sizes and central column densities; the massive tail of the distribution should be within reach of future blind H I surveys. Local Group RELHICS (LGRs) have some properties consistent with observed Ultra Compact High Velocity Clouds (UCHVCs) but the sheer number of the latter suggests that most UCHVCs are not RELHICS. Our results suggest that LGRs (i) should typically be beyond 500 kpc from the Milky Way or M31; (ii) have positive Galactocentric radial velocities; (iii) H I sizes not exceeding 1 kpc, and (iv) should be nearly round. The detection and characterization of RELHICS would offer a unique probe of the small-scale clustering of CDM