The origin of the enhanced metallicity of satellite galaxies

Observations of galaxies in the local Universe have shown that both the ionized gas and the stars of satellites are more metal-rich than of equally massive centrals. To gain insight into the connection between this metallicity enhancement and other differences between centrals and satellites, such as their star formation rates, gas content, and growth history, we study the metallicities of >3600 galaxies with Mstar > 1010  M⊙ in the cosmological hydrodynamical EAGLE 100 Mpc ‘Reference’ simulation, including ∼1500 in the vicinity of galaxy groups and clusters (M200 ≥ 1013 M⊙). The simulation predicts excess gas and stellar metallicities in satellites consistent with observations, except for stellar metallicities at Mstar ≲ 1010.2 M⊙ where the predicted excess is smaller than observed. The exact magnitude of the effect depends on galaxy selection, aperture, and on whether the metallicity is weighted by stellar mass or luminosity. The stellar metallicity excess in clusters is also sensitive to the efficiency scaling of star formation feedback. We identify stripping of low-metallicity gas from the galaxy outskirts, as well as suppression of metal-poor inflows towards the galaxy centre, as key drivers of the enhancement of gas metallicity. Stellar metallicities in satellites are higher than in the field as a direct consequence of the more metal-rich star-forming gas, whereas stripping of stars and suppressed stellar mass growth, as well as differences in accreted versus in situ star formation between satellites and the field, are of secondary importance

Star formation quenching in simulated group and cluster galaxies: When, how, and why?

Star formation is observed to be suppressed in group and cluster galaxies compared to the field. To gain insight into the quenching process, we have analysed ~2000 galaxies formed in the GIMIC suite of cosmological hydrodynamical simulations. The time of quenching varies from ~2 Gyr before accretion (first crossing of r200,c) to >4 Gyr after, depending on satellite and host mass. Once begun, quenching is rapid (>~ 500 Myr) in low-mass galaxies (M* < 10^10 M_Sun), but significantly more protracted for more massive satellites. The simulations predict a substantial role of outflows driven by ram pressure -- but not tidal forces -- in removing the star-forming interstellar matter (ISM) from satellite galaxies, especially dwarfs (M* ~ 10^9 M_Sun) where they account for nearly two thirds of ISM loss in both groups and clusters. Immediately before quenching is complete, this fraction rises to ~80% even for Milky Way analogues (M* ~ 10^10.5 M_Sun) in groups (M_host ~ 10^13.5 M_Sun). We show that (i) ISM stripping was significantly more effective at early times than at z = 0; (ii) approximately half the gas is stripped from `galactic fountains' and half directly from the star forming disk; (iii) galaxies undergoing stripping experience ram pressure up to ~100 times the average at a given group/cluster-centric radius, because they are preferentially located in overdense ICM regions. Remarkably, stripping causes at most half the loss of the extended gas haloes surrounding our simulated satellites. These results contrast sharply with the current picture of strangulation -- removal of the ISM through star formation after stripping of the hot halo -- being the dominant mechanism quenching group and cluster satellites

The imprint of dark subhaloes on the circumgalactic medium

The standard model of cosmology, the LCDM model, robustly predicts the existence of a multitude of dark matter 'subhaloes' around galaxies like the Milky Way. A wide variety of observations have been proposed to look for the gravitational effects such subhaloes would induce in observable matter. Most of these approaches pertain to the stellar or cool gaseous phases of matter. Here we propose a new approach, which is to search for the perturbations that such dark subhaloes would source in the warm/hot circumgalactic medium (CGM) around normal galaxies. With a combination of analytic theory, carefully-controlled high-resolution idealised simulations, and full cosmological hydrodynamical simulations, we calculate the expected signal and how it depends on important physical parameters (subhalo mass, CGM temperature, and relative velocity). We find that dark subhaloes enhance the local CGM pressure, density, and temperature, in order of decreasing magnitude of the effect. For the pressure, the fluctuations can vary in magnitude from tens of percent (for subhaloes with M_sub=10^10 Msun) to a few percent (for subhaloes with M_sub=10^8 Msun), although this depends strongly on the CGM temperature. The subhaloes also induce fluctuations in the velocity field ranging in magnitude from a few km/s up to 25 km/s. We propose that X-ray, Sunyaev-Zel'dovich effect, radio dispersion measure, and quasar absorption line observations can be used to measure these fluctuations and place constraints on the abundance and distribution of dark subhaloes, thereby placing constraints on the nature of dark matter

Imprint of baryons and massive neutrinos on velocity statistics

We explore the impact of baryonic effects (namely stellar and AGN feedback) on the moments of pairwise velocity using the Illustris-TNG, EAGLE, cosmo-OWLS, and BAHAMAS suites of cosmological hydrodynamical simulations. The assumption that the mean pairwise velocity of the gas component follows that of the dark matter is studied here at small separations, and we find that even at pair separations of 10-20 $h^{-1}\mathrm{Mpc}$ there is a 4-5% velocity bias. At smaller separations, it gets larger with strength varying depending on the subgrid prescription. By isolating different physical processes, our findings suggest that the large scale velocity bias is mainly driven by stellar rather than AGN feedback. If unaccounted for, this velocity offset could possibly bias cosmological constraints from the kinetic Sunyaev-Zel'dovich effect in future cosmic microwave background (CMB) surveys. Furthermore, we examine how the first and the second moment of the pairwise velocity are affected by both the baryonic and the neutrino free-streaming effects for both the matter and gas components. For both moments, we were able to disentangle the effects of baryonic processes from those of massive neutrinos; and below pair separations of 20 $h^{-1}\mathrm{Mpc}$, we find that these moments of the pairwise velocity decrease with increasing neutrino mass. Our work thus paves a way in which the pairwise velocity statistics can be utilised to constrain the summed mass of neutrinos from future CMB surveys and peculiar velocity surveys

Exploring the effects of galaxy formation on matter clustering through a library of simulation power spectra

Upcoming weak lensing surveys require a detailed theoretical understanding of the matter power spectrum in order to derive accurate and precise cosmological parameter values. While galaxy formation is known to play an important role, its precise effects are currently unknown. We present a set of 92 matter power spectra from the OWLS, cosmo-OWLS and BAHAMAS simulation suites, including different $\Lambda$CDM cosmologies, neutrino masses, subgrid prescriptions and AGN feedback strengths. We conduct a detailed investigation of the dependence of the relative difference between the total matter power spectra in hydrodynamical and collisionless simulations on the effectiveness of stellar and AGN feedback, cosmology and redshift. The strength of AGN feedback can greatly affect the power on a range of scales, while a lack of stellar feedback can greatly increase the effectiveness of AGN feedback on large scales. We also examine differences in the initial conditions of hydrodynamic and N-body simulations that can lead to a ~1% discrepancy in the large-scale power, and furthermore show our results to be insensitive to cosmic variance. We present an empirical model capable of predicting the effect of galaxy formation on the matter power spectrum at z=0 to within 1% for k<1 h/Mpc, given only the mean baryon fraction in galaxy groups. Differences in group baryon fractions can also explain the quantitative disagreement between predictions from the literature. All total and dark matter only power spectra in this library will be made publicly available at powerlib.strw.leidenuniv.nl

A detection of wobbling brightest cluster galaxies within massive galaxy clusters

A striking signal of dark matter beyond the standard model is the existence of cores in the centre of galaxy clusters. Recent simulations predict that a brightest cluster galaxy (BCG) inside a cored galaxy cluster will exhibit residual wobbling due to previous major mergers, long after the relaxation of the overall cluster. This phenomenon is absent with standard cold dark matter where a cuspy density profile keeps a BCG tightly bound at the centre. We test this hypothesis using cosmological simulations and deep observations of 10 galaxy clusters acting as strong gravitational lenses. Modelling the BCG wobble as a simple harmonic oscillator, we measure the wobble amplitude, Aw, in the BAHAMAS suite of cosmological hydrodynamical simulations, finding an upper limit for the cold dark matter paradigm of Aw < 2 kpc at the 95 per cent confidence limit. We carry out the same test on the data finding a non-zero amplitude of Aw=11.82+7.3−3.0 kpc, with the observations dis-favouring Aw = 0 at the 3σ confidence level. This detection of BCG wobbling is evidence for a dark matter core at the heart of galaxy clusters. It also shows that strong lensing models of clusters cannot assume that the BCG is exactly coincident with the large-scale halo. While our small sample of galaxy clusters already indicates a non-zero Aw, with larger surveys, e.g. Euclid, we will be able to not only confirm the effect but also to use it to determine whether or not the wobbling finds its origin in new fundamental physics or astrophysical process

Weak lensing minima and peaks: Cosmological constraints and the impact of baryons

We present a novel statistic to extract cosmological information in weak lensing data: the lensing minima. We also investigate the effect of baryons on the cosmological constraints from peak and minimum counts. Using the \texttt{MassiveNuS} simulations, we find that lensing minima are sensitive to non-Gaussian cosmological information and are complementary to the lensing power spectrum and peak counts. For an LSST-like survey, we obtain $95\%$ credible intervals from a combination of lensing minima and peaks that are significantly stronger than from the power spectrum alone, by $44\%$, $11\%$, and $63\%$ for the neutrino mass sum $\sum m_\nu$, matter density $\Omega_m$, and amplitude of fluctuation $A_s$, respectively. We explore the effect of baryonic processes on lensing minima and peaks using the hydrodynamical simulations \texttt{BAHAMAS} and \texttt{Osato15}. We find that ignoring baryonic effects would lead to strong ($\approx 4 \sigma$) biases in inferences from peak counts, but negligible ($\approx 0.5 \sigma$) for minimum counts, suggesting lensing minima are a potentially more robust tool against baryonic effects. Finally, we demonstrate that the biases can in principle be mitigated without significantly degrading cosmological constraints when we model and marginalize the baryonic effects.UK Science and Technology Facilities Council (grant number ST/N000927/1)

On the relationship between cooling flows and bubbles

A common feature of the X-ray bubbles observed in Chandra images of some cooling flow clusters is that they appear to be surrounded by bright, cool shells. Temperature maps of a few nearby luminous clusters reveal that the shells consist of the coolest gas in the clusters—much cooler than the surrounding medium. Using simple models, we study the effects of this cool emission on the inferred cooling flow properties of clusters. We find that the introduction of bubbles into model clusters that do not have cooling flows results in temperature and surface brightness profiles that resemble those seen in nearby cooling flow clusters. They also approximately reproduce the recent XMM-Newton and Chandra observations of a high minimum temperature of ~1-3 keV. Hence, bubbles, if present, must be taken into account when inferring the physical properties of the intracluster medium. In the case of some clusters, bubbles may account entirely for these observed features, calling into question their designation as clusters with cooling flows. However, since not all nearby cooling flow clusters show bubble-like features, we suggest that there may be a diverse range of physical phenomena that give rise to the same observed features

Testing Sunyaev-Zel'dovich measurements of the hot gas content of dark matter haloes using synthetic skies

The thermal Sunyaev-Zel'dovich (tSZ) effect offers a means of probing the hot gas in and around massive galaxies and galaxy groups and clusters, which is thought to constitute a large fraction of the baryon content of the Universe. The Planck collaboration recently performed a stacking analysis of a large sample of `locally brightest galaxies' (LBGs) and, surprisingly, inferred an approximately self-similar relation between the tSZ flux and halo mass. At face value, this implies that the hot gas mass fraction is independent of halo mass, a result which is in apparent conflict with resolved X-ray observations. We test the robustness of the inferred trend using synthetic tSZ maps generated from cosmological hydrodynamical simulations and using the same tools and assumptions applied in the Planck study. We show that, while the detection and the estimate of the `total' flux (within $5 r_{500}$) is reasonably robust, the inferred flux originating from within $r_{500}$ (i.e. the limiting radius to which X-ray observations typically probe) is highly sensitive to the assumed pressure distribution of the gas. Using our most realistic simulations with AGN feedback, that reproduce a wide variety of X-ray and optical properties of groups and clusters, we estimate that the derived tSZ flux within $r_{500}$ is biased high by up to to an order of magnitude for haloes with masses $M_{500} \sim 10^{13}$ M$_{\odot}$. Moreover, we show that the AGN simulations are consistent with the total tSZ flux-mass relation observed with Planck, whereas a self-similar model is ruled out

The scatter and evolution of the global hot gas properties of simulated galaxy cluster populations

We use the cosmo-OWLS suite of cosmological hydrodynamical simulations to investigate the scatter and evolution of the global hot gas properties of large simulated populations of galaxy groups and clusters. Our aim is to compare the predictions of different physical models and to explore the extent to which commonly-adopted assumptions in observational analyses (e.g. self-similar evolution) are violated. We examine the relations between (true) halo mass and the X-ray temperature, X-ray luminosity, gas mass, Sunyaev-Zel'dovich (SZ) flux, the X-ray analogue of the SZ flux ($Y_X$) and the hydrostatic mass. For the most realistic models, which include AGN feedback, the slopes of the various mass-observable relations deviate substantially from the self-similar ones, particularly at late times and for low-mass clusters. The amplitude of the mass-temperature relation shows negative evolution with respect to the self-similar prediction (i.e. slower than the prediction) for all models, driven by an increase in non-thermal pressure support at higher redshifts. The AGN models predict strong positive evolution of the gas mass fractions at low halo masses. The SZ flux and $Y_X$ show positive evolution with respect to self-similarity at low mass but negative evolution at high mass. The scatter about the relations is well approximated by log-normal distributions, with widths that depend mildly on halo mass. The scatter decreases significantly with increasing redshift. The exception is the hydrostatic mass-halo mass relation, for which the scatter increases with redshift. Finally, we discuss the relative merits of various hot gas-based mass proxies