Source author record

Ian G. McCarthy

Ian G. McCarthy appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph.HE astro-ph

Catalog footprint

What is connected

60works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

High-resolution synthetic UV-submm images for Milky Way-mass simulated galaxies from the ARTEMIS project

We present redshift-zero synthetic dust-aware observations for the 45 Milky Way-mass simulated galaxies of the ARTEMIS project, calculated with the SKIRT radiative transfer code. The post-processing procedure includes components for star-forming regions, stellar sources, and diffuse dust. We produce and publicly release realistic high-resolution images for 50 commonly-used broadband filters from ultraviolet to sub-millimetre wavelengths and for 18 different viewing angles. We compare the simulated ARTEMIS galaxies to observed galaxies in the DustPedia database with similar stellar mass and star formation rate, and to synthetic observations of the simulated galaxies of the Auriga project produced in previous work using a similar post-processing technique. In all cases, global galaxy properties are derived using SED fitting. We find that, similar to Auriga, the post-processed ARTEMIS galaxies generally reproduce the observed scaling relations for global fluxes and physical properties, although dust extinction at FUV/UV wavelengths is underestimated and representative dust temperatures are lower than observed. At a resolved scale, we compare multi-wavelength non-parametric morphological properties of selected disc galaxies across the data sets. We find that the ARTEMIS galaxies largely reproduce the observed morphological trends as a function of wavelength, although they appear to be more clumpy and less symmetrical than observed. We note that the ARTEMIS and Auriga galaxies occupy adjacent regions in the specific star formation versus stellar mass plane, so that the synthetic observation data sets supplement each other.

preprint2022arXiv

Intrinsic alignments of the extended radio continuum emission of galaxies in the EAGLE simulations

We present measurements of the intrinsic alignments (IAs) of the star-forming gas of galaxies in the EAGLE simulations. Radio continuum imaging of this gas enables cosmic shear measurements complementary to optical surveys. We measure the orientation of star-forming gas with respect to the direction to, and orientation of, neighbouring galaxies. Star-forming gas exhibits a preferentially radial orientation-direction alignment that is a decreasing function of galaxy pair separation, but remains significant to $\gtrsim 1$ Mpc at $z=0$. The alignment is qualitatively similar to that exhibited by the stars, but is weaker at fixed separation. Pairs of galaxies hosted by more massive subhaloes exhibit stronger alignment at fixed separation, but the strong alignment of close pairs is dominated by ${\sim}L^\star$ galaxies and their satellites. At fixed comoving separation, the radial alignment is stronger at higher redshift. The orientation-orientation alignment is consistent with random at all separations, despite subhaloes exhibiting preferential parallel minor axis alignment. The weaker IA of star-forming gas than for stars stems from the former's tendency to be less well aligned with the dark matter structure of galaxies than the latter, and implies that the systematic uncertainty due to IA may be less severe in radio continuum weak lensing surveys than in optical counterparts. Alignment models equating the orientation of star-forming gas discs to that of stellar discs or the DM structure of host subhaloes will therefore overestimate the impact of IAs on radio continuum cosmic shear measurements.

preprint2022arXiv

Merger-induced galaxy transformations in the ARTEMIS simulations

Using the ARTEMIS set of 45 high-resolution cosmological simulations, we investigate a range of merger-induced dynamical transformations of Milky Way-like galaxies. We first identify populations of accreted stars on highly radial orbits, similar to the 'Gaia Sausage' in the Milky Way. We show that $\approx1/3$ of the ARTEMIS galaxies contain a similar feature, and confirm that they usually comprise stellar debris from the most massive accreted satellite. Selecting 15 galaxies with discs at the present-day, we study their changes around the times of the GS-like mergers. Dark matter haloes of many of these exhibit global changes in shape and orientation, with almost half becoming significantly more spherical when the mergers occur. Focusing on the galaxies themselves, we find that 4/15 have stellar discs which experience large changes in the orientation of their angular momentum (AM) axes, at rates of up to $\sim60$ degrees Gyr$^{-1}$. By calculating the orbital angular momentum axes of the satellites before they are accreted, we show that there is a tendency for the disc's AM to become more aligned with this axis after the merger. We also investigate the origin of in situ retrograde stars, analogous to the 'Splash' in the Milky Way. Tracing them back to earlier snapshots, we demonstrate that they were often disrupted onto their extreme orbits by multiple early mergers. We also find that the total mass of these stars outside the central regions positively correlates with the total accreted stellar mass.

preprint2022arXiv

Quenching of satellite galaxies of Milky Way analogues: reconciling theory and observations

The vast majority of low-mass satellite galaxies around the Milky Way and M31 appear virtually devoid of cool gas and show no signs of recent or ongoing star formation. Cosmological simulations demonstrate that such quenching is expected and is due to the harsh environmental conditions that satellites face when joining the Local Group (LG). However, recent observations of Milky Way analogues in the SAGA survey present a very different picture, showing the majority of observed satellites to be actively forming stars, calling into question the realism of current simulations and the typicality of the LG. Here we use the ARTEMIS suite of high-resolution cosmological hydrodynamical simulations to carry out a careful comparison with observations of dwarf satellites in the LG, SAGA, and the Local Volume (LV) survey. We show that differences between SAGA and the LG and LV surveys, as well as between SAGA and the ARTEMIS simulations, can be strongly reduced by considering differences in the host mass distributions and (more importantly) observational selection effects, specifically that low-mass satellites which have only recently been accreted are more likely to be star-forming, have a higher optical surface brightness, and are therefore more likely to be included in the SAGA survey. This picture is confirmed using data from the deeper LV survey, which shows pronounced quenching at low masses, in accordance with the predictions of LCDM-based simulations.

preprint2022arXiv

The morphology of star-forming gas and its alignment with galaxies and dark matter haloes in the EAGLE simulations

We present measurements of the morphology of star-forming gas in galaxies from the EAGLE simulations, and its alignment relative to stars and dark matter (DM). Imaging of such gas in the radio continuum enables weak lensing experiments that complement traditional optical approaches. Star-forming gas is typically more flattened than its associated stars and DM, particularly for present-day subhaloes of total mass $\sim$$10^{ 12-12.5} \mathrm{M_{ \odot}}$, which preferentially host star-forming galaxies with rotationally-supported stellar discs. Such systems have oblate, spheroidal star-forming gas distributions, but in both less- and more-massive subhaloes the distributions tend to be prolate, and its morphology correlates positively and significantly with that of its host galaxy's stars, both in terms of sphericity and triaxiality. The minor axis of star-forming gas most commonly aligns with the minor axis of its host subhalo's DM, but often aligns more closely with one of the other two principal axes of the DM distribution in prolate subhaloes. Star-forming gas aligns with DM less strongly than is the case for stars, but its morphological minor axis aligns closely with its kinematic axis, affording a route to observational identification of the unsheared morphological axis. The projected ellipticities of star-forming gas in EAGLE are consistent with shapes inferred from high-fidelity radio continuum images, and they exhibit greater shape noise than is the case for images of the stars, owing to the greater characteristic flattening of star-forming gas with respect to stars.

preprint2021arXiv

Constraining AGN feedback model with SZ profile

Relativistic jets from AGN have a wide range of impacts on galaxy groups and clusters and are key for understanding their formation and physical properties. However, this non-gravitational process is not well understood. Galaxy groups with shallow gravitational potentials are ideal laboratories to study and constrain the AGN feedback model. We studied hot gas in ~66,000 SDSS LRG halos with an average halo mass of 3 x 10^13 Msun using the Planck tSZ map. We have detected their average tSZ radial profile at ~17 sigma and compared it with the cosmo-OWLS cosmological hydrodynamical simulations with different AGN feedback models. The best agreement has been obtained for the AGN 8.0 model in the simulations. We have also compared our measured tSZ profile with the prediction from the universal pressure profile assuming the self-similar relation and found them consistent if the model accounts for the clustering of neighboring haloes via a two-halo term.

preprint2021arXiv

Galaxy Velocity Bias in Cosmological Simulations: Towards Percent-level Calibration

Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the ratio of galaxy and DM velocity dispersions -- as a function of redshift, host halo mass, and galaxy stellar mass threshold (Mstarsat), for massive halos (M200c > 1e13.5 msun) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of b_v increasing with M200c and decreasing with redshift and Mstarsat. The ensemble-estimated theoretical uncertainty in b_v is 2-3% but becomes percent-level when considering only the three highest resolution simulations. We update the mass-richness normalization previously estimated by Farahi et al. (2016) for an SDSS redMaPPer cluster sample. The improved accuracy of our b_v estimates reduces the mass normalization uncertainty from 22% to 8%, demonstrating that dynamical estimation techniques can be competitive with weak lensing in calibrating population mean masses. We discuss necessary steps for further improving this precision. Our estimates for b_v(M200c, Mstarsat, z) are made publicly available.

preprint2020arXiv

Connecting the structure of dark matter haloes to the primordial power spectrum

A large body of work based on collisionless cosmological N-body simulations going back over two decades has advanced the idea that collapsed dark matter haloes have simple and approximately universal forms for their mass density and pseudo-phase space density (PPSD) distributions. However, a general consensus on the physical origin of these results has not yet been reached. In the present study, we explore to what extent the apparent universality of these forms holds when we vary the initial conditions (i.e., the primordial power spectrum of density fluctuations) away from the standard CMB-normalised case, but still within the context of LCDM with a fixed expansion history. Using simulations that vary the initial amplitude and shape, we show that the structure of dark matter haloes retains a clear memory of the initial conditions. Specifically, increasing (lowering) the amplitude of fluctuations increases (decreases) the concentration of haloes and, if pushed far enough, the density profiles deviate strongly from the NFW form that is a good approximation for the CMB-normalised case. Although, an Einasto form works well. Rather than being universal, the slope of the PPSD (or pseudo-entropy) profile steepens (flattens) with increasing (decreasing) power spectrum amplitude and can exhibit a strong halo mass dependence. Our results therefore indicate that the previously identified universality of the structure of dark matter haloes is mostly a consequence of adopting a narrow range of (CMB-normalised) initial conditions for the simulations. Our new suite provides a useful test-bench against which physical models for the origin of halo structure can be validated.

preprint2020arXiv

Environment from cross-correlations: connecting hot gas and the quenching of galaxies

The observable properties of galaxies depend on both internal processes and the external environment. In terms of the environmental role, we still do not have a clear picture of the processes driving the transformation of galaxies. The use of proxies for environment (e.g., host halo mass, distance to the N^th nearest neighbour, etc.), as opposed to the real physical conditions (e.g., hot gas density) may bear some responsibility for this. Here we propose a new method that directly links galaxies to their local environments, by using spatial cross-correlations of galaxy catalogues with maps from large-scale structure surveys (e.g., thermal Sunyaev-Zel'dovich [tSZ] effect, diffuse X-ray emission, weak lensing of galaxies or the CMB). We focus here on the quenching of galaxies and its link to local hot gas properties. Maps of galaxy overdensity and quenched fraction excess are constructed from volume-limited SDSS catalogs, which are cross-correlated with tSZ effect and X-ray maps from Planck and ROSAT, respectively. Strong signals out to Mpc scales are detected for most cross-correlations and are compared to predictions from the EAGLE and BAHAMAS cosmological hydrodynamical simulations. The simulations successfully reproduce many, but not all, of the observed power spectra, with an indication that environmental quenching may be too efficient in the simulations. We demonstrate that the cross-correlations are sensitive to both the internal (e.g., AGN and stellar feedback) and external processes (e.g., ram pressure stripping, harassment, strangulation, etc.) responsible for quenching. The methods outlined in this paper can be adapted to other observables and, with upcoming surveys, will provide a stringent test of physical models for environmental transformation.

preprint2020arXiv

Exploring extensions to the standard cosmological model and the impact of baryons on small scales

It has been claimed that the standard model of cosmology (LCDM) cannot easily account for a number of observations on relatively small scales, motivating extensions to the standard model. Here we introduce a new suite of cosmological simulations that systematically explores three plausible extensions: warm dark matter, self-interacting dark matter, and a running of the scalar spectral index of density fluctuations. Current observational constraints are used to specify the additional parameters that come with these extensions. We examine a large range of observable metrics on small scales, including the halo mass function, density and circular velocity profiles, the abundance of satellite subhaloes, and halo concentrations. For any given metric, significant degeneracies can be present between the extensions. In detail, however, the different extensions have quantitatively distinct mass and radial dependencies, suggesting that a multi-probe approach over a range of scales can be used to break the degeneracies. We also demonstrate that the relative effects on the radial density profiles in the different extensions (compared to the standard model) are converged down to significantly smaller radii than are the absolute profiles. We compare the derived cosmological trends with the impact of baryonic physics using the EAGLE and ARTEMIS simulations. Significant degeneracies are also present between baryonic physics and cosmological variations (with both having similar magnitude effects on some observables). Given the inherent uncertainties both in the modelling of galaxy formation physics and extensions to LCDM, a systematic and simultaneous exploration of both is strongly warranted.

preprint2020arXiv

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.

preprint2020arXiv

Observable tests of self-interacting dark matter in galaxy clusters: cosmological simulations with SIDM and baryons

We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations. As is well known, SIDM haloes are systematically less dense in their centres, and rounder, than CDM haloes. Here we find that that these changes are not reflected in the distribution of gas or stars within galaxy clusters, or in their X-ray luminosities. However, gravitational lensing observables can discriminate between DM models, and we present a menu of tests that future surveys could use to measure the SIDM interaction strength. We ray-trace our simulated galaxy clusters to produce strong lensing maps. Including baryons boosts the lensing strength of clusters that produce no critical curves in SIDM-only simulations. Comparing the Einstein radii of our simulated clusters with those observed in the CLASH survey, we find that at velocities around 1000 km/s an SIDM cross-section of sigma/m > 1 cm^2/g is likely incompatible with observed cluster lensing.

preprint2020arXiv

Stellar Property Statistics of Massive Halos from Cosmological Hydrodynamics Simulations: Common Kernel Shapes

We study stellar property statistics, including satellite galaxy occupation, of massive halo populations realized by three cosmological hydrodynamics simulations: BAHAMAS + MACSIS, TNG300 of the IllustrisTNG suite, and Magneticum Pathfinder. The simulations incorporate independent sub-grid methods for astrophysical processes with spatial resolutions ranging from $1.5$ to $6$ kpc, and each generates samples of $1000$ or more halos with $M_{\rm halo}> 10^{13.5} M_{\odot}$ at redshift $z=0$. Applying localized, linear regression (LLR), we extract halo mass-conditioned statistics (normalizations, slopes, and intrinsic covariance) for a three-element stellar property vector consisting of: i) $N_{sat}$, the number of satellite galaxies with stellar mass, $M_{\star, \rm sat} > 10^{10} M_{\odot}$ within radius $R_{200c}$ of the halo; ii) $M_{\star,\rm tot}$, the total stellar mass within that radius, and; iii) $M_{\star,\rm BCG}$, the gravitationally-bound stellar mass of the central galaxy within a $100 \, \rm kpc$ radius. Scaling parameters for the three properties with halo mass show mild differences among the simulations, in part due to numerical resolution, but there is qualitative agreement on property correlations, with halos having smaller than average central galaxies tending to also have smaller total stellar mass and a larger number of satellite galaxies. Marginalizing over total halo mass, we find the satellite galaxy kernel, $p(\ln N_{sat}\,|\,M_{\rm halo},z)$ to be consistently skewed left, with skewness parameter $γ= -0.91 \pm 0.02$, while that of $\ln M_{\star,\rm tot}$ is closer to log-normal, in all three simulations. The highest resolution simulations find $γ\simeq -0.8$ for the $z=0$ shape of $p(\ln M_{\star,\rm BCG}\,|\,M_{\rm halo},z)$ and also that the fractional scatter in total stellar mass is below $10\%$ in halos more massive than $10^{14.3} M_{\odot}$.

preprint2020arXiv

The BAHAMAS project: Effects of a running scalar spectral index on large-scale structure

Recent analyses of the cosmic microwave background (CMB) and the Lyman-alpha forest indicate a mild preference for a deviation from a power law primordial matter power spectrum (a so-called negative `running'). We use an extension to the BAHAMAS suite of cosmological hydrodynamic simulations to explore the effects that a running scalar spectral index has on large-scale structure (LSS), using Planck CMB constraints to initialize the simulations. We focus on 5 key statistics: i) the non-linear matter power spectrum ii) the halo mass function; iii) the halo two-point auto correlation function; iv) total mass halo density profiles; and v) the halo concentration-mass relation. In terms of the matter power spectrum, we find that a running scalar spectral index affects all k-scales examined in this study, with a negative (positive) running leading to an amplification (suppression) of power. These effects should be easily detectable with upcoming surveys such as LSST and Euclid. In the mass range sampled, a positive running leads to an increase in the mass of galaxy groups and clusters, with the favoured negative running leading to a decrease in mass of lower-mass (M <~ 10^13 M_solar) halos, but an increase for the most massive (M >~ 10^13 M_solar) halos. Changes in the mass are generally confined to 5-10% which, while not insignificant, cannot by itself reconcile the claimed tension between the primary CMB and cluster number counts. We find that running does not significantly affect the shapes of density profiles of matched halos, changing only their amplitude. Finally, we demonstrate that the observed effects on LSS due to a running scalar spectral index are separable from those of baryonic effects to typically a few percent precision.

preprint2020arXiv

The BAHAMAS project: Effects of dynamical dark energy on large-scale structure

In this work we consider the impact of spatially-uniform but time-varying dark energy (or `dynamical dark energy', DDE) on large-scale structure in a spatially flat universe, using large cosmological hydrodynamical simulations that form part of the BAHAMAS project. As DDE changes the expansion history of the universe, it impacts the growth of structure. We explore variations in DDE that are constrained to be consistent with the cosmic microwave background. We find that DDE can affect the clustering of matter and haloes at the ~10% level (suppressing it for so-called `freezing' models, while enhancing it for `thawing' models), which should be distinguishable with upcoming large-scale structure surveys. DDE cosmologies can also enhance or suppress the halo mass function (with respect to LCDM) over a wide range of halo masses. The internal properties of haloes are minimally affected by changes in DDE, however. Finally, we show that the impact of baryons and associated feedback processes is largely independent of the change in cosmology and that these processes can be modelled separately to typically better than a few percent accuracy

preprint2020arXiv

The GOGREEN survey: Post-infall environmental quenching fails to predict the observed age difference between quiescent field and cluster galaxies at z>1

We study the star formation histories (SFHs) and mass-weighted ages of 331 UVJ-selected quiescent galaxies in 11 galaxy clusters and in the field at 1<z<1.5 from the Gemini Observations of Galaxies in Rich Early ENvironments (GOGREEN) survey. We determine the SFHs of individual galaxies by simultaneously fitting rest-frame optical spectroscopy and broadband photometry to stellar population models. We confirm that the SFHs are consistent with more massive galaxies having on average earlier formation times. Comparing galaxies found in massive clusters with those in the field, we find galaxies with $M_\ast<10^{11.3}$ M$_{\odot}$ in the field have more extended SFHs. From the SFHs we calculate the mass-weighted ages, and compare age distributions of galaxies between the two environments, at fixed mass. We constrain the difference in mass-weighted ages between field and cluster galaxies to $0.31_{^{-0.33}}^{_{+0.51}}$ Gyr, in the sense that cluster galaxies are older. We place this result in the context of two simple quenching models and show that neither environmental quenching based on time since infall (without pre-processing) nor a difference in formation times alone can reproduce both the average age difference and relative quenched fractions. This is distinctly different from local clusters, for which the majority of the quenched population is consistent with having been environmentally quenched upon infall. Our results suggest that quenched population in galaxy clusters at z>1 has been driven by different physical processes than those at play at z=0.

preprint2020arXiv

Weak lensing Analysis of X-Ray-selected XXL Galaxy Groups and Clusters with Subaru HSC Data

We present a weak-lensing analysis of X-ray galaxy groups and clusters selected from the XMM-XXL survey using the first-year data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. Our joint weak-lensing and X-ray analysis focuses on 136 spectroscopically confirmed X-ray-selected systems at 0.031 < z < 1.033 detected in the 25sqdeg XXL-N region. We characterize the mass distributions of individual clusters and establish the concentration-mass (c-M) relation for the XXL sample, by accounting for selection bias and statistical effects, and marginalizing over the remaining mass calibration uncertainty. We find the mass-trend parameter of the c-M relation to be β= -0.07 \pm 0.28 and the normalization to be c200 = 4.8 \pm 1.0 (stat) \pm 0.8 (syst) at M200=10^{14}Msun/h and z = 0.3. We find no statistical evidence for redshift evolution. Our weak-lensing results are in excellent agreement with dark-matter-only c-M relations calibrated for recent LCDM cosmologies. The level of intrinsic scatter in c200 is constrained as σ(\ln[c200]) < 24% (99.7% CL), which is smaller than predicted for the full population of LCDM halos. This is likely caused in part by the X-ray selection bias in terms of the relaxation state. We determine the temperature-mass (Tx-M500) relation for a subset of 105 XXL clusters that have both measured HSC lensing masses and X-ray temperatures. The resulting Tx-M500 relation is consistent with the self-similar prediction. Our Tx-M500 relation agrees with the XXL DR1 results at group scales, but has a slightly steeper mass trend, implying a smaller mass scale in the cluster regime. The overall offset in the Tx-M500 relation is at the $1.5σ$ level, corresponding to a mean mass offset of (34\pm 20)%. We also provide bias-corrected, weak-lensing-calibrated M200 and M500 mass estimates of individual XXL clusters based on their measured X-ray temperatures.

preprint2019arXiv

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 $Λ$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.

preprint2019arXiv

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_ν$, matter density $Ω_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 σ$) biases in inferences from peak counts, but negligible ($\approx 0.5 σ$) 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.

preprint2016arXiv

Dissecting the thermal Sunyaev-Zeldovich-gravitational lensing cross-correlation with hydrodynamical simulations

We use the cosmo-OWLS suite of cosmological hydrodynamical simulations, which includes different galactic feedback models, to predict the cross-correlation signal between weak gravitational lensing and the thermal Sunyaev-Zeldovich (tSZ) $y$-parameter. The predictions are compared to the recent detection reported by van Waerbeke and collaborators. The simulations reproduce the weak lensing-tSZ cross-correlation, $ξ_{yκ}(θ)$, well. The uncertainty arising from different possible feedback models appears to be important on small scales only ($θ\lesssim 10$ arcmin), while the amplitude of the correlation on all scales is sensitive to cosmological parameters that control the growth rate of structure (such as $σ_8$, $Ω_m$ and $Ω_b$). This study confirms our previous claim (in Ma et al.) that a significant proportion of the signal originates from the diffuse gas component in low-mass ($M_{\rm{halo}} \lesssim 10^{14} M_{\odot}$) clusters as well as from the region beyond the virial radius. We estimate that approximately 20$\%$ of the detected signal comes from low-mass clusters, which corresponds to about 30$\%$ of the baryon density of the Universe. The simulations also suggest that more than half of the baryons in the Universe are in the form of diffuse gas outside halos ($\gtrsim 5$ times the virial radius) which is not hot or dense enough to produce a significant tSZ signal or be observed by X-ray experiments. Finally, we show that future high-resolution tSZ-lensing cross-correlation observations will serve as a powerful tool for discriminating between different galactic feedback models.

preprint2016arXiv

nIFTY galaxy cluster simulations III: The Similarity & Diversity of Galaxies & Subhaloes

We examine subhaloes and galaxies residing in a simulated LCDM galaxy cluster ($M^{\rm crit}_{200}=1.1\times10^{15}M_\odot/h$) produced by hydrodynamical codes ranging from classic Smooth Particle Hydrodynamics (SPH), newer SPH codes, adaptive and moving mesh codes. These codes use subgrid models to capture galaxy formation physics. We compare how well these codes reproduce the same subhaloes/galaxies in gravity only, non-radiative hydrodynamics and full feedback physics runs by looking at the overall subhalo/galaxy distribution and on an individual objects basis. We find the subhalo population is reproduced to within $\lesssim10\%$ for both dark matter only and non-radiative runs, with individual objects showing code-to-code scatter of $\lesssim0.1$ dex, although the gas in non-radiative simulations shows significant scatter. Including feedback physics significantly increases the diversity. Subhalo mass and $V_{max}$ distributions vary by $\approx20\%$. The galaxy populations also show striking code-to-code variations. Although the Tully-Fisher relation is similar in almost all codes, the number of galaxies with $10^{9}M_\odot/h\lesssim M_*\lesssim 10^{12}M_\odot/h$ can differ by a factor of 4. Individual galaxies show code-to-code scatter of $\sim0.5$ dex in stellar mass. Moreover, strong systematic differences exist, with some codes producing galaxies $70\%$ smaller than others. The diversity partially arises from the inclusion/absence of AGN feedback. Our results combined with our companion papers demonstrate that subgrid physics is not just subject to fine-tuning, but the complexity of building galaxies in all environments remains a challenge. We argue even basic galaxy properties, such as the stellar mass to halo mass, should be treated with errors bars of $\sim0.2-0.4$ dex.

preprint2016arXiv

nIFTy Galaxy Cluster simulations IV: Quantifying the Influence of Baryons on Halo Properties

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.

preprint2016arXiv

nIFTy galaxy cluster simulations V: Investigation of the Cluster Infall Region

We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $Λ$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a dark-matter only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution $m_{\text{DM}}=9.01\times 10^8h^{-1}\text{M}_{\odot}$ and $m_{\text{gas}}=1.9\times 10^8h^{-1}\text{M}_{\odot}$ respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60% of haloes in the infall region with mass ~$10^{12.5} - 10^{14} h^{-1}\text{M}_{\odot}$, including 2-3 group-sized haloes ($> 10^{13}h^{-1}\text{M}_{\odot}$). However, we find that only ~10% of objects in the infall region are subhaloes residing in haloes, which may suggest that there is not much ongoing preprocessing occurring in the infall region at z=0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in subgrid schemes and calibration procedures that each model uses. Models that do not include AGN feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude.

preprint2016arXiv

The EAGLE simulations of galaxy formation: public release of halo and galaxy catalogues

We present the public data release of halo and galaxy catalogues extracted from the EAGLE suite of cosmological hydrodynamical simulations of galaxy formation. These simulations were performed with an enhanced version of the GADGET code that includes a modified hydrodynamics solver, time-step limiter and subgrid treatments of baryonic physics, such as stellar mass loss, element-by-element radiative cooling, star formation and feedback from star formation and black hole accretion. The simulation suite includes runs performed in volumes ranging from 25 to 100 comoving megaparsecs per side, with numerical resolution chosen to marginally resolve the Jeans mass of the gas at the star formation threshold. The free parameters of the subgrid models for feedback are calibrated to the redshift z=0 galaxy stellar mass function, galaxy sizes and black hole mass - stellar mass relation. The simulations have been shown to match a wide range of observations for present-day and higher-redshift galaxies. The raw particle data have been used to link galaxies across redshifts by creating merger trees. The indexing of the tree produces a simple way to connect a galaxy at one redshift to its progenitors at higher redshift and to identify its descendants at lower redshift. In this paper we present a relational database which we are making available for general use. A large number of properties of haloes and galaxies and their merger trees are stored in the database, including stellar masses, star formation rates, metallicities, photometric measurements and mock gri images. Complex queries can be created to explore the evolution of more than 10^5 galaxies, examples of which are provided in appendix. (abridged)

preprint2016arXiv

The EAGLE simulations: atomic hydrogen associated with galaxies

We examine the properties of atomic hydrogen (HI) associated with galaxies in the EAGLE simulations of galaxy formation. EAGLE's feedback parameters were calibrated to reproduce the stellar mass function and galaxy sizes at $z=0.1$, and we assess whether this calibration also yields realistic HI properties. We estimate the self-shielding density with a fitting function calibrated using radiation transport simulations, and correct for molecular hydrogen with empirical or theoretical relations. The `standard-resolution' simulations systematically underestimate HI column densities, leading to an HI deficiency in low-mass ($M_\star < 10^{10}M_\odot$) galaxies and poor reproduction of the observed HI mass function. These shortcomings are largely absent from EAGLE simulations featuring a factor of 8 (2) better mass (spatial) resolution, within which the HI mass of galaxies evolves more mildly from $z=1$ to $0$ than in the standard-resolution simulations. The largest-volume simulation reproduces the observed clustering of HI systems, and its dependence on HI-richness. At fixed $M_\star$, galaxies acquire more HI in simulations with stronger feedback, as they become associated with more massive haloes and higher infall rates. They acquire less HI in simulations with a greater star formation efficiency, since the star formation and feedback necessary to balance the infall rate is produced by smaller gas reservoirs. The simulations indicate that the HI of present-day galaxies was acquired primarily by the smooth accretion of ionized, intergalactic gas at $z\simeq1$, which later self-shields, and that only a small fraction is contributed by the reincorporation of gas previously heated strongly by feedback. HI reservoirs are highly dynamic: over $40$ percent of HI associated with $z=0.1$ galaxies is converted to stars or ejected by $z=0$.

preprint2016arXiv

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 M_star > 10^10 M_sun in the cosmological hydrodynamical EAGLE 100 Mpc `Reference' simulation, including ~1500 in the vicinity of galaxy groups and clusters (M_200 >= 10^13 M_sun). The simulation predicts excess gas and stellar metallicities in satellites consistent with observations, except for stellar metallicities at M_star <~ 10^10.2 M_sun 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 vs. in-situ star formation between satellites and the field, are of secondary importance.

preprint2016arXiv

The redshift evolution of massive galaxy clusters in the MACSIS simulations

We present the MAssive ClusterS and Intercluster Structures (MACSIS) project, a suite of 390 clusters simulated with baryonic physics that yields realistic massive galaxy clusters capable of matching a wide range of observed properties. MACSIS extends the recent BAHAMAS simulation to higher masses, enabling robust predictions for the redshift evolution of cluster properties and an assessment of the effect of selecting only the hottest systems. We study the observable-mass scaling relations and the X-ray luminosity-temperature relation over the complete observed cluster mass range. As expected, we find the slope of these scaling relations and the evolution of their normalization with redshift departs significantly from the self-similar predictions. However, for a sample of hot clusters with core-excised temperatures $k_{\rm{B}}T\geq5\,\rm{keV}$ the normalization and slope of the observable-mass relations and their evolution are significantly closer to self-similar. The exception is the temperature-mass relation, for which the increased importance of non-thermal pressure support and biased X-ray temperatures leads to a greater departure from self-similarity in the hottest systems. As a consequence, these also affect the slope and evolution of the normalization in the luminosity-temperature relation. The median hot gas profiles show good agreement with observational data at $z=0$ and $z=1$, with their evolution again departing significantly from the self-similar prediction. However, selecting a hot sample of clusters yields profiles that evolve significantly closer to the self-similar prediction. In conclusion, our results show that understanding the selection function is vital for robust calibration of cluster properties with mass and redshift.

preprint2016arXiv

The XMM Cluster Survey: evolution of the velocity dispersion -- temperature relation over half a Hubble time

We measure the evolution of the velocity dispersion--temperature ($σ_{\rm v}$--$T_{\rm X}$) relation up to $z = 1$ using a sample of 38 galaxy clusters drawn from the \textit{XMM} Cluster Survey. This work improves upon previous studies by the use of a homogeneous cluster sample and in terms of the number of high redshift clusters included. We present here new redshift and velocity dispersion measurements for 12 $z > 0.5$ clusters observed with the GMOS instruments on the Gemini telescopes. Using an orthogonal regression method, we find that the slope of the relation is steeper than that expected if clusters were self-similar, and that the evolution of the normalisation is slightly negative, but not significantly different from zero ($σ_{\rm v} \propto T^{0.86 \pm 0.14} E(z)^{-0.37 \pm 0.33}$). We verify our results by applying our methods to cosmological hydrodynamical simulations. The lack of evolution seen in our data is consistent with simulations that include both feedback and radiative cooling.

preprint2015arXiv

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.

preprint2015arXiv

nIFTy galaxy cluster simulations I: dark matter & non-radiative models

We have simulated the formation of a galaxy cluster in a $Λ$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods.

preprint2015arXiv

nIFTy galaxy cluster simulations II: radiative models

We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $Λ$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction.

preprint2015arXiv

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.

preprint2015arXiv

The alignment and shape of dark matter, stellar, and hot gas distributions in the EAGLE and cosmo-OWLS simulations

We report the alignment and shape of dark matter, stellar, and hot gas distributions in the EAGLE and cosmo-OWLS simulations. The combination of these state-of-the-art hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($11 < log_{10}(M_{200}/ [h^{-1}M_\odot]) < 15$), a wide radial range ($-2.3 < log_{10}(r/[h^{-1}Mpc ]) < 1.3$) and redshifts $0 < z < 1$. The shape parameters of the dark matter, stellar and hot gas distributions follow qualitatively similar trends: they become more aspherical (and triaxial) with increasing halo mass, radius and redshift. We measure the misalignment of the baryonic components (hot gas and stars) of galaxies with their host halo as a function of halo mass, radius, redshift, and galaxy type (centrals vs satellites and early- vs late-type). Overall, galaxies align well with the local distribution of the total (mostly dark) matter. However, the stellar distributions on galactic scales exhibit a median misalignment of about 45-50 degrees with respect to their host haloes. This misalignment is reduced to 25-30 degrees in the most massive haloes ($13 < log_{10}(M_{200}/ [h^{-1}M_\odot ]) < 15$). Half of the disc galaxies in the EAGLE simulations have a misalignment angle with respect to their host haloes larger than 40 degrees. We present fitting functions and tabulated values for the probability distribution of galaxy-halo misalignment to enable a straightforward inclusion of our results into models of galaxy formations based on purely collisionless N-body simulations.

preprint2015arXiv

The EAGLE simulations of galaxy formation: calibration of subgrid physics and model variations

We present results from thirteen cosmological simulations that explore the parameter space of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) simulation project. Four of the simulations follow the evolution of a periodic cube L = 50 cMpc on a side, and each employs a different subgrid model of the energetic feedback associated with star formation. The relevant parameters were adjusted so that the simulations each reproduce the observed galaxy stellar mass function at z = 0.1. Three of the simulations fail to form disc galaxies as extended as observed, and we show analytically that this is a consequence of numerical radiative losses that reduce the efficiency of stellar feedback in high-density gas. Such losses are greatly reduced in the fourth simulation - the EAGLE reference model - by injecting more energy in higher density gas. This model produces galaxies with the observed size distribution, and also reproduces many galaxy scaling relations. In the remaining nine simulations, a single parameter or process of the reference model was varied at a time. We find that the properties of galaxies with stellar mass <~ M* (the "knee" of the galaxy stellar mass function) are largely governed by feedback associated with star formation, while those of more massive galaxies are also controlled by feedback from accretion onto their central black holes. Both processes must be efficient in order to reproduce the observed galaxy population. In general, simulations that have been calibrated to reproduce the low-redshift galaxy stellar mass function will still not form realistic galaxies, but the additional requirement that galaxy sizes be acceptable leads to agreement with a large range of observables.

preprint2015arXiv

The EAGLE simulations of galaxy formation: the importance of the hydrodynamics scheme

We present results from a subset of simulations from the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) suite in which the formulation of the hydrodynamics scheme is varied. We compare simulations that use the same subgrid models without re-calibration of the parameters but employing the standard GADGET flavour of smoothed particle hydrodynamics (SPH) instead of the more recent state-of-the-art ANARCHY formulation of SPH that was used in the fiducial EAGLE runs. We find that the properties of most galaxies, including their masses and sizes, are not significantly affected by the details of the hydrodynamics solver. However, the star formation rates of the most massive objects are affected by the lack of phase mixing due to spurious surface tension in the simulation using standard SPH. This affects the efficiency with which AGN activity can quench star formation in these galaxies and it also leads to differences in the intragroup medium that affect the X-ray emission from these objects. The differences that can be attributed to the hydrodynamics solver are, however, likely to be less important at lower resolution. We also find that the use of a time step limiter is important for achieving the feedback efficiency required to match observations of the low-mass end of the galaxy stellar mass function.

preprint2015arXiv

The evolution of galaxy metallicity scaling relations in cosmological hydrodynamical simulations

The evolution of the metal content of galaxies and its relations to other global properties [such as total stellar mass (M*), circular velocity, star formation rate (SFR), halo mass, etc.] provides important constraints on models of galaxy formation. Here we examine the evolution of metallicity scaling relations of simulated galaxies in the Galaxies-Intergalactic Medium Interaction Calculation suite of cosmological simulations. We make comparisons to observations of the correlation of gas-phase abundances with M* (the mass-metallicity relation, MZR), as well as with both M* and SFR or gas mass fraction (the so-called 3D fundamental metallicity relations, FMRs). The simulated galaxies follow the observed local MZR and FMRs over an order of magnitude in M*, but overpredict the metallicity of massive galaxies (log M* > 10.5), plausibly due to inefficient feedback in this regime. We discuss the origin of the MZR and FMRs in the context of galactic outflows and gas accretion. We examine the evolution of mass-metallicity relations defined using different elements that probe the three enrichment channels (SNII, SNIa, and AGB stars). Relations based on elements produced mainly by SNII evolve weakly, whereas those based on elements produced preferentially in SNIa/AGB exhibit stronger evolution, due to the longer timescales associated with these channels. Finally, we compare the relations of central and satellite galaxies, finding systematically higher metallicities for satellites, as observed. We show this is due to the removal of the metal poor gas reservoir that normally surrounds galaxies and acts to dilute their gas-phase metallicity (via cooling/accretion onto the disk), but is lost due to ram pressure stripping for satellites.

preprint2015arXiv

The XXL Survey IV. Mass-temperature relation of the bright cluster sample

The XXL survey is the largest survey carried out by XMM-Newton. Covering an area of 50deg$^2$, the survey contains $\sim450$ galaxy clusters out to a redshift $\sim$2 and to an X-ray flux limit of $\sim5\times10^{-15}erg\,s^{-1}cm^{-2}$. This paper is part of the first release of XXL results focussed on the bright cluster sample. We investigate the scaling relation between weak-lensing mass and X-ray temperature for the brightest clusters in XXL. The scaling relation is used to estimate the mass of all 100 clusters in XXL-100-GC. Based on a subsample of 38 objects that lie within the intersection of the northern XXL field and the publicly available CFHTLenS catalog, we derive the $M_{WL}$ of each system with careful considerations of the systematics. The clusters lie at $0.1<z<0.6$ and span a range of $ T\simeq1-5keV$. We combine our sample with 58 clusters from the literature, increasing the range out to 10keV. To date, this is the largest sample of clusters with $M_{WL}$ measurements that has been used to study the mass-temperature relation. The fit ($M\propto T^b$) to the XXL clusters returns a slope $b=1.78^{+0.37}_{-0.32}$ and intrinsic scatter $σ_{\ln M|T}\simeq0.53$; the scatter is dominated by disturbed clusters. The fit to the combined sample of 96 clusters is in tension with self-similarity, $b=1.67\pm0.12$ and $σ_{\ln M|T}\simeq0.41$. Overall our results demonstrate the feasibility of ground-based weak-lensing scaling relation studies down to cool systems of $\sim1keV$ temperature and highlight that the current data and samples are a limit to our statistical precision. As such we are unable to determine whether the validity of hydrostatic equilibrium is a function of halo mass. An enlarged sample of cool systems, deeper weak-lensing data, and robust modelling of the selection function will help to explore these issues further.

preprint2014arXiv

Bent by baryons: the low mass galaxy-halo relation

The relation between galaxies and dark matter halos is of vital importance for evaluating theoretical predictions of structure formation and galaxy formation physics. We show that the widely used method of abundance matching based on dark matter only simulations fails at the low mass end because two of its underlying assumptions are broken: only a small fraction of low mass (below 10^9.5 solar masses) halos host a visible galaxy, and halos grow at a lower rate due to the effect of baryons. In this regime, reliance on dark matter only simulations for abundance matching is neither accurate nor self-consistent. We find that the reported discrepancy between observational estimates of the halo masses of dwarf galaxies and the values predicted by abundance matching does not point to a failure of LCDM, but simply to a failure to account for baryonic effects. Our results also imply that the Local Group contains only a few hundred observable galaxies in contrast with the thousands of faint dwarfs that abundance matching would suggest. We show how relations derived from abundance matching can be corrected, so that they can be used self-consistently to calibrate models of galaxy formation.

preprint2014arXiv

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.

preprint2014arXiv

The impact of baryonic processes on the two-point correlation functions of galaxies, subhaloes and matter

The observed clustering of galaxies and the cross-correlation of galaxies and mass provide important constraints on both cosmology and models of galaxy formation. Even though the dissipation and feedback processes associated with galaxy formation are thought to affect the distribution of matter, essentially all models used to predict clustering data are based on collisionless simulations. Here, we use large hydrodynamical simulations to investigate how galaxy formation affects the autocorrelation functions of galaxies and subhaloes, as well as their cross-correlation with matter. We show that the changes due to the inclusion of baryons are not limited to small scales and are even present in samples selected by subhalo mass. Samples selected by subhalo mass cluster ~10% more strongly in a baryonic run on scales r > 1Mpc/h, and this difference increases for smaller separations. While the inclusion of baryons boosts the clustering at fixed subhalo mass on all scales, the sign of the effect on the cross-correlation of subhaloes with matter can vary with radius. We show that the large-scale effects are due to the change in subhalo mass caused by the strong feedback associated with galaxy formation and may therefore not affect samples selected by number density. However, on scales r < r_vir significant differences remain after accounting for the change in subhalo mass. We conclude that predictions for galaxy-galaxy and galaxy-mass clustering from models based on collisionless simulations will have errors greater than 10% on sub-Mpc scales, unless the simulation results are modified to correctly account for the effects of baryons on the distributions of mass and satellites.

preprint2014arXiv

The impact of galaxy formation on the total mass, mass profile and abundance of haloes

We use cosmological hydrodynamical simulations to investigate how the inclusion of physical processes relevant to galaxy formation (star formation, metal-line cooling, stellar winds, supernovae and feedback from Active Galactic Nuclei, AGN) change the properties of haloes, over four orders of magnitude in mass. We find that gas expulsion and the associated dark matter (DM) expansion induced by supernova-driven winds are important for haloes with masses M200 < 10^13 Msun, lowering their masses by up to 20% relative to a DM-only model. AGN feedback, which is required to prevent overcooling, has a significant impact on halo masses all the way up to cluster scales (M200 ~ 10^15 Msun). Baryonic physics changes the total mass profiles of haloes out to several times the virial radius, a modification that cannot be captured by a change in the halo concentration. The decrease in the total halo mass causes a decrease in the halo mass function of about 20%. This effect can have important consequences for abundance matching technique as well as for most semi-analytic models of galaxy formation. We provide analytic fitting formulae, derived from simulations that reproduce the observed baryon fractions, to correct halo masses and mass functions from DM-only simulations. The effect of baryonic physics (AGN feedback in particular) on cluster number counts is about as large as changing the cosmology from WMAP7 to Planck, even when a moderately high mass limit of M500 ~ 10^14 Msun is adopted. Thus, for precision cosmology the effects of baryons must be accounted for.

preprint2014arXiv

The thermal Sunyaev Zel'dovich effect power spectrum in light of Planck

(Abridged) The amplitude of the thermal Sunyaev Zel'dovich effect (tSZ) power spectrum is extremely sensitive to the abundance of galaxy clusters and therefore to fundamental cosmological parameters that control their growth, such as sigma_8 and Omega_m. Here we explore the sensitivity of the tSZ power spectrum to important non-gravitational ('sub-grid') physics by employing the cosmo-OWLS suite of large-volume cosmological hydrodynamical simulations, run in both the Planck and WMAP7 best-fit cosmologies. On intermediate and small angular scales (ell > ~1000, or theta < ~10 arcmin), accessible with the South Pole Telescope and the Atacama Cosmology Telescope, the predicted tSZ power spectrum is highly model dependent, with AGN feedback having a particularly large effect. However, at large scales, observable with the Planck telescope, the effects of sub-grid physics are minor. Comparing the simulations with observations, we find a significant amplitude offset on all measured angular scales (including large scales), if the Planck best-fit cosmology is assumed by the simulations. This is shown to be a generic result for all current tSZ models. By contrast, if the WMAP7 cosmology is adopted, there is full consistency with the Planck power spectrum measurements on large scales and agreement at the 2 sigma level with the SPT/ACT measurements at intermediate scales for our fiducial AGN model, which Le Brun et al. (2014) have shown reproduces the 'resolved' properties of the local cluster population remarkably well. These findings strongly suggest that there are significantly fewer massive galaxy clusters than expected for the Planck best-fit cosmology, which is consistent with recent measurements of the tSZ number counts. Our findings therefore pose a significant challenge to the cosmological parameter values preferred (and/or the model adopted) by the Planck primary CMB analyses.

preprint2014arXiv

Towards a realistic population of simulated galaxy groups and clusters

We present a new suite of large-volume cosmological hydrodynamical simulations called cosmo-OWLS. They form an extension to the OverWhelmingly Large Simulations (OWLS) project, and have been designed to help improve our understanding of cluster astrophysics and non-linear structure formation, which are now the limiting systematic errors when using clusters as cosmological probes. Starting from identical initial conditions in either the Planck or WMAP7 cosmologies, we systematically vary the most important `sub-grid' physics, including feedback from supernovae and active galactic nuclei (AGN). We compare the properties of the simulated galaxy groups and clusters to a wide range of observational data, such as X-ray luminosity and temperature, gas mass fractions, entropy and density profiles, Sunyaev-Zel'dovich flux, I-band mass-to-light ratio, dominance of the brightest cluster galaxy, and central massive black hole (BH) masses, by producing synthetic observations and mimicking observational analysis techniques. These comparisons demonstrate that some AGN feedback models can produce a realistic population of galaxy groups and clusters, broadly reproducing both the median trend and, for the first time, the scatter in physical properties over approximately two decades in mass ($10^{13} \lesssim M_{500} \lesssim 10^{15}~\textrm{M}_{\odot}$) and 1.5 decades in radius ($0.05 \lesssim r/r_{500} \lesssim 1.5$). However, in other models, the AGN feedback is too violent (even though they reproduce the observed BH scaling relations), implying calibration of the models is required. The production of realistic populations of simulated groups and clusters, as well as models that bracket the observations, opens the door to the creation of synthetic surveys for assisting the astrophysical and cosmological interpretation of cluster surveys, as well as quantifying the impact of selection effects.

preprint2013arXiv

Enriching the hot circumgalactic medium

Models of galaxy formation in a CDM universe predict that massive galaxies are surrounded by a hot, quasi-hydrostatic circumgalactic corona of slowly cooling gas, predominantly accreted from the IGM. This prediction is borne out by the cosmological hydrodynamical simulations of Crain et al., which reproduce 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 galaxy formation models. We show here that, while the hot CGM of galaxies formed in the simulations is 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 synthesised 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 SNe-driven winds and convection over several gigayears, establishing a strong negative radial metallicity gradient. Whilst metal ions synthesised by stars are necessary to produce the X-ray emissivity that enables the hot CGM of isolated galaxies to be detected, 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.

preprint2013arXiv

Why does the environmental influence on group and cluster galaxies extend beyond the virial radius?

In the local Universe, galaxies in groups and clusters contain less gas and are less likely to be forming stars than their field counterparts. This effect is not limited to the central group/cluster regions, but is shown by recent observations to persist out to several virial radii. To gain insight into the extent and cause of this large-scale environmental influence, we use a suite of high-resolution cosmological hydrodynamic simulations to analyse galaxies around simulated groups and clusters of a wide range of mass (log M/M_sun = [13.0, 15.2]). In qualitative agreement with the observations, we find a systematic depletion of both hot and cold gas and a decline in the star forming fraction of galaxies as far out as ~ 5 r200 from the host centre. While a substantial fraction of these galaxies are on highly elliptical orbits and are not infalling for the first time (~ 50 per cent at 2 r200, independent of host mass) or are affected by `pre-processing' (less than 10 per cent of galaxies around groups, increasing to ~ 50 per cent around a massive cluster), even a combination of these indirect mechanisms does not fully account for the environmental influence, particularly in the case of the hot gas content. Direct ram pressure interaction with an extended gas `halo' surrounding groups and clusters is shown to be sufficiently strong to strip the hot gas atmospheres of infalling galaxies out to ~ 5 r200. We show that this influence is highly anisotropic, with ram pressure along filaments enhanced by up to a factor of 100 despite significant co-flow of gas and galaxies.

preprint2012arXiv

Mock weak lensing analysis of simulated galaxy clusters: bias and scatter in mass and concentration

(Abridged) We quantify the bias and scatter in galaxy cluster masses and concentrations derived from an idealised mock weak gravitational lensing (WL) survey, and their effect on the cluster mass-concentration relation. For this, we simulate WL distortions on a population of background galaxies due to a large (~3000) sample of galaxy cluster haloes extracted from the Millennium Simulation at z~0.2. This study takes into account the influence of shape noise, cluster substructure and asphericity as well as correlated large-scale structure, but not uncorrelated large-scale structure along the line of sight and observational effects. We find a small, but non-negligble, negative median bias in both mass and concentration at a level of ~5%, the exact value depending both on cluster mass and radial survey range. Both the mass and concentration derived from WL show considerable scatter about their true values. This scatter has, even for the highest mass clusters of M200 > 10^14.8 M_sun, a level of ~30% and ~20% for concentration and mass respectively and increases strongly with decreasing cluster mass. For a typical survey analysing 30 galaxies per arcmin^2 over a radial range from 30" to 15' from the cluster centre, the derived M200-c relation has a slope and normalisation too low compared to the underlying true (3D) relation by ~40% and ~15% respectively. The scatter and bias in mass are shown to reflect a departure at large radii of the true WL shear/matter distribution of the simulated clusters from the NFW profile adopted in modelling the mock observations. Orientation of the triaxial cluster haloes dominates the concentration scatter (except at low masses, where galaxy shape noise becomes dominant), while the bias in c is mostly due to substructure within the virial radius.

preprint2012arXiv

Quantifying the effect of baryon physics on weak lensing tomography

We use matter power spectra from cosmological hydrodynamic simulations to quantify the effect of baryon physics on the weak gravitational lensing shear signal. The simulations consider a number of processes, such as radiative cooling, star formation, supernovae and feedback from active galactic nuclei (AGN). Van Daalen et al. (2011) used the same simulations to show that baryon physics, in particular the strong feedback that is required to solve the overcooling problem, modifies the matter power spectrum on scales relevant for cosmological weak lensing studies. As a result, the use of power spectra from dark matter simulations can lead to significant biases in the inferred cosmological parameters. We show that the typical biases are much larger than the precision with which future missions aim to constrain the dark energy equation of state, w_0. For instance, the simulation with AGN feedback, which reproduces X-ray and optical properties of groups of galaxies, gives rise to a ~40% bias in w_0. We demonstrate that the modification of the power spectrum is dominated by groups and clusters of galaxies, the effect of which can be modelled. We consider an approach based on the popular halo model and show that simple modifications can capture the main features of baryonic feedback. Despite its simplicity, we find that our model, when calibrated on the simulations, is able to reduce the bias in w_0 to a level comparable to the size of the statistical uncertainties for a Euclid-like mission. While observations of the gas and stellar fractions as a function of halo mass can be used to calibrate the model, hydrodynamic simulations will likely still be needed to extend the observed scaling relations down to halo masses of 10 ^12 M_sun/h.

preprint2012arXiv

Rotation rates, sizes, and star formation efficiencies of a representative population of simulated disc galaxies

We examine the rotation rates, sizes, and star formation (SF) efficiencies of a representative population of simulated disc galaxies extracted from the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamic simulations. These simulations include efficient, but energetically feasible supernova feedback, but have not been tuned in any way to produce 'realistic' disc galaxies. Yet, they generate a large number of discs, without requiring extremely high resolution. Over the wide galaxy stellar mass range, 9.0 < log10[Mstar (Msun)] < 10.5, the simulations reproduce the observed Tully-Fisher relation, the rotation curves of disc galaxies in bins of stellar mass, the mass-size relation of disc galaxies, the optical rotation to virial circular velocity ratio (Vopt/Vvir), and the SF efficiencies of disc galaxies as inferred from stacked weak lensing and stacked satellite kinematics observations. They also reproduce the specific star formation rates of ~L* galaxies but predict too low levels of star formation for low-mass galaxies, which is plausibly due to the finite resolution of the simulations. At higher stellar masses, log10[Mstar (\Msun)] > 10.6, the simulated galaxies are too concentrated and have too high SF efficiencies. We conjecture that this shortcoming reflects the neglect of feedback from accreting supermassive black holes in these simulations. We conclude that it is possible to generate a representative population of disc galaxies that reproduces many of the observed trends of local disc galaxies using standard numerical hydrodynamic techniques and a plausible implementation of the "subgrid" astrophysical processes thought to be relevant to galaxy formation.

preprint2012arXiv

The competition between confinement and ram pressure and its implications for galaxies in groups and clusters

Ram pressure stripping of the hot gas that surrounds normal galaxies as they fall into groups and clusters (also referred to as `strangulation' or `starvation') is generally thought to shut down star formation on a time scale of a few Gyr. However, it has recently been suggested, on the basis of X-ray-optical scaling relations of galaxies in the field and the group/cluster environment, that confinement pressure by the intra-cluster medium can actually lead to an increase in the mass of hot gas surrounding these galaxies. We investigate the competition between pressure confinement and ram pressure stripping for satellite galaxies in orbit about galaxy groups and clusters using simple analytic models and detailed cosmological hydrodynamic simulations. It is found that, independent of host mass, ram pressure is generally dominant over confinement pressure --- only ~16 % of galaxies find themselves in the reverse situation. Furthermore, these galaxies have, on average, less hot gas than ram-pressure dominated ones, contrary to simple expectations. This is explained by the fact that the small number of galaxies which are confinement dominated are typically at first or second apocentre and have therefore already been maximally affected by ram pressure stripping around first pericentre. Our results are shown to be insensitive to host halo mass; we argue that the same is true for uncertain sub-grid processes, such as feedback.

preprint2012arXiv

The XMM Cluster Survey: The interplay between the brightest cluster galaxy and the intra-cluster medium via AGN feedback

Using a sample of 123 X-ray clusters and groups drawn from the XMM-Cluster Survey first data release, we investigate the interplay between the brightest cluster galaxy (BCG), its black hole, and the intra-cluster/group medium (ICM). It appears that for groups and clusters with a BCG likely to host significant AGN feedback, gas cooling dominates in those with Tx > 2 keV while AGN feedback dominates below. This may be understood through the sub-unity exponent found in the scaling relation we derive between the BCG mass and cluster mass over the halo mass range 10^13 < M500 < 10^15Msol and the lack of correlation between radio luminosity and cluster mass, such that BCG AGN in groups can have relatively more energetic influence on the ICM. The Lx - Tx relation for systems with the most massive BCGs, or those with BCGs co-located with the peak of the ICM emission, is steeper than that for those with the least massive and most offset, which instead follows self-similarity. This is evidence that a combination of central gas cooling and powerful, well fuelled AGN causes the departure of the ICM from pure gravitational heating, with the steepened relation crossing self-similarity at Tx = 2 keV. Importantly, regardless of their black hole mass, BCGs are more likely to host radio-loud AGN if they are in a massive cluster (Tx > 2 keV) and again co-located with an effective fuel supply of dense, cooling gas. This demonstrates that the most massive black holes appear to know more about their host cluster than they do about their host galaxy. The results lead us to propose a physically motivated, empirical definition of 'cluster' and 'group', delineated at 2 keV.

preprint2011arXiv

Cosmological simulations of the formation of the stellar haloes around disc galaxies

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the formation of stellar spheroids of Milky Way-mass disc galaxies. The simulations contain accurate treatments of metal-dependent radiative cooling, star formation, supernova feedback, and chemodynamics, and the large volumes that have been simulated yield an unprecedentedly large sample of ~400 simulated L_* disc galaxies. The simulated galaxies are surrounded by low-mass, low-surface brightness stellar haloes that extend out to ~100 kpc and beyond. The diffuse stellar distributions bear a remarkable resemblance to those observed around the Milky Way, M31 and other nearby galaxies, in terms of mass density, surface brightness, and metallicity profiles. We show that in situ star formation typically dominates the stellar spheroids by mass at radii of r < 30 kpc, whereas accretion of stars dominates at larger radii and this change in origin induces a change in slope of the surface brightness and metallicity profiles, which is also present in the observational data. The system-to-system scatter in the in situ mass fractions of the spheroid, however, is large and spans over a factor of 4. Consequently, there is a large degree of scatter in the shape and normalisation of the spheroid density profile within r < 30 kpc (e.g., when fit by a spherical powerlaw profile the indices range from -2.6 to -3.4). We show that the in situ mass fraction of the spheroid is linked to the formation epoch of the system. Dynamically older systems have, on average, larger contributions from in situ star formation, although there is significant system-to-system scatter in this relationship. Thus, in situ star formation likely represents the solution to the longstanding failure of pure accretion-based models to reproduce the observed properties of the inner spheroid.

preprint2011arXiv

Global structure and kinematics of stellar haloes in cosmological hydrodynamic simulations

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the global structure and kinematics of stellar spheroids of Milky Way mass disc galaxies. Font et al. have recently demonstrated that these simulations are able to successfully reproduce the satellite luminosity functions and the metallicity and surface brightness profiles of the spheroids of the Milky Way and M31. A key to the success of the simulations is a significant contribution to the spheroid from stars that formed in situ. While the outer halo is dominated by accreted stars, stars formed in the main progenitor of the galaxy dominate at r < ~30 kpc. In the present study we show that this component was primarily formed in a proto-disc at high redshift and was subsequently liberated from the disc by dynamical heating associated with mass accretion. As a consequence of its origin, the in situ component of the spheroid has different kinematics (namely net prograde rotation with respect to the disc) than that of the spheroid component built from the disruption of satellites. In addition, the in situ component has a flattened distribution, that is due in part to its rotation. We make comparisons with measurements of the shape and kinematics of local galaxies, including the Milky Way and M31, and stacked observations of more distant galaxies. We find that the simulated disc galaxies have spheroids of the correct shape (oblate with a median axis ratio of ~0.6 at radii of < ~30 kpc, but note there is significant system-to-system scatter in this quantity) and that the kinematics show evidence for two components (due to in situ vs. accreted), as observed. Our findings therefore add considerable weight to the importance of dissipative processes in the formation of stellar haloes and to the notion of a 'dual stellar halo'.

preprint2011arXiv

The population of Milky Way satellites in the LambdaCDM cosmology

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations.

preprint2010arXiv

Gas expulsion by quasar-driven winds as a solution to the over-cooling problem in galaxy groups and clusters

Galaxy groups are not scaled down versions of massive galaxy clusters - the hot gas in groups (known as the intragroup medium, IGrM for short) is, on average, less dense than the intracluster medium, implying that one or more non-gravitational processes (e.g., radiative cooling, star formation, and/or feedback) has had a relatively larger effect on groups. In the present study, we compare a number of cosmological hydrodynamic simulations that form part of the OverWhelmingly Large Simulations project to isolate and quantify the effects of cooling and feedback from supernovae (SNe) and active galactic nuclei (AGN) on the gas. This is achieved by comparing Lagrangian thermal histories of the gas in the different runs, which were all started from identical initial conditions. While radiative cooling, star formation, and SN feedback are all necessary ingredients, only runs that also include AGN feedback are able to successfully reproduce the optical and X-ray properties of groups and low-mass clusters. We isolate how, when, and exactly what gas is heated by AGN. Interestingly, we find that the gas that constitutes the present-day IGrM is that which was not strongly heated by AGN. Instead, the low median density/high median entropy of the gas in present-day groups is achieved by the ejection of lower entropy gas from low-mass progenitor galaxies at high redshift (primarily 2 < z < 4). This corresponds to the epoch when supermassive black holes accreted most of their mass, typically at a rate that is close to the Eddington limit (i.e., when the black holes are in a `quasar mode').

preprint2010arXiv

The case for AGN feedback in galaxy groups

[Abridged] The relatively recent insight that energy input from supermassive black holes (BHs) can have a substantial effect on the star formation rates (SFRs) of galaxies motivates us to examine its effects on the scale of galaxy groups. At present, groups contain most of the galaxies and a significant fraction of the overall baryon content of the universe. To explore the effects of BH feedback on groups, we analyse two high resolution cosmological hydro simulations from the OverWhelmingly Large Simulations project. While both include galactic winds driven by supernovae, only one includes feedback from BHs. We compare the properties of the simulated groups to a wide range of observational data, including hot gas radial profiles and gas mass fractions (fgas), luminosity-mass-temperature (L-M-T) scaling relations, K-band luminosity of the group and its central brightest galaxy (CBG), SFRs and ages of the CBG, and gas/stellar metallicities. Both runs yield entropy profiles similar to the data, while the run without AGN feedback yields highly peaked temperature profiles, in discord with the observations. Energy input from BHs significantly reduces fgas for groups with masses less than ~10^14 Msun, yielding fgas-T and L-T relations that are in agreement with the data. The run without AGN feedback suffers from the well known overcooling problem; the resulting K-band luminosities are much larger than observed. By contrast, the run that includes BH feedback yields K-band luminosities and CBG SFRs and ages in agreement with current estimates. Both runs yield very similar gas-phase metallicities that match X-ray data, but they predict very different stellar metallicities. Based on the above, galaxy groups provide a compelling case that BH feedback is a crucial ingredient in the formation of massive galaxies.

preprint2010arXiv

The similarity of observed X-ray coronae associated with L* disc and elliptical galaxies

The existence of hot, X-ray luminous gaseous coronae surrounding present day L* galaxies is a generic prediction of galaxy formation theory in the cold dark matter cosmogony. While extended X-ray emission has been known to exist around elliptical galaxies for a long time, diffuse extra-planar emission has only recently been detected around disc galaxies. We compile samples of elliptical and disc galaxies that have Chandra and XMM-Newton measurements, and compare the scaling of the coronal X-ray luminosity (L_X) with both the K-band luminosity (L_K) and the coronal X-ray temperature (T_X). The X-ray flux measurements are corrected for non-thermal point source contamination by spatial excision and spectral subtraction for resolved and unresolved sources respectively. We find that the properties of the extended X-ray emission from galaxies of different morphological types are similar: for both elliptical and disc galaxies, the L_X - L_K and L_X - T_X relations have similar slope, normalisation and scatter. The observed universality of coronal X-ray properties suggests that the bulk of this emission originates from gas that has been accreted, shock-heated and compressed during the assembly of the galaxy and that outflows triggered by stellar processes make only a minor contribution to the X-ray emission. This reservoir of cooling gas is a potential source of morphological transformation; it provides a fresh supply of material for discs to grow around galaxies of all morphological types.

preprint2010arXiv

X-ray coronae in simulations of disc galaxy formation

The existence of X-ray luminous gaseous coronae around massive disc galaxies is a long-standing prediction of galaxy formation theory in the cold dark matter cosmogony. This prediction has garnered little observational support, with non-detections commonplace and detections for only a relatively small number of galaxies which are much less luminous than expected. We investigate the coronal properties of a large sample of bright, disc-dominated galaxies extracted from the GIMIC suite of cosmological hydrodynamic simulations recently presented by Crain et al. Remarkably, the simulations reproduce the observed scalings of X-ray luminosity with K-band luminosity and star formation rate and, when account is taken of the density structure of the halo, with disc rotation velocity as well. Most of the star formation in the simulated galaxies (which have realistic stellar mass fractions) is fuelled by gas cooling from a quasi-hydrostatic hot corona. However, these coronae are more diffuse, and of a lower luminosity, than predicted by the analytic models of White & Frenk because of a substantial increase in entropy at z ~ 1-3. Both the removal of low entropy gas by star formation and energy injection from supernovae contribute to this increase in entropy, but the latter is dominant for halo masses M_200 <~ 10^(12.5) Msun. Only a small fraction of the mass of the hot gas is outflowing as a wind but, because of its high density and metallicity, it contributes disproportionally to the X-ray emission. The bulk of the X-ray emission, however, comes from the diffuse quasi-hydrostatic corona which supplies the fuel for ongoing star formation in discs today. Future deep X-ray observations with high spectral resolution (e.g. with NeXT/ASTRO-H or IXO) should be able to map the velocity structure of the hot gas and test this fundamental prediction of current galaxy formation theory.

preprint2009arXiv

Probing the cosmic web: inter-cluster filament detection using gravitational lensing

The problem of detecting dark matter filaments in the cosmic web is considered. Weak lensing is an ideal probe of dark matter, and therefore forms the basis of particularly promising detection methods. We consider and develop a number of weak lensing techniques that could be used to detect filaments in individual or stacked cluster fields, and apply them to synthetic lensing data sets in the fields of clusters from the Millennium Simulation. These techniques are multipole moments of the shear and convergence, mass reconstruction, and parameterized fits to filament mass profiles using a Markov Chain Monte Carlo approach. In particular, two new filament detection techniques are explored (multipole shear filters and Markov Chain Monte Carlo mass profile fits), and we outline the quality of data required to be able to identify and quantify filament profiles. We also consider the effects of large scale structure on filament detection. We conclude that using these techniques, there will be realistic prospects of detecting filaments in data from future space-based missions. The methods presented in this paper will be of great use in the identification of dark matter filaments in future surveys.

preprint2009arXiv

The accretion of galaxies into groups and clusters

We use the galaxy stellar mass and halo merger tree information from the semi-analytic model galaxy catalogue of Font et al. (2009) to examine the accretion of galaxies into a large sample of groups and clusters, covering a wide range in halo mass (10E12.9 to 10E15.3 Msun/h), and selected from each of four redshift epochs (z=0, 0.5, 1.0 and 1.5). We find that clusters at all examined redshifts have accreted a significant fraction of their final galaxy populations through galaxy groups. A 10E14.5 Msun/h mass cluster at z=0 has, on average, accreted ~ 40% of its galaxies (Mstellar > 10E9 Msun/h) from halos with masses greater than 10E13 Msun/h. Further, the galaxies which are accreted through groups are more massive, on average, than galaxies accreted through smaller halos or from the field population. We find that at a given epoch, the fraction of galaxies accreted from isolated environments is independent of the final cluster or group mass. In contrast, we find that observing a cluster of the same halo mass at each redshift epoch implies different accretion rates of isolated galaxies, from 5-6 % per Gyr at z=0 to 15% per Gyr at z=1.5. We find that combining the existence of a Butcher Oemler effect at z=0.5 and the observations that galaxies within groups display significant environmental effects with galaxy accretion histories justifies striking conclusions. Namely, that the dominant environmental process must begin to occur in halos of 10E12 -- 10E13 Msun/h, and act over timescales of > 2 Gyrs. This argues in favor of a mechanism like "strangulation", in which the hot halo of a galaxy is stripped upon infalling into a more massive halo . This simple model predicts that by z=1.5 galaxy groups and clusters will display little to no environmental effects.

preprint2008arXiv

The Colours of Satellite Galaxies in Groups and Clusters

Current models of galaxy formation predict satellite galaxies in groups and clusters that are redder than observed. We investigate the effect on the colours of satellite galaxies produced by the ram pressure stripping of their hot gaseous atmospheres as the satellites orbit within their parent halo. We incorporate a model of the stripping process based on detailed hydrodynamic simulations within the Durham semi-analytic model of galaxy formation. The simulations show that the environment in groups and clusters is less aggressive than previously assumed. The main uncertainty in the model is the treatment of gas expelled by supernovae. With reasonable assumptions for the stripping of this material, we find that satellite galaxies are able to retain a significant fraction of their hot gas for several Gigayears, thereby replenishing their reservoirs of cold, star forming gas and remaining blue for a relatively long period of time. A bimodal distribution of galaxy colours, similar to that observed in SDSS data, is established and the colours of the satellite galaxies are in good agreement with the data. In addition, our model naturally accounts for the observed dependence of satellite colours on environment, from small groups to high mass clusters.

Ian G. McCarthy

What is connected

Connect this record

See the researcher in context

Building this map preview

60 published item(s)

High-resolution synthetic UV-submm images for Milky Way-mass simulated galaxies from the ARTEMIS project

Intrinsic alignments of the extended radio continuum emission of galaxies in the EAGLE simulations

Merger-induced galaxy transformations in the ARTEMIS simulations

Quenching of satellite galaxies of Milky Way analogues: reconciling theory and observations

The morphology of star-forming gas and its alignment with galaxies and dark matter haloes in the EAGLE simulations

Constraining AGN feedback model with SZ profile

Galaxy Velocity Bias in Cosmological Simulations: Towards Percent-level Calibration

Connecting the structure of dark matter haloes to the primordial power spectrum

Environment from cross-correlations: connecting hot gas and the quenching of galaxies

Exploring extensions to the standard cosmological model and the impact of baryons on small scales

Imprint of baryons and massive neutrinos on velocity statistics

Observable tests of self-interacting dark matter in galaxy clusters: cosmological simulations with SIDM and baryons

Stellar Property Statistics of Massive Halos from Cosmological Hydrodynamics Simulations: Common Kernel Shapes

The BAHAMAS project: Effects of a running scalar spectral index on large-scale structure

The BAHAMAS project: Effects of dynamical dark energy on large-scale structure

The GOGREEN survey: Post-infall environmental quenching fails to predict the observed age difference between quiescent field and cluster galaxies at z>1

Weak lensing Analysis of X-Ray-selected XXL Galaxy Groups and Clusters with Subaru HSC Data

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

Weak Lensing Minima and Peaks: Cosmological Constraints and the Impact of Baryons

Dissecting the thermal Sunyaev-Zeldovich-gravitational lensing cross-correlation with hydrodynamical simulations

nIFTY galaxy cluster simulations III: The Similarity & Diversity of Galaxies & Subhaloes

nIFTy Galaxy Cluster simulations IV: Quantifying the Influence of Baryons on Halo Properties

nIFTy galaxy cluster simulations V: Investigation of the Cluster Infall Region

The EAGLE simulations of galaxy formation: public release of halo and galaxy catalogues

The EAGLE simulations: atomic hydrogen associated with galaxies

The origin of the enhanced metallicity of satellite galaxies

The redshift evolution of massive galaxy clusters in the MACSIS simulations

The XMM Cluster Survey: evolution of the velocity dispersion -- temperature relation over half a Hubble time

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

nIFTy galaxy cluster simulations I: dark matter & non-radiative models

nIFTy galaxy cluster simulations II: radiative models

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

The alignment and shape of dark matter, stellar, and hot gas distributions in the EAGLE and cosmo-OWLS simulations

The EAGLE simulations of galaxy formation: calibration of subgrid physics and model variations

The EAGLE simulations of galaxy formation: the importance of the hydrodynamics scheme

The evolution of galaxy metallicity scaling relations in cosmological hydrodynamical simulations

The XXL Survey IV. Mass-temperature relation of the bright cluster sample

Bent by baryons: the low mass galaxy-halo relation

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

The impact of baryonic processes on the two-point correlation functions of galaxies, subhaloes and matter

The impact of galaxy formation on the total mass, mass profile and abundance of haloes

The thermal Sunyaev Zel'dovich effect power spectrum in light of Planck

Towards a realistic population of simulated galaxy groups and clusters

Enriching the hot circumgalactic medium

Why does the environmental influence on group and cluster galaxies extend beyond the virial radius?

Mock weak lensing analysis of simulated galaxy clusters: bias and scatter in mass and concentration

Quantifying the effect of baryon physics on weak lensing tomography

Rotation rates, sizes, and star formation efficiencies of a representative population of simulated disc galaxies

The competition between confinement and ram pressure and its implications for galaxies in groups and clusters

The XMM Cluster Survey: The interplay between the brightest cluster galaxy and the intra-cluster medium via AGN feedback

Cosmological simulations of the formation of the stellar haloes around disc galaxies

Global structure and kinematics of stellar haloes in cosmological hydrodynamic simulations

The population of Milky Way satellites in the LambdaCDM cosmology

Gas expulsion by quasar-driven winds as a solution to the over-cooling problem in galaxy groups and clusters

The case for AGN feedback in galaxy groups

The similarity of observed X-ray coronae associated with L* disc and elliptical galaxies

X-ray coronae in simulations of disc galaxy formation

Probing the cosmic web: inter-cluster filament detection using gravitational lensing

The accretion of galaxies into groups and clusters

The Colours of Satellite Galaxies in Groups and Clusters